[yt-users] ionization states in ramses data sets

[Geoffrey So]

Hi,

according to
http://yt-project.org/docs/2.6/analyzing/creating_derived_fields.html?

what worked for me for enzo data is the following

def _HINeutralFraction(field, data):
    return data['HI_Density']/(data['HI_Density']+data['HII_Density'])

add_field('HINeutralFraction',function=_HINeutralFraction,
          units=r"\rho_\mathrm{HI}/\rho_\mathrm{H}")

unless there's something strange with ramses dataset, the above should work.

From
G.S.


On Thu, Aug 14, 2014 at 10:12 AM, Professor Maller <
amaller at citytech.cuny.edu> wrote:

> Hi,
>
> Does anyone have code to create derived fields of ionization states, e.g.
> HI, for ramses data sets?  If that's a no, does anyone have an example of a
> derived field where you have to read in a table and interpolate the field
> values.
>
> Thanks,
>
> Ari
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20140814/14caf322/attachment.htm>