User: knarf Date: 2011/05/18 09:07 AM

Removed: /examples/tov/ tov_pert.rpar

Modified: /examples/tov/ tov.rpar

Log: combine rpars

File Changes:

Directory: /examples/tov/ =========================

File [modified]: tov.rpar Delta lines: +3 -2 =================================================================== --- examples/tov/tov.rpar 2011-05-18 13:46:36 UTC (rev 87) +++ examples/tov/tov.rpar 2011-05-18 14:07:46 UTC (rev 88) @@ -8,7 +8,7 @@ my $K_TOV = 100; my $K_change = ""; if ($type eq "_pert") { - $K_TOV = 99; + $K_TOV = 90; $K_change = "TOVSolver::TOV_Transform_to_K = 100"; }

@@ -191,7 +191,8 @@ ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero"

-TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 +TOVSolver::Mass[0] = 1.4 TOVSolver::TOV_Gamma[0] = 2.0 TOVSolver::TOV_K[0] = $K_TOV $K_change

File [removed]: tov_pert.rpar Delta lines: +0 -242 =================================================================== --- examples/tov/tov_pert.rpar 2011-05-18 13:46:36 UTC (rev 87) +++ examples/tov/tov_pert.rpar 2011-05-18 14:07:46 UTC (rev 88) @@ -1,242 +0,0 @@ -#!/usr/bin/perl -w -use strict; - -foreach my $res ((8, 4, 2)) { - - my $out_every = 256/$res; - my $par = <<PAR_EOF; -## mclachlan tov_static -ActiveThorns = "Time MoL" -ActiveThorns = "Coordbase CartGrid3d Boundary StaticConformal" -ActiveThorns = "SymBase ADMBase TmunuBase HydroBase InitBase ADMCoupling ADMMacros" -ActiveThorns = "IOUtil" -#ActiveThorns = "Formaline" -ActiveThorns = "SpaceMask CoordGauge Constants LocalReduce aeilocalinterp LoopControl" -ActiveThorns = "Carpet CarpetLib CarpetReduce CarpetRegrid2 CarpetInterp" -ActiveThorns = "CarpetIOASCII CarpetIOScalar CarpetIOHDF5 CarpetIOBasic" -ActiveThorns = "SystemStatistics" - -ActiveThorns = "ADMConstraints" - -# grid parameters -CartGrid3D::type = "coordbase" -CartGrid3D::domain = "full" -CartGrid3D::avoid_origin = "no" -CoordBase::xmin = 0.0 -CoordBase::ymin = 0.0 -CoordBase::zmin = 0.0 -CoordBase::xmax = 240.0 -CoordBase::ymax = 240.0 -CoordBase::zmax = 240.0 -CoordBase::dx = $res -CoordBase::dy = $res -CoordBase::dz = $res - -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 -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 0 -CoordBase::boundary_shiftout_z_upper = 0 - - - - -ActiveThorns = "ReflectionSymmetry" - -ReflectionSymmetry::reflection_x = "yes" -ReflectionSymmetry::reflection_y = "yes" -ReflectionSymmetry::reflection_z = "yes" -ReflectionSymmetry::avoid_origin_x = "no" -ReflectionSymmetry::avoid_origin_y = "no" -ReflectionSymmetry::avoid_origin_z = "no" - - - -TmunuBase::stress_energy_storage = yes -TmunuBase::stress_energy_at_RHS = yes -TmunuBase::timelevels = 1 -TmunuBase::prolongation_type = none - -HydroBase::timelevels = 3 - -ADMMacros::spatial_order = 4 - -ADMBase::metric_type = "physical" - -ADMConstraints::bound = "static" -ADMConstraints::constraints_timelevels = 3 -ADMConstraints::constraints_persist = yes - -SpaceMask::use_mask = "yes" - -Cactus::terminate = "time" -Cactus::cctk_final_time = 1000 - -Carpet::domain_from_coordbase = "yes" -Carpet::enable_all_storage = no -Carpet::use_buffer_zones = "yes" - -Carpet::poison_new_timelevels = "yes" -Carpet::check_for_poison = "no" -Carpet::poison_value = 113 - -Carpet::init_3_timelevels = no -Carpet::init_fill_timelevels = "yes" - -CarpetLib::poison_new_memory = "yes" -CarpetLib::poison_value = 114 - -# system specific Carpet paramters -Carpet::max_refinement_levels = 10 -driver::ghost_size = 3 -Carpet::prolongation_order_space = 3 -Carpet::prolongation_order_time = 2 - -CarpetRegrid2::regrid_every = 0 -CarpetRegrid2::num_centres = 1 - -CarpetRegrid2::num_levels_1 = 5 -CarpetRegrid2::radius_1[1] =120.0 -CarpetRegrid2::radius_1[2] = 60.0 -CarpetRegrid2::radius_1[3] = 30.0 -CarpetRegrid2::radius_1[4] = 15.0 - - -time::dtfac = 0.25 -MoL::ODE_Method = "rk4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -# check all physical variables for NaNs -ActiveThorns = "NaNChecker" -#NaNChecker::check_every = 1 -#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" -#NaNChecker::check_vars = "ADMBase::metric ADMBase::lapse ADMBase::shift HydroBase::rho HydroBase::eps HydroBase::press HydroBase::vel" - -## Lapse Condition: \partial_t alpha = - alpha K -## Shift Condition: \partial_t beta^i = 0 - -# Hydro paramters -ActiveThorns = "EOS_Omni" -ActiveThorns = "GRHydro" - -HydroBase::evolution_method = "GRHydro" - -GRHydro::riemann_solver = "Marquina" -GRHydro::GRHydro_eos_type = "Polytype" -GRHydro::GRHydro_eos_table = "2D_Polytrope" -GRHydro::recon_method = "ppm" -GRHydro::GRHydro_stencil = 3 -GRHydro::bound = "none" -GRHydro::rho_abs_min = 1.e-10 -#GRHydro::GRHydro = 18 # Tmunu(10), rho,press,eps,w_lorentz,vel, tau -#GRHydro::GRHydro = 10 # gij(6), alpha, beta(3) - - -ActiveThorns = "GenericFD NewRad" -ActiveThorns = "ML_BSSN ML_BSSN_Helper" -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" - -TmunuBase::support_old_CalcTmunu_mechanism = "no" - -ML_BSSN::timelevels = 3 - -ML_BSSN::harmonicN = 1 # 1+log -ML_BSSN::harmonicF = 1.0 # 1+log -ML_BSSN::LapseACoeff = 1.0 -ML_BSSN::ShiftBCoeff = 1.0 -ML_BSSN::ShiftGammaCoeff = 0.0 -ML_BSSN::AlphaDriver = 0.0 -ML_BSSN::BetaDriver = 0.0 -ML_BSSN::LapseAdvectionCoeff = 0.0 -ML_BSSN::ShiftAdvectionCoeff = 0.0 - -ML_BSSN::MinimumLapse = 1.0e-8 - -ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" -ML_BSSN::my_rhs_boundary_condition = "NewRad" - -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" - - - -# init parameters -InitBase::initial_data_setup_method = "init_some_levels" - -ActiveThorns = "TOVSolver" - -ADMBase::initial_data = "tov" -ADMBase::initial_lapse = "tov" -ADMBase::initial_shift = "tov" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -#TOVSolver::TOV_Rho_Central[0] = 1.28e-3 -TOVSolver::TOV_Mass[0] = 1.4 -TOVSolver::TOV_Gamma[0] = 2.0 -TOVSolver::TOV_K[0] = 100.0 - -# I/O - -IO::out_dir = $parfile - -IOBasic::outInfo_every = $out_every -IOBasic::outInfo_vars = "HydroBase::rho ADMBase::lapse" - -IOScalar::outScalar_every = $out_every -IOScalar::one_file_per_group = yes -IOScalar::outScalar_vars = " - HydroBase::rho - HydroBase::press - HydroBase::eps - HydroBase::vel - ADMBase::lapse - ADMBase::metric - ADMBase::curv - ADMConstraints::ham - ADMConstraints::momentum - SystemStatistics::process_memory -" - -IOASCII::out1D_every = $out_every -IOASCII::one_file_per_group = yes -IOASCII::output_symmetry_points = no -IOASCII::out3D_ghosts = no -IOASCII::out3D_outer_ghosts = no -IOASCII::out1D_vars = " - HydroBase::rho - HydroBase::press - HydroBase::eps - HydroBase::vel - ADMBase::lapse - ADMBase::metric - ADMBase::curv - ADMConstraints::ham - ADMConstraints::momentum -" -PAR_EOF - - open(OUT, ">tov_${res}.par"); - print OUT "$par"; - close(OUT); -} -