diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-19 23:16:09 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-09-19 23:16:09 +0000 |
| commit | 246882fdfc1f1bd00ea2969193aa665ea6513019 (patch) | |
| tree | 988845c113729ffc1353bcca77ce6468eef2fb7d /src/GRHydro_ENOReconstruct_drv.F90 | |
| parent | bab01622dcacbc47d2c4a5a89549374bb863c9bd (diff) | |
fix non-multipatch support for MHD
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@276 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_ENOReconstruct_drv.F90')
| -rw-r--r-- | src/GRHydro_ENOReconstruct_drv.F90 | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/src/GRHydro_ENOReconstruct_drv.F90 b/src/GRHydro_ENOReconstruct_drv.F90 index 1c2a28a..a46de80 100644 --- a/src/GRHydro_ENOReconstruct_drv.F90 +++ b/src/GRHydro_ENOReconstruct_drv.F90 @@ -20,9 +20,9 @@ #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) -#define Bvecx(i,j,k) lBvec(i,j,k,1) -#define Bvecy(i,j,k) lBvec(i,j,k,2) -#define Bvecz(i,j,k) lBvec(i,j,k,3) +#define Bvecx(i,j,k) Bprim(i,j,k,1) +#define Bvecy(i,j,k) Bprim(i,j,k,2) +#define Bvecz(i,j,k) Bprim(i,j,k,3) #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) @@ -67,12 +67,15 @@ subroutine GRHydro_ENOReconstruct_drv(CCTK_ARGUMENTS) ! save memory when MP is not used CCTK_INT :: GRHydro_UseGeneralCoordinates TARGET vel, lvel - CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup + TARGET Bvec, lBvec + CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup, Bprim if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then vup => lvel + Bprim => lBvec else vup => vel + Bprim => Bvec end if allocate(trivial_rp(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)),STAT=ierr) |