<div dir="ltr">Hi John,<div><br></div><div>Yes, sorry it is vertex centered, I forgot to mention that. The data comes from a staggered grid however, (the velocities and pressures defined on separate grids); currently I'm just interpolating everything to the same grid. In principle I could interpolate everything to a cell centered grid (which is the native grid for the pressure), but I'd have to think about how to handle boundary conditions in that case. Does YT handle staggered grids in a natural way? </div><div><br></div><div>Thanks again,</div><div><br></div><div>Lukas</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 2, 2016 at 3:49 PM, John Zuhone <span dir="ltr"><<a href="mailto:jzuhone@gmail.com" target="_blank">jzuhone@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 style="word-wrap:break-word">Hi Lukas,<div><br></div><div>So I was actually able to convert your file to binary VTK and get started with it in yt, but I hit a roadblock. </div><div><br></div><div>Your data is “vertex-centered” in our parlance, but we don’t as of yet support vertex-centered data. </div><div><br></div><div>If you can convert to cell-centered data, we can probably support you out of the box, but for now vertex-centered data is not yet implemented. </div><div><br></div><div>For example, here is what the header for an Athena VTK file looks like:</div><div><br></div><div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span># vtk DataFile Version 3.0</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>PRIMITIVE vars at time= 2.700111e+03, level= 0, domain= 0</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>BINARY</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>DATASET STRUCTURED_POINTS</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>DIMENSIONS 257 257 257</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>ORIGIN -2.000000e+00 -2.000000e+00 -2.000000e+00</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>SPACING 1.562500e-02 1.562500e-02 1.562500e-02</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>CELL_DATA 16777216</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>SCALARS density float</span></div><div style="margin:0px;font-size:11px;line-height:normal;font-family:Menlo"><span>LOOKUP_TABLE default</span></div><div><span><br></span></div><div><span>Let us know…</span></div><div><span><br></span></div><div><span>John</span></div><div><span><br></span></div><div><blockquote type="cite"><div><div class="h5"><div>On May 16, 2016, at 10:04 PM, Lukas Bystricky <<a href="mailto:lukasbystricky@gmail.com" target="_blank">lukasbystricky@gmail.com</a>> wrote:</div><br></div></div><div><div><div class="h5"><div dir="ltr">Hi all,<div><br></div><div>I'm trying to load in some data from an ASCII VTK file (attached for reference) that I created following the guidelines here: <a href="http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf" target="_blank">http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf</a>. I can load this file into Visit and I was hoping that to load it into yt would be as simple as :</div><div><br></div><div>ds = yt.load("lid_driven_cavity_20.vtk")</div><div><br></div><div>When I try that however I get the error message:</div><div><div><br></div><div>Traceback (most recent call last):</div><div> File "visualization.py", line 3, in <module></div><div> ds = yt.load("lid_driven_cavity_20.vtk")</div><div> File "C:\Users\Dave\Anaconda2\lib\site-packages\yt\convenience.py", line 84, in load</div><div> return candidates[0](*args, **kwargs)</div><div> File "C:\Users\Dave\Anaconda2\lib\site-packages\yt\frontends\athena\data_structures.py", line 470, in __init__</div><div> Dataset.__init__(self, filename, dataset_type, units_override=units_override)</div><div> File "C:\Users\Dave\Anaconda2\lib\site-packages\yt\data_objects\static_output.py", line 190, in __init__</div><div> self._parse_parameter_file()</div><div> File "C:\Users\Dave\Anaconda2\lib\site-packages\yt\frontends\athena\data_structures.py", line 528, in _parse_parameter_file</div><div> self.domain_left_edge = grid['left_edge']</div><div>KeyError: 'left_edge'</div><div><br></div><div>Is it possible to load in such a file to yt? If not, since I'm writing my own data file, how could I edit it to make it readable by yt?</div></div><div><br></div><div>Thanks in advance,</div><div><br></div><div>Lukas</div></div>
</div></div><span><lid_driven_cavity_20.vtk></span>_______________________________________________<br>yt-users mailing list<br><a href="mailto:yt-users@lists.spacepope.org" target="_blank">yt-users@lists.spacepope.org</a><br><a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br></div></blockquote></div><br></div></div><br>_______________________________________________<br>
yt-users mailing list<br>
<a href="mailto:yt-users@lists.spacepope.org">yt-users@lists.spacepope.org</a><br>
<a href="http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org" rel="noreferrer" target="_blank">http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org</a><br>
<br></blockquote></div><br></div>