Found a set of parameter files for planar waves with periodic

Found a set of parameter files for planar waves with periodic
boundaries.  Was more complicated than I suspected.  Embarrassing.

darcs-hash:20010316214231-f6438-34a413431040d7a968ca0f97d5ec10875168749f.gz

6 files changed, 88 insertions, 66 deletions

diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par
index 7954e1fe0..e5ff3f880 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par
@@ -7,22 +7,26 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
 Cactus::cctk_itlast		= 120
- 
-Time::dtfac			= 0.5
 
-driver::global_nx		= 33
-driver::global_ny		= 45
-driver::global_nz		= 61
+Time::dtfac			= 0.5
 
-grid::type			= byspacing
-grid::dxyz			= 0.3
+driver::global_nx		= 32
+driver::global_ny		= 46
+driver::global_nz		= 62
+
+grid::type			= byrange
+grid::xmin			= -4.8
+grid::ymin			= -6.6
+grid::zmin			= -9.0
+grid::xmax			=  4.5
+grid::ymax			=  6.9
+grid::zmax			=  9.3
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic"
 
@@ -38,6 +42,6 @@ IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.21505376344086	# 2 / ((33-2) * 0.3)
-IDScalarWave::ky		= 0.155038759689922	# 2 / ((45-2) * 0.3)
-IDScalarWave::kz		= 0.225988700564972	# 4 / ((61-2) * 0.3)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par
index 86e6b91c6..a639af040 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par
@@ -7,22 +7,26 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
 Cactus::cctk_itlast		= 60
- 
+
 Time::dtfac			= 0.5
 
 driver::global_nx		= 17
-driver::global_ny		= 23
-driver::global_nz		= 31
-
-grid::type			= byspacing
-grid::dxyz			= 0.6
+driver::global_ny		= 24
+driver::global_nz		= 32
+
+grid::type			= byrange
+grid::xmin			= -4.8
+grid::ymin			= -6.6
+grid::zmin			= -9.0
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic_coarse"
 
@@ -32,12 +36,12 @@ IOBasic::outinfo_every		= 1 # 10
 #IOBasic::outScalar_every	= 2
 #IOBasic::outScalar_vars		= "wavetoy::phi"
 
-IOASCII::out1D_every		= 2 # 1 # 2
+IOASCII::out1D_every		= 1 # 2
 IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.222222222222222	# 2 / ((17-2) * 0.6)
-IDScalarWave::ky		= 0.158730158730159	# 2 / ((23-2) * 0.6)
-IDScalarWave::kz		= 0.229885057471264	# 4 / ((31-2) * 0.6)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par
index e16bc1bbc..871f26f6f 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par
@@ -7,24 +7,28 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl2.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
-Cactus::cctk_itlast		= 60
+Cactus::cctk_itlast		= 120
  
 Time::dtfac			= 0.5
 
 driver::global_nx		= 17
-driver::global_ny		= 23
-driver::global_nz		= 31
+driver::global_ny		= 24
+driver::global_nz		= 32
 
 Carpet::max_refinement_levels	= 2
 
-grid::type			= byspacing
-grid::dxyz			= 0.6
+grid::type			= byrange
+grid::xmin			= -4.8
+grid::ymin			= -6.6
+grid::zmin			= -9.0
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic_coarse_rl2"
 
@@ -40,6 +44,6 @@ IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.222222222222222	# 2 / ((17-2) * 0.6)
-IDScalarWave::ky		= 0.158730158730159	# 2 / ((23-2) * 0.6)
-IDScalarWave::kz		= 0.229885057471264	# 4 / ((31-2) * 0.6)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par
index 6167a146d..f7464f676 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par
@@ -7,24 +7,28 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_coarse_rl3.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
-Cactus::cctk_itlast		= 60
+Cactus::cctk_itlast		= 240
  
 Time::dtfac			= 0.5
 
 driver::global_nx		= 17
-driver::global_ny		= 23
-driver::global_nz		= 31
+driver::global_ny		= 24
+driver::global_nz		= 32
 
 Carpet::max_refinement_levels	= 3
 
