diff options
-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. |