Ignore:
Timestamp:
2018-05-04T07:57:05+02:00 (2 years ago)
Author:
gm
Message:

dev_merge_2017: add a TKE parameter to control the activation of nn_etau below sea-ice

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r9440 r9546  
    7575   INTEGER  ::      nn_htau   ! type of tke profile of penetration (=0/1) 
    7676   REAL(wp) ::      rn_efr    ! fraction of TKE surface value which penetrates in the ocean 
     77   REAL(wp) ::      rn_eice   ! =0 ON below sea-ice, =4 OFF when ice fraction > 1/4    
    7778   LOGICAL  ::   ln_lc     ! Langmuir cells (LC) as a source term of TKE or not 
    7879   REAL(wp) ::      rn_lc     ! coef to compute vertical velocity of Langmuir cells 
     
    298299                  !                                           ! TKE Langmuir circulation source term 
    299300                  en(ji,jj,jk) = en(ji,jj,jk) + rdt * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * ( zwlc * zwlc * zwlc )   & 
    300                      &                              / zhlc(ji,jj) * wmask(ji,jj,jk) 
    301 !!gm                     &                          / zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     301                     &                              / zhlc(ji,jj) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    302302               END DO 
    303303            END DO 
     
    397397       
    398398      IF( nn_etau == 1 ) THEN           !* penetration below the mixed layer (rn_efr fraction) 
    399          DO jk = 2, jpkm1 
     399         DO jk = 2, jpkm1                       ! rn_eice =0 ON below sea-ice, =4 OFF when ice fraction > 0.25 
    400400            DO jj = 2, jpjm1 
    401401               DO ji = fs_2, fs_jpim1   ! vector opt. 
    402402                  en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    403                      &                                 * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) 
    404 !!gm                     &                             * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    405 !!gm    multiplication by surface tmask useless (already includes in MAX( 0, 1-4*fr_i ) 
     403                     &                                 * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    406404               END DO 
    407405            END DO 
     
    412410               jk = nmln(ji,jj) 
    413411               en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    414                   &                                 * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) 
    415 !!gm                  &                             * MAX(0.,1._wp - 4.*fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     412                  &                                 * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) )  * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    416413            END DO 
    417414         END DO 
     
    426423                  zdif = rhftau_scl * MAX( 0._wp, zdif + rhftau_add )  ! apply some modifications... 
    427424                  en(ji,jj,jk) = en(ji,jj,jk) + zbbrau * zdif * EXP( -pdepw(ji,jj,jk) / htau(ji,jj) )   & 
    428                      &                        * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * wmask(ji,jj,jk) 
    429 !!gm                     &                    * MAX(0.,1._wp - 4.*fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
     425                     &                        * MAX(0.,1._wp - rn_eice *fr_i(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) 
    430426               END DO 
    431427            END DO 
     
    683679         WRITE(numout,*) '          type of tke penetration profile            nn_htau   = ', nn_htau 
    684680         WRITE(numout,*) '          fraction of TKE that penetrates            rn_efr    = ', rn_efr 
     681         WRITE(numout,*) '          below sea-ice:  =0 ON                      rn_eice   = ', rn_eice 
     682         WRITE(numout,*) '          =4 OFF when ice fraction > 1/4   ' 
    685683         IF( ln_drg ) THEN 
    686684            WRITE(numout,*) 
Note: See TracChangeset for help on using the changeset viewer.