Added a routine for ensuring that the grid set up in the parameter file

is cactus compatible. If the new parameter "allow_grid_resize" is set,
then cartoon will now actually modify the grid that was specified by the
user. Though not really a nice solution, it is potentially less annoying
than having the user try to figure out appropriate grid dimensions
themselves.
The old test suite for cartoon (test_cartoon.par) was not a very good
example, as it actually set up an less-than-optimal grid (with bitant
specified in the wrong plane), and could only run on 2-processors. So
I've removed it and replaced it with a pair of new test suites which
demonstrate cartoon in "full" and "bitant" mode using the new grid
resizing mechanism.
Documentation as such does not exist yet, but will in the fullness of
time. For the time being, please refer to the sample parameter files
in the test directory to see how to set things up.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Cartoon2D/trunk@26 eec4d7dc-71c2-46d6-addf-10296150bf52

1 files changed, 90 insertions, 0 deletions

diff --git a/test/test_cartoon_1.par b/test/test_cartoon_1.par
new file mode 100644
index 0000000..08895e4
--- /dev/null
+++ b/test/test_cartoon_1.par
@@ -0,0 +1,90 @@
+#==============================================================================
+# Cartoon2D, using a bitant grid
+#==============================================================================
+
+ActiveThorns = "cartoon2d idanalyticbh einstein adm_bssn admconstraints ahfinder pughinterp boundary time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic"
+
+#------------------------------------------------------------------------------
+# Run parameters for Cartoon
+#------------------------------------------------------------------------------
+
+# Cartoon
+cartoon2d::cartoon_active		= "yes"
+cartoon2d::verbose              	= "yes"
+cartoon2d::stencil              	= "yes"
+cartoon2d::order                	= 3
+cartoon2d::allow_grid_resize		= "yes"
+
+# This is for 2 processors, you have to
+# change manually by processor numbers.
+
+#driver::processor_topology       	= "manual"
+#driver::processor_topology_3d_x  	= 1
+#driver::processor_topology_3d_y  	= 1
+#driver::processor_topology_3d_z  	= 2
+
+# Grid
+
+driver::global_nx                       = 16
+driver::global_ny                       = 3
+driver::global_nz                       = 16
+
+driver::ghost_size_x                    = 2
+driver::ghost_size_y                    = 1
+driver::ghost_size_z                    = 2
+
+grid::avoid_originy                     = "no"
+
+grid::type                              = "byspacing"
+grid::domain				= "bitant"
+grid::bitant_plane			= "xy"
+grid::dxyz                              = 0.2
+
+time::dtfac                             = 0.25
+
+cactus::cctk_itlast                     = 1
+
+# MODEL 
+
+einstein::initial_data                  = "schwarzschild"
+
+#--- Evolution
+einstein::evolution_system 		= "adm_bssn"
+adm_bssn::method 			= "leapfrog"
+adm_bssn::bound  			= "flat"
+adm_bssn::cartoon                	= "yes"
+admconstraints::cartoon			= "yes"
+
+#--- gauge
+einstein::slicing                       = "1+log"
+
+#--- AHFinder.
+
+ahfinder::ahf_active   			= "yes"
+ahfinder::ahf_cartoon			= "yes"
+ahfinder::ahf_gaussout 			= "yes"
+
+ahfinder::ahf_phi    			= "no"
+ahfinder::ahf_octant 			= "high"
+
+ahfinder::ahf_lmax   			= 0
+
+# Interpolation.
+
+ahfinder::interpolation_order 		= 1
+
+#------------------------------------------------------------------------------
+# Output
+#------------------------------------------------------------------------------
+IO::outdir                              = "test_cartoon_1"
+IO::out_fileinfo   			= "none"
+
+IOBasic::outScalar_every                = 1
+IOBasic::outScalar_vars                 = "admconstraints::hamiltonian admconstraints::momentum einstein::alp einstein::grr einstein::metric einstein::trace_of_K einstein::betax einstein::betay einstein::betaz adm_bssn::adm_bs_gxx adm_bssn::adm_bs_phi adm_bssn::ADM_BSSN_K adm_bssn::ADM_BSSN_gamma"
+
+IOASCII::out1D_every                    = 1
+IOASCII::out1D_vars                     = "admconstraints::hamiltonian admconstraints::momentum einstein::alp einstein::grr einstein::metric einstein::trace_of_K einstein::betax einstein::betay einstein::betaz adm_bssn::adm_bs_gxx adm_bssn::adm_bs_phi adm_bssn::ADM_BSSN_K adm_bssn::ADM_BSSN_gamma"
+
+IOBasic::outInfo_every                  = 1
+IOBasic::outInfo_vars                   = "einstein::alp adm_bssn::adm_bs_dtalp"
+#==============================================================================