Loading...
Searching...
No Matches
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
Opm::DeckKeyword
Definition DeckKeyword.hpp:36
Opm::DeckRecord
Definition DeckRecord.hpp:32
Opm
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition Exceptions.hpp:30