diff --git a/opm/core/props/rock/RockFromDeck.cpp b/opm/core/props/rock/RockFromDeck.cpp index 10ee7bf75..50e6982a2 100644 --- a/opm/core/props/rock/RockFromDeck.cpp +++ b/opm/core/props/rock/RockFromDeck.cpp @@ -121,8 +121,10 @@ namespace Opm // K(i,j) = (*tensor[kmap[kix]])[glob]; permeability_[off + (i + dim*j)] = (*tensor[kmap[kix]])[glob]; } + // K(i,i) = std::max(K(i,i), perm_threshold); - permeability_[off + 3*i + i] = std::max(permeability_[off + 3*i + i], perm_threshold); + double& kii = permeability_[off + i*(dim + 1)]; + kii = std::max(kii, perm_threshold); } permfield_valid_[c] = std::vector::value_type(1);