<div dir="ltr"><div class="gmail_default" style=""><div class="gmail_default" style="">Hi,</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">I have been doing runs trying to simulate the evaporation of interstellar clouds via thermal conduction.  The runs have been working but the clouds are evaporating a bit faster than I expected.  This led me to look into the implementation of flux limitation.  Looking at diff_advanceTherm there is a call to Eos_getAbarZbar, which returns the value of Ye, but it's unclear what value will be returned for a single species run.  I have been setting eos_singleSpeciesA and that seems to work as expected.  I know in multispecies runs Ye = zbar/abar, so I wonder if its 1/abar.  Can anyone help me out on this?</div><div class="gmail_default" style=""><br></div><div class="gmail_default" style="">Regards,</div><div class="gmail_default" style="">Jon</div></div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">________________________________________________________<br>Jonathan D. Slavin                 Harvard-Smithsonian CfA<br><a href="mailto:jslavin@cfa.harvard.edu" target="_blank">jslavin@cfa.harvard.edu</a>       60 Garden Street, MS 83<br>phone: (617) 496-7981       Cambridge, MA 02138-1516<br>cell: (781) 363-0035             USA<br>________________________________________________________<br><br></div></div></div></div>
</div>