[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 495)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Wed Mar 27 20:46:39 CDT 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 at 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
More information about the Commits
mailing list