
<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi,<div><br></div><div>There is a bug in hy_uhd_getFaceFlux relating to the "EOSforRiemann" flag. There is a variable eosMf that defaults to 1.0 and is never set afterwards, when it should be set to the cell averages. This means that as soon as there is any mixing this routine is not using the correct mass fractions when calling the equation of state.</div><div><br></div><div>The fix is easy, simply put this line above all 6 Eos calls (2 for each direction), replacing the indices with the appropriate ones that match that particular location in the code:</div><div>eosMf = 0.5*(U(SPECIES_BEGIN:SPECIES_END,i,j,k) + U(SPECIES_BEGIN:SPECIES_END,i,j,k+1))</div><div><br></div><div>Regards,</div><div>- James<br><div>

<div style="font-size: 12px; "><br class="Apple-interchange-newline">-- <br>James Guillochon<br>Department of Astronomy & Astrophysics<br>University of California, Santa Cruz<font class="Apple-style-span" color="#144FAE"><span class="Apple-style-span" style="text-decoration: underline; "><br></span></font></div><div style="font-size: 12px; "><a href="mailto:jfg@ucolick.org">jfg@ucolick.org</a></div>

</div>

<br></div></body></html>