AMR boundaries

Add support for using boundaries with AMR. Seems to work with symmetry
boundaries but gets a bit confused if you try refining outer
boundaries. Requires that you setup the domain using CoordBase.

darcs-hash:20050211135704-58c7f-6ed6ce1449de17599dd8b23abc909f5ba313d048.gz

2 files changed, 93 insertions, 5 deletions

diff --git a/CarpetDev/CarpetAdaptiveRegrid/par/AMR1.par b/CarpetDev/CarpetAdaptiveRegrid/par/AMR1.par
index 359c35d88..588d5de82 100644
--- a/CarpetDev/CarpetAdaptiveRegrid/par/AMR1.par
+++ b/CarpetDev/CarpetAdaptiveRegrid/par/AMR1.par
@@ -13,15 +13,25 @@ mol::initial_data_is_crap = "yes"
 ##carpet::veryverbose = "yes"
 
 grid::domain = "full"
-grid::type = "byspacing"
+grid::type = "coordbase"
 grid::avoid_origin = "no"
-driver::global_nx = 51
-driver::global_ny = 51
-driver::global_nz = 51
-grid::dxyz = 0.02
+
+coordbase::xmin = -0.5
+coordbase::ymin = -0.5
+coordbase::zmin = -0.5
+coordbase::xmax =  0.5
+coordbase::ymax =  0.5
+coordbase::zmax =  0.5
+
+coordbase::dx = 0.02
+coordbase::dy = 0.02
+coordbase::dz = 0.02
+
 driver::ghost_size = 1
 time::dtfac = 0.5
 
+carpet::domain_from_coordbase = yes
+
 carpet::max_refinement_levels = 2
 carpet::init_3_timelevels = yes
 carpetadaptiveregrid::max_error = 1.e-6
@@ -33,6 +43,8 @@ carpetadaptiveregrid::regrid_every = 1
 
 #carpetadaptiveregrid::verbose = yes
 #carpetadaptiveregrid::veryverbose = yes
+#carpet::verbose = yes
+#carpet::veryverbose = yes
 
 cactus::terminate = "time"
 cactus::cctk_final_time = 0.5
diff --git a/CarpetDev/CarpetAdaptiveRegrid/par/AMR2.par b/CarpetDev/CarpetAdaptiveRegrid/par/AMR2.par
new file mode 100644
index 000000000..e53527333
--- /dev/null
+++ b/CarpetDev/CarpetAdaptiveRegrid/par/AMR2.par
@@ -0,0 +1,76 @@
+
+ActiveThorns = "coordbase SymBase NaNChecker carpetReduce CartGrid3D carpet carpetlib carpetadaptiveregrid Boundary IOBasic IOUtil carpetIOASCII IDWaveMoL carpetSlab WaveMoL Time MoL reflectionsymmetry"
+
+
+IDWaveMoL::initial_data = "gaussian"
+
+wavemol::bound = "radiation"
+
+mol::initial_data_is_crap = "yes"
+
+##carpet::adaptive_stepsize = "yes"
+
+##carpet::veryverbose = "yes"
+
+grid::domain = "full"
+grid::type = "coordbase"
+grid::avoid_origin = "no"
+
+coordbase::xmin = -0.5
+coordbase::ymin = -0.5
+coordbase::zmin =  0.0
+coordbase::xmax =  0.5
+coordbase::ymax =  0.5
+coordbase::zmax =  0.5
+
+coordbase::dx = 0.02
+coordbase::dy = 0.02
+coordbase::dz = 0.02
+
+reflectionsymmetry::reflection_z = yes
+reflectionsymmetry::avoid_origin_z = no
+coordbase::boundary_size_z_lower = 2
+coordbase::boundary_shiftout_z_lower = 1
+
+carpet::prolongation_order_space = 3
+carpet::prolongation_order_time = 2
+driver::ghost_size = 2
+time::dtfac = 0.5
+
+carpet::domain_from_coordbase = yes
+
+carpet::max_refinement_levels = 2
+carpet::init_3_timelevels = yes
+carpetadaptiveregrid::max_error = 1.e-6
+carpetadaptiveregrid::pad = 2
+carpetadaptiveregrid::min_width = 6
+carpetadaptiveregrid::error_var = "mol::errorestimate[0]"
+#carpetadaptiveregrid::error_var = "wavemol::phi"
+carpetadaptiveregrid::regrid_every = 1
+
+#carpetadaptiveregrid::verbose = yes
+#carpetadaptiveregrid::veryverbose = yes
+#carpet::verbose = yes
+#carpet::veryverbose = yes
+
+cactus::terminate = "time"
+cactus::cctk_final_time = 0.5
+
+iobasic::outScalar_every = 1
+iobasic::outScalar_vars = "wavemol::phi mol::errorestimate"
+
+iobasic::outInfo_every = 1
+iobasic::outInfo_vars = "wavemol::phi mol::errorestimate[0]"
+
+ioascii::out1D_every = 1
+ioascii::out1D_vars = "wavemol::scalarevolvemol_scalar mol::errorestimate"
+
+IO::out_dir = $.
+
+#mol::ode_method = "Generic"
+#mol::mol_intermediate_steps = 4
+#mol::mol_num_scratch_levels = 3
+mol::ode_method = "RK45"
+mol::mol_intermediate_steps = 6
+mol::mol_num_scratch_levels = 6
+mol::adaptive_stepsize = "no"