<div dir="ltr"><div dir="ltr"><div>Hello,</div><div><br></div><div>I am following the Parma group&#39;s par file to simulate binary neutron stars merging into a black hole. In particular I adopt the following lines into my parfile but it looks like everything is &quot;excised&quot; at the very first iteration (the maximum density suddenly becomes ~10e-11, basically the atmosphere density settings).</div><div><br></div><div>Could someone please guide me how this should be done properly? or what parameter settings do I need to set as well?</div><div><br></div><div>With Regards,</div><div>Ken</div><div><br></div><div><br></div><div>#------------------------------------------------------------------------------</div><div># AHFinder</div><div><br></div><div>ActiveThorns = &quot;AHFinderDirect&quot;</div><div><br></div><div># This is later steered</div><div>AHFinderDirect::find_every = 0</div><div><br></div><div>AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no</div><div><br></div><div>AHFinderDirect::move_origins            = yes</div><div>AHFinderDirect::reshape_while_moving    = yes</div><div>AHFinderDirect::predict_origin_movement = yes</div><div><br></div><div>AHFinderDirect::geometry_interpolator_name = &quot;Lagrange polynomial interpolation&quot;</div><div>AHFinderDirect::geometry_interpolator_pars = &quot;order=4&quot;</div><div>AHFinderDirect::surface_interpolator_name  = &quot;Lagrange polynomial interpolation&quot;</div><div>AHFinderDirect::surface_interpolator_pars  = &quot;order=4&quot;</div><div><br></div><div>AHFinderDirect::output_h_every = 128</div><div><br></div><div>AHFinderDirect::N_horizons = 1</div><div><br></div><div>AHFinderDirect::which_surface_to_store_info          [1] = 0</div><div>AHFinderDirect::reset_horizon_after_not_finding      [1] = no</div><div>AHFinderDirect::initial_guess__coord_sphere__radius  [1] = 1.3528</div><div><br></div><div># Hydroexcision</div><div>ActiveThorns = &quot;SetMask_SphericalSurface&quot;</div><div><br></div><div>SetMask_SphericalSurface::SetMask_SurfaceIndex[0] = 0</div><div>SetMask_SphericalSurface::SetMask_RadiusFactor[0] = 0.6</div><div>SetMask_SphericalSurface::SetMask_TrueShape[0]    = &quot;false&quot;</div></div></div>