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
>
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org