<div dir="ltr"><div><div>Hi, everyone:<br><br></div>I am still having problems outputting 2d VTK files. I have tried modifying the function vtk_out_hex() in postnek/trap.f (see attached file) by adding an extra conditional for when mz == 1 (which is true for 2d). In this case, I have simply copied the original code, but omitted the k loop and set k = 1.</div><div><br></div><div>When I do this, postnek will output the cell information for 2d simulations (CELLS), although it is still outputting the 8 points per cell (3d) since I am not sure which 4 of the 8 points might be correct for the 2d case. I have tried visualizing the 2d VTK file using various combinations of the 4 points: (1,2,3,4), (5,6,7,8), (1,3,5,7), etc., but the results still do not look correct.</div><div><br></div><div>I have attached the modified version of vtk_out_hex() that I have written. I have also attached the VTK output that it produces: eddy_uv_mod.vtk.</div><div><br></div><div>I'm afraid I am a bit stuck here without better knowledge of the underlying data structures in Nek. Any help or nudges in the correct direction are appreciated!</div><div><br></div><div>Thanks,</div><div>Taiyo<br></div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Dec 6, 2017 at 1:47 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</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"><div><div><div><div>Hi, all:<br><br></div>I have run into problems using postnek/postex to convert simulation output to VTK format for 2d simulations (3d simulations work fine). I have tested this using the eddy_uv example from the NekExamples directory. When using postex, I use the option "DUMP vtk ALL", which produces the attached VTK output. It seems that the grid point position information is incorrect.<br><br>Also, I have looked through the postnek functions in trap.f and it seems that VTK cell information (CELLS, CELL_TYPES) is omitted for 2d cases. Is there a way to output this information for 2d results?<br><br></div>Perhaps there is something I am overlooking? Any help is greatly appreciated.<br><br></div>Thanks,<br></div>Taiyo<br></div>
<br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div><div><div><div><div>Taiyo J. Wilson<br></div></div>Research Engineer<br></div>Illinois Rocstar LLC<br></div>108 Hessel Blvd.<br></div>Champaign, IL 61820<br></div><a href="mailto:twilson@illinoisrocstar.com" target="_blank">twilson@illinoisrocstar.com</a><br></div><a href="http://www.illinoisrocstar.com" target="_blank">www.illinoisrocstar.com</a><br></div></div></div></div>
</div>