[yt-users] Star particles error

[Stephen Skory]

Fernando and Matt,

Fernando sent me one of his datasets and I think I know what's going
on. The dataset I'm looking at (g1e10s1e11dm1e12_0088) has 1449 grids,
but only 21 of them have (star) particles in them. Normally, when yt
first opens an enzo dataset and there is no .yt file in the directory,
it opens 20 semi-randomly chosen grids to see what fields exist in the
dataset. The assumption is that if these 20 grids are well-spaced
enough over all the levels, it should see all the fields in the
collective set of grids. From this, it builds a list of fields
(pf.h.field_list) and derived fields (pf.h.derived_field_list).

If I force yt to check all the grids (by modifiying
yt/frontends/enzo/data_structures.py) for fields on Fernando's
dataset, the particle fields work, and I got Fernando's SFR script at
the very top to work (modulo a bug I discovered in the cylinder object
volume calculation. I just pushed a fix and it should be in the
development copy of yt soon. Perhaps the same time as a fix for this
field list issue is also made available).

Matt - what do you think is the best solution for this? Add a ytcfg
parameter so users with datasets like this can simply tell yt to check
all grids for fields?

--
Stephen Skory
s at skory.us
http://stephenskory.com/
510.621.3687 (google voice)