# 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"
 
if (einstein_register_slicing)
{
   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"
}

if (einstein_register_slicing)
{
   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: trace_of_K,detofg
} "Compute the trace of the extrinsic curvature"

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

schedule curv_carttosphere at CCTK_ANALYSIS
{
  STORAGE: spherical_curv
  LANG: Fortran
  TRIGGERS: spherical_curv
} "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"
#}