[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 532)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/05/29 11:35 AM

Modified:
 /trunk/src/
  GRHydro_PPM.F90

Log:
 GRHydro: fixed critical bug in temperature reconstruction in ePPM
 
 From: Christian David Ott <cott at zwicky.(none)>

File Changes:

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

File [modified]: GRHydro_PPM.F90
Delta lines: +7 -6
===================================================================
--- trunk/src/GRHydro_PPM.F90	2013-05-29 16:35:12 UTC (rev 531)
+++ trunk/src/GRHydro_PPM.F90	2013-05-29 16:35:16 UTC (rev 532)
@@ -1158,13 +1158,14 @@
             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.