- Timestamp:
- 2021-09-08T16:07:02+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DIA/diahth.F90
r13497 r15234 334 334 SUBROUTINE dia_hth_htc( Kmm, pdep, pt, phtc ) 335 335 ! 336 INTEGER , INTENT(in) :: Kmm ! ocean time level index337 REAL(wp), INTENT(in) :: pdep ! depth over the heat content338 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) ::pt339 REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) ::phtc340 ! 341 INTEGER :: ji, jj, jk, ik342 REAL(wp), DIMENSION(jpi,jpj) :: zthick343 INTEGER , DIMENSION(jpi,jpj) :: ilevel336 INTEGER , INTENT(in) :: Kmm ! ocean time level index 337 REAL(wp), INTENT(in) :: pdep ! depth over the heat content 338 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in) :: pt 339 REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: phtc 340 ! 341 INTEGER :: ji, jj, jk, ik 342 REAL(wp), DIMENSION(jpi,jpj) :: zthick 343 INTEGER , DIMENSION(jpi,jpj) :: ilevel 344 344 345 345 346 346 ! surface boundary condition 347 347 348 IF( .NOT. ln_linssh ) THEN ; zthick(:,:) = 0._wp ; phtc(:,:) = 0._wp348 IF( .NOT. ln_linssh ) THEN ; zthick(:,:) = 0._wp ; phtc(:,:) = 0._wp 349 349 ELSE ; zthick(:,:) = ssh(:,:,Kmm) ; phtc(:,:) = pt(:,:,1) * ssh(:,:,Kmm) * tmask(:,:,1) 350 350 ENDIF 351 351 ! 352 352 ilevel(:,:) = 1 353 DO_3D( 1, 1, 1, 1, 2, jpkm1 )354 IF( ( gdep t(ji,jj,jk,Kmm) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN355 ilevel(ji,jj) = jk 353 DO_3D( 1, 1, 1, 1, 1, jpkm1 ) 354 IF( ( gdepw(ji,jj,jk+1,Kmm) < pdep ) .AND. ( tmask(ji,jj,jk) == 1 ) ) THEN 355 ilevel(ji,jj) = jk+1 356 356 zthick(ji,jj) = zthick(ji,jj) + e3t(ji,jj,jk,Kmm) 357 357 phtc (ji,jj) = phtc (ji,jj) + e3t(ji,jj,jk,Kmm) * pt(ji,jj,jk) … … 361 361 DO_2D( 1, 1, 1, 1 ) 362 362 ik = ilevel(ji,jj) 363 zthick(ji,jj) = pdep - zthick(ji,jj) ! remaining thickness to reach depht pdep364 phtc(ji,jj) = phtc(ji,jj) &365 & + pt (ji,jj,ik+1) * MIN( e3t(ji,jj,ik+1,Kmm), zthick(ji,jj) ) &366 * tmask(ji,jj,ik+1)363 IF( tmask(ji,jj,ik) == 1 ) THEN 364 zthick(ji,jj) = MIN ( gdepw(ji,jj,ik+1,Kmm), pdep ) - zthick(ji,jj) ! remaining thickness to reach dephw pdep 365 phtc(ji,jj) = phtc(ji,jj) + pt(ji,jj,ik) * zthick(ji,jj) 366 ENDIF 367 367 END_2D 368 !369 368 ! 370 369 END SUBROUTINE dia_hth_htc
Note: See TracChangeset
for help on using the changeset viewer.