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 3798 – NEMO

Changeset 3798


Ignore:
Timestamp:
2013-02-11T15:16:58+01:00 (11 years ago)
Author:
cbricaud
Message:

add same modifications than TKE to have any impact of evd,tmx,ddm on the turbulent closure solution ; see ticket 954

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_MERGE_2012/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfgls.F90

    r3680 r3798  
    4444   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   mxln    !: now mixing length 
    4545   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   zwall   !: wall function 
     46   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   avt_k   ! not enhanced Kz 
     47   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   avm_k   ! not enhanced Kz 
     48   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   avmu_k  ! not enhanced Kz 
     49   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   avmv_k  ! not enhanced Kz 
    4650   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   ustars2 !: Squared surface velocity scale at T-points 
    4751   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   ustarb2 !: Squared bottom  velocity scale at T-points 
     
    119123      !!---------------------------------------------------------------------- 
    120124      ALLOCATE( en(jpi,jpj,jpk),  mxln(jpi,jpj,jpk), zwall(jpi,jpj,jpk) ,     & 
     125         &      avt_k (jpi,jpj,jpk) , avm_k (jpi,jpj,jpk),                    & 
     126         &      avmu_k(jpi,jpj,jpk) , avmv_k(jpi,jpj,jpk),                    & 
    121127         &      ustars2(jpi,jpj), ustarb2(jpi,jpj)                      , STAT= zdf_gls_alloc ) 
    122128         ! 
     
    159165 
    160166      ustars2 = 0._wp   ;   ustarb2 = 0._wp   ;   psi  = 0._wp   ;   zwall_psi = 0._wp 
     167 
     168      IF( kt /= nit000 ) THEN   ! restore before value to compute tke 
     169         avt (:,:,:) = avt_k (:,:,:) 
     170         avm (:,:,:) = avm_k (:,:,:) 
     171         avmu(:,:,:) = avmu_k(:,:,:) 
     172         avmv(:,:,:) = avmv_k(:,:,:)  
     173      ENDIF 
    161174 
    162175      ! Compute surface and bottom friction at T-points 
     
    882895      ENDIF 
    883896      ! 
     897      avt_k (:,:,:) = avt (:,:,:) 
     898      avm_k (:,:,:) = avm (:,:,:) 
     899      avmu_k(:,:,:) = avmu(:,:,:) 
     900      avmv_k(:,:,:) = avmv(:,:,:) 
     901      ! 
    884902      CALL wrk_dealloc( jpi,jpj, zdep, zflxs, zhsro ) 
    885903      CALL wrk_dealloc( jpi,jpj,jpk, eb, mxlb, shear, eps, zwall_psi, z_elem_a, z_elem_b, z_elem_c, psi ) 
     
    12441262         !                                   ! ------------------- 
    12451263         IF(lwp) WRITE(numout,*) '---- gls-rst ----' 
    1246          CALL iom_rstput( kt, nitrst, numrow, 'en'   , en    ) 
    1247          CALL iom_rstput( kt, nitrst, numrow, 'avt'  , avt   ) 
    1248          CALL iom_rstput( kt, nitrst, numrow, 'avm'  , avm   ) 
    1249          CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu  ) 
    1250          CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv  ) 
    1251          CALL iom_rstput( kt, nitrst, numrow, 'mxln' , mxln  ) 
     1264         CALL iom_rstput( kt, nitrst, numrow, 'en'   , en     ) 
     1265         CALL iom_rstput( kt, nitrst, numrow, 'avt'  , avt_k  ) 
     1266         CALL iom_rstput( kt, nitrst, numrow, 'avm'  , avm_k  ) 
     1267         CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu_k ) 
     1268         CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv_k ) 
     1269         CALL iom_rstput( kt, nitrst, numrow, 'mxln' , mxln   ) 
    12521270         ! 
    12531271      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.