schedule BaikalX_BSSN_to_ADM in ODESolvers_RHS after BaikalX_T4UU before BaikalX_Ricci

schedule BaikalX_BSSN_to_ADM in ODESolvers_RHS after BaikalX_zero_rhss after BaikalX_T4UU before BaikalX_Ricci

  WRITES: ADMBase::dtshift(everywhere)

  WRITES: BaikalX::evol_variables_rhs(interior)

  WRITES: BaikalX::evol_variables_rhs(everywhere)
  # this applies boundary conditions and since ODESolvers_RHS is called after ODESolver's
  # ADD, this is the equivalent of MoL_PostStep
  SYNC: BaikalX::evol_variables

schedule enforce_detgammabar_constraint in MoL_PostStep before BC_Update

schedule enforce_detgammabar_constraint in ODESolvers_RHS before BaikalX_zero_rhss

# Next update ADM quantities
schedule BaikalX_BSSN_to_ADM in MoL_PostStep after BaikalX_ApplyBCs before ADMBase_SetADMVars
{
  LANG: C
  READS: BaikalX::evol_variables(everywhere)
  WRITES: ADMBase::metric(everywhere)
  WRITES: ADMBase::curv(everywhere)
  WRITES: ADMBASE::shift(everywhere)
  WRITES: ADMBASE::dtshift(everywhere)
  WRITES: ADMBASE::lapse(everywhere)
} "Perform BSSN-to-ADM conversion. Useful for diagnostics."

#include "Symmetry.h"