diff --git a/applications/ebos/eclwellmanager.hh b/applications/ebos/eclwellmanager.hh index a46813823..2c7bad5eb 100644 --- a/applications/ebos/eclwellmanager.hh +++ b/applications/ebos/eclwellmanager.hh @@ -651,7 +651,15 @@ protected: const Opm::Completion* completion = compInfo.first; std::shared_ptr eclWell = compInfo.second; - eclWell->setRadius(elemCtx, dofIdx, 0.5*completion->getDiameter()); + // the catch is a hack for a ideosyncrasy of opm-parser with regard to + // defaults handling: if the deck did not specify a radius for the + // completion, there seems to be no other way to detect this except for + // catching the exception + try { + eclWell->setRadius(elemCtx, dofIdx, 0.5*completion->getDiameter()); + } + catch (const std::logic_error& e) + {} // overwrite the automatically computed effective // permeability by the one specified in the deck. Note: this