Update the example parameter file in thorn IDConstraintViolate to

using McLachlan.  (Note: the test case should also be updated.)


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDConstraintViolate/trunk@14 9ea3c693-0e4f-0410-9b72-c04bdc00281e

3 files changed, 245 insertions, 29 deletions

diff --git a/interface.ccl b/interface.ccl
index 5893745..e5c5d58 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -2,4 +2,4 @@
 # $Header$
 
 implements: idconstraintviolate
-inherits: admbase, grid
\ No newline at end of file
+inherits: admbase, grid
diff --git a/par/Gaussian.par b/par/Gaussian.par
index d4feba5..3df4924 100644
--- a/par/Gaussian.par
+++ b/par/Gaussian.par
@@ -1,35 +1,251 @@
+Cactus::cctk_run_title = "Constraint violating Gaussian initial data"
 
-!DESC "Constraint violating Gaussian initial data"
+Cactus::cctk_full_warnings         = yes
+Cactus::highlight_warning_messages = no
 
-ActiveThorns = "SpaceMask Boundary ADMBase ADM ADMMacros ADMAnalysis ADMConstraints CoordGauge StaticConformal CartGrid3D IOBasic IOUtil IOASCII IDConstraintViolate PsiKadelia Pugh PughReduce PughSlab Time"
- 
-#------------------------------------------------------------------------------
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 10.0
 
-driver::global_nsize       = 40
 
-cactus::cctk_itlast        = 1000
 
-IOBasic::outInfo_vars      = "admbase::alp admbase::gxx admconstraints::ham"
-IOBasic::outInfo_every     = 5
-IOASCII::out1D_vars        = "admbase::alp admbase::gxx admanalysis::grr psikadelia::psi0re admconstraints::ham admconstraints::momentum"
-IOASCII::out1D_every       =  5
+ActiveThorns = "IOUtil"
 
-#------------------------------------------------------------------------------
+IO::out_dir = $parfile
 
-IO::out_dir                      = "Gaussian"
 
-grid::type                       = "box"
 
-time::dtfac                      = 0.25 
- 
-admbase::metric_type             = "physical"
-admbase::initial_data            = "constraint violating gaussian"  
-admbase::evolution_method        = "ADM"     
-admbase::lapse_evolution_method  = "1+log"
-admbase::initial_lapse	         = "one"
+ActiveThorns = "Fortran"
 
-adm::bound                       = "radiative"      
-adm::method                      = "stagleap"
+ActiveThorns = "GSL"
 
-#------------------------------------------------------------------------------
+ActiveThorns = "GenericFD"
 
