Refactor inflow tests and flux assignments.

This commit is contained in:
Bård Skaflestad 2012-02-10 11:40:21 +01:00
parent deb868c233
commit 8a2d216a83

View File

@ -76,8 +76,10 @@ namespace Opm
{
cell = cell_index;
s0 = tm.saturation_[cell];
influx = tm.source_[cell] > 0 ? -tm.source_[cell] : 0.0;
outflux = tm.source_[cell] <= 0 ? -tm.source_[cell] : 0.0;
double dflux = -tm.source_[cell];
bool src_is_inflow = dflux < 0.0;
influx = src_is_inflow ? dflux : 0.0;
outflux = !src_is_inflow ? dflux : 0.0;
dtpv = tm.dt_/tm.porevolume_[cell];
for (int i = tm.grid_.cell_facepos[cell]; i < tm.grid_.cell_facepos[cell+1]; ++i) {