Changeset 15643
- Timestamp:
- 2022-01-14T14:13:55+01:00 (3 years ago)
- Location:
- NEMO/branches/UKMO/NEMO_4.0.4_GO6_mixing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/UKMO/NEMO_4.0.4_GO6_mixing/cfgs/SHARED/namelist_ref
r14075 r15643 1076 1076 ! = 0 constant 10 m length scale 1077 1077 ! = 1 0.5m at the equator to 30m poleward of 40 degrees 1078 rn_htau_scaling = 1.0 ! scaling factor to apply to depth of TKE penetration 1078 1079 nn_eice = 1 ! attenutaion of langmuir & surface wave breaking under ice 1079 1080 ! ! = 0 no impact of ice cover on langmuir & surface wave breaking -
NEMO/branches/UKMO/NEMO_4.0.4_GO6_mixing/src/OCE/ZDF/zdftke.F90
r14099 r15643 81 81 INTEGER :: nn_etau ! type of depth penetration of surface tke (=0/1/2/3) 82 82 INTEGER :: nn_htau ! type of tke profile of penetration (=0/1) 83 REAL(wp) :: rn_htau_scaling ! a acaling factor to apply to the penetration of TKE 83 84 REAL(wp) :: rn_efr ! fraction of TKE surface value which penetrates in the ocean 84 85 LOGICAL :: ln_lc ! Langmuir cells (LC) as a source term of TKE or not … … 720 721 & rn_mxl0 , nn_mxlice, rn_mxlice, & 721 722 & nn_pdl , ln_lc , rn_lc, & 723 & rn_htau_scaling , & 722 724 & nn_etau , nn_htau , rn_efr , nn_eice 723 725 !!---------------------------------------------------------------------- … … 766 768 WRITE(numout,*) ' test param. to add tke induced by wind nn_etau = ', nn_etau 767 769 WRITE(numout,*) ' type of tke penetration profile nn_htau = ', nn_htau 770 WRITE(numout,*) ' scaling factor for tke penetration depth rn_htau_scaling = ', rn_htau_scaling 768 771 WRITE(numout,*) ' fraction of TKE that penetrates rn_efr = ', rn_efr 769 772 WRITE(numout,*) ' langmuir & surface wave breaking under ice nn_eice = ', nn_eice … … 817 820 SELECT CASE( nn_htau ) ! Choice of the depth of penetration 818 821 CASE( 0 ) ! constant depth penetration (here 10 meters) 819 htau(:,:) = 10._wp822 htau(:,:) = rn_htau_scaling*10._wp 820 823 CASE( 1 ) ! F(latitude) : 0.5m to 30m poleward of 40 degrees 821 htau(:,:) = MAX( 0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) ) )824 htau(:,:) = MAX( 0.5_wp, rn_htau_scaling*MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) ) ) 822 825 CASE( 4 ) ! F(latitude) : 0.5m to 10m/30m poleward of 13/40 degrees north/south 823 826 DO jj = 2, jpjm1 824 827 DO ji = fs_2, fs_jpim1 ! vector opt. 825 828 IF( gphit(ji,jj) <= 0._wp ) THEN 826 htau(ji,jj) = MAX( 0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) )829 htau(ji,jj) = MAX( 0.5_wp, rn_htau_scaling*MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) ) 827 830 ELSE 828 htau(ji,jj) = MAX( 0.5_wp, MIN( 10._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) )831 htau(ji,jj) = MAX( 0.5_wp, rn_htau_scaling*MIN( 10._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) ) ) 829 832 ENDIF 830 833 END DO
Note: See TracChangeset
for help on using the changeset viewer.