Bugfix in calculation of adsorbed polymer for balance reports.

This commit is contained in:
Atgeirr Flø Rasmussen
2012-03-26 15:33:38 +02:00
parent c070f0ad4d
commit 294c8b03be
3 changed files with 9 additions and 4 deletions

View File

@@ -940,7 +940,7 @@ main(int argc, char** argv)
// Report volume balances.
Opm::computeSaturatedVol(porevol, state.saturation(), satvol);
polymass = Opm::computePolymerMass(porevol, state.saturation(), state.concentration(), polyprop.deadPoreVol());
polymass_adsorbed = Opm::computePolymerAdsorbed(polyprop, porevol, state.cmax());
polymass_adsorbed = Opm::computePolymerAdsorbed(*props, polyprop, porevol, state.cmax());
Opm::computeInjectedProduced(*props, polyprop, state.saturation(), state.concentration(),
src, simtimer.currentStepLength(), inflow_c,
injected, produced, polyinj, polyprod);