ECLPvtCurveCollection.hpp

Go to the documentation of this file.

/*
  Copyright 2017 Statoil ASA.
 
  This file is part of the Open Porous Media Project (OPM).
 
  OPM is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  OPM is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.
 
  You should have received a copy of the GNU General Public License
  along with OPM.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef OPM_ECLPVTCURVECOLLECTION_HEADER_INCLUDED
#define OPM_ECLPVTCURVECOLLECTION_HEADER_INCLUDED
 
#include <opm/utility/ECLGraph.hpp>
#include <opm/utility/ECLPhaseIndex.hpp>
#include <opm/utility/ECLPropertyUnitConversion.hpp>
#include <opm/utility/ECLPvtCommon.hpp>
#include <opm/utility/ECLPvtGas.hpp>
#include <opm/utility/ECLPvtOil.hpp>
#include <opm/utility/ECLUnitHandling.hpp>
 
#include <memory>
#include <vector>
 
 
namespace Opm {
    class ECLInitFileData;
} // Opm
 
namespace Opm { namespace ECLPVT {
 
    class ECLPvtCurveCollection
    {
    public:
        ECLPvtCurveCollection(const ECLGraph&        G,
                              const ECLInitFileData& init);
 
        void setOutputUnits(std::unique_ptr<const ECLUnits::UnitSystem> usys);
 
        std::vector<PVTGraph>
        getPvtCurve(const RawCurve      curve,
                    const ECLPhaseIndex phase,
                    const int           pvtRegionID) const;
 
        std::vector<double>
        getDynamicPropertySI(const RawCurve             property,
                             const ECLPhaseIndex        phase,
                             const int                  pvtRegionID,
                             const std::vector<double>& phasePress,
                             const std::vector<double>& mixRatio
                                 = std::vector<double>()) const;
 
        std::vector<double>
        getDynamicPropertyNative(const RawCurve      property,
                                 const ECLPhaseIndex phase,
                                 const int           pvtRegionID,
                                 std::vector<double> phasePress, // Mutable copy
                                 std::vector<double> mixRatio    // Mutable copy
                                 = std::vector<double>()) const;
 
    private:
        std::shared_ptr<Gas> gas_; // shared => default special member funcs.
 
        std::shared_ptr<Oil> oil_;
 
        std::shared_ptr<const ECLUnits::UnitSystem> usys_native_;
 
        std::shared_ptr<const ECLUnits::UnitSystem> usys_internal_;
 
        std::shared_ptr<const ECLUnits::UnitSystem> usys_output_{nullptr};
 
        bool isValidRequest(const ECLPhaseIndex phase,
                            const int           pvtRegionID) const;
 
        std::vector<PVTGraph>
        convertToOutputUnits(std::vector<PVTGraph>&& graph,
                             const RawCurve          curve,
                             const ECLPhaseIndex     phase) const;
    };
 
}} // Opm::ECLPVT
 
#endif // OPM_ECLPVTCURVECOLLECTION_HEADER_INCLUDED