Add the option to use user defined regions for the elliptic smoothing

procedure. Works by setting the admbase variables to zero in these
regions before doing the smoothing. Needs some testing.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/NoExcision/trunk@12 4ec1db94-0e4f-0410-ada3-8bed251432c9

1 files changed, 7 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 54d402b..824d112 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -20,6 +20,13 @@ if (smooth_regions) {
       LANG: Fortran
     } "Smooth regions"
   } else {
+    
+    if (use_user_regions) {
+      SCHEDULE NoExcision_Set_Zero IN ADMBase_PostInitial BEFORE NoExcision_CGSmoothing
+      {
+        LANG: Fortran
+      } "set variablse to zero in used defined regions"
+    }
 
     SCHEDULE GROUP NoExcision_CGSmoothing IN ADMBase_PostInitial
     {