diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:35:16 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-05-29 16:35:16 +0000 |
| commit | a8d6c1d56668352c89da49bc82fc6a06d688052d (patch) | |
| tree | bf7dfe2cdaaec222067c7fb10e12a464236b1817 | |
| parent | f1d2591498f08a6912fa6cb81466683552173874 (diff) | |
GRHydro: fixed critical bug in temperature reconstruction in ePPM
From: Christian David Ott <cott@zwicky.(none)>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@532 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_PPM.F90 | 13 |
1 files changed, 7 insertions, 6 deletions
diff --git a/src/GRHydro_PPM.F90 b/src/GRHydro_PPM.F90 index b53ded9..2f65bcb 100644 --- a/src/GRHydro_PPM.F90 +++ b/src/GRHydro_PPM.F90 @@ -1158,13 +1158,14 @@ subroutine SimplePPM_temperature_1d(& tempminus(i+1) = tempplus(i) enddo - - !! Finally compute pressure gradient needed for flattening and shock detection - do i = 2, nx-1 - dpress(i) = press(i+1) - press(i-1) - end do - + endif + + !! Finally compute pressure gradient needed for flattening and shock detection + do i = 2, nx-1 + dpress(i) = press(i+1) - press(i-1) + end do + endif !!$ Zone flattening. See appendix of C&W, p. 197-8. |