OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

remove computeBhpAtThpLimitProd() output from .DBG file

Browse Source
This commit is contained in:
Tor Harald Sandve 2022-09-12 09:29:37 +02:00
parent f36ac3eeed
commit 001d53858e
2 changed files with 24 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
opm/simulators/wells/WellInterfaceGeneric.cpp
View File

@ -459,11 +459,13 @@ bhpMax(const std::function<double(const double)>& fflo,
double high = maxPerfPress + 1.0 * unit::barsa; double high = maxPerfPress + 1.0 * unit::barsa;
double f_low = fflo(low); double f_low = fflo(low);
double f_high = fflo(high); double f_high = fflo(high);
deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() + if constexpr (extraBhpAtThpLimitProdOutput) {
" low = " + std::to_string(low) + deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() +
" high = " + std::to_string(high) + " low = " + std::to_string(low) +
" f(low) = " + std::to_string(f_low) + " high = " + std::to_string(high) +
" f(high) = " + std::to_string(f_high)); " f(low) = " + std::to_string(f_low) +
" f(high) = " + std::to_string(f_high));
}
int adjustments = 0; int adjustments = 0;
const int max_adjustments = 10; const int max_adjustments = 10;
const double adjust_amount = 5.0 * unit::barsa; const double adjust_amount = 5.0 * unit::barsa;
@ -512,12 +514,14 @@ bhpMax(const std::function<double(const double)>& fflo,
return std::nullopt; return std::nullopt;
} }
} }
deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() + if constexpr (extraBhpAtThpLimitProdOutput) {
" low = " + std::to_string(low) + deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() +
" high = " + std::to_string(high) + " low = " + std::to_string(low) +
" f(low) = " + std::to_string(f_low) + " high = " + std::to_string(high) +
" f(high) = " + std::to_string(f_high) + " f(low) = " + std::to_string(f_low) +
" bhp_max = " + std::to_string(bhp_max)); " f(high) = " + std::to_string(f_high) +
" bhp_max = " + std::to_string(bhp_max));
}
return bhp_max; return bhp_max;
} }
@ -537,11 +541,13 @@ bisectBracket(const std::function<double(const double)>& eq,
double eq_high = eq(high); double eq_high = eq(high);
double eq_low = eq(low); double eq_low = eq(low);
const double eq_bhplimit = eq_low; const double eq_bhplimit = eq_low;
deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() + if constexpr (extraBhpAtThpLimitProdOutput) {
" low = " + std::to_string(low) + deferred_logger.debug("computeBhpAtThpLimitProd(): well = " + this->name() +
" high = " + std::to_string(high) + " low = " + std::to_string(low) +
" eq(low) = " + std::to_string(eq_low) + " high = " + std::to_string(high) +
" eq(high) = " + std::to_string(eq_high)); " eq(low) = " + std::to_string(eq_low) +
" eq(high) = " + std::to_string(eq_high));
}
if (eq_low * eq_high > 0.0) { if (eq_low * eq_high > 0.0) {
// Failed to bracket the zero. // Failed to bracket the zero.
// If this is due to having two solutions, bisect until bracketed. // If this is due to having two solutions, bisect until bracketed.

2
opm/simulators/wells/WellInterfaceGeneric.hpp
View File

@ -51,6 +51,8 @@ class Schedule;
class WellInterfaceGeneric { class WellInterfaceGeneric {
public: public:
static constexpr bool extraBhpAtThpLimitProdOutput = false;
WellInterfaceGeneric(const Well& well, WellInterfaceGeneric(const Well& well,
const ParallelWellInfo& parallel_well_info, const ParallelWellInfo& parallel_well_info,
const int time_step, const int time_step,