[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