In the first place, I do not understand the difference between 3DRectMesh
and 3DCoRectMesh.

I picked the following coordinate for testing:

Nz,Ny,Nx = [80, 90, 100]
z0,y0,x0 = [-3.,-2.,-1.]
lz,ly,lx = [Nz*3,Ny*2,Nx]
dz,dy,dx = [lz/(Nz-1), ly/(Ny-1), lx/(Nx-1)]

I found I have to use the following lines in the xdmf to get the correct
grid to be shown in Paraview for 3DRectMesh and 3DCoRectMesh, respectively:

<!--For 3DRectMesh -->
      <Topology TopologyType="3DRectMesh" Dimensions=" 80 90 100"/>
      <Geometry GeometryType="Origin_DXDYDZ">
        <DataItem Name="Origin" DataType="Float" Dimensions="3"
Format="XML"> -1 -2 -3</DataItem>
        <DataItem Name="DxDyDz" DataType="Float" Dimensions="3"
Format="XML"> 1 2 3</DataItem>
      </Geometry>

<!--For 3DCoRectMesh -->
      <Topology TopologyType="3DCoRectMesh" Dimensions=" 80 90 100"/>
      <Geometry GeometryType="Origin_DXDYDZ">
        <DataItem Name="Origin" DataType="Float" Dimensions="3"
Format="XML"> -3 -2 -1</DataItem>
        <DataItem Name="DxDyDz" DataType="Float" Dimensions="3"
Format="XML"> 3 2 1</DataItem>
      </Geometry>

Note that, the Dimensions for Topology in both cases are in zyx order,
while the Origin/DxDyDz arrays have different orders in the two cases:

3DRectMesh requires Origin=[x0,y0,z0], DX=[dx,dy,dz], i.e., xyz order
3DCoRectMesh requires Origin=[z0,y0,x0], DX=[dz,dy,dx], i.e., zyx order.

Am I writing the correct xdmf codes? How would you specify these data?
Thanks.
_______________________________________________
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

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/paraview

Reply via email to