I accidentally committed more changes than intended.
The only other relevant change is to make Y_e_con_rhs public instead of private. This is in line with the other rhs variables, and is needed for refluxing.
-erik
On Wed, Sep 12, 2012 at 2:00 PM, schnetter@cct.lsu.edu wrote:
User: eschnett Date: 2012/09/12 01:00 PM
Modified: /trunk/ interface.ccl /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_TVDReconstruct.F90, GRHydro_Tmunu.F90, GRHydro_TmunuM.F90
Log: Correct format specifier
File Changes:
Directory: /trunk/src/
File [modified]: GRHydro_Con2Prim.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Con2Prim.F90 2012-09-04 01:20:40 UTC (rev 420) +++ trunk/src/GRHydro_Con2Prim.F90 2012-09-12 18:00:56 UTC (rev 421) @@ -859,7 +859,7 @@ if(nf.gt.nfudgemax) then if(GRHydro_c2p_reset_eps_tau_hot_eos.ne.1) then call CCTK_WARN(failinfomode,"EOS error in c2p 0: injected heat too many times")
write(warnline,"(i8.4i5,1P10E15.6)")cctk_iteration,GRHydro_reflevel,ii,jj,kk,x,y,z
write(warnline,"(i8,4i5,1P10E15.6)")cctk_iteration,GRHydro_reflevel,ii,jj,kk,x,y,z call CCTK_WARN(failinfomode,warnline) write(warnline,"(1P10E15.6)") rho,epsilon,temp,ye call CCTK_WARN(failinfomode,warnline)
File [modified]: GRHydro_TVDReconstruct.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_TVDReconstruct.F90 2012-09-04 01:20:40 UTC (rev 420) +++ trunk/src/GRHydro_TVDReconstruct.F90 2012-09-12 18:00:56 UTC (rev 421) @@ -115,7 +115,7 @@ CCTK_REAL,intent(IN)::a_in,b_in CCTK_REAL::minmod_result
minmod_result=0.5D0*(sign(1.0D0,a_in)+sign(1.0D0,b_in))*min(abs(a_in),abs(b_in))
minmod_result=0.5D0*(sign(1.0d0,a_in)+sign(1.0d0,b_in))*min(abs(a_in),abs(b_in)) end function minmod_func
end subroutine tvdreconstruct
File [modified]: GRHydro_Tmunu.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Tmunu.F90 2012-09-04 01:20:40 UTC (rev 420) +++ trunk/src/GRHydro_Tmunu.F90 2012-09-12 18:00:56 UTC (rev 421) @@ -105,7 +105,7 @@ !!$ Calculate the specific relativistic enthalpy times rho times the !!$ square of the lorentz factor.
rhoenthalpy = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0D0 +eps(i,j,k)) + press(i,j,k))
rhoenthalpy = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0d0 +eps(i,j,k)) + press(i,j,k))
!!$ Calculate lower components of 4-velocity (without the Lorent factor).
File [modified]: GRHydro_TmunuM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_TmunuM.F90 2012-09-04 01:20:40 UTC (rev 420) +++ trunk/src/GRHydro_TmunuM.F90 2012-09-12 18:00:56 UTC (rev 421) @@ -122,7 +122,7 @@ !!$ Calculate the specific relativistic enthalpy times rho + the mag. field contribution times the !!$ square of the lorentz factor.
rhohstarw2 = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0D0 +eps(i,j,k)) + press(i,j,k) + b2)
rhohstarw2 = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0d0 +eps(i,j,k)) + press(i,j,k) + b2) pstar = press(i,j,k)+0.5d0*b2
!!$ Calculate lower components of 4-velocity (without the Lorent factor).
Directory: /trunk/
File [modified]: interface.ccl Delta lines: +1 -1 =================================================================== --- trunk/interface.ccl 2012-09-04 01:20:40 UTC (rev 420) +++ trunk/interface.ccl 2012-09-12 18:00:56 UTC (rev 421) @@ -611,8 +611,8 @@ fs_alpha5 } "Maximum characteristic speeds for the flux splitting"
+public: real Y_e_con_rhs TYPE=GF tags='Prolongation="None" checkpoint="no"' "RHS for the electron fraction" -public: real Y_e_con_flux TYPE=GF tags='Prolongation="None" checkpoint="no"' "Flux for the electron fraction" private: real Y_e_plus TYPE=GF tags='Prolongation="None" checkpoint="no"' "Plus state for the electron fraction"
Commits mailing list Commits@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/commits