<div dir="ltr"><div><div><div><div><div><div><div>Hello, <br><br></div>I noticed a bug in the EoS subroutine of Flash 4.3. In the file eos_mgamma.F90 (source/physics/Eos/EosMain/Multigamma/eos_mgamma) line 377 the code is<br><br><br> eosData(det+ilo:det+ihi) = eos_gasConstant / eosData(abar+ilo:abar+ihi)* & <br> (eosData(gamc+ilo:gamc+ihi) - 1.0)<br> <br><br></div>Again, in line 429, <br><br> eosData(c_v+ilo:c_v+ihi) = eosData(det+ilo:det+ihi)<br><br></div>This
does not give the correct value of CV. There seems to be missing
bracket in calculation of eosData(det+ilo:det+ihi) which can be
corrected by writing code as:<br> <br> eosData(det+ilo:det+ihi) = eos_gasConstant / (eosData(abar+ilo:abar+ihi)* & <br> (eosData(gamc+ilo:gamc+ihi) - 1.0))<br><br></div>Thanks<br><br></div>Regards,<br></div>Jeevan Dahal<br></div>University of Missouri-Columbia</div>