[yt-users] Making slices of Athena SMR data

John ZuHone jzuhone at gmail.com
Mon Aug 5 14:22:12 PDT 2013


Hi Eve,

Sorry that you're having trouble!

Just wanted to check real quick before anything else--which version of yt are you using? Also, are you using the stable branch or the development branch? Support for Athena SMR exists in yt versions 2.5.3 and later. 

Best,

John

On Aug 5, 2013, at 5:13 PM, Eve Lee <elee at cita.utoronto.ca> wrote:

> Hello,
> 
> I'm trying to make a slice plot of Athena MPI SMR data using TimeSeriesData, and I find that the code only make slices out of the base level and blanks out deeper level regions.
> 
> Following is how I read in the Athena vtk files:
> 
> ts = TimeSeriesData.from_filenames('id0/kh*.vtk', parallel=True)
> my_storage={}
> for sto,pf in ts.piter(storage=my_storage):
>   slc = SlicePlot(pf, 'z', 'Density')
>   slc.set_zlim('Density', 0.0, 3.0)
>   num = (pf.h.hierarchy_filename).split('.')[1]
>   slc.save('./slice/'+num)
> 
> all my vtk files are in id0, ..., idn where n=number of processors I used, and I'm only using 1 more level of refinement.
> 
> Thanks in advance!
> 
> Best Regards,
> 
> Eve
> 
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org




More information about the yt-users mailing list