# Schedule definitions for thorn ADMBase
# $Header$

STORAGE: shift_state, metric[1], curv[1], lapse[1]

SCHEDULE ADMBase_ParamCheck at CCTK_PARAMCHECK
{
  LANG: C
  OPTIONS: global
} "Check consistency of parameters"

# Schedule group for setting initial data

SCHEDULE group ADMBase_InitialData at CCTK_INITIAL
{
} "Schedule group for calculating ADM initial data"

SCHEDULE group ADMBase_InitialGauge at CCTK_INITIAL after ADMBase_InitialData
{
} "Schedule group for the ADM initial gauge condition"

SCHEDULE group ADMBase_PostInitial at CCTK_INITIAL after ADMBase_InitialGauge
{
} "Schedule group for modifying the ADM initial data, such as e.g. adding noise"

# Initial data for the spatial slices

if(CCTK_Equals(initial_data, "Cartesian Minkowski"))
{
  SCHEDULE ADMBase_CartesianMinkowski in ADMBase_InitialData
  {
    LANG: C
  } "Set the metric and extrinsic curvature to cartesian minkowski values"
}

# Initial data for the lapse

if(CCTK_Equals(initial_lapse, "one"))
{
  SCHEDULE ADMBase_LapseOne in ADMBase_InitialGauge
  {
    LANG: C
  } "Set the lapse to 1 at all points"
}

# Do we have storage for the shift ?
if (!CCTK_Equals(initial_shift,"none")) 
{
  STORAGE: shift[1]

  SCHEDULE ADMBase_SetShiftStateOn at CCTK_BASEGRID
  {
    LANG: C
  } "Set the shift_state variable to 1"
}
else
{
  SCHEDULE ADMBase_SetShiftStateOff at CCTK_BASEGRID
  {
    LANG: C
  } "Set the shift_state variable to 0"
}

# Initial data for the shift

if(CCTK_Equals(initial_shift, "zero"))
{
  SCHEDULE ADMBase_ShiftZero in ADMBase_InitialGauge
  {
    LANG: C
  } "Set the shift to 0 at all points"
}

# FIXME: Remove this when symmetries done better 
schedule Einstein_InitSymBound at CCTK_WRAGH
{
  LANG: C
  OPTIONS: global
} "Set up GF symmetries"