New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 12622 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdftke.F90 – NEMO

Ignore:
Timestamp:
2020-03-27T20:55:44+01:00 (4 years ago)
Author:
techene
Message:

all: add e3 substitute (sometimes it requires to add ze3t/u/v/w) and limit precompiled files lines to about 130 character

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/ZDF/zdftke.F90

    r12616 r12622  
    261261         ! 
    262262         !                        !* total energy produce by LC : cumulative sum over jk 
    263          zpelc(:,:,1) =  MAX( rn2b(:,:,1), 0._wp ) * gdepw(:,:,1,Kmm) * e3w(:,:,1,Kmm) 
     263         zpelc(:,:,1) = MAX( rn2b(:,:,1), 0._wp ) * gdepw(:,:,1,Kmm) * e3w(:,:,1,Kmm) 
    264264         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) 
    266267         END DO 
    267268         !                        !* finite Langmuir Circulation depth 
     
    317318         !                                   ! eddy coefficient (ensure numerical stability) 
    318319         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)  ) 
    320322         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)  ) 
    322325         ! 
    323326         zd_up(ji,jj,jk) = zzd_up            ! Matrix (zdiag, zd_up, zd_lw) 
     
    483486      CASE ( 2 )           ! |dk[xml]| bounded by e3t : 
    484487         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) ) 
    486490         END_3D 
    487491         DO_3DS_00_00( jpkm1, 2, -1 ) 
     
    493497      CASE ( 3 )           ! lup and ldown, |dk[xml]| bounded by e3t : 
    494498         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) ) 
    496501         END_3D 
    497502         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) ) 
    499505         END_3D 
    500506         DO_3D_00_00( 2, jpkm1 ) 
Note: See TracChangeset for help on using the changeset viewer.