Changeset 7295


Ignore:
Timestamp:
2016-11-21T17:14:56+01:00 (5 years ago)
Author:
till
Message:

reverting the changes made before (zmax_scale and zpdlr, in zdftke.F90)
sorry

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_nn_etau_revision/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

    r7294 r7295  
    714714                  !                                          ! local Richardson number 
    715715                  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 )  ) 
    718717!!gm and even better with the use of the "true" ri_crit=0.22222...  (this change the results!) 
    719718!!gm              zpdlr = MAX(  0.1_wp,  ri_crit / MAX( ri_crit , zri )  ) 
     
    832831            htau(:,:) = MAX(  0.5_wp, MIN( 15._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) )   ) 
    833832         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._wp 
     833            zmax_scale = 1._wp                     ! Amount to scale the maximum value by. 
    835834            zmin = 0.5_wp                          ! The minimum value 
    836835            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.