Changeset 6748 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90

Timestamp:

2016-06-28T11:53:56+02:00 (8 years ago)

Author:

mocavero

Message:

GYRE hybrid parallelization

File:

• branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90 (modified) (6 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90

@@ -77 +77 @@
       IF( l_trddyn )   THEN         ! Save ua and va trends
          CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
+!$OMP PARALLEL WORKSHARE 
          ztrdu(:,:,:) = ua(:,:,:) 
          ztrdv(:,:,:) = va(:,:,:) 
+!$OMP END PARALLEL WORKSHARE 
       ENDIF
       
       DO jk = 2, jpkm1              ! Vertical momentum advection at level w and u- and v- vertical
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpj                   ! vertical fluxes 
             DO ji = fs_2, jpi             ! vector opt.
@@ -87 +90 @@
             END DO
          END DO
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1                 ! vertical momentum advection at w-point
             DO ji = fs_2, fs_jpim1        ! vector opt.
@@ -97 +101 @@
       ! Surface and bottom advective fluxes set to zero
       IF ( ln_isfcav ) THEN
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1           ! vector opt.
@@ -106 +111 @@
          END DO
       ELSE
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1        
             DO ji = fs_2, fs_jpim1           ! vector opt.
@@ -116 +122 @@
       END IF
 
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji, zua, zva)
       DO jk = 1, jpkm1              ! Vertical momentum advection at u- and v-points
          DO jj = 2, jpjm1
@@ -130 +137 @@
 
       IF( l_trddyn ) THEN           ! save the vertical advection trends for diagnostic
+!$OMP PARALLEL WORKSHARE 
          ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
          ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
+!$OMP END PARALLEL WORKSHARE 
          CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt )
          CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv )