Changeset 1601 for trunk/NEMO/OPA_SRC/ZDF/zdftke.F90
- Timestamp:
- 2009-08-11T12:09:19+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/ZDF/zdftke.F90
r1537 r1601 63 63 #endif 64 64 65 ! !!! ** Namelist nam _tke **65 ! !!! ** Namelist namzdf_tke ** 66 66 LOGICAL :: ln_mxl0 = .FALSE. ! mixing length scale surface value as function of wind stress or not 67 67 INTEGER :: nn_mxl = 2 ! type of mixing length (=0/1/2/3) … … 616 616 !! viscosity when using a tke turbulent closure scheme 617 617 !! 618 !! ** Method : Read the nam _tke namelist and check the parameters618 !! ** Method : Read the namzdf_tke namelist and check the parameters 619 619 !! called at the first timestep (nit000) 620 620 !! 621 !! ** input : Namlist nam _tke621 !! ** input : Namlist namzdf_tke 622 622 !! 623 623 !! ** Action : Increase by 1 the nstop flag is setting problem encounter … … 625 625 INTEGER :: ji, jj, jk ! dummy loop indices 626 626 !! 627 NAMELIST/nam _tke/ rn_ediff, rn_ediss , rn_ebb, rn_emin, &628 & rn_emin0, rn_bshear, nn_mxl, ln_mxl0, &629 & rn_lmin , rn_lmin0 , nn_pdl, nn_etau, &630 & nn_htau , rn_efr , ln_lc , rn_lc631 !!---------------------------------------------------------------------- 632 633 REWIND ( numnam ) !* Read Namelist nam _tke : Turbulente Kinetic Energy634 READ ( numnam, nam _tke )627 NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb, rn_emin, & 628 & rn_emin0, rn_bshear, nn_mxl, ln_mxl0, & 629 & rn_lmin , rn_lmin0 , nn_pdl, nn_etau, & 630 & nn_htau , rn_efr , ln_lc , rn_lc 631 !!---------------------------------------------------------------------- 632 633 REWIND ( numnam ) !* Read Namelist namzdf_tke : Turbulente Kinetic Energy 634 READ ( numnam, namzdf_tke ) 635 635 636 636 ri_cri = 2. / ( 2. + rn_ediss / rn_ediff ) ! resulting critical Richardson number … … 640 640 WRITE(numout,*) 'zdf_tke : tke turbulent closure scheme - initialisation' 641 641 WRITE(numout,*) '~~~~~~~~' 642 WRITE(numout,*) ' Namelist nam_tke : set tke mixing parameters'643 WRITE(numout,*) ' coef. to compute avtrn_ediff = ', rn_ediff644 WRITE(numout,*) ' Kolmogoroff dissipation coef.rn_ediss = ', rn_ediss645 WRITE(numout,*) ' tke surface input coef.rn_ebb = ', rn_ebb646 WRITE(numout,*) ' minimum value of tkern_emin = ', rn_emin647 WRITE(numout,*) ' surface minimum value of tkern_emin0 = ', rn_emin0648 WRITE(numout,*) ' background shear (>0)rn_bshear= ', rn_bshear649 WRITE(numout,*) ' mixing length typenn_mxl = ', nn_mxl650 WRITE(numout,*) ' prandl number flagnn_pdl = ', nn_pdl651 WRITE(numout,*) ' surface mixing length = F(stress) or notln_mxl0 = ', ln_mxl0652 WRITE(numout,*) ' surface mixing length minimum valuern_lmin0 = ', rn_lmin0653 WRITE(numout,*) ' interior mixing length minimum valuern_lmin0 = ', rn_lmin0654 WRITE(numout,*) ' test param. to add tke induced by windnn_etau = ', nn_etau655 WRITE(numout,*) ' flag for computation of exp. tke profilenn_htau = ', nn_htau656 WRITE(numout,*) ' % of rn_emin0 which pene. the thermoclinern_efr = ', rn_efr657 WRITE(numout,*) ' flag to take into acc. Langmuir circ.ln_lc = ', ln_lc658 WRITE(numout,*) ' coef to compute verticla velocity of LCrn_lc = ', rn_lc642 WRITE(numout,*) ' Namelist namzdf_tke : set tke mixing parameters' 643 WRITE(numout,*) ' coef. to compute avt rn_ediff = ', rn_ediff 644 WRITE(numout,*) ' Kolmogoroff dissipation coef. rn_ediss = ', rn_ediss 645 WRITE(numout,*) ' tke surface input coef. rn_ebb = ', rn_ebb 646 WRITE(numout,*) ' minimum value of tke rn_emin = ', rn_emin 647 WRITE(numout,*) ' surface minimum value of tke rn_emin0 = ', rn_emin0 648 WRITE(numout,*) ' background shear (>0) rn_bshear= ', rn_bshear 649 WRITE(numout,*) ' mixing length type nn_mxl = ', nn_mxl 650 WRITE(numout,*) ' prandl number flag nn_pdl = ', nn_pdl 651 WRITE(numout,*) ' surface mixing length = F(stress) or not ln_mxl0 = ', ln_mxl0 652 WRITE(numout,*) ' surface mixing length minimum value rn_lmin0 = ', rn_lmin0 653 WRITE(numout,*) ' interior mixing length minimum value rn_lmin0 = ', rn_lmin0 654 WRITE(numout,*) ' test param. to add tke induced by wind nn_etau = ', nn_etau 655 WRITE(numout,*) ' flag for computation of exp. tke profile nn_htau = ', nn_htau 656 WRITE(numout,*) ' % of rn_emin0 which pene. the thermocline rn_efr = ', rn_efr 657 WRITE(numout,*) ' flag to take into acc. Langmuir circ. ln_lc = ', ln_lc 658 WRITE(numout,*) ' coef to compute verticla velocity of LC rn_lc = ', rn_lc 659 659 WRITE(numout,*) 660 WRITE(numout,*) ' critical Richardson nb with your choice of coefs.= ', ri_cri660 WRITE(numout,*) ' critical Richardson nb with your parameters ri_cri = ', ri_cri 661 661 ENDIF 662 662 … … 671 671 ! !* depth of penetration of surface tke 672 672 IF( nn_etau /= 0 ) THEN 673 SELECT CASE( nn_htau ) ! Choice of the depth of penetration673 SELECT CASE( nn_htau ) ! Choice of the depth of penetration 674 674 CASE( 0 ) ! constant depth penetration (here 10 meters) 675 675 htau(:,:) = 10.e0 … … 760 760 ! ! ------------------- 761 761 IF(lwp) WRITE(numout,*) '---- tke-rst ----' 762 CALL iom_rstput( kt, nitrst, numrow, 'en' , en 763 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt 764 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm 765 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu 766 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv 767 CALL iom_rstput( kt, nitrst, numrow, 'dissl', dissl 762 CALL iom_rstput( kt, nitrst, numrow, 'en' , en ) 763 CALL iom_rstput( kt, nitrst, numrow, 'avt' , avt ) 764 CALL iom_rstput( kt, nitrst, numrow, 'avm' , avm ) 765 CALL iom_rstput( kt, nitrst, numrow, 'avmu' , avmu ) 766 CALL iom_rstput( kt, nitrst, numrow, 'avmv' , avmv ) 767 CALL iom_rstput( kt, nitrst, numrow, 'dissl', dissl ) 768 768 ! 769 769 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.