[tyskies]

Hello everyone,

I'm having difficulties to access the y coordinate of the wall patches of my test case.
Here is the code I use:

    const fvPatchList& patches = this->mesh_.boundary();
    forAll(patches, patchi)
    {
        const fvPatch& curPatch = patches[patchi];
        if (isA<wallFvPatch>(curPatch))
        {
          forAll(curPatch, facei)
          {
              label celli = curPatch.faceCells()[facei];
              const double y = curPatch.Cf()[facei].y(); // 0.0003;
              Info << "\n\n" << y << endl;
              epsilon_[celli] = 2.0 * nu_[celli] * k_[celli] / y / y;
              // epsilon_[celli] = 2.0 * nu_[celli] * k_[celli] / mag(y) / mag(y);
           }
       }
    }

Despite the fact that no errors are detected during the compilation, the simulation crashes before it completes the first step.

I found that the simulation crashes at the line where epsilon_[celli] is evaluated and that it is because y is equal to 0...

Do you have any idea why my code isn't working?

Thanks for your support!