Changeset 7037 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.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/LIM_SRC_2/limdyn_2.F90 (modified) (5 diffs)

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

@@ -78 +78 @@
          !
          ! Mean ice and snow thicknesses.          
+!$OMP PARALLEL WORKSHARE
          hsnm(:,:)  = ( 1.0 - frld(:,:) ) * hsnif(:,:)
          hicm(:,:)  = ( 1.0 - frld(:,:) ) * hicif(:,:)
+!$OMP END PARALLEL WORKSHARE
          !
          !                                     ! Rheology (ice dynamics)
@@ -170 +172 @@
          SELECT CASE( cp_ice_msh )           ! ice-ocean relative velocity at u- & v-pts
          CASE( 'C' )                               ! EVP : C-grid ice dynamics
+!$OMP PARALLEL WORKSHARE
             zu_io(:,:) = u_ice(:,:) - ssu_m(:,:)           ! ice-ocean & ice velocity at ocean velocity points
             zv_io(:,:) = v_ice(:,:) - ssv_m(:,:)
+!$OMP END PARALLEL WORKSHARE
          CASE( 'I' )                               ! VP  : B-grid ice dynamics (I-point) 
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
             DO jj = 1, jpjm1                               ! u_ice v_ice at I-point ; ssu_m, ssv_m at U- & V-points
                DO ji = 1, jpim1   ! NO vector opt.         !
@@ -183 +188 @@
          ! frictional velocity at T-point
          zcoef = 0.5_wp * cw
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1
             DO ji = 2, jpim1   ! NO vector opt. because of zu_io
@@ -193 +199 @@
          !
          zcoef = SQRT( 0.5 ) / rau0
+!$OMP PARALLEL DO schedule(static) private(jj, ji)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -278 +285 @@
       pstarh = pstar / 2.0
       !
+!$OMP PARALLEL WORKSHARE
       ahiu(:,:) = ahi0 * umask(:,:,1)            ! Ice eddy Diffusivity coefficients.
       ahiv(:,:) = ahi0 * vmask(:,:,1)
+!$OMP END PARALLEL WORKSHARE
       !
    END SUBROUTINE lim_dyn_init_2