GRHydro_InitData: adapt rotor parfile to current code

From: Roland Haas <roland.haas@physics.gatech.edu>

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@141 ac85fae7-cede-4708-beff-ae01c7fa1c26

1 files changed, 60 insertions, 27 deletions

diff --git a/par/rotor.par b/par/rotor.par
index 476bc08..d2caef3 100644
--- a/par/rotor.par
+++ b/par/rotor.par
@@ -1,15 +1,16 @@
 ##
-ActiveThorns = "time MoL cartgrid3d carpetioascii ioutil Fortran boundary  hydrobase grhydro grhydro_initdata nanchecker eos_omni CoordBase  LoopControl CarpetLib Carpet CarpetReduce spacemask ADMCoupling ADMMacros ADMBASE STATICCONFORMAL TMUNUBASE SymBase localreduce CarpetIOBasic initbase "
+ActiveThorns = "time MoL cartgrid3d carpetioascii ioutil boundary  hydrobase grhydro grhydro_initdata nanchecker eos_omni CoordBase  LoopControl CarpetLib Carpet CarpetReduce spacemask ADMCoupling ADMMacros ADMBASE STATICCONFORMAL TMUNUBASE SymBase localreduce CarpetIOBasic initbase "
 ActiveThorns = "CarpetIOHDF5"
-#This is a hack. ghost_size should be the same as grhydro::grhydro_stencil.
 
-driver::ghost_size=3
-grhydro::grhydro_stencil=3
+driver::ghost_size=2       # for TVD
+grhydro::grhydro_stencil=2
 
-time::dtfac = 0.25
+Time::dtfac = 0.25         # 0.5 see p 409 (0.25 from David Neilsen, Eric WHirschmann, and R Steven Millward arXiv:gr-qc/0512147v1)
 
-methodoflines::ODE_Method = "rk2"
-methodoflines::MoL_Intermediate_Steps=2
+# rk2 dies quite quickly
+methodoflines::ODE_Method = "rk3"
+methodoflines::MoL_Intermediate_Steps=3
+#methodoflines::MoL_Num_Scratch_Levels=1
 
 ADMBase::initial_shift = "zero"
 
@@ -19,11 +20,16 @@ HydroBase::timelevels=2
 HydroBase::evolution_method="grhydro"
 HydroBase::Bvec_evolution_method="grhydro"
 
+# rotor problem from DelZanna
+# needs to be before EOS_Omni
+grhydro_initdata::rotor_use_smoothing = no
+
 SpaceMask::use_mask="yes"
-EOS_Omni::gl_gamma=2.0
+EOS_Omni::gl_gamma=5./3. # relativistic rotor problem of DelZanna & Buccianti & Londrille AaA_400_397_2003
+# try to mimic the core fluid
+EOS_Omni::poly_k = "grhydro_initdata::rotor_pressin/grhydro_initdata::rotor_rhoin**EOS_Omni::gl_gamma"
+EOS_Omni::poly_gamma="EOS_Omni::gl_gamma"
 grhydro::recon_method="tvd"
-grhydro::method_type="RSA FV"
-grhydro::recon_vars="primitive"
 grhydro::riemann_solver="HLLE"
 grhydro::tvd_limiter="vanleerMC"
 grhydro::bound="flat"
@@ -31,29 +37,60 @@ grhydro::grhydro_eos_type = "General"
 grhydro::grhydro_eos_table="Ideal_Fluid"
 grhydro::Grhydro_MaxNumConstrainedVars=13
 grhydro::GRHydro_MaxNumEvolvedVars=10
+grhydro::GRHydro_MaxNumSandRVars=0 # ADMBase is not evolved
+grhydro::rho_abs_min = 1e-14 # need some reasonable value here, even though no density ever gets this low. This is about the same as GRHydro's default value from rho_central and rho_rel_min
 
 grhydro::clean_divergence = "yes"
-grhydro::kap_dc = 1.0
+grhydro::transport_constraints = "no"
+grhydro::decouple_normal_Bfield = "yes"
+grhydro::kap_dc = 5.0
 grhydro::track_divB = "yes"
 
-grid::type = "BySpacing"
+grid::type = "coordbase"
 grid::domain = "full"
-grid::dxyz = 0.004
 
-driver::global_nx = 250
-driver::global_ny = 250
-driver::global_nz = 10
+CoordBase::domainsize = "minmax"
+CoordBase::xmin = 0.
+CoordBase::ymin = 0.
+CoordBase::zmin = -0.0125
+CoordBase::xmax = 1.0
+CoordBase::ymax = 1.0
+CoordBase::zmax = 0.0125
+
+CoordBase::boundary_size_x_lower = carpet::ghost_size
+CoordBase::boundary_size_y_lower = carpet::ghost_size
+CoordBase::boundary_size_z_lower = carpet::ghost_size
+CoordBase::boundary_size_x_upper = carpet::ghost_size
+CoordBase::boundary_size_y_upper = carpet::ghost_size
+CoordBase::boundary_size_z_upper = carpet::ghost_size
+
+CoordBase::spacing = "numcells"
+CoordBase::ncells_x = 400
+CoordBase::ncells_y = 400
+CoordBase::ncells_z = 10
+
+carpet::domain_from_coordbase = "yes"
 
 Cactus::terminate       = "time"
-Cactus::cctk_final_time = 0.4
-#Cactus::terminate = "iteration"
-#cactus::cctk_itlast = 2
+Cactus::cctk_final_time = 0.4 # see DelZanna p 410
 
 IO::out_dir = $parfile
-CarpetIOBasic::outInfo_vars="hydrobase::rho Hydrobase::bvec GRHydro::bcons"
+CarpetIOBasic::outInfo_vars="
+  hydrobase::rho
+  hydrobase::press
+  hydrobase::eps
+  hydrobase::vel[0] hydrobase::vel[1]
+  hydrobase::w_lorentz
+  GRHydro::bcons[0] GRHydro::bcons[1]
+  GRHydro::divB GRHydro::psidc
+"
 CarpetIOBasic::outInfo_every=1
-CarpetIOASCII::out1D_every = 20
-CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps  HydroBase::vel grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec grhydro::Bcons"
+CarpetIOASCII::out1D_every = 10
+CarpetIOASCII::out1D_z    = "no"
+CarpetIOASCII::out1D_d    = "no"
+CarpetIOASCII::out1D_xline_y = 0.5
+CarpetIOASCII::out1D_yline_x = 0.5
+CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps  HydroBase::vel grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec grhydro::Bcons grhydro::divB grhydro::psidc"
 
 CarpetIOHDF5::one_file_per_group             = "yes"
 CarpetIOHDF5::open_one_input_file_at_a_time  = "yes"
@@ -63,9 +100,6 @@ CarpetIOHDF5::out2D_xz                       = "no"
 CarpetIOHDF5::out2D_yz                       = "no"
 CarpetIOHDF5::out2D_xyplane_z                = 0.0
 CarpetIOHDF5::out2D_vars      = "
-#  CarpetReduce::weight
-#  GRHydro::GRHydro_tracers
-#  Grid::coordinates
   GRHydro::Bcons
   GrHydro::dens 
   GrHydro::divB
@@ -77,6 +111,5 @@ CarpetIOHDF5::out2D_vars      = "
   HydroBase::rho
   HydroBase::vel
   HydroBase::w_lorentz
+  GRHydro::psidc
 "
-
-