<div dir="ltr"><div><div>Hi yt-users,<br><br></div>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.  <br><br>The relevant parts of my code look like this:<br><br>print(ds.field_list)<br><br></div><div>#only select the refined region of the box<br></div><div>ads = ds.r[0.42426:0.599523,0.450995:0.65,0.423461:0.590858]<br>adv = dv.r[0.373047:0.628906,0.369141:0.625,0.369141:0.662109]<br>ads.set_field_parameter('omega_baryon',0.046)<br>adv.set_field_parameter('omega_baryon',0.046)<br><br>plot = yt.PhasePlot(ads,"baryon_overdensity","temperature",["HI_Density"])<br><br>plot.save()<br><br>plot = yt.PhasePlot(adv,"baryon_overdensity","temperature",["HI_Density"])<br><br>plot.save()<br><br></div><div>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.  <br><br></div><div>Can anyone help me get cgs units in this phaseplot?<br><br></div><div>Thanks!<br><br></div><div>Stephanie<br><br></div><div><div><div><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div><div><font face="garamond,serif">--<br>Dr. Stephanie Tonnesen<br></font></div><font face="garamond,serif">Alvin E. Nashman Postdoctoral Fellow<br></font></div><font face="garamond,serif">Carnegie Observatories, Pasadena, CA<br></font></div><font face="garamond,serif"><a href="mailto:stonnes@gmail.com" target="_blank">stonnes@gmail.com</a><br></font></div></div></div>
</div></div></div></div>