Changeset 12196


Ignore:
Timestamp:
2019-12-11T18:13:20+01:00 (7 months ago)
Author:
smasson
Message:

trunk: limit bottom slope to ½*e3/e1 in ldf_slp, see #1180

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/LDF/ldfslp.F90

    r10425 r12196  
    207207               zbu = MIN(  zbu, - z1_slpmax * ABS( zau ) , -7.e+3_wp/e3u_n(ji,jj,jk)* ABS( zau )  ) 
    208208               zbv = MIN(  zbv, - z1_slpmax * ABS( zav ) , -7.e+3_wp/e3v_n(ji,jj,jk)* ABS( zav )  ) 
     209               !                                      ! Fred Dupont: add a correction for bottom partial steps: 
     210               !                                      !              max slope = 1/2 * e3 / e1 
     211               IF (ln_zps .AND. jk==mbku(ji,jj)) & 
     212                  zbu = MIN(  zbu, - z1_slpmax * ABS( zau ) , - 2._wp * e1u(ji,jj) / e3u_n(ji,jj,jk)* ABS( zau )  ) 
     213               IF (ln_zps .AND. jk==mbkv(ji,jj)) & 
     214                  zbv = MIN(  zbv, - z1_slpmax * ABS( zav ) , - 2._wp * e2v(ji,jj) / e3v_n(ji,jj,jk)* ABS( zav )  ) 
    209215               !                                      ! uslp and vslp output in zwz and zww, resp. 
    210216               zfi = MAX( omlmask(ji,jj,jk), omlmask(ji+1,jj,jk) ) 
Note: See TracChangeset for help on using the changeset viewer.