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

27 Mar 2013

User: rhaas
Date: 2013/03/27 08:46 PM
Modified:
 /trunk/src/
  GRHydro_Boundaries.F90, GRHydro_Con2Prim.F90, GRHydro_Differences.F90, GRHydro_ENOReconstruct.F90, GRHydro_Eigenproblem.F90, GRHydro_EigenproblemM.F90, GRHydro_Eigenproblem_Marquina.F90, GRHydro_EoSChangeGamma.F90, GRHydro_Flux.F90, GRHydro_FluxSplit.F90, GRHydro_HLLE.F90, GRHydro_HLLEM.F90, GRHydro_HLLE_AM.F90, GRHydro_MP5Reconstruct.F90, GRHydro_Marquina.F90, GRHydro_PPM.F90, GRHydro_PPMM.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_Particle.F90, GRHydro_Prim2Con.F90, GRHydro_Prim2ConAM.F90, GRHydro_Prim2ConM.F90, GRHydro_Reconstruct.F90, GRHydro_RiemannSolve.F90, GRHydro_RiemannSolveAM.F90, GRHydro_RiemannSolveM.F90, GRHydro_RoeSolver.F90, GRHydro_WENOReconstruct.F90
Log:
 GRHydro: remove unused variables
From: Roland Haas rhaas@tapir.caltech.edu
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_Boundaries.F90
Delta lines: +0 -3
===================================================================

--- trunk/src/GRHydro_Boundaries.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Boundaries.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -376,7 +376,6 @@
integer, dimension(3) :: sw
   integer :: ierr = 0
-  integer :: i,j,k
   CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates
CCTK_INT, parameter :: faces=CCTK_ALL_FACES
@@ -515,7 +514,6 @@
integer, dimension(3) :: sw
   integer :: ierr = 0
-  integer :: i,j,k
   CCTK_INT :: GRHydro_UseGeneralCoordinates, general_coordinates
CCTK_INT, parameter :: faces=CCTK_ALL_FACES
@@ -651,7 +649,6 @@
integer, dimension(3) :: sw
   integer :: ierr = 0
-  integer :: i,j,k
CCTK_INT, parameter :: faces=CCTK_ALL_FACES
   CCTK_INT, parameter :: ione=1
File [modified]: GRHydro_Con2Prim.F90
Delta lines: +8 -9
===================================================================
--- trunk/src/GRHydro_Con2Prim.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Con2Prim.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -54,7 +54,6 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
-  integer :: CCTK_MyProc
   integer :: i, j, k, itracer, nx, ny, nz
   CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, det, pmin, epsmin, dummy1, dummy2
   logical :: epsnegative
@@ -416,11 +415,11 @@
        press, uxx, uxy, uxz, uyy, uyz, uzz, det, w_lorentz, x, &
        y, z, r, GRHydro_rho_min
-  CCTK_REAL s2, c0, c1, c2, c3, c4, f, df, ftol, v2, w, vlowx, vlowy, vlowz
+  CCTK_REAL s2, f, df, vlowx, vlowy, vlowz
   CCTK_INT count, i, handle, GRHydro_reflevel
   CCTK_REAL GRHydro_C2P_failed, dummy1, dummy2
-  CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, epsold, epsnew, w2, &
-       w2mhalf, temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin
+  CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, &
+            temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin
   character(len=200) warnline
   logical epsnegative, mustbisect
@@ -778,11 +777,11 @@
        press, uxx, uxy, uxz, uyy, uyz, uzz, det, w_lorentz, x, &
        y, z, r, GRHydro_rho_min
   CCTK_REAL temp, ye
-  CCTK_REAL s2, c0, c1, c2, c3, c4, f, df, ftol, v2, w, vlowx, vlowy, vlowz
+  CCTK_REAL s2, f, df, vlowx, vlowy, vlowz
   CCTK_INT cctk_iteration, ii,jj,kk,count, i, handle, GRHydro_reflevel
   CCTK_REAL GRHydro_C2P_failed
