call Flame_rhJumpReactive ( real(inout) :: eosData_u(:), real(in) :: qbar_u, real(out) :: eosData_b(:), real(out) :: qbar_b, integer(in) :: eos_mode )
Calculate state of burned (and unburned) material by applying Rankine-Hugoniot jump condition. This version is for reactive ash (NSE), so that all information including the energy release (for a constant pressure burn) is calculated. Unburned state is found from calling Eos on eosData_u with mode eos_mode qbar values are in MeV/Baryon
eosData_u - thermodynamic state data for unburned material qbar_u - average binding energy per baryon (MeV) of unburned material eosData_b - thermodynamic state data for burned material qbar_u - average binding energy per baryon (MeV) of burned material eos_mode - equation of state mode for calculating complete unburned state.
See Flame_interface.F90 for possible updates