intialize RHS to zero in slow evolution sector

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

1 files changed, 26 insertions, 0 deletions

diff --git a/src/InitialCopy.c b/src/InitialCopy.c
index 650f248..3810ab6 100644
--- a/src/InitialCopy.c
+++ b/src/InitialCopy.c
@@ -590,6 +590,32 @@ void MoL_InitRHS(CCTK_ARGUMENTS)
     MoL_LinearCombination(cctkGH, RHSVariableIndex[var], rl, tl, 0.0,
                           NULL, NULL, NULL, 0);
   }
+  
+  for (var = 0; var < MoLNumEvolvedVariablesSlow; var++)
+  {
+    StorageOn = CCTK_QueryGroupStorageI(cctkGH,
+                                       CCTK_GroupIndexFromVarI(RHSVariableIndexSlow[var]));
+    
+    if (StorageOn < 0)
+    {
+      CCTK_VWarn(1,__LINE__,__FILE__,CCTK_THORNSTRING,"Warning for index %i", 
+                 RHSVariableIndexSlow[var]);
+      CCTK_WARN(0, "The index passed does not correspond to a GF.");
+    }
+    else if (StorageOn == 0) {
+#ifdef MOLDEBUG
+      printf("Aargh! Vars %d var %d index %d name %s\n",
+             MoLNumEvolvedVariablesSlow, var, RHSVariableIndexSlow[var],  
+                 CCTK_VarName(RHSVariableIndexSlow[var]));
+#endif
+      CCTK_VWarn(1,__LINE__,__FILE__,CCTK_THORNSTRING,"Warning for GF %s", 
+                 CCTK_VarName(RHSVariableIndexSlow[var]));
+      CCTK_WARN(0, "The grid function does not have storage assigned.");
+    }
+    
+    MoL_LinearCombination(cctkGH, RHSVariableIndexSlow[var], rl, tl, 0.0,
+                          NULL, NULL, NULL, 0);
+  }
     
   for (var = 0; var < MoLNumEvolvedArrayVariables; var++)
   {