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

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

@@ -185 +185 @@
       IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'ldf_tra_init: failed to allocate arrays')
       !
+!$OMP PARALLEL WORKSHARE
       ahtu(:,:,jpk) = 0._wp                           ! last level always 0  
       ahtv(:,:,jpk) = 0._wp
+!$OMP END PARALLEL WORKSHARE
       !
       !                                               ! value of eddy mixing coef.
@@ -201 +203 @@
          CASE(   0  )      !==  constant  ==!
             IF(lwp) WRITE(numout,*) '          tracer mixing coef. = constant = ', rn_aht_0
+!$OMP PARALLEL WORKSHARE
             ahtu(:,:,:) = zah0 * umask(:,:,:)
             ahtv(:,:,:) = zah0 * vmask(:,:,:)
+!$OMP END PARALLEL WORKSHARE
             !
          CASE(  10  )      !==  fixed profile  ==!
             IF(lwp) WRITE(numout,*) '          tracer mixing coef. = F( depth )'
+!$OMP PARALLEL WORKSHARE
             ahtu(:,:,1) = zah0 * umask(:,:,1)                      ! constant surface value
             ahtv(:,:,1) = zah0 * vmask(:,:,1)
+!$OMP END PARALLEL WORKSHARE
             CALL ldf_c1d( 'TRA', r1_4, ahtu(:,:,1), ahtv(:,:,1), ahtu, ahtv )
             !
@@ -216 +222 @@
             CALL iom_get ( inum, jpdom_data, 'ahtv_2D', ahtv(:,:,1) )
             CALL iom_close( inum )
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 2, jpkm1
                ahtu(:,:,jk) = ahtu(:,:,1) * umask(:,:,jk)
@@ -245 +252 @@
             CALL iom_get ( inum, jpdom_data, 'ahtv_3D', ahtv )
             CALL iom_close( inum )
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 1, jpkm1
                ahtu(:,:,jk) = ahtu(:,:,jk) * umask(:,:,jk)
@@ -268 +276 @@
          !
          IF( ln_traldf_blp .AND. .NOT. l_ldftra_time ) THEN
+!$OMP PARALLEL WORKSHARE
             ahtu(:,:,:) = SQRT( ahtu(:,:,:) )
             ahtv(:,:,:) = SQRT( ahtv(:,:,:) )
+!$OMP END PARALLEL WORKSHARE
          ENDIF
          !
@@ -422 +432 @@
          CASE(   0  )      !==  constant  ==!
             IF(lwp) WRITE(numout,*) '          eddy induced velocity coef. = constant = ', rn_aeiv_0
+!$OMP PARALLEL WORKSHARE
             aeiu(:,:,:) = rn_aeiv_0
             aeiv(:,:,:) = rn_aeiv_0
+!$OMP END PARALLEL WORKSHARE
             !
          CASE(  10  )      !==  fixed profile  ==!
             IF(lwp) WRITE(numout,*) '          eddy induced velocity coef. = F( depth )'
+!$OMP PARALLEL WORKSHARE
             aeiu(:,:,1) = rn_aeiv_0                                ! constant surface value
             aeiv(:,:,1) = rn_aeiv_0
+!$OMP END PARALLEL WORKSHARE
             CALL ldf_c1d( 'TRA', r1_4, aeiu(:,:,1), aeiv(:,:,1), aeiu, aeiv )
             !
@@ -437 +451 @@
             CALL iom_get  ( inum, jpdom_data, 'aeiv', aeiv(:,:,1) )
             CALL iom_close( inum )
+!$OMP PARALLEL DO schedule(static) private(jk)
             DO jk = 2, jpk
                aeiu(:,:,jk) = aeiu(:,:,1)