<div dir="ltr">I am currently trying the Function itemtype approach. In all of the examples for the function itemtype it assumes that arrays are all of the same size. In my case I have three 1D arrays of different sizes: x(nx), y(ny), and z(nz). The operation I need to do to convert to cylindrical coordinates is:<div>
<br></div><div>do k=1,nz</div><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 20, 2013 at 10:57 AM, David E DeMarle <span dir="ltr"><<a href="mailto:dave.demarle@kitware.com" target="_blank">dave.demarle@kitware.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">See if XDMF's Function itemtype will do the trick. (If so please consider posting an example to the <a href="http://xdmf.org" target="_blank">xdmf.org</a> wiki).<div>
<br><div><div>Otherwise you'll have to apply a filter to do the geometry transformation after you load the RectMesh into ParaView.</div>
</div><div><br></div></div></div><div class="gmail_extra"><div class="im"><br clear="all"><div>David E DeMarle<br>Kitware, Inc.<br>R&D Engineer<br>21 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: <a href="tel:518-881-4909" value="+15188814909" target="_blank">518-881-4909</a></div>
<br><br></div><div><div class="h5"><div class="gmail_quote">On Wed, Mar 20, 2013 at 10:52 AM, Nikolaos Beratlis <span dir="ltr"><<a href="mailto:nikos.beratlis@gmail.com" target="_blank">nikos.beratlis@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">Hi David,<div><br></div><div>I looked at your examples and I noticed that you changed the GeometryType to XYZ in order to plot a curvilinear grid. So it can't be done by specifying three separate 1D arrays instead (GeometryType=VXVYVZ)? My only concern is storage as I use grids with sizes of 2000x2000x2000. Writing that to a file will be huge compared to writing only three 1D arrays as I do for the cartesian case. Or is there a way to convert a grid read in cartesian coordinates to cylindrical the way it's done in Tecplot?</div>
<div><br></div><div>Thank you.</div><div><br></div><div>Nikos</div></div><div><div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Mar 20, 2013 at 10:18 AM, David E DeMarle <span dir="ltr"><<a href="mailto:dave.demarle@kitware.com" target="_blank">dave.demarle@kitware.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">XDMF's 3DRectMesh == VTK'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><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're example will work.</span></div><div>
<span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Examples of both (produced in ParaView via Data Object Generator (RG1 | SG1)->Write as XDMF) are attached.</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><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&D Engineer<br>21 Corporate Drive<br>Clifton Park, NY 12065-8662<br>Phone: <a href="tel:518-881-4909" value="+15188814909" target="_blank">518-881-4909</a></div>
<br><br><div class="gmail_quote"><div><div>On Tue, Mar 19, 2013 at 11:04 PM, Nikolaos Beratlis <span dir="ltr"><<a href="mailto:nikos.beratlis@gmail.com" target="_blank">nikos.beratlis@gmail.com</a>></span> wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
<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 "3DRectMesh" and GeometryType="VXVYVZ" (see attached image):<div>
<br></div><div><div><?xml version="1.0" ?></div><div><!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []></div><div><Xdmf Version="2.0"></div><div><Domain></div><div><Grid Name="mesh" GridType="Uniform"></div>
<div><Topology TopologyType="3DRectMesh" NumberOfElements="10 4 6"/></div><div><Geometry GeometryType="VXVYVZ"></div><div><DataItem Dimensions="6" NumberType="Float" Precision="4" Format="XML"></div>
<div>0.0 1.0 2.0 3.0 4.0 5.0</div><div></DataItem></div><div><DataItem Dimensions="4" NumberType="Float" Precision="4" Format="XML"></div><div>0.0 1.5708 3.1415 4.7124</div>
<div></DataItem></div><div><DataItem Dimensions="10" NumberType="Float" Precision="4" Format="XML"></div><div>0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0</div><div></DataItem></div>
<div></Geometry></div><div></Grid></div><div></Domain></div><div></Xdmf></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><?xml version="1.0" ?></div><div><!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []></div><div><Xdmf Version="2.0"></div><div><Domain></div>
<div><Grid Name="mesh" GridType="Uniform"></div><div><Topology TopologyType="3DSMesh" NumberOfElements="10 4 6"/></div><div><Geometry GeometryType="VXVYVZ"></div>
<div><DataItem Dimensions="6" NumberType="Float" Precision="4" Format="XML"></div><div>0.0 1.0 2.0 3.0 4.0 5.0</div><div></DataItem></div><div><DataItem Dimensions="4" NumberType="Float" Precision="4" Format="XML"></div>
<div>0.0 1.5708 3.1415 4.7124</div><div></DataItem></div><div><DataItem Dimensions="10" NumberType="Float" Precision="4" Format="XML"></div><div>0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0</div>
<div></DataItem></div><div></Geometry></div><div></Grid></div><div></Domain></div><div></Xdmf></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></div></div>_______________________________________________<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>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div>
</blockquote></div><br></div>