The current version of the code has changed the temperature from computational units to Kelvin or solves the energy equation instead of temperature. This note is to check the parameter and equations changes need to do to be consistent comparing the old ways of blog:bliu01092014 and blog:bliu01142014 . So the term \gamma_{7} won't show in the thermal diffusion equations.

1. subroutine getdEdT in EOS.f90

This subroutine will return the value of dE/dT which has a unit as C.U. Energy/Kelvin

2. source and stencil_fix(0) in diffusion.f90 will be modified as

    source = ((ndiff-1.0)/(ndiff+1.0))*dot_product(T(0:2*nDim),stencil_fixed(0:2*nDim)) - T(0)*GetdEdT(Info%q(i,j,k,:))
    stencil_fixed(0) = stencil_fixed(0) - GetdEdT(Info%q(i,j,k,:)) 

comparing with the old way

!    source = ((ndiff-1.0)/(ndiff+1.0))*dot_product(T(0:2*nDim),stencil_fixed(0:2*nDim)) - T(0)*gamma7*Info%q(i,j,k,1)
!    stencil_fixed(0) = stencil_fixed(0) - gamma7*Info%q(i,j,k,1)