Update of hysteresis model is now performed prior

to each time step, thus also capturing the initial state.
This commit is contained in:
osae 2014-08-12 19:53:38 +02:00
parent 8e679382dd
commit 6eea5e22a1

View File

@ -292,6 +292,9 @@ namespace Opm
// Max oil saturation
props_.updateSatOilMax(state.saturation());
// Hysteresis
props_.updateSatHyst(state.saturation(), allcells_);
// Solve pressure equation.
// if (check_well_controls_) {
// computeFractionalFlow(props_, allcells_,
@ -339,9 +342,6 @@ namespace Opm
computePorevolume(AutoDiffGrid::numCells(grid_), AutoDiffGrid::beginCellVolumes(grid_), props_.porosity(), *rock_comp_props_, state.pressure(), porevol);
}
// Hysteresis
props_.updateSatHyst(state.saturation(), allcells_);
sreport.total_time = step_timer.secsSinceStart();
if (output_) {
sreport.reportParam(tstep_os);