id summary reporter owner description type status priority milestone component version severity resolution keywords cc
1963 CMIP6 / OMIP diagnostics of KE lateral/vertical dissipation and EKE to PE conversion when using EIV gm systeam "In OMIP/ CMIP6 diag., the dissipation of Kinetic Energy due to horizontal and vertical mixing is demanded.
It is not possible to simply use the build-in KE-trend diagnostics to compute them.
Indeed, in KE-trend diagnostics we simply evaluate Un*U-mixing trend. Such a term represent the combination of both diffusive and dissipation trends of KE (see NEMO doc appendix C). In case of a laplacian mixing with a coefficient K, the dissipation trends of KE is -K*[grad(U)]**2, NOT u.div[K*grad(U)]
In NEMO, I modify the following modules:
M NEMO/OPA_SRC/DYN/dynldf_bilap.F90
M NEMO/OPA_SRC/DYN/dynldf_lap.F90
M NEMO/OPA_SRC/DYN/dynzdf_imp.F90
so that they output the vertical sum of KE dissipation trend as demanded by OMIP. There is no need to activate the KE-trend diag (ln_KE_trd=.true.) to obtains theses fields (which slightly speed up the code)
In addition, I also correct a misplaced mixing coefficient in Smagorinski case in dynldf_bilap.F90 (an error which concerns only the v3.6, not the trunk)
A second point concerns the conversion of eddy kinetic energy (EKE) into potential energy due to EIV parameterization (Gent & McWilliams param)
No error here, but misplaced diagnostics. We don't need anymore to activated ln_KE_trd=T to obtain this diag.
I also change the name as it is a conversion from '''E'''KE to PE, not form KE...
To do that I modify the following modules:
M NEMO/OPA_SRC/TRA/traadv_eiv.F90
M NEMO/OPA_SRC/TRD/trd_oce.F90
M NEMO/OPA_SRC/TRD/trdken.F90
I introduce new variables in field_def.xml file which output the required fields:
{{{
}}}
NB this also concerns the '''trunk''' version
I will update the v3.6_stable version tomorrow (on Monday October the 16th) as well as the trunk.
Gurvan" Defect closed normal OCE v3.6 minor fixed CMIP6 diagnostics