- Timestamp:
- 2020-07-03T20:46:53+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/ZDF/zdftke.F90
r13247 r13248 98 98 !! * Substitutions 99 99 # include "do_loop_substitute.h90" 100 # include "domzgr_substitute.h90" 100 101 !!---------------------------------------------------------------------- 101 102 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 267 268 zpelc(:,:,1) = MAX( rn2b(:,:,1), 0._wp ) * gdepw(:,:,1,Kmm) * e3w(:,:,1,Kmm) 268 269 DO jk = 2, jpk 269 zpelc(:,:,jk) = zpelc(:,:,jk-1) + MAX( rn2b(:,:,jk), 0._wp ) * gdepw(:,:,jk,Kmm) * e3w(:,:,jk,Kmm) 270 zpelc(:,:,jk) = zpelc(:,:,jk-1) + & 271 & MAX( rn2b(:,:,jk), 0._wp ) * gdepw(:,:,jk,Kmm) * e3w(:,:,jk,Kmm) 270 272 END DO 271 273 ! !* finite Langmuir Circulation depth … … 325 327 ! ! eddy coefficient (ensure numerical stability) 326 328 zzd_up = zcof * MAX( p_avm(ji,jj,jk+1) + p_avm(ji,jj,jk ) , 2.e-5_wp ) & ! upper diagonal 327 & / ( e3t(ji,jj,jk ,Kmm) * e3w(ji,jj,jk ,Kmm) ) 329 & / ( e3t(ji,jj,jk ,Kmm) & 330 & * e3w(ji,jj,jk ,Kmm) ) 328 331 zzd_lw = zcof * MAX( p_avm(ji,jj,jk ) + p_avm(ji,jj,jk-1) , 2.e-5_wp ) & ! lower diagonal 329 & / ( e3t(ji,jj,jk-1,Kmm) * e3w(ji,jj,jk ,Kmm) ) 332 & / ( e3t(ji,jj,jk-1,Kmm) & 333 & * e3w(ji,jj,jk ,Kmm) ) 330 334 ! 331 335 zd_up(ji,jj,jk) = zzd_up ! Matrix (zdiag, zd_up, zd_lw) … … 515 519 & gdepw(ji,jj,mbkt(ji,jj)+1,Kmm) - gdepw(ji,jj,jk,Kmm) ) 516 520 ! wmask prevent zmxlm = 0 if jk = mikt(ji,jj) 517 zmxlm(ji,jj,jk) = zemxl * wmask(ji,jj,jk) + MIN( zmxlm(ji,jj,jk) , e3w(ji,jj,jk,Kmm) ) * (1 - wmask(ji,jj,jk)) 518 zmxld(ji,jj,jk) = zemxl * wmask(ji,jj,jk) + MIN( zmxlm(ji,jj,jk) , e3w(ji,jj,jk,Kmm) ) * (1 - wmask(ji,jj,jk)) 521 zmxlm(ji,jj,jk) = zemxl * wmask(ji,jj,jk) & 522 & + MIN( zmxlm(ji,jj,jk) , e3w(ji,jj,jk,Kmm) ) * (1 - wmask(ji,jj,jk)) 523 zmxld(ji,jj,jk) = zemxl * wmask(ji,jj,jk) & 524 & + MIN( zmxlm(ji,jj,jk) , e3w(ji,jj,jk,Kmm) ) * (1 - wmask(ji,jj,jk)) 519 525 END_3D 520 526 ! … … 528 534 CASE ( 2 ) ! |dk[xml]| bounded by e3t : 529 535 DO_3D_00_00( 2, jpkm1 ) 530 zmxlm(ji,jj,jk) = MIN( zmxlm(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 536 zmxlm(ji,jj,jk) = & 537 & MIN( zmxlm(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 531 538 END_3D 532 539 DO_3DS_00_00( jpkm1, 2, -1 ) … … 538 545 CASE ( 3 ) ! lup and ldown, |dk[xml]| bounded by e3t : 539 546 DO_3D_00_00( 2, jpkm1 ) 540 zmxld(ji,jj,jk) = MIN( zmxld(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 547 zmxld(ji,jj,jk) = & 548 & MIN( zmxld(ji,jj,jk-1) + e3t(ji,jj,jk-1,Kmm), zmxlm(ji,jj,jk) ) 541 549 END_3D 542 550 DO_3DS_00_00( jpkm1, 2, -1 ) 543 zmxlm(ji,jj,jk) = MIN( zmxlm(ji,jj,jk+1) + e3t(ji,jj,jk+1,Kmm), zmxlm(ji,jj,jk) ) 551 zmxlm(ji,jj,jk) = & 552 & MIN( zmxlm(ji,jj,jk+1) + e3t(ji,jj,jk+1,Kmm), zmxlm(ji,jj,jk) ) 544 553 END_3D 545 554 DO_3D_00_00( 2, jpkm1 )
Note: See TracChangeset
for help on using the changeset viewer.