From ca9ef5432b81a9678d2b8003465cd66421d7e87a Mon Sep 17 00:00:00 2001 From: schnetter <> Date: Tue, 8 Jul 2003 21:12:00 +0000 Subject: Make the width of the boundary a Cactus parameter "bound_width". darcs-hash:20030708211245-07bb3-04251d647d052d857072f63961daafaec3d05c60.gz --- CarpetExtra/FOWaveToyF77/param.ccl | 7 ++++++- CarpetExtra/FOWaveToyF77/src/FOWaveToy.F77 | 17 +++++------------ 2 files changed, 11 insertions(+), 13 deletions(-) (limited to 'CarpetExtra/FOWaveToyF77') 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") -- cgit v1.2.3