FromWSEG.hpp

/*
  Copyright 2020 Equinor 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 FROM_WSEG_HPP
#define FROM_WSEG_HPP
 
#include <opm/input/eclipse/Deck/DeckKeyword.hpp>
 
#include <map>
#include <string>
#include <vector>
 
namespace Opm {
 
template<class ICD>
std::map<std::string, std::vector<std::pair<int, ICD> > >
fromWSEG(const DeckKeyword& wseg) {
         std::map<std::string, std::vector<std::pair<int, ICD> > > res;
 
    for (const DeckRecord &record : wseg) {
        const std::string well_name = record.getItem("WELL").getTrimmedString(0);
 
        const int start_segment = record.getItem("SEGMENT1").get<int>(0);
        const int end_segment = record.getItem("SEGMENT2").get<int>(0);
 
        if (start_segment < 2 || end_segment < 2 || end_segment < start_segment) {
            const std::string message = "Segment numbers " + std::to_string(start_segment) + " and "
                + std::to_string(end_segment) + " specified in WSEGSICD for well " +
                well_name
                + " are illegal ";
            throw std::invalid_argument(message);
        }
 
        const ICD spiral_icd(record);
        for (int seg = start_segment; seg <= end_segment; seg++) {
            res[well_name].push_back(std::make_pair(seg, spiral_icd));
        }
    }
    return res;
}
 
}
 
#endif