OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

EclipseWriter: make it not crash if no wells are present

Browse Source
This commit is contained in:
Andreas Lauser 2014-07-15 13:15:34 +02:00
parent 89cbfac621
commit c6d3994722
1 changed files with 13 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/core/io/eclipse/EclipseWriter.cpp
View File

@ -682,17 +682,25 @@ protected:
{ {
// convert m^3/s of injected fluid to m^3/d of produced fluid // convert m^3/s of injected fluid to m^3/d of produced fluid
const double convFactor = Opm::unit::convert::to(1., Opm::unit::day); const double convFactor = Opm::unit::convert::to(1., Opm::unit::day);
const double value = sign_ * wellState.wellRates()[index_] * convFactor; double value = 0;
if (wellState.wellRates().size() > 0) {
assert(int(wellState.wellRates().size()) > index_);
value = sign_ * wellState.wellRates()[index_] * convFactor;
}
return value; return value;
} }
double bhp(const WellState& wellstate) double bhp(const WellState& wellstate)
{ {
if (wellstate.bhp().size() > 0) {
// Note that 'index_' is used here even though it is meant // Note that 'index_' is used here even though it is meant
// to give a (well,phase) pair. // to give a (well,phase) pair.
const int num_phases = wellstate.wellRates().size() / wellstate.bhp().size(); const int num_phases = wellstate.wellRates().size() / wellstate.bhp().size();
return wellstate.bhp()[index_/num_phases]; return wellstate.bhp()[index_/num_phases];
} }
return 0.0;
}
}; };
/// Monitors the rate given by a well. /// Monitors the rate given by a well.