diff options
Diffstat (limited to 'src/InitialCopy.c')
| -rw-r--r-- | src/InitialCopy.c | 26 |
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++) { |