- Timestamp:
- 2020-03-26T11:15:02+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traldf_iso.F90
r12590 r12606 254 254 ! 255 255 DO_2D_00_00 256 pt_rhs(ji,jj,jk,jn) = pt_rhs(ji,jj,jk,jn) + zsign * ( zftu(ji,jj,jk) - zftu(ji-1,jj,jk)&257 & 258 & 256 pt_rhs(ji,jj,jk,jn) = pt_rhs(ji,jj,jk,jn) & 257 & + zsign * ( zftu(ji,jj,jk) - zftu(ji-1,jj,jk) + zftv(ji,jj,jk) - zftv(ji,jj-1,jk) ) & 258 & * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 259 259 END_2D 260 260 END DO ! End of slab … … 301 301 CASE( 1 ) ! 1st pass : eddy coef = ah_wslp2 302 302 DO_3D_00_00( 2, jpkm1 ) 303 ztfw(ji,jj,jk) = ztfw(ji,jj,jk)&304 & + ah_wslp2(ji,jj,jk)* e1e2t(ji,jj) &305 & 303 ztfw(ji,jj,jk) = & 304 & ztfw(ji,jj,jk) + ah_wslp2(ji,jj,jk) * e1e2t(ji,jj) & 305 & * ( pt(ji,jj,jk-1,jn) - pt(ji,jj,jk,jn) ) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk) 306 306 END_3D 307 307 CASE( 2 ) ! 2nd pass : eddy flux = ah_wslp2 and akz applied on pt and pt2 gradients, resp. 308 308 DO_3D_00_00( 2, jpkm1 ) 309 ztfw(ji,jj,jk) = ztfw(ji,jj,jk) + e1e2t(ji,jj) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk) 309 ztfw(ji,jj,jk) = ztfw(ji,jj,jk) + e1e2t(ji,jj) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk) & 310 310 & * ( ah_wslp2(ji,jj,jk) * ( pt (ji,jj,jk-1,jn) - pt (ji,jj,jk,jn) ) & 311 311 & + akz(ji,jj,jk) * ( pt2(ji,jj,jk-1,jn) - pt2(ji,jj,jk,jn) ) ) … … 315 315 ! 316 316 DO_3D_00_00( 1, jpkm1 ) 317 pt_rhs(ji,jj,jk,jn) = pt_rhs(ji,jj,jk,jn) + zsign * ( ztfw (ji,jj,jk) - ztfw(ji,jj,jk+1) ) &318 & * r1_e1e2t(ji,jj)/ e3t(ji,jj,jk,Kmm)317 pt_rhs(ji,jj,jk,jn) = pt_rhs(ji,jj,jk,jn) + zsign * ( ztfw (ji,jj,jk) - ztfw(ji,jj,jk+1) ) * r1_e1e2t(ji,jj) & 318 & / e3t(ji,jj,jk,Kmm) 319 319 END_3D 320 320 !
Note: See TracChangeset
for help on using the changeset viewer.