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 15033 for NEMO/trunk/src/OCE/DYN – NEMO

Ignore:
Timestamp:
2021-06-21T12:24:45+02:00 (3 years ago)
Author:
smasson
Message:

trunk: suppress jpim1 et jpjm1, #2699

Location:
NEMO/trunk/src/OCE/DYN
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/OCE/DYN/dynldf_lap_blp.F90

    r14834 r15033  
    114114            ! 
    115115            DO_2D( iij-1, iij, iij-1, iij ) 
    116                !                                      ! ahm * e3 * curl  (computed from 1 to jpim1/jpjm1) 
     116               !                                      ! ahm * e3 * curl  (warning: computed for ji-1,jj-1) 
    117117               zcur(ji-1,jj-1) = ahmf(ji-1,jj-1,jk) * e3f(ji-1,jj-1,jk) * r1_e1e2f(ji-1,jj-1)       &   ! ahmf already * by fmask 
    118118                  &     * (  e2v(ji  ,jj-1) * pv(ji  ,jj-1,jk) - e2v(ji-1,jj-1) * pv(ji-1,jj-1,jk)  & 
    119119                  &        - e1u(ji-1,jj  ) * pu(ji-1,jj  ,jk) + e1u(ji-1,jj-1) * pu(ji-1,jj-1,jk)  ) 
    120                !                                      ! ahm * div        (computed from 2 to jpi/jpj) 
     120               !                                      ! ahm * div        (warning: computed for ji,jj) 
    121121               zdiv(ji,jj)     = ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kbb)               &   ! ahmt already * by tmask 
    122122                  &     * (  e2u(ji,jj)*e3u(ji,jj,jk,Kbb) * pu(ji,jj,jk) - e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kbb) * pu(ji-1,jj,jk)  & 

  • NEMO/trunk/src/OCE/DYN/dynldf_lap_blp_lf.F90

    r14834 r15033  
    101101         ! 
    102102         DO_3D( iij-1, iij-1, iij-1, iij-1, 1, jpkm1 )                           ! Horizontal slab 
    103             !                                      ! ahm * e3 * curl  (computed from 1 to jpim1/jpjm1) 
     103            !                                      ! ahm * e3 * curl  (warning: computed for ji-1,jj-1) 
    104104            zcur     = ahmf(ji,jj,jk) * e3f(ji,jj,jk) * r1_e1e2f(ji,jj)               &   ! ahmf already * by fmask    
    105105               &       * ( e2v(ji+1,jj) * pv(ji+1,jj,jk) - e2v(ji,jj) * pv(ji,jj,jk)  & 
    … …  
    111111               &       * ( e2v(ji,jj) * pv(ji,jj,jk) - e2v(ji-1,jj) * pv(ji-1,jj,jk)  & 
    112112               &       - e1u(ji-1,jj+1) * pu(ji-1,jj+1,jk) + e1u(ji-1,jj) * pu(ji-1,jj,jk) ) 
    113             !                                      ! ahm * div        (computed from 2 to jpi/jpj) 
     113            !                                      ! ahm * div        (warning: computed for ji,jj) 
    114114            zdiv     = ahmt(ji,jj,jk) * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kbb)               &   ! ahmt already * by tmask 
    115115               &     * ( e2u(ji,jj)*e3u(ji,jj,jk,Kbb) * pu(ji,jj,jk) - e2u(ji-1,jj)*e3u(ji-1,jj,jk,Kbb) * pu(ji-1,jj,jk)  & 
Note: See TracChangeset for help on using the changeset viewer.