+ActiveThorns = "LoopControl"
+
+
+
+ActiveThorns = "InitBase"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetReduce"
+
+Carpet::domain_from_coordbase = yes
+
+driver::ghost_size = 3
+
+Carpet::init_fill_timelevels = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 512
+NaNChecker::action_if_found = "just warn"
+NaNChecker::check_vars      = "
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+"
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase SymBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmin = -1.0
+CoordBase::ymin = -1.0
+CoordBase::zmin = -1.0
+CoordBase::xmax = +1.0
+CoordBase::ymax = +1.0
+CoordBase::zmax = +1.0
+CoordBase::dx   =  0.05
+CoordBase::dy   =  0.05
+CoordBase::dz   =  0.05
+
+CoordBase::boundary_size_x_lower = 3
+CoordBase::boundary_size_y_lower = 3
+CoordBase::boundary_size_z_lower = 3
+CoordBase::boundary_size_x_upper = 3
+CoordBase::boundary_size_y_upper = 3
+CoordBase::boundary_size_z_upper = 3
+
+GenericFD::boundary_width = 3
+
+CartGrid3D::type = "coordbase"
+
+
+
+ActiveThorns = "SphericalSurface"
+
+
+
+ActiveThorns = "MoL Time"
+
+MoL::ODE_Method             = "RK4"
+MoL::MoL_Intermediate_Steps = 4
+MoL::MoL_Num_Scratch_Levels = 1
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"
+
+
+
+ActiveThorns = "IDConstraintViolate"
+
+ADMBase::initial_data    = "constraint violating gaussian"
+ADMBase::initial_lapse   = "one"
+ADMBase::initial_shift   = "zero"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
+
+
+
+ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad"
+
+ADMBase::evolution_method         = "ML_BSSN"
+ADMBase::lapse_evolution_method   = "ML_BSSN"
+ADMBase::shift_evolution_method   = "ML_BSSN"
+ADMBase::dtlapse_evolution_method = "ML_BSSN"
+ADMBase::dtshift_evolution_method = "ML_BSSN"
+
+ML_BSSN::timelevels = 3
+
+ML_BSSN::harmonicF           = 2.0    # 1+log
+ML_BSSN::harmonicN           = 1      # 1+log
+ML_BSSN::ShiftGammaCoeff     = 0.75
+ML_BSSN::LapseACoeff         = 1.0
+ML_BSSN::ShiftBCoeff         = 1.0
+ML_BSSN::AlphaDriver         = 1.0
+ML_BSSN::BetaDriver          = 1.0
+ML_BSSN::LapseAdvectionCoeff = 0.0
+ML_BSSN::ShiftAdvectionCoeff = 0.0
+
+ML_BSSN::MinimumLapse        = 1.0e-8
+
+ML_BSSN::dt_lapse_shift_method = "noLapseShiftAdvection"
+
+ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
+ML_BSSN::my_rhs_boundary_condition     = "NewRad"
+Boundary::radpower                     = 2
+
+ML_BSSN::ML_log_confac_bound = "none"
+ML_BSSN::ML_metric_bound     = "none"
+ML_BSSN::ML_Gamma_bound      = "none"
+ML_BSSN::ML_trace_curv_bound = "none"
+ML_BSSN::ML_curv_bound       = "none"
+ML_BSSN::ML_lapse_bound      = "none"
+ML_BSSN::ML_dtlapse_bound    = "none"
+ML_BSSN::ML_shift_bound      = "none"
+ML_BSSN::ML_dtshift_bound    = "none"
+
+
+
+ActiveThorns = "Dissipation"
+
+Dissipation::order = 5
+Dissipation::vars  = "
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
+"
+
+
+
+ActiveThorns = "ML_ADMConstraints"
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        Carpet::physical_time_per_hour
+        ML_ADMConstraints::H
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::one_file_per_group = yes
+
+IOScalar::outScalar_every = 1
+IOScalar::outScalar_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ML_ADMConstraints::ML_Ham
+        ML_ADMConstraints::ML_mom
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::output_symmetry_points = no
+IOASCII::out3D_ghosts           = no
+
+IOASCII::out0D_every = 1
+IOASCII::out0D_vars  = "
+        Carpet::timing
+"
+
+IOASCII::out1D_every = 10
+IOASCII::out1D_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ML_ADMConstraints::ML_Ham
+        ML_ADMConstraints::ML_mom
+"
+
+
+
+ActiveThorns = "Formaline"
+
+
+
+ActiveThorns = "TimerReport"
+
+TimerReport::out_every                  = 10
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers_together = yes
+TimerReport::output_all_timers_readable = yes
+TimerReport::n_top_timers               = 20
diff --git a/param.ccl b/param.ccl
index 9b8bde3..9a7fa73 100644
--- a/param.ccl
+++ b/param.ccl
@@ -7,19 +7,19 @@ uses KEYWORD metric_type
 
 extends KEYWORD initial_data
 {
-  "constraint violating gaussian" :: "Constraint violating gaussian"
+  "constraint violating gaussian" :: "Constraint violating Gaussian"
 }
 
 private:
 
-REAL radius "The radial position of the gaussian waves"
+REAL radius "The radial position of the Gaussian wave"
 {
- 0:* :: "Positive"
+  0:* :: "Positive"
 } 0.0
 
-REAL sigma "The sigma for the gaussian wave"
+REAL sigma "The sigma for the Gaussian wave"
 {
- 0:* :: "Positive"
+  (0:* :: "Positive"
 } 0.1
 
 REAL amplitude "The amplitude of the waves"