Changeset 7698 for trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90

Timestamp:

2017-02-18T10:02:03+01:00 (7 years ago)

Author:

mocavero

Message:

update trunk with OpenMP parallelization

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynldf_lap_blp.F90

@@ -75 +75 @@
       !
       !                                                ! ===============
+!$OMP PARALLEL
       DO jk = 1, jpkm1                                 ! Horizontal slab
          !                                             ! ===============
+!$OMP DO schedule(static) private(jj, ji)
          DO jj = 2, jpj
             DO ji = fs_2, jpi   ! vector opt.
@@ -93 +95 @@
          END DO  
          !
+!$OMP DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1                             ! - curl( curl) + grad( div )
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -106 +109 @@
          !                                             ! ===============
       END DO                                           !   End of slab
+!$OMP END PARALLEL
       !                                                ! ===============
       CALL wrk_dealloc( jpi, jpj, zcur, zdiv ) 
@@ -128 +132 @@
       !!----------------------------------------------------------------------
       INTEGER                         , INTENT(in   ) ::   kt         ! ocean time-step index
+      INTEGER                                         ::   jk, jj, ji
       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in   ) ::   pub, pvb   ! before velocity fields
       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pua, pva   ! momentum trend
@@ -144 +149 @@
       ENDIF
       !
-      zulap(:,:,:) = 0._wp
-      zvlap(:,:,:) = 0._wp
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
+      DO jk = 1, jpk
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               zulap(ji,jj,jk) = 0._wp
+               zvlap(ji,jj,jk) = 0._wp
+            END DO
+         END DO
+      END DO
       !
       CALL dyn_ldf_lap( kt, pub, pvb, zulap, zvlap, 1 )   ! rotated laplacian applied to ptb (output in zlap)