diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-19 23:15:08 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-19 23:15:08 +0000 |
| commit | bab01622dcacbc47d2c4a5a89549374bb863c9bd (patch) | |
| tree | 0552e722be64101eabbb4d70b71a65264d61df16 /src/GRHydro_ENOReconstruct_drv.F90 | |
| parent | 17f782f2a241db7f2ed8795cb4a1af7ed7ec28e5 (diff) | |
fix non-multipatch velocities in ENO reconstruction
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@275 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_ENOReconstruct_drv.F90')
| -rw-r--r-- | src/GRHydro_ENOReconstruct_drv.F90 | 17 |
1 files changed, 14 insertions, 3 deletions
diff --git a/src/GRHydro_ENOReconstruct_drv.F90 b/src/GRHydro_ENOReconstruct_drv.F90 index 247b398..1c2a28a 100644 --- a/src/GRHydro_ENOReconstruct_drv.F90 +++ b/src/GRHydro_ENOReconstruct_drv.F90 @@ -14,9 +14,9 @@ #include "SpaceMask.h" -#define velx(i,j,k) lvel(i,j,k,1) -#define vely(i,j,k) lvel(i,j,k,2) -#define velz(i,j,k) lvel(i,j,k,3) +#define velx(i,j,k) vup(i,j,k,1) +#define vely(i,j,k) vup(i,j,k,2) +#define velz(i,j,k) vup(i,j,k,3) #define sx(i,j,k) scon(i,j,k,1) #define sy(i,j,k) scon(i,j,k,2) #define sz(i,j,k) scon(i,j,k,3) @@ -64,6 +64,17 @@ subroutine GRHydro_ENOReconstruct_drv(CCTK_ARGUMENTS) CCTK_INT :: ierr + ! save memory when MP is not used + CCTK_INT :: GRHydro_UseGeneralCoordinates + TARGET vel, lvel + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + + if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then + vup => lvel + else + vup => vel + end if + allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) if (ierr .ne. 0) then call CCTK_WARN(0, "Allocation problems with trivial_rp") |