New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 6457 for branches/UKMO/dev_r5518_GC3p0_package/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90 – NEMO

Ignore:
Timestamp:
2016-04-08T11:26:47+02:00 (9 years ago)
Author:
dancopsey
Message:

Merged in changeset 6238 of dev_r5021_nn_etau_revision

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GC3p0_package/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r6450 r6457  
    9393   REAL(wp) ::   rhftau_scl = 1.0_wp       ! scale factor applied to HF part of taum  (nn_etau=3) 
    9494 
     95   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   en             !: now turbulent kinetic energy   [m2/s2] 
     96   REAL(wp)        , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   htau           ! depth of tke penetration (nn_htau) 
    9597   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   e_niw          !: TKE budget- near-inertial waves term 
    96    REAL(wp)        , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   htau           ! depth of tke penetration (nn_htau) 
    9798   REAL(wp)        , ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   efr            ! surface boundary condition for nn_etau = 4 
    9899   REAL(wp)        , ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   dissl          ! now mixing lenght of dissipation 
     
    118119      !!---------------------------------------------------------------------- 
    119120      ALLOCATE(                                                                    & 
     121         &      efr  (jpi,jpj)     , e_niw(jpi,jpj,jpk) ,                         &       
    120122#if defined key_c1d 
    121123         &      e_dis(jpi,jpj,jpk) , e_mix(jpi,jpj,jpk) ,                          & 
    122124         &      e_pdl(jpi,jpj,jpk) , e_ric(jpi,jpj,jpk) ,                          & 
    123125#endif 
    124          &      htau  (jpi,jpj)    , dissl(jpi,jpj,jpk) ,                          & 
    125          &      efr  (jpi,jpj)     , e_niw(jpi,jpj,jpk) ,                          & 
    126          &      STAT= zdf_tke_alloc      ) 
     126         &      en    (jpi,jpj,jpk) , htau  (jpi,jpj)    , dissl(jpi,jpj,jpk) ,     &  
     127         &      avt_k (jpi,jpj,jpk) , avm_k (jpi,jpj,jpk),                          & 
     128         &      avmu_k(jpi,jpj,jpk) , avmv_k(jpi,jpj,jpk), STAT= zdf_tke_alloc      ) 
    127129         ! 
    128130      IF( lk_mpp             )   CALL mpp_sum ( zdf_tke_alloc ) 
     
    842844      ENDIF 
    843845       
    844       IF( nn_etau == 2  .OR. ( nn_etau /= 0 .AND. nn_htau == 2 ) )   CALL zdf_mxl( nit000 - 1 )      ! Initialization of nmln and hmlp 
     846      IF( nn_etau == 2  )   CALL zdf_mxl( nit000 )      ! Initialization of nmln  
     847 
     848      IF( nn_etau /= 0 .and. nn_htau == 2 ) THEN 
     849          ierr = zdf_mxl_alloc() 
     850          nmln(:,:) = nlb10           ! Initialization of nmln 
     851      ENDIF 
    845852 
    846853      !                               !* depth of penetration of surface tke 
Note: See TracChangeset for help on using the changeset viewer.