- Timestamp:
- 2020-03-27T20:55:44+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdfosm.F90
r12616 r12622 596 596 zwb_ent(ji,jj) = 0._wp 597 597 ENDIF 598 inhml = & 599 & MAX( INT( zari * zhbl(ji,jj) / e3t(ji,jj,ibld(ji,jj),Kmm) ), & 600 & 1 ) 598 inhml = MAX( INT( zari * zhbl(ji,jj) & 599 & / e3t(ji,jj,ibld(ji,jj),Kmm) ), 1 ) 601 600 imld(ji,jj) = MAX( ibld(ji,jj) - inhml, 1) 602 601 zhml(ji,jj) = gdepw(ji,jj,imld(ji,jj),Kmm) … … 612 611 IF ( zdb_bl(ji,jj) > 0._wp ) THEN 613 612 ! pycnocline thickness set by stratification - use same relationship as for neutral conditions. 614 zari = MIN( 4.5 * ( zvstr(ji,jj)**2 ) & 615 & / ( zdb_bl(ji,jj) * zhbl(ji,jj) ) + 0.01 , 0.2 ) 616 inhml = & 617 & MAX( INT( zari * zhbl(ji,jj) / e3t(ji,jj,ibld(ji,jj),Kmm) ), & 618 & 1 ) 613 zari = MIN( 4.5 * ( zvstr(ji,jj)**2 ) & 614 & / ( zdb_bl(ji,jj) * zhbl(ji,jj) ) + 0.01 , 0.2 ) 615 inhml = MAX( INT( zari * zhbl(ji,jj) & 616 & / e3t(ji,jj,ibld(ji,jj),Kmm) ), 1 ) 619 617 imld(ji,jj) = MAX( ibld(ji,jj) - inhml, 1) 620 618 zhml(ji,jj) = gdepw(ji,jj,imld(ji,jj),Kmm)
Note: See TracChangeset
for help on using the changeset viewer.