21#ifndef OPM_TABLE_MANAGER_HPP
22#define OPM_TABLE_MANAGER_HPP
28#include <opm/input/eclipse/EclipseState/Tables/DenT.hpp>
29#include <opm/input/eclipse/EclipseState/Tables/JouleThomson.hpp>
30#include <opm/input/eclipse/EclipseState/Tables/PvtgTable.hpp>
31#include <opm/input/eclipse/EclipseState/Tables/PvtgwTable.hpp>
32#include <opm/input/eclipse/EclipseState/Tables/PvtgwoTable.hpp>
33#include <opm/input/eclipse/EclipseState/Tables/PvtoTable.hpp>
34#include <opm/input/eclipse/EclipseState/Tables/PvtsolTable.hpp>
35#include <opm/input/eclipse/EclipseState/Tables/RocktabTable.hpp>
36#include <opm/input/eclipse/EclipseState/Tables/Rock2dTable.hpp>
37#include <opm/input/eclipse/EclipseState/Tables/Rock2dtrTable.hpp>
38#include <opm/input/eclipse/EclipseState/Tables/PlyshlogTable.hpp>
39#include <opm/input/eclipse/EclipseState/Tables/PvtwsaltTable.hpp>
40#include <opm/input/eclipse/EclipseState/Tables/RwgsaltTable.hpp>
41#include <opm/input/eclipse/EclipseState/Tables/BrineDensityTable.hpp>
42#include <opm/input/eclipse/EclipseState/Tables/SolventDensityTable.hpp>
43#include <opm/input/eclipse/EclipseState/Tables/StandardCond.hpp>
44#include <opm/input/eclipse/EclipseState/Tables/FlatTable.hpp>
45#include <opm/input/eclipse/EclipseState/Tables/SorwmisTable.hpp>
46#include <opm/input/eclipse/EclipseState/Tables/SgcwmisTable.hpp>
47#include <opm/input/eclipse/EclipseState/Tables/MiscTable.hpp>
48#include <opm/input/eclipse/EclipseState/Tables/PmiscTable.hpp>
49#include <opm/input/eclipse/EclipseState/Tables/MsfnTable.hpp>
50#include <opm/input/eclipse/EclipseState/Tables/JFunc.hpp>
51#include <opm/input/eclipse/EclipseState/Tables/Tabdims.hpp>
52#include <opm/input/eclipse/EclipseState/Tables/TableContainer.hpp>
53#include <opm/input/eclipse/EclipseState/Tables/Aqudims.hpp>
54#include <opm/input/eclipse/EclipseState/Tables/PlymwinjTable.hpp>
55#include <opm/input/eclipse/EclipseState/Tables/SkprwatTable.hpp>
56#include <opm/input/eclipse/EclipseState/Tables/SkprpolyTable.hpp>
57#include <opm/input/eclipse/EclipseState/Tables/Eqldims.hpp>
58#include <opm/input/eclipse/EclipseState/Tables/Regdims.hpp>
59#include <opm/input/eclipse/EclipseState/Tables/TLMixpar.hpp>
72 const TableContainer& getTables(
const std::string& tableName )
const;
73 const TableContainer& operator[](
const std::string& tableName)
const;
74 bool hasTables(
const std::string& tableName )
const;
76 const Tabdims& getTabdims()
const;
77 const Eqldims& getEqldims()
const;
78 const Aqudims& getAqudims()
const;
79 const Regdims& getRegdims()
const;
84 size_t numFIPRegions()
const;
135 const JFunc& getJFunc()
const;
137 const std::vector<PvtgTable>& getPvtgTables()
const;
138 const std::vector<PvtgwTable>& getPvtgwTables()
const;
139 const std::vector<PvtgwoTable>& getPvtgwoTables()
const;
140 const std::vector<PvtoTable>& getPvtoTables()
const;
141 const std::vector<PvtsolTable>& getPvtsolTables()
const;
142 const std::vector<Rock2dTable>& getRock2dTables()
const;
143 const std::vector<Rock2dtrTable>& getRock2dtrTables()
const;
147 const DenT& WatDenT()
const;
148 const DenT& GasDenT()
const;
149 const DenT& OilDenT()
const;
154 std::size_t gas_comp_index()
const;
156 const std::vector<PvtwsaltTable>& getPvtwSaltTables()
const;
157 const std::vector<RwgsaltTable>& getRwgSaltTables()
const;
158 const std::vector<BrineDensityTable>& getBrineDensityTables()
const;
159 const std::vector<SolventDensityTable>& getSolventDensityTables()
const;
173 const std::map<int, PlymwinjTable>& getPlymwinjTables()
const;
174 const std::map<int, SkprwatTable>& getSkprwatTables()
const;
175 const std::map<int, SkprpolyTable>& getSkprpolyTables()
const;
176 const std::map<std::string, TableContainer>& getSimpleTables()
const;
193 double rtemp()
const;
195 double salinity()
const;
199 template<
class Serializer>
202 auto simpleTables = m_simpleTables;
203 auto split = splitSimpleTable(simpleTables);
204 serializer(simpleTables);
205 serializer(split.plyshMax);
206 serializer(split.plyshMap);
207 serializer(split.rockMax);
208 serializer(split.rockMap);
209 serializer(m_pvtgTables);
210 serializer(m_pvtgwTables);
211 serializer(m_pvtgwoTables);
212 serializer(m_pvtoTables);
213 serializer(m_pvtsolTables);
214 serializer(m_rock2dTables);
215 serializer(m_rock2dtrTables);
216 serializer(m_pvtwTable);
217 serializer(m_pvcdoTable);
218 serializer(m_densityTable);
219 serializer(m_diffCoeffTable);
220 serializer(m_plyvmhTable);
221 serializer(m_rockTable);
222 serializer(m_plmixparTable);
223 serializer(m_shrateTable);
224 serializer(m_stone1exTable);
225 serializer(m_viscrefTable);
226 serializer(m_watdentTable);
227 serializer(m_sgofletTable);
228 serializer(m_swofletTable);
229 serializer(m_pvtwsaltTables);
230 serializer(m_rwgsaltTables);
231 serializer(m_bdensityTables);
232 serializer(m_sdensityTables);
233 serializer(m_plymwinjTables);
234 serializer(m_skprwatTables);
235 serializer(m_skprpolyTables);
236 serializer(m_tabdims);
237 serializer(m_regdims);
238 serializer(m_eqldims);
239 serializer(m_aqudims);
240 serializer(hasImptvd);
241 serializer(hasEnptvd);
242 serializer(hasEqlnum);
243 serializer(hasShrate);
252 serializer(m_gas_comp_index);
254 serializer(m_salinity);
255 serializer(m_tlmixpar);
257 m_simpleTables = simpleTables;
258 if (split.plyshMax > 0) {
260 for (
const auto& it : split.plyshMap) {
261 container.addTable(it.first, it.second);
263 m_simpleTables.insert(std::make_pair(
"PLYSHLOG", container));
265 if (split.rockMax > 0) {
266 TableContainer container(split.rockMax);
267 for (
const auto& it : split.rockMap) {
268 container.addTable(it.first, it.second);
270 m_simpleTables.insert(std::make_pair(
"ROCKTAB", container));
276 TableContainer& forceGetTables(
const std::string& tableName ,
size_t numTables);
278 void complainAboutAmbiguousKeyword(
const Deck& deck,
const std::string& keywordName);
280 void addTables(
const std::string& tableName ,
size_t numTables);
281 void initSimpleTables(
const Deck& deck);
282 void initRTempTables(
const Deck& deck);
283 void initDims(
const Deck& deck);
284 void initRocktabTables(
const Deck& deck);
285 void initGasvisctTables(
const Deck& deck);
287 void initPlymaxTables(
const Deck& deck);
288 void initPlyrockTables(
const Deck& deck);
289 void initPlyshlogTables(
const Deck& deck);
291 void initPlymwinjTables(
const Deck& deck);
292 void initSkprwatTables(
const Deck& deck);
293 void initSkprpolyTables(
const Deck& deck);
297 template <
class TableType>
298 void initRockTables(
const Deck& deck,
const std::string& keywordName, std::vector<TableType>& rocktable );
300 template <
class TableType>
301 void initPvtwsaltTables(
const Deck& deck, std::vector<TableType>& pvtwtables );
303 template <
class TableType>
304 void initRwgsaltTables(
const Deck& deck, std::vector<TableType>& rwgtables );
306 template <
class TableType>
307 void initBrineTables(
const Deck& deck, std::vector<TableType>& brinetables );
309 void initSolventTables(
const Deck& deck, std::vector<SolventDensityTable>& solventtables);
314 template <
class TableType>
315 void initSimpleTableContainerWithJFunc(
const Deck& deck,
316 const std::string& keywordName,
317 const std::string& tableName,
320 template <
class TableType>
321 void initSimpleTableContainer(
const Deck& deck,
322 const std::string& keywordName,
323 const std::string& tableName,
326 template <
class TableType>
327 void initSimpleTableContainer(
const Deck& deck,
328 const std::string& keywordName,
331 template <
class TableType>
332 void initSimpleTableContainerWithJFunc(
const Deck& deck,
333 const std::string& keywordName,
336 template <
class TableType>
337 void initSimpleTable(
const Deck& deck,
338 const std::string& keywordName,
339 std::vector<TableType>& tableVector);
341 template <
class TableType>
342 void initFullTables(
const Deck& deck,
343 const std::string& keywordName,
344 std::vector<TableType>& tableVector);
346 void checkPVTOMonotonicity(
const Deck& deck)
const;
348 void logPVTOMonotonicityFailure(
const Deck& deck,
349 const std::size_t tableID,
350 const std::vector<PvtoTable::FlippedFVF>& flipped_Bo)
const;
352 std::map<std::string , TableContainer> m_simpleTables;
353 std::vector<PvtgTable> m_pvtgTables;
354 std::vector<PvtgwTable> m_pvtgwTables;
355 std::vector<PvtgwoTable> m_pvtgwoTables;
356 std::vector<PvtoTable> m_pvtoTables;
357 std::vector<PvtsolTable> m_pvtsolTables;
358 std::vector<Rock2dTable> m_rock2dTables;
359 std::vector<Rock2dtrTable> m_rock2dtrTables;
360 PvtwTable m_pvtwTable;
361 PvcdoTable m_pvcdoTable;
362 DensityTable m_densityTable;
363 DiffCoeffTable m_diffCoeffTable;
364 PlyvmhTable m_plyvmhTable;
365 RockTable m_rockTable;
366 PlmixparTable m_plmixparTable;
367 ShrateTable m_shrateTable;
368 Stone1exTable m_stone1exTable;
369 ViscrefTable m_viscrefTable;
370 WatdentTable m_watdentTable;
371 SgofletTable m_sgofletTable;
372 SwofletTable m_swofletTable;
373 std::vector<PvtwsaltTable> m_pvtwsaltTables;
374 std::vector<RwgsaltTable> m_rwgsaltTables;
375 std::vector<BrineDensityTable> m_bdensityTables;
376 std::vector<SolventDensityTable> m_sdensityTables;
377 std::map<int, PlymwinjTable> m_plymwinjTables;
378 std::map<int, SkprwatTable> m_skprwatTables;
379 std::map<int, SkprpolyTable> m_skprpolyTables;
387 bool hasImptvd =
false;
388 bool hasEnptvd =
false;
389 bool hasEqlnum =
false;
390 bool hasShrate =
false;
391 std::optional<JFunc> jfunc;
400 std::size_t m_gas_comp_index = 77;
401 double m_rtemp {288.7056};
402 double m_salinity {0.0};
404 struct SplitSimpleTables {
407 std::map<size_t, std::shared_ptr<PlyshlogTable>> plyshMap;
408 std::map<size_t, std::shared_ptr<RocktabTable>> rockMap;
411 SplitSimpleTables splitSimpleTable(std::map<std::string,TableContainer>& simpleTables);
Definition: Aqudims.hpp:34
Definition: Eqldims.hpp:32
Definition: JouleThomson.hpp:30
Definition: Regdims.hpp:36
Class for (de-)serializing.
Definition: Serializer.hpp:75
bool isSerializing() const
Returns true if we are currently doing a serialization operation.
Definition: Serializer.hpp:174
Definition: TLMixpar.hpp:55
Definition: Tabdims.hpp:36
Definition: TableContainer.hpp:31
Definition: TableManager.hpp:65
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:130
Definition: FlatTable.hpp:185
Definition: FlatTable.hpp:321
Definition: FlatTable.hpp:355
Definition: FlatTable.hpp:296
Definition: FlatTable.hpp:223
Definition: FlatTable.hpp:259
Definition: FlatTable.hpp:575
Definition: FlatTable.hpp:380
Definition: StandardCond.hpp:24
Definition: FlatTable.hpp:405
Definition: FlatTable.hpp:565
Definition: FlatTable.hpp:461
Definition: FlatTable.hpp:492