Changeset 13905 for NEMO/branches/2020/dev_r13508_HPC-09_communications_cleanup/src/OCE/ZDF/zdftke.F90
- Timestamp:
- 2020-11-28T14:29:12+01:00 (3 years ago)
- Location:
- NEMO/branches/2020/dev_r13508_HPC-09_communications_cleanup
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13508_HPC-09_communications_cleanup
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@135 07sette10 ^/utils/CI/sette@13559 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r13508_HPC-09_communications_cleanup/src/OCE/ZDF/zdftke.F90
r13497 r13905 52 52 #endif 53 53 ! 54 #if defined key_si355 USE ice, ONLY: hm_i, h_i56 #endif57 #if defined key_cice58 USE sbc_ice, ONLY: h_i59 #endif60 54 USE in_out_manager ! I/O manager 61 55 USE iom ! I/O manager library … … 382 376 !!gm BUG : in the exp remove the depth of ssh !!! 383 377 !!gm i.e. use gde3w in argument (gdepw(:,:,:,Kmm)) 384 385 378 ! 379 ! penetration is partly switched off below sea-ice if nn_eice/=0 380 ! 386 381 IF( nn_etau == 1 ) THEN !* penetration below the mixed layer (rn_efr fraction) 387 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) ! nn_eice=0 : ON below sea-ice ; nn_eice>0 : partly OFF382 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 388 383 en(ji,jj,jk) = en(ji,jj,jk) + rn_efr * en(ji,jj,1) * EXP( -gdepw(ji,jj,jk,Kmm) / htau(ji,jj) ) & 389 384 & * MAX( 0._wp, 1._wp - zice_fra(ji,jj) ) * wmask(ji,jj,jk) * tmask(ji,jj,1) … … 396 391 END_2D 397 392 ELSEIF( nn_etau == 3 ) THEN !* penetration belox the mixed layer (HF variability) 398 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) ! nn_eice=0 : ON below sea-ice ; nn_eice>0 : partly OFF393 DO_3D( 0, 0, 0, 0, 2, jpkm1 ) 399 394 ztx2 = utau(ji-1,jj ) + utau(ji,jj) 400 395 zty2 = vtau(ji ,jj-1) + vtau(ji,jj) … … 655 650 WRITE(numout,*) ' mixing length type nn_mxl = ', nn_mxl 656 651 WRITE(numout,*) ' surface mixing length = F(stress) or not ln_mxl0 = ', ln_mxl0 657 IF( ln_mxl0 ) THEN658 WRITE(numout,*) ' type of scaling under sea-ice nn_mxlice = ', nn_mxlice659 IF( nn_mxlice == 1 ) &660 WRITE(numout,*) ' ice thickness when scaling under sea-ice rn_mxlice = ', rn_mxlice661 ENDIF662 652 WRITE(numout,*) ' surface mixing length minimum value rn_mxl0 = ', rn_mxl0 663 653 IF( ln_mxl0 ) THEN … … 688 678 CALL ctl_stop( 'zdf_tke_init: wrong value for nn_eice, should be 0,1,2, or 3') 689 679 END SELECT 690 IF( .NOT.ln_drg_OFF ) THEN691 WRITE(numout,*)692 WRITE(numout,*) ' Namelist namdrg_top/_bot: used values:'693 WRITE(numout,*) ' top ocean cavity roughness (m) rn_z0(_top)= ', r_z0_top694 WRITE(numout,*) ' Bottom seafloor roughness (m) rn_z0(_bot)= ', r_z0_bot695 ENDIF696 680 WRITE(numout,*) 697 681 WRITE(numout,*) ' ==>>> critical Richardson nb with your parameters ri_cri = ', ri_cri
Note: See TracChangeset
for help on using the changeset viewer.