# Schedule definitions for thorn Einstein
# $Header$

STORAGE: metric, curv, lapse
COMMUNICATION: metric, curv, lapse

if (!CCTK_Equals(shift,"none")) 
{
  STORAGE: shift
}

if (use_conformal) {
  STORAGE: confac
}

if (use_conformal_derivs) {
  STORAGE: confac_1derivs,confac_2derivs
}

schedule Einstein_InitSymBound at CCTK_BASEGRID
{
  LANG: Fortran
} "Set up symmetries for Einstein GFs"
 
schedule Einstein_ActivateSlicing at CCTK_BASEGRID
{
  LANG: C
} "Initialize slicing, setup priorities for mixed slicings"

schedule InitialEinstein at CCTK_INITIAL 
{
    LANG: Fortran
} "Initialisation for Einstein methods"

if (CCTK_Equals(initial_data,"flat")) 
{
  schedule InitialFlat at CCTK_INITIAL after InitialEinstein
  {
    LANG: Fortran
  } "Flat initial data"
}

if (CCTK_Equals(initial_lapse,"one"))
{
  schedule LapseOne at CCTK_INITIAL
  {
    LANG: Fortran
  } "Set initial lapse to one"
}

if (CCTK_Equals(initial_lapse,"gaussian"))
{
  schedule LapseGaussian at CCTK_INITIAL
  {
    LANG: Fortran
  } "Set initial lapse to a gaussian"
}

schedule Einstein_SetNextSlicing at CCTK_PRESTEP
{
  LANG: C
} "Identify the slicing for the next iteration"

schedule evaltrK at CCTK_ANALYSIS
{
  STORAGE: trace_of_K,detofg
  LANG: Fortran
  TRIGGERS: trK,detg
} "Compute the trace of the extrinsic curvature"

schedule metric_carttosphere at CCTK_ANALYSIS
{
  STORAGE: spherical_metric
  LANG: Fortran
  TRIGGERS: grr,grp,grq,gpp,gqp,gqq
} "Calculate the spherical metric in r,theta(q), phi(p)"

schedule curv_carttosphere at CCTK_ANALYSIS
{
  STORAGE: spherical_curv
  LANG: Fortran
  TRIGGERS: hrr,hrp,hrq,hpp,hqp,hqq
} "Calculate the spherical ex. curvature in r, theta(q), phi(p)"

#if (use_courant) 
#{
#  schedule CalcCourant at CCTK_PRESTEP
#  {
#   LANG: C
#  } "Calculate the timestep base on the Courant condition"
#}