[Commits] [svn:einsteintoolkit] Paper_EinsteinToolkit_2010/examples/tov/ (Rev. 88)
knarf at cct.lsu.edu
knarf at cct.lsu.edu
Wed May 18 09:07:46 CDT 2011
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);
-}
-
More information about the Commits
mailing list