Commit ee8ee427 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Removed const from level_0 and level_N (OpenMP restriction).

parent 86a77401
Pipeline #4341 passed with stages
in 15 minutes
......@@ -227,8 +227,8 @@ vert_gen_weights3d1d(bool expol, size_t gridsize, int nlev1, const Varray<T> &xl
const auto lup = levelDirUp(nlev1, lev1p2[0].data());
const auto ldown = levelDirDown(nlev1, lev1p2[0].data());
if (!lup && !ldown) cdoAbort("Non monotonic zaxis!");
const double level_0 = lup ? -1.e33 : 1.e33;
const double level_N = lup ? 1.e33 : -1.e33;
double level_0 = lup ? -1.e33 : 1.e33;
double level_N = lup ? 1.e33 : -1.e33;
#ifdef _OPENMP
#pragma omp parallel for default(none) shared(expol, gridsize, nlev1, nlev2, lev2, level_0, level_N, lev1p2, xlev1, xlev_idx1, xlev_idx2, xlev_wgt1, xlev_wgt2, lev_idx1, lev_idx2, lev_wgt1, lev_wgt2)
......
......@@ -60,8 +60,8 @@ vert_gen_weights3d(bool expol, size_t gridsize, int nlev1, Varray<float> &xlev1,
const auto lup = levelDirUp(nlev1, lev1p2[0].data());
const auto ldown = levelDirDown(nlev1, lev1p2[0].data());
if (!lup && !ldown) cdoAbort("Non monotonic zaxis!");
const double level_0 = lup ? -1.e33 : 1.e33;
const double level_N = lup ? 1.e33 : -1.e33;
double level_0 = lup ? -1.e33 : 1.e33;
double level_N = lup ? 1.e33 : -1.e33;
#ifdef _OPENMP
#pragma omp parallel for default(none) shared(expol, gridsize, nlev1, nlev2, lev2, xlev2, level_0, level_N, lev1p2, xlev1, xlev_idx1, xlev_idx2, xlev_wgt1, xlev_wgt2, lev_idx1, lev_idx2, lev_wgt1, lev_wgt2)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment