diff options
Diffstat (limited to 'src/GRHydro_PPM.F90')
-rw-r--r-- | src/GRHydro_PPM.F90 | 18 |
1 files changed, 0 insertions, 18 deletions
diff --git a/src/GRHydro_PPM.F90 b/src/GRHydro_PPM.F90 index 9a8a3ce..cc42787 100644 --- a/src/GRHydro_PPM.F90 +++ b/src/GRHydro_PPM.F90 @@ -64,11 +64,6 @@ subroutine SimplePPM_1d(handle,poly,& DECLARE_CCTK_PARAMETERS -#ifdef _EOS_BASE_INC_ -#undef _EOS_BASE_INC_ -#endif -#include "EOS_Base.inc" - CCTK_INT :: handle,poly,nx CCTK_REAL :: dx CCTK_REAL, dimension(nx) :: rho,velx,vely,velz,eps @@ -109,19 +104,6 @@ subroutine SimplePPM_1d(handle,poly,& trivial_rp = .true. -!!$ Loop to set the pressure and eps -!!$ poly=0 -!!$ if (poly .eq. 0) then -!!$ do i = 1, nx -!!$ press(i) = EOS_Pressure(handle, rho(i), eps(i)) -!!$ end do -!!$ else -!!$ do i = 1, nx -!!$ press(i) = EOS_Pressure(handle, rho(i), eps(i)) -!!$ eps(i) = EOS_SpecificIntEnergy(handle, rho(i), press(i)) -!!$ end do -!!$ end if - !!$ Average slopes delta_m(a). See (1.7) of Colella and Woodward, p.178 !!$ This is the expression for an even grid. |