[FLASH-USERS] bug in eos_mgamma.F90 file

[jeevan dahal]

Hello,

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


         eosData(det+ilo:det+ihi) = eos_gasConstant /
eosData(abar+ilo:abar+ihi)* &
                (eosData(gamc+ilo:gamc+ihi) - 1.0)


Again, in line 429,

 eosData(c_v+ilo:c_v+ihi) = eosData(det+ilo:det+ihi)

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:

 eosData(det+ilo:det+ihi) = eos_gasConstant / (eosData(abar+ilo:abar+ihi)*
&
                (eosData(gamc+ilo:gamc+ihi) - 1.0))

Thanks

Regards,
Jeevan Dahal
University of Missouri-Columbia
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20151203/ad281e3b/attachment.htm>