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 12606 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traldf_iso.F90 – NEMO

Ignore:
Timestamp:
2020-03-26T11:15:02+01:00 (4 years ago)
Author:
techene
Message:

all: add e3 substitute and limit precompiled files lines to about 130 character, OCE/TRA/traisf.F90: remove ONLY : e3t, r1_e1e2t

File:
1 edited

Legend:

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

    r12590 r12606  
    254254            ! 
    255255            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                   &                                                 + zftv(ji,jj,jk) - zftv(ji,jj-1,jk)  )   & 
    258                   &                                             * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
     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) 
    259259            END_2D 
    260260         END DO                                        !   End of slab 
     
    301301            CASE(  1  )                            ! 1st pass : eddy coef = ah_wslp2 
    302302               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                      &           * ( pt(ji,jj,jk-1,jn) - pt(ji,jj,jk,jn) ) / e3w(ji,jj,jk,Kmm) * wmask(ji,jj,jk) 
     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) 
    306306               END_3D 
    307307            CASE(  2  )                         ! 2nd pass : eddy flux = ah_wslp2 and akz applied on pt  and pt2 gradients, resp. 
    308308               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)   & 
    310310                     &                            * (  ah_wslp2(ji,jj,jk) * ( pt (ji,jj,jk-1,jn) - pt (ji,jj,jk,jn) )   & 
    311311                     &                            +         akz(ji,jj,jk) * ( pt2(ji,jj,jk-1,jn) - pt2(ji,jj,jk,jn) )   ) 
     
    315315         ! 
    316316         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) 
    319319         END_3D 
    320320         ! 
Note: See TracChangeset for help on using the changeset viewer.