| Version 1 (modified by davestorkey, 3 years ago) (diff) |
|---|
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.
Summary
| Action | ENHANCE-01_davestorkey_fix_3D_momentum_trends |
|---|---|
| PI(S) | Dave Storkey |
| Digest | Fix 3D momentum trends diagnostics to work properly with split-explicit timestepping |
| Dependencies | None |
| Branch | source:/NEMO/branches/{YEAR}/dev_r{REV}_{ACTION_NAME} |
| Previewer(s) | Gurvan Madec |
| Reviewer(s) | TBD |
| Ticket | #2344 |
Description
There are various issues with the 3D momentum trends diagnostics (ln_dyn_trd=T) when using the split-explicit solution of the momentum equation (dynspg_ts - default from NEMO 4.0 onwards):
- The "SPG" trend is not the trend due to the surface pressure gradient. It is the change in the velocities due to the solution of the barotropic part of the equations in dyn_spg_ts.
- More generally, for trends which are allowed to vary over the barotropic subcycling (PVO, BFR), there is a small inaccuracy because currently this variation is not taken into account.
- At various points in the code the barotropic part of the velocities is subtracted or added back in (eg. at the beginning of dyn_spg_ts, or at the beginning of dyn_zdf for the implicit bottom friction option). This subtraction/addition of the barotropic component is sometimes incorrectly included in the trend associated with that module.
These problems mean that it is currently not possible to close the momentum budget, as already noted in tickets #1584 and #1984. The 3D momentum trends are used to calculate various other diagnostics such as 2D vorticity trends (ln_vor_trd=T) or kinetic energy trends (ln_KE_trd=T) so the problems also affect these derived diagnostics.
Implementation
- Momentum trends diagnostics should correspond to terms in the primitive equations (rather than NEMO modules).
- We should have a small number of 3D trends which are sufficient to close the budget: a) Each 3D trend should include the barotropic and baroclinic part. (We can output the barotropic part separately if we want). Where the barotropic part of the trend varies over the subcycling this part can be calculated accurately using the time-averaging weights already available in dynspg_ts.F90. b) The horizontal pressure gradient trend (HPG) should include the surface pressure gradient (SPG). The SPG trend can be output as an extra diagnostics not required to close the budget. c) The vertical diffusion trend (ZDF) already includes the contribution of the wind stress, top stress (in ice shelf cavities) and implicit bottom stress, but it should also include the contribution from the explicit bottom stress terms. The explicit bottom stress can still be output as an extra diagnostic not required to close the budget.
- We should make sure that the subtraction or addition of the barotropic component of velocities or velocity trends is not included in the momentum trend diagnostics. This particularly applies to the ZDF trend. (This is similar to a fix proposed but not implemented by Christoph in ticket #1584).
Documentation updates
...
Preview
...
Tests
...
Review
...
