[Paraview] XDMF and structured grid in cylindrical coordinates (curvilinear)

[David E DeMarle]

XDMF's 3DRectMesh == VTK's vtkRectilinearGrid, in which you case you
specify an X array, Y array and a Z array.

3DSMesh == vtkStructuredGrid, in which case you specify X,Y,Z for every
node. Once you do that you're example will work.

Examples of both (produced in ParaView via Data Object Generator (RG1 |
SG1)->Write as XDMF) are attached.



David E DeMarle
Kitware, Inc.
R&D Engineer
21 Corporate Drive
Clifton Park, NY 12065-8662
Phone: 518-881-4909


On Tue, Mar 19, 2013 at 11:04 PM, Nikolaos Beratlis <
nikos.beratlis at gmail.com> wrote:

> I am using a 3D structured orthogonal grid in cylindrical coordinates. I
> can read the grid as a cartesian grid using the following XDMF file by
> specifying TopologyType "3DRectMesh" and GeometryType="VXVYVZ" (see
> attached image):
>
> <?xml version="1.0" ?>
> <!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
> <Xdmf Version="2.0">
> <Domain>
> <Grid Name="mesh" GridType="Uniform">
> <Topology TopologyType="3DRectMesh" NumberOfElements="10 4 6"/>
> <Geometry GeometryType="VXVYVZ">
> <DataItem Dimensions="6" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.0 2.0 3.0 4.0 5.0
> </DataItem>
> <DataItem Dimensions="4" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.5708 3.1415 4.7124
> </DataItem>
> <DataItem Dimensions="10" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0
> </DataItem>
> </Geometry>
> </Grid>
> </Domain>
> </Xdmf>
>
> When I try to read the grid as 3DSMesh which I think corresponds to
> curvilinear Paraview crashes.
>
> <?xml version="1.0" ?>
> <!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
> <Xdmf Version="2.0">
> <Domain>
> <Grid Name="mesh" GridType="Uniform">
> <Topology TopologyType="3DSMesh" NumberOfElements="10 4 6"/>
> <Geometry GeometryType="VXVYVZ">
> <DataItem Dimensions="6" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.0 2.0 3.0 4.0 5.0
> </DataItem>
> <DataItem Dimensions="4" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.5708 3.1415 4.7124
> </DataItem>
> <DataItem Dimensions="10" NumberType="Float" Precision="4" Format="XML">
> 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0
> </DataItem>
> </Geometry>
> </Grid>
> </Domain>
> </Xdmf>
>
>  Is there any way to read a 3D cylindrical grid using three 1D arrays for
> the coordinates in each direction?
>
> Thank you,
>
> Nikos
>
>
>
> _______________________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Please keep messages on-topic and check the ParaView Wiki at:
> http://paraview.org/Wiki/ParaView
>
> Follow this link to subscribe/unsubscribe:
> http://www.paraview.org/mailman/listinfo/paraview
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20130320/fd76a37c/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testrg.xmf
Type: application/octet-stream
Size: 2211 bytes
Desc: not available
URL: <http://www.paraview.org/pipermail/paraview/attachments/20130320/fd76a37c/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: testsg.xmf
Type: application/octet-stream
Size: 2064 bytes
Desc: not available
URL: <http://www.paraview.org/pipermail/paraview/attachments/20130320/fd76a37c/attachment-0003.obj>