diff options
Diffstat (limited to 'schedule.ccl')
-rw-r--r-- | schedule.ccl | 60 |
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 { |