Changeset 7698 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90

Timestamp:

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

Author:

mocavero

Message:

update trunk with OpenMP parallelization

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90

@@ -69 +69 @@
       CALL wrk_alloc( jpi, jpj, jpk, zco3, zcaldiss, zhinit, zhi, zco3sat )
       !
-      zco3    (:,:,:) = 0.
-      zcaldiss(:,:,:) = 0.
-      zhinit(:,:,:)   = hi(:,:,:) * 1000. / ( rhop(:,:,:) + rtrn )
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
+       DO jk = 1, jpk
+          DO jj = 1, jpj
+             DO ji = 1, jpi
+                zco3    (ji,jj,jk) = 0.
+                zcaldiss(ji,jj,jk) = 0.
+                zhinit(ji,jj,jk)   = hi(ji,jj,jk) * 1000. / ( rhop(ji,jj,jk) + rtrn )
+             END DO
+          END DO
+      END DO
       !     -------------------------------------------
       !     COMPUTE [CO3--] and [H+] CONCENTRATIONS
@@ -78 +85 @@
       CALL solve_at_general(zhinit, zhi)
 
+!$OMP PARALLEL
+!$OMP DO schedule(static) private(jk, jj, ji)
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -94 +103 @@
       !     ---------------------------------------------------------
 
+!$OMP DO schedule(static) private(jk,jj,ji,zcalcon,zfact,zomegaca,zexcess0,zexcess,zdispot)
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -124 +134 @@
          END DO
       END DO
+!$OMP END PARALLEL
       !