Hello everyone,

I am studying how to use Einstein_toolkit, I want to try to add atmosphere to a BBH evolution.

I add the below things to BBHLowRes.par in the par dir.

However, it reports "nan" error soon after the beginning of evolution.

I do not know how to correct it.

Could anyone give me some suggestion?

ActiveThorns = "GRHydro"

ActiveThorns = "GRHydro_InitData"

ActiveThorns = "HydroBase"

ActiveThorns = "EOS_Omni"

HydroBase::initial_hydro = "only_atmo"

HydroBase::evolution_method = "GRHydro"

HydroBase::timelevels = 3

SpaceMask::use_mask = "yes"

GRHydro::GRHydro_countmax = 100

GRHydro::riemann_solver = "HLLE" # Marquina is currently not supported by MP

GRHydro::recon_method = "ppm"

GRHydro::GRHydro_stencil = 3

GRHydro::bound = "flat"

GRHydro::rho_abs_min = 1.0e-6

GRHydro::GRHydro_atmo_tolerance = 0.

GRHydro::c2p_reset_pressure = "yes"

GRHydro::GRHydro_eos_type = "General"

GRHydro::GRHydro_eos_table = "Ideal_Fluid"

GRHydro::GRHydro_MaxNumSandRVars = 0

GRHydro::use_enhanced_ppm = "yes"

GRHydro::sync_conserved_only = "yes"

GRHydro::reconstruct_Wv = "yes"

GRHydro::c2p_resort_to_bisection = "yes"

GRHydro::use_cxx_code = "yes"

TmunuBase::stress_energy_storage = "yes"

TmunuBase::stress_energy_at_RHS = "yes"

TmunuBase::timelevels = 1

TmunuBase::prolongation_type = "none"

TmunuBase::support_old_CalcTmunu_mechanism = "no"

EOS_Omni::gl_gamma = 2

EOS_Omni::gl_k=100

EOS_Omni::poly_k = 100

Yours,

Zhi-Chao, Zhao

Institute of High Energy Physics, CAS, China.