aboutsummaryrefslogtreecommitdiff
path: root/src/GRHydro_UpdateMask.F90
diff options
context:
space:
mode:
authorrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2014-04-15 19:49:55 +0000
committerrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2014-04-15 19:49:55 +0000
commite2e9ee82855cad7b1d32d83ec0a599b0555d031d (patch)
tree5bb9c592ce14dab89871db6bb6ed93181393a0f5 /src/GRHydro_UpdateMask.F90
parente98f136dc52b3c92a8627c3ce6de44e2becc394e (diff)
GRHydro: clean up handling of Avec evolution and boundaries
Remove almost empty Analysis_Init routine, apply boundary conditions consistently. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@635 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_UpdateMask.F90')
-rw-r--r--src/GRHydro_UpdateMask.F906
1 files changed, 6 insertions, 0 deletions
diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90
index aa7d340..9dd8a3e 100644
--- a/src/GRHydro_UpdateMask.F90
+++ b/src/GRHydro_UpdateMask.F90
@@ -65,6 +65,9 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS)
if (evolve_Bvec) then
Bconsrhs(i,j,k,:) = 0.0d0
endif
+ if (evolve_Avec) then
+ Avecrhs(i,j,k,:) = 0.0d0
+ endif
! TODO: Need to set Avecrhs and Aphirhs to zero for centered / or staggered case if vector potential is evolved
@@ -93,6 +96,9 @@ subroutine GRHydroUpdateAtmosphereMask(CCTK_ARGUMENTS)
if (evolve_Bvec) then
Bconsrhs(i,j,k,:) = 0.0d0
endif
+ if (evolve_Avec) then
+ Avecrhs(i,j,k,:) = 0.0d0
+ endif
! TODO: Need to set Avecrhs and Aphirhs to zero for centered / or staggered case if vector potential is evolved