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 15020 for NEMO/branches/2021/ticket2680_C1D_PAPA/src/OCE/BDY/bdyini.F90 – NEMO

Ignore:
Timestamp:
2021-06-18T15:21:42+02:00 (3 years ago)
Author:
gsamson
Message:

merge trunk into branch (#2680)

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/ticket2680_C1D_PAPA/src/OCE/BDY/bdyini.F90

    r14433 r15020  
    4444   INTEGER, DIMENSION(jp_nseg) ::   jpjnob, jpindt, jpinft, npckgn   ! 
    4545   INTEGER, DIMENSION(jp_nseg) ::   jpjsob, jpisdt, jpisft, npckgs   ! 
     46    
     47   !! * Substitutions 
     48#  include "do_loop_substitute.h90" 
    4649   !!---------------------------------------------------------------------- 
    4750   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    … …  
    140143      !! ** Input   :  bdy_init.nc, input file for unstructured open boundaries 
    141144      !!----------------------------------------------------------------------       
     145      INTEGER  ::   ji, jj                                 ! dummy loop indices 
    142146      INTEGER  ::   ib_bdy, ii, ij, igrd, ib, ir, iseg     ! dummy loop indices 
    143147      INTEGER  ::   icount, icountr, icountr0, ibr_max     ! local integers 
    … …  
    630634      ! For the flagu/flagv calculation below we require a version of fmask without 
    631635      ! the land boundary condition (shlat) included: 
    632       DO ij = 1, jpjm1 
    633          DO ii = 1, jpim1 
    634             zfmask(ii,ij) =  ztmask(ii,ij  ) * ztmask(ii+1,ij  )   & 
    635                &           * ztmask(ii,ij+1) * ztmask(ii+1,ij+1) 
    636          END DO 
    637       END DO 
     636      DO_2D( 0, 0, 0, 0 ) 
     637         zfmask(ji,jj) =  ztmask(ji,jj  ) * ztmask(ji+1,jj  )   & 
     638            &           * ztmask(ji,jj+1) * ztmask(ji+1,jj+1) 
     639      END_2D 
    638640      CALL lbc_lnk( 'bdyini', zfmask, 'F', 1.0_wp ) 
    639641 
    … …  
    646648 
    647649      ! Derive mask on U and V grid from mask on T grid 
    648       DO ij = 1, jpjm1 
    649          DO ii = 1, jpim1 
    650             bdyumask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii+1,ij  ) 
    651             bdyvmask(ii,ij) = bdytmask(ii,ij) * bdytmask(ii  ,ij+1)   
    652          END DO 
    653       END DO 
     650      DO_2D( 0, 0, 0, 0 ) 
     651            bdyumask(ji,jj) = bdytmask(ji,jj) * bdytmask(ji+1,jj  ) 
     652            bdyvmask(ji,jj) = bdytmask(ji,jj) * bdytmask(ji  ,jj+1)   
     653      END_2D 
    654654      CALL lbc_lnk( 'bdyini', bdyumask, 'U', 1.0_wp , bdyvmask, 'V', 1.0_wp )   ! Lateral boundary cond.  
    655655 
    … …  
    687687 
    688688      ! Recompute zfmask 
    689       DO ij = 1, jpjm1 
    690          DO ii = 1, jpim1 
    691             zfmask(ii,ij) =  ztmask(ii,ij  ) * ztmask(ii+1,ij  )   & 
    692                &           * ztmask(ii,ij+1) * ztmask(ii+1,ij+1) 
    693          END DO 
    694       END DO 
     689      DO_2D( 0, 0, 0, 0 ) 
     690         zfmask(ji,jj) =  ztmask(ji,jj  ) * ztmask(ji+1,jj  )   & 
     691            &           * ztmask(ji,jj+1) * ztmask(ji+1,jj+1) 
     692      END_2D 
    695693      CALL lbc_lnk( 'bdyini', zfmask, 'F', 1.0_wp ) 
    696694 
Note: See TracChangeset for help on using the changeset viewer.