diff options
Diffstat (limited to 'src/Initialisation.c')
-rw-r--r-- | src/Initialisation.c | 23 |
1 files changed, 22 insertions, 1 deletions
diff --git a/src/Initialisation.c b/src/Initialisation.c index 9aa31ec..762b402 100644 --- a/src/Initialisation.c +++ b/src/Initialisation.c @@ -20,7 +20,14 @@ void HydroBase_Zero (CCTK_ARGUMENTS) w_lorentz[i] = 1.0; eps[i ] = 0.0; } - + + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 1) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + Abar[i ] = 0.0; + } + } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 1) { #pragma omp parallel for for (int i=0; i<np; ++i) { @@ -55,6 +62,13 @@ void HydroBase_Zero (CCTK_ARGUMENTS) } } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 2) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + Abar_p[i ] = 0.0; + } + } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 2) { #pragma omp parallel for for (int i=0; i<np; ++i) { @@ -86,6 +100,13 @@ void HydroBase_Zero (CCTK_ARGUMENTS) "Too many active time levels for HydroBase variables"); } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::Abar") >= 3) { +#pragma omp parallel for + for (int i=0; i<np; ++i) { + Abar_p_p[i ] = 0.0; + } + } + if (CCTK_ActiveTimeLevels(cctkGH, "HydroBase::temperature") >= 3) { #pragma omp parallel for for (int i=0; i<np; ++i) { |