diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-29 17:17:07 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-11-29 17:17:07 +0000 |
commit | 3392abc93c456d7faf1e1c7683bca7460466a82c (patch) | |
tree | 1180abd23e76d4cecfe1d08170d1b7d34b1c06c1 /src/GRHydro_HLLEM.F90 | |
parent | b6d504ec102733365346aa056744aaf0055eb008 (diff) |
expand region where fluxes are computed for constraint transport
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@307 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_HLLEM.F90')
-rw-r--r-- | src/GRHydro_HLLEM.F90 | 8 |
1 files changed, 5 insertions, 3 deletions
diff --git a/src/GRHydro_HLLEM.F90 b/src/GRHydro_HLLEM.F90 index 7e9c541..6f1791e 100644 --- a/src/GRHydro_HLLEM.F90 +++ b/src/GRHydro_HLLEM.F90 @@ -74,9 +74,11 @@ subroutine GRHydro_HLLEM(CCTK_ARGUMENTS) call CCTK_WARN(0, "Flux direction not x,y,z") end if - do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil - do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil - do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + ! constraint transport needs to be able to average fluxes in the directions + ! other that flux_direction + do k = GRHydro_stencil, cctk_lsh(3) - GRHydro_stencil + transport_constraints*(1-zoffset) + do j = GRHydro_stencil, cctk_lsh(2) - GRHydro_stencil + transport_constraints*(1-yoffset) + do i = GRHydro_stencil, cctk_lsh(1) - GRHydro_stencil + transport_constraints*(1-xoffset) f1 = 0.d0 lamminus = 0.d0 |