21 #ifndef OPM_TABLE_MANAGER_HPP
22 #define OPM_TABLE_MANAGER_HPP
28 #include <opm/parser/eclipse/EclipseState/Tables/DenT.hpp>
29 #include <opm/parser/eclipse/EclipseState/Tables/PvtgTable.hpp>
30 #include <opm/parser/eclipse/EclipseState/Tables/PvtgwTable.hpp>
31 #include <opm/parser/eclipse/EclipseState/Tables/PvtgwoTable.hpp>
32 #include <opm/parser/eclipse/EclipseState/Tables/PvtoTable.hpp>
33 #include <opm/parser/eclipse/EclipseState/Tables/PvtsolTable.hpp>
34 #include <opm/parser/eclipse/EclipseState/Tables/RocktabTable.hpp>
35 #include <opm/parser/eclipse/EclipseState/Tables/Rock2dTable.hpp>
36 #include <opm/parser/eclipse/EclipseState/Tables/Rock2dtrTable.hpp>
37 #include <opm/parser/eclipse/EclipseState/Tables/PlyshlogTable.hpp>
38 #include <opm/parser/eclipse/EclipseState/Tables/PvtwsaltTable.hpp>
39 #include <opm/parser/eclipse/EclipseState/Tables/RwgsaltTable.hpp>
40 #include <opm/parser/eclipse/EclipseState/Tables/BrineDensityTable.hpp>
41 #include <opm/parser/eclipse/EclipseState/Tables/SolventDensityTable.hpp>
42 #include <opm/parser/eclipse/EclipseState/Tables/StandardCond.hpp>
43 #include <opm/parser/eclipse/EclipseState/Tables/FlatTable.hpp>
44 #include <opm/parser/eclipse/EclipseState/Tables/SorwmisTable.hpp>
45 #include <opm/parser/eclipse/EclipseState/Tables/SgcwmisTable.hpp>
46 #include <opm/parser/eclipse/EclipseState/Tables/MiscTable.hpp>
47 #include <opm/parser/eclipse/EclipseState/Tables/PmiscTable.hpp>
48 #include <opm/parser/eclipse/EclipseState/Tables/MsfnTable.hpp>
49 #include <opm/parser/eclipse/EclipseState/Tables/JFunc.hpp>
50 #include <opm/parser/eclipse/EclipseState/Tables/Tabdims.hpp>
51 #include <opm/parser/eclipse/EclipseState/Tables/TableContainer.hpp>
52 #include <opm/parser/eclipse/EclipseState/Tables/Aqudims.hpp>
53 #include <opm/parser/eclipse/EclipseState/Tables/PlymwinjTable.hpp>
54 #include <opm/parser/eclipse/EclipseState/Tables/SkprwatTable.hpp>
55 #include <opm/parser/eclipse/EclipseState/Tables/SkprpolyTable.hpp>
56 #include <opm/parser/eclipse/EclipseState/Tables/Eqldims.hpp>
57 #include <opm/parser/eclipse/EclipseState/Tables/Regdims.hpp>
58 #include <opm/parser/eclipse/EclipseState/Tables/TLMixpar.hpp>
71 const TableContainer& getTables(
const std::string& tableName )
const;
72 const TableContainer& operator[](
const std::string& tableName)
const;
73 bool hasTables(
const std::string& tableName )
const;
75 const Tabdims& getTabdims()
const;
76 const Eqldims& getEqldims()
const;
77 const Aqudims& getAqudims()
const;
78 const Regdims& getRegdims()
const;
83 size_t numFIPRegions()
const;
132 const JFunc& getJFunc()
const;
134 const std::vector<PvtgTable>& getPvtgTables()
const;
135 const std::vector<PvtgwTable>& getPvtgwTables()
const;
136 const std::vector<PvtgwoTable>& getPvtgwoTables()
const;
137 const std::vector<PvtoTable>& getPvtoTables()
const;
138 const std::vector<PvtsolTable>& getPvtsolTables()
const;
139 const std::vector<Rock2dTable>& getRock2dTables()
const;
140 const std::vector<Rock2dtrTable>& getRock2dtrTables()
const;
144 const DenT& WatDenT()
const;
145 const DenT& GasDenT()
const;
146 const DenT& OilDenT()
const;
148 std::size_t gas_comp_index()
const;
150 const std::vector<PvtwsaltTable>& getPvtwSaltTables()
const;
151 const std::vector<RwgsaltTable>& getRwgSaltTables()
const;
152 const std::vector<BrineDensityTable>& getBrineDensityTables()
const;
153 const std::vector<SolventDensityTable>& getSolventDensityTables()
const;
165 const std::map<int, PlymwinjTable>& getPlymwinjTables()
const;
166 const std::map<int, SkprwatTable>& getSkprwatTables()
const;
167 const std::map<int, SkprpolyTable>& getSkprpolyTables()
const;
168 const std::map<std::string, TableContainer>& getSimpleTables()
const;
185 double rtemp()
const;
187 double salinity()
const;
191 template<
class Serializer>
194 auto simpleTables = m_simpleTables;
195 auto split = splitSimpleTable(simpleTables);
196 serializer.map(simpleTables);
197 serializer(split.plyshMax);
198 serializer.map(split.plyshMap);
199 serializer(split.rockMax);
200 serializer.map(split.rockMap);
201 serializer.vector(m_pvtgTables);
202 serializer.vector(m_pvtgwTables);
203 serializer.vector(m_pvtgwoTables);
204 serializer.vector(m_pvtoTables);
205 serializer.vector(m_pvtsolTables);
206 serializer.vector(m_rock2dTables);
207 serializer.vector(m_rock2dtrTables);
208 m_pvtwTable.serializeOp(serializer);
209 m_pvcdoTable.serializeOp(serializer);
210 m_densityTable.serializeOp(serializer);
211 m_diffCoeffTable.serializeOp(serializer);
212 m_plyvmhTable.serializeOp(serializer);
213 m_rockTable.serializeOp(serializer);
214 m_plmixparTable.serializeOp(serializer);
215 m_shrateTable.serializeOp(serializer);
216 m_stone1exTable.serializeOp(serializer);
217 m_viscrefTable.serializeOp(serializer);
218 m_watdentTable.serializeOp(serializer);
219 serializer.vector(m_pvtwsaltTables);
220 serializer.vector(m_rwgsaltTables);
221 serializer.vector(m_bdensityTables);
222 serializer.vector(m_sdensityTables);
223 serializer.map(m_plymwinjTables);
224 serializer.map(m_skprwatTables);
225 serializer.map(m_skprpolyTables);
226 m_tabdims.serializeOp(serializer);
227 m_regdims.serializeOp(serializer);
228 m_eqldims.serializeOp(serializer);
229 m_aqudims.serializeOp(serializer);
230 serializer(hasImptvd);
231 serializer(hasEnptvd);
232 serializer(hasEqlnum);
233 serializer(hasShrate);
235 oilDenT.serializeOp(serializer);
236 gasDenT.serializeOp(serializer);
237 watDenT.serializeOp(serializer);
238 stcond.serializeOp(serializer);
239 serializer(m_gas_comp_index);
241 serializer(m_salinity);
242 m_tlmixpar.serializeOp(serializer);
243 if (!serializer.isSerializing()) {
244 m_simpleTables = simpleTables;
245 if (split.plyshMax > 0) {
247 for (
const auto& it : split.plyshMap) {
248 container.addTable(it.first, it.second);
250 m_simpleTables.insert(std::make_pair(
"PLYSHLOG", container));
252 if (split.rockMax > 0) {
253 TableContainer container(split.rockMax);
254 for (
const auto& it : split.rockMap) {
255 container.addTable(it.first, it.second);
257 m_simpleTables.insert(std::make_pair(
"ROCKTAB", container));
263 TableContainer& forceGetTables(
const std::string& tableName ,
size_t numTables);
265 void complainAboutAmbiguousKeyword(
const Deck& deck,
const std::string& keywordName);
267 void addTables(
const std::string& tableName ,
size_t numTables);
268 void initSimpleTables(
const Deck& deck);
269 void initRTempTables(
const Deck& deck);
270 void initDims(
const Deck& deck);
271 void initRocktabTables(
const Deck& deck);
272 void initGasvisctTables(
const Deck& deck);
274 void initPlymaxTables(
const Deck& deck);
275 void initPlyrockTables(
const Deck& deck);
276 void initPlyshlogTables(
const Deck& deck);
278 void initPlymwinjTables(
const Deck& deck);
279 void initSkprwatTables(
const Deck& deck);
280 void initSkprpolyTables(
const Deck& deck);
284 template <
class TableType>
285 void initRockTables(
const Deck& deck,
const std::string& keywordName, std::vector<TableType>& rocktable );
287 template <
class TableType>
288 void initPvtwsaltTables(
const Deck& deck, std::vector<TableType>& pvtwtables );
290 template <
class TableType>
291 void initRwgsaltTables(
const Deck& deck, std::vector<TableType>& rwgtables );
293 template <
class TableType>
294 void initBrineTables(
const Deck& deck, std::vector<TableType>& brinetables );
296 void initSolventTables(
const Deck& deck, std::vector<SolventDensityTable>& solventtables);
301 template <
class TableType>
302 void initSimpleTableContainerWithJFunc(
const Deck& deck,
303 const std::string& keywordName,
304 const std::string& tableName,
307 template <
class TableType>
308 void initSimpleTableContainer(
const Deck& deck,
309 const std::string& keywordName,
310 const std::string& tableName,
313 template <
class TableType>
314 void initSimpleTableContainer(
const Deck& deck,
315 const std::string& keywordName,
318 template <
class TableType>
319 void initSimpleTableContainerWithJFunc(
const Deck& deck,
320 const std::string& keywordName,
323 template <
class TableType>
324 void initSimpleTable(
const Deck& deck,
325 const std::string& keywordName,
326 std::vector<TableType>& tableVector);
328 template <
class TableType>
329 void initFullTables(
const Deck& deck,
330 const std::string& keywordName,
331 std::vector<TableType>& tableVector);
333 void checkPVTOMonotonicity(
const Deck& deck)
const;
335 void logPVTOMonotonicityFailure(
const Deck& deck,
336 const std::size_t tableID,
337 const std::vector<PvtoTable::FlippedFVF>& flipped_Bo)
const;
339 std::map<std::string , TableContainer> m_simpleTables;
340 std::vector<PvtgTable> m_pvtgTables;
341 std::vector<PvtgwTable> m_pvtgwTables;
342 std::vector<PvtgwoTable> m_pvtgwoTables;
343 std::vector<PvtoTable> m_pvtoTables;
344 std::vector<PvtsolTable> m_pvtsolTables;
345 std::vector<Rock2dTable> m_rock2dTables;
346 std::vector<Rock2dtrTable> m_rock2dtrTables;
347 PvtwTable m_pvtwTable;
348 PvcdoTable m_pvcdoTable;
349 DensityTable m_densityTable;
350 DiffCoeffTable m_diffCoeffTable;
351 PlyvmhTable m_plyvmhTable;
352 RockTable m_rockTable;
353 PlmixparTable m_plmixparTable;
354 ShrateTable m_shrateTable;
355 Stone1exTable m_stone1exTable;
356 ViscrefTable m_viscrefTable;
357 WatdentTable m_watdentTable;
358 std::vector<PvtwsaltTable> m_pvtwsaltTables;
359 std::vector<RwgsaltTable> m_rwgsaltTables;
360 std::vector<BrineDensityTable> m_bdensityTables;
361 std::vector<SolventDensityTable> m_sdensityTables;
362 std::map<int, PlymwinjTable> m_plymwinjTables;
363 std::map<int, SkprwatTable> m_skprwatTables;
364 std::map<int, SkprpolyTable> m_skprpolyTables;
372 bool hasImptvd =
false;
373 bool hasEnptvd =
false;
374 bool hasEqlnum =
false;
375 bool hasShrate =
false;
376 std::optional<JFunc> jfunc;
382 std::size_t m_gas_comp_index = 77;
386 struct SplitSimpleTables {
389 std::map<size_t, std::shared_ptr<PlyshlogTable>> plyshMap;
390 std::map<size_t, std::shared_ptr<RocktabTable>> rockMap;
393 SplitSimpleTables splitSimpleTable(std::map<std::string,TableContainer>& simpleTables);
Definition: Aqudims.hpp:34
Definition: Eqldims.hpp:32
Definition: Regdims.hpp:36
Definition: Serializer.hpp:38
Definition: TLMixpar.hpp:55
Definition: Tabdims.hpp:36
Definition: TableContainer.hpp:31
Definition: TableManager.hpp:64
bool useEqlnum() const
deck has keyword "EQLNUM" — Equilibriation region numbers
bool useJFunc() const
deck has keyword "JFUNC" — Use Leverett's J Function for capillary pressure
bool useEnptvd() const
deck has keyword "ENPTVD" — Saturation end-point versus depth tables
bool useImptvd() const
deck has keyword "IMPTVD" — Imbition end-point versus depth tables
bool useShrate() const
deck has keyword "SHRATE"
This class implements a small container which holds the transmissibility mulitpliers for all the face...
Definition: Exceptions.hpp:29
Definition: FlatTable.hpp:45
Definition: FlatTable.hpp:91
Definition: FlatTable.hpp:218
Definition: FlatTable.hpp:252
Definition: FlatTable.hpp:193
Definition: FlatTable.hpp:128
Definition: FlatTable.hpp:156
Definition: FlatTable.hpp:277
Definition: StandardCond.hpp:24
Definition: FlatTable.hpp:302
Definition: FlatTable.hpp:358
Definition: FlatTable.hpp:389