Changeset 15231 for NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diahth.F90
- Timestamp:
- 2021-09-08T09:58:57+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diahth.F90
r12276 r15231 365 365 ! 366 366 ilevel(:,:) = 1 367 DO jk = 2, jpkm1367 DO jk = 1, jpkm1 368 368 DO jj = 1, jpj 369 369 DO ji = 1, jpi 370 IF( ( gdep t_n(ji,jj,jk) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN371 ilevel(ji,jj) = jk372 373 370 IF( ( gdepw_n(ji,jj,jk+1) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN 371 ilevel(ji,jj) = jk+1 372 zthick(ji,jj) = zthick(ji,jj) + e3t_n(ji,jj,jk) 373 phtc (ji,jj) = phtc (ji,jj) + e3t_n(ji,jj,jk) * ptn(ji,jj,jk) 374 374 ENDIF 375 375 ENDDO … … 380 380 DO ji = 1, jpi 381 381 ik = ilevel(ji,jj) 382 zthick(ji,jj) = pdep - zthick(ji,jj) ! remaining thickness to reach depht pdep 383 phtc(ji,jj) = phtc(ji,jj) + ptn(ji,jj,ik+1) * MIN( e3t_n(ji,jj,ik+1), zthick(ji,jj) ) & 384 * tmask(ji,jj,ik+1) 382 IF( tmask(ji,jj,ik) == 1 ) THEN 383 zthick(ji,jj) = MIN ( gdepw_n(ji,jj,ik+1), pdep ) - zthick(ji,jj) ! remaining thickness to reach dephw pdep 384 phtc(ji,jj) = phtc(ji,jj) + ptn(ji,jj,ik) * zthick(ji,jj) 385 ENDIF 385 386 END DO 386 387 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.