From 2512639ec0be67086fadcfee57549f2087ca90bb Mon Sep 17 00:00:00 2001 From: allen Date: Sat, 16 Mar 2002 17:23:58 +0000 Subject: Added static boundary conditions git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/WaveToyC/trunk@67 90eea020-d82d-4da5-bf6e-4ee79ff7632f --- par/wavetoyc_static.par | 45 +++++++++++++++++++++++++++++++++++++++++++++ param.ccl | 1 + src/WaveToy.c | 4 ++++ 3 files changed, 50 insertions(+) create mode 100644 par/wavetoyc_static.par 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", -- cgit v1.2.3