User: eschnett
Date: 2010/06/30 08:07 PM
Modified:
/trunk/
interface.ccl, param.ccl
/trunk/par/
Gaussian.par
Log:
Update the example parameter file in thorn IDConstraintViolate to
using McLachlan. (Note: the test case should also be updated.)
File Changes:
Directory: /trunk/par/
======================
File [modified]: Gaussian.par
Delta lines: +239 -23
===================================================================
--- trunk/par/Gaussian.par 2010-06-17 19:59:27 UTC (rev 13)
+++ trunk/par/Gaussian.par 2010-07-01 01:07:10 UTC (rev 14)
@@ -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
Directory: /trunk/
==================
File [modified]: interface.ccl
Delta lines: +1 -1
===================================================================
--- trunk/interface.ccl 2010-06-17 19:59:27 UTC (rev 13)
+++ trunk/interface.ccl 2010-07-01 01:07:10 UTC (rev 14)
@@ -2,4 +2,4 @@
# $Header$
implements: idconstraintviolate
-inherits: admbase, grid
\ No newline at end of file
+inherits: admbase, grid
File [modified]: param.ccl
Delta lines: +5 -5
===================================================================
--- trunk/param.ccl 2010-06-17 19:59:27 UTC (rev 13)
+++ trunk/param.ccl 2010-07-01 01:07:10 UTC (rev 14)
@@ -7,19 +7,19 @@
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"
