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