Initial revision

darcs-hash:20010301114010-f6438-12fb8a9ffcc80e86c0a97e37b5b0dae0dbc59b79.gz

33 files changed, 1581 insertions, 0 deletions

diff --git a/CarpetExtra/WaveToyF77/par/convtest.gnuplot b/CarpetExtra/WaveToyF77/par/convtest.gnuplot
new file mode 100644
index 000000000..acd6d1526
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/convtest.gnuplot
@@ -0,0 +1,25 @@
+# initialise
+kx=0.222222222222222
+ky=0.151515151515152
+kz=0.222222222222222
+omega=sqrt(kx**2 + ky**2 + kz**2)
+
+
+
+# without refinement
+
+# solutions
+p "wavetoyf77_periodic/phi.zl" u ($8-1):($1==time&&$3==0?$9:0/0) w lp, "wavetoyf77_periodic_coarse/phi.zl" u (2*($8-1)):(2*$1==time&&$3==0?$9:0/0) w lp, cos((kz*(0.3*x-9) + omega*(0.15*time))*pi) w l
+
+# errors
+p "wavetoyf77_periodic/phi.zl" u ($8-1):($1==time&&$3==0?$9-cos((kz*(0.3*($8-1)-9) + omega*(0.15*time))*pi):0/0) w lp, "wavetoyf77_periodic_coarse/phi.zl" u (2*$8):(2*$1==time&&$3==0?$9-cos((kz*(0.3*(2*($8-1))-9) + omega*(0.15*time))*pi):0/0) w lp
+
+
+
+# with refinement
+
+# solutions
+p "wavetoyf77_periodic_rl2/phi.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "wavetoyf77_periodic_coarse_rl2/phi.zl" u (2*($8-2)):(2*$1==time&&$3==0?$9:0/0) w lp, cos((kz*(0.15*x-9) + omega*(0.075*time))*pi) w l
+
+# errors
+p "wavetoyf77_periodic_rl2/phi.zl" u ($8-2):($1==time&&$3==0?$9-cos((kz*(0.15*($8-2)-9) + omega*(0.075*time))*pi):0/0) w lp, "wavetoyf77_periodic_coarse_rl2/phi.zl" u (2*($8-2)):(2*$1==time&&$3==0?$9-cos((kz*(0.15*(2*($8-2))-9) + omega*(0.075*time))*pi):0/0) w lp
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77.par b/CarpetExtra/WaveToyF77/par/wavetoyf77.par
new file mode 100644
index 000000000..823062ddb
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77.par
@@ -0,0 +1,40 @@
+# /*@@
+#   @file      wavetoyf77.par
+#   @date      2001-01-07
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating minimal evolution
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77.par,v 1.8 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::verbose			= no
+
+Time::dtfac			= 0.5
+
+Carpet::global_nx		= 21
+Carpet::global_ny		= 21
+Carpet::global_nz		= 21
+
+grid::type			= byrange
+grid::xyzmin			= -1.0
+grid::xyzmax			= +1.0
+
+grid::mode			= none
+
+WaveToy::bound			= zero
+
+IO::out_dir			= "wavetoyf77"
+
+IOBasic::outInfo_every		= 1
+#IOBasic::outInfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 1
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+#IOBasic::outScalar_style	= gnuplot
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_carpet.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_carpet.par
new file mode 100644
index 000000000..d03856cba
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_carpet.par
@@ -0,0 +1,39 @@
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= octant
+grid::dxyz			= 0.3
+ 
+IO::out_dir			= "wavetoyf77_carpet"
+
+IOBasic::outinfo_every		= 10
+#IOBasic::outinfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 2
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+IOFlexIO::out3D_every		= 10
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+
+WaveToyF77::bound		= radiation
+
+IOFlexIO::indir3D		= "wavetoyf77_pugh"
+IOFlexIO::in3D_vars		= "wavetoy::phi"
+IOFlexIO::in3D_format		= IEEE
+IOFlexIO::in3D_extension	= "_3d.ieee"
+
+#IOFlexIO::indir3D		= "wavetoyf77_rad"
+#IOFlexIO::in3D_vars		= "wavetoy::phi"
+#IOFlexIO::in3D_format		= HDF5
+#IOFlexIO::in3D_extension	= ".h5"
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_excision.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision.par
new file mode 100644
index 000000000..fb011661c
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision.par
@@ -0,0 +1,38 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_excision.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::terminate		= time
+Cactus::cctk_final_time		= 10.0
+ 
+Time::dtfac			= 0.5
+
+driver::global_nsize		= 35
+
+Carpet::max_refinement_levels	= 4
+
+grid::type			= byspacing
+grid::dxyz			= 0.0625
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_excision"
+
+IOBasic::outinfo_every		= 8
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+IOBasic::outScalar_style	= gnuplot
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi"
+
+WaveToyF77::bound		= radiation
+WaveToyF77::excision_bound	= 1/r
+WaveToyF77::excision_radius	= 0.25
+
+IDScalarWave::initial_data	= 1/r
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl2.par
new file mode 100644
index 000000000..304e3ad2f
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl2.par
@@ -0,0 +1,40 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl2.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::terminate		= time
+Cactus::cctk_final_time		= 10.0
+ 
+Time::dtfac			= 0.5
+
+driver::global_nsize		= 35
+
+Carpet::max_refinement_levels	= 4
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byspacing
+grid::dxyz			= 0.0625
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_excision_rl2"
+
+IOBasic::outinfo_every		= 8
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+IOBasic::outScalar_style	= gnuplot
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi"
+
+WaveToyF77::bound		= radiation
+WaveToyF77::excision_bound	= 1/r
+WaveToyF77::excision_radius	= 0.25
+
+IDScalarWave::initial_data	= 1/r
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl3.par
new file mode 100644
index 000000000..153afb55d
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl3.par
@@ -0,0 +1,40 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_excision_rl3.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::terminate		= time
+Cactus::cctk_final_time		= 10.0
+ 
+Time::dtfac			= 0.5
+
+driver::global_nsize		= 35
+
+Carpet::max_refinement_levels	= 4
+
+CarpetRegrid::refinement_levels	= 3
+
+grid::type			= byspacing
+grid::dxyz			= 0.0625
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_excision_rl3"
+
+IOBasic::outinfo_every		= 8
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+IOBasic::outScalar_style	= gnuplot
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi"
+
+WaveToyF77::bound		= radiation
+WaveToyF77::excision_bound	= 1/r
+WaveToyF77::excision_radius	= 0.25
+
+IDScalarWave::initial_data	= 1/r
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_flat.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_flat.par
new file mode 100644
index 000000000..6f765c6f9
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_flat.par
@@ -0,0 +1,41 @@
+# /*@@
+#   @file      wavetoy_flat.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating flat boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_flat.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 500
+
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= octant
+grid::dxyz			= 0.3
+
+IO::out_dir			= "wavetoyf77_flat"
+
+IOBasic::outInfo_every		= 50
+#IOBasic::outInfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 10
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 10
+IOASCII::out1D_vars		= "wavetoy::phi "
+
+WaveToyF77::bound		= flat
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_none.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_none.par
new file mode 100644
index 000000000..7eca77dd9
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_none.par
@@ -0,0 +1,41 @@
+# /*@@
+#   @file      wavetoy_none.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating no boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_none.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 500
+
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= octant
+grid::dxyz			= 0.3
+
+IO::out_dir			= "wavetoyf77_none"
+
+IOBasic::outInfo_every		= 50
+#IOBasic::outInfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 10
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 10
+IOASCII::out1D_vars		= "wavetoy::phi "
+
+WaveToyF77::bound		= none
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par
new file mode 100644
index 000000000..d9637f1c3
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par
@@ -0,0 +1,57 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par,v 1.11 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 120
+
+Time::dtfac			= 0.5
+
+driver::global_nx		= 34
+driver::global_ny		= 48
+driver::global_nz		= 64
+
+driver::ghost_size		= 2
+
+grid::type			= byrange
+grid::xmin			= -5.1
+grid::ymin			= -6.9
+grid::zmin			= -9.3
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 24
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par
new file mode 100644
index 000000000..8a2f3f1f2
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par
@@ -0,0 +1,58 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par,v 1.12 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+#Carpet::poison_new_timelevels	= yes
+#Carpet::check_for_poison	= yes
+#Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 60
+
+Time::dtfac			= 0.5
+
+driver::global_nx		= 19
+driver::global_ny		= 26
+driver::global_nz		= 34
+
+driver::ghost_size		= 2
+
+grid::type			= byrange
+grid::xmin			= -5.4
+grid::ymin			= -7.2
+grid::zmin			= -9.6
+grid::xmax			=  5.4
+grid::ymax			=  7.8
+grid::zmax			= 10.2
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_coarse"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 12
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par
new file mode 100644
index 000000000..eca2f4784
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par
@@ -0,0 +1,64 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par,v 1.15 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 19
+driver::global_ny		= 26
+driver::global_nz		= 34
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 2
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byrange
+grid::xmin			= -5.4
+grid::ymin			= -7.2
+grid::zmin			= -9.6
+grid::xmax			=  5.4
+grid::ymax			=  7.8
+grid::zmax			= 10.2
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_coarse_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 24
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par
new file mode 100644
index 000000000..00dd3b4d7
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par
@@ -0,0 +1,64 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par,v 1.13 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 19
+driver::global_ny		= 26
+driver::global_nz		= 34
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 3
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+CarpetRegrid::refinement_levels	= 3
+
+grid::type			= byrange
+grid::xmin			= -5.4
+grid::ymin			= -7.2
+grid::zmin			= -9.6
+grid::xmax			=  5.4
+grid::ymax			=  7.8
+grid::zmax			= 10.2
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_coarse_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+IOFlexIO::out3D_every		= 48
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine.par
new file mode 100644
index 000000000..a6c784c72
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine.par
@@ -0,0 +1,57 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 64
+driver::global_ny		= 92
+driver::global_nz		= 124
+
+driver::ghost_size		= 2
+
+grid::type			= byrange
+grid::xmin			= -4.95
+grid::ymin			= -6.75
+grid::zmin			= -9.15
+grid::xmax			=  4.5
+grid::ymax			=  6.9
+grid::zmax			=  9.3
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_fine"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 #2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 48
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl2.par
new file mode 100644
index 000000000..06dd57676
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl2.par
@@ -0,0 +1,63 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl2.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 480
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 64
+driver::global_ny		= 92
+driver::global_nz		= 124
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 2
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byrange
+grid::xmin			= -4.95
+grid::ymin			= -6.75
+grid::zmin			= -9.15
+grid::xmax			=  4.5
+grid::ymax			=  6.9
+grid::zmax			=  9.3
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_fine_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 96
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl3.par
new file mode 100644
index 000000000..d7ebd76ea
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl3.par
@@ -0,0 +1,61 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_fine_rl3.par,v 1.4 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 960
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 64
+driver::global_ny		= 92
+driver::global_nz		= 124
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 3
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+grid::type			= byrange
+grid::xmin			= -4.95
+grid::ymin			= -6.75
+grid::zmin			= -9.15
+grid::xmax			=  4.5
+grid::ymax			=  6.9
+grid::zmax			=  9.3
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_fine_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+IOFlexIO::out3D_every		= 192
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par
new file mode 100644
index 000000000..52d23d1f7
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par
@@ -0,0 +1,63 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par,v 1.13 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 34
+driver::global_ny		= 48
+driver::global_nz		= 64
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 2
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byrange
+grid::xmin			= -5.1
+grid::ymin			= -6.9
+grid::zmin			= -9.3
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+#IOFlexIO::out3D_every		= 48
+#IOFlexIO::out3D_vars		= "wavetoy::phi"
+#IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par
new file mode 100644
index 000000000..d959e8d9d
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par
@@ -0,0 +1,61 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par,v 1.12 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Carpet::poison_new_timelevels	= yes
+Carpet::check_for_poison	= yes
+Carpet::checksum_timelevels	= yes
+
+Cactus::cctk_itlast		= 480
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 34
+driver::global_ny		= 48
+driver::global_nz		= 64
+
+driver::ghost_size		= 2
+
+Carpet::max_refinement_levels	= 3
+Carpet::prolongation_order_space= 3
+Carpet::prolongation_order_time	= 2
+
+grid::type			= byrange
+grid::xmin			= -5.1
+grid::ymin			= -6.9
+grid::zmin			= -9.3
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
+grid::mode			= periodic
+
+IO::out_dir			= "wavetoyf77_periodic_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+IOFlexIO::out3D_every		= 96
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= plane
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_pugh.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_pugh.par
new file mode 100644
index 000000000..5e0f28ff4
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_pugh.par
@@ -0,0 +1,35 @@
+ActiveThorns = "boundary iobasic idscalarwave time wavetoyf77 pugh pughslab pughreduce cartgrid3d ioutil ioascii ioflexio flexio"
+ 
+time::dtfac = 0.5
+
+idscalarwave::initial_data = "gaussian"
+idscalarwave::sigma = 2.8
+idscalarwave::radius = 0
+
+wavetoyf77::bound = "radiation"
+
+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 = 120
+
+IOBasic::outScalar_every = 2
+IOBasic::outScalar_vars = "wavetoy::phi"
+
+IOASCII::out1D_every = 2
+IOASCII::out1D_z = "no"
+IOASCII::out1D_y = "no"
+IOASCII::out1D_vars = "wavetoy::phi"
+
+IOBasic::outinfo_every = 10
+IOBasic::outinfo_vars = "wavetoy::phi"
+
+IOFlexIO::out3D_every = 1
+IOFlexIO::out3D_vars = "wavetoy::phi"
+ 
+IO::out_dir = "wavetoyf77_pugh"
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad.par
new file mode 100644
index 000000000..5f7b87c17
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad.par
@@ -0,0 +1,47 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad.par,v 1.9 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= octant
+grid::dxyz			= 0.3
+ 
+IO::out_dir			= "wavetoyf77_rad"
+
+IOBasic::outinfo_every		= 10
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+IOASCII::separate_grids		= no
+
+IOFlexIO::out3D_every		= 10
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_automatic.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_automatic.par
new file mode 100644
index 000000000..66df240d8
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_automatic.par
@@ -0,0 +1,59 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_automatic.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= none
+grid::dxyz			= 0.6
+
+Carpet::verbose			= yes
+Carpet::max_refinement_levels	= 3
+
+CarpetRegrid::refined_regions	= automatic
+CarpetRegrid::refinement_levels	= 3
+CarpetRegrid::regrid_every	= 40
+CarpetRegrid::minwidth		= 8
+CarpetRegrid::minfraction	= 0.75
+CarpetRegrid::errorvar		= "wavetoy::phi"
+CarpetRegrid::maxerror		= 0.1
+
+IO::out_dir			= "wavetoyf77_rad_automatic"
+
+IOBasic::outinfo_every		= 40
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+IOASCII::separate_grids		= yes
+IOASCII::separate_components	= yes
+
+IOFlexIO::out3D_every		= 40
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_centre.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_centre.par
new file mode 100644
index 000000000..28d0b193b
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_centre.par
@@ -0,0 +1,54 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_centre.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= none
+grid::dxyz			= 0.6
+
+Carpet::verbose			= yes
+Carpet::max_refinement_levels	= 3
+
+CarpetRegrid::refined_regions	= centre
+CarpetRegrid::refinement_levels	= 3
+CarpetRegrid::regrid_every	= 40
+
+IO::out_dir			= "wavetoyf77_rad_centre"
+
+IOBasic::outinfo_every		= 40
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+IOASCII::separate_grids		= no
+
+IOFlexIO::out3D_every		= 40
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full.par
new file mode 100644
index 000000000..4b5ff3046
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full.par
@@ -0,0 +1,41 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 33
+driver::global_ny		= 45
+driver::global_nz		= 61
+
+grid::type			= byspacing
+grid::dxyz			= 0.3
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 1
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse.par
new file mode 100644
index 000000000..2b3041897
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse.par
@@ -0,0 +1,41 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse.par,v 1.7 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 60
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 17
+driver::global_ny		= 23
+driver::global_nz		= 31
+
+grid::type			= byspacing
+grid::dxyz			= 0.6
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_coarse"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl2.par
new file mode 100644
index 000000000..34fe7d655
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl2.par
@@ -0,0 +1,45 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl2.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 17
+driver::global_ny		= 23
+driver::global_nz		= 31
+
+Carpet::max_refinement_levels	= 2
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byspacing
+grid::dxyz			= 0.6
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_coarse_rl2"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl3.par
new file mode 100644
index 000000000..3dd715e08
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl3.par
@@ -0,0 +1,45 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_coarse_rl3.par,v 1.6 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 17
+driver::global_ny		= 23
+driver::global_nz		= 31
+
+Carpet::max_refinement_levels	= 3
+
+CarpetRegrid::refinement_levels	= 3
+
+grid::type			= byspacing
+grid::dxyz			= 0.6
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_coarse_rl3"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse.par
new file mode 100644
index 000000000..194ce0039
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse.par
@@ -0,0 +1,41 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse.par,v 1.4 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 30
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 9
+driver::global_ny		= 12
+driver::global_nz		= 16
+
+grid::type			= byspacing
+grid::dxyz			= 1.2
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_realcoarse"
+
+IOBasic::outinfo_every		= 1 # 10
+#IOBasic::outinfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 2
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl2.par
new file mode 100644
index 000000000..19fa632ca
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl2.par
@@ -0,0 +1,43 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl2.par,v 1.4 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 30
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 9
+driver::global_ny		= 12
+driver::global_nz		= 16
+
+Carpet::max_refinement_levels	= 2
+
+grid::type			= byspacing
+grid::dxyz			= 1.2
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_realcoarse_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+#IOBasic::outinfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 2
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl3.par
new file mode 100644
index 000000000..451503e37
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl3.par
@@ -0,0 +1,43 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_realcoarse_rl3.par,v 1.4 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 30
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 9
+driver::global_ny		= 12
+driver::global_nz		= 16
+
+Carpet::max_refinement_levels	= 3
+
+grid::type			= byspacing
+grid::dxyz			= 1.2
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_realcoarse_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+#IOBasic::outinfo_vars		= "wavetoy::phi"
+
+#IOBasic::outScalar_every	= 2
+#IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl2.par
new file mode 100644
index 000000000..d283afe38
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl2.par
@@ -0,0 +1,45 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl2.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 33
+driver::global_ny		= 45
+driver::global_nz		= 61
+
+Carpet::max_refinement_levels	= 2
+
+CarpetRegrid::refinement_levels	= 2
+
+grid::type			= byspacing
+grid::dxyz			= 0.3
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_rl2"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl3.par
new file mode 100644
index 000000000..11474bbd9
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl3.par
@@ -0,0 +1,45 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_full_rl3.par,v 1.3 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 33
+driver::global_ny		= 45
+driver::global_nz		= 61
+
+Carpet::max_refinement_levels	= 3
+
+CarpetRegrid::refinement_levels	= 3
+
+grid::type			= byspacing
+grid::dxyz			= 0.3
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_rad_full_rl3"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_none.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_none.par
new file mode 100644
index 000000000..e2c07db00
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_none.par
@@ -0,0 +1,53 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_rad_none.par,v 1.5 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO FlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 240
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 30
+driver::global_ny		= 30
+driver::global_nz		= 30
+
+grid::type			= byspacing
+grid::mode			= none
+grid::dxyz			= 0.6
+
+Carpet::max_refinement_levels	= 3
+
+CarpetRegrid::refined_regions	= none
+CarpetRegrid::refinement_levels	= 3
+CarpetRegrid::regrid_every	= 40
+
+IO::out_dir			= "wavetoyf77_rad_none"
+
+IOBasic::outinfo_every		= 40
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 8
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 8
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+IOASCII::separate_grids		= no
+
+IOFlexIO::out3D_every		= 40
+IOFlexIO::out3D_vars		= "wavetoy::phi"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+WaveToyF77::bound		= radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse.par
new file mode 100644
index 000000000..66121b671
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse.par
@@ -0,0 +1,44 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse.par,v 1.3 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 60
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 17
+driver::global_ny		= 23
+driver::global_nz		= 31
+
+Carpet::base_extents		= "[ ([1,1,1]:[15,21,14]:[1,1,1]), ([1,1,18]:[15,21,29]:[1,1,1]) ]"
+Carpet::base_bboxes		= "[ [[1,1],[1,1],[1,1]], [[1,1],[1,1],[1,1]] ]"
+
+grid::type			= byspacing
+grid::dxyz			= 0.6
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_zero_ell_coarse"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= zero # radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse_rl2.par
new file mode 100644
index 000000000..bf73eafc9
--- /dev/null
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse_rl2.par
@@ -0,0 +1,48 @@
+# /*@@
+#   @file      wavetoyf77_rad.par
+#   @date      2001-03-06
+#   @author    Erik Schnetter
+#   @desc
+#   Wavetoy parameter file demonstrating radiation boundaries in octant mode
+#   @enddesc
+# @@*/
+#
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_zero_ell_coarse_rl2.par,v 1.3 2004/03/23 12:00:37 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDScalarWave WaveToyF77"
+
+Cactus::cctk_itlast		= 120
+ 
+Time::dtfac			= 0.5
+
+driver::global_nx		= 17
+driver::global_ny		= 23
+driver::global_nz		= 31
+
+Carpet::base_extents		= "[ ([2,2,2]:[30,42,28]:[2,2,2]), ([2,2,36]:[30,42,58]:[2,2,2]) ]"
+Carpet::base_bboxes		= "[ [[1,1],[1,1],[1,1]], [[1,1],[1,1],[1,1]] ]"
+
+Carpet::max_refinement_levels	= 2
+
+CarpetRegrid::refinement_levels	= 1
+
+grid::type			= byspacing
+grid::dxyz			= 0.6
+grid::avoid_origin		= no
+
+IO::out_dir			= "wavetoyf77_zero_ell_coarse_rl2"
+
+IOBasic::outinfo_every		= 1
+IOBasic::outinfo_vars		= "wavetoy::phi"
+
+IOBasic::outScalar_every	= 2
+IOBasic::outScalar_vars		= "wavetoy::phi"
+
+IOASCII::out1D_every		= 1
+IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
+
+WaveToyF77::bound		= zero # radiation
+
+IDScalarWave::initial_data	= gaussian
+IDScalarWave::sigma		= 2.8
+IDScalarWave::radius		= 0