Changeset 15014 for NEMO/trunk/src/OCE/BDY
- Timestamp:
- 2021-06-17T19:02:04+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/BDY/bdyini.F90
r14433 r15014 44 44 INTEGER, DIMENSION(jp_nseg) :: jpjnob, jpindt, jpinft, npckgn ! 45 45 INTEGER, DIMENSION(jp_nseg) :: jpjsob, jpisdt, jpisft, npckgs ! 46 47 !! * Substitutions 48 # include "do_loop_substitute.h90" 46 49 !!---------------------------------------------------------------------- 47 50 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 140 143 !! ** Input : bdy_init.nc, input file for unstructured open boundaries 141 144 !!---------------------------------------------------------------------- 145 INTEGER :: ji, jj ! dummy loop indices 142 146 INTEGER :: ib_bdy, ii, ij, igrd, ib, ir, iseg ! dummy loop indices 143 147 INTEGER :: icount, icountr, icountr0, ibr_max ! local integers … … 630 634 ! For the flagu/flagv calculation below we require a version of fmask without 631 635 ! 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 638 640 CALL lbc_lnk( 'bdyini', zfmask, 'F', 1.0_wp ) 639 641 … … 646 648 647 649 ! 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 654 654 CALL lbc_lnk( 'bdyini', bdyumask, 'U', 1.0_wp , bdyvmask, 'V', 1.0_wp ) ! Lateral boundary cond. 655 655 … … 687 687 688 688 ! 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 695 693 CALL lbc_lnk( 'bdyini', zfmask, 'F', 1.0_wp ) 696 694
Note: See TracChangeset
for help on using the changeset viewer.