Thank you very much Stephen. I'll be waiting for Matt's reply. Also, thanks for solving the bug using SFR. Actually it was the next question I'll ask once I solve the "creation-time" problem. Regards, Fernando. On Jan 3, 2012, at 8:50 PM, Stephen Skory wrote:
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@skory.us http://stephenskory.com/ 510.621.3687 (google voice) _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org