Skip to content

Commit

Permalink
Merge branch 'feature/insitu-concat' of https://github.com/NOAA-EMC/G…
Browse files Browse the repository at this point in the history
…DASApp into feature/insitu-concat
  • Loading branch information
AndrewEichmann-NOAA committed Feb 21, 2025
2 parents 3876928 + 055a80f commit 2efb0c3
Showing 1 changed file with 16 additions and 0 deletions.
16 changes: 16 additions & 0 deletions utils/soca/gdas_ens_handler.h
Original file line number Diff line number Diff line change
Expand Up @@ -170,6 +170,22 @@ namespace gdasapp {
// Set variables to zero if specified in the configuration
postProcIncr.setToZero(incr);

// Optionally apply inflation
if (fullConfig.has("ensemble inflation.value")) {
const double inflation = fullConfig.getDouble("ensemble inflation.value");
mom6_incr *= inflation;
oops::Log::info() << "incr after scalar inflation " << i << ":"
<< mom6_incr << std::endl;
}
if (fullConfig.has("ensemble inflation.field")) {
soca::Increment weight(geomOut, mom6_incr.variables(), mom6_incr.validTime());
const eckit::LocalConfiguration weightConf(fullConfig, "ensemble inflation.field");
weight.read(weightConf);
mom6_incr.schur_product_with(weight);
oops::Log::info() << "incr after field inflation " << i << ":"
<< mom6_incr << std::endl;
}

// Save the increments used to initialize the ensemble forecast
result = postProcIncr.save(mom6_incr, i+1);

Expand Down

0 comments on commit 2efb0c3

Please sign in to comment.