Initialize well solvent fraction

The well solvent fraction is set to zero as default.
This commit is contained in:
Tor Harald Sandve 2015-09-17 15:59:39 +02:00
parent 18613b9c7b
commit 23f4f48df5

View File

@ -97,9 +97,11 @@ namespace Opm
const Wells* wells)
{
// compute solvent inflow
const int nw = wells->number_of_wells;
std::vector<double> perfcells_fraction(wells->well_connpos[nw]);
std::fill(perfcells_fraction.begin(), perfcells_fraction.end(), 0.0);
if (deck_->hasKeyword("WSOLVENT")) {
const int nw = wells->number_of_wells;
std::vector<double> perfcells_fraction(wells->well_connpos[nw]);
size_t currentStep = timer.currentStepNum();
ScheduleConstPtr schedule = BaseType::eclipse_state_->getSchedule();
@ -135,11 +137,10 @@ namespace Opm
}
}
}
well_state.solventFraction() = perfcells_fraction;
}
well_state.solventFraction() = perfcells_fraction;
}
} // namespace Opm
#endif // OPM_SIMULATORFULLYIMPLICITBLACKOILSOLVENT_IMPL_HEADER_INCLUDED