Ignore:
Timestamp:
2017-02-18T10:02:03+01:00 (4 years ago)
Author:
mocavero
Message:

update trunk with OpenMP parallelization

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/USR/usrdef_istate.F90

    r6923 r7698  
    5555      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~   Ocean at rest, with an horizontally uniform T and S profiles' 
    5656      ! 
    57       pu  (:,:,:) = 0._wp        ! ocean at rest 
    58       pv  (:,:,:) = 0._wp 
    59       pssh(:,:)   = 0._wp 
     57!$OMP PARALLEL 
     58!$OMP DO schedule(static) private(jk,jj,ji) 
     59      DO jk = 1, jpk 
     60         DO jj = 1, jpj 
     61            DO ji = 1, jpi 
     62               pu  (ji,jj,jk) = 0._wp        ! ocean at rest 
     63               pv  (ji,jj,jk) = 0._wp 
     64            END DO 
     65         END DO 
     66      END DO 
     67!$OMP END DO NOWAIT 
     68!$OMP DO schedule(static) private(jj,ji) 
     69      DO jj = 1, jpj 
     70         DO ji = 1, jpi 
     71            pssh(ji,jj)   = 0._wp 
     72         END DO 
     73      END DO 
     74!$OMP END DO NOWAIT 
    6075      ! 
     76!$OMP DO schedule(static) private(jk,jj,ji) 
    6177      DO jk = 1, jpk             ! horizontally uniform T & S profiles 
    6278         DO jj = 1, jpj 
     
    7995         END DO 
    8096      END DO 
     97!$OMP END PARALLEL 
    8198      !    
    8299   END SUBROUTINE usr_def_istate 
Note: See TracChangeset for help on using the changeset viewer.