[FLASH-USERS] A patch for a corrected version of hy_uhd_getFaceFlux

[dongwook at flash.uchicago.edu]

Dear FLASH users,

As James Guillochon kindly reported a bug in hy_uhd_getFaceFlux.F90
couple of days ago, I'd like to provide a corrected routine of it.

The bug was related to obtaining wrong states for several variables
including temperatures, ABAR, ZBAR and also mass fraction when calling
EOS.

The change only affects when "EOSforRiemann=.true." in the unsplit
hydro and MHD solvers. With EOSforRiemann=.true., the unsplit
solvers calculate gamc & game based on Riemann states of density and
pressure by calling extra EOS.

A default setup for EOSforRiemann is .false., in which gamc and game
are getting reconstructed from cell center values.

In order to use this new corrected routine, please replace your
old one (in FLASH3.3):

source/physics/Hydro/HydroMain/unsplit/hy_uhd_getFaceFlux.F90

with this new one.

Best regards,
Dongwook
-------------- next part --------------
A non-text attachment was scrubbed...
Name: hy_uhd_getFaceFlux.F90
Type: application/octet-stream
Size: 33664 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20110210/c6c239de/attachment-0001.obj>