<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">XDMF&#39;s 3DRectMesh == VTK&#39;s vtkRectilinearGrid, in which you case you specify an X array, Y array and a Z array.</span><br><div><span style="font-family:arial,sans-serif;font-size:13px"><br>

</span></div><div style><span style="font-family:arial,sans-serif;font-size:13px">3DSMesh == vtkStructuredGrid, in which case you specify X,Y,Z for every node. Once you do that you&#39;re example will work.</span></div><div style>

<span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:13px">Examples of both (produced in ParaView via Data Object Generator (RG1 | SG1)-&gt;Write as XDMF) are attached.</span></div>

<div style><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div style><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div></div><div class="gmail_extra"><br clear="all"><div>

David E DeMarle<br>Kitware, Inc.<br>R&amp;D Engineer<br>21 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: 518-881-4909</div>
<br><br><div class="gmail_quote">On Tue, Mar 19, 2013 at 11:04 PM, Nikolaos Beratlis <span dir="ltr">&lt;<a href="mailto:nikos.beratlis@gmail.com" target="_blank">nikos.beratlis@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div dir="ltr">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 &quot;3DRectMesh&quot; and GeometryType=&quot;VXVYVZ&quot; (see attached image):<div>


<br></div><div><div>&lt;?xml version=&quot;1.0&quot; ?&gt;</div><div>&lt;!DOCTYPE Xdmf SYSTEM &quot;Xdmf.dtd&quot; []&gt;</div><div>&lt;Xdmf Version=&quot;2.0&quot;&gt;</div><div>&lt;Domain&gt;</div><div>&lt;Grid Name=&quot;mesh&quot; GridType=&quot;Uniform&quot;&gt;</div>


<div>&lt;Topology TopologyType=&quot;3DRectMesh&quot; NumberOfElements=&quot;10 4 6&quot;/&gt;</div><div>&lt;Geometry GeometryType=&quot;VXVYVZ&quot;&gt;</div><div>&lt;DataItem Dimensions=&quot;6&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div>


<div>0.0 1.0 2.0 3.0 4.0 5.0</div><div>&lt;/DataItem&gt;</div><div>&lt;DataItem Dimensions=&quot;4&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div><div>0.0 1.5708 3.1415 4.7124</div>


<div>&lt;/DataItem&gt;</div><div>&lt;DataItem Dimensions=&quot;10&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div><div>0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0</div><div>&lt;/DataItem&gt;</div>


<div>&lt;/Geometry&gt;</div><div>&lt;/Grid&gt;</div><div>&lt;/Domain&gt;</div><div>&lt;/Xdmf&gt;</div></div><div><br></div><div>When I try to read the grid as 3DSMesh which I think corresponds to curvilinear Paraview crashes.</div>


<div><br></div><div><div><div>&lt;?xml version=&quot;1.0&quot; ?&gt;</div><div>&lt;!DOCTYPE Xdmf SYSTEM &quot;Xdmf.dtd&quot; []&gt;</div><div>&lt;Xdmf Version=&quot;2.0&quot;&gt;</div><div>&lt;Domain&gt;</div>
<div>&lt;Grid Name=&quot;mesh&quot; GridType=&quot;Uniform&quot;&gt;</div><div>&lt;Topology TopologyType=&quot;3DSMesh&quot; NumberOfElements=&quot;10 4 6&quot;/&gt;</div><div>&lt;Geometry GeometryType=&quot;VXVYVZ&quot;&gt;</div>


<div>&lt;DataItem Dimensions=&quot;6&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div><div>0.0 1.0 2.0 3.0 4.0 5.0</div><div>&lt;/DataItem&gt;</div><div>&lt;DataItem Dimensions=&quot;4&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div>


<div>0.0 1.5708 3.1415 4.7124</div><div>&lt;/DataItem&gt;</div><div>&lt;DataItem Dimensions=&quot;10&quot; NumberType=&quot;Float&quot; Precision=&quot;4&quot; Format=&quot;XML&quot;&gt;</div><div>0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0</div>


<div>&lt;/DataItem&gt;</div><div>&lt;/Geometry&gt;</div><div>&lt;/Grid&gt;</div><div>&lt;/Domain&gt;</div><div>&lt;/Xdmf&gt;</div></div><div><br></div></div><div> Is there any way to read a 3D cylindrical grid using three 1D arrays for the coordinates in each direction?</div>


<div><br></div><div>Thank you,</div><div><br></div><div>Nikos</div><div><br></div><div><br></div></div>
<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.paraview.org/mailman/listinfo/paraview" target="_blank">http://www.paraview.org/mailman/listinfo/paraview</a><br>
<br></blockquote></div><br></div>