-  CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, epsold, epsnew, w2, &
-       w2mhalf, temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin
+  CCTK_REAL udens, usx, usy, usz, utau, pold, pnew, &
+            temp1, drhobydpress, depsbydpress, dpressbydeps, dpressbydrho, pmin, epsmin
   CCTK_REAL epsminl,pminl,plow,tmp, dummy1, dummy2
   CCTK_REAL local_perc_ptol
@@ -1248,7 +1247,7 @@
     w_lorentz = 1.d0
   end if
-51 press = pnew
+  press = pnew
   vlowx = usx / ( (rho + rho*epsilon + press) * w_lorentz**2)
   vlowy = usy / ( (rho + rho*epsilon + press) * w_lorentz**2)
   vlowz = usz / ( (rho + rho*epsilon + press) * w_lorentz**2)
@@ -1666,7 +1665,7 @@
CCTK_REAL s2, f, df, vlowx, vlowy, vlowz
   CCTK_INT count, handle, GRHydro_reflevel
-  CCTK_REAL udens, usx, usy, usz, rhoold, rhonew, epsold, epsnew, &
+  CCTK_REAL udens, usx, usy, usz, rhoold, rhonew, &
        enthalpy, denthalpy, sqrtdet, invsqrtdet, invfac, GRHydro_C2P_failed, dummy1, dummy2
   character(len=200) warnline
File [modified]: GRHydro_Differences.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_Differences.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Differences.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -27,7 +27,7 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i, j, k, nx, ny, nz, kk, ll
+  CCTK_INT :: i, j, k, nx, ny, nz
   CCTK_REAL :: idx, idy, idz
CCTK_REAL :: dr_eps = 1.d-2
File [modified]: GRHydro_ENOReconstruct.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_ENOReconstruct.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_ENOReconstruct.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -40,7 +40,7 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
-  CCTK_INT :: i, j, k, l, m, q, allocstat
+  CCTK_INT :: i, j, l, m, q, allocstat
   CCTK_REAL :: denominator, numerator, numerator_product
if(.not.coeffs_allocated) then
File [modified]: GRHydro_Eigenproblem.F90
Delta lines: +1 -2
===================================================================
--- trunk/src/GRHydro_Eigenproblem.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Eigenproblem.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -145,7 +145,6 @@
   CCTK_REAL vlowx,vlowy,vlowz,v2,w
   CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta
   CCTK_INT handle,ii,jj,kk
-  CCTK_REAL dpdrho,dpdeps,press
! begin EOS Omni vars
   integer :: n,keytemp,anyerr,keyerr(1)
@@ -598,7 +597,7 @@
   CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta
   CCTK_REAL cs2,one,two
   CCTK_REAL vlowx,vlowy,vlowz,v2,w
-  CCTK_REAL press,dpdrho,dpdeps,enthalpy,kappa
+  CCTK_REAL press,dpdeps,enthalpy,kappa
   CCTK_REAL axp,axm,vxp,vxm,cxx,cxy,cxz,gam,xsi,dlt
   CCTK_INT handle
File [modified]: GRHydro_EigenproblemM.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_EigenproblemM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_EigenproblemM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -155,7 +155,7 @@
   CCTK_REAL vlowx,vlowy,vlowz,v2,w
   CCTK_REAL lam1,lam2,lam3,lamm,lamp,lamm_nobeta,lamp_nobeta
   CCTK_INT handle, ii,jj,kk
-  CCTK_REAL dpdrho,dpdeps,press
+  CCTK_REAL press
CCTK_REAL Bvcxlow,Bvcylow,Bvczlow,Bvc2,rhohstar,va2
   CCTK_REAL Bdotv,b2
File [modified]: GRHydro_Eigenproblem_Marquina.F90
Delta lines: +10 -14
===================================================================
--- trunk/src/GRHydro_Eigenproblem_Marquina.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Eigenproblem_Marquina.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -49,8 +49,8 @@
DECLARE_CCTK_PARAMETERS
-  CCTK_REAL rhor,velxr,velyr,velzr,epsr,w_lorentzr
-  CCTK_REAL rhol,velxl,velyl,velzl,epsl,w_lorentzl
+  CCTK_REAL rhor,velxr,velyr,velzr,epsr
+  CCTK_REAL rhol,velxl,velyl,velzl,epsl
   CCTK_REAL densl,sxl,syl,szl,taul
   CCTK_REAL densr,sxr,syr,szr,taur
