New URL for NEMO forge!

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#1119 (Reproducibility using enveloping bathymetry) – NEMO

Opened 11 years ago

Closed 7 years ago

#1119 closed Bug (wontfix)

Reproducibility using enveloping bathymetry

Reported by: jamesharle Owned by: jchanut
Priority: low Milestone:
Component: OCE Version: trunk
Severity: Keywords: DOM s_sigma


Currently in the iterative loop on ln 1173 of domzgr.F90 (R3909 trunk) the land points are included in the bathymetric smoothing. This may not be an issue if a square processor map domain is retained i.e. jpni*jpnj = jpnij. However, if jpni*jpnj >jpnij then the number of land cells that participate in the smoothing could vary depending on the jpni and jpnj chosen. An example of the AMM12 domain is shown in the attached pdf + a full summary of the issues.

Commit History (1)


Solution to ticket #1119: reproducibility using enveloping bathymetry.

Attachments (1)

domzgr_final.pdf (430.7 KB) - added by jamesharle 11 years ago.
Summary of enveloping bathymetry issue + solution

Download all attachments as: .zip

Change History (9)

Changed 11 years ago by jamesharle

Summary of enveloping bathymetry issue + solution

comment:1 Changed 11 years ago by jamesharle

  • Resolution set to fixed
  • Status changed from new to closed

bugfix done in r3926
tested for hbatt and mbathy reproducibility in AMM12 8x4, 4x8 and with removal of land processors.

comment:2 Changed 11 years ago by jamesharle

  • Resolution fixed deleted
  • Status changed from closed to reopened

Although this fix solves several of the issues, there still remains one problem to fix. A bug in the new code was also found to produce NaNs? when scaling factors are calculated. Ticket reopened and old code re-inserted.

comment:3 Changed 11 years ago by jchanut

Why not assuming that we read enveloping bathymetry in a netcdf file (and move the computation in an external pre-processing program) ? This would at least be useful for Agrif in order to ensure bathymetry matching at open boundaries.

comment:4 Changed 11 years ago by clevy

  • Owner changed from NEMO team to jchanut
  • Status changed from reopened to new

comment:5 Changed 11 years ago by rfurner

Hi All,

Apologies for the confusion on this, its partly my fault.  James has come up with a new fix which entirely resolves the problem.  However, it requires new boundaries to be created for the AMM12 configuration.  In the interest of input clarity I thought it best to therefore add this bug fix as part of the merge, so the 3.5 Amm12 inputs remain valid for all of 3.5 and a new set of inputs will be available for 3.6. This isnt ideal from a developers perspective, but I think the alternative is so user unfriendly (having various inputs for different revisions of the same release) that its justified here...

I hope thats ok,


comment:6 Changed 8 years ago by nicolasmartin

  • Keywords DOM added; domzgr removed

comment:7 Changed 8 years ago by nicolasmartin

  • Keywords s_sigma added; sigma removed

comment:8 Changed 7 years ago by clevy

  • Resolution set to wontfix
  • Status changed from new to closed

Suppressing tickets opened on the trunk since more than 2 years... cleaning started!

Note: See TracTickets for help on using tickets.