Make the width of the boundary a Cactus parameter "bound_width".

darcs-hash:20030708211245-07bb3-04251d647d052d857072f63961daafaec3d05c60.gz

2 files changed, 11 insertions, 13 deletions

diff --git a/CarpetExtra/FOWaveToyF77/param.ccl b/CarpetExtra/FOWaveToyF77/param.ccl
index 7b2263a38..9204f6329 100644
--- a/CarpetExtra/FOWaveToyF77/param.ccl
+++ b/CarpetExtra/FOWaveToyF77/param.ccl
@@ -1,5 +1,5 @@
 # Parameter definitions for thorn WaveToyF77
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/FOWaveToyF77/param.ccl,v 1.4 2003/06/27 15:54:58 schnetter Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/FOWaveToyF77/param.ccl,v 1.5 2003/07/08 23:12:45 schnetter Exp $
 
 private:
  
@@ -8,6 +8,11 @@ STRING bound "Type of boundary condition to use"
   ".*" :: "Must be a registered boundary condition"
 } "Zero"
 
+INT bound_width "Number of boundary points"
+{
+  1:* :: ""
+} 1
+
 KEYWORD excision_bound "Type of excision boundary condition to use"
 {
   "1/r"  :: "use 1/r Dirichlet data"
diff --git a/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77 b/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
index 9e6db95bd..829961283 100644
--- a/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
+++ b/CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77
@@ -42,9 +42,6 @@ c     Declare variables in argument list
       DECLARE_CCTK_ARGUMENTS
       DECLARE_CCTK_FUNCTIONS
       DECLARE_CCTK_PARAMETERS
-      
-      INTEGER   bndwidth
-      PARAMETER (bndwidth=2)
 
       INTEGER   i,j,k
       INTEGER   istart(3), iend(3)
@@ -69,12 +66,12 @@ c     -----------------
          if (cctk_bbox(2*d-1).eq.0) then
             istart(d) = 1+cctk_nghostzones(d)
          else
-            istart(d) = 1+bndwidth
+            istart(d) = 1+bound_width
          end if
          if (cctk_bbox(2*d).eq.0) then
             iend(d) = cctk_lsh(d)-cctk_nghostzones(d)
          else
-            iend(d) = cctk_lsh(d)-bndwidth
+            iend(d) = cctk_lsh(d)-bound_width
          end if
       end do
       
@@ -134,12 +131,8 @@ c     ----------------
       CCTK_INT Boundary_SelectGroupForBC 
       
 c     Local declarations
-      INTEGER   bndwidth
-      PARAMETER (bndwidth=2)
-
       integer table
-      save table
-      data table /-1/
+      data    table /-1/
       
       character fbound*100
       CCTK_INT  fboundlen
@@ -196,10 +189,10 @@ c     -----------------------------------
       call CCTK_FortranString (fboundlen, bound, fbound)
       if (fboundlen.lt.0) call CCTK_WARN (0, "internal error")
       
-      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bndwidth,
+      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bound_width,
      $     table, "fowavetoy::scalarevolve", fbound)
       if (ierr.ne.0) call CCTK_WARN (0, "internal error")
-      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bndwidth,
+      ierr = Boundary_SelectGroupForBC (cctkGH, CCTK_ALL_FACES, bound_width,
      $     table, "fowavetoy::scalarevolve_derivs", fbound)
       if (ierr.ne.0) call CCTK_WARN (0, "internal error")