Hi Geoffrey, Stephen, So there isn't currently a comoving density field. I think it can be added like so: def _ComovingDensity(field, data): ef = (1.0 + data.pf["CosmologyCurrentRedshift"])**3.0 return data["Density"]/ef add_field("ComovingDensity", function=_ComovingDensity) it won't need a conversion function because Density is already in CGS. Geoffrey, I think I'd prefer not to have a flag to swap things in and out -- globally turning on "comoving units" would likely end up having to add a lot of machinery to handle all of the comoving conversions. I think it's okay just to give them when requested. (Which, for what it's worth, is what Enzo does in the parameter file output.) -Matt On Sun, Feb 21, 2010 at 3:13 AM, <gso@physics.ucsd.edu> wrote:
Is it possible to leave the default in proper units, but add a new feature using a flag that will set the data in either the co-moving, or proper units? (or is it much easier for the user to multiply it by the right factors of 1+CurrentRedshift?)
From G.S.
Hi all,
I would like to make a phase plot that uses co-moving density, rather than proper density, something like this:
http://yt.enzotools.org/doc/cookbook/recipes.html#simple-phase
do I have to define a new field, or is there a built-in comoving density field (that I can't find)?
Thanks! _______________________________________________________ sskory@physics.ucsd.edu o__ Stephen Skory http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student ________________________________(_)_\(_)_______________
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org