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@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.