diff options
author | allen <allen@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 1999-05-13 18:03:42 +0000 |
---|---|---|
committer | allen <allen@b7a48df3-cbbf-4440-997f-b4b717c9f7fc> | 1999-05-13 18:03:42 +0000 |
commit | 0d84029ae8ed58a9188c0dc8b01782b7a1e0015d (patch) | |
tree | 0dc21affccbe96424912103e06806894319406db /src | |
parent | d9f8a66b546b7c0c3e6477f03a677301126e7d9f (diff) |
Added boundary conditions (Flat, just to look tidy and be consistent with thw
old version)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/ADMConstraints/trunk@5 b7a48df3-cbbf-4440-997f-b4b717c9f7fc
Diffstat (limited to 'src')
-rw-r--r-- | src/ADMConstraints.F | 12 |
1 files changed, 7 insertions, 5 deletions
diff --git a/src/ADMConstraints.F b/src/ADMConstraints.F index d91cdf7..53c3e96 100644 --- a/src/ADMConstraints.F +++ b/src/ADMConstraints.F @@ -58,6 +58,10 @@ c Perhaps this and others should go into cctk.h integer :: i,j,k,i1,i2,j1,j2,k1,k2 + +c Stencil width used for calculating constraints (needed for outer boundary condition) + integer, dimension(3),parameter :: sw = 1 + REAL :: dx,dy,dz REAL :: pi,det,uxx,uxy,uxz,uyy,uyz,uzz @@ -173,12 +177,10 @@ c ---------------------------------- #include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_undefine.h" -c Synchronize and apply boundaries conditions +c Synchronize and apply flat boundary conditions -c call onebound(ARENA, ADM_ham, 1d0, 1d0, 1d0) -c call onebound(ARENA, ADM_momx,-1d0, 1d0, 1d0) -c call onebound(ARENA, ADM_momy, 1d0,-1d0, 1d0) -c call onebound(ARENA, ADM_momz, 1d0,-1d0, 1d0) + call ApplySym(GH,"admconstraints::admconstraints") + call ApplyFlatBC(GH,sw,"admconstraints::admconstraints") if (constraint_communication.eq.1) then call CCTK_SyncGroup(GH,"adm::ADM_constraints") |