- Timestamp:
- 2019-05-16T15:23:56+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0_mirror_text_diagnostics/src/OCE/ZDF/zdftke.F90
r10968 r10986 699 699 WRITE(numout,*) ' ==>>> critical Richardson nb with your parameters ri_cri = ', ri_cri 700 700 WRITE(numout,*) 701 IF(lflush) CALL FLUSH(numout) 701 702 ENDIF 702 703 ! … … 704 705 rn_emin = 1.e-10_wp ! specific values of rn_emin & rmxl_min are used 705 706 rmxl_min = 1.e-03_wp ! associated avt minimum = molecular salt diffusivity (10^-9 m2/s) 706 IF(lwp) WRITE(numout,*) ' ==>>> Internal wave-driven mixing case: force rn_emin = 1.e-10 and rmxl_min = 1.e-3' 707 IF(lwp) THEN 708 WRITE(numout,*) ' ==>>> Internal wave-driven mixing case: force rn_emin = 1.e-10 and rmxl_min = 1.e-3' 709 IF(lflush) CALL FLUSH(numout) 710 ENDIF 707 711 ELSE ! standard case : associated avt minimum = molecular viscosity (10^-6 m2/s) 708 712 rmxl_min = 1.e-6_wp / ( rn_ediff * SQRT( rn_emin ) ) ! resulting minimum length to recover molecular viscosity 709 IF(lwp) WRITE(numout,*) ' ==>>> minimum mixing length with your parameters rmxl_min = ', rmxl_min 713 IF(lwp) THEN 714 WRITE(numout,*) ' ==>>> minimum mixing length with your parameters rmxl_min = ', rmxl_min 715 IF(lflush) CALL FLUSH(numout) 716 ENDIF 710 717 ENDIF 711 718 ! … … 720 727 ! 721 728 IF( ln_mxl0 ) THEN 722 IF(lwp) WRITE(numout,*) 723 IF(lwp) WRITE(numout,*) ' ==>>> use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min' 729 IF(lwp) THEN 730 WRITE(numout,*) 731 WRITE(numout,*) ' ==>>> use a surface mixing length = F(stress) : set rn_mxl0 = rmxl_min' 732 IF(lflush) CALL FLUSH(numout) 733 ENDIF 724 734 rn_mxl0 = rmxl_min 725 735 ENDIF … … 781 791 CALL iom_get( numror, jpdom_autoglo, 'dissl', dissl, ldxios = lrxios ) 782 792 ELSE ! start TKE from rest 783 IF(lwp) WRITE(numout,*) 784 IF(lwp) WRITE(numout,*) ' ==>>> previous run without TKE scheme, set en to background values' 793 IF(lwp) THEN 794 WRITE(numout,*) 795 WRITE(numout,*) ' ==>>> previous run without TKE scheme, set en to background values' 796 IF(lflush) CALL FLUSH(numout) 797 ENDIF 785 798 en (:,:,:) = rn_emin * wmask(:,:,:) 786 799 dissl(:,:,:) = 1.e-12_wp … … 788 801 ENDIF 789 802 ELSE !* Start from rest 790 IF(lwp) WRITE(numout,*) 791 IF(lwp) WRITE(numout,*) ' ==>>> start from rest: set en to the background value' 803 IF(lwp) THEN 804 WRITE(numout,*) 805 WRITE(numout,*) ' ==>>> start from rest: set en to the background value' 806 IF(lflush) CALL FLUSH(numout) 807 ENDIF 792 808 en (:,:,:) = rn_emin * wmask(:,:,:) 793 809 dissl(:,:,:) = 1.e-12_wp … … 797 813 ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN ! Create restart file 798 814 ! ! ------------------- 799 IF(lwp .AND. nprint > 0) WRITE(numout,*) '---- tke_rst ----' 815 IF(lwp .AND. nprint > 0) THEN 816 WRITE(numout,*) '---- tke_rst ----' 817 IF(lflush) CALL FLUSH(numout) 818 ENDIF 800 819 IF( lwxios ) CALL iom_swap( cwxios_context ) 801 820 CALL iom_rstput( kt, nitrst, numrow, 'en' , en , ldxios = lwxios )
Note: See TracChangeset
for help on using the changeset viewer.