diff options
Diffstat (limited to 'src/GRHydro_PPM.F90')
| -rw-r--r-- | src/GRHydro_PPM.F90 | 20 |
1 files changed, 2 insertions, 18 deletions
diff --git a/src/GRHydro_PPM.F90 b/src/GRHydro_PPM.F90 index f71a0e1..892d453 100644 --- a/src/GRHydro_PPM.F90 +++ b/src/GRHydro_PPM.F90 @@ -73,13 +73,6 @@ subroutine SimplePPM_1d(handle,poly,& CCTK_REAL, dimension(nx) :: rho,velx,vely,velz,eps CCTK_REAL, dimension(nx) :: rhominus,velxminus,velyminus,velzminus,epsminus CCTK_REAL, dimension(nx) :: rhoplus,velxplus,velyplus,velzplus,epsplus - CCTK_REAL, dimension(nx) :: rhominusl,velxminusl,velyminusl,velzminusl - CCTK_REAL, dimension(nx) :: epsminusl - CCTK_REAL, dimension(nx) :: rhoplusl,velxplusl,velyplusl,velzplusl,epsplusl - CCTK_REAL, dimension(nx) :: rhominusr,velxminusr,velyminusr,velzminusr - CCTK_REAL, dimension(nx) :: epsminusr - CCTK_REAL, dimension(nx) :: rhoplusr,velxplusr,velyplusr,velzplusr,epsplusr - CCTK_REAL, dimension(nx) :: atmosphere_mask CCTK_INT :: i,s CCTK_REAL, dimension(nx) :: drho,dvelx,dvely,dvelz,deps @@ -98,7 +91,6 @@ subroutine SimplePPM_1d(handle,poly,& CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, det CCTK_REAL, dimension(5) :: lam - CCTK_REAL :: dupw, dloc, delta CCTK_REAL :: agxx, agxy, agxz, agyy, agyz, agzz CCTK_REAL, dimension(nx) :: xwind, l_ev_l, l_ev_r @@ -1046,17 +1038,12 @@ subroutine SimplePPM_temperature_1d(& CCTK_REAL, dimension(nx) :: temperature, velx CCTK_REAL, dimension(nx) :: tempminus CCTK_REAL, dimension(nx) :: tempplus - CCTK_REAL, dimension(nx) :: tempminusl - CCTK_REAL, dimension(nx) :: tempplusl - CCTK_REAL, dimension(nx) :: tempminusr - CCTK_REAL, dimension(nx) :: tempplusr - CCTK_REAL, dimension(nx) :: atmosphere_mask CCTK_INT :: i,s CCTK_REAL, dimension(nx) :: dtemp CCTK_REAL, dimension(nx) :: dmtemp CCTK_REAL, dimension(nx) :: press,dpress,tilde_flatten - CCTK_REAL :: dpress2,dvel,w,flatten,eta,etatilde + CCTK_REAL :: dpress2,dvel,w,flatten CCTK_INT, dimension(nx) :: hydro_excision_mask @@ -1569,7 +1556,6 @@ subroutine SimplePPM_tracer_1d(nx,dx,rho,velx,vely,velz, & CCTK_REAL, dimension(nx) :: press,dpress,tilde_flatten CCTK_REAL, dimension(nx,number_of_tracers) :: dmtracer, dtracer, tracerflat!, d2tracer CCTK_REAL :: dpress2,w,flatten,dvel - CCTK_REAL :: eta, etatilde !!$ Average slopes delta_m(a). See (1.7) of Colella and Woodward, p.178 !!$ This is the expression for an even grid. @@ -1777,14 +1763,12 @@ subroutine SimplePPM_ye_1d(nx,dx,rho,velx,vely,velz, & CCTK_REAL :: dx CCTK_REAL, dimension(nx) :: rho,velx,vely,velz CCTK_REAL, dimension(nx) :: Y_e,Y_e_minus,Y_e_plus - CCTK_REAL :: Y_eflatomega CCTK_INT :: i,s CCTK_REAL, dimension(nx) :: press,dpress,tilde_flatten - CCTK_REAL, dimension(nx) :: dmY_e, dY_e, Y_eflat!, d2tracer + CCTK_REAL, dimension(nx) :: dmY_e, dY_e CCTK_REAL :: dpress2,w,flatten,dvel - CCTK_REAL :: eta, etatilde CCTK_REAL :: D2a, D2aL, D2aR, D2aC, D2aLim, rhi, daplus, daminus, D3a, D3aLL, D3aL, D3aR, D2aLL, D2aRR, D3aMin, D3aMax |