Dear James,<br><br>Thanks for the correction. I will make this change in the routine.<br>I will also send out a follow-up email replying to your second email shortly.<br><br>Best regards,<br>Dongwook<br><br><div class="gmail_quote">
On Thu, Feb 3, 2011 at 5:32 PM, James Guillochon <span dir="ltr"><<a href="mailto:jfg@ucolick.org">jfg@ucolick.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div style="word-wrap: break-word;">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>-- <br>James Guillochon<br>Department of Astronomy & Astrophysics<br>University of California, Santa Cruz<font color="#144fae"><span style="text-decoration: underline;"><br></span></font></div>
<div style="font-size: 12px;"><a href="mailto:jfg@ucolick.org" target="_blank">jfg@ucolick.org</a></div>
</div>
<br></div></div></blockquote></div><br>