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

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

@@ -83 +83 @@
       IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends
          CALL wrk_alloc( jpi,jpj,jpk,   ztrdu, ztrdv ) 
+!$OMP PARALLEL WORKSHARE
          ztrdu(:,:,:) = ua(:,:,:)
          ztrdv(:,:,:) = va(:,:,:)
+!$OMP END PARALLEL WORKSHARE
       ENDIF
       !
@@ -91 +93 @@
          .OR.  nn_ice_embd == 2  ) THEN                                      ! embedded sea-ice
          !
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -100 +103 @@
          IF( ln_apr_dyn .AND. .NOT.ln_dynspg_ts ) THEN   !==  Atmospheric pressure gradient (added later in time-split case) ==!
             zg_2 = grav * 0.5
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
             DO jj = 2, jpjm1                          ! gradient of Patm using inverse barometer ssh
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -115 +119 @@
             CALL upd_tide( kt )                      ! update tide potential
             !
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
             DO jj = 2, jpjm1                         ! add tide potential forcing
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -129 +134 @@
             zgrau0r     = - grav * r1_rau0
             zpice(:,:) = (  zintp * snwice_mass(:,:) + ( 1.- zintp ) * snwice_mass_b(:,:)  ) * zgrau0r
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
             DO jj = 2, jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -139 +145 @@
          ENDIF
          !
+!$OMP PARALLEL DO schedule(static) private(jk, jj, ji)
          DO jk = 1, jpkm1                    !== Add all terms to the general trend
             DO jj = 2, jpjm1
@@ -158 +165 @@
       !                    
       IF( l_trddyn )   THEN                  ! save the surface pressure gradient trends for further diagnostics
+!$OMP PARALLEL WORKSHARE
          ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
          ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
+!$OMP END PARALLEL WORKSHARE
          CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt )
          CALL wrk_dealloc( jpi,jpj,jpk,   ztrdu, ztrdv )