diff options
author | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-02-16 23:18:45 +0000 |
---|---|---|
committer | reisswig <reisswig@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-02-16 23:18:45 +0000 |
commit | 52d7f2ffa9c499f1b35b24771b889302bc61f50f (patch) | |
tree | e4c7950aced0242782671ecc1d81c1a796c87c13 | |
parent | 0a5f2247d543da932fe3a1459be8e44a369dc8e0 (diff) |
Cosmetic change to ePPM3.
Improved handling of cs2 < 0.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@314 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r-- | src/GRHydro_Eigenproblem.F90 | 2 | ||||
-rw-r--r-- | src/GRHydro_PPM.F90 | 10 |
2 files changed, 6 insertions, 6 deletions
diff --git a/src/GRHydro_Eigenproblem.F90 b/src/GRHydro_Eigenproblem.F90 index 605f372..1bbcf35 100644 --- a/src/GRHydro_Eigenproblem.F90 +++ b/src/GRHydro_Eigenproblem.F90 @@ -87,7 +87,7 @@ subroutine eigenvalues(handle,rho,velx,vely,velz,eps, & 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 diff --git a/src/GRHydro_PPM.F90 b/src/GRHydro_PPM.F90 index aedeff7..832cd38 100644 --- a/src/GRHydro_PPM.F90 +++ b/src/GRHydro_PPM.F90 @@ -667,16 +667,16 @@ trivial_rp = .true. 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 @@ trivial_rp = .true. endif &&\ trivial_rp(i-1) = .false. &&\ trivial_rp(i) = .false. &&\ - endif + endif |