VFPProdProperties::ADB VFPProdProperties::bhp(const std::vector<int>& table_id, const Wells& wells, const ADB& qs, const ADB& thp_arg, const ADB& alq) const { const int nw = wells.number_of_wells; //Short-hands for water / oil / gas phases //TODO enable support for two-phase. assert(wells.number_of_phases == 3); const ADB& w = subset(qs, Span(nw, 1, BlackoilPhases::Aqua*nw)); const ADB& o = subset(qs, Span(nw, 1, BlackoilPhases::Liquid*nw)); const ADB& g = subset(qs, Span(nw, 1, BlackoilPhases::Vapour*nw)); return bhp(table_id, w, o, g, thp_arg, alq); }
/// The number of wells present. int numWells() const { return bhp().size(); }