Index: interface.ccl =================================================================== --- interface.ccl (revision 13) +++ interface.ccl (working copy) @@ -2,4 +2,4 @@ # $Header$ implements: idconstraintviolate -inherits: admbase, grid \ No newline at end of file +inherits: admbase, grid Index: param.ccl =================================================================== --- param.ccl (revision 13) +++ param.ccl (working copy) @@ -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" Index: par/Gaussian.par =================================================================== --- par/Gaussian.par (revision 13) +++ par/Gaussian.par (working copy) @@ -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