Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/nemo_v3_6_STABLE/NEMOGCM/NEMO/OPA_SRC/TRD/trdken.F90

    r7564 r8627  
    195195                    CALL ken_p2k( kt , zke ) 
    196196                      CALL iom_put( "ketrd_convP2K", zke )     ! conversion -rau*g*w 
    197         CASE( jpdyn_eivke ) 
    198             ! CMIP6 diagnostic tknebto = tendency of KE from 
    199             ! parameterized mesoscale eddy advection 
    200             ! = vertical_integral( k (N S)^2 ) rho dz 
    201             ! rho = reference density 
    202             ! S = isoneutral slope. 
    203             ! Most terms are on W grid so work on this grid 
    204 #ifdef key_traldf_eiv 
    205             CALL wrk_alloc( jpi, jpj, zke2d ) 
    206             zke2d(:,:) = 0._wp 
    207             DO jk = 1,jpk 
    208                DO ji = 1,jpi 
    209                   DO jj = 1,jpj 
    210                      zke2d(ji,jj) = zke2d(ji,jj) +  rau0 * fsaeiw(ji, jj, jk)               & 
    211                           &                      * ( wslpi(ji, jj, jk) * wslpi(ji,jj,jk)    & 
    212                           &                      +   wslpj(ji, jj, jk) * wslpj(ji,jj,jk) )  & 
    213                           &                      *   rn2(ji,jj,jk) * fse3w(ji, jj, jk) 
    214                   ENDDO 
    215                ENDDO 
    216             ENDDO 
    217             CALL iom_put("ketrd_eiv", zke2d) 
    218             CALL wrk_dealloc( jpi, jpj, zke2d ) 
    219 #endif 
     197!!gm moved in traadv_eiv ===>>> diag becomes accessible without ln_trdtra=T 
     198!        CASE( jpdyn_eivke ) 
     199!            ! CMIP6 diagnostic tknebto = tendency of EKE from 
     200!            ! parameterized mesoscale eddy advection 
     201!            ! = vertical_integral( k (N S)^2 ) rho dz 
     202!            ! rho = reference density 
     203!            ! S = isoneutral slope. 
     204!            ! Most terms are on W grid so work on this grid 
     205!#ifdef key_traldf_eiv 
     206!            CALL wrk_alloc( jpi, jpj, zke2d ) 
     207!            zke2d(:,:) = 0._wp 
     208!            DO jk = 1,jpk 
     209!               DO ji = 1,jpi 
     210!                  DO jj = 1,jpj 
     211!                     zke2d(ji,jj) = zke2d(ji,jj) +  rau0 * fsaeiw(ji, jj, jk)               & 
     212!                          &                      * ( wslpi(ji, jj, jk) * wslpi(ji,jj,jk)    & 
     213!                          &                      +   wslpj(ji, jj, jk) * wslpj(ji,jj,jk) )  & 
     214!                          &                      *   rn2(ji,jj,jk) * fse3w(ji, jj, jk) 
     215!                  ENDDO 
     216!               ENDDO 
     217!            ENDDO 
     218!            CALL iom_put("ketrd_eiv", zke2d) 
     219!            CALL wrk_dealloc( jpi, jpj, zke2d ) 
     220!#endif 
     221!!gm end 
    220222         ! 
    221223      END SELECT 
Note: See TracChangeset for help on using the changeset viewer.