Changeset 11480 for NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90
- Timestamp:
- 2019-08-29T11:23:25+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/ZDF/zdfosm.F90
r10955 r11480 489 489 490 490 zhbl_t(:,:) = hbl(:,:) + (zdhdt(:,:) - ww(ji,jj,ibld(ji,jj)))* rn_rdt ! certainly need wb here, so subtract it 491 zhbl_t(:,:) = MIN(zhbl_t(:,:), ht _n(:,:))491 zhbl_t(:,:) = MIN(zhbl_t(:,:), ht(:,:)) 492 492 zdhdt(:,:) = MIN(zdhdt(:,:), (zhbl_t(:,:) - hbl(:,:))/rn_rdt + ww(ji,jj,ibld(ji,jj))) ! adjustment to represent limiting by ocean bottom 493 493 … … 525 525 526 526 zhbl_s = zhbl_s + MIN( - zwb_ent(ji,jj) / zdb * rn_rdt / FLOAT(ibld(ji,jj)-imld(ji,jj) ), e3w(ji,jj,jk,Kmm) ) 527 zhbl_s = MIN(zhbl_s, ht _n(ji,jj))527 zhbl_s = MIN(zhbl_s, ht(ji,jj)) 528 528 529 529 IF ( zhbl_s >= gdepw(ji,jj,jm+1,Kmm) ) jm = jm + 1 … … 546 546 & * zwstrl(ji,jj)**3 / hbli(ji,jj) ) / zdb * e3w(ji,jj,jk,Kmm) / zdhdt(ji,jj) ! ALMG to investigate whether need to include ww here 547 547 548 zhbl_s = MIN(zhbl_s, ht _n(ji,jj))548 zhbl_s = MIN(zhbl_s, ht(ji,jj)) 549 549 IF ( zhbl_s >= gdepw(ji,jj,jm,Kmm) ) jm = jm + 1 550 550 END DO
Note: See TracChangeset
for help on using the changeset viewer.