Changeset 7295
- Timestamp:
- 2016-11-21T17:14:56+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_nn_etau_revision/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90
r7294 r7295 714 714 ! ! local Richardson number 715 715 zri = MAX( rn2b(ji,jj,jk), 0._wp ) * zcoef / (zdku + zdkv + rn_bshear ) 716 ! below a mod suggested by Colin Jones (implementation by TK) 717 zpdlr = MAX( 0.05_wp, 1.0 / MAX( 1.0 , ( 1.0 + 10.0*zri ) ) ) 716 zpdlr = MAX( 0.1_wp, 0.2 / MAX( 0.2 , zri ) ) 718 717 !!gm and even better with the use of the "true" ri_crit=0.22222... (this change the results!) 719 718 !!gm zpdlr = MAX( 0.1_wp, ri_crit / MAX( ri_crit , zri ) ) … … 832 831 htau(:,:) = MAX( 0.5_wp, MIN( 15._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) ) ) 833 832 CASE( 4 ) ! F(latitude) : single sine curve per hemisphere 834 zmax_scale = 2._wp/3._wp ! Amount to scale the maximum value by. Originally 1._wp833 zmax_scale = 1._wp ! Amount to scale the maximum value by. 835 834 zmin = 0.5_wp ! The minimum value 836 835 zmax_h_n = 10._wp * zmax_scale ! The hard maximum value (cutoff imposed on sine curve)
Note: See TracChangeset
for help on using the changeset viewer.