Changeset 7040


Ignore:
Timestamp:
2016-10-19T12:54:45+02:00 (4 years ago)
Author:
timgraham
Message:

Added tnkebto diagnostic = KE trend from mesoscale eddy advection

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  
    782782      <field id="ketrd_convP2K" long_name="ke-trend: conversion (potential to kinetic)"      unit="W/s^3"                        /> 
    783783      <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"                        />    
    784785 
    785786      <!-- variables available with ln_PE_trd --> 
  • branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90

    r6910 r7040  
    240240       CALL dia_ptr_ohst_components( jp_sal, 'eiv', z3d ) 
    241241    ENDIF 
     242 
     243    IF( ln_trddyn .AND. ln_KE_trd ) CALL trd_dyn(u_eiv, v_eiv, jptrd_eivke, kt ) 
    242244# endif   
    243245 
  • branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trd_oce.F90

    r6427 r7040  
    5656   ! 
    5757   !                                                  !!!* Momentum trends indices 
    58    INTEGER, PUBLIC, PARAMETER ::   jptot_dyn  = 15     !: Total trend nb: change it when adding/removing one indice below 
     58   INTEGER, PUBLIC, PARAMETER ::   jptot_dyn  = 16     !: Total trend nb: change it when adding/removing one indice below 
    5959   !                               ===============     !   
    6060   INTEGER, PUBLIC, PARAMETER ::   jpdyn_hpg  =  1     !: hydrostatic pressure gradient  
     
    7373   INTEGER, PUBLIC, PARAMETER ::   jpdyn_spgflt  = 14  !: filter contribution to surface pressure gradient (spg_flt) 
    7474   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 
    7576   ! 
    7677   !!---------------------------------------------------------------------- 
  • branches/UKMO/v3_6_extra_CMIP6_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90

    r6427 r7040  
    2727   USE lib_mpp        ! MPP library 
    2828   USE wrk_nemo       ! Memory allocation 
     29   USE ldfslp         ! Isopycnal slopes 
    2930 
    3031   IMPLICIT NONE 
     
    4243#  include "domzgr_substitute.h90" 
    4344#  include "vectopt_loop_substitute.h90" 
     45#  include "ldfeiv_substitute.h90" 
     46 
    4447   !!---------------------------------------------------------------------- 
    4548   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    102105      ENDIF 
    103106      ! 
    104       zke(:,:,jpk) = 0._wp 
    105       zke(1,:, : ) = 0._wp 
    106107      zke(:,1, : ) = 0._wp 
    107108      DO jk = 1, jpkm1 
     
    192193                    CALL ken_p2k( kt , zke ) 
    193194                      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) 
    194214         ! 
    195215      END SELECT 
Note: See TracChangeset for help on using the changeset viewer.