- Timestamp:
- 2020-07-03T19:15:31+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/ZDF/zdftke.F90
r13015 r13247 310 310 DO_3D_00_00( 2, jpkm1 ) 311 311 ! ! local Richardson number 312 zri = MAX( rn2b(ji,jj,jk), 0._wp ) * p_avm(ji,jj,jk) / ( p_sh2(ji,jj,jk) + rn_bshear ) 312 IF (rn2b(ji,jj,jk) <= 0.0_wp) then 313 zri = 0.0_wp 314 ELSE 315 zri = rn2b(ji,jj,jk) * p_avm(ji,jj,jk) / ( p_sh2(ji,jj,jk) + rn_bshear ) 316 ENDIF 313 317 ! ! inverse of Prandtl number 314 318 apdlr(ji,jj,jk) = MAX( 0.1_wp, ri_cri / MAX( ri_cri , zri ) )
Note: See TracChangeset
for help on using the changeset viewer.