Opened 6 years ago

Closed 3 years ago

#1350 closed Enhancement (fixed)

dev CNRS-0 --- new version of tranpc (non-penetrative convection) based local stability

Reported by: gm Owned by: nemo
Priority: low Milestone:
Component: OCE Version: release-3.6
Severity: Keywords:
Cc: Review:
MP ready?:


npc, the non-penetrative convection algorithm, has been designed to restore the static stability of the water column at each model time step (Madec et al. JPO 1991). Its implementation in OPA was using potential density to calculate the static stability, a limitation of the algorithm that make it useless when using a full equation of state. In addition, the passage from OPA8 (Fortran 77) to OPA9 (NEMO, Fortran 90) introduces bugs in the algorithm that haven't been fixed (see #1092, #1326).

Here, a new version of the algorithm has been developed. It is based on the same principle as the original one, but is now able to handle any equation of state (i.e. working on locally computed Brunt-Vaisala frequency). To that purpose, the algorithm is now using alpha and beta (thermal expansion and haline contraction coeff.) to compute the stability of the water column instead of gradient of potential density.

The algorithm has been implemented by Laurent Brodeau (brodeau@…).
He intensively tested in ORCA2-LIM2 and ORCA1-LIM2 configurations.

The new algorithm has 2 main properties :

1- it is designed for any equation of state
2- it is working !

Note that since the algorithm is using alpha and beta, it is usable only in the new TEOS10 branche (see #1294).

Commit History (1)


#1350 : new version of tranpc.F90 based local static stability

Change History (1)

comment:1 Changed 3 years ago by lovato

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

The developments of this ticket are included in the NEMO V3.6_stable (r5550) and trunk (r5386).

Note: See TracTickets for help on using tickets.