- Timestamp:
- 2017-12-26T17:32:56+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r9104 r9169 678 678 WRITE(numout,*) ' type of tke penetration profile nn_htau = ', nn_htau 679 679 WRITE(numout,*) ' fraction of TKE that penetrates rn_efr = ', rn_efr 680 WRITE(numout,*)681 680 IF( ln_drg ) THEN 681 WRITE(numout,*) 682 682 WRITE(numout,*) ' Namelist namdrg_top/_bot: used values:' 683 683 WRITE(numout,*) ' top ocean cavity roughness (m) rn_z0(_top)= ', r_z0_top … … 685 685 ENDIF 686 686 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 689 688 WRITE(numout,*) 690 689 ENDIF … … 693 692 rn_emin = 1.e-10_wp ! specific values of rn_emin & rmxl_min are used 694 693 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' 696 695 ELSE ! standard case : associated avt minimum = molecular viscosity (10^-6 m2/s) 697 696 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_min697 IF(lwp) WRITE(numout,*) ' ==>> minimum mixing length with your parameters rmxl_min = ', rmxl_min 699 698 ENDIF 700 699 ! … … 709 708 ! 710 709 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' 712 712 rn_mxl0 = rmxl_min 713 713 ENDIF … … 763 763 CALL iom_get( numror, jpdom_autoglo, 'dissl', dissl ) 764 764 ELSE ! start TKE from rest 765 IF(lwp) WRITE(numout,*) 765 766 IF(lwp) WRITE(numout,*) ' ==>> previous run without TKE scheme, set en to background values' 766 767 en (:,:,:) = rn_emin * wmask(:,:,:) … … 769 770 ENDIF 770 771 ELSE !* Start from rest 772 IF(lwp) WRITE(numout,*) 771 773 IF(lwp) WRITE(numout,*) ' ==>> start from rest: set en to the background value' 772 774 en (:,:,:) = rn_emin * wmask(:,:,:) … … 777 779 ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file 778 780 ! ! ------------------- 779 IF(lwp) WRITE(numout,*) '---- tke -rst ----'781 IF(lwp) WRITE(numout,*) '---- tke_rst ----' 780 782 CALL iom_rstput( kt, nitrst, numrow, 'en' , en ) 781 783 CALL iom_rstput( kt, nitrst, numrow, 'avt_k', avt_k )
Note: See TracChangeset
for help on using the changeset viewer.