Initial revision

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

9 files changed, 847 insertions, 0 deletions

diff --git a/CarpetExtra/WaveToyMoL/par/sw1d.gnu b/CarpetExtra/WaveToyMoL/par/sw1d.gnu
new file mode 100644
index 000000000..834756543
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d.gnu
@@ -0,0 +1,55 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d.gnu,v 1.2 2003/06/30 17:28:51 schnetter Exp $
+
+set grid
+
+dt=0.05
+f(x)=cos(2*pi*(x+t*dt))
+
+t=200
+
+
+p [0:1] "sw1d_1l_0020/phi.xl" i t u 10:13 w l, "sw1d_1l_0040/phi.xl" i t u 10:13 w l, "sw1d_1l_0080/phi.xl" i t u 10:13 w l, f(x)
+
+p [0:1] "sw1d_1l_0020/phi.xl" i t u 10:($13-f($10)) w l, "sw1d_1l_0040/phi.xl" i t u 10:(4*($13-f($10))) w l, "sw1d_1l_0080/phi.xl" i t u 10:(16*($13-f($10))) w l
+
+
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:13 w l, f(x)
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:($13-f($10)) w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:(4*($13-f($10))) w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:(16*($13-f($10))) w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:(64*($13-f($10))) w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:(256*($13-f($10))) w l
+
+
+
+
+
+set grid
+
+
+dt=0.1
+A=1.0/16.0
+
+f(x)=A*t*dt*sin(2*pi*x)*cos(2*pi*dt*t)
+
+
+
+t=400
+
+p [0:1] "sw1d_1l_0020/phierror.xl" i t u 10:13 w lp, "sw1d_1l_0040/phierror.xl" i 2*t u 10:(4*$13) w lp, "sw1d_1l_0080/phierror.xl" i 4*t u 10:(16*$13) w lp
+
+
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 3*t ev :::2::2 u 10:13 w l, "sw1d_2l_0040/phi.xl" i 3*2*t ev :::2::2 u 10:13 w l, "sw1d_2l_0080/phi.xl" i 3*4*t ev :::2::2 u 10:13 w l, "sw1d_2l_0160/phi.xl" i 3*8*t ev :::2::2 u 10:13 w l, "sw1d_2l_0320/phi.xl" i 3*16*t ev :::2::2 u 10:13 w l, cos(2*pi*x)
+p [0:1] "sw1d_2l_20/phierror.xl" i 3*t u 10:13 w l, "sw1d_2l_40/phierror.xl" i 3*2*t u 10:(4*$13) w l, "sw1d_2l_80/phierror.xl" i 3*4*t u 10:(16*$13) w l, "sw1d_2l_160/phierror.xl" i 3*8*t u 10:(64*$13) w l, "sw1d_2l_320/phierror.xl" i 3*16*t u 10:(256*$13) w l
+
+p [0:1] "sw1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "sw1d_2l_0320/phi.xl" i 2*t u 10:13 w l, cos(2*pi*x)
+p [0:1] "sw1d_2l_0020/phierror.xl" i 2*t u 10:13 w l, "sw1d_2l_0040/phierror.xl" i 2*t u 10:(4*$13) w l, "sw1d_2l_0080/phierror.xl" i 2*t u 10:(16*$13) w l, "sw1d_2l_0160/phierror.xl" i 2*t u 10:(64*$13) w l, "sw1d_2l_0320/phierror.xl" i 2*t u 10:(256*$13) w l
+p [0:1] "sw1d_2l_0020/phierror.xl" i 2*t u 10:($13-f($10)/20) w l, "sw1d_2l_0040/phierror.xl" i 2*t u 10:(4*($13-f($10)/40)) w l, "sw1d_2l_0080/phierror.xl" i 2*t u 10:(16*($13-f($10)/80)) w l, "sw1d_2l_0160/phierror.xl" i 2*t u 10:(64*($13-f($10)/160)) w l, "sw1d_2l_0320/phierror.xl" i 2*t u 10:(256*($13-f($10)/320)) w l
+
+
+
+p [:] "< tail +6 sw1d_1l_0020/phierror_norm2.xg" t "unigrid error, dx=1/20" w l, "< tail +6 sw1d_1l_0040/phierror_norm2.xg" u 1:(4*$2) t "4 * unigrid error, dx=1/40" w l, "< tail +6 sw1d_2l_0020/phierror_norm2.xg" t "FMR error, dx=1/20" w l, "< tail +6 sw1d_2l_0040/phierror_norm2.xg" u 1:(4*$2) t "4 * FMR error, dx=1/40" w l, "< tail +6 sw1d_2l_0080/phierror_norm2.xg" u 1:(16*$2) t "16 * FMR error, dx=1/80" w l, "< tail +6 sw1d_2l_0160/phierror_norm2.xg" u 1:(64*$2) t "64 * FMR error, dx=1/160" w l, "< tail +6 sw1d_2l_0320/phierror_norm2.xg" u 1:(256*$2) t "256 * FMR error, dx=1/320" w l
+, "< tail +6 sw1d_1l_0080/phierror_norm2.xg" u 1:(16*$2) t "16 * unigrid error, dx=1/80" w l
+
+p [:5] "< tail +6 sw1d_1l_20/phierror_norm2.xg" t "unigrid error, dx=1/20" w l, "< tail +6 sw1d_1l_40/phierror_norm2.xg" u 1:(4*$2) t "4 * unigrid error, dx=1/40" w l, "< tail +6 sw1d_1l_80/phierror_norm2.xg" u 1:(16*$2) t "16 * unigrid error, dx=1/80" w l, "< tail +6 sw1d_2l_20/phierror_norm2.xg" ev 2 t "FMR error, dx=1/20" w l, "< tail +6 sw1d_2l_40/phierror_norm2.xg" ev 2 u 1:(4*$2) t "4 * FMR error, dx=1/40" w l, "< tail +6 sw1d_2l_80/phierror_norm2.xg" ev 2 u 1:(16*$2) t "16 * FMR error, dx=1/80" w l, "< tail +6 sw1d_2l_160/phierror_norm2.xg" ev 2 u 1:(64*$2) t "64 * FMR error, dx=1/160" w l, "< tail +6 sw1d_2l_320/phierror_norm2.xg" ev 2 u 1:(256*$2) t "256 * FMR error, dx=1/320" w l
+
+p [:5] "< tail +6 sw1d_1l_20/psierror_norm2.xg" w l, "< tail +6 sw1d_1l_40/psierror_norm2.xg" u 1:(4*$2) w l, "< tail +6 sw1d_1l_80/psierror_norm2.xg" u 1:(16*$2) w l, "< tail +6 sw1d_2l_20/psierror_norm2.xg" ev 2 w l, "< tail +6 sw1d_2l_40/psierror_norm2.xg" ev 2 u 1:(4*$2) w l, "< tail +6 sw1d_2l_80/psierror_norm2.xg" ev 2 u 1:(16*$2) w l, "< tail +6 sw1d_2l_160/psierror_norm2.xg" ev 2 u 1:(64*$2) w l
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par
new file mode 100644
index 000000000..7cd2fb167
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0020.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([28,4,4]:[68,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 16
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.1
+grid::ymin = -0.15
+grid::zmin = -0.15
+grid::xmax = +1.05
+grid::ymax = +0.15
+grid::zmax = +0.15
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par
new file mode 100644
index 000000000..9c9f93248
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0040.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([48,4,4]:[128,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.05
+grid::zmin = -0.05
+grid::xmax = +1.025
+grid::ymax = +0.05
+grid::zmax = +0.05
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par
new file mode 100644
index 000000000..c8a040c7c
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_1l_0080.par,v 1.4 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 1
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([88,4,4]:[248,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.025
+grid::zmin = -0.025
+grid::xmax = +1.0125
+grid::ymax = +0.025
+grid::zmax = +0.025
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_1l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par
new file mode 100644
index 000000000..13f2de3ee
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0020.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([30,4,4]:[66,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 16
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.1
+grid::ymin = -0.1
+grid::zmin = -0.1
+grid::xmax = +1.05
+grid::ymax = +0.1
+grid::zmax = +0.1
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par
new file mode 100644
index 000000000..5f215164d
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0040.par,v 1.6 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([50,4,4]:[126,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "CoordBase LocalInterp Slab SymBase"
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.05
+grid::zmin = -0.05
+grid::xmax = +1.025
+grid::ymax = +0.05
+grid::zmax = +0.05
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par
new file mode 100644
index 000000000..5bfcdf1cc
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0080.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([90,4,4]:[246,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.025
+grid::zmin = -0.025
+grid::xmax = +1.0125
+grid::ymax = +0.025
+grid::zmax = +0.025
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par
new file mode 100644
index 000000000..c47ead004
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0160.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 25600
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 164
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([170,4,4]:[486,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 128
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.0125
+grid::ymin = -0.0125
+grid::zmin = -0.0125
+grid::xmax = +1.00625
+grid::ymax = +0.0125
+grid::zmax = +0.0125
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0160"
+
+IOBasic::outInfo_every = 128
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 128
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 128
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 128
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
diff --git a/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par
new file mode 100644
index 000000000..4836c42b9
--- /dev/null
+++ b/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyMoL/par/sw1d_2l_0320.par,v 1.5 2004/07/07 16:15:31 tradke Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 51200
+
+
+
+ActiveThorns = "IOUtil Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 324
+driver::global_ny = 5
+driver::global_nz = 5
+driver::ghost_size = 2
+
+#Carpet::verbose = yes
+#Carpet::veryverbose = yes
+Carpet::max_refinement_levels = 3
+Carpet::prolongation_order_space = 3
+Carpet::prolongation_order_time = 2
+#Carpet::init_3_timelevels = yes
+Carpet::init_each_timelevel = yes
+
+#CarpetLib::output_bboxes = yes
+
+CarpetRegrid::refinement_levels = 2
+CarpetRegrid::refined_regions = manual-gridpoint-list
+CarpetRegrid::gridpoints = "[[ ([330,4,4]:[966,12,12]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 256
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Boundary CoordBase SymBase CartGrid3D Periodic"
+
+grid::domain = full
+grid::type = byrange
+grid::xmin = -0.00625
+grid::ymin = -0.00625
+grid::zmin = -0.00625
+grid::xmax = +1.003125
+grid::ymax = +0.00625
+grid::zmax = +0.00625
+
+Periodic::periodic = yes
+
+
+
+ActiveThorns = "Time MoL"
+
+Time::dtfac = 0.25
+
+MoL::ODE_Method = ICN
+
+#MoL::ODE_Method = RK2
+#MoL::MoL_Intermediate_Steps = 2
+
+
+
+ActiveThorns = "IDScalarWaveMoL WaveToyMoL"
+
+IDScalarWaveMoL::wave_number[0] = 1.0
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "WaveToyMoL::scalarevolve"
+
+
+
+ActiveThorns = "IOBasic CarpetIOASCII"
+
+IO::out_dir = "sw1d_2l_0320"
+
+IOBasic::outInfo_every = 256
+IOBasic::outInfo_vars = "WaveToyMoL::phi IDScalarWaveMoL::phierror"
+ 
+IOBasic::outScalar_every = 256
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 256
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"
+
+#IOASCII::out3D_every = 256
+#IOASCII::out3D_vars = "WaveToyMoL::scalarevolve IDScalarWaveMoL::scalarevolveerror"