@@ -63,8 +63,8 @@
!!$  LOCAL VARS
-  CCTK_REAL du(5),aa(5,5),qdiff1,qdiff2,qdiff3,qdiff4,qdiff5
-  integer i,j,k,l
+  CCTK_REAL du(5),aa(5,5)
+  integer i,j,k
   CCTK_REAL paug(5,6),tmp1,tmp2,sump,summ,f_du(5)
   CCTK_REAL leivec1l(5),leivec2l(5),leivec3l(5),leivecpl(5),leivecml(5)
   CCTK_REAL leivec1r(5),leivec2r(5),leivec3r(5),leivecpr(5),leivecmr(5)
@@ -78,8 +78,6 @@
   CCTK_REAL vlowxr,vlowyr,vlowzr,v2r,wr
   CCTK_REAL vlowxl,vlowyl,vlowzl,v2l,wl
-  CCTK_REAL lamp_nobeta,lamm_nobeta
-
   CCTK_REAL pressl,dpdrhol,dpdepsl,enthalpyl,kappal
   CCTK_REAL pressr,dpdrhor,dpdepsr,enthalpyr,kappar
   CCTK_REAL axpl,axml,vxpl,vxml,xsil,dltl
@@ -774,8 +772,8 @@
DECLARE_CCTK_PARAMETERS
-  CCTK_REAL rhor,velxr,velyr,velzr,epsr,w_lorentzr
-  CCTK_REAL rhol,velxl,velyl,velzl,epsl,w_lorentzl
+  CCTK_REAL rhor,velxr,velyr,velzr,epsr
+  CCTK_REAL rhol,velxl,velyl,velzl,epsl
   CCTK_REAL densl,sxl,syl,szl,taul
   CCTK_REAL densr,sxr,syr,szr,taur
   CCTK_REAL templ,tempr,ye_l,ye_r
@@ -789,8 +787,8 @@
!!$  LOCAL VARS
-  CCTK_REAL du(5),aa(5,5),qdiff1,qdiff2,qdiff3,qdiff4,qdiff5
-  integer i,j,k,l
+  CCTK_REAL du(5),aa(5,5)
+  integer i,j,k
   CCTK_REAL paug(5,6),tmp1,tmp2,sump,summ,f_du(5)
   CCTK_REAL leivec1l(5),leivec2l(5),leivec3l(5),leivecpl(5),leivecml(5)
   CCTK_REAL leivec1r(5),leivec2r(5),leivec3r(5),leivecpr(5),leivecmr(5)
@@ -804,10 +802,8 @@
   CCTK_REAL vlowxr,vlowyr,vlowzr,v2r,wr
   CCTK_REAL vlowxl,vlowyl,vlowzl,v2l,wl
-  CCTK_REAL lamp_nobeta,lamm_nobeta
-
-  CCTK_REAL pressl,dpdrhol,dpdepsl,enthalpyl,kappal
-  CCTK_REAL pressr,dpdrhor,dpdepsr,enthalpyr,kappar
+  CCTK_REAL pressl,dpdepsl,enthalpyl,kappal
+  CCTK_REAL pressr,dpdepsr,enthalpyr,kappar
   CCTK_REAL axpl,axml,vxpl,vxml,xsil,dltl
   CCTK_REAL axpr,axmr,vxpr,vxmr,xsir,dltr
   CCTK_REAL cxx,cxy,cxz,gam,vxa,vxb
File [modified]: GRHydro_EoSChangeGamma.F90
Delta lines: +0 -4
===================================================================
--- trunk/src/GRHydro_EoSChangeGamma.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_EoSChangeGamma.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -291,10 +291,6 @@
CCTK_REAL :: local_Gamma, local_k, eos_k_initial
-  CCTK_REAL, dimension(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)) :: Q
-
-  character(len=100) infoline
-
   ! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
   CCTK_REAL, DIMENSION(:,:,:), POINTER :: g11, g12, g13, g22, g23, g33
File [modified]: GRHydro_Flux.F90
Delta lines: +0 -1
===================================================================
--- trunk/src/GRHydro_Flux.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Flux.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -36,7 +36,6 @@
   CCTK_REAL :: dens,sx,sy,sz,tau,velx,vely,velz,det
   CCTK_REAL :: densf,sxf,syf,szf,tauf
   CCTK_REAL :: alp,beta,press 
-  integer :: i,j,k
   integer :: xoffset,yoffset,zoffset
densf = dens * ( (velx*xoffset + vely*yoffset + &
File [modified]: GRHydro_FluxSplit.F90
Delta lines: +2 -2
===================================================================
--- trunk/src/GRHydro_FluxSplit.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_FluxSplit.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -40,7 +40,7 @@
   DECLARE_CCTK_FUNCTIONS
integer :: nx, ny, nz, i, j, k, ierr, max_handle
-  CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, det, psi4pt, beta
+  CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, det, beta
   CCTK_REAL, dimension(5) :: lambda
   CCTK_REAL :: alpha1_local, alpha2_local, alpha3_local, alpha4_local, &
        alpha5_local
@@ -287,7 +287,7 @@
   integer :: nx, ny, nz, i, j, k
CCTK_REAL, dimension(:), allocatable :: upper, det, dummy
-  CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz, psi4pt
+  CCTK_REAL :: uxx, uxy, uxz, uyy, uyz, uzz
nx = cctk_lsh(1)
   ny = cctk_lsh(2)
File [modified]: GRHydro_HLLE.F90
Delta lines: +4 -4
===================================================================
--- trunk/src/GRHydro_HLLE.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_HLLE.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -166,9 +166,9 @@
   CCTK_REAL, dimension(5) :: qdiff
   CCTK_REAL ::  charmin, charmax, charpm,avg_alp,avg_det, etabar
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
! sign requires its arguments to be of identical KIND
   CCTK_REAL, parameter :: one = 1d0
@@ -634,9 +634,9 @@
   CCTK_REAL, dimension(number_of_tracers) :: qdiff
   CCTK_REAL ::  charmin, charmax, charpm,avg_alp,avg_det
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
File [modified]: GRHydro_HLLEM.F90
Delta lines: +4 -6
===================================================================
--- trunk/src/GRHydro_HLLEM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_HLLEM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -53,8 +53,7 @@
   CCTK_REAL, dimension(5) :: lamminus,lamplus
   CCTK_REAL ::  charmin, charmax, charpm,chartop,avg_alp,avg_det, sdet
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r, &
-       cs2_p, cs2_m, dpdeps_p, dpdeps_m
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
   CCTK_REAL :: rhoenth_p, rhoenth_m, avg_betax, avg_betay, avg_betaz
   CCTK_REAL :: vxtp,vytp,vztp,vxtm,vytm,vztm,ab0p,ab0m,b2p,b2m,bdotvp,bdotvm
   CCTK_REAL :: wp,wm,v2p,v2m,bxlowp,bxlowm,bylowp,bylowm,bzlowp,bzlowm,vA2m,vA2p
@@ -66,7 +65,7 @@
   CCTK_REAL :: psidcp, psidcm, psidcf, psidcdiff, psidcfp, psidcfm
   CCTK_REAL :: charmax_dc, charmin_dc, charpm_dc
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
   CCTK_REAL :: xtemp
! save memory when MP is not used
@@ -724,11 +723,10 @@
   CCTK_REAL, dimension(3) :: mag_p, mag_m
   CCTK_REAL ::  charmin, charmax, charpm,avg_alp,avg_det
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r, &
-       cs2_p, cs2_m, dpdeps_p, dpdeps_m
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
   CCTK_REAL :: b2p,b2m,vA2m,vA2p
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
File [modified]: GRHydro_HLLE_AM.F90
Delta lines: +4 -6
===================================================================
--- trunk/src/GRHydro_HLLE_AM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_HLLE_AM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -53,8 +53,7 @@
   CCTK_REAL, dimension(5) :: lamminus,lamplus
   CCTK_REAL ::  charmin, charmax, charpm,chartop,avg_alp,avg_det, sdet
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r, &
-       cs2_p, cs2_m, dpdeps_p, dpdeps_m
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
   CCTK_REAL :: rhoenth_p, rhoenth_m, avg_betax, avg_betay, avg_betaz
   CCTK_REAL :: vxtp,vytp,vztp,vxtm,vytm,vztm,ab0p,ab0m,b2p,b2m,bdotvp,bdotvm
   CCTK_REAL :: wp,wm,v2p,v2m,bxlowp,bxlowm,bylowp,bylowm,bzlowp,bzlowm,vA2m,vA2p
@@ -64,7 +63,7 @@
   CCTK_REAL :: psidcp, psidcm, psidcf, psidcdiff, psidcfp, psidcfm
   CCTK_REAL :: charmax_dc, charmin_dc, charpm_dc
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
   CCTK_REAL :: xtemp
! save memory when MP is not used
@@ -681,11 +680,10 @@
   CCTK_REAL, dimension(3) :: mag_p, mag_m
   CCTK_REAL ::  charmin, charmax, charpm,avg_alp,avg_det
   CCTK_REAL :: gxxh, gxyh, gxzh, gyyh, gyzh, gzzh, uxxh, uxyh, &
-       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh, alp_l, alp_r, &
-       cs2_p, cs2_m, dpdeps_p, dpdeps_m
+       uxzh, uyyh, uyzh, uzzh, avg_beta, usendh
   CCTK_REAL :: b2p,b2m,vA2m,vA2p
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
File [modified]: GRHydro_MP5Reconstruct.F90
Delta lines: +2 -2
===================================================================
--- trunk/src/GRHydro_MP5Reconstruct.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_MP5Reconstruct.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -40,7 +40,7 @@
DECLARE_CCTK_PARAMETERS
-  CCTK_INT :: nx, i, j, k, r
+  CCTK_INT :: nx, i
   CCTK_REAL, dimension(nx) :: v, vplus, vminus
CCTK_INT, dimension(nx) :: hydro_excision_mask
@@ -95,7 +95,7 @@
 #define MINMOD4(w,x,y,z) \
      0.125d0*( sign(one,w)+sign(one,x) )*abs( (sign(one,w)+sign(one,y)) * (sign(one,w)+sign(one,z)) )*min(abs(w), abs(x), abs(y), abs(z))
-#define MP5(am2, am1, a, ap1, ap2, arecon) \    
+#define MP5(am2, am1, a, ap1, ap2, arecon) \
       vl = (2.0d0*am2 - 13.0d0*am1 + 47.0d0*a + 27.0d0*ap1 - 3.0d0*ap2)/60.0d0  &&\
       vmp = a + MINMOD( ap1-a, mp5_alpha*(a-am1) )  &&\
       if ((vl-a)*(vl-vmp) .le. mp5_eps) then &&\
File [modified]: GRHydro_Marquina.F90
Delta lines: +1 -2
===================================================================
--- trunk/src/GRHydro_Marquina.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Marquina.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -45,13 +45,12 @@
          consp,consm_i,fplus,fminus,f_marquina,primp,primm_i
     CCTK_REAL :: avg_alp,avg_beta,gxxh,gxyh,gxzh,gyyh,gyzh,gzzh, &
          avg_det,uxxh,uxyh,uxzh,uyyh,uyzh,uzzh,&
-         pressp,pressm_i, &
          tmp_w_lorentzp, tmp_w_lorentzm_i, w_lorentzp,w_lorentzm_i, usendh, psi4h
     integer :: m
     integer :: i,j,k
     integer :: keytemp
-    CCTK_INT :: type_bits, trivial, not_trivial
+    CCTK_INT :: type_bits, trivial
if(evolve_temper.eq.1.and.reconstruct_temper.eq.1) then
        keytemp = 1
File [modified]: GRHydro_PPM.F90
Delta lines: +2 -18
===================================================================
--- trunk/src/GRHydro_PPM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_PPM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -73,13 +73,6 @@
   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 @@
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 @@
   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 @@
   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 @@
   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
File [modified]: GRHydro_PPMM.F90
Delta lines: +0 -1
===================================================================
--- trunk/src/GRHydro_PPMM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_PPMM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -73,7 +73,6 @@
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
File [modified]: GRHydro_PPMMReconstruct_drv.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_PPMMReconstruct_drv.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_PPMMReconstruct_drv.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -60,7 +60,7 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  integer :: nx, ny, nz, i, j, k, itracer
+  integer :: nx, ny, nz, i, j, k
logical, dimension(:,:,:), allocatable :: trivial_rp
File [modified]: GRHydro_PPMReconstruct_drv.F90
Delta lines: +1 -1
===================================================================
--- trunk/src/GRHydro_PPMReconstruct_drv.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_PPMReconstruct_drv.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -53,7 +53,7 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  integer :: nx, ny, nz, i, j, k, itracer
+  integer :: nx, ny, nz, i, j, k
logical, dimension(:,:,:), allocatable :: trivial_rp
File [modified]: GRHydro_Particle.F90
Delta lines: +1 -2
===================================================================
--- trunk/src/GRHydro_Particle.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Particle.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -36,14 +36,13 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i, l
+  CCTK_INT :: i
   CCTK_INT :: interp_handle, table_handle, status, coord_system_handle
character(len=200) :: particle_interp
   character(len=128) :: warn_message
   CCTK_INT :: particle_interp_len
   character(len=7) :: particle_order
-  character(len=15) :: vname
CCTK_INT, dimension(1) :: lsh
   CCTK_POINTER, dimension(3) :: interp_coords
File [modified]: GRHydro_Prim2Con.F90
Delta lines: +0 -1
===================================================================
--- trunk/src/GRHydro_Prim2Con.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Prim2Con.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -441,7 +441,6 @@
   CCTK_INT :: i, j, k
   CCTK_INT :: keytemp
   CCTK_REAL :: det
-  character(len=512) :: warnline
! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
File [modified]: GRHydro_Prim2ConAM.F90
Delta lines: +1 -2
===================================================================
--- trunk/src/GRHydro_Prim2ConAM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Prim2ConAM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -438,7 +438,6 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
-  CCTK_REAL :: xtemp(1)
   CCTK_INT :: i, j, k
   CCTK_REAL :: det
   CCTK_REAL :: maxtau0
@@ -631,7 +630,7 @@
   CCTK_REAL :: gxx, gxy, gxz, gyy, gyz, gzz, det
   CCTK_REAL :: ddens, dsx, dsy, dsz, dtau, &
        drho, dvelx, dvely, dvelz, deps, dpress, dBvcx, dBvcy, dBvcz, &
-       w_tmp, w, vlowx, vlowy, vlowz, sqrtdet   
+       w_tmp, w, vlowx, vlowy, vlowz
   CCTK_INT :: handle
   CCTK_REAL :: Bdotv,ab0,b2,blowx,blowy,blowz
   character(len=256) NaN_WarnLine
File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +1 -2
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Prim2ConM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -455,7 +455,6 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
-  CCTK_REAL :: xtemp(1)
   CCTK_INT :: i, j, k
   CCTK_REAL :: det
   CCTK_REAL :: maxtau0
@@ -654,7 +653,7 @@
   CCTK_REAL :: gxx, gxy, gxz, gyy, gyz, gzz, det
   CCTK_REAL :: ddens, dsx, dsy, dsz, dtau, dBconsx, dBconsy, dBconsz, &
        drho, dvelx, dvely, dvelz, deps, dpress, dBvcx, dBvcy, dBvcz, &
-       w_tmp, w, vlowx, vlowy, vlowz, sqrtdet   
+       w_tmp, w, vlowx, vlowy, vlowz
   CCTK_INT :: handle
   CCTK_REAL :: Bdotv,ab0,b2,blowx,blowy,blowz
   character(len=256) NaN_WarnLine
File [modified]: GRHydro_Reconstruct.F90
Delta lines: +0 -1
===================================================================
--- trunk/src/GRHydro_Reconstruct.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_Reconstruct.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -49,7 +49,6 @@
   CCTK_REAL :: local_min_tracer, dummy1, dummy2
   CCTK_INT :: type_bits, not_trivial
   CCTK_REAL :: agxx, agxy, agxz, agyy, agyz, agzz, w
-  character*256 :: warnline
! save memory when MP is not used
   CCTK_INT :: GRHydro_UseGeneralCoordinates
File [modified]: GRHydro_RiemannSolve.F90
Delta lines: +0 -4
===================================================================
--- trunk/src/GRHydro_RiemannSolve.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_RiemannSolve.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -36,8 +36,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE")) then
call GRHydro_HLLE(CCTK_PASS_FTOF)
@@ -93,8 +91,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE")) then
call GRHydro_HLLE(CCTK_PASS_FTOF)
File [modified]: GRHydro_RiemannSolveAM.F90
Delta lines: +0 -4
===================================================================
--- trunk/src/GRHydro_RiemannSolveAM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_RiemannSolveAM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -36,8 +36,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE").or.CCTK_EQUALS(riemann_solver,"LLF")) then
call GRHydro_HLLE_AM(CCTK_PASS_FTOF)
@@ -100,8 +98,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE")) then
call GRHydro_HLLE_AM(CCTK_PASS_FTOF)
File [modified]: GRHydro_RiemannSolveM.F90
Delta lines: +0 -4
===================================================================
--- trunk/src/GRHydro_RiemannSolveM.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_RiemannSolveM.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -36,8 +36,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE").or.CCTK_EQUALS(riemann_solver,"LLF")) then
call GRHydro_HLLEM(CCTK_PASS_FTOF)
@@ -100,8 +98,6 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
-  CCTK_INT :: i,j,k
-  
   if (CCTK_EQUALS(riemann_solver,"HLLE")) then
call GRHydro_HLLEM(CCTK_PASS_FTOF)
File [modified]: GRHydro_RoeSolver.F90
Delta lines: +4 -4
===================================================================
--- trunk/src/GRHydro_RoeSolver.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_RoeSolver.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -38,15 +38,15 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
   CCTK_REAL, dimension(5) :: roeflux,roeave,qdiff,consp,consm_i,&
-       fplus,fminus,f_roe,primp,primm_i,consh
+       fplus,fminus,f_roe,primp,primm_i
   CCTK_REAL :: avg_alp,avg_beta,gxxh,gxyh,gxzh,gyyh,gyzh,gzzh, &
        avg_det,uxxh,uxyh,uxzh,uyyh,uyzh,uzzh,&
-       rhoave, velxave, velyave, velzave, pressave, epsave, &
-       w_lorentzave, usendh, alp_l, alp_r, psi4h
+       rhoave, velxave, velyave, velzave, epsave, &
+       w_lorentzave, usendh
   integer :: m
   integer :: i,j,k
-  CCTK_INT :: type_bits, trivial, not_trivial
+  CCTK_INT :: type_bits, trivial
if (flux_direction == 1) then
     call SpaceMask_GetTypeBits(type_bits, "Hydro_RiemannProblemX")
File [modified]: GRHydro_WENOReconstruct.F90
Delta lines: +3 -4
===================================================================
--- trunk/src/GRHydro_WENOReconstruct.F90	2013-03-28 01:46:34 UTC (rev 494)
+++ trunk/src/GRHydro_WENOReconstruct.F90	2013-03-28 01:46:38 UTC (rev 495)
@@ -40,8 +40,7 @@
   DECLARE_CCTK_ARGUMENTS
   DECLARE_CCTK_PARAMETERS
-  CCTK_INT :: i, j, k, l, m, q, allocstat
-  CCTK_REAL :: denominator, numerator, numerator_product
+  CCTK_INT :: allocstat
if(.not.coeffs_allocated) then
      ! Right now we hardcode to 5th order
@@ -156,7 +155,7 @@
DECLARE_CCTK_PARAMETERS
-  CCTK_INT :: order, nx, i, j, k, r
+  CCTK_INT :: order, nx, i, j
   CCTK_REAL, dimension(nx) :: v, vplus, vminus
CCTK_INT, dimension(nx) :: hydro_excision_mask
@@ -164,7 +163,7 @@
   logical, dimension(nx) :: excise
   logical :: normal_weno
-  CCTK_REAL :: large = 1.d10, gamma1, gamma2, gamma3, beta1, beta2, beta3, vnorm, betanorm
+  CCTK_REAL :: beta1, beta2, beta3, vnorm, betanorm
   CCTK_REAL :: wplus1, wplus2, wplus3, wbarplus1, wbarplus2, wbarplus3
   CCTK_REAL :: wminus1, wminus2, wminus3, wbarminus1, wbarminus2, wbarminus3

    