Reinstate correct porevolume usage.

The porevolume at the start of the timestep should be provided to
the solver.
This commit is contained in:
Atgeirr Flø Rasmussen 2013-03-15 15:10:16 +01:00
parent 63c7380d04
commit e42cfc1303

View File

@ -543,7 +543,7 @@ namespace Opm
double injected[2] = { 0.0 };
double produced[2] = { 0.0 };
for (int tr_substep = 0; tr_substep < num_transport_substeps_; ++tr_substep) {
tsolver_->solve(&porevol[0], &transport_src[0], stepsize, state);
tsolver_->solve(&initial_porevol[0], &transport_src[0], stepsize, state);
double substep_injected[2] = { 0.0 };
double substep_produced[2] = { 0.0 };