From a8d6c1d56668352c89da49bc82fc6a06d688052d Mon Sep 17 00:00:00 2001
From: rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>
Date: Wed, 29 May 2013 16:35:16 +0000
Subject: 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
---
 src/GRHydro_PPM.F90 | 13 +++++++------
 1 file 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.
-- 
cgit v1.2.3