Add explanatory comment.

This commit is contained in:
Atgeirr Flø Rasmussen 2012-10-02 15:46:33 +02:00
parent 5acbf00b14
commit a9783eefc7

View File

@ -77,8 +77,12 @@ namespace Opm
std::vector<double> prod_surfvol(np);
for (int c = 0; c < num_cells; ++c) {
if (transport_src[c] > 0.0) {
// Inflowing transport source is a surface volume flux
// for the first phase.
injected[0] += transport_src[c]*dt;
} else if (transport_src[c] < 0.0) {
// Outflowing transport source is a total reservoir
// volume flux.
const double flux = -transport_src[c]*dt;
const double* sat = &s[np*c];
props.relperm(1, sat, &c, &mob[0], 0);