I accidentally committed more changes than intended.<div><br></div><div>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.</div>
<div><br></div><div>-erik<br><br><div class="gmail_quote">On Wed, Sep 12, 2012 at 2:00 PM, <span dir="ltr"><<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
User: eschnett<br>
Date: 2012/09/12 01:00 PM<br>
<br>
Modified:<br>
/trunk/<br>
interface.ccl<br>
/trunk/src/<br>
GRHydro_Con2Prim.F90, GRHydro_TVDReconstruct.F90, GRHydro_Tmunu.F90, GRHydro_TmunuM.F90<br>
<br>
Log:<br>
Correct format specifier<br>
<br>
File Changes:<br>
<br>
Directory: /trunk/src/<br>
======================<br>
<br>
File [modified]: GRHydro_Con2Prim.F90<br>
Delta lines: +1 -1<br>
===================================================================<br>
--- trunk/src/GRHydro_Con2Prim.F90 2012-09-04 01:20:40 UTC (rev 420)<br>
+++ trunk/src/GRHydro_Con2Prim.F90 2012-09-12 18:00:56 UTC (rev 421)<br>
@@ -859,7 +859,7 @@<br>
if(nf.gt.nfudgemax) then<br>
if(GRHydro_c2p_reset_eps_tau_hot_eos.ne.1) then<br>
call CCTK_WARN(failinfomode,"EOS error in c2p 0: injected heat too many times")<br>
- write(warnline,"(i8.4i5,1P10E15.6)") cctk_iteration,GRHydro_reflevel,ii,jj,kk,x,y,z<br>
+ write(warnline,"(i8,4i5,1P10E15.6)") cctk_iteration,GRHydro_reflevel,ii,jj,kk,x,y,z<br>
call CCTK_WARN(failinfomode,warnline)<br>
write(warnline,"(1P10E15.6)") rho,epsilon,temp,ye<br>
call CCTK_WARN(failinfomode,warnline)<br>
<br>
File [modified]: GRHydro_TVDReconstruct.F90<br>
Delta lines: +1 -1<br>
===================================================================<br>
--- trunk/src/GRHydro_TVDReconstruct.F90 2012-09-04 01:20:40 UTC (rev 420)<br>
+++ trunk/src/GRHydro_TVDReconstruct.F90 2012-09-12 18:00:56 UTC (rev 421)<br>
@@ -115,7 +115,7 @@<br>
CCTK_REAL,intent(IN)::a_in,b_in<br>
CCTK_REAL::minmod_result<br>
<br>
- minmod_result=0.5D0*(sign(1.0D0,a_in)+sign(1.0D0,b_in))*min(abs(a_in),abs(b_in))<br>
+ minmod_result=0.5D0*(sign(1.0d0,a_in)+sign(1.0d0,b_in))*min(abs(a_in),abs(b_in))<br>
end function minmod_func<br>
<br>
end subroutine tvdreconstruct<br>
<br>
File [modified]: GRHydro_Tmunu.F90<br>
Delta lines: +1 -1<br>
===================================================================<br>
--- trunk/src/GRHydro_Tmunu.F90 2012-09-04 01:20:40 UTC (rev 420)<br>
+++ trunk/src/GRHydro_Tmunu.F90 2012-09-12 18:00:56 UTC (rev 421)<br>
@@ -105,7 +105,7 @@<br>
!!$ Calculate the specific relativistic enthalpy times rho times the<br>
!!$ square of the lorentz factor.<br>
<br>
- rhoenthalpy = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0D0 + eps(i,j,k)) + press(i,j,k))<br>
+ rhoenthalpy = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0d0 + eps(i,j,k)) + press(i,j,k))<br>
<br>
!!$ Calculate lower components of 4-velocity (without the Lorent factor).<br>
<br>
<br>
File [modified]: GRHydro_TmunuM.F90<br>
Delta lines: +1 -1<br>
===================================================================<br>
--- trunk/src/GRHydro_TmunuM.F90 2012-09-04 01:20:40 UTC (rev 420)<br>
+++ trunk/src/GRHydro_TmunuM.F90 2012-09-12 18:00:56 UTC (rev 421)<br>
@@ -122,7 +122,7 @@<br>
!!$ Calculate the specific relativistic enthalpy times rho + the mag. field contribution times the<br>
!!$ square of the lorentz factor.<br>
<br>
- rhohstarw2 = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0D0 + eps(i,j,k)) + press(i,j,k) + b2)<br>
+ rhohstarw2 = w_lorentz(i,j,k)**2*(rho(i,j,k)*(1.0d0 + eps(i,j,k)) + press(i,j,k) + b2)<br>
pstar = press(i,j,k)+0.5d0*b2<br>
<br>
!!$ Calculate lower components of 4-velocity (without the Lorent factor).<br>
<br>
Directory: /trunk/<br>
==================<br>
<br>
File [modified]: interface.ccl<br>
Delta lines: +1 -1<br>
===================================================================<br>
--- trunk/interface.ccl 2012-09-04 01:20:40 UTC (rev 420)<br>
+++ trunk/interface.ccl 2012-09-12 18:00:56 UTC (rev 421)<br>
@@ -611,8 +611,8 @@<br>
fs_alpha5<br>
} "Maximum characteristic speeds for the flux splitting"<br>
<br>
+public:<br>
real Y_e_con_rhs TYPE=GF tags='Prolongation="None" checkpoint="no"' "RHS for the electron fraction"<br>
-public:<br>
real Y_e_con_flux TYPE=GF tags='Prolongation="None" checkpoint="no"' "Flux for the electron fraction"<br>
private:<br>
real Y_e_plus TYPE=GF tags='Prolongation="None" checkpoint="no"' "Plus state for the electron fraction"<br>
<br>
_______________________________________________<br>
Commits mailing list<br>
<a href="mailto:Commits@einsteintoolkit.org">Commits@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/commits" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/commits</a><br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu" target="_blank">schnetter@cct.lsu.edu</a>><br><a href="http://www.perimeterinstitute.ca/personal/eschnetter/" target="_blank">http://www.perimeterinstitute.ca/personal/eschnetter/</a><br>
</div>