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.
2020WP/ASINTER-02_emanuelaclementi_Waves – NEMO

Name and subject of the action

Last edition: Wikinfo(changed_ts)? by Wikinfo(changed_by)?

The PI is responsible to closely follow the progress of the action, and especially to contact NEMO project manager if the delay on preview (or review) are longer than the 2 weeks expected.

  1. Summary
  2. Preview
  3. Tests
  4. Review


Action ASINTER-02_emanuelaclementi_Waves
PI(S) Emanuela Clementi
Digest New ocean-wave mixing compliant with bulks and OBL mixing
Dependencies #2155
Branch source:/NEMO/branches/{YEAR}/dev_r{REV}_{ACTION_NAME}
Previewer(s) Gurvan Madec, Rachid Benshila?
Reviewer(s) Rachid Benshila, Gurvan Madec
Ticket #2339


This action is an activity of IMMERSE Task5.2 and focuses on the upgrade the ocean-wave coupling representation.

The final aim is to include an additional wave coupling development by mid 2020, which consists in modified
vertical mixing due to breaking waves: the energy flux explicitly evaluated by a wave model
(through the dissipation term in the wave energy balance equation) will be used as a new source of vertical mixing.

The first step is to include recent works on enhanced ocean mixing due to breaking waves (Couvelard et al., 2019)
and outcomes from WAVE2NEMO CMEMS service evolution project.

In addition, new wave developments detailed in Ticket #2155
should be revised to be compliant with bulks and OBL mixing in NEMO reference version.


A branch has been created including wave induced terms following Couvelard et al., 2019.

Namelist, ZDF, DYN and SBC routines are impacted:

  • Namelist modifications to move all related wave parameters in namsbc_wave (except ln_wave) and include new parameters.
  • Update sbcmod.F90, sbcwave.F90 , sbccpl.F90, sbccpl.F90, cpl_oasis3.F90 to include additional wave parameters and coupling:
    • Chranock coefficient from wave model
    • wind stress corrected by wave intake
    • TKE surface BC from wave model
    • Bernoulli head, waves' inuced pressure
    • Breivik 2016 profile
    • Vortex force activation
  • Modify sbcblk_algo_ecmwf.F90 to add Charnok coupling
  • Modify dynvor.F90, dynzad.F90 to add vortex focre term:
  • Modify dynspg.F90 to account for depth-independent Bernoulli head
  • Modifications to tke.F90:
    • Langmuir mixing parameterization including the Stokes Drift contribution using Craik-Leibovich velocity scale (cpl_sdrftx) or Surface Stokes drift deduced from surface stress.
    • Modified shear production term in the TKE equation due to the Stokes Drift
    • Add 2 cases to evaluate en, zdiag, zd_lw, zd_up: Dirichlet Boundary Condition using surface TKE flux from waves and without wave coupling
    • Modifications to allow for Neumann boundary condition
    • Mixing length: add evaluation if coupling with waves
    • Modifiy TKE due to surface and internal wave breaking
    • In tke_avn include the mixing length evaluated if waves data are available

Documentation updates

Error: Failed to load processor box
No macro or processor named 'box' found



Error: Failed to load processor box
No macro or processor named 'box' found



Error: Failed to load processor box
No macro or processor named 'box' found

SETTE tests:

Current code is : NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves @ r13852 ( last change @ r13852 )

SETTE validation report generated for :

NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves @ r13852+ (last changed revision)

on ifort_zeus_xios arch file

!!---------------1st pass------------------!!


WGYRE_PISCES_ST run.stat restartability passed : 13852+
WGYRE_PISCES_ST tracer.stat restartability passed : 13852+
WORCA2_ICE_PISCES_ST run.stat restartability passed : 13852+
WORCA2_ICE_PISCES_ST tracer.stat restartability passed : 13852+
WORCA2_OFF_PISCES_ST tracer.stat restartability passed : 13852+
WAMM12_ST run.stat restartability passed : 13852+
WORCA2_SAS_ICE_ST run.stat restartability passed : 13852+
WAGRIF_DEMO_ST run.stat restartability passed : 13852+
WWED025_ST run.stat restartability passed : 13852+
WISOMIP+_ST run.stat restartability passed : 13852+
WOVERFLOW_ST run.stat restartability passed : 13852+
WLOCK_EXCHANGE_ST run.stat restartability passed : 13852+
WVORTEX_ST run.stat restartability passed : 13852+
WICE_AGRIF_ST run.stat restartability passed : 13852+


WGYRE_PISCES_ST run.stat reproducibility passed : 13852+
WGYRE_PISCES_ST tracer.stat reproducibility passed : 13852+
WORCA2_ICE_PISCES_ST run.stat reproducibility passed : 13852+
WORCA2_ICE_PISCES_ST tracer.stat reproducibility passed : 13852+
WORCA2_OFF_PISCES_ST tracer.stat reproducibility passed : 13852+
WAMM12_ST run.stat reproducibility passed : 13852+
WORCA2_SAS_ICE_ST run.stat reproducibility passed : 13852+
WORCA2_ICE_OBS_ST run.stat reproducibility passed : 13852+
WAGRIF_DEMO_ST run.stat reproducibility passed : 13852+
WWED025_ST run.stat reproducibility passed : 13852+
WISOMIP+_ST run.stat reproducibility passed : 13852+
WVORTEX_ST run.stat reproducibility passed : 13852+
WICE_AGRIF_ST run.stat reproducibility passed : 13852+

!----agrif check----!

ORCA2 AGRIF vs ORCA2 NOAGRIF run.stat unchanged - passed : 13852+ 13852+

Validation tests:

The branch validation has been performed in both forced and coupled modes based on the ORCA2_ICE_PISCES reference configuration, for 20 simulation days (see attachment).

Test case

  • Test case name: Not defined yet. Draft name: NEMO-WAVE coupling
  • Test case functionality: This test case should check the correct NEMO-WAVE processes implementation and consists of an adiabatic case for the validation of the Generalized Lagrangian Mean implementation for the wave-current interaction.
  • Test case setup: The adiabatic test case consists of a steady monochromatic wave shoaling from 4 to 6 m depth on a slope without bottom friction. A symmetric bottom sloping is used to allow for periodic boundary conditions in the x-direction. Open boundary conditions are defined in y-directions. The characteristics of the monochromatic waves are a significant wave height of 1.02 m, a wave period of 5.24s and a wave direction propagation in the x-direction. Wave breaking and wave induced mixing are not taken into account. The horizontal resolution is 10 meters in x-direction and 50m in y-direction. The NEMO configuration has a high vertical resolution using 60 levels.
  • Test case verification value: Validation of the Generalized Lagrangian Mean implementation
  • Status of the test case as for now: Started


Error: Failed to load processor box
No macro or processor named 'box' found
  • Is the proposed methodology now implemented? : YES
  • Are the code changes in agreement with the flowchart defined at preview step? : YES
  • Are the code changes in agreement with list of routines and variables as proposed at preview step? YES
  • Is the in-line documentation accurate and sufficient? YES
  • Do the code changes comply with NEMO coding standards? YES
  • Is the development documented with sufficient details for others to understand the impact of the change? YES, to me it's part of Couvelard et al
  • Is the project literature (manual, guide, web, …) now updated or completed following the proposed summary in preview section? Can't say, I can't compiled NEMO documentation for a while now

In summary the branch is suitable for merging

Last modified 3 years ago Last modified on 2020-12-01T23:26:19+01:00

Attachments (1)

Download all attachments as: .zip