Changeset 6748 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DYN/dynbfr.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/dynbfr.F90 (modified) (3 diffs)

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

@@ -65 +65 @@
         IF( l_trddyn ) THEN      ! trends: store the input trends
            CALL wrk_alloc( jpi,jpj,jpk,   ztrdu, ztrdv )
+!$OMP PARALLEL WORKSHARE
            ztrdu(:,:,:) = ua(:,:,:)
            ztrdv(:,:,:) = va(:,:,:)
+!$OMP END PARALLEL WORKSHARE
         ENDIF
 
-
+!$OMP PARALLEL DO schedule(static) private(jj, ji, ikbu, ikbv)
         DO jj = 2, jpjm1
            DO ji = 2, jpim1
@@ -82 +84 @@
         !
         IF( ln_isfcav ) THEN        ! ocean cavities
+!$OMP PARALLEL DO schedule(static) private(jj, ji, ikbu, ikbv)
            DO jj = 2, jpjm1
               DO ji = 2, jpim1
@@ -99 +102 @@
         !
         IF( l_trddyn ) THEN      ! trends: send trends to trddyn for further diagnostics
+!$OMP PARALLEL WORKSHARE
            ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
            ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
+!$OMP END PARALLEL WORKSHARE
            CALL trd_dyn( ztrdu(:,:,:), ztrdv(:,:,:), jpdyn_bfr, kt )
            CALL wrk_dealloc( jpi,jpj,jpk,   ztrdu, ztrdv )