Changeset 12279 for NEMO/trunk/src/OCE/TRA
- Timestamp:
- 2019-12-20T14:36:15+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/TRA/tranxt.F90
r10425 r12279 267 267 INTEGER :: ji, jj, jk, jn ! dummy loop indices 268 268 REAL(wp) :: zfact, zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d ! local scalar 269 REAL(wp) :: zfact2, ze3t_b, ze3t_n, ze3t_a, ze3t_f, ze3t_d 269 REAL(wp) :: zfact2, ze3t_b, ze3t_n, ze3t_a, ze3t_f, ze3t_d, zscale ! - - 270 270 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) :: ztrd_atf 271 271 !!---------------------------------------------------------------------- … … 312 312 ztc_f = ztc_n + atfp * ztc_d 313 313 ! 314 zscale = zfact2 * e3t_n(ji,jj,jk) / ( ht_n(ji,jj) + 1._wp - ssmask(ji,jj) ) 315 ze3t_f = ze3t_f - zscale * ( emp_b(ji,jj) - emp(ji,jj) ) 316 IF ( ll_rnf ) ze3t_f = ze3t_f - zscale * ( rnf_b(ji,jj) - rnf(ji,jj) ) 317 IF ( ll_isf ) ze3t_f = ze3t_f - zscale * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) ) 318 314 319 IF( jk == mikt(ji,jj) ) THEN ! first level 315 ze3t_f = ze3t_f - zfact2 * ( (emp_b(ji,jj) - emp(ji,jj) ) &316 & + (fwfisf_b(ji,jj) - fwfisf(ji,jj)) )317 320 ztc_f = ztc_f - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) ) 318 321 ENDIF 319 IF( ln_rnf_depth ) THEN320 ! Rivers are not just at the surface must go down to nk_rnf(ji,jj)321 IF( mikt(ji,jj) <=jk .and. jk <= nk_rnf(ji,jj) ) THEN322 ze3t_f = ze3t_f - zfact2 * ( - (rnf_b(ji,jj) - rnf(ji,jj) ) ) &323 & * ( e3t_n(ji,jj,jk) / h_rnf(ji,jj) )324 ENDIF325 ELSE326 IF( jk == mikt(ji,jj) ) THEN ! first level327 ze3t_f = ze3t_f - zfact2 * ( - (rnf_b(ji,jj) - rnf(ji,jj) ) )328 ENDIF329 ENDIF330 331 322 ! 332 323 ! solar penetration (temperature only)
Note: See TracChangeset
for help on using the changeset viewer.