OilfieldToolsNet/opm-common
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Templated internalizing of WSEGSICD keyword

Browse Source
This commit is contained in:
Joakim Hove
2020-06-11 09:32:24 +02:00
parent 4c100d0bd5
commit 73a244ceda
2 changed files with 33 additions and 26 deletions
Download Patch File Download Diff File Expand all files Collapse all files

32
opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.hpp
View File

@@ -27,12 +27,11 @@
#include <string>
#include <opm/parser/eclipse/EclipseState/Schedule/MSW/icd.hpp>
#include <opm/parser/eclipse/Deck/DeckRecord.hpp>
#include <opm/parser/eclipse/Deck/DeckKeyword.hpp>
namespace Opm {
class DeckRecord;
class DeckKeyword;
class SICD {
public:
@@ -91,6 +90,33 @@ namespace Opm {
serializer(m_scaling_factor);
}
template<class ICD>
static 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("SEG1").get<int>(0);
const int end_segment = record.getItem("SEG2").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;
}
private:
double m_strength;
double m_length;

27
src/opm/parser/eclipse/EclipseState/Schedule/MSW/SICD.cpp
View File

@@ -101,32 +101,13 @@ namespace Opm {
return result;
}
std::map<std::string, std::vector<std::pair<int, SICD> > >
SICD::fromWSEGSICD(const DeckKeyword& wsegsicd)
{
std::map<std::string, std::vector<std::pair<int, SICD> > > res;
for (const DeckRecord &record : wsegsicd) {
const std::string well_name = record.getItem("WELL").getTrimmedString(0);
const int start_segment = record.getItem("SEG1").get<int>(0);
const int end_segment = record.getItem("SEG2").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 SICD 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;
return SICD::fromWSEG<SICD>(wsegsicd);
}
double SICD::maxAbsoluteRate() const {