Initial revision

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

26 files changed, 1852 insertions, 0 deletions

diff --git a/CarpetExtra/SpaceTimeToy/par/convtest.gnuplot b/CarpetExtra/SpaceTimeToy/par/convtest.gnuplot
new file mode 100644
index 000000000..a1042c939
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/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 "spacetimetoy_periodic/phi.zl" u ($8-1):($1==time&&$3==0?$9:0/0) w lp, "spacetimetoy_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 "spacetimetoy_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, "spacetimetoy_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 "spacetimetoy_periodic_rl2/phi.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "spacetimetoy_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 "spacetimetoy_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, "spacetimetoy_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/SpaceTimeToy/par/convtest_double.gnuplot b/CarpetExtra/SpaceTimeToy/par/convtest_double.gnuplot
new file mode 100644
index 000000000..660430e24
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/convtest_double.gnuplot
@@ -0,0 +1,36 @@
+# initialise
+kx=0.222222222222222
+ky=0.151515151515152
+kz=0.222222222222222
+omega=sqrt(kx**2 + ky**2 + kz**2)
+
+
+
+# compare to uncoupled systems
+
+# coarse solutions
+p "spacetimetoy_periodic_coarse/phi.zl" u (2*($8-1)):(2*$1==time&&$3==0?$9:0/0) w lp, "hydrotoy_periodic_coarse/u.zl" u (2*($8-1)):(2*$1==time&&$3==0?$9:0/0) w lp, "doubletoy_periodic_coarse/phi.zl" u (2*($8-1)):(2*$1==time&&$3==0?$9:0/0) w lp, "doubletoy_periodic_coarse/u.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
+
+# fine solutions
+p "spacetimetoy_periodic/phi.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "hydrotoy_periodic/u.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "doubletoy_periodic/phi.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "doubletoy_periodic/u.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, cos((kz*(0.3*x-9) + omega*(0.15*time))*pi) w l
+
+
+
+# without refinement
+
+# solutions
+
+p "doubletoy_periodic/phi.zl" u ($8-1):($1==time&&$3==0?$9:0/0) w lp, "doubletoy_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 "doubletoy_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, "doubletoy_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 "doubletoy_periodic_rl2/phi.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "doubletoy_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 "doubletoy_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, "doubletoy_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/SpaceTimeToy/par/doubletoy_periodic.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic.par
new file mode 100644
index 000000000..7210fc79a
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic.par
@@ -0,0 +1,66 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic.par,v 1.7 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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			= "doubletoy_periodic"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 24
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse.par
new file mode 100644
index 000000000..2985234ef
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse.par
@@ -0,0 +1,66 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic_coarse.par,v 1.7 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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			= "doubletoy_periodic_coarse"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 12
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl2.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl2.par
new file mode 100644
index 000000000..55534b4a5
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl2.par
@@ -0,0 +1,72 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic_coarse_rl2.par,v 1.8 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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			= "doubletoy_periodic_coarse_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 24
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl3.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl3.par
new file mode 100644
index 000000000..2b7442bce
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_coarse_rl3.par
@@ -0,0 +1,72 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic_coarse_rl3.par,v 1.8 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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			= "doubletoy_periodic_coarse_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 48
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl2.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl2.par
new file mode 100644
index 000000000..2611c62c4
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl2.par
@@ -0,0 +1,72 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic_rl2.par,v 1.8 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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			= "doubletoy_periodic_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 48
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl3.par b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl3.par
new file mode 100644
index 000000000..5cbcc1290
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/doubletoy_periodic_rl3.par
@@ -0,0 +1,72 @@
+# /*@@
+#   @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/SpaceTimeToy/par/doubletoy_periodic_rl3.par,v 1.8 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy IDHydroToy HydroToy"
+
+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
+
+CarpetRegrid::refinement_levels	= 3
+
+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			= "doubletoy_periodic_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 96
+IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
+
+HydroToy::bound			= radiation
+
+IDHydroToy::initial_data	= plane
+IDHydroToy::kx			= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDHydroToy::ky			= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDHydroToy::kz			= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu
new file mode 100644
index 000000000..b78bdf751
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu
@@ -0,0 +1,30 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.gnu,v 1.1 2003/06/18 18:24:29 schnetter Exp $
+
+set grid
+set style data linespoints
+
+kx=0.222222222222222
+ky=0.151515151515152
+kz=0.222222222222222
+omega=sqrt(kx**2+ky**2+kz**2)
+dt=0.3
+phi(x)=cos(pi*(kx*x+omega*dt*t))
+psi(x)=-pi*omega*sin(pi*(kx*x+omega*dt*t))
+
+t=0
+
+
+
+# unigrid function
+p [-4.5:4.5] "spacetimetoy_periodic_coarse/phi.xl" i t u 10:13, "spacetimetoy_periodic/phi.xl" i 2*t u 10:13, "spacetimetoy_periodic_fine/phi.xl" i 4*t u 10:13, phi(x)
+
+# unigrid error
+p [-4.5:4.5] "spacetimetoy_periodic_coarse/phi.xl" i t u 10:($13-phi($10)), "spacetimetoy_periodic/phi.xl" i 2*t u 10:(4*($13-phi($10))), "spacetimetoy_periodic_fine/phi.xl" i 4*t u 10:(16*($13-phi($10)))
+
+
+
+# refinement function
+p [-4.5:4.5] "spacetimetoy_periodic_coarse_rl2/phi.xl" i 3*t u 10:13, "spacetimetoy_periodic_rl2/phi.xl" i 6*t u 10:13, "spacetimetoy_periodic_fine_rl2/phi.xl" i 12*t u 10:13, phi(x)
+
+# refinement error
+p [-4.5:4.5] "spacetimetoy_periodic_coarse_rl2/phi.xl" i 3*t u 10:($13-phi($10)), "spacetimetoy_periodic_rl2/phi.xl" i 6*t u 10:(4*($13-phi($10))), "spacetimetoy_periodic_fine_rl2/phi.xl" i 12*t u 10:(16*($13-phi($10)))
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.par
new file mode 100644
index 000000000..6f157a979
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic.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/SpaceTimeToy/par/spacetimetoy_periodic.par,v 1.9 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+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			= "spacetimetoy_periodic"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 24
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse.par
new file mode 100644
index 000000000..ca078c0d4
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse.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/SpaceTimeToy/par/spacetimetoy_periodic_coarse.par,v 1.11 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+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			= "spacetimetoy_periodic_coarse"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 12
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl2.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl2.par
new file mode 100644
index 000000000..c76217944
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl2.par
@@ -0,0 +1,67 @@
+# /*@@
+#   @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/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl2.par,v 1.11 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+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			= "spacetimetoy_periodic_coarse_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 24
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl3.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl3.par
new file mode 100644
index 000000000..0b155bece
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl3.par
@@ -0,0 +1,67 @@
+# /*@@
+#   @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/SpaceTimeToy/par/spacetimetoy_periodic_coarse_rl3.par,v 1.12 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+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			= "spacetimetoy_periodic_coarse_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 48
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine.par
new file mode 100644
index 000000000..9bd3fb198
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine.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/SpaceTimeToy/par/spacetimetoy_periodic_fine.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+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			= "spacetimetoy_periodic_fine"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 48
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine_rl2.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine_rl2.par
new file mode 100644
index 000000000..ca19dbbcc
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_fine_rl2.par
@@ -0,0 +1,67 @@
+# /*@@
+#   @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/SpaceTimeToy/par/spacetimetoy_periodic_fine_rl2.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+#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			= "spacetimetoy_periodic_fine_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 96
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl2.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl2.par
new file mode 100644
index 000000000..8e4c45437
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl2.par
@@ -0,0 +1,67 @@
+# /*@@
+#   @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/SpaceTimeToy/par/spacetimetoy_periodic_rl2.par,v 1.11 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+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			= "spacetimetoy_periodic_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 48
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl3.par b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl3.par
new file mode 100644
index 000000000..54c264ed8
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/spacetimetoy_periodic_rl3.par
@@ -0,0 +1,67 @@
+# /*@@
+#   @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/SpaceTimeToy/par/spacetimetoy_periodic_rl3.par,v 1.10 2004/03/23 12:00:15 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Slab Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab IDSpaceTimeToy SpaceTimeToy Hydrotoy"
+
+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
+
+CarpetRegrid::refinement_levels	= 3
+
+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			= "spacetimetoy_periodic_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "spacetimetoy::spacetimeevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "spacetimetoy::spacetimeevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "spacetimetoy::spacetimeevolve grid::coordinates"
+
+#IOFlexIO::out3D_every		= 96
+#IOFlexIO::out3D_vars		= "spacetimetoy::spacetimeevolve"
+#IOFlexIO::out3D_format		= IEEE
+#IOFlexIO::out3D_extension	= ".vcamr"
+
+SpaceTimeToy::icn_iterations	= 2
+SpaceTimeToy::bound		= radiation
+
+HydroToy::hydrotoy_active	= no
+
+IDSpaceTimeToy::initial_data	= plane
+IDSpaceTimeToy::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDSpaceTimeToy::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDSpaceTimeToy::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d.gnu b/CarpetExtra/SpaceTimeToy/par/st1d.gnu
new file mode 100644
index 000000000..1c7eee34b
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d.gnu
@@ -0,0 +1,31 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d.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] "st1d_1l_0020/phi.xl" i t u 10:13 w l, "st1d_1l_0040/phi.xl" i t u 10:13 w l, "st1d_1l_0080/phi.xl" i t u 10:13 w l, f(x)
+
+p [0:1] "st1d_1l_0020/phi.xl" i t u 10:($13-f($10)) w l, "st1d_1l_0040/phi.xl" i t u 10:(4*($13-f($10))) w l, "st1d_1l_0080/phi.xl" i t u 10:(16*($13-f($10))) w l
+
+
+
+p [0:1] "st1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0320/phi.xl" i 2*t u 10:13 w l, f(x)
+
+p [0:1] "st1d_2l_0020/phi.xl" i 2*t u 10:($13-f($10)) w l, "st1d_2l_0040/phi.xl" i 2*t u 10:(4*($13-f($10))) w l, "st1d_2l_0080/phi.xl" i 2*t u 10:(16*($13-f($10))) w l, "st1d_2l_0160/phi.xl" i 2*t u 10:(64*($13-f($10))) w l, "st1d_2l_0320/phi.xl" i 2*t u 10:(256*($13-f($10))) w l
+
+
+
+set grid
+
+dt=0.05
+
+p [0:1] "st1d_1l_0020/phi.xl" i t u 10:13 w l, "st1d_1l_0040/phi.xl" i t u 10:13 w l, "st1d_1l_0080/phi.xl" i t u 10:13 w l, cos (2*pi*(x+dt*t))
+p [0:1] "st1d_2l_0020/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0040/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0080/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0160/phi.xl" i 2*t u 10:13 w l, "st1d_2l_0320/phi.xl" i 2*t u 10:13 w l, cos (2*pi*(x+dt*t))
+
+p [0:1] "st1d_1l_0020/phi.xl" i t u 10:($13-cos(2*pi*($10+dt*t))) w l, "st1d_1l_0040/phi.xl" i t u 10:(4*($13-cos(2*pi*($10+dt*t)))) w l, "st1d_1l_0080/phi.xl" i t u 10:(16*($13-cos(2*pi*($10+dt*t)))) w l
+p [0:1] "st1d_2l_0020/phi.xl" i 2*t u 10:($13-cos(2*pi*($10+dt*t))) w l, "st1d_2l_0040/phi.xl" i 2*t u 10:(4*($13-cos(2*pi*($10+dt*t)))) w l, "st1d_2l_0080/phi.xl" i 2*t u 10:(16*($13-cos(2*pi*($10+dt*t)))) w l, "st1d_2l_0160/phi.xl" i 2*t u 10:(64*($13-cos(2*pi*($10+dt*t)))) w l, "st1d_2l_0320/phi.xl" i 2*t u 10:(256*($13-cos(2*pi*($10+dt*t)))) w l
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_1l_0020.par b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0020.par
new file mode 100644
index 000000000..96da55975
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_1l_0020.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 7
+driver::global_nz = 7
+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,6,6]:[68,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+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
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_1l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_1l_0040.par b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0040.par
new file mode 100644
index 000000000..ca8d2fa85
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_1l_0040.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 7
+driver::global_nz = 7
+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,6,6]:[128,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.075
+grid::zmin = -0.075
+grid::xmax = +1.025
+grid::ymax = +0.075
+grid::zmax = +0.075
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_1l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_1l_0080.par b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0080.par
new file mode 100644
index 000000000..a2535989b
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_1l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_1l_0080.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 128
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 7
+driver::global_nz = 7
+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,6,6]:[248,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.0375
+grid::zmin = -0.0375
+grid::xmax = +1.0125
+grid::ymax = +0.0375
+grid::zmax = +0.0375
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_1l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_2l_0020.par b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0020.par
new file mode 100644
index 000000000..0d4263611
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0020.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_2l_0020.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 3200
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 32
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 24
+driver::global_ny = 7
+driver::global_nz = 7
+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 = "[[ ([28,6,6]:[68,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+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
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_2l_0020"
+
+IOBasic::outInfo_every = 16
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 16
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 16
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 16
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_2l_0040.par b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0040.par
new file mode 100644
index 000000000..90261f6aa
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0040.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_2l_0040.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 6400
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 64
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 44
+driver::global_ny = 7
+driver::global_nz = 7
+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 = "[[ ([48,6,6]:[128,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.05
+grid::ymin = -0.075
+grid::zmin = -0.075
+grid::xmax = +1.025
+grid::ymax = +0.075
+grid::zmax = +0.075
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_2l_0040"
+
+IOBasic::outInfo_every = 32
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 32
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 32
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 32
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_2l_0080.par b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0080.par
new file mode 100644
index 000000000..ab4ff726c
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0080.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_2l_0080.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 12800
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 128
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 84
+driver::global_ny = 7
+driver::global_nz = 7
+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 = "[[ ([88,6,6]:[248,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.025
+grid::ymin = -0.0375
+grid::zmin = -0.0375
+grid::xmax = +1.0125
+grid::ymax = +0.0375
+grid::zmax = +0.0375
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_2l_0080"
+
+IOBasic::outInfo_every = 64
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 64
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 64
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 64
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_2l_0160.par b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0160.par
new file mode 100644
index 000000000..cb222029e
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0160.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_2l_0160.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 25600
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 256
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 164
+driver::global_ny = 7
+driver::global_nz = 7
+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 = "[[ ([168,6,6]:[488,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.0125
+grid::ymin = -0.01875
+grid::zmin = -0.01875
+grid::xmax = +1.00625
+grid::ymax = +0.01875
+grid::zmax = +0.01875
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_2l_0160"
+
+IOBasic::outInfo_every = 128
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 128
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 128
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 128
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"
diff --git a/CarpetExtra/SpaceTimeToy/par/st1d_2l_0320.par b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0320.par
new file mode 100644
index 000000000..054046275
--- /dev/null
+++ b/CarpetExtra/SpaceTimeToy/par/st1d_2l_0320.par
@@ -0,0 +1,99 @@
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/SpaceTimeToy/par/st1d_2l_0320.par,v 1.2 2004/03/23 12:00:15 schnetter Exp $
+
+Cactus::cctk_initial_time = 0
+Cactus::cctk_itlast = 51200
+
+
+
+ActiveThorns = "LocalInterp NaNChecker Slab"
+
+NaNChecker::check_every = 512
+NaNChecker::check_vars = "all"
+NaNChecker::action_if_found = terminate
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetRegrid CarpetInterp CarpetReduce CarpetSlab"
+
+driver::global_nx = 324
+driver::global_ny = 7
+driver::global_nz = 7
+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 = "[[ ([328,6,6]:[968,16,16]:[2,2,2]) ]]"
+CarpetRegrid::outerbounds = "[[ [[0,0],[1,1],[1,1]] ]]"
+
+
+
+ActiveThorns = "Boundary Cart3d"
+
+grid::mode = periodic
+grid::type = byrange
+grid::xmin = -0.00625
+grid::ymin = -0.009375
+grid::zmin = -0.009375
+grid::xmax = +1.003125
+grid::ymax = +0.009375
+grid::zmax = +0.009375
+
+
+
+ActiveThorns = "Time"
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "IDSpaceTimeToy SpaceTimeToy HydroToy"
+
+IDSpaceTimeToy::initial_data = plane
+IDSpaceTimeToy::kx = 2.0
+IDSpaceTimeToy::ky = 0.0
+IDSpaceTimeToy::kz = 0.0
+
+SpaceTimeToy::icn_iterations = 2
+
+HydroToy::hydrotoy_active = no
+
+
+
+#ActiveThorns = "Dissipation"
+#
+#Dissipation::epsdis = 0.1
+#Dissipation::timelevel = "n+1"
+#Dissipation::vars = "SpaceTimeToy::spacetimeevolve"
+
+
+
+ActiveThorns = "IOUtil IOBasic CarpetIOASCII"
+
+IO::out_dir = "st1d_2l_0320"
+
+IOBasic::outInfo_every = 256
+IOBasic::outInfo_vars = "SpaceTimeToy::phi"
+ 
+IOBasic::outScalar_every = 256
+IOBasic::outScalar_reductions = "minimum maximum norm1 norm2 norm_inf"
+IOBasic::outScalar_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::output_all_timelevels = yes
+
+IOASCII::out1D_every = 256
+IO::out_yline_x = 0.5
+IO::out_zline_x = 0.5
+IOASCII::out1D_vars = "SpacetimeToy::spacetimeevolve"
+
+#IOASCII::out3D_every = 256
+#IOASCII::out3D_vars = "SpacetimeToy::spacetimeevolve"