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

branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

@@ -136 +136 @@
       IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'ldf_dyn_init: failed to allocate arrays')
       !
+!$OMP PARALLEL WORKSHARE
       ahmt(:,:,jpk) = 0._wp                           ! last level always 0  
       ahmf(:,:,jpk) = 0._wp
+!$OMP END PARALLEL WORKSHARE
       !
       !                                               ! value of eddy mixing coef.
@@ -154 +156 @@
          CASE(   0  )      !==  constant  ==!
             IF(lwp) WRITE(numout,*) '          momentum mixing coef. = constant '
+!$OMP PARALLEL WORKSHARE
             ahmt(:,:,:) = zah0 * tmask(:,:,:)
             ahmf(:,:,:) = zah0 * fmask(:,:,:)
+!$OMP END PARALLEL WORKSHARE
             !
          CASE(  10  )      !==  fixed profile  ==!
             IF(lwp) WRITE(numout,*) '          momentum mixing coef. = F( depth )'
+!$OMP PARALLEL WORKSHARE
             ahmt(:,:,1) = zah0 * tmask(:,:,1)                      ! constant surface value
             ahmf(:,:,1) = zah0 * fmask(:,:,1)
+!$OMP END PARALLEL WORKSHARE
             CALL ldf_c1d( 'DYN', r1_4, ahmt(:,:,1), ahmf(:,:,1), ahmt, ahmf )
             !
@@ -172 +178 @@
 !!              do we introduce a scaling by the max value of the array, and then multiply by zah0 ????
 !!              better:  check that the max is <=1  i.e. it is a shape from 0 to 1, not a coef that has physical dimension
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 2, jpkm1
                ahmt(:,:,jk) = ahmt(:,:,1) * tmask(:,:,jk)
@@ -190 +197 @@
 !!gm Question : info for LAP or BLP case  to take into account the SQRT in the bilaplacian case ????
 !!              do we introduce a scaling by the max value of the array, and then multiply by zah0 ????
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 1, jpkm1
                ahmt(:,:,jk) = ahmt(:,:,jk) * tmask(:,:,jk)