- Timestamp:
- 2016-10-19T12:54:45+02:00 (8 years ago)
- Location:
- branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/CONFIG/SHARED/field_def.xml
r7011 r7040 782 782 <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)" unit="W/s^3" /> 783 783 <field id="KE" long_name="kinetic energy: u(n)*u(n+1)/2" unit="W/s^2" /> 784 <field id="ketrd_eiv" long_name="Kinetic energy trend from mesoscale eddy advection" unit="W/s^2" /> 784 785 785 786 <!-- variables available with ln_PE_trd --> -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90
r6910 r7040 240 240 CALL dia_ptr_ohst_components( jp_sal, 'eiv', z3d ) 241 241 ENDIF 242 243 IF( ln_trddyn .AND. ln_KE_trd ) CALL trd_dyn(u_eiv, v_eiv, jptrd_eivke, kt ) 242 244 # endif 243 245 -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trd_oce.F90
r6427 r7040 56 56 ! 57 57 ! !!!* Momentum trends indices 58 INTEGER, PUBLIC, PARAMETER :: jptot_dyn = 1 5!: Total trend nb: change it when adding/removing one indice below58 INTEGER, PUBLIC, PARAMETER :: jptot_dyn = 16 !: Total trend nb: change it when adding/removing one indice below 59 59 ! =============== ! 60 60 INTEGER, PUBLIC, PARAMETER :: jpdyn_hpg = 1 !: hydrostatic pressure gradient … … 73 73 INTEGER, PUBLIC, PARAMETER :: jpdyn_spgflt = 14 !: filter contribution to surface pressure gradient (spg_flt) 74 74 INTEGER, PUBLIC, PARAMETER :: jpdyn_spgexp = 15 !: explicit contribution to surface pressure gradient (spg_flt) 75 INTEGER, PUBLIC, PARAMETER :: jpdyn_eivke = 16 !: K.E trend from Gent McWilliams scheme 75 76 ! 76 77 !!---------------------------------------------------------------------- -
branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90
r6427 r7040 27 27 USE lib_mpp ! MPP library 28 28 USE wrk_nemo ! Memory allocation 29 USE ldfslp ! Isopycnal slopes 29 30 30 31 IMPLICIT NONE … … 42 43 # include "domzgr_substitute.h90" 43 44 # include "vectopt_loop_substitute.h90" 45 # include "ldfeiv_substitute.h90" 46 44 47 !!---------------------------------------------------------------------- 45 48 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 102 105 ENDIF 103 106 ! 104 zke(:,:,jpk) = 0._wp105 zke(1,:, : ) = 0._wp106 107 zke(:,1, : ) = 0._wp 107 108 DO jk = 1, jpkm1 … … 192 193 CALL ken_p2k( kt , zke ) 193 194 CALL iom_put( "ketrd_convP2K", zke ) ! conversion -rau*g*w 195 CASE( jpdyn_eivke ) 196 ! CMIP6 diagnostic tknebto = tendency of KE from 197 ! parameterized mesoscale eddy advection 198 ! = vertical_integral( k (N S)^2 ) rho dz 199 ! rho = reference density 200 ! S = isoneutral slope. 201 ! Most terms are on W grid so work on this grid 202 zke2d(:,:) = 0._wp 203 DO jk = 1,jpk 204 DO ji = 1,jpi 205 DO jj = 1,jpj 206 zke2d(ji,jj) = zke2d(ji,jj) + rau0 * fsaeiw(ji, jj, jk) & 207 & * ( wslpi(ji, jj, jk) * wslpi(ji,jj,jk) & 208 & + wslpj(ji, jj, jk) * wslpj(ji,jj,jk) ) & 209 & * bn2(ji,jj,jk) * fse3w(ji, jj, jk) 210 ENDDO 211 ENDDO 212 ENDDO 213 iom_put("ketrd_eiv", zke2d) 194 214 ! 195 215 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.