[Commits] [svn:einsteintoolkit] Paper_EinsteinToolkit_2010/ (Rev. 87)

knarf at cct.lsu.edu knarf at cct.lsu.edu
Wed May 18 08:46:36 CDT 2011


User: knarf
Date: 2011/05/18 08:46 AM

Added:
 /examples/tov/
  tov_pert.rpar

Log:
 parameter file for perturbed TOV star - note, this so far uses uncommited code within the TOVSolver

Directory Changes:

Directory: /svn:executable/
===========================

   + *

File Changes:

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

File [added]: tov_pert.rpar
Delta lines: +242 -0
===================================================================
--- examples/tov/tov_pert.rpar	                        (rev 0)
+++ examples/tov/tov_pert.rpar	2011-05-18 13:46:36 UTC (rev 87)
@@ -0,0 +1,242 @@
+#!/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);
+}
+



Property changes on: examples/tov/tov_pert.rpar
___________________________________________________________________



More information about the Commits mailing list