<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div dir="auto" style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Hi,<div class=""><br class=""></div><div class="">I’m looking at the Ionize module and suspect there might be a bug in <font face="Menlo" class="">IonizeMain/Eqi/nei.F90</font> (perhaps this file is out of date). I think the line</div><div class=""><br class=""></div><div class=""><font face="Menlo" class="">xout(jj2) = xout2(i)/selnel</font></div><div class=""><br class=""></div><div class="">should be</div><div class=""><br class=""></div><div class=""><font face="Menlo" class="">xout(jj2) = xout2(i) * selnel(jj2)</font></div><div class=""><br class=""></div><div class="">Even from the view of code grammar, <font face="Menlo" class="">selnel</font> is defined as a single real number, but it refers to an array <font face="Menlo" class="">ion_xfrac</font> when <font face="Menlo" class="">nei.F90</font> is called from <font face="Menlo" class="">Ioinze.F90</font>. In addition, population fraction (<span style="font-family: Menlo;" class="">xout2</span>) needs to be converted into mass fraction (<span style="font-family: Menlo;" class="">xout</span>), so the operation here should be multiplication instead of division.</div><div class=""><br class=""></div><div class="">And I have an additional question: if I understood this correctly, the electron population fraction is set to constant in the code (by electron to proton ratio). But even if we can ignore the contribution of metals, helium might not be fully ionized and may make a significant change to electron density. Then why is it valid that electron density remain unchanged throughout ionization?</div><div class=""><br class=""></div><div class="">Thanks in advance!</div><div class=""><br class=""><div class="">
<div class=""><div class="">Best wishes,</div><div class="">--</div><div class="">Suoqing JI</div><div class="">Ph.D Candidate</div><div class="">Department of Physics</div><div class="">University of California, Santa Barbara</div><div class=""><a href="http://web.physics.ucsb.edu/~suoqing" class="">http://web.physics.ucsb.edu/~suoqing</a></div></div>
</div>
<br class=""></div></div></div></div></body></html>