# Schedule definitions for thorn ADMBase
# $Header$

STORAGE: shift_state, dtlapse_state, dtshift_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"
}

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

  SCHEDULE ADMBase_SetDtlapseStateOn at CCTK_BASEGRID
  {
    LANG: C
  } "Set the dtlapse_state variable to 1"
}
else
{
  SCHEDULE ADMBase_SetDtlapseStateOff at CCTK_BASEGRID
  {
    LANG: C
  } "Set the dtlase_state variable to 0"
}

# Initial data for dtlapse

if(CCTK_Equals(initial_dtlapse, "zero"))
{
  SCHEDULE ADMBase_DtlapseZero in ADMBase_InitialGauge
  {
    LANG: C
  } "Set dtlapse to 0 at all points"
}

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

  SCHEDULE ADMBase_SetDtshiftStateOn at CCTK_BASEGRID
  {
    LANG: C
  } "Set the dtshift_state variable to 1"
}
else
{
  SCHEDULE ADMBase_SetDtshiftStateOff at CCTK_BASEGRID
  {
    LANG: C
  } "Set the dtlase_state variable to 0"
}

# Initial data for dtshift

if(CCTK_Equals(initial_dtshift, "zero"))
{
  SCHEDULE ADMBase_DtshiftZero in ADMBase_InitialGauge
  {
    LANG: C
  } "Set dtshift 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"