[yt-users] HI_Density units

[Stephanie Tonnesen]

Hi yt-users,

I am trying to get HI_Density in units of g/cm^3 but can only seem to get
it in terms of code mass/code length^3.

The relevant parts of my code look like this:

print(ds.field_list)

#only select the refined region of the box
ads = ds.r[0.42426:0.599523,0.450995:0.65,0.423461:0.590858]
adv = dv.r[0.373047:0.628906,0.369141:0.625,0.369141:0.662109]
ads.set_field_parameter('omega_baryon',0.046)
adv.set_field_parameter('omega_baryon',0.046)

plot = yt.PhasePlot(ads,"baryon_overdensity","temperature",["HI_Density"])

plot.save()

plot = yt.PhasePlot(adv,"baryon_overdensity","temperature",["HI_Density"])

plot.save()

Then I see that one of my fields is ('enzo','HI_Density') and I get a
phaseplot that looks good except HI_Density is in code mass/code length^3.
If I try "hi_density" I get that there is no field ('stars','hi_density').
If I try ('gas','hi_density') I get that there is no field
('stars',('gas','hi_density').  If I try "gas","hi_density" then I get that
there is no field 'stars','gas'.  I also tried to write
["HI_Density"].in_cgs() and that didn't work.  [("HI_Density").in_cgs()]
also didn't work.

Can anyone help me get cgs units in this phaseplot?

Thanks!

Stephanie

--
Dr. Stephanie Tonnesen
Alvin E. Nashman Postdoctoral Fellow
Carnegie Observatories, Pasadena, CA
stonnes at gmail.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160818/b9ebaa3d/attachment.htm>