Added static boundary conditions

git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyC/trunk@67 90eea020-d82d-4da5-bf6e-4ee79ff7632f

3 files changed, 50 insertions, 0 deletions

diff --git a/par/wavetoyc_static.par b/par/wavetoyc_static.par
new file mode 100644
index 0000000..939a92c
--- /dev/null
+++ b/par/wavetoyc_static.par
@@ -0,0 +1,45 @@
+# wavetoyc_static.par - wavetoy evolution with static boundaries
+# $Id$
+#
+# /*@@
+#   @file      wavetoyc_flat.par
+#   @date      Sat 16th March 2002
+#   @author    Gabrielle Allen
+#   @desc
+#   Wavetoy parameter file demonstrating static boundaries in Octant
+#   @enddesc
+# @@*/          
+
+ActiveThorns = "idscalarwavec time boundary wavetoyc pugh pughreduce pughslab cartgrid3d ioutil ioascii iobasic"
+ 
+time::dtfac = 0.5
+
+idscalarwave::initial_data = "gaussian"
+idscalarwave::sigma = 4.8
+idscalarwave::radius = 0
+
+wavetoy::bound = "static"
+
+grid::type = "BySpacing"
+grid::domain = "octant"
+grid::dxyz = 0.3
+
+driver::global_nx = 30
+driver::global_ny = 30
+driver::global_nz = 30
+
+cactus::cctk_itlast = 500
+
+IOBasic::outScalar_every = 10
+IOBasic::outScalar_vars = "wavetoy::phi"
+
+IOASCII::out1D_every = 10
+IOASCII::out1D_z = "no"
+IOASCII::out1D_y = "no"
+IOASCII::out1D_vars = "wavetoy::phi "
+
+IOBasic::outInfo_every = 50
+IOBasic::outInfo_vars = "wavetoy::phi"
+ 
+IO::outdir = "wavetoyc_static"
+
diff --git a/param.ccl b/param.ccl
index 1ca884b..ffb6720 100644
--- a/param.ccl
+++ b/param.ccl
@@ -7,6 +7,7 @@ KEYWORD bound "Type of boundary condition to use"
 {
   "none"      :: "No boundary condition"
   "flat"      :: "Flat boundary condition"
+  "static"    :: "Static boundary condition"
   "radiation" :: "Radiation boundary condition"
 } "none"
 
diff --git a/src/WaveToy.c b/src/WaveToy.c
index 4cd365b..e959cea 100644
--- a/src/WaveToy.c
+++ b/src/WaveToy.c
@@ -134,6 +134,10 @@ void WaveToyC_Boundaries(CCTK_ARGUMENTS)
   {
     ierr = BndFlatVN(cctkGH,sw,"wavetoy::phi");
   }
+  else if (CCTK_Equals(bound,"static"))
+  {
+    ierr = BndStaticVN(cctkGH,sw,"wavetoy::phi");
+  }
   else if (CCTK_Equals(bound,"radiation")) 
   {
     ierr = BndRadiativeVN(cctkGH,sw,0.0,1.0,"wavetoy::phi",