Skip to content

Commit

Permalink
Fix to the halo issue in the variance partitioning code (#1171)
Browse files Browse the repository at this point in the history
- fixes #1170
  • Loading branch information
guillaumevernieres authored Jun 13, 2024
1 parent 9ffce97 commit f5bbec0
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions utils/soca/gdas_soca_diagb.h
Original file line number Diff line number Diff line change
Expand Up @@ -300,12 +300,12 @@ namespace gdasapp {
}

// Update the layer thickness halo
xbFs["hocn"].haloExchange();
nodeColumns.haloExchange(xbFs["hocn"]);

// Loop through variables
for (auto & var : configD.socaVars.variables()) {
// Update the halo
xbFs[var].haloExchange();
nodeColumns.haloExchange(xbFs[var]);

// Skip the layer thickness variable
if (var == "hocn") {
Expand Down Expand Up @@ -361,7 +361,7 @@ namespace gdasapp {
// Horizontal averaging
for (int iter = 0; iter < configD.niterHoriz; ++iter) {
// Update the halo points
bkgErrFs[var].haloExchange();
nodeColumns.haloExchange(bkgErrFs[var]);
auto stdDevBkg = atlas::array::make_view<double, 2>(bkgErrFs[var]);

// Loops through nodes and levels
Expand Down

0 comments on commit f5bbec0

Please sign in to comment.