[yt-users] Finding the max density from a checkpoint dump?

Nathan Goldbaum nathan12343 at gmail.com
Tue May 28 20:32:46 PDT 2013


Hi Joseph,

Excuse me - I saw in your paste that your yt python sources live in a yt-2.4 egg, but apparently that has nothing to do with the current version.

Looking at the code, it seems that somehow the 'x', 'y' and 'z' fields are different shapes for some of the grids in your data.  You can see for yourself what's going on by turning on the %debug magic in IPython and rerunning the command that crashes.  You'll be dumped into an ipdb prompt in the context where IndexError is raised.  You can see how 'x', 'y', and 'z' are defined by inspecting _coordX, _coordY, and _coordZ in $YT_HG/yt/data_objects/universal_fields.py.

I've never tried opening a checkpoint data dump in yt, although I suspect it's not done very often.  In the past I've used Enzo's subcycle-based data outputs to get intermediate interpolated data dumps in situations where I can't simply lower dtDatDump. There are many other enzo users on this list so perhaps one of them have had luck with this in the past.  

Sorry I can't be of more help here.

Cheers,

Nathan

On May 28, 2013, at 6:45 PM, Joseph Smidt <josephsmidt at gmail.com> wrote:

> Hey everyone,
> 
>     My Enzo simulation has refined to enough levels that I need to make checkpoint dumps mid-timestep and do YT analysis from there.  I can read in the checkpoint dump but when I go to find the max density I get this error:  http://pastebin.com/2mNVxNyW
> 
>    Anybody have any idea what is going on and how to fix this?  Also, since this dump is in the middle of the W-cycle, I believe the time on the finest grid is different then the root grid.  How do I find the time on the finest grid for this checkpoint dump?  Thanks.
> 
> -- 
> ------------------------------------------------------------------------
> Joseph Smidt <josephsmidt at gmail.com>
> 
> Theoretical Division
> P.O. Box 1663, Mail Stop B283
> Los Alamos, NM 87545
> Office: 505-665-9752
> Fax:    505-667-1931
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20130528/d3c106f2/attachment.html>


More information about the yt-users mailing list