User: rhaas Date: 2013/12/02 03:50 PM
Modified: /trunk/src/ GRHydro_ENOReconstruct_drv.F90, GRHydro_MP5Reconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv_opt.F90, GRHydro_ReconstructPoly.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_WENOReconstruct_drv.F90
Log: GRHydro: use ash instead of lsh when allocating arrays
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -61,7 +61,7 @@
CCTK_INT :: ierr
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if @@ -69,9 +69,9 @@
!!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr)
if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm")
File [modified]: GRHydro_MP5Reconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -65,15 +65,15 @@ CCTK_INT :: ierr
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
!!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr)
if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm")
File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -61,15 +61,15 @@
CCTK_INT :: ierr
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
!!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr)
call SpaceMask_GetTypeBits(type_bitsx, "Hydro_RiemannProblemX") call SpaceMask_GetStateBits(trivialx, "Hydro_RiemannProblemX", &
File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -57,7 +57,7 @@
logical :: apply_enhanced_ppm
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
File [modified]: GRHydro_PPMReconstruct_drv_opt.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMReconstruct_drv_opt.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_PPMReconstruct_drv_opt.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -90,7 +90,7 @@ pvup = loc(vel) end if
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: +8 -8 =================================================================== --- trunk/src/GRHydro_ReconstructPoly.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_ReconstructPoly.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -74,18 +74,18 @@
CCTK_INT :: GRHydro_UseGeneralCoordinates
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
- allocate(psi4(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetax(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetay(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - lbetaz(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(psi4(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetax(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetay(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + lbetaz(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr)
if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with lbeta")
File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -84,7 +84,7 @@ CCTK_REAL, DIMENSION(:,:,:,:), ALLOCATABLE :: wvel ! vel*w_lorentz CCTK_REAL :: agxx,agxy,agxz,agyy,agyz,agzz ! metric components
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: +4 -4 =================================================================== --- trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-12-02 21:50:25 UTC (rev 593) +++ trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-12-02 21:50:34 UTC (rev 594) @@ -125,15 +125,15 @@ ! CCTK_REAL, dimension(ny) :: vxY,vyY,vzY ! CCTK_REAL, dimension(nz) :: vxZ,vyZ,vzZ
- allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(trivial_rp(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") end if
!!$ The dum variables are used as dummies if MHD on but divergence cleaning off - allocate(dum(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dump(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),& - dumm(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) + allocate(dum(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dump(cctk_ash(1),cctk_ash(2),cctk_ash(3)),& + dumm(cctk_ash(1),cctk_ash(2),cctk_ash(3)),STAT=ierr)
if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with dum dump dumm")