User: rhaas Date: 2012/08/09 01:26 AM

Modified: /trunk/par/ rotor.par

Log: GRHydro_InitData: adapt rotor parfile to current code

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

File Changes:

Directory: /trunk/par/ ======================

File [modified]: rotor.par Delta lines: +60 -27 =================================================================== --- trunk/par/rotor.par 2012-08-09 06:26:40 UTC (rev 140) +++ trunk/par/rotor.par 2012-08-09 06:26:42 UTC (rev 141) @@ -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::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_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_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 @@ HydroBase::rho HydroBase::vel HydroBase::w_lorentz + GRHydro::psidc " - -