Changeset 6508 for branches/2016/dev_v3_6_STABLE_OMP/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

Timestamp:

2016-05-04T16:23:46+02:00 (8 years ago)

Author:

dkuts

Message:

First version of OMP changes, partly ported from previous branch

File:

• branches/2016/dev_v3_6_STABLE_OMP/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90 (modified) (10 diffs)

branches/2016/dev_v3_6_STABLE_OMP/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

@@ -142 +142 @@
          !
          IF( ln_dynadv_vec ) THEN      ! vector form : applied on velocity
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
             DO jk = 1, jpkm1
                DO jj = 2, jpjm1
@@ -152 +153 @@
             !
          ELSE                          ! flux form : applied on thickness weighted velocity
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
             DO jk = 1, jpkm1
                DO jj = 2, jpjm1
@@ -175 +177 @@
             END DO 
          END DO 
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
          DO jk = 1, jpkm1              ! unweighted time stepping 
             DO jj = 2, jpjm1
@@ -186 +189 @@
          IF( l_trddyn )   THEN                      ! temporary save of spg trends
             CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv )
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
             DO jk = 1, jpkm1              ! unweighted time stepping 
                DO jj = 2, jpjm1
@@ -210 +214 @@
       ! compute the next vertically averaged velocity (effect of the additional force not included)
       ! ---------------------------------------------
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -217 +222 @@
       END DO
       DO jk = 2, jpkm1                     ! vertical sum
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -224 +230 @@
          END DO
       END DO
-
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
       DO jj = 2, jpjm1                     ! transport: multiplied by the horizontal scale factor
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -238 +244 @@
       ! Right hand side of the elliptic equation and first guess
       ! --------------------------------------------------------
+!$OMP PARALLEL DO schedule(static) private(jj, ji, zgcb,zbtd)
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -296 +303 @@
       ! Transport divergence gradient multiplied by z2dt
       ! --------------------------------------------====
+!$OMP PARALLEL DO schedule(static) private(jj, ji, ztdgu, ztdgv)
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -354 +362 @@
       !                       trend, the leap-frog time stepping will not
       !                       be done in dynnxt.F90 routine)
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
       DO jk = 1, jpkm1
          DO jj = 2, jpjm1