aboutsummaryrefslogtreecommitdiff
path: root/schedule.ccl
diff options
context:
space:
mode:
Diffstat (limited to 'schedule.ccl')
-rw-r--r--schedule.ccl60
1 files changed, 58 insertions, 2 deletions
diff --git a/schedule.ccl b/schedule.ccl
index a72ade7..200f405 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -1,7 +1,7 @@
# Schedule definitions for thorn ADMBase
# $Header$
-STORAGE: shift_state, metric[1], curv[1], lapse[1]
+STORAGE: shift_state, dtlapse_state, dtshift_state, metric[1], curv[1], lapse[1]
SCHEDULE ADMBase_ParamCheck at CCTK_PARAMCHECK
{
@@ -44,7 +44,7 @@ if(CCTK_Equals(initial_lapse, "one"))
}
# Do we have storage for the shift ?
-if (!CCTK_Equals(initial_shift,"none"))
+if (!CCTK_Equals(initial_shift, "none"))
{
STORAGE: shift[1]
@@ -71,6 +71,62 @@ if(CCTK_Equals(initial_shift, "zero"))
} "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
{