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 13469 for NEMO/branches/2020/temporary_r4_trunk/src/ICE/iceitd.F90 – NEMO

Ignore:
Timestamp:
2020-09-15T12:49:18+02:00 (4 years ago)
Author:
smasson
Message:

r4_trunk: first change of DO loops for routines to be merged, see #2523

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/temporary_r4_trunk/src/ICE/iceitd.F90

    r13466 r13469  
    9797      ! 
    9898      npti = 0   ;   nptidx(:) = 0 
    99       DO jj = 1, jpj 
    100          DO ji = 1, jpi 
    101             IF ( at_i(ji,jj) > epsi10 ) THEN 
    102                npti = npti + 1 
    103                nptidx( npti ) = (jj - 1) * jpi + ji 
    104             ENDIF 
    105          END DO 
    106       END DO 
     99      DO_2D_11_11 
     100         IF ( at_i(ji,jj) > epsi10 ) THEN 
     101            npti = npti + 1 
     102            nptidx( npti ) = (jj - 1) * jpi + ji 
     103         ENDIF 
     104      END_2D 
    107105       
    108106      !----------------------------------------------------------------------------------------------- 
    … …  
    606604         !                    !--------------------------------------- 
    607605         npti = 0   ;   nptidx(:) = 0 
    608          DO jj = 1, jpj 
    609             DO ji = 1, jpi 
    610                IF( a_i(ji,jj,jl) > 0._wp .AND. v_i(ji,jj,jl) > (a_i(ji,jj,jl) * hi_max(jl)) ) THEN 
    611                   npti = npti + 1 
    612                   nptidx( npti ) = (jj - 1) * jpi + ji                   
    613                ENDIF 
    614             END DO 
    615          END DO 
     606         DO_2D_11_11 
     607            IF( a_i(ji,jj,jl) > 0._wp .AND. v_i(ji,jj,jl) > (a_i(ji,jj,jl) * hi_max(jl)) ) THEN 
     608               npti = npti + 1 
     609               nptidx( npti ) = (jj - 1) * jpi + ji                   
     610            ENDIF 
     611         END_2D 
    616612         ! 
    617613!!clem   CALL tab_2d_1d( npti, nptidx(1:npti), h_i_1d(1:npti), h_i(:,:,jl) ) 
    … …  
    647643         !                    !----------------------------------------- 
    648644         npti = 0 ; nptidx(:) = 0 
    649          DO jj = 1, jpj 
    650             DO ji = 1, jpi 
    651                IF( a_i(ji,jj,jl+1) > 0._wp .AND. v_i(ji,jj,jl+1) <= (a_i(ji,jj,jl+1) * hi_max(jl)) ) THEN 
    652                   npti = npti + 1 
    653                   nptidx( npti ) = (jj - 1) * jpi + ji                   
    654                ENDIF 
    655             END DO 
    656          END DO 
     645         DO_2D_11_11 
     646            IF( a_i(ji,jj,jl+1) > 0._wp .AND. v_i(ji,jj,jl+1) <= (a_i(ji,jj,jl+1) * hi_max(jl)) ) THEN 
     647               npti = npti + 1 
     648               nptidx( npti ) = (jj - 1) * jpi + ji                   
     649            ENDIF 
     650         END_2D 
    657651         ! 
    658652         CALL tab_2d_1d( npti, nptidx(1:npti), a_i_1d(1:npti), a_i(:,:,jl+1) ) ! jl+1 is ok 
Note: See TracChangeset for help on using the changeset viewer.