EquilWrapper.hpp

Go to the documentation of this file.

/*
  Copyright (C) 2014 by Andreas Lauser

  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_PARSER_EQUIL_WRAPPER_HPP
#define OPM_PARSER_EQUIL_WRAPPER_HPP

#include <opm/parser/eclipse/Deck/Deck.hpp>

#include <vector>
#include <algorithm>

namespace Opm {
    class EquilWrapper {
    public:
        EquilWrapper(Opm::DeckKeywordConstPtr keyword)
            : m_keyword(keyword)
        {
        }

        int numRegions() const
        { return m_keyword->size(); }

        double datumDepth(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(0)->getSIDouble(0); }

        double datumDepthPressure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(1)->getSIDouble(0); }

        double waterOilContactDepth(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(2)->getSIDouble(0); }

        double waterOilContactCapillaryPressure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(3)->getSIDouble(0); }

        double gasOilContactDepth(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(4)->getSIDouble(0); }

        double gasOilContactCapillaryPressure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(5)->getSIDouble(0); }

        int liveOilInitProceedure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(6)->getInt(0); }

        int wetGasInitProceedure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(7)->getInt(0); }

        int initializationTargetAccuracy(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(8)->getInt(0); }

        int compositionalInitializationProceedure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(9)->getInt(0); }

        bool useSaturationPressure(int regionIdx) const
        { return m_keyword->getRecord(regionIdx)->getItem(10)->getInt(0) == 0; }

    private:
        Opm::DeckKeywordConstPtr m_keyword;
    };
}

#endif  // OPM_PARSER_EQUIL_WRAPPER_HPP