Ignore:
Timestamp:
2017-03-03T12:46:59+01:00 (4 years ago)
Author:
mocavero
Message:

Reverting trunk to remove OpenMP

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r7698 r7753  
    105105      !! ** Purpose :      passive tracers inventories at initialsation phase 
    106106      !!---------------------------------------------------------------------- 
    107       INTEGER ::  jk, jn, jj, ji    ! dummy loop indices 
     107      INTEGER ::  jk, jn    ! dummy loop indices 
    108108      CHARACTER (len=25) :: charout 
    109109      !!---------------------------------------------------------------------- 
    110110      !                                                              ! masked grid volume 
    111 !$OMP PARALLEL 
    112 !$OMP DO schedule(static) private(jk,jj,ji) 
    113111      DO jk = 1, jpk 
    114          DO jj = 1, jpj 
    115             DO ji = 1, jpi 
    116                cvol(ji,jj,jk) = e1e2t(ji,jj) * e3t_n(ji,jj,jk) * tmask(ji,jj,jk) 
    117             END DO 
    118          END DO 
     112         cvol(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk) 
    119113      END DO 
    120       ! 
    121 !$OMP DO schedule(static) private(jn) 
    122       DO jn = 1, jptra 
    123          trai(jn) = 0._wp                                               ! initial content of all tracers 
    124       END DO 
    125 !$OMP END PARALLEL 
    126114      !                                                              ! total volume of the ocean  
    127115      areatot = glob_sum( cvol(:,:,:) ) 
    128116      ! 
     117      trai(:) = 0._wp                                                   ! initial content of all tracers 
    129118      DO jn = 1, jptra 
    130119         trai(jn) = trai(jn) + glob_sum( trn(:,:,:,jn) * cvol(:,:,:)   ) 
     
    231220      USE trcdta          ! initialisation from files 
    232221      ! 
    233       INTEGER :: jn, jl, jk, jj, ji   ! dummy loop indices 
     222      INTEGER :: jn, jl   ! dummy loop indices 
    234223      !!---------------------------------------------------------------------- 
    235224      ! 
     
    265254        ENDIF 
    266255        ! 
    267 !$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji) 
    268         DO jn = 1, jptra 
    269            DO jk = 1, jpk 
    270               DO jj = 1, jpj 
    271                  DO ji = 1, jpi 
    272                     trb(ji,jj,jk,jn) = trn(ji,jj,jk,jn) 
    273                  END DO 
    274               END DO 
    275            END DO 
    276         END DO 
     256        trb(:,:,:,:) = trn(:,:,:,:) 
    277257        !  
    278258      ENDIF 
    279259  
    280 !$OMP PARALLEL DO schedule(static) private(jn,jk,jj,ji) 
    281       DO jn = 1, jptra 
    282          DO jk = 1, jpk 
    283             DO jj = 1, jpj 
    284                DO ji = 1, jpi 
    285                   tra(ji,jj,jk,jn) = 0._wp 
    286                END DO 
    287             END DO 
    288          END DO 
    289       END DO 
     260      tra(:,:,:,:) = 0._wp 
    290261      !                                                         ! Partial top/bottom cell: GRADh(trn) 
    291262   END SUBROUTINE trc_ini_state 
Note: See TracChangeset for help on using the changeset viewer.