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 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90 – NEMO

Ignore:
Timestamp:
2017-12-26T17:32:56+01:00 (6 years ago)
Author:
gm
Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

File:
1 edited

Legend:

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

    r9104 r9169  
    678678         WRITE(numout,*) '          type of tke penetration profile            nn_htau   = ', nn_htau 
    679679         WRITE(numout,*) '          fraction of TKE that penetrates            rn_efr    = ', rn_efr 
    680          WRITE(numout,*) 
    681680         IF( ln_drg ) THEN 
     681            WRITE(numout,*) 
    682682            WRITE(numout,*) '   Namelist namdrg_top/_bot:   used values:' 
    683683            WRITE(numout,*) '      top    ocean cavity roughness (m)          rn_z0(_top)= ', r_z0_top 
     
    685685         ENDIF 
    686686         WRITE(numout,*) 
    687          WRITE(numout,*) 
    688          WRITE(numout,*) '   ==>> critical Richardson nb with your parameters  ri_cri = ', ri_cri 
     687         WRITE(numout,*) '   ==>>   critical Richardson nb with your parameters  ri_cri = ', ri_cri 
    689688         WRITE(numout,*) 
    690689      ENDIF 
     
    693692         rn_emin  = 1.e-10_wp             ! specific values of rn_emin & rmxl_min are used 
    694693         rmxl_min = 1.e-03_wp             ! associated avt minimum = molecular salt diffusivity (10^-9 m2/s) 
    695          IF(lwp) WRITE(numout,*) '      Internal wave-driven mixing case:   force   rn_emin = 1.e-10 and rmxl_min = 1.e-3 ' 
     694         IF(lwp) WRITE(numout,*) '   ==>>   Internal wave-driven mixing case:   force   rn_emin = 1.e-10 and rmxl_min = 1.e-3' 
    696695      ELSE                          ! standard case : associated avt minimum = molecular viscosity (10^-6 m2/s) 
    697696         rmxl_min = 1.e-6_wp / ( rn_ediff * SQRT( rn_emin ) )    ! resulting minimum length to recover molecular viscosity 
    698          IF(lwp) WRITE(numout,*) '      minimum mixing length with your parameters rmxl_min = ', rmxl_min 
     697         IF(lwp) WRITE(numout,*) '   ==>>   minimum mixing length with your parameters rmxl_min = ', rmxl_min 
    699698      ENDIF 
    700699      ! 
     
    709708      ! 
    710709      IF( ln_mxl0 ) THEN 
    711          IF(lwp) WRITE(numout,*) '   use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min' 
     710         IF(lwp) WRITE(numout,*) 
     711         IF(lwp) WRITE(numout,*) '   ==>>   use a surface mixing length = F(stress) :   set rn_mxl0 = rmxl_min' 
    712712         rn_mxl0 = rmxl_min 
    713713      ENDIF 
     
    763763               CALL iom_get( numror, jpdom_autoglo, 'dissl', dissl ) 
    764764            ELSE                                          ! start TKE from rest 
     765               IF(lwp) WRITE(numout,*) 
    765766               IF(lwp) WRITE(numout,*) '   ==>>   previous run without TKE scheme, set en to background values' 
    766767               en   (:,:,:) = rn_emin * wmask(:,:,:) 
     
    769770            ENDIF 
    770771         ELSE                                   !* Start from rest 
     772            IF(lwp) WRITE(numout,*) 
    771773            IF(lwp) WRITE(numout,*) '   ==>>   start from rest: set en to the background value' 
    772774            en   (:,:,:) = rn_emin * wmask(:,:,:) 
     
    777779      ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN   ! Create restart file 
    778780         !                                   ! ------------------- 
    779          IF(lwp) WRITE(numout,*) '---- tke-rst ----' 
     781         IF(lwp) WRITE(numout,*) '---- tke_rst ----' 
    780782         CALL iom_rstput( kt, nitrst, numrow, 'en'   , en    ) 
    781783         CALL iom_rstput( kt, nitrst, numrow, 'avt_k', avt_k ) 
Note: See TracChangeset for help on using the changeset viewer.