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

rhaas at tapir.caltech.edu rhaas at tapir.caltech.edu
Mon Dec 2 15:50:34 CST 2013


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")



More information about the Commits mailing list