User: reisswig Date: 2012/02/16 05:18 PM

Modified: /trunk/src/ GRHydro_Eigenproblem.F90, GRHydro_PPM.F90

Log: Cosmetic change to ePPM3.

Improved handling of cs2 < 0.

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2012-02-09 01:03:16 UTC (rev 313) +++ trunk/src/GRHydro_Eigenproblem.F90 2012-02-16 23:18:45 UTC (rev 314) @@ -87,7 +87,7 @@

if(cs2.lt.0.0d0) then if (abs(cs2) .gt. 1.0d-4) then - write(warnline,'(a50,6g16.7)') 'rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho + write(warnline,'(a50,6g16.7)') 'abs(cs2), rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho call CCTK_WARN(1,warnline) call CCTK_WARN(0,"cs2 < 0! Check speed of sound calculation!") else

File [modified]: GRHydro_PPM.F90 Delta lines: +5 -5 =================================================================== --- trunk/src/GRHydro_PPM.F90 2012-02-09 01:03:16 UTC (rev 313) +++ trunk/src/GRHydro_PPM.F90 2012-02-16 23:18:45 UTC (rev 314) @@ -667,16 +667,16 @@ aplus(i) = a(i) &&\ endif &&\ endif &&\ - endif &&\ + endif &&\ trivial_rp(i-1) = .false. &&\ trivial_rp(i) = .false. &&\ else &&\ if (abs(daplus) .le. abs(a(i)) .and. abs(daminus) .le. abs(a(i))) then &&\ if (abs(daplus) .ge. 2.0d0*abs(daminus)) then &&\ - aplus(i) = a(i) - 2.0d0*(aminus(i) - a(i)) &&\ + aplus(i) = a(i) + 2.0d0*daminus &&\ else if (abs(daminus) .ge. 2.0d0*abs(daplus)) then &&\ - aminus(i) = a(i) - 2.0d0*(aplus(i) - a(i)) &&\ - end if &&\ + aminus(i) = a(i) - 2.0d0*daplus &&\ + end if &&\ else &&\ if (abs(daplus) .ge. 2.0d0*abs(daminus)) then &&\ aplus(i) = a(i) &&\ @@ -686,7 +686,7 @@ endif &&\ trivial_rp(i-1) = .false. &&\ trivial_rp(i) = .false. &&\ - endif + endif