Changeset 7037 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

Timestamp:

2016-10-18T15:32:04+02:00 (8 years ago)

Author:

mocavero

Message:

ORCA2_LIM_PISCES hybrid version update

File:

• branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90 (modified) (6 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

@@ -109 +109 @@
          ELSE                                         ! No restart or restart not found: Euler forward time stepping
            zfact = 1._wp
+!$OMP PARALLEL WORKSHARE
            sbc_trc_b(:,:,:) = 0._wp
+!$OMP END PARALLEL WORKSHARE
          ENDIF
       ELSE                                         ! Swap of forcing fields
          IF( ln_top_euler ) THEN
             zfact = 1._wp
+!$OMP PARALLEL WORKSHARE
             sbc_trc_b(:,:,:) = 0._wp
+!$OMP END PARALLEL WORKSHARE
          ELSE
             zfact = 0.5_wp
+!$OMP PARALLEL WORKSHARE
             sbc_trc_b(:,:,:) = sbc_trc(:,:,:)
+!$OMP END PARALLEL WORKSHARE
          ENDIF
          !
@@ -127 +133 @@
       !
       IF( .NOT. lk_offline .AND. .NOT.ln_linssh ) THEN  ! online coupling with vvl
+!$OMP PARALLEL WORKSHARE
          zsfx(:,:) = 0._wp
+!$OMP END PARALLEL WORKSHARE
       ELSE                                      ! online coupling free surface or offline with free surface
+!$OMP PARALLEL WORKSHARE
          zsfx(:,:) = emp(:,:)
+!$OMP END PARALLEL WORKSHARE
       ENDIF
 
@@ -135 +145 @@
       DO jn = 1, jptra
          !
-         IF( l_trdtrc )   ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends
-
+         IF( l_trdtrc ) THEN
+!$OMP PARALLEL WORKSHARE
+            ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends
+!$OMP END PARALLEL WORKSHARE
+         END IF
          IF ( nn_ice_tr == -1 ) THEN  ! No tracers in sea ice (null concentration in sea ice)
 
+!$OMP PARALLEL DO schedule(static) private(jj, ji) 
             DO jj = 2, jpj
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -146 +160 @@
 
          ELSE
-
+!$OMP PARALLEL DO schedule(static) private(jj,ji,zse3t,zftra,zcd,ztfx,zdtra,zratio)
             DO jj = 2, jpj
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -170 +184 @@
          CALL lbc_lnk( sbc_trc(:,:,jn), 'T', 1. )
          !                                       Concentration dilution effect on tracers due to evaporation & precipitation 
+!$OMP PARALLEL DO schedule(static) private(jj,ji,zse3t) 
          DO jj = 2, jpj
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -178 +193 @@
          !
          IF( l_trdtrc ) THEN
+!$OMP PARALLEL WORKSHARE
             ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:)
+!$OMP END PARALLEL WORKSHARE
             CALL trd_tra( kt, 'TRC', jn, jptra_nsr, ztrtrd )
          END IF