[Users] [Commits] [svn:einsteintoolkit] GRHydro/trunk/ (Rev. 421)
Erik Schnetter
schnetter at cct.lsu.edu
Wed Sep 12 13:07:42 CDT 2012
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 at 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 at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/commits
>
--
Erik Schnetter <schnetter at cct.lsu.edu>
http://www.perimeterinstitute.ca/personal/eschnetter/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20120912/9ba51567/attachment.html
More information about the Users
mailing list