Initialise the RHS with zero before scheduling the physics routines.

This is necessary because MoL switches off boundary prolongation when
there is mesh refinement, leaving certain grid points undefined.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@54 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b

1 files changed, 9 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 83832c5..8d71041 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -444,6 +444,15 @@ if (CCTK_Equals(ODE_Method,"ICN-avg"))
   } "Averages the time levels for the averaging ICN method"
 }
 
+##############################################
+###  The time step is initialised to zero  ###
+##############################################
+
+schedule MoL_InitRHS IN MoL_Step BEFORE MoL_CalcRHS
+{
+  LANG: C
+} "Initialise the RHS functions"
+
 #####################################################
 ###  The group where all the physics takes place  ###
 #####################################################