diff --git a/ebos/ecltransmissibility.hh b/ebos/ecltransmissibility.hh index c9cd29c83..e824701f9 100644 --- a/ebos/ecltransmissibility.hh +++ b/ebos/ecltransmissibility.hh @@ -759,22 +759,25 @@ private: // over several processes.) const auto& fp = vanguard_.eclState().fieldProps(); if (fp.has_double("PERMX")) { - const std::vector& permxData = fp.get_global_double("PERMX"); + const std::vector& permxData = fp.get_double("PERMX"); - std::vector permyData(permxData); + std::vector permyData; if (fp.has_double("PERMY")) - permyData = fp.get_global_double("PERMY"); + permyData = fp.get_double("PERMY"); + else + permyData = permxData; - std::vector permzData(permxData); + std::vector permzData; if (fp.has_double("PERMZ")) - permzData = fp.get_global_double("PERMZ"); + permzData = fp.get_double("PERMZ"); + else + permzData = permxData; for (size_t dofIdx = 0; dofIdx < numElem; ++ dofIdx) { - unsigned cartesianElemIdx = vanguard_.cartesianIndex(dofIdx); permeability_[dofIdx] = 0.0; - permeability_[dofIdx][0][0] = permxData[cartesianElemIdx]; - permeability_[dofIdx][1][1] = permyData[cartesianElemIdx]; - permeability_[dofIdx][2][2] = permzData[cartesianElemIdx]; + permeability_[dofIdx][0][0] = permxData[dofIdx]; + permeability_[dofIdx][1][1] = permyData[dofIdx]; + permeability_[dofIdx][2][2] = permzData[dofIdx]; } // for now we don't care about non-diagonal entries