On Wed, Oct 7, 2015 at 2:03 PM, Cynthia Saad (Student) <crs07@mail.aub.edu> wrote:

Hello all, I am using enzo for 3d mhd collapse. I need to check conservation of divB=0. So i turned on BAnyl in the parameter file. Now i need to use yt for B_Div. And i'm having trouble doing so. I tried using phaseplot, but it says could not find field B_Div.

Can you try using the field name ('enzo', 'DivB'). If BAnyl is turned on, Enzo writes the divergence to disk as 'DivB'. You can also check what on-disk fields are recognized for your dataset by looking at ds.field_list.

One more question how may i load several data files in the same .py script?

Simply call the load function multiple times in the same script. You can have multiple datasets loaded at the same time. You might also be interested in the DatasetSeries class, which lets you load a time series or whole simulation of datasets. The easiest way to create a DatasetSeries is by passing a glob pattern to load, like in the following example: time_series = yt.load('DD????/DD????)' for ds in time_series: print ds 'ds' in this script is an individual simulation output. time_series is a DatasetSeries instance. hope that helps, Nathan

Thank you so much for taking the time to reply! Your input is really needed and mch appreciated.

Cynthia Saad _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org