Changeset 7508 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90

Timestamp:

2016-12-19T13:15:59+01:00 (8 years ago)

Author:

mocavero

Message:

changes on code duplication and workshare construct

File:

• branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90 (modified) (2 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/lib_fortran.F90

@@ -103 +103 @@
       !!-----------------------------------------------------------------------
       REAL(wp), INTENT(in), DIMENSION(:,:,:) ::   ptab          ! input 3D array
-      REAL(wp)                               ::   glob_sum_3d   ! global masked sum
-      !!
-      INTEGER :: jk
+      REAL(wp)                               ::   glob_sum_3d, a, b   ! global masked sum
+      !!
+      INTEGER :: jk, jj, ji
       INTEGER :: ijpk ! local variable: size of the 3d dimension of ptab
       !!-----------------------------------------------------------------------
@@ -112 +112 @@
       !
       glob_sum_3d = 0.e0
+!$OMP PARALLEL
+!$OMP DO schedule(static) private(jk,jj,ji,b) REDUCTION(+:a,glob_sum_3d)
       DO jk = 1, ijpk
-         glob_sum_3d = glob_sum_3d + SUM( ptab(:,:,jk)*tmask_i(:,:) )
-      END DO
+      a = 0.e0
+      b = 0.e0
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               b = ptab(ji,jj,jk)*tmask_i(ji,jj)
+               a = a + b
+           END DO
+         END DO
+         glob_sum_3d = glob_sum_3d + a
+      END DO
+!$OMP END PARALLEL
+
       IF( lk_mpp )   CALL mpp_sum( glob_sum_3d )
       !