Dear all,
I was expecting the two following snippets to be semantically equivalent
for two-phase water+oil flow:
auto& well_perfrates = well_state.perfRates();for (int w = 0; w <
wells->number_of_wells; ++w) {
double well_rate = 0.;
for (int perf = wells->well_connpos[w]; perf <
wells->well_connpos[w+1]; ++perf) {
well_rate += perf_rate;
}
cout << well_rate << '\t';}
-------------------------------------------------------------------------------------------------------------
auto& well_rates = well_state.wellRates();for (int w = 0; w <
wells->number_of_wells; ++w) {
double water_rate = well_rates[w];
double oil_rate = well_rates[w+1];
cout << water_rate + oil_rate << '\t';}
Why does the later print 0 everywhere? WellState::wellrates_ is not
implemented yet? Did I misunderstand something?
Sincerely,
Júlio.
_______________________________________________
Opm mailing list
[email protected]
http://www.opm-project.org/mailman/listinfo/opm
