- Timestamp:
- 2020-03-27T20:55:44+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdftke.F90
r12616 r12622 261 261 ! 262 262 ! !* total energy produce by LC : cumulative sum over jk 263 zpelc(:,:,1) = 263 zpelc(:,:,1) = MAX( rn2b(:,:,1), 0._wp ) * gdepw(:,:,1,Kmm) * e3w(:,:,1,Kmm) 264 264 DO jk = 2, jpk 265 zpelc(:,:,jk) = zpelc(:,:,jk-1) + MAX( rn2b(:,:,jk), 0._wp ) * gdepw(:,:,jk,Kmm) * e3w(:,:,jk,Kmm) 265 zpelc(:,:,jk) = zpelc(:,:,jk-1) + & 266 & MAX( rn2b(:,:,jk), 0._wp ) * gdepw(:,:,jk,Kmm) * e3w(:,:,jk,Kmm) 266 267 END DO 267 268 ! !* finite Langmuir Circulation depth … … 317 318 ! ! eddy coefficient (ensure numerical stability) 318 319 zzd_up = zcof * MAX( p_avm(ji,jj,jk+1) + p_avm(ji,jj,jk ) , 2.e-5_wp ) & ! upper diagonal 319 & / ( e3t(ji,jj,jk ,Kmm) * e3w(ji,jj,jk ,Kmm) ) 320 & / ( e3t(ji,jj,jk ,Kmm) & 321 & * e3w(ji,jj,jk ,Kmm) ) 320 322 zzd_lw = zcof * MAX( p_avm(ji,jj,jk ) + p_avm(ji,jj,jk-1) , 2.e-5_wp ) & ! lower diagonal 321 & / ( e3t(ji,jj,jk-1,Kmm) * e3w(ji,jj,jk ,Kmm) ) 323 & / ( e3t(ji,jj,jk-1,Kmm) & 324 & * e3w(ji,jj,jk ,Kmm) ) 322 325 ! 323 326 zd_up(ji,jj,jk) = zzd_up ! Matrix (zdiag, zd_up, zd_lw) … … 483 486 CASE ( 2 ) ! |dk[xml]| bounded by e3t : 484 487 DO_3D_00_00( 2, jpkm1 ) 485 zmxlm(ji,jj,jk) = MIN( zmxlm(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 488 zmxlm(ji,jj,jk) = & 489 & MIN( zmxlm(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 486 490 END_3D 487 491 DO_3DS_00_00( jpkm1, 2, -1 ) … … 493 497 CASE ( 3 ) ! lup and ldown, |dk[xml]| bounded by e3t : 494 498 DO_3D_00_00( 2, jpkm1 ) 495 zmxld(ji,jj,jk) = MIN( zmxld(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 499 zmxld(ji,jj,jk) = & 500 & MIN( zmxld(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 496 501 END_3D 497 502 DO_3DS_00_00( jpkm1, 2, -1 ) 498 zmxlm(ji,jj,jk) = MIN( zmxlm(ji,jj,jk+1) + e3t(ji,jj,jk+1,Kmm), zmxlm(ji,jj,jk) ) 503 zmxlm(ji,jj,jk) = & 504 & MIN( zmxlm(ji,jj,jk+1) + e3t(ji,jj,jk+1,Kmm), zmxlm(ji,jj,jk) ) 499 505 END_3D 500 506 DO_3D_00_00( 2, jpkm1 )
Note: See TracChangeset
for help on using the changeset viewer.