2009 Stream 3 : New features =

Last edited Timestamp?

---

S3.1 : from the community

It is of paramount importance to incorporate features developed in the NEMO community and identified as priorities by the Developer’s Committee in order to improve the system and encourage further developments. The 7 actions are:

S3.1-a — Tidal mixing scheme

Status: zdftmx.F90 added in the trunk and validated in orca2-lim (NEMO-paris), see ticket #425 . It includes both a bottom intensified vertical mixing in area of internal tide generation (Simmons OM 2005, Bessière PhD 2007) and a intensification throughout the water column in the Indonesian area (Koch-Larrouy et al. 2007, 2008a,b).
Pending: on-line comments and minor optimisation (see detail in ticket #457), bug correction in the call (see ticket #462) , and paper documentation #458

S3.1-b — Three bands light penetration + ocean colour

Status: The light penetration developed by Lengaigne et al. (CD07, GRL08) has been introduced in traqsr.F90 in the trunk. It include an optional coupling to PISCES or the use of ocean color data. It has been validated in orca2-lim-pisces (NEMO-paris), see ticket #428. Note that new namelist parameters have been introduced to control it.
Pending: paper documentation, see ticket #459

S3.1-c — Diurnal cycle on solar radiation

Status: Bernie et al. (2007) distribution of the daily solar radiation onto the diurnal cycle has been implemented in an on-line module and validated in orca05 simulation on NEMO v2 by S. Masson and co.
Pending: its incorporation in v3. To be done by NEMO-Paris

S3.1-d Conservative Leap-Frog + Asselin filter

Status: ongoing work (NEMO-Paris). See detail 2009WP/2009Stream3/LeapFrogC

S3.1-e — Observations operator code and increments application

Status: to be done by NEMO-Exeter.

S3.1-f — Ice-shelf melting parameterisation

Status: See Wang & Beckmann (AG 2007), PhD of P. Mathiot (2008). To be done by NEMO-Paris. The key issue here is that the existing code does not include the volume flux associated with the ice-shelf melting. A major change have to be performed prior to its incorporation in the reference version.

S3.1-g — Parameterisation of katabatic winds

Status: See PhD of P. Mathiot (2008). to be done by NEMO-Paris.

---

S3.2 : from NEMO team (short-term)

It consists of the short-term developments associated with the priorities of each consortium members that will be incorporated into the system, as they can be useful for the user community. The actions of the 8 sub-stream are:

S3.2-a — New forcing

Tidal potential, atmospheric pressure) and river input scheme (continuation & end)

S3.2-b — Revisit of mass and salt fluxes between ice, ocean and air

continuation & end.

S3.2-c — Embedding sea-ice inside the ocean

continuation & end.

S3.2-d — New advection schemes

PPM, Quickest and Prather.

S3.2-e — Lateral ocean physics

isoneutral operator, skew flux, Visbeck coefficients, slope limitation, Aiki formulation. (continuation & end)

S3.2-f — Ice improvements

momentum forcing, snow-ice formation, EVP rheology in LIM2.

S3.2-g — Improved s-coordinate option

hydrostatic pressure gradient computation, POLCOM-style s-coordinate. (continuation & end)

S3.2-h — Vertical physics

improvements of TKE scheme, addition of a Generic length scale model

--- zdftke2.F90 : validation in global condiguration, update of the paper documentation.
Status: zdftke2.F90 has been provided by NEMO-Toulouse, with its paper documentation, and put in release v3.1 (except the doc).
Problems have been found with in orca2-lim in both forced and coupled mode (NEMO-paris).
Pending issues:
(1) add the paper documentation in the trunk (DOC).
(2) bug corrections:

correct the prandl number computation in order to avoid a floating point exception, i.e. replace in tke_mlmc routine:

                  zri = MAX( rn2(ji,jj,jk), 0. ) / (zcoef * (zdku + zdkv + 1.e-20) / avmt(ji,jj,jk)) ! local Richardson number

by

                  zri = MAX( rn2(ji,jj,jk), 0. ) * avmt(ji,jj,jk) / (zcoef * (zdku + zdkv + 1.e-20) )   ! local Richardson number

partial step topography: add in step a missing call to zpshde when zdftke2 is used with ln_dynhpg_imp=FALSE.

(3) change the restart strategy to allow coupled simulation
In the current version, the restart of tke2 requires to compute en at nitend+1. In coupled mode, the forcing of the next coupled period cannot be known. The restart must be change switch from (en , dissl and avt at nitend+1) to (en, dissl, avt, avm at nitend) (ongoing)

--- zdfgls.F90 : addition of the Generic Length scale scheme (Warner et al. Ocean Modelling, 8, 81-113, 2005).
Status: on-going work (NEMO-Toulouse).

---

S3.3 : from NEMO team (long-term)

It consists of the short-term developments associated with the priorities of each consortium members that will be incorporated into the system, as they can be useful for the user community. The actions of the 8 sub-stream are:

S3.3-a — Revisit of the time stepping in TOP

(2 years, start in 2009)

S3.3-b — Restore the variable volume option (key_vvl)

(2 years, started in 2008)
Status: ongoing work (NEMO-Paris). See detail 2009WP/2009Stream3/VVL

S3.3-c — Toward a generalisation of the split-explicit free surface

(2 years, start in 2009)

S3.3-d — Inclusion of the effect of waves in the marginal ice zone

(2 years, start in 2009)

S3.3-e — Interface to the OASIS 4 coupler

(2-3 years, ongoing, started in 2008)

S3.3-f — AGRIF and sea-ice

(from 0 order to LIM2 and 3 fully AGRIFed) (3 years, started in 2008)

S3.3-g — Development of TAM

(ANR VODA), with a major contribution from the LJK team (3 years, started in 2008)

S3.3-h — High parallelisation on new plateforms

(GPU graphic cards, OpenMP etc …, start in 2009)

---