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

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

@@ -91 +91 @@
       IF( l_trddyn ) THEN                    ! Temporary saving of ua and va trends (l_trddyn)
          CALL wrk_alloc( jpi,jpj,jpk, ztrdu, ztrdv )
+!$OMP PARALLEL WORKSHARE
          ztrdu(:,:,:) = ua(:,:,:)
          ztrdv(:,:,:) = va(:,:,:)
+!$OMP END PARALLEL WORKSHARE
       ENDIF
       !
@@ -105 +107 @@
       !
       IF( l_trddyn ) THEN      ! save the hydrostatic pressure gradient trends for momentum trend diagnostics
+!$OMP PARALLEL WORKSHARE
          ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
          ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
+!$OMP END PARALLEL WORKSHARE
          CALL trd_dyn( ztrdu, ztrdv, jpdyn_hpg, kt )
          CALL wrk_dealloc( jpi,jpj,jpk, ztrdu, ztrdv )
@@ -238 +242 @@
             END DO
          END DO
+!$OMP PARALLEL WORKSHARE
          riceload(:,:)=ziceload(:,:)  ! need to be saved for diaar5
+!$OMP END PARALLEL WORKSHARE
 
          CALL wrk_dealloc( jpi,jpj, 2,  ztstop) 
@@ -282 +288 @@
 
       ! Surface value
+!$OMP PARALLEL DO private(ji,jj, zcoef1)
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -297 +304 @@
       ! interior value (2=<jk=<jpkm1)
       DO jk = 2, jpkm1
+!$OMP PARALLEL DO private(ji,jj, zcoef1)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.