Initial revision

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

7 files changed, 397 insertions, 0 deletions

diff --git a/CarpetExtra/HydroToy/par/convtest.gnuplot b/CarpetExtra/HydroToy/par/convtest.gnuplot
new file mode 100644
index 000000000..15eb3716b
--- /dev/null
+++ b/CarpetExtra/HydroToy/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 "hydrotoy_periodic/u.zl" u ($8-1):($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, cos((kz*(0.3*x-9) + omega*(0.15*time))*pi) w l
+
+# errors
+p "hydrotoy_periodic/u.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, "hydrotoy_periodic_coarse/u.zl" u (2*($8-1)):(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 "hydrotoy_periodic_rl2/u.zl" u ($8-2):($1==time&&$3==0?$9:0/0) w lp, "hydrotoy_periodic_coarse_rl2/u.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 "hydrotoy_periodic_rl2/u.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, "hydrotoy_periodic_coarse_rl2/u.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/HydroToy/par/hydrotoy_periodic.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic.par
new file mode 100644
index 000000000..19da43136
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_periodic.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/HydroToy/par/hydrotoy_periodic.par,v 1.8 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab 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			= "hydrotoy_periodic"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 24
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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/HydroToy/par/hydrotoy_periodic_coarse.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic_coarse.par
new file mode 100644
index 000000000..7a27a685e
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_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/HydroToy/par/hydrotoy_periodic_coarse.par,v 1.7 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab 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			= "hydrotoy_periodic_coarse"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 12
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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/HydroToy/par/hydrotoy_periodic_coarse_rl2.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic_coarse_rl2.par
new file mode 100644
index 000000000..fe9ac250a
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_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/HydroToy/par/hydrotoy_periodic_coarse_rl2.par,v 1.8 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab 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			= "hydrotoy_periodic_coarse_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 24
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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/HydroToy/par/hydrotoy_periodic_coarse_rl3.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic_coarse_rl3.par
new file mode 100644
index 000000000..4229b977f
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_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/HydroToy/par/hydrotoy_periodic_coarse_rl3.par,v 1.8 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab 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			= "hydrotoy_periodic_coarse_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 48
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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/HydroToy/par/hydrotoy_periodic_rl2.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic_rl2.par
new file mode 100644
index 000000000..f9612d03c
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_periodic_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/HydroToy/par/hydrotoy_periodic_rl2.par,v 1.8 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduc CarpetRegrid CarpetSlab 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			= "hydrotoy_periodic_rl2"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 48
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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/HydroToy/par/hydrotoy_periodic_rl3.par b/CarpetExtra/HydroToy/par/hydrotoy_periodic_rl3.par
new file mode 100644
index 000000000..a05ac4719
--- /dev/null
+++ b/CarpetExtra/HydroToy/par/hydrotoy_periodic_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/HydroToy/par/hydrotoy_periodic_rl3.par,v 1.8 2004/03/23 12:00:04 schnetter Exp $
+
+ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetIOFlexIO CarpetLib CarpetReduce CarpetRegrid CarpetSlab 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	= 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			= "hydrotoy_periodic_rl3"
+
+IOBasic::outinfo_every		= 1 # 10
+IOBasic::outinfo_vars		= "hydrotoy::hydroevolve"
+
+IOBasic::outScalar_every	= 1 # 2
+IOBasic::outScalar_vars		= "hydrotoy::hydroevolve"
+
+IOASCII::out1D_every		= 1 # 2
+IOASCII::out1D_vars		= "hydrotoy::hydroevolve grid::coordinates"
+
+IOFlexIO::out3D_every		= 96
+IOFlexIO::out3D_vars		= "hydrotoy::hydroevolve"
+IOFlexIO::out3D_format		= IEEE
+IOFlexIO::out3D_extension	= ".vcamr"
+
+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)