OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-23 16:00:01 -06:00
Code Issues Packages Projects Releases Wiki Activity
e8bf79d774
opm-simulators/ebos/ecloutputblackoilmodule.hh

2183 lines
114 KiB
C++
Raw Normal View History

Added corrections to cancel unintended changes to code
2019-04-23 01:56:06 -05:00
// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
re-add the emacs and vim modelines to all source files this makes things easier and IMHO these two lines do not cause any disturbance.
2015-06-18 06:43:59 -05:00
// vi: set et ts=4 sw=4 sts=4:
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
/*
This file is part of the Open Porous Media project (OPM).
fix
2019-10-30 11:51:49 -05:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
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 2 of the License, or
(at your option) any later version.
fix
2019-10-30 11:51:49 -05:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
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.
fix
2019-10-30 11:51:49 -05:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
You should have received a copy of the GNU General Public License
along with OPM. If not, see <http://www.gnu.org/licenses/>.
fix
2019-10-30 11:51:49 -05:00
clean up the licensing preable of source files the in-file lists of authors has been removed in favor of a global list of authors in the LICENSE file. this is done because (a) maintaining a list of authors at the beginning of a file is a major pain in the a**, (b) the list of authors was not accurate in about 85% of all cases where more than one person was involved and (c) this list is not legally binding in any way (the copyright is at the person who authored a given change, if these lists had any legal relevance, one could "aquire" the copyright of the module by forking it and removing the lists...) the only exception of this is the eWoms fork of dune-istl's solvers.hh file. This is beneficial because the authors of that file do not appear in the global list. Further, carrying the fork of that file is required because we would like to use a reasonable convergence criterion for the linear solver. (the solvers from dune-istl do neither support user-defined convergence criteria not do the developers want support for it. (my patch was rejected a few years ago.))
2016-03-14 07:21:47 -05:00
Consult the COPYING file in the top-level source directory of this
module for the precise wording of the license and the list of
copyright holders.
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
*/
/*!
* \file
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
* \copydoc Opm::EclOutputBlackOilModule
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
*/
#ifndef EWOMS_ECL_OUTPUT_BLACK_OIL_MODULE_HH
#define EWOMS_ECL_OUTPUT_BLACK_OIL_MODULE_HH
changed: ewoms/models/blackoil -> opm/models/blackoil
2019-09-16 04:22:14 -05:00
#include <opm/models/blackoil/blackoilproperties.hh>
make all headers fully autonomous again
2018-03-12 08:31:36 -05:00
changed: ewoms/common -> opm/models/utils
2019-09-16 03:58:20 -05:00
#include <opm/models/utils/propertysystem.hh>
#include <opm/models/utils/parametersystem.hh>
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
adapt to the move of code from opm-common to opm-material
2018-02-01 07:40:01 -06:00
#include <opm/material/common/Valgrind.hpp>
adjustments for eclipse support in opm-common
2018-03-05 06:00:51 -06:00
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
#include <opm/parser/eclipse/Units/Units.hpp>
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
#include <opm/parser/eclipse/EclipseState/SummaryConfig/SummaryConfig.hpp>
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
#include <opm/output/data/Cells.hpp>
#include <opm/output/eclipse/EclipseIO.hpp>
make all headers fully autonomous again
2018-03-12 08:31:36 -05:00
#include <opm/common/OpmLog/OpmLog.hpp>
fix most pedantic compiler warnings in the basic infrastructure i.e., using clang 3.8 to compile the test suite with the following flags: ``` -Weverything -Wno-documentation -Wno-documentation-unknown-command -Wno-c++98-compat -Wno-c++98-compat-pedantic -Wno-undef -Wno-padded -Wno-global-constructors -Wno-exit-time-destructors -Wno-weak-vtables -Wno-float-equal ``` should not produce any warnings anymore. In my opinion the only flag which would produce beneficial warnings is -Wdocumentation. This has not been fixed in this patch because writing documentation is left for another day (or, more likely, year). note that this patch consists of a heavy dose of the OPM_UNUSED macro and plenty of static_casts (to fix signedness issues). Fixing the singedness issues were quite a nightmare and the fact that the Dune API is quite inconsistent in that regard was not exactly helpful. :/ Finally this patch includes quite a few formatting changes (e.g., all occurences of 'T &t' should be changed to `T& t`) and some fixes for minor issues which I've found during the excercise. I've made sure that all unit tests the test suite still pass successfully and I've made sure that flow_ebos still works for Norne and that it did not regress w.r.t. performance. (Note that this patch does not fix compiler warnings triggered `ebos` and `flow_ebos` but only those caused by the basic infrastructure or the unit tests.) v2: fix the warnings that occur if the dune-localfunctions module is not available. thanks to [at]atgeirr for testing. v3: fix dune 2.3 build issue
2016-11-07 08:14:07 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
#include <dune/common/fvector.hh>
#include <type_traits>
do explicit put properties into the the Ewoms::Properties namespace anymore instead, do it implicitly by using the BEGIN_PROPERTIES and END_PROPERTIES macros.
2018-06-14 09:06:05 -05:00
BEGIN_PROPERTIES
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
// create new type tag for the Ecl-output
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
NEW_TYPE_TAG(EclOutputBlackOil);
ebos: add possibility to disable the FIP output even if requested in the deck this significantly reduces the noise level of debug files if the FIP values are not of interest.
2018-07-26 05:31:32 -05:00
NEW_PROP_TAG(ForceDisableFluidInPlaceOutput);
SET_BOOL_PROP(EclOutputBlackOil, ForceDisableFluidInPlaceOutput, false);
do explicit put properties into the the Ewoms::Properties namespace anymore instead, do it implicitly by using the BEGIN_PROPERTIES and END_PROPERTIES macros.
2018-06-14 09:06:05 -05:00
END_PROPERTIES
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
namespace Opm {
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
// forward declaration
template <class TypeTag>
class EcfvDiscretization;
/*!
doxygen: add a separate group for the ECL black-oil simulator application note: the doxygen groups are quite a bit behind the curve and should be overhauled soon. (e.g. now there's not only the vertex centered finite volume space discretization anymore...)
2014-12-22 12:19:03 -06:00
* \ingroup EclBlackOilSimulator
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
*
* \brief Output module for the results black oil model writing in
* ECL binary format.
*/
template <class TypeTag>
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
class EclOutputBlackOilModule
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
{
typedef typename GET_PROP_TYPE(TypeTag, Simulator) Simulator;
typedef typename GET_PROP_TYPE(TypeTag, Discretization) Discretization;
typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
make ebos AD aware
2015-05-21 09:18:45 -05:00
typedef typename GET_PROP_TYPE(TypeTag, Evaluation) Evaluation;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
typedef typename GET_PROP_TYPE(TypeTag, ElementContext) ElementContext;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
typedef typename GET_PROP_TYPE(TypeTag, MaterialLaw) MaterialLaw;
typedef typename GET_PROP_TYPE(TypeTag, MaterialLawParams) MaterialLawParams;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
typedef typename GET_PROP_TYPE(TypeTag, FluidSystem) FluidSystem;
Fix parallel restart and pav in water zone
2018-01-16 08:44:01 -06:00
typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<0>::Iterator ElementIterator;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
enum { numPhases = FluidSystem::numPhases };
enum { oilPhaseIdx = FluidSystem::oilPhaseIdx };
enum { gasPhaseIdx = FluidSystem::gasPhaseIdx };
enum { waterPhaseIdx = FluidSystem::waterPhaseIdx };
enum { gasCompIdx = FluidSystem::gasCompIdx };
adapt to the new PVT API of the black oil fluid system
2015-02-03 06:58:18 -06:00
enum { oilCompIdx = FluidSystem::oilCompIdx };
ebos: only write out temperature in thermal runs
2018-11-16 07:51:07 -06:00
enum { enableEnergy = GET_PROP_VALUE(TypeTag, EnableEnergy) };
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
typedef std::vector<Scalar> ScalarBuffer;
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
typedef std::vector<std::string> StringBuffer;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
struct FipDataType
{
enum FipId
{
Clean up PR review
2018-01-16 07:24:21 -06:00
WaterInPlace = 0, //WIP
OilInPlace = 1, //OIP
GasInPlace = 2, //GIP
OilInPlaceInLiquidPhase = 3, //OIPL
OilInPlaceInGasPhase = 4, //OIPG
GasInPlaceInLiquidPhase = 5, //GIPL
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
GasInPlaceInGasPhase = 6, //GIPG
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
PoreVolume = 7, //PV
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
};
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
static const int numFipValues = PoreVolume + 1 ;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
};
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
struct WellProdDataType
{
enum WPId
{
WellLocationi = 0, //WLi
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
WellLocationj = 1, //WLj
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
OilRate = 2, //OR
WaterRate = 3, //WR
GasRate = 4, //GR
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
FluidResVol = 5, //FRV
WaterCut = 6, //WC
GasOilRatio = 7, //GOR
WatGasRatio = 8, //WGR
BHP = 9, //BHP
THP = 10, //THP
SteadyStatePI = 11, //SteadyStatePI
WellName = 0, //WName
CTRLMode = 1, //CTRL
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
};
static const int numWPValues = 12;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
static const int numWPNames = 2;
};
struct WellInjDataType
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
enum WIId
{
WellLocationi = 0, //WLi
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
WellLocationj = 1, //WLj
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
OilRate = 2, //OR
WaterRate = 3, //WR
GasRate = 4, //GR
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
FluidResVol = 5, //FRV
BHP = 6, //BHP
THP = 7, //THP
SteadyStateII = 8, //SteadyStateII
WellName = 0, //WName
CTRLModeOil = 1, //CTRLo
CTRLModeWat = 2, //CTRLw
CTRLModeGas = 3, //CTRLg
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
};
static const int numWIValues = 9;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
static const int numWINames = 4;
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
};
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
struct WellCumDataType
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
enum WCId
{
WellLocationi = 0, //WLi
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
WellLocationj = 1, //WLj
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
OilProd = 2, //OP
WaterProd = 3, //WP
GasProd = 4, //GP
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
FluidResVolProd = 5, //FRVP
OilInj = 6, //OI
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
WaterInj = 7, //WI
GasInj = 8, //GI
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
FluidResVolInj = 9, //FRVI
WellName = 0, //WName
WellType = 1, //WType
WellCTRL = 2, //WCTRL
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
};
static const int numWCValues = 10;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
static const int numWCNames = 3;
added groups to reports
2019-10-30 08:09:51 -05:00
};
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
public:
Minor clean-up based on review
2018-02-01 06:16:37 -06:00
template<class CollectDataToIORankType>
EclOutputBlackOilModule(const Simulator& simulator, const CollectDataToIORankType& collectToIORank)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
: simulator_(simulator)
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
{
createLocalFipnum_();
Minor clean-up based on review
2018-02-01 06:16:37 -06:00
// Summary output is for all steps
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::SummaryConfig summaryConfig = simulator_.vanguard().summaryConfig();
Minor clean-up based on review
2018-02-01 06:16:37 -06:00
// Initialize block output
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
for (const auto& node: summaryConfig) {
adjust for SummaryNode -> SummaryConfigNode rename
2020-03-13 09:56:33 -05:00
if (node.category() == SummaryConfigNode::Category::Block) {
Blackoil Output: Chase API Change in SummaryNode Class In particular, the .type() function is renamed to .category(), and it no longer returns a LibECL type. Similarly, the .num() function has been renamed to .number().
2019-09-27 08:56:55 -05:00
if (collectToIORank.isGlobalIdxOnThisRank(node.number() - 1)) {
std::pair<std::string, int> key = std::make_pair(node.keyword(), node.number());
Split the PackUnPack class into PackUnPackBlockData, PackUnPackCellData and PackUnPackWellData
2018-02-09 05:57:34 -06:00
blockData_[key] = 0.0;
Minor clean-up based on review
2018-02-01 06:16:37 -06:00
}
}
}
ebos: add possibility to disable the FIP output even if requested in the deck this significantly reduces the noise level of debug files if the FIP values are not of interest.
2018-07-26 05:31:32 -05:00
forceDisableFipOutput_ = EWOMS_GET_PARAM(TypeTag, bool, ForceDisableFluidInPlaceOutput);
}
/*!
* \brief Register all run-time parameters for the Vtk output module.
*/
static void registerParameters()
{
EWOMS_REGISTER_PARAM(TypeTag, bool, ForceDisableFluidInPlaceOutput,
"Do not print fluid-in-place values after each report step even if requested by the deck.");
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
/*!
* \brief Allocate memory for the scalar fields we would like to
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
* write to ECL output files
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
*/
Always allocate buffers when restarting
2019-10-18 02:07:06 -05:00
void allocBuffers(unsigned bufferSize, unsigned reportStepNum, const bool substep, const bool log, const bool isRestart)
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
{
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
if (!std::is_same<Discretization, Opm::EcfvDiscretization<TypeTag> >::value)
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
return;
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
// Summary output is for all steps
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::SummaryConfig summaryConfig = simulator_.vanguard().summaryConfig();
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
// Only output RESTART_AUXILIARY asked for by the user.
RestartConfig has moved from EclipseState to Schedule
2020-02-17 01:39:21 -06:00
const Opm::RestartConfig& restartConfig = simulator_.vanguard().schedule().restart();
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
std::map<std::string, int> rstKeywords = restartConfig.getRestartKeywords(reportStepNum);
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
for (auto& [keyword, should_write] : rstKeywords) {
if (this->isOutputCreationDirective_(keyword)) {
// 'BASIC', 'FREQ' and similar. Don't attempt to create
// cell-based output for these keywords and don't warn about
// not being able to create such cell-based result vectors.
should_write = 0;
}
else {
should_write = restartConfig.getKeyword(keyword, reportStepNum);
}
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
}
outputFipRestart_ = false;
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
computeFip_ = false;
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
// Fluid in place
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
if (!substep || summaryConfig.require3DField(fipEnumToString_(i))) {
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (rstKeywords["FIP"] > 0) {
rstKeywords["FIP"] = 0;
outputFipRestart_ = true;
}
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
fip_[i].resize(bufferSize, 0.0);
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
computeFip_ = true;
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
else
fip_[i].clear();
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
}
if (!substep || summaryConfig.hasKeyword("FPR") || summaryConfig.hasKeyword("FPRP") || summaryConfig.hasKeyword("RPR")) {
Some more minor clean-up
2018-01-29 05:23:23 -06:00
fip_[FipDataType::PoreVolume].resize(bufferSize, 0.0);
Minor clean-up after review
2018-01-29 01:43:42 -06:00
hydrocarbonPoreVolume_.resize(bufferSize, 0.0);
pressureTimesPoreVolume_.resize(bufferSize, 0.0);
pressureTimesHydrocarbonVolume_.resize(bufferSize, 0.0);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
}
else {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
hydrocarbonPoreVolume_.clear();
pressureTimesPoreVolume_.clear();
pressureTimesHydrocarbonVolume_.clear();
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
// Well RFT data
if (!substep) {
Use RFTConfig() object to query for RFT active wells
2019-04-04 00:27:39 -05:00
const auto& schedule = simulator_.vanguard().schedule();
const auto& rft_config = schedule.rftConfig();
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
for (const auto& well: schedule.getWells(reportStepNum)) {
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
// don't bother with wells not on this process
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const auto& defunctWellNames = simulator_.vanguard().defunctWellNames();
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
if (defunctWellNames.find(well.name()) != defunctWellNames.end()) {
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
continue;
}
Use RFTConfig() object to query for RFT active wells
2019-04-04 00:27:39 -05:00
if (!rft_config.active(reportStepNum))
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
continue;
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
for (const auto& connection: well.getConnections()) {
rename completion to connection
2018-06-28 08:49:45 -05:00
const size_t i = size_t(connection.getI());
const size_t j = size_t(connection.getJ());
const size_t k = size_t(connection.getK());
changed: avoid use of EclipseGrid outputting RFT data
2020-02-25 05:33:15 -06:00
const size_t index = simulator_.vanguard().eclState().gridDims().getGlobalIndex(i, j, k);
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
rename completion to connection
2018-06-28 08:49:45 -05:00
oilConnectionPressures_.emplace(std::make_pair(index, 0.0));
waterConnectionSaturations_.emplace(std::make_pair(index, 0.0));
gasConnectionSaturations_.emplace(std::make_pair(index, 0.0));
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
}
}
use `enableEnergy` consistently thanks to [at]totto82 for noticing.
2018-12-17 04:37:07 -06:00
// always allocate memory for temperature
temperature_.resize(bufferSize, 0.0);
Always allocate buffers when restarting
2019-10-18 02:07:06 -05:00
// field data should be allocated
// 1) when we want to restart
// 2) when it is ask for by the user via restartConfig
// 3) when it is not a substep
if (!isRestart && (!restartConfig.getWriteRestartFile(reportStepNum, log) || substep))
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
return;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
// always output saturation of active phases
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx))
continue;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
saturation_[phaseIdx].resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
// and oil pressure
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
oilPressure_.resize(bufferSize, 0.0);
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
rstKeywords["PRES"] = 0;
rstKeywords["PRESSURE"] = 0;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
if (FluidSystem::phaseIsActive(oilPhaseIdx))
rstKeywords["SOIL"] = 0;
if (FluidSystem::phaseIsActive(gasPhaseIdx))
rstKeywords["SGAS"] = 0;
if (FluidSystem::phaseIsActive(waterPhaseIdx))
rstKeywords["SWAT"] = 0;
if (FluidSystem::enableDissolvedGas()) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
rs_.resize(bufferSize, 0.0);
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
rstKeywords["RS"] = 0;
}
if (FluidSystem::enableVaporizedOil()) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
rv_.resize(bufferSize, 0.0);
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
rstKeywords["RV"] = 0;
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (GET_PROP_VALUE(TypeTag, EnableSolvent))
sSol_.resize(bufferSize, 0.0);
if (GET_PROP_VALUE(TypeTag, EnablePolymer))
cPolymer_.resize(bufferSize, 0.0);
Use foam module.
2019-04-26 03:38:37 -05:00
if (GET_PROP_VALUE(TypeTag, EnableFoam))
cFoam_.resize(bufferSize, 0.0);
renamed saltwater-->brine and removed keywords from missingFeatures
2019-12-03 10:48:17 -06:00
if (GET_PROP_VALUE(TypeTag, EnableBrine))
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
cSalt_.resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (simulator_.problem().vapparsActive())
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
soMax_.resize(bufferSize, 0.0);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (simulator_.problem().materialLawManager()->enableHysteresis()) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
pcSwMdcOw_.resize(bufferSize, 0.0);
krnSwMdcOw_.resize(bufferSize, 0.0);
pcSwMdcGo_.resize(bufferSize, 0.0);
krnSwMdcGo_.resize(bufferSize, 0.0);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
if (simulator_.vanguard().eclState().fieldProps().has_double("SWATINIT")) {
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
ppcw_.resize(bufferSize, 0.0);
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
rstKeywords["PPCW"] = 0;
}
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (FluidSystem::enableDissolvedGas() && rstKeywords["RSSAT"] > 0) {
rstKeywords["RSSAT"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
gasDissolutionFactor_.resize(bufferSize, 0.0);
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
}
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (FluidSystem::enableVaporizedOil() && rstKeywords["RVSAT"] > 0) {
rstKeywords["RVSAT"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
oilVaporizationFactor_.resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(waterPhaseIdx) && rstKeywords["BW"] > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["BW"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
invB_[waterPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(oilPhaseIdx) && rstKeywords["BO"] > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["BO"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
invB_[oilPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(gasPhaseIdx) && rstKeywords["BG"] > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["BG"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
invB_[gasPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
if (rstKeywords["DEN"] > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
rstKeywords["DEN"] = 0;
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx))
continue;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
density_[phaseIdx].resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
const bool hasVWAT = (rstKeywords["VISC"] > 0) || (rstKeywords["VWAT"] > 0);
const bool hasVOIL = (rstKeywords["VISC"] > 0) || (rstKeywords["VOIL"] > 0);
const bool hasVGAS = (rstKeywords["VISC"] > 0) || (rstKeywords["VGAS"] > 0);
rstKeywords["VISC"] = 0;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(waterPhaseIdx) && hasVWAT) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["VWAT"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
viscosity_[waterPhaseIdx].resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(oilPhaseIdx) && hasVOIL > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["VOIL"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
viscosity_[oilPhaseIdx].resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(gasPhaseIdx) && hasVGAS > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["VGAS"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
viscosity_[gasPhaseIdx].resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(waterPhaseIdx) && rstKeywords["KRW"] > 0) {
Added corrections to cancel unintended changes to code
2019-04-23 01:56:06 -05:00
rstKeywords["KRW"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
relativePermeability_[waterPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(oilPhaseIdx) && rstKeywords["KRO"] > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["KRO"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
relativePermeability_[oilPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(gasPhaseIdx) && rstKeywords["KRG"] > 0) {
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
rstKeywords["KRG"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
relativePermeability_[gasPhaseIdx].resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (rstKeywords["PBPD"] > 0) {
rstKeywords["PBPD"] = 0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
bubblePointPressure_.resize(bufferSize, 0.0);
dewPointPressure_.resize(bufferSize, 0.0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
// tracers
const int numTracers = simulator_.problem().tracerModel().numTracers();
if (numTracers > 0){
tracerConcentrations_.resize(numTracers);
for (int tracerIdx = 0; tracerIdx < numTracers; ++tracerIdx)
{
tracerConcentrations_[tracerIdx].resize(bufferSize, 0.0);
}
}
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
// ROCKC
if (rstKeywords["ROCKC"] > 0) {
rstKeywords["ROCKC"] = 0;
rockCompPorvMultiplier_.resize(bufferSize, 0.0);
rockCompTransMultiplier_.resize(bufferSize, 0.0);
swMax_.resize(bufferSize, 0.0);
minimumOilPressure_.resize(bufferSize, 0.0);
overburdenPressure_.resize(bufferSize, 0.0);
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
//Warn for any unhandled keyword
if (log) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
for (auto& keyValue: rstKeywords) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (keyValue.second > 0) {
std::string logstring = "Keyword '";
logstring.append(keyValue.first);
logstring.append("' is unhandled for output to file.");
Opm::OpmLog::warning("Unhandled output keyword", logstring);
}
}
}
failedCellsPb_.clear();
failedCellsPd_.clear();
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
// Not supported in flow legacy
if (false)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
saturatedOilFormationVolumeFactor_.resize(bufferSize, 0.0);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (false)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
oilSaturationPressure_.resize(bufferSize, 0.0);
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
/*!
* \brief Modify the internal buffers according to the intensive quanties relevant
* for an element
*/
fix most pedantic compiler warnings in the basic infrastructure i.e., using clang 3.8 to compile the test suite with the following flags: ``` -Weverything -Wno-documentation -Wno-documentation-unknown-command -Wno-c++98-compat -Wno-c++98-compat-pedantic -Wno-undef -Wno-padded -Wno-global-constructors -Wno-exit-time-destructors -Wno-weak-vtables -Wno-float-equal ``` should not produce any warnings anymore. In my opinion the only flag which would produce beneficial warnings is -Wdocumentation. This has not been fixed in this patch because writing documentation is left for another day (or, more likely, year). note that this patch consists of a heavy dose of the OPM_UNUSED macro and plenty of static_casts (to fix signedness issues). Fixing the singedness issues were quite a nightmare and the fact that the Dune API is quite inconsistent in that regard was not exactly helpful. :/ Finally this patch includes quite a few formatting changes (e.g., all occurences of 'T &t' should be changed to `T& t`) and some fixes for minor issues which I've found during the excercise. I've made sure that all unit tests the test suite still pass successfully and I've made sure that flow_ebos still works for Norne and that it did not regress w.r.t. performance. (Note that this patch does not fix compiler warnings triggered `ebos` and `flow_ebos` but only those caused by the basic infrastructure or the unit tests.) v2: fix the warnings that occur if the dune-localfunctions module is not available. thanks to [at]atgeirr for testing. v3: fix dune 2.3 build issue
2016-11-07 08:14:07 -06:00
void processElement(const ElementContext& elemCtx)
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
{
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
if (!std::is_same<Discretization, Opm::EcfvDiscretization<TypeTag> >::value)
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
return;
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
const auto& problem = elemCtx.simulator().problem();
make most indices unsigned (instead of using 'int'.) This triggered quite a few compiler warnings which are also dealt-with by this patch.
2015-11-18 04:54:35 -06:00
for (unsigned dofIdx = 0; dofIdx < elemCtx.numPrimaryDof(/*timeIdx=*/0); ++dofIdx) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
const auto& intQuants = elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0);
const auto& fs = intQuants.fluidState();
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
adapt to the recent blackoil API changes of opm-material
2016-01-04 08:32:55 -06:00
typedef typename std::remove_const<typename std::remove_reference<decltype(fs)>::type>::type FluidState;
make most indices unsigned (instead of using 'int'.) This triggered quite a few compiler warnings which are also dealt-with by this patch.
2015-11-18 04:54:35 -06:00
unsigned globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, /*timeIdx=*/0);
ECL restart file output: rename 'regionIdx' to 'pvtRegionIdx'
2015-12-31 06:20:41 -06:00
unsigned pvtRegionIdx = elemCtx.primaryVars(dofIdx, /*timeIdx=*/0).pvtRegionIndex();
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (saturation_[phaseIdx].size() == 0)
continue;
implement two-phase blackoil simulations
2016-11-15 11:00:48 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
saturation_[phaseIdx][globalDofIdx] = Opm::getValue(fs.saturation(phaseIdx));
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Opm::Valgrind::CheckDefined(saturation_[phaseIdx][globalDofIdx]);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (oilPressure_.size() > 0) {
Enable single-phase runs.
2019-10-09 08:24:23 -05:00
if (FluidSystem::phaseIsActive(oilPhaseIdx)) {
oilPressure_[globalDofIdx] = Opm::getValue(fs.pressure(oilPhaseIdx));
}else{
// put pressure in oil pressure for output
if (FluidSystem::phaseIsActive(waterPhaseIdx)) {
oilPressure_[globalDofIdx] = Opm::getValue(fs.pressure(waterPhaseIdx));
} else {
oilPressure_[globalDofIdx] = Opm::getValue(fs.pressure(gasPhaseIdx));
}
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Opm::Valgrind::CheckDefined(oilPressure_[globalDofIdx]);
Add output of temperature
2018-01-31 04:10:37 -06:00
}
implement two-phase blackoil simulations
2016-11-15 11:00:48 -06:00
use `enableEnergy` consistently thanks to [at]totto82 for noticing.
2018-12-17 04:37:07 -06:00
if (enableEnergy) {
Add output of temperature
2018-01-31 04:10:37 -06:00
temperature_[globalDofIdx] = Opm::getValue(fs.temperature(oilPhaseIdx));
Opm::Valgrind::CheckDefined(temperature_[globalDofIdx]);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (gasDissolutionFactor_.size() > 0) {
move the maximum oil saturation stuff into the problem this also simplifies creating a more proper selection logic for the DR[SV]DT keywords or VAPPARS.
2017-12-19 11:00:14 -06:00
Scalar SoMax = elemCtx.problem().maxOilSaturation(globalDofIdx);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
gasDissolutionFactor_[globalDofIdx] =
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
FluidSystem::template saturatedDissolutionFactor<FluidState, Scalar>(fs, oilPhaseIdx, pvtRegionIdx, SoMax);
adapt to the move of the valgrind client requests into the Opm namespace
2017-02-09 11:25:44 -06:00
Opm::Valgrind::CheckDefined(gasDissolutionFactor_[globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (oilVaporizationFactor_.size() > 0) {
move the maximum oil saturation stuff into the problem this also simplifies creating a more proper selection logic for the DR[SV]DT keywords or VAPPARS.
2017-12-19 11:00:14 -06:00
Scalar SoMax = elemCtx.problem().maxOilSaturation(globalDofIdx);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
oilVaporizationFactor_[globalDofIdx] =
FluidSystem::template saturatedDissolutionFactor<FluidState, Scalar>(fs, gasPhaseIdx, pvtRegionIdx, SoMax);
Opm::Valgrind::CheckDefined(oilVaporizationFactor_[globalDofIdx]);
implement two-phase blackoil simulations
2016-11-15 11:00:48 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (gasFormationVolumeFactor_.size() > 0) {
Revert "rename {gil,gas}FormationVolumeFactor to \1FormationFactor" This reverts commit 3e247b059f20b1d04f47bd832719a99df8b9ac26. After a more thorough investigation, the cannonical name of these quantities turned out to be "* formation volume factor"...
2015-01-26 04:55:37 -06:00
gasFormationVolumeFactor_[globalDofIdx] =
adapt to the recent API changes in opm-material the change for the inverse formation volume factors (instead of the FVFs) is needed to keep the build happy, using the new getRv_()/getRs_() infrastructure might slightly improve performance. (the emphasis is *slightly* because it only changes something if BlackOilFluidState is used for initialization and initialization is not time critical in the first place.)
2016-01-29 04:57:41 -06:00
1.0/FluidSystem::template inverseFormationVolumeFactor<FluidState, Scalar>(fs, gasPhaseIdx, pvtRegionIdx);
adapt to the move of the valgrind client requests into the Opm namespace
2017-02-09 11:25:44 -06:00
Opm::Valgrind::CheckDefined(gasFormationVolumeFactor_[globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (saturatedOilFormationVolumeFactor_.size() > 0) {
saturatedOilFormationVolumeFactor_[globalDofIdx] =
1.0/FluidSystem::template saturatedInverseFormationVolumeFactor<FluidState, Scalar>(fs, oilPhaseIdx, pvtRegionIdx);
Opm::Valgrind::CheckDefined(saturatedOilFormationVolumeFactor_[globalDofIdx]);
}
if (oilSaturationPressure_.size() > 0) {
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
oilSaturationPressure_[globalDofIdx] =
ECL restart file output: rename 'regionIdx' to 'pvtRegionIdx'
2015-12-31 06:20:41 -06:00
FluidSystem::template saturationPressure<FluidState, Scalar>(fs, oilPhaseIdx, pvtRegionIdx);
adapt to the move of the valgrind client requests into the Opm namespace
2017-02-09 11:25:44 -06:00
Opm::Valgrind::CheckDefined(oilSaturationPressure_[globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
if (rs_.size()) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
rs_[globalDofIdx] = Opm::getValue(fs.Rs());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Opm::Valgrind::CheckDefined(rs_[globalDofIdx]);
}
if (rv_.size()) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
rv_[globalDofIdx] = Opm::getValue(fs.Rv());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Opm::Valgrind::CheckDefined(rv_[globalDofIdx]);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (invB_[phaseIdx].size() == 0)
continue;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
invB_[phaseIdx][globalDofIdx] = Opm::getValue(fs.invB(phaseIdx));
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Opm::Valgrind::CheckDefined(invB_[phaseIdx][globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (density_[phaseIdx].size() == 0)
continue;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
density_[phaseIdx][globalDofIdx] = Opm::getValue(fs.density(phaseIdx));
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Opm::Valgrind::CheckDefined(density_[phaseIdx][globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (viscosity_[phaseIdx].size() == 0)
continue;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
viscosity_[phaseIdx][globalDofIdx] = Opm::getValue(fs.viscosity(phaseIdx));
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Opm::Valgrind::CheckDefined(viscosity_[phaseIdx][globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (relativePermeability_[phaseIdx].size() == 0)
continue;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
relativePermeability_[phaseIdx][globalDofIdx] = Opm::getValue(intQuants.relativePermeability(phaseIdx));
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
Opm::Valgrind::CheckDefined(relativePermeability_[phaseIdx][globalDofIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
if (sSol_.size() > 0) {
sSol_[globalDofIdx] = intQuants.solventSaturation().value();
}
if (cPolymer_.size() > 0) {
cPolymer_[globalDofIdx] = intQuants.polymerConcentration().value();
}
Use foam module.
2019-04-26 03:38:37 -05:00
if (cFoam_.size() > 0) {
cFoam_[globalDofIdx] = intQuants.foamConcentration().value();
}
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
if (cSalt_.size() > 0) {
renamed saltwater-->brine and removed keywords from missingFeatures
2019-12-03 10:48:17 -06:00
cSalt_[globalDofIdx] = fs.saltConcentration().value();
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (bubblePointPressure_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
try {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
bubblePointPressure_[globalDofIdx] = Opm::getValue(FluidSystem::bubblePointPressure(fs, intQuants.pvtRegionIndex()));
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
fix a few "unused variable" warnings that occur on more uncommon compilers i.e. on ICC. the compiler is right in producing warnings here, though.
2018-03-03 06:06:58 -06:00
catch (const Opm::NumericalIssue&) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const auto cartesianIdx = elemCtx.simulator().vanguard().grid().globalCell()[globalDofIdx];
failedCellsPb_.push_back(cartesianIdx);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (dewPointPressure_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
try {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
dewPointPressure_[globalDofIdx] = Opm::getValue(FluidSystem::dewPointPressure(fs, intQuants.pvtRegionIndex()));
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
fix a few "unused variable" warnings that occur on more uncommon compilers i.e. on ICC. the compiler is right in producing warnings here, though.
2018-03-03 06:06:58 -06:00
catch (const Opm::NumericalIssue&) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const auto cartesianIdx = elemCtx.simulator().vanguard().grid().globalCell()[globalDofIdx];
failedCellsPd_.push_back(cartesianIdx);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
}
if (soMax_.size() > 0)
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
soMax_[globalDofIdx] =
std::max(Opm::getValue(fs.saturation(oilPhaseIdx)),
problem.maxOilSaturation(globalDofIdx));
if (swMax_.size() > 0)
swMax_[globalDofIdx] =
std::max(Opm::getValue(fs.saturation(waterPhaseIdx)),
problem.maxWaterSaturation(globalDofIdx));
if (minimumOilPressure_.size() > 0)
minimumOilPressure_[globalDofIdx] =
std::min(Opm::getValue(fs.pressure(oilPhaseIdx)),
problem.minOilPressure(globalDofIdx));
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
if (overburdenPressure_.size() > 0)
overburdenPressure_[globalDofIdx] = problem.overburdenPressure(globalDofIdx);
if (rockCompPorvMultiplier_.size() > 0)
rockCompPorvMultiplier_[globalDofIdx] = problem.template rockCompPoroMultiplier<Scalar>(intQuants, globalDofIdx);
if (rockCompTransMultiplier_.size() > 0)
rockCompTransMultiplier_[globalDofIdx] = problem.template rockCompTransMultiplier<Scalar>(intQuants, globalDofIdx);
const auto& matLawManager = problem.materialLawManager();
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (matLawManager->enableHysteresis()) {
if (pcSwMdcOw_.size() > 0 && krnSwMdcOw_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
matLawManager->oilWaterHysteresisParams(
pcSwMdcOw_[globalDofIdx],
krnSwMdcOw_[globalDofIdx],
globalDofIdx);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
if (pcSwMdcGo_.size() > 0 && krnSwMdcGo_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
matLawManager->gasOilHysteresisParams(
pcSwMdcGo_[globalDofIdx],
krnSwMdcGo_[globalDofIdx],
globalDofIdx);
}
}
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
Simulate with only active cells (#2213) Use FieldProps implementation for 3D properties
2020-01-13 08:46:50 -06:00
if (ppcw_.size() > 0) {
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
ppcw_[globalDofIdx] = matLawManager->oilWaterScaledEpsInfoDrainage(globalDofIdx).maxPcow;
Simulate with only active cells (#2213) Use FieldProps implementation for 3D properties
2020-01-13 08:46:50 -06:00
//printf("ppcw_[%d] = %lg\n", globalDofIdx, ppcw_[globalDofIdx]);
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
// hack to make the intial output of rs and rv Ecl compatible.
// For cells with swat == 1 Ecl outputs; rs = rsSat and rv=rvSat, in all but the initial step
// where it outputs rs and rv values calculated by the initialization. To be compatible we overwrite
// rs and rv with the values computed in the initially.
// Volume factors, densities and viscosities need to be recalculated with the updated rs and rv values.
// This can be removed when ebos has 100% controll over output
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (elemCtx.simulator().episodeIndex() < 0 && FluidSystem::phaseIsActive(oilPhaseIdx) && FluidSystem::phaseIsActive(gasPhaseIdx)) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
const auto& fsInitial = problem.initialFluidState(globalDofIdx);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
// use initial rs and rv values
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (rv_.size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
rv_[globalDofIdx] = fsInitial.Rv();
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (rs_.size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
rs_[globalDofIdx] = fsInitial.Rs();
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
// re-compute the volume factors, viscosities and densities if asked for
if (density_[oilPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
density_[oilPhaseIdx][globalDofIdx] = FluidSystem::density(fsInitial,
oilPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (density_[gasPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
density_[gasPhaseIdx][globalDofIdx] = FluidSystem::density(fsInitial,
gasPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (invB_[oilPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
invB_[oilPhaseIdx][globalDofIdx] = FluidSystem::inverseFormationVolumeFactor(fsInitial,
oilPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (invB_[gasPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
invB_[gasPhaseIdx][globalDofIdx] = FluidSystem::inverseFormationVolumeFactor(fsInitial,
gasPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (viscosity_[oilPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
viscosity_[oilPhaseIdx][globalDofIdx] = FluidSystem::viscosity(fsInitial,
oilPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (viscosity_[gasPhaseIdx].size() > 0)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
viscosity_[gasPhaseIdx][globalDofIdx] = FluidSystem::viscosity(fsInitial,
gasPhaseIdx,
intQuants.pvtRegionIndex());
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Minor clean-up after review
2018-01-29 01:43:42 -06:00
// Add fluid in Place values
updateFluidInPlace_(elemCtx, dofIdx);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Split the PackUnPack class into PackUnPackBlockData, PackUnPackCellData and PackUnPackWellData
2018-02-09 05:57:34 -06:00
// Adding block data
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const auto cartesianIdx = elemCtx.simulator().vanguard().grid().globalCell()[globalDofIdx];
for (auto& val: blockData_) {
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
const auto& key = val.first;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
int cartesianIdxBlock = key.second - 1;
if (cartesianIdx == cartesianIdxBlock) {
if (key.first == "BWSAT")
Minor clean-up after review
2018-01-29 01:43:42 -06:00
val.second = Opm::getValue(fs.saturation(waterPhaseIdx));
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
else if (key.first == "BGSAT")
Minor clean-up after review
2018-01-29 01:43:42 -06:00
val.second = Opm::getValue(fs.saturation(gasPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BOSAT")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = 1. - Opm::getValue(fs.saturation(gasPhaseIdx)) - Opm::getValue(fs.saturation(waterPhaseIdx));
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
else if (key.first == "BPR")
Minor clean-up after review
2018-01-29 01:43:42 -06:00
val.second = Opm::getValue(fs.pressure(oilPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BWKR" || key.first == "BKRW")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(intQuants.relativePermeability(waterPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BGKR" || key.first == "BKRG")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(intQuants.relativePermeability(gasPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BOKR" || key.first == "BKRO")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(intQuants.relativePermeability(oilPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BWPC")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(fs.pressure(oilPhaseIdx)) - Opm::getValue(fs.pressure(waterPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BGPC")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(fs.pressure(gasPhaseIdx)) - Opm::getValue(fs.pressure(oilPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BVWAT" || key.first == "BWVIS")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(fs.viscosity(waterPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BVGAS" || key.first == "BGVIS")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(fs.viscosity(gasPhaseIdx));
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
else if (key.first == "BVOIL" || key.first == "BOVIS")
changes to provide Summary data output of additional block variables Provide output data for: BOSAT BWKR BOKR BKRO BGKR BKRG BKRW BWPC BGPC BVWAT BWVIS BVGAS BGVIS BVOIL BOVIS
2019-04-12 09:48:27 -05:00
val.second = Opm::getValue(fs.viscosity(oilPhaseIdx));
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
else {
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
std::string logstring = "Keyword '";
logstring.append(key.first);
logstring.append("' is unhandled for output to file.");
Opm::OpmLog::warning("Unhandled output keyword", logstring);
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
}
}
}
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
// Adding Well RFT data
rename completion to connection
2018-06-28 08:49:45 -05:00
if (oilConnectionPressures_.count(cartesianIdx) > 0) {
oilConnectionPressures_[cartesianIdx] = Opm::getValue(fs.pressure(oilPhaseIdx));
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
rename completion to connection
2018-06-28 08:49:45 -05:00
if (waterConnectionSaturations_.count(cartesianIdx) > 0) {
waterConnectionSaturations_[cartesianIdx] = Opm::getValue(fs.saturation(waterPhaseIdx));
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
rename completion to connection
2018-06-28 08:49:45 -05:00
if (gasConnectionSaturations_.count(cartesianIdx) > 0) {
gasConnectionSaturations_[cartesianIdx] = Opm::getValue(fs.saturation(gasPhaseIdx));
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
// tracers
const auto& tracerModel = simulator_.problem().tracerModel();
if (tracerConcentrations_.size()>0) {
for (int tracerIdx = 0; tracerIdx < tracerModel.numTracers(); tracerIdx++){
if (tracerConcentrations_[tracerIdx].size() == 0)
continue;
tracerConcentrations_[tracerIdx][globalDofIdx] = tracerModel.tracerConcentration(tracerIdx, globalDofIdx);
}
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
}
void outputErrorLog()
{
const size_t maxNumCellsFaillog = 20;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
int pbSize = failedCellsPb_.size(), pdSize = failedCellsPd_.size();
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
std::vector<int> displPb, displPd, recvLenPb, recvLenPd;
const auto& comm = simulator_.gridView().comm();
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (isIORank_()) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
displPb.resize(comm.size()+1, 0);
displPd.resize(comm.size()+1, 0);
recvLenPb.resize(comm.size());
recvLenPd.resize(comm.size());
}
comm.gather(&pbSize, recvLenPb.data(), 1, 0);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
comm.gather(&pdSize, recvLenPd.data(), 1, 0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
std::partial_sum(recvLenPb.begin(), recvLenPb.end(), displPb.begin()+1);
std::partial_sum(recvLenPd.begin(), recvLenPd.end(), displPd.begin()+1);
std::vector<int> globalFailedCellsPb, globalFailedCellsPd;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (isIORank_()) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
globalFailedCellsPb.resize(displPb.back());
globalFailedCellsPd.resize(displPd.back());
}
comm.gatherv(failedCellsPb_.data(), static_cast<int>(failedCellsPb_.size()),
globalFailedCellsPb.data(), recvLenPb.data(),
displPb.data(), 0);
comm.gatherv(failedCellsPd_.data(), static_cast<int>(failedCellsPd_.size()),
globalFailedCellsPd.data(), recvLenPd.data(),
displPd.data(), 0);
std::sort(globalFailedCellsPb.begin(), globalFailedCellsPb.end());
std::sort(globalFailedCellsPd.begin(), globalFailedCellsPd.end());
if (globalFailedCellsPb.size() > 0) {
std::stringstream errlog;
errlog << "Finding the bubble point pressure failed for " << globalFailedCellsPb.size() << " cells [";
errlog << globalFailedCellsPb[0];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const size_t maxElems = std::min(maxNumCellsFaillog, globalFailedCellsPb.size());
for (size_t i = 1; i < maxElems; ++i) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
errlog << ", " << globalFailedCellsPb[i];
}
if (globalFailedCellsPb.size() > maxNumCellsFaillog) {
errlog << ", ...";
}
errlog << "]";
Opm::OpmLog::warning("Bubble point numerical problem", errlog.str());
}
if (globalFailedCellsPd.size() > 0) {
std::stringstream errlog;
errlog << "Finding the dew point pressure failed for " << globalFailedCellsPd.size() << " cells [";
errlog << globalFailedCellsPd[0];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const size_t maxElems = std::min(maxNumCellsFaillog, globalFailedCellsPd.size());
for (size_t i = 1; i < maxElems; ++i) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
errlog << ", " << globalFailedCellsPd[i];
}
if (globalFailedCellsPd.size() > maxNumCellsFaillog) {
errlog << ", ...";
}
errlog << "]";
Opm::OpmLog::warning("Dew point numerical problem", errlog.str());
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
}
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
void addRftDataToWells(Opm::data::Wells& wellDatas, size_t reportStepNum)
make all headers fully autonomous again
2018-03-12 08:31:36 -05:00
{
Use RFTConfig() object to query for RFT active wells
2019-04-04 00:27:39 -05:00
const auto& schedule = simulator_.vanguard().schedule();
const auto& rft_config = schedule.rftConfig();
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
for (const auto& well: schedule.getWells(reportStepNum)) {
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
// don't bother with wells not on this process
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const auto& defunctWellNames = simulator_.vanguard().defunctWellNames();
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
if (defunctWellNames.find(well.name()) != defunctWellNames.end()) {
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
continue;
}
//add data infrastructure for shut wells
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
if (!wellDatas.count(well.name())) {
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
Opm::data::Well wellData;
Use RFTConfig() object to query for RFT active wells
2019-04-04 00:27:39 -05:00
if (!rft_config.active(reportStepNum))
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
continue;
Use RFTConfig() object to query for RFT active wells
2019-04-04 00:27:39 -05:00
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
wellData.connections.resize(well.getConnections().size());
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
size_t count = 0;
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
for (const auto& connection: well.getConnections()) {
rename completion to connection
2018-06-28 08:49:45 -05:00
const size_t i = size_t(connection.getI());
const size_t j = size_t(connection.getJ());
const size_t k = size_t(connection.getK());
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
changed: avoid use of EclipseGrid outputting RFT data
2020-02-25 05:33:15 -06:00
const size_t index = simulator_.vanguard().eclState().gridDims().getGlobalIndex(i, j, k);
rename completion to connection
2018-06-28 08:49:45 -05:00
auto& connectionData = wellData.connections[count];
connectionData.index = index;
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
count++;
}
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
wellDatas.emplace(std::make_pair(well.name(), wellData));
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
Use new well implementation Well2 from opm-common
2019-05-02 05:51:25 -05:00
Opm::data::Well& wellData = wellDatas.at(well.name());
rename completion to connection
2018-06-28 08:49:45 -05:00
for (auto& connectionData: wellData.connections) {
const auto index = connectionData.index;
if (oilConnectionPressures_.count(index) > 0)
connectionData.cell_pressure = oilConnectionPressures_.at(index);
if (waterConnectionSaturations_.count(index) > 0)
connectionData.cell_saturation_water = waterConnectionSaturations_.at(index);
if (gasConnectionSaturations_.count(index) > 0)
connectionData.cell_saturation_gas = gasConnectionSaturations_.at(index);
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
}
rename completion to connection
2018-06-28 08:49:45 -05:00
oilConnectionPressures_.clear();
waterConnectionSaturations_.clear();
gasConnectionSaturations_.clear();
Pass RFT output via data::well With this no cell data is used by the output code in opm-output
2018-02-07 07:30:43 -06:00
}
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
/*!
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
* \brief Move all buffers to data::Solution.
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
*/
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
void assignToSolution(Opm::data::Solution& sol)
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
{
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
if (!std::is_same<Discretization, Opm::EcfvDiscretization<TypeTag>>::value)
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
return;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (oilPressure_.size() > 0) {
sol.insert("PRESSURE", Opm::UnitSystem::measure::pressure, std::move(oilPressure_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ebos: only write out temperature in thermal runs
2018-11-16 07:51:07 -06:00
if (enableEnergy) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
sol.insert("TEMP", Opm::UnitSystem::measure::temperature, std::move(temperature_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ECL output: use the same units as the input deck for ECL output this makes it possible to compare ECL restart and summary files produced by ebos directly with the ones stemming from Eclips. (But be aware that VTK output files from ebos are still all-SI!)
2015-01-25 10:27:08 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(waterPhaseIdx) && saturation_[waterPhaseIdx].size() > 0) {
sol.insert("SWAT", Opm::UnitSystem::measure::identity, std::move(saturation_[waterPhaseIdx]), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (FluidSystem::phaseIsActive(gasPhaseIdx) && saturation_[gasPhaseIdx].size() > 0) {
sol.insert("SGAS", Opm::UnitSystem::measure::identity, std::move(saturation_[gasPhaseIdx]), Opm::data::TargetType::RESTART_SOLUTION);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
if (ppcw_.size() > 0) {
sol.insert ("PPCW", Opm::UnitSystem::measure::pressure, std::move(ppcw_), Opm::data::TargetType::RESTART_SOLUTION);
Store capillary pressure in SOLUTION section
2018-09-04 02:07:18 -05:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (gasDissolutionFactor_.size() > 0) {
sol.insert("RSSAT", Opm::UnitSystem::measure::gas_oil_ratio, std::move(gasDissolutionFactor_), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (oilVaporizationFactor_.size() > 0) {
sol.insert("RVSAT", Opm::UnitSystem::measure::oil_gas_ratio, std::move(oilVaporizationFactor_), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (rs_.size() > 0) {
sol.insert("RS", Opm::UnitSystem::measure::gas_oil_ratio, std::move(rs_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (rv_.size() > 0) {
sol.insert("RV", Opm::UnitSystem::measure::oil_gas_ratio, std::move(rv_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (invB_[waterPhaseIdx].size() > 0) {
sol.insert("1OVERBW", Opm::UnitSystem::measure::water_inverse_formation_volume_factor, std::move(invB_[waterPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (invB_[oilPhaseIdx].size() > 0) {
sol.insert("1OVERBO", Opm::UnitSystem::measure::oil_inverse_formation_volume_factor, std::move(invB_[oilPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (invB_[gasPhaseIdx].size() > 0) {
sol.insert("1OVERBG", Opm::UnitSystem::measure::gas_inverse_formation_volume_factor, std::move(invB_[gasPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (density_[waterPhaseIdx].size() > 0) {
sol.insert("WAT_DEN", Opm::UnitSystem::measure::density, std::move(density_[waterPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (density_[oilPhaseIdx].size() > 0) {
sol.insert("OIL_DEN", Opm::UnitSystem::measure::density, std::move(density_[oilPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (density_[gasPhaseIdx].size() > 0) {
sol.insert("GAS_DEN", Opm::UnitSystem::measure::density, std::move(density_[gasPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (viscosity_[waterPhaseIdx].size() > 0) {
sol.insert("WAT_VISC", Opm::UnitSystem::measure::viscosity, std::move(viscosity_[waterPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (viscosity_[oilPhaseIdx].size() > 0) {
sol.insert("OIL_VISC", Opm::UnitSystem::measure::viscosity, std::move(viscosity_[oilPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (viscosity_[gasPhaseIdx].size() > 0) {
sol.insert("GAS_VISC", Opm::UnitSystem::measure::viscosity, std::move(viscosity_[gasPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
if (relativePermeability_[waterPhaseIdx].size() > 0) {
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
sol.insert("WATKR", Opm::UnitSystem::measure::identity, std::move(relativePermeability_[waterPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (relativePermeability_[oilPhaseIdx].size() > 0) {
sol.insert("OILKR", Opm::UnitSystem::measure::identity, std::move(relativePermeability_[oilPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (relativePermeability_[gasPhaseIdx].size() > 0) {
sol.insert("GASKR", Opm::UnitSystem::measure::identity, std::move(relativePermeability_[gasPhaseIdx]), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
if (pcSwMdcOw_.size() > 0)
sol.insert ("PCSWM_OW", Opm::UnitSystem::measure::identity, std::move(pcSwMdcOw_), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (krnSwMdcOw_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("KRNSW_OW", Opm::UnitSystem::measure::identity, std::move(krnSwMdcOw_), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (pcSwMdcGo_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("PCSWM_GO", Opm::UnitSystem::measure::identity, std::move(pcSwMdcGo_), Opm::data::TargetType::RESTART_AUXILIARY);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (krnSwMdcGo_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("KRNSW_GO", Opm::UnitSystem::measure::identity, std::move(krnSwMdcGo_), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (soMax_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("SOMAX", Opm::UnitSystem::measure::identity, std::move(soMax_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (sSol_.size() > 0)
change ecl output name from SSOL to SSOLVENT
2018-05-09 07:52:43 -05:00
sol.insert ("SSOLVENT", Opm::UnitSystem::measure::identity, std::move(sSol_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (cPolymer_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("POLYMER", Opm::UnitSystem::measure::identity, std::move(cPolymer_), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Use foam module.
2019-04-26 03:38:37 -05:00
if (cFoam_.size() > 0)
sol.insert ("FOAM", Opm::UnitSystem::measure::identity, std::move(cFoam_), Opm::data::TargetType::RESTART_SOLUTION);
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
if (cSalt_.size() > 0)
sol.insert ("SALT", Opm::UnitSystem::measure::identity, std::move(cSalt_), Opm::data::TargetType::RESTART_SOLUTION);
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (dewPointPressure_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("PDEW", Opm::UnitSystem::measure::pressure, std::move(dewPointPressure_), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (bubblePointPressure_.size() > 0)
Move instead of copy data to solution vector
2018-01-09 04:28:20 -06:00
sol.insert ("PBUB", Opm::UnitSystem::measure::pressure, std::move(bubblePointPressure_), Opm::data::TargetType::RESTART_AUXILIARY);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
if (swMax_.size() > 0)
sol.insert ("SWMAX", Opm::UnitSystem::measure::identity, std::move(swMax_), Opm::data::TargetType::RESTART_SOLUTION);
if (minimumOilPressure_.size() > 0)
sol.insert ("PRESROCC", Opm::UnitSystem::measure::pressure, std::move(minimumOilPressure_), Opm::data::TargetType::RESTART_SOLUTION);
if (overburdenPressure_.size() > 0)
sol.insert ("PRES_OVB", Opm::UnitSystem::measure::pressure, std::move(overburdenPressure_), Opm::data::TargetType::RESTART_SOLUTION);
if (rockCompPorvMultiplier_.size() > 0)
sol.insert ("PORV_RC", Opm::UnitSystem::measure::identity, std::move(rockCompPorvMultiplier_), Opm::data::TargetType::RESTART_SOLUTION);
if (rockCompTransMultiplier_.size() > 0)
sol.insert ("TMULT_RC", Opm::UnitSystem::measure::identity, std::move(rockCompTransMultiplier_), Opm::data::TargetType::RESTART_SOLUTION);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
// Fluid in place
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (outputFipRestart_ && fip_[i].size() > 0) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
sol.insert(fipEnumToString_(i),
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Opm::UnitSystem::measure::volume,
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
fip_[i],
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Opm::data::TargetType::SUMMARY);
}
}
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
// tracers
const auto& tracerModel = simulator_.problem().tracerModel();
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
if (tracerConcentrations_.size() > 0) {
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
for (int tracerIdx = 0; tracerIdx<tracerModel.numTracers(); tracerIdx++){
std::string tmp = tracerModel.tracerName(tracerIdx) + "F";
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
sol.insert(tmp, Opm::UnitSystem::measure::identity, std::move(tracerConcentrations_[tracerIdx]), Opm::data::TargetType::RESTART_SOLUTION);
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
}
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
Clean up PR review
2018-01-16 07:24:21 -06:00
// write Fluid In Place to output log
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
void outputFipLog(std::map<std::string, double>& miscSummaryData, std::map<std::string, std::vector<double>>& regionData, const bool substep)
{
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
const auto& comm = simulator_.gridView().comm();
Cater for empty fipnum vectors in output code. Due to loadbalancing this vector might be empty under some circumstances. Hence the old code would dereference the end iterator which is a bug.
2019-03-22 09:14:08 -05:00
auto maxElement = std::max_element(fipnum_.begin(), fipnum_.end());
size_t ntFip = 0;
if ( maxElement != fipnum_.end() ) {
ntFip = *maxElement;
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
ntFip = comm.max(ntFip);
// sum values over each region
Some more minor clean-up
2018-01-29 05:23:23 -06:00
ScalarBuffer regionFipValues[FipDataType::numFipValues];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
for (int i = 0; i < FipDataType::numFipValues; i++) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
regionFipValues[i] = computeFipForRegions_(fip_[i], fipnum_, ntFip);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
if (isIORank_() && origRegionValues_[i].empty())
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
origRegionValues_[i] = regionFipValues[i];
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
// sum all region values to compute the field total
std::vector<int> fieldNum(ntFip, 1);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
ScalarBuffer fieldFipValues(FipDataType::numFipValues, 0.0);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
bool comunicateSum = false; // the regionValues are already summed over all ranks.
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
const ScalarBuffer& tmp = computeFipForRegions_(regionFipValues[i], fieldNum, 1, comunicateSum);
fieldFipValues[i] = tmp[0];
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// compute the hydrocarbon averaged pressure over the regions.
Minor clean-up after review
2018-01-29 01:43:42 -06:00
ScalarBuffer regPressurePv = computeFipForRegions_(pressureTimesPoreVolume_, fipnum_, ntFip);
ScalarBuffer regPvHydrocarbon = computeFipForRegions_(hydrocarbonPoreVolume_, fipnum_, ntFip);
ScalarBuffer regPressurePvHydrocarbon = computeFipForRegions_(pressureTimesHydrocarbonVolume_, fipnum_, ntFip);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
Minor clean-up after review
2018-01-29 01:43:42 -06:00
ScalarBuffer fieldPressurePv = computeFipForRegions_(regPressurePv, fieldNum, 1, comunicateSum);
ScalarBuffer fieldPvHydrocarbon = computeFipForRegions_(regPvHydrocarbon, fieldNum, 1, comunicateSum);
ScalarBuffer fieldPressurePvHydrocarbon = computeFipForRegions_(regPressurePvHydrocarbon, fieldNum, 1, comunicateSum);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// output on io rank
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
// the original Fip values are stored on the first step
// TODO: Store initial Fip in the init file and restore them
// and use them here.
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::SummaryConfig summaryConfig = simulator_.vanguard().summaryConfig();
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (isIORank_()) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// Field summary output
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
std::string key = "F" + fipEnumToString_(i);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (summaryConfig.hasKeyword(key))
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
miscSummaryData[key] = fieldFipValues[i];
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (summaryConfig.hasKeyword("FOE") && !origTotalValues_.empty())
Some more minor clean-up
2018-01-29 05:23:23 -06:00
miscSummaryData["FOE"] = fieldFipValues[FipDataType::OilInPlace] / origTotalValues_[FipDataType::OilInPlace];
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (summaryConfig.hasKeyword("FPR"))
Some more minor clean-up
2018-01-29 05:23:23 -06:00
miscSummaryData["FPR"] = pressureAverage_(fieldPressurePvHydrocarbon[0], fieldPvHydrocarbon[0], fieldPressurePv[0], fieldFipValues[FipDataType::PoreVolume], true);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (summaryConfig.hasKeyword("FPRP"))
Fix FPRP output
2019-05-27 04:24:55 -05:00
miscSummaryData["FPRP"] = pressureAverage_(fieldPressurePvHydrocarbon[0], fieldPvHydrocarbon[0], fieldPressurePv[0], fieldFipValues[FipDataType::PoreVolume], false);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// Region summary output
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
std::string key = "R" + fipEnumToString_(i);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (summaryConfig.hasKeyword(key))
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
regionData[key] = regionFipValues[i];
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
}
if (summaryConfig.hasKeyword("RPR"))
Some more minor clean-up
2018-01-29 05:23:23 -06:00
regionData["RPR"] = pressureAverage_(regPressurePvHydrocarbon, regPvHydrocarbon, regPressurePv, regionFipValues[FipDataType::PoreVolume], true);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (summaryConfig.hasKeyword("RPRP"))
Some more minor clean-up
2018-01-29 05:23:23 -06:00
regionData["RPRP"] = pressureAverage_(regPressurePvHydrocarbon, regPvHydrocarbon, regPressurePv, regionFipValues[FipDataType::PoreVolume], false);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// Output to log
if (!substep) {
Minor clean-up after review
2018-01-29 01:43:42 -06:00
fipUnitConvert_(fieldFipValues);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (origTotalValues_.empty())
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
origTotalValues_ = fieldFipValues;
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
Some more minor clean-up
2018-01-29 05:23:23 -06:00
Scalar fieldHydroCarbonPoreVolumeAveragedPressure = pressureAverage_(fieldPressurePvHydrocarbon[0], fieldPvHydrocarbon[0], fieldPressurePv[0], fieldFipValues[FipDataType::PoreVolume], true);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
pressureUnitConvert_(fieldHydroCarbonPoreVolumeAveragedPressure);
outputRegionFluidInPlace_(origTotalValues_, fieldFipValues, fieldHydroCarbonPoreVolumeAveragedPressure, 0);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
for (size_t reg = 0; reg < ntFip; ++reg) {
ScalarBuffer tmpO(FipDataType::numFipValues, 0.0);
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
tmpO[i] = origRegionValues_[i][reg];
}
Minor clean-up after review
2018-01-29 01:43:42 -06:00
fipUnitConvert_(tmpO);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
ScalarBuffer tmp(FipDataType::numFipValues, 0.0);
Some more minor clean-up
2018-01-29 05:23:23 -06:00
for (int i = 0; i<FipDataType::numFipValues; i++) {
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
tmp[i] = regionFipValues[i][reg];
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
}
Minor clean-up after review
2018-01-29 01:43:42 -06:00
fipUnitConvert_(tmp);
Some more minor clean-up
2018-01-29 05:23:23 -06:00
Scalar regHydroCarbonPoreVolumeAveragedPressure = pressureAverage_(regPressurePvHydrocarbon[reg], regPvHydrocarbon[reg], regPressurePv[reg], regionFipValues[FipDataType::PoreVolume][reg], true);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
pressureUnitConvert_(regHydroCarbonPoreVolumeAveragedPressure);
outputRegionFluidInPlace_(tmpO, tmp, regHydroCarbonPoreVolumeAveragedPressure, reg + 1);
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
}
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
fix
2019-10-30 11:51:49 -05:00
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
// write production report to output
void outputProdLog(size_t reportStepNum, const bool substep, bool forceDisableProdOutput)
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if (!substep) {
ScalarBuffer tmp_values(WellProdDataType::numWPValues, 0.0);
StringBuffer tmp_names(WellProdDataType::numWPNames, "");
outputProductionReport_(tmp_values, tmp_names, forceDisableProdOutput);
const auto& st = simulator_.vanguard().summaryState();
const auto& schedule = simulator_.vanguard().schedule();
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
for (const auto& gname: schedule.groupNames()) {
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto gName = static_cast<std::string>(gname);
auto get = [&st, &gName](const std::string& vector)
{
const auto key = vector + ':' + gName;
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
return st.has(key) ? st.get(key) : 0.0;
};
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_names[0] = gname;
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_values[2] = get("GOPR"); //WellProdDataType::OilRate
tmp_values[3] = get("GWPR"); //WellProdDataType::WaterRate
tmp_values[4] = get("GGPR"); //WellProdDataType::GasRate
tmp_values[5] = get("GVPR"); //WellProdDataType::FluidResVol
tmp_values[6] = get("GWCT"); //WellProdDataType::WaterCut
tmp_values[7] = get("GGOR"); //WellProdDataType::GasOilRatio
tmp_values[8] = get("GWPR")/get("GGPR"); //WellProdDataType::WaterGasRatio
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
outputProductionReport_(tmp_values, tmp_names, forceDisableProdOutput);
}
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
for (const auto& wname: schedule.wellNames(reportStepNum)) {
// don't bother with wells not on this process
const auto& defunctWellNames = simulator_.vanguard().defunctWellNames();
if (defunctWellNames.find(wname) != defunctWellNames.end()) {
continue;
}
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
const auto& well = schedule.getWell(wname, reportStepNum);
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
// Ignore injector wells
if (well.isInjector()){
continue;
}
tmp_names[0] = wname;//WellProdDataType::WellName
auto wName = static_cast<std::string>(wname);
auto get = [&st, &wName](const std::string& vector)
{
const auto key = vector + ':' + wName;
return st.has(key) ? st.get(key) : 0.0;
};
const auto& controls = well.productionControls(st);
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
using CMode = ::Opm::Well::ProducerCMode;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto fctl = [](const auto wmctl) -> std::string
{
switch (wmctl) {
case CMode::ORAT: return "ORAT";
case CMode::WRAT: return "WRAT";
case CMode::GRAT: return "GRAT";
case CMode::LRAT: return "LRAT";
case CMode::RESV: return "RESV";
case CMode::THP: return "THP";
case CMode::BHP: return "BHP";
case CMode::CRAT: return "CRate";
case CMode::GRUP: return "GRUP";
default:
{
return "none";
}
}
};
tmp_names[1] = fctl(controls.cmode); //WellProdDataType::CTRLMode
tmp_values[0] = well.getHeadI() + 1;//WellProdDataType::WellLocationi
tmp_values[1] = well.getHeadJ() + 1;//WellProdDataType::WellLocationj
tmp_values[2] = get("WOPR"); //WellProdDataType::OilRate
tmp_values[3] = get("WWPR"); //WellProdDataType::WaterRate
tmp_values[4] = get("WGPR"); //WellProdDataType::GasRate
tmp_values[5] = get("WVPR"); //WellProdDataType::FluidResVol
tmp_values[6] = get("WWCT"); //WellProdDataType::WaterCut
tmp_values[7] = get("WGOR"); //WellProdDataType::GasOilRatio
tmp_values[8] = get("WWPR")/get("WGPR"); //WellProdDataType::WaterGasRatio
tmp_values[9] = get("WBHP"); //WellProdDataType::BHP
tmp_values[10] = get("WTHP"); //WellProdDataType::THP
//tmp_values[11] = 0; //WellProdDataType::SteadyStatePI //
outputProductionReport_(tmp_values, tmp_names, forceDisableProdOutput);
}
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
// write injection report to output
void outputInjLog(size_t reportStepNum, const bool substep, bool forceDisableInjOutput)
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if (!substep) {
ScalarBuffer tmp_values(WellInjDataType::numWIValues, 0.0);
StringBuffer tmp_names(WellInjDataType::numWINames, "");
outputInjectionReport_(tmp_values, tmp_names, forceDisableInjOutput);
const auto& st = simulator_.vanguard().summaryState();
const auto& schedule = simulator_.vanguard().schedule();
for (const auto& gname: schedule.groupNames()) {
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto gName = static_cast<std::string>(gname);
auto get = [&st, &gName](const std::string& vector)
{
const auto key = vector + ':' + gName;
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
return st.has(key) ? st.get(key) : 0.0;
};
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_names[0] = gname;
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_values[2] = get("GOIR");//WellInjDataType::OilRate
tmp_values[3] = get("GWIR"); //WellInjDataType::WaterRate
tmp_values[4] = get("GGIR"); //WellInjDataType::GasRate
tmp_values[5] = get("GVIR");//WellInjDataType::FluidResVol
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
outputInjectionReport_(tmp_values, tmp_names, forceDisableInjOutput);
}
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
for (const auto& wname: schedule.wellNames(reportStepNum)) {
// don't bother with wells not on this process
const auto& defunctWellNames = simulator_.vanguard().defunctWellNames();
if (defunctWellNames.find(wname) != defunctWellNames.end()) {
continue;
}
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
const auto& well = schedule.getWell(wname, reportStepNum);
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
// Ignore Producer wells
if (well.isProducer()){
continue;
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_names[0] = wname; //WellInjDataType::WellName
auto wName = static_cast<std::string>(wname);
auto get = [&st, &wName](const std::string& vector)
{
const auto key = vector + ':' + wName;
return st.has(key) ? st.get(key) : 0.0;
};
const auto& controls = well.injectionControls(st);
const auto ctlMode = controls.cmode;
const auto injType = controls.injector_type;
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
using CMode = ::Opm::Well::InjectorCMode;
InjectorType enum is moved out from Well class
2020-03-03 02:08:02 -06:00
using WType = ::Opm::InjectorType;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto ftype = [](const auto wtype) -> std::string
{
switch (wtype) {
case WType::OIL: return "Oil";
case WType::WATER: return "Wat";
case WType::GAS: return "Gas";
case WType::MULTI: return "Multi";
default:
{
return "";
}
}
};
auto fctl = [](const auto wmctl) -> std::string
{
switch (wmctl) {
case CMode::RATE: return "RATE";
case CMode::RESV: return "RESV";
case CMode::THP: return "THP";
case CMode::BHP: return "BHP";
case CMode::GRUP: return "GRUP";
default:
{
return "";
}
}
};
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
const auto flowtype = ftype(injType);
const auto flowctl = fctl(ctlMode);
if(flowtype == "Oil") //WellInjDataType::CTRLModeOil
{
if (flowctl == "RATE"){ tmp_names[1] = "ORAT"; }
else { tmp_names[1] = flowctl; }
}
else if (flowtype == "Wat") //WellInjDataType::CTRLModeWat
{
if (flowctl == "RATE"){ tmp_names[3] = "WRAT"; }
else { tmp_names[2] = flowctl; }
}
else if (flowtype == "Gas") //WellInjDataType::CTRLModeGas
{
if (flowctl == "RATE"){ tmp_names[3] = "GRAT"; }
else { tmp_names[3] = flowctl; }
}
tmp_values[0] = well.getHeadI() + 1; //WellInjDataType::wellLocationi
tmp_values[1] = well.getHeadJ() + 1; //WellInjDataType::wellLocationj
tmp_values[2] = get("WOIR"); //WellInjDataType::OilRate
tmp_values[3] = get("WWIR"); //WellInjDataType::WaterRate
tmp_values[4] = get("WGIR"); //WellInjDataType::GasRate
tmp_values[5] = get("WVIR");//WellInjDataType::FluidResVol
tmp_values[6] = get("WBHP"); //WellInjDataType::BHP
tmp_values[7] = get("WTHP"); //WellInjDataType::THP
//tmp_values[8] = 0; //WellInjDataType::SteadyStateII
outputInjectionReport_(tmp_values, tmp_names, forceDisableInjOutput);
}
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
// write cumulative production and injection reports to output
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
void outputCumLog(size_t reportStepNum, const bool substep, bool forceDisableCumOutput)
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if (!substep) {
ScalarBuffer tmp_values(WellCumDataType::numWCValues, 0.0);
StringBuffer tmp_names(WellCumDataType::numWCNames, "");
outputCumulativeReport_(tmp_values, tmp_names, forceDisableCumOutput);
const auto& st = simulator_.vanguard().summaryState();
const auto& schedule = simulator_.vanguard().schedule();
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
for (const auto& gname: schedule.groupNames()) {
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto gName = static_cast<std::string>(gname);
auto get = [&st, &gName](const std::string& vector)
{
const auto key = vector + ':' + gName;
return st.has(key) ? st.get(key) : 0.0;
};
tmp_names[0] = gname;
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_values[2] = get("GOPT"); //WellCumDataType::OilProd
tmp_values[3] = get("GWPT"); //WellCumDataType::WaterProd
tmp_values[4] = get("GGPT"); //WellCumDataType::GasProd
tmp_values[5] = get("GVPT");//WellCumDataType::FluidResVolProd
Update ecloutputblackoilmodule.hh
2019-11-06 06:25:04 -06:00
tmp_values[6] = get("GOIT"); //WellCumDataType::OilInj
tmp_values[7] = get("GWIT"); //WellCumDataType::WaterInj
tmp_values[8] = get("GGIT"); //WellCumDataType::GasInj
tmp_values[9] = get("GVIT");//WellCumDataType::FluidResVolInj
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
outputCumulativeReport_(tmp_values, tmp_names, forceDisableCumOutput);
}
added groups to reports
2019-10-30 08:09:51 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
for (const auto& wname : schedule.wellNames(reportStepNum)) {
// don't bother with wells not on this process
const auto& defunctWellNames = simulator_.vanguard().defunctWellNames();
if (defunctWellNames.find(wname) != defunctWellNames.end()) {
continue;
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
const auto& well = schedule.getWell(wname, reportStepNum);
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
tmp_names[0] = wname; //WellCumDataType::WellName
auto wName = static_cast<std::string>(wname);
auto get = [&st, &wName](const std::string& vector)
{
const auto key = vector + ':' + wName;
return st.has(key) ? st.get(key) : 0.0;
};
if (well.isInjector()) {
const auto& controls = well.injectionControls(st);
const auto ctlMode = controls.cmode;
const auto injType = controls.injector_type;
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
using CMode = ::Opm::Well::InjectorCMode;
InjectorType enum is moved out from Well class
2020-03-03 02:08:02 -06:00
using WType = ::Opm::InjectorType;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto ftype = [](const auto wtype) -> std::string
{
switch (wtype) {
case WType::OIL: return "Oil";
case WType::WATER: return "Wat";
case WType::GAS: return "Gas";
case WType::MULTI: return "Multi";
default:
{
return "";
}
}
};
auto fctl = [](const auto wmctl) -> std::string
{
switch (wmctl) {
case CMode::RATE: return "RATE";
case CMode::RESV: return "RESV";
case CMode::THP: return "THP";
case CMode::BHP: return "BHP";
case CMode::GRUP: return "GRUP";
default:
{
return "";
}
}
};
tmp_names[1] = "INJ"; //WellCumDataType::WellType
const auto flowctl = fctl(ctlMode);
if (flowctl == "RATE") //WellCumDataType::WellCTRL
{
const auto flowtype = ftype(injType);
if(flowtype == "Oil"){ tmp_names[2] = "ORAT"; }
else if(flowtype == "Wat"){ tmp_names[2] = "WRAT"; }
else if(flowtype == "Gas"){ tmp_names[2] = "GRAT"; }
}
else
{
tmp_names[2] = flowctl;
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
}
else if (well.isProducer()) {
const auto& controls = well.productionControls(st);
Rename Well2 -> Well and Group2 -> Group
2019-11-13 16:16:11 -06:00
using CMode = ::Opm::Well::ProducerCMode;
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
auto fctl = [](const auto wmctl) -> std::string
{
switch (wmctl) {
case CMode::ORAT: return "ORAT";
case CMode::WRAT: return "WRAT";
case CMode::GRAT: return "GRAT";
case CMode::LRAT: return "LRAT";
case CMode::RESV: return "RESV";
case CMode::THP: return "THP";
case CMode::BHP: return "BHP";
case CMode::CRAT: return "CRAT";
case CMode::GRUP: return "GRUP";
default:
{
return "none";
}
}
};
tmp_names[1] = "PROD"; //WellProdDataType::CTRLMode
tmp_names[2] = fctl(controls.cmode); //WellProdDataType::CTRLMode
}
tmp_values[0] = well.getHeadI() + 1; //WellCumDataType::wellLocationi
tmp_values[1] = well.getHeadJ() + 1; //WellCumDataType::wellLocationj
tmp_values[2] = get("WOPT"); //WellCumDataType::OilProd
tmp_values[3] = get("WWPT"); //WellCumDataType::WaterProd
tmp_values[4] = get("WGPT"); //WellCumDataType::GasProd
tmp_values[5] = get("WVPT");//WellCumDataType::FluidResVolProd
Update ecloutputblackoilmodule.hh
2019-11-06 06:25:04 -06:00
tmp_values[6] = get("WOIT"); //WellCumDataType::OilInj
tmp_values[7] = get("WWIT"); //WellCumDataType::WaterInj
tmp_values[8] = get("WGIT"); //WellCumDataType::GasInj
tmp_values[9] = get("WVIT");//WellCumDataType::FluidResVolInj
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
outputCumulativeReport_(tmp_values, tmp_names, forceDisableCumOutput);
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
}
fix
2019-10-30 11:51:49 -05:00
ebos: replace tabs with spaces and remove trailing white space for some reason this keeps creping back in...
2019-06-03 04:11:41 -05:00
void setRestart(const Opm::data::Solution& sol, unsigned elemIdx, unsigned globalDofIndex)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
{
Scalar so = 1.0;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (saturation_[waterPhaseIdx].size() > 0 && sol.has("SWAT")) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
saturation_[waterPhaseIdx][elemIdx] = sol.data("SWAT")[globalDofIndex];
so -= sol.data("SWAT")[globalDofIndex];
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (saturation_[gasPhaseIdx].size() > 0 && sol.has("SGAS")) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
saturation_[gasPhaseIdx][elemIdx] = sol.data("SGAS")[globalDofIndex];
so -= sol.data("SGAS")[globalDofIndex];
}
Fix solvent restart
2019-09-25 05:28:58 -05:00
if (sSol_.size() > 0) {
// keep the SSOL option for backward compatibility
// should be removed after 10.2018 release
if (sol.has("SSOL"))
sSol_[elemIdx] = sol.data("SSOL")[globalDofIndex];
else if (sol.has("SSOLVENT"))
sSol_[elemIdx] = sol.data("SSOLVENT")[globalDofIndex];
so -= sSol_[elemIdx];
}
Make the restart code more robust. Don't assume all needed data is given in the restart file. i.e. make it possible to run a solvent restart run from a standard blackoil run.
2018-02-19 01:49:02 -06:00
assert(saturation_[oilPhaseIdx].size() > 0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
saturation_[oilPhaseIdx][elemIdx] = so;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (oilPressure_.size() > 0 && sol.has("PRESSURE"))
oilPressure_[elemIdx] = sol.data("PRESSURE")[globalDofIndex];
use `enableEnergy` consistently thanks to [at]totto82 for noticing.
2018-12-17 04:37:07 -06:00
if (enableEnergy && sol.has("TEMP"))
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
temperature_[elemIdx] = sol.data("TEMP")[globalDofIndex];
if (rs_.size() > 0 && sol.has("RS"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
rs_[elemIdx] = sol.data("RS")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (rv_.size() > 0 && sol.has("RV"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
rv_[elemIdx] = sol.data("RV")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (cPolymer_.size() > 0 && sol.has("POLYMER"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
cPolymer_[elemIdx] = sol.data("POLYMER")[globalDofIndex];
Use foam module.
2019-04-26 03:38:37 -05:00
if (cFoam_.size() > 0 && sol.has("FOAM"))
cFoam_[elemIdx] = sol.data("FOAM")[globalDofIndex];
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
if (cSalt_.size() > 0 && sol.has("SALT"))
cSalt_[elemIdx] = sol.data("SALT")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (soMax_.size() > 0 && sol.has("SOMAX"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
soMax_[elemIdx] = sol.data("SOMAX")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (pcSwMdcOw_.size() > 0 &&sol.has("PCSWM_OW"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
pcSwMdcOw_[elemIdx] = sol.data("PCSWM_OW")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (krnSwMdcOw_.size() > 0 && sol.has("KRNSW_OW"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
krnSwMdcOw_[elemIdx] = sol.data("KRNSW_OW")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (pcSwMdcGo_.size() > 0 && sol.has("PCSWM_GO"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
pcSwMdcGo_[elemIdx] = sol.data("PCSWM_GO")[globalDofIndex];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (krnSwMdcGo_.size() > 0 && sol.has("KRNSW_GO"))
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
krnSwMdcGo_[elemIdx] = sol.data("KRNSW_GO")[globalDofIndex];
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
if (ppcw_.size() > 0 && sol.has("PPCW"))
ppcw_[elemIdx] = sol.data("PPCW")[globalDofIndex];
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
template <class FluidState>
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
void assignToFluidState(FluidState& fs, unsigned elemIdx) const
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
{
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++ phaseIdx) {
if (saturation_[phaseIdx].size() == 0)
continue;
ECL output: use the same units as the input deck for ECL output this makes it possible to compare ECL restart and summary files produced by ebos directly with the ones stemming from Eclips. (But be aware that VTK output files from ebos are still all-SI!)
2015-01-25 10:27:08 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
fs.setSaturation(phaseIdx, saturation_[phaseIdx][elemIdx]);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (oilPressure_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
// this assumes that capillary pressures only depend on the phase saturations
// and possibly on temperature. (this is always the case for ECL problems.)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
Dune::FieldVector< Scalar, numPhases > pc(0);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
const MaterialLawParams& matParams = simulator_.problem().materialLawParams(elemIdx);
MaterialLaw::capillaryPressures(pc, matParams, fs);
Opm::Valgrind::CheckDefined(oilPressure_[elemIdx]);
Opm::Valgrind::CheckDefined(pc);
assert(FluidSystem::phaseIsActive(oilPhaseIdx));
for (unsigned phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx) {
if (!FluidSystem::phaseIsActive(phaseIdx))
continue;
fs.setPressure(phaseIdx, oilPressure_[elemIdx] + (pc[phaseIdx] - pc[oilPhaseIdx]));
}
ECL output: use the same units as the input deck for ECL output this makes it possible to compare ECL restart and summary files produced by ebos directly with the ones stemming from Eclips. (But be aware that VTK output files from ebos are still all-SI!)
2015-01-25 10:27:08 -06:00
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
use `enableEnergy` consistently thanks to [at]totto82 for noticing.
2018-12-17 04:37:07 -06:00
if (enableEnergy)
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
fs.setTemperature(temperature_[elemIdx]);
if (rs_.size() > 0)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
fs.setRs(rs_[elemIdx]);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
if (rv_.size() > 0)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
fs.setRv(rv_[elemIdx]);
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
}
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
void initHysteresisParams(Simulator& simulator, unsigned elemIdx) const
{
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (soMax_.size() > 0)
Fix bug in setting the max oil saturation after restart
2018-02-19 01:46:27 -06:00
simulator.problem().setMaxOilSaturation(elemIdx, soMax_[elemIdx]);
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
if (simulator.problem().materialLawManager()->enableHysteresis()) {
auto matLawManager = simulator.problem().materialLawManager();
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
if (pcSwMdcOw_.size() > 0 && krnSwMdcOw_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
matLawManager->setOilWaterHysteresisParams(
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
pcSwMdcOw_[elemIdx],
krnSwMdcOw_[elemIdx],
elemIdx);
}
if (pcSwMdcGo_.size() > 0 && krnSwMdcGo_.size() > 0) {
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
matLawManager->setGasOilHysteresisParams(
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
pcSwMdcGo_[elemIdx],
krnSwMdcGo_[elemIdx],
elemIdx);
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
Simulate with only active cells (#2213) Use FieldProps implementation for 3D properties
2020-01-13 08:46:50 -06:00
if (simulator_.vanguard().eclState().fieldProps().has_double("SWATINIT")) {
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
auto oilWaterScaledEpsInfoDrainage = simulator.problem().materialLawManager()->oilWaterScaledEpsInfoDrainagePointerReferenceHack(elemIdx);
oilWaterScaledEpsInfoDrainage->maxPcow = ppcw_[elemIdx];
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
}
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
Scalar getSolventSaturation(unsigned elemIdx) const
{
Improve out-of-range check for extra components.
2019-08-07 05:03:43 -05:00
if (sSol_.size() > elemIdx)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
return sSol_[elemIdx];
return 0;
}
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
Scalar getPolymerConcentration(unsigned elemIdx) const
{
Improve out-of-range check for extra components.
2019-08-07 05:03:43 -05:00
if (cPolymer_.size() > elemIdx)
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
return cPolymer_[elemIdx];
return 0;
}
Use foam module.
2019-04-26 03:38:37 -05:00
Scalar getFoamConcentration(unsigned elemIdx) const
{
Improve out-of-range check for extra components.
2019-08-07 05:03:43 -05:00
if (cFoam_.size() > elemIdx)
Use foam module.
2019-04-26 03:38:37 -05:00
return cFoam_[elemIdx];
return 0;
}
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
Scalar getSaltConcentration(unsigned elemIdx) const
{
if (cSalt_.size() > elemIdx)
return cSalt_[elemIdx];
return 0;
}
fix some minor coding style issues in Ewoms::EclOutputBlackOilModule
2019-02-01 10:33:30 -06:00
const std::map<std::pair<std::string, int>, double>& getBlockData()
{ return blockData_; }
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
private:
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
Only compute auxiliary values if asked for Ask the restartConfig if restart files are written. Solution variables are always written since they may be needed by the summary files Remove output of saturated data to align with legacy code.
2018-01-09 07:01:30 -06:00
bool isIORank_() const
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
{
const auto& comm = simulator_.gridView().comm();
return comm.rank() == 0;
}
Minor clean-up after review
2018-01-29 01:43:42 -06:00
void updateFluidInPlace_(const ElementContext& elemCtx, unsigned dofIdx)
{
const auto& intQuants = elemCtx.intensiveQuantities(dofIdx, /*timeIdx=*/0);
const auto& fs = intQuants.fluidState();
unsigned globalDofIdx = elemCtx.globalSpaceIndex(dofIdx, /*timeIdx=*/0);
// Fluid in Place calculations
// calculate the pore volume of the current cell. Note that the porosity
// returned by the intensive quantities is defined as the ratio of pore
// space to total cell volume and includes all pressure dependent (->
// rock compressibility) and static modifiers (MULTPV, MULTREGP, NTG,
// PORV, MINPV and friends). Also note that because of this, the porosity
// returned by the intensive quantities can be outside of the physical
// range [0, 1] in pathetic cases.
const double pv =
elemCtx.simulator().model().dofTotalVolume(globalDofIdx)
* intQuants.porosity().value();
if (pressureTimesHydrocarbonVolume_.size() > 0 && pressureTimesPoreVolume_.size() > 0) {
assert(hydrocarbonPoreVolume_.size() == pressureTimesHydrocarbonVolume_.size());
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
assert(fip_[FipDataType::PoreVolume].size() == pressureTimesPoreVolume_.size());
Minor clean-up after review
2018-01-29 01:43:42 -06:00
Some more minor clean-up
2018-01-29 05:23:23 -06:00
fip_[FipDataType::PoreVolume][globalDofIdx] = pv;
Minor clean-up after review
2018-01-29 01:43:42 -06:00
Scalar hydrocarbon = 0.0;
if (FluidSystem::phaseIsActive(oilPhaseIdx))
hydrocarbon += Opm::getValue(fs.saturation(oilPhaseIdx));
if (FluidSystem::phaseIsActive(gasPhaseIdx))
hydrocarbon += Opm::getValue(fs.saturation(gasPhaseIdx));
hydrocarbonPoreVolume_[globalDofIdx] = pv * hydrocarbon;
if (FluidSystem::phaseIsActive(oilPhaseIdx)) {
pressureTimesPoreVolume_[globalDofIdx] = Opm::getValue(fs.pressure(oilPhaseIdx)) * pv;
pressureTimesHydrocarbonVolume_[globalDofIdx] = pressureTimesPoreVolume_[globalDofIdx] * hydrocarbon;
}
}
if (computeFip_) {
Scalar fip[FluidSystem::numPhases];
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
for (unsigned phaseIdx = 0; phaseIdx < FluidSystem::numPhases; ++phaseIdx) {
fip[phaseIdx] = 0.0;
if (!FluidSystem::phaseIsActive(phaseIdx))
Minor clean-up after review
2018-01-29 01:43:42 -06:00
continue;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
const double b = Opm::getValue(fs.invB(phaseIdx));
const double s = Opm::getValue(fs.saturation(phaseIdx));
fip[phaseIdx] = b * s * pv;
Minor clean-up after review
2018-01-29 01:43:42 -06:00
}
Some more minor clean-up
2018-01-29 05:23:23 -06:00
if (FluidSystem::phaseIsActive(oilPhaseIdx) && fip_[FipDataType::OilInPlace].size() > 0)
fip_[FipDataType::OilInPlace][globalDofIdx] = fip[oilPhaseIdx];
if (FluidSystem::phaseIsActive(gasPhaseIdx) && fip_[FipDataType::GasInPlace].size() > 0)
fip_[FipDataType::GasInPlace][globalDofIdx] = fip[gasPhaseIdx];
if (FluidSystem::phaseIsActive(waterPhaseIdx) && fip_[FipDataType::WaterInPlace].size() > 0)
fip_[FipDataType::WaterInPlace][globalDofIdx] = fip[waterPhaseIdx];
Minor clean-up after review
2018-01-29 01:43:42 -06:00
Some more minor clean-up
2018-01-29 05:23:23 -06:00
// Store the pure oil and gas Fip
if (FluidSystem::phaseIsActive(oilPhaseIdx) && fip_[FipDataType::OilInPlaceInLiquidPhase].size() > 0)
fip_[FipDataType::OilInPlaceInLiquidPhase][globalDofIdx] = fip[oilPhaseIdx];
Minor clean-up after review
2018-01-29 01:43:42 -06:00
Some more minor clean-up
2018-01-29 05:23:23 -06:00
if (FluidSystem::phaseIsActive(gasPhaseIdx) && fip_[FipDataType::GasInPlaceInGasPhase].size() > 0)
fip_[FipDataType::GasInPlaceInGasPhase][globalDofIdx] = fip[gasPhaseIdx];
Minor clean-up after review
2018-01-29 01:43:42 -06:00
if (FluidSystem::phaseIsActive(oilPhaseIdx) && FluidSystem::phaseIsActive(gasPhaseIdx)) {
// Gas dissolved in oil and vaporized oil
Scalar gasInPlaceLiquid = Opm::getValue(fs.Rs()) * fip[oilPhaseIdx];
Scalar oilInPlaceGas = Opm::getValue(fs.Rv()) * fip[gasPhaseIdx];
fix typo in fip output
2019-05-22 08:39:57 -05:00
if (fip_[FipDataType::GasInPlaceInLiquidPhase].size() > 0)
Some more minor clean-up
2018-01-29 05:23:23 -06:00
fip_[FipDataType::GasInPlaceInLiquidPhase][globalDofIdx] = gasInPlaceLiquid;
if (fip_[FipDataType::OilInPlaceInGasPhase].size() > 0)
fip_[FipDataType::OilInPlaceInGasPhase][globalDofIdx] = oilInPlaceGas;
// Add dissolved gas and vaporized oil to total Fip
if (fip_[FipDataType::OilInPlace].size() > 0)
fip_[FipDataType::OilInPlace][globalDofIdx] += oilInPlaceGas;
if (fip_[FipDataType::GasInPlace].size() > 0)
fip_[FipDataType::GasInPlace][globalDofIdx] += gasInPlaceLiquid;
Minor clean-up after review
2018-01-29 01:43:42 -06:00
}
}
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
void createLocalFipnum_()
{
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const auto& gridView = simulator_.vanguard().gridView();
use compressed field property for fipnum
2020-02-25 07:27:52 -06:00
if (simulator_.vanguard().eclState().fieldProps().has_int("FIPNUM"))
fipnum_ = simulator_.vanguard().eclState().fieldProps().get_int("FIPNUM");
else
fipnum_.resize(gridView.size(/*codim=*/0), 1);
ElementContext elemCtx(simulator_);
ElementIterator elemIt = simulator_.gridView().template begin</*codim=*/0>();
const ElementIterator& elemEndIt = simulator_.gridView().template end</*codim=*/0>();
size_t elemIdx = 0;
for (; elemIt != elemEndIt; ++elemIt, ++elemIdx) {
const Element& elem = *elemIt;
if (elem.partitionType() != Dune::InteriorEntity)
fipnum_[elemIdx] = 0;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
}
// Sum Fip values over regions.
Minor clean-up after review
2018-01-29 01:43:42 -06:00
ScalarBuffer computeFipForRegions_(const ScalarBuffer& fip, std::vector<int>& regionId, size_t maxNumberOfRegions, bool commSum = true)
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
{
ScalarBuffer totals(maxNumberOfRegions, 0.0);
Minor clean-up based on review
2018-02-01 06:16:37 -06:00
if (fip.empty())
return totals;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
assert(regionId.size() == fip.size());
for (size_t j = 0; j < regionId.size(); ++j) {
const int regionIdx = regionId[j] - 1;
// the cell is not attributed to any region. ignore it!
if (regionIdx < 0)
continue;
assert(regionIdx < static_cast<int>(maxNumberOfRegions));
totals[regionIdx] += fip[j];
}
if (commSum) {
const auto& comm = simulator_.gridView().comm();
for (size_t i = 0; i < maxNumberOfRegions; ++i)
totals[i] = comm.sum(totals[i]);
}
return totals;
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
ScalarBuffer pressureAverage_(const ScalarBuffer& pressurePvHydrocarbon, const ScalarBuffer& pvHydrocarbon, const ScalarBuffer& pressurePv, const ScalarBuffer& pv, bool hydrocarbon)
{
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
size_t size = pressurePvHydrocarbon.size();
assert(pvHydrocarbon.size() == size);
assert(pressurePv.size() == size);
assert(pv.size() == size);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
ScalarBuffer fraction(size, 0.0);
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
for (size_t i = 0; i < size; ++i) {
fraction[i] = pressureAverage_(pressurePvHydrocarbon[i], pvHydrocarbon[i], pressurePv[i], pv[i], hydrocarbon);
}
return fraction;
}
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
Scalar pressureAverage_(const Scalar& pressurePvHydrocarbon, const Scalar& pvHydrocarbon, const Scalar& pressurePv, const Scalar& pv, bool hydrocarbon)
{
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
if (pvHydrocarbon > 1e-10 && hydrocarbon)
return pressurePvHydrocarbon / pvHydrocarbon;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
return pressurePv / pv;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
Minor clean-up after review
2018-01-29 01:43:42 -06:00
void fipUnitConvert_(ScalarBuffer& fip)
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
{
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
Some more minor clean-up
2018-01-29 05:23:23 -06:00
fip[FipDataType::WaterInPlace] = Opm::unit::convert::to(fip[FipDataType::WaterInPlace], Opm::unit::stb);
fip[FipDataType::OilInPlace] = Opm::unit::convert::to(fip[FipDataType::OilInPlace], Opm::unit::stb);
fip[FipDataType::OilInPlaceInLiquidPhase] = Opm::unit::convert::to(fip[FipDataType::OilInPlaceInLiquidPhase], Opm::unit::stb);
fip[FipDataType::OilInPlaceInGasPhase] = Opm::unit::convert::to(fip[FipDataType::OilInPlaceInGasPhase], Opm::unit::stb);
fip[FipDataType::GasInPlace] = Opm::unit::convert::to(fip[FipDataType::GasInPlace], 1000*Opm::unit::cubic(Opm::unit::feet));
fip[FipDataType::GasInPlaceInLiquidPhase] = Opm::unit::convert::to(fip[FipDataType::GasInPlaceInLiquidPhase], 1000*Opm::unit::cubic(Opm::unit::feet));
fip[FipDataType::GasInPlaceInGasPhase] = Opm::unit::convert::to(fip[FipDataType::GasInPlaceInGasPhase], 1000*Opm::unit::cubic(Opm::unit::feet));
fip[FipDataType::PoreVolume] = Opm::unit::convert::to(fip[FipDataType::PoreVolume], Opm::unit::stb);
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
convert FIP values to LAB scale
2019-01-21 07:55:40 -06:00
else if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_LAB) {
Scalar scc = Opm::unit::cubic(Opm::prefix::centi * Opm::unit::meter); //standard cubic cm.
fip[FipDataType::WaterInPlace] = Opm::unit::convert::to(fip[FipDataType::WaterInPlace], scc);
fip[FipDataType::OilInPlace] = Opm::unit::convert::to(fip[FipDataType::OilInPlace], scc);
fip[FipDataType::OilInPlaceInLiquidPhase] = Opm::unit::convert::to(fip[FipDataType::OilInPlaceInLiquidPhase], scc);
fip[FipDataType::OilInPlaceInGasPhase] = Opm::unit::convert::to(fip[FipDataType::OilInPlaceInGasPhase], scc);
fip[FipDataType::GasInPlace] = Opm::unit::convert::to(fip[FipDataType::GasInPlace], scc);
fip[FipDataType::GasInPlaceInLiquidPhase] = Opm::unit::convert::to(fip[FipDataType::GasInPlaceInLiquidPhase], scc);
fip[FipDataType::GasInPlaceInGasPhase] = Opm::unit::convert::to(fip[FipDataType::GasInPlaceInGasPhase], scc);
fip[FipDataType::PoreVolume] = Opm::unit::convert::to(fip[FipDataType::PoreVolume], scc);
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
else if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
// nothing to do
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
else {
adapt to the move of code from opm-common to opm-material
2018-02-01 07:40:01 -06:00
throw std::runtime_error("Unsupported unit type for fluid in place output.");
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
}
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
void pressureUnitConvert_(Scalar& pav)
{
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
pav = Opm::unit::convert::to(pav, Opm::unit::psia);
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
}
else if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
pav = Opm::unit::convert::to(pav, Opm::unit::barsa);
}
convert FIP values to LAB scale
2019-01-21 07:55:40 -06:00
else if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_LAB) {
pav = Opm::unit::convert::to(pav, Opm::unit::atm);
}
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
else {
adapt to the move of code from opm-common to opm-material
2018-02-01 07:40:01 -06:00
throw std::runtime_error("Unsupported unit type for fluid in place output.");
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
}
}
void outputRegionFluidInPlace_(const ScalarBuffer& oip, const ScalarBuffer& cip, const Scalar& pav, const int reg)
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
{
ebos: add possibility to disable the FIP output even if requested in the deck this significantly reduces the noise level of debug files if the FIP values are not of interest.
2018-07-26 05:31:32 -05:00
if (forceDisableFipOutput_)
return;
Avoid output of trivial fipnum regions
2018-11-06 02:02:30 -06:00
// don't output FIPNUM report if the region has no porv.
if (cip[FipDataType::PoreVolume] == 0)
return;
rename the "grid manager" to "vanguard" IMO the term "vanguard" expresses better what these classes are supposed to do: level the ground for the cavalry. Normally this simply means to create and distribute a grid object, but it can become quite a bit more complicated, as exemplified by the vanguard classes of ebos..
2018-02-01 09:26:58 -06:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
std::ostringstream ss;
if (!reg) {
ss << " ===================================================\n"
<< " : Field Totals :\n";
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
}
else {
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
ss << " ===================================================\n"
<< " : FIPNUM report region "
<< std::setw(2) << reg << " :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
ss << " : PAV =" << std::setw(14) << pav << " BARSA :\n"
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
<< std::fixed << std::setprecision(0)
Some more minor clean-up
2018-01-29 05:23:23 -06:00
<< " : PORV =" << std::setw(14) << cip[FipDataType::PoreVolume] << " RM3 :\n";
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
if (!reg) {
ss << " : Pressure is weighted by hydrocarbon pore volume :\n"
<< " : Porv volumes are taken at reference conditions :\n";
}
ss << " :--------------- Oil SM3 ---------------:-- Wat SM3 --:--------------- Gas SM3 ---------------:\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
Pass region and field FIP to opm-output
2018-01-17 08:24:24 -06:00
ss << " : PAV =" << std::setw(14) << pav << " PSIA :\n"
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
<< std::fixed << std::setprecision(0)
Some more minor clean-up
2018-01-29 05:23:23 -06:00
<< " : PORV =" << std::setw(14) << cip[FipDataType::PoreVolume] << " RB :\n";
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
if (!reg) {
ss << " : Pressure is weighted by hydrocarbon pore volume :\n"
<< " : Pore volumes are taken at reference conditions :\n";
}
ss << " :--------------- Oil STB ---------------:-- Wat STB --:--------------- Gas MSCF ---------------:\n";
}
ss << " : Liquid Vapour Total : Total : Free Dissolved Total :" << "\n"
<< ":------------------------:------------------------------------------:----------------:------------------------------------------:" << "\n"
Some more minor clean-up
2018-01-29 05:23:23 -06:00
<< ":Currently in place :" << std::setw(14) << cip[FipDataType::OilInPlaceInLiquidPhase] << std::setw(14) << cip[FipDataType::OilInPlaceInGasPhase] << std::setw(14) << cip[FipDataType::OilInPlace] << ":"
<< std::setw(13) << cip[FipDataType::WaterInPlace] << " :" << std::setw(14) << (cip[FipDataType::GasInPlaceInGasPhase]) << std::setw(14) << cip[FipDataType::GasInPlaceInLiquidPhase] << std::setw(14) << cip[FipDataType::GasInPlace] << ":\n"
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
<< ":------------------------:------------------------------------------:----------------:------------------------------------------:\n"
Some more minor clean-up
2018-01-29 05:23:23 -06:00
<< ":Originally in place :" << std::setw(14) << oip[FipDataType::OilInPlaceInLiquidPhase] << std::setw(14) << oip[FipDataType::OilInPlaceInGasPhase] << std::setw(14) << oip[FipDataType::OilInPlace] << ":"
<< std::setw(13) << oip[FipDataType::WaterInPlace] << " :" << std::setw(14) << oip[FipDataType::GasInPlaceInGasPhase] << std::setw(14) << oip[FipDataType::GasInPlaceInLiquidPhase] << std::setw(14) << oip[FipDataType::GasInPlace] << ":\n"
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
<< ":========================:==========================================:================:==========================================:\n";
Opm::OpmLog::note(ss.str());
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
void outputProductionReport_(const ScalarBuffer& wellProd, const StringBuffer& wellProdNames, const bool forceDisableProdOutput)
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if(forceDisableProdOutput)
return;
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
std::ostringstream ss;
if (wellProdNames[WellProdDataType::WellName].empty()) {
ss << "======================================================= PRODUCTION REPORT =======================================================\n"//=================== \n"
<< ": WELL : LOCATION :CTRL: OIL : WATER : GAS : FLUID : WATER : GAS/OIL : WAT/GAS : BHP OR : THP OR :\n"// STEADY-ST PI :\n"
<< ": NAME : (I,J,K) :MODE: RATE : RATE : RATE : RES.VOL. : CUT : RATIO : RATIO : CON.PR.: BLK.PR.:\n";// OR POTN OF PREF. PH:\n";
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
ss << ": : : : SCM/DAY : SCM/DAY : SCM/DAY : RCM/DAY : SCM/SCM : SCM/SCM : SCM/SCM : BARSA : BARSA :\n";// :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
ss << ": : : : STB/DAY : STB/DAY : MSCF/DAY : RB/DAY : : MSCF/STB : STB/MSCF : PSIA : PSIA :\n";// :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_LAB) {
ss << ": : : : SCC/HR : SCC/HR : SCC/HR : RCC : SCC/SCC : SCC/SCC : SCC/SCC : ATMA : ATMA :\n";// :\n";
}
ss << "=================================================================================================================================\n";//=================== \n";
}
else {
added groups to reports
2019-10-30 08:09:51 -05:00
if (wellProd[WellProdDataType::WellLocationi] < 1) {
ss << std::right << std::fixed << ":" << std::setw (8) << wellProdNames[WellProdDataType::WellName] << ":" << std::setprecision(0) << std::setw(11) << "" << ":" << std::setw(4) << wellProdNames[WellProdDataType::CTRLMode] << ":" << std::setprecision(1) << std::setw(11) << wellProd[WellProdDataType::OilRate] << ":" << std::setw(11) << wellProd[WellProdDataType::WaterRate] << ":" << std::setw(11)<< wellProd[WellProdDataType::GasRate] << ":" << std::setw(11) << wellProd[WellProdDataType::FluidResVol] << std::setprecision(3) << ":" << std::setw(11) << wellProd[WellProdDataType::WaterCut] << std::setprecision(2) << ":" << std::setw(10) << wellProd[WellProdDataType::GasOilRatio] << std::setprecision(4) << ":" << std::setw(12) << wellProd[WellProdDataType::WatGasRatio] << std::setprecision(1) << ":" << std::setw(8) << "" << ":" << std::setw(8) << "" << ": \n";//wellProd[WellProdDataType::SteadyStatePI] << std::setw(10) << "\n"
}
else {
ss << std::right << std::fixed << ":" << std::setw (8) << wellProdNames[WellProdDataType::WellName] << ":" << std::setprecision(0) << std::setw(5) << wellProd[WellProdDataType::WellLocationi] << "," << std::setw(5) << wellProd[WellProdDataType::WellLocationj] << ":" << std::setw(4) << wellProdNames[WellProdDataType::CTRLMode] << ":" << std::setprecision(1) << std::setw(11) << wellProd[WellProdDataType::OilRate] << ":" << std::setw(11) << wellProd[WellProdDataType::WaterRate] << ":" << std::setw(11)<< wellProd[WellProdDataType::GasRate] << ":" << std::setw(11) << wellProd[WellProdDataType::FluidResVol] << std::setprecision(3) << ":" << std::setw(11) << wellProd[WellProdDataType::WaterCut] << std::setprecision(2) << ":" << std::setw(10) << wellProd[WellProdDataType::GasOilRatio] << std::setprecision(4) << ":" << std::setw(12) << wellProd[WellProdDataType::WatGasRatio] << std::setprecision(1) << ":" << std::setw(8) << wellProd[WellProdDataType::BHP] << ":" << std::setw(8) << wellProd[WellProdDataType::THP] << ": \n";//wellProd[WellProdDataType::SteadyStatePI] << std::setw(10) << "\n"
}
fix
2019-10-30 11:51:49 -05:00
ss << ":"<< std::setfill ('-') << std::setw (9) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (5) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (12) << ":" << std::setfill ('-') << std::setw (11) << ":" << std::setfill ('-') << std::setw (13) << ":" << std::setfill ('-') << std::setw (9) << ":" << std::setfill ('-') << std::setw (9) << ":" << "\n";
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
Opm::OpmLog::note(ss.str());
}
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
void outputInjectionReport_(const ScalarBuffer& wellInj, const StringBuffer& wellInjNames, const bool forceDisableInjOutput)
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if(forceDisableInjOutput)
return;
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
std::ostringstream ss;
if (wellInjNames[WellInjDataType::WellName].empty()) {
ss << "=================================================== INJECTION REPORT ========================================\n"//===================== \n"
<< ": WELL : LOCATION : CTRL : CTRL : CTRL : OIL : WATER : GAS : FLUID : BHP OR : THP OR :\n"// STEADY-ST II :\n"
<< ": NAME : (I,J,K) : MODE : MODE : MODE : RATE : RATE : RATE : RES.VOL. : CON.PR.: BLK.PR.:\n";// OR POTENTIAL :\n";
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
ss << ": : : OIL : WAT : GAS : SCM/DAY : SCM/DAY : SCM/DAY : RCM/DAY : BARSA : BARSA :\n";// :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
ss << ": : : OIL : WAT : GAS : STB/DAY : STB/DAY : MSCF/DAY : RB/DAY : PSIA : PSIA :\n";// :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_LAB) {
ss << ": : : OIL : WAT : GAS : SCC/HR : SCC/HR : SCC/HR : RCC/HR : ATMA : ATMA :\n";// :\n";
}
ss << "==============================================================================================================\n";//===================== \n";
}
else {
added groups to reports
2019-10-30 08:09:51 -05:00
if (wellInj[WellInjDataType::WellLocationi] < 1) {
ss << std::right << std::fixed << std::setprecision(0) << ":" << std::setw (8) << wellInjNames[WellInjDataType::WellName] << ":" << std::setw(11) << "" << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeOil] << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeWat] << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeGas] << ":" << std::setprecision(1) << std::setw(11) << wellInj[WellInjDataType::OilRate] << ":" << std::setw(11) << wellInj[WellInjDataType::WaterRate] << ":" << std::setw(11)<< wellInj[WellInjDataType::GasRate] << ":" << std::setw(11) << wellInj[WellInjDataType::FluidResVol] << ":" << std::setw(8)<< "" << ":" << std::setw(8)<< "" << ": \n";//wellInj[WellInjDataType::SteadyStateII] << std::setw(10) << "\n"
}
else {
ss << std::right << std::fixed << std::setprecision(0) << ":" << std::setw (8) << wellInjNames[WellInjDataType::WellName] << ":" << std::setw(5) << wellInj[WellInjDataType::WellLocationi] << "," << std::setw(5) << wellInj[WellInjDataType::WellLocationj] << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeOil] << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeWat] << ":" << std::setw(6) << wellInjNames[WellInjDataType::CTRLModeGas] << ":" << std::setprecision(1) << std::setw(11) << wellInj[WellInjDataType::OilRate] << ":" << std::setw(11) << wellInj[WellInjDataType::WaterRate] << ":" << std::setw(11)<< wellInj[WellInjDataType::GasRate] << ":" << std::setw(11) << wellInj[WellInjDataType::FluidResVol] << ":" << std::setw(8)<< wellInj[WellInjDataType::BHP] << ":" << std::setw(8)<< wellInj[WellInjDataType::THP] << ": \n";//wellInj[WellInjDataType::SteadyStateII] << std::setw(10) << "\n"
}
ss << ":--------:-----------:------:------:------:------------:----------:-----------:-----------:--------:--------: \n";//--------------------:\n";
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
Opm::OpmLog::note(ss.str());
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
void outputCumulativeReport_(const ScalarBuffer& wellCum, const StringBuffer& wellCumNames, const bool forceDisableCumOutput)
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
if(forceDisableCumOutput)
return;
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
const Opm::UnitSystem& units = simulator_.vanguard().eclState().getUnits();
std::ostringstream ss;
if (wellCumNames[WellCumDataType::WellName].empty()) {
ss << "=================================================== CUMULATIVE PRODUCTION/INJECTION REPORT =========================================\n"
<< ": WELL : LOCATION : WELL :CTRL: OIL : WATER : GAS : Prod : OIL : WATER : GAS : INJ :\n"
<< ": NAME : (I,J,K) : TYPE :MODE: PROD : PROD : PROD : RES.VOL. : INJ : INJ : INJ : RES.VOL. :\n";
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_METRIC) {
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
ss << ": : : : : MSCM : MSCM : MMSCM : MRCM : MSCM : MSCM : MMSCM : MRCM :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_FIELD) {
ss << ": : : : : MSTB : MSTB : MMSCF : MRB : MSTB : MSTB : MMSCF : MRB :\n";
}
if (units.getType() == Opm::UnitSystem::UnitType::UNIT_TYPE_LAB) {
ss << ": : : : : MSCC : MSCC : MMSCC : MRCC : MSCC : MSCC : MMSCC : MRCC :\n";
}
ss << "====================================================================================================================================\n";
}
else {
added groups to reports
2019-10-30 08:09:51 -05:00
if (wellCum[WellCumDataType::WellLocationi] < 1) {
Update ecloutputblackoilmodule.hh
2019-11-06 06:25:04 -06:00
ss << std::right << std::fixed << std::setprecision(0) << ":" << std::setw (8) << wellCumNames[WellCumDataType::WellName] << ":" << std::setw(11) << "" << ":" << std::setw(8) << wellCumNames[WellCumDataType::WellType] << ":" << std::setw(4) << wellCumNames[WellCumDataType::WellCTRL] << ":" << std::setprecision(1) << std::setw(11) << wellCum[WellCumDataType::OilProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::WaterProd]/1000 << ":" << std::setw(11)<< wellCum[WellCumDataType::GasProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::FluidResVolProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::OilInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::WaterInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::GasInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::FluidResVolInj]/1000 << ": \n";
added groups to reports
2019-10-30 08:09:51 -05:00
}
else {
Update ecloutputblackoilmodule.hh
2019-11-06 06:25:04 -06:00
ss << std::right << std::fixed << std::setprecision(0) << ":" << std::setw (8) << wellCumNames[WellCumDataType::WellName] << ":" << std::setw(5) << wellCum[WellCumDataType::WellLocationi] << "," << std::setw(5) << wellCum[WellCumDataType::WellLocationj] << ":" << std::setw(8) << wellCumNames[WellCumDataType::WellType] << ":" << std::setw(4) << wellCumNames[WellCumDataType::WellCTRL] << ":" << std::setprecision(1) << std::setw(11) << wellCum[WellCumDataType::OilProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::WaterProd]/1000 << ":" << std::setw(11)<< wellCum[WellCumDataType::GasProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::FluidResVolProd]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::OilInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::WaterInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::GasInj]/1000 << ":" << std::setw(11) << wellCum[WellCumDataType::FluidResVolInj]/1000 << ": \n";
added groups to reports
2019-10-30 08:09:51 -05:00
}
ss << ":--------:-----------:--------:----:------------:----------:-----------:-----------:------------:----------:-----------:-----------: \n";
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
Opm::OpmLog::note(ss.str());
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
std::string fipEnumToString_(int i)
{
Some more minor clean-up
2018-01-29 05:23:23 -06:00
typedef typename FipDataType::FipId FipId;
ecloutputblackoilmodule.hh: adapt coding style to the rest of eWoms the style was mostly fine, but there were a few leftovers from the pre-eWoms epoch of that code.
2018-04-30 10:17:22 -05:00
switch(static_cast<FipId>(i))
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
{
Some more minor clean-up
2018-01-29 05:23:23 -06:00
case FipDataType::WaterInPlace: return "WIP";
case FipDataType::OilInPlace: return "OIP";
case FipDataType::GasInPlace: return "GIP";
case FipDataType::OilInPlaceInLiquidPhase: return "OIPL";
case FipDataType::OilInPlaceInGasPhase: return "OIPG";
case FipDataType::GasInPlaceInLiquidPhase: return "GIPL";
case FipDataType::GasInPlaceInGasPhase: return "GIPG";
case FipDataType::PoreVolume: return "PV";
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
}
return "ERROR";
}
fix
2019-10-30 09:27:41 -05:00
std::string WPEnumToString_(int i)
{
typedef typename WellProdDataType::WPId WPId;
switch(static_cast<WPId>(i))
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
case WellProdDataType::WellName: return "WName";
case WellProdDataType::WellLocationi: return "WLi";
case WellProdDataType::WellLocationj: return "WLj";
case WellProdDataType::CTRLMode: return "CTRL";
case WellProdDataType::OilRate: return "OR";
case WellProdDataType::WaterRate: return "WR";
case WellProdDataType::GasRate: return "GR";
case WellProdDataType::FluidResVol: return "FRV";
case WellProdDataType::WaterCut: return "WC";
case WellProdDataType::GasOilRatio: return "GOR";
case WellProdDataType::WatGasRatio: return "WGR";
case WellProdDataType::BHP: return "BHP";
case WellProdDataType::THP: return "THP";
case WellProdDataType::SteadyStatePI: return "SteadyStatePI";
fix
2019-10-30 09:27:41 -05:00
}
return "ERROR";
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
}
std::string WIEnumToString_(int i)
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
typedef typename WellInjDataType::WIId WIId;
switch(static_cast<WIId>(i))
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
case WellInjDataType::WellName: return "WName";
case WellInjDataType::WellLocationi: return "WLi";
case WellInjDataType::WellLocationj: return "WLj";
case WellInjDataType::CTRLModeOil: return "CTRLo";
case WellInjDataType::CTRLModeWat: return "CTRLw";
case WellInjDataType::CTRLModeGas: return "CTRLg";
case WellInjDataType::OilRate: return "OR";
case WellInjDataType::WaterRate: return "WR";
case WellInjDataType::GasRate: return "GR";
case WellInjDataType::FluidResVol: return "FRV";
case WellInjDataType::BHP: return "BHP";
case WellInjDataType::THP: return "THP";
case WellInjDataType::SteadyStateII: return "SteadyStateII";
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
return "ERROR";
}
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
std::string WCEnumToString_(int i)
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
{
typedef typename WellCumDataType::WCId WCId;
switch(static_cast<WCId>(i))
{
çonverted tabs to spaces
2019-11-01 06:56:20 -05:00
case WellCumDataType::WellName: return "WName";
case WellCumDataType::WellLocationi: return "WLi";
case WellCumDataType::WellLocationj: return "WLj";
case WellCumDataType::WellType: return "WType";
case WellCumDataType::WellCTRL: return "WCTRL";
case WellCumDataType::OilProd: return "OP";
case WellCumDataType::WaterProd: return "WP";
case WellCumDataType::GasProd: return "GP";
case WellCumDataType::FluidResVolProd: return "FRVP";
case WellCumDataType::OilInj: return "OI";
case WellCumDataType::WaterInj: return "WI";
case WellCumDataType::GasInj: return "GI";
case WellCumDataType::FluidResVolInj: return "FRVI";
changed as revised by OPMUsers
2019-10-26 06:13:44 -05:00
}
return "ERROR";
}
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
Cell Data: Don't Warn About Supported Output Directives This commit ensures that we don't issue misleading warning about output creation directives such as BASIC FREQ RESTART that could be, and typically are, put into keywords like RPTRST and RPTSCHED. Previously, we would interpret these directives as names of cell vectors and we would get very confusing diagnostic messages of the form Keyword 'BASIC' is unhandled for output to file which suggests that 'BASIC' is unsupported despite the output system supporting the directive for quite some time. While here, also apply the same reasoning to the vectors we always output such as SWAT, SGAS, or PRESSURE.
2020-03-25 05:15:35 -05:00
bool isOutputCreationDirective_(const std::string& keyword) const
{
return (keyword == "BASIC") || (keyword == "FREQ")
|| (keyword == "RESTART") // From RPTSCHED
|| (keyword == "SAVE") || (keyword == "SFREQ"); // Not really supported
}
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
const Simulator& simulator_;
Some clean-up Dont' output rs, rv, somax and hysteresis if they are not enabled. We no longer output the same as legacy.
2018-01-22 06:38:50 -06:00
bool outputFipRestart_;
Fix slow code Avoid checking summaryConfig at every cell Add temporary fix around issue with RPTSCHED in restartConfig
2018-01-23 02:46:56 -06:00
bool computeFip_;
ebos: add possibility to disable the FIP output even if requested in the deck this significantly reduces the noise level of debug files if the FIP values are not of interest.
2018-07-26 05:31:32 -05:00
bool forceDisableFipOutput_;
Pass block values to OPM-output
2018-01-19 08:45:42 -06:00
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ScalarBuffer saturation_[numPhases];
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ScalarBuffer oilPressure_;
ScalarBuffer temperature_;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ScalarBuffer gasDissolutionFactor_;
implement two-phase blackoil simulations
2016-11-15 11:00:48 -06:00
ScalarBuffer oilVaporizationFactor_;
Revert "rename {gil,gas}FormationVolumeFactor to \1FormationFactor" This reverts commit 3e247b059f20b1d04f47bd832719a99df8b9ac26. After a more thorough investigation, the cannonical name of these quantities turned out to be "* formation volume factor"...
2015-01-26 04:55:37 -06:00
ScalarBuffer gasFormationVolumeFactor_;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ScalarBuffer saturatedOilFormationVolumeFactor_;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
ScalarBuffer oilSaturationPressure_;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ScalarBuffer rs_;
ScalarBuffer rv_;
ScalarBuffer invB_[numPhases];
ScalarBuffer density_[numPhases];
ScalarBuffer viscosity_[numPhases];
ScalarBuffer relativePermeability_[numPhases];
ScalarBuffer sSol_;
ScalarBuffer cPolymer_;
Use foam module.
2019-04-26 03:38:37 -05:00
ScalarBuffer cFoam_;
addd simple salt/brine implementation
2019-11-07 02:39:42 -06:00
ScalarBuffer cSalt_;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ScalarBuffer soMax_;
ScalarBuffer pcSwMdcOw_;
ScalarBuffer krnSwMdcOw_;
ScalarBuffer pcSwMdcGo_;
ScalarBuffer krnSwMdcGo_;
Output and read PPCW when SWATINIT is enabled
2018-09-07 03:19:00 -05:00
ScalarBuffer ppcw_;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
ScalarBuffer bubblePointPressure_;
ScalarBuffer dewPointPressure_;
Add support for water induced compaction using ROCKCOMP, ROCK2D, ROCK2DTR, ROCKWNOD and OVERBURD
2019-03-12 09:51:41 -05:00
ScalarBuffer rockCompPorvMultiplier_;
ScalarBuffer rockCompTransMultiplier_;
ScalarBuffer swMax_;
ScalarBuffer overburdenPressure_;
ScalarBuffer minimumOilPressure_;
Add ecl output from opm-output The new eclwriter output and restart using the eclIO from opm-output All tests in opm-simulator pass and MPI restart works. Standard initialization is done prior to restart in order to compute correct initial thpressure values. This is not necessary if thpressures are written to the restart file TODO: Some trivial fields are written out in order mimic legacy code, this should be cleaned up TODO: Output of wells, FIP, NNC is still done in opm-simulators. This should be moved later.
2017-12-07 05:38:48 -06:00
std::vector<int> failedCellsPb_;
std::vector<int> failedCellsPd_;
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
std::vector<int> fipnum_;
Some more minor clean-up
2018-01-29 05:23:23 -06:00
ScalarBuffer fip_[FipDataType::numFipValues];
Add FIP, INIT and EGRID output Writes INIT and EGRID files initially Adds Fip to summary output if required. Output Fip to log (.PRT) if Opm-logger is set
2018-01-12 08:32:43 -06:00
ScalarBuffer origTotalValues_;
Some more minor clean-up
2018-01-29 05:23:23 -06:00
ScalarBuffer origRegionValues_[FipDataType::numFipValues];
Minor clean-up after review
2018-01-29 01:43:42 -06:00
ScalarBuffer hydrocarbonPoreVolume_;
ScalarBuffer pressureTimesPoreVolume_;
ScalarBuffer pressureTimesHydrocarbonVolume_;
Split the PackUnPack class into PackUnPackBlockData, PackUnPackCellData and PackUnPackWellData
2018-02-09 05:57:34 -06:00
std::map<std::pair<std::string, int>, double> blockData_;
rename completion to connection
2018-06-28 08:49:45 -05:00
std::map<size_t, Scalar> oilConnectionPressures_;
std::map<size_t, Scalar> waterConnectionSaturations_;
std::map<size_t, Scalar> gasConnectionSaturations_;
Add ecl tracer model reads tracer input from deck, solves tracer equation fully implicit as a post processing step in endTimeStep tested on a simple modified SPE1CASE1 deck and compared with eclipse TODO: restart and parallel
2018-11-14 06:10:01 -06:00
std::vector<ScalarBuffer> tracerConcentrations_;
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
};
changed: namespace Ewoms -> namespace Opm
2019-09-05 10:04:39 -05:00
} // namespace Opm
move everything which is ECL specific to applications/ebos this helps to keep the core blackoil model code lean and mean and it is also less confusing for newbies because the ECL blackoil simulator is not a "test" anymore. in case somebody wonders, "ebos" stands for "&eWoms &Black-&Oil &Simulator". I picked this name because it is short, a syllable, has not been taken by anything else (as far as I know) and "descriptive" names are rare for programs anyway: everyone who does not yet know about 'git' or 'emacs' and tells me that based on their names they must be a source-code managment system and an editor gets a crate of beer sponsored by me!
2014-11-28 05:58:48 -06:00
#endif
Reference in New Issue Copy Permalink