-grid::type			= byspacing
-grid::dxyz			= 0.6
+grid::type			= byrange
+grid::xmin			= -4.8
+grid::ymin			= -6.6
+grid::zmin			= -9.0
+grid::xmax			=  4.8
+grid::ymax			=  7.2
+grid::zmax			=  9.6
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic_coarse_rl3"
 
@@ -40,6 +44,6 @@ IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.222222222222222	# 2 / ((17-2) * 0.6)
-IDScalarWave::ky		= 0.158730158730159	# 2 / ((23-2) * 0.6)
-IDScalarWave::kz		= 0.229885057471264	# 4 / ((31-2) * 0.6)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par
index f5ab34e5d..b33487e77 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par
@@ -7,24 +7,26 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl2.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
-Cactus::cctk_itlast		= 120
+Cactus::cctk_itlast		= 240
  
 Time::dtfac			= 0.5
 
-driver::global_nx		= 33
-driver::global_ny		= 45
-driver::global_nz		= 61
+driver::global_nx		= 32
+driver::global_ny		= 46
+driver::global_nz		= 62
 
 Carpet::max_refinement_levels	= 2
 
-grid::type			= byspacing
-grid::dxyz			= 0.3
+grid::type			= byrange
+grid::xyzmin			= -0.3
+grid::xmax			=  9.0
+grid::ymax			= 13.2
+grid::zmax			= 18.0
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic_rl2"
 
@@ -40,6 +42,6 @@ IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.21505376344086	# 2 / ((33-2) * 0.3)
-IDScalarWave::ky		= 0.155038759689922	# 2 / ((45-2) * 0.3)
-IDScalarWave::kz		= 0.225988700564972	# 4 / ((61-2) * 0.3)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)
diff --git a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par
index 9a089ab13..ff5560cec 100644
--- a/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par
+++ b/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par
@@ -7,24 +7,28 @@
 #   @enddesc
 # @@*/
 #
-# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par,v 1.2 2001/03/16 21:32:18 eschnett Exp $
+# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/WaveToyF77/par/wavetoyf77_periodic_rl3.par,v 1.3 2001/03/16 22:42:31 eschnett Exp $
 
 ActiveThorns = "Boundary IOBasic IOUtil Time Cart3d Carpet CarpetIOASCII CarpetLib CarpetSlab IDScalarWave WaveToyF77"
 
-Cactus::cctk_itlast		= 120
+Cactus::cctk_itlast		= 480
  
 Time::dtfac			= 0.5
 
-driver::global_nx		= 33
-driver::global_ny		= 45
-driver::global_nz		= 61
+driver::global_nx		= 32
+driver::global_ny		= 46
+driver::global_nz		= 62
 
 Carpet::max_refinement_levels	= 3
 
-grid::type			= byspacing
-grid::dxyz			= 0.3
+grid::type			= byrange
+grid::xmin			= -4.8
+grid::ymin			= -6.6
+grid::zmin			= -9.0
+grid::xmax			=  4.5
+grid::ymax			=  6.9
+grid::zmax			=  9.3
 grid::mode			= periodic
-grid::avoid_origin		= no
 
 IO::outdir			= "wavetoyf77_periodic_rl3"
 
@@ -40,6 +44,6 @@ IOASCII::out1D_vars		= "wavetoy::phi grid::coordinates"
 WaveToyF77::bound		= radiation
 
 IDScalarWave::initial_data	= plane
-IDScalarWave::kx		= 0.21505376344086	# 2 / ((33-2) * 0.3)
-IDScalarWave::ky		= 0.155038759689922	# 2 / ((45-2) * 0.3)
-IDScalarWave::kz		= 0.225988700564972	# 4 / ((61-2) * 0.3)
+IDScalarWave::kx		= 0.222222222222222	# 2 / ((32-2) * 0.3)
+IDScalarWave::ky		= 0.151515151515152	# 2 / ((46-2) * 0.3)
+IDScalarWave::kz		= 0.222222222222222	# 4 / ((62-2) * 0.3)