<div dir="ltr"><div>Hello everyone, </div><div><br></div><div>I am studying how to use Einstein_toolkit, I want to try to add atmosphere to a BBH evolution.</div><div>I add the below things to BBHLowRes.par in the par dir.</div><div>However, it reports &quot;nan&quot; error soon after the beginning of evolution.</div><div>I do not know how to correct it. </div><div>Could anyone give me some suggestion?</div><div><br></div><div>&quot;<br></div><div><div>ActiveThorns = &quot;GRHydro&quot;</div><div>ActiveThorns = &quot;GRHydro_InitData&quot;</div><div>ActiveThorns = &quot;HydroBase&quot;</div><div>ActiveThorns = &quot;EOS_Omni&quot;</div><div><br></div><div>HydroBase::initial_hydro         = &quot;only_atmo&quot;</div><div>HydroBase::evolution_method      = &quot;GRHydro&quot;</div><div>HydroBase::timelevels            = 3</div><div><br></div><div>SpaceMask::use_mask = &quot;yes&quot;</div><div>GRHydro::GRHydro_countmax       = 100</div><div>GRHydro::riemann_solver            = &quot;HLLE&quot;   # Marquina is currently not supported by MP</div><div>GRHydro::recon_method              = &quot;ppm&quot;</div><div>GRHydro::GRHydro_stencil            = 3</div><div>GRHydro::bound                     = &quot;flat&quot;</div><div>GRHydro::rho_abs_min               = 1.0e-6</div><div>GRHydro::GRHydro_atmo_tolerance    = 0.</div><div>GRHydro::c2p_reset_pressure        = &quot;yes&quot;</div><div>GRHydro::GRHydro_eos_type           = &quot;General&quot;</div><div>GRHydro::GRHydro_eos_table          = &quot;Ideal_Fluid&quot;</div><div>GRHydro::GRHydro_MaxNumSandRVars = 0</div><div><br></div><div>GRHydro::use_enhanced_ppm            = &quot;yes&quot;</div><div>GRHydro::sync_conserved_only     = &quot;yes&quot;<br></div><div>GRHydro::reconstruct_Wv          = &quot;yes&quot;</div><div>GRHydro::c2p_resort_to_bisection = &quot;yes&quot;</div><div>GRHydro::use_cxx_code            = &quot;yes&quot;</div><div><br></div><div>TmunuBase::stress_energy_storage = &quot;yes&quot;</div><div>TmunuBase::stress_energy_at_RHS  = &quot;yes&quot;</div><div>TmunuBase::timelevels            =  1</div><div>TmunuBase::prolongation_type     = &quot;none&quot;</div><div>TmunuBase::support_old_CalcTmunu_mechanism = &quot;no&quot;</div><div><br></div><div>EOS_Omni::gl_gamma = 2</div><div>EOS_Omni::gl_k=100</div><div>EOS_Omni::poly_k = 100</div></div><div>&quot;</div><div>Yours,</div><div>Zhi-Chao, Zhao</div><div>Institute of High Energy Physics, CAS, China.</div><div><br></div></div>