New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7508 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90 – NEMO

Ignore:
Timestamp:
2016-12-19T13:15:59+01:00 (8 years ago)
Author:
mocavero
Message:

changes on code duplication and workshare construct

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfmxl.F90

    r7037 r7508  
    9999      ! w-level of the mixing and mixed layers 
    100100!$OMP PARALLEL 
    101 !$OMP WORKSHARE 
    102       nmln(:,:)  = nlb10               ! Initialization to the number of w ocean point 
    103       hmlp(:,:)  = 0._wp               ! here hmlp used as a dummy variable, integrating vertically N^2 
    104 !$OMP END WORKSHARE 
     101!$OMP DO schedule(static) private(jj, ji) 
     102      DO jj = 1, jpj 
     103         DO ji = 1, jpi 
     104            nmln(ji,jj)  = nlb10               ! Initialization to the number of w ocean point 
     105            hmlp(ji,jj)  = 0._wp               ! here hmlp used as a dummy variable, integrating vertically N^2 
     106         END DO 
     107      END DO 
    105108      DO jk = nlb10, jpkm1 
    106109!$OMP DO schedule(static) private(jj, ji, ikt) 
     
    115118      ! 
    116119      ! w-level of the turbocline and mixing layer (iom_use) 
    117 !$OMP WORKSHARE 
    118       imld(:,:) = mbkt(:,:) + 1        ! Initialization to the number of w ocean point 
    119 !$OMP END WORKSHARE 
     120!$OMP DO schedule(static) private(jj, ji) 
     121      DO jj = 1, jpj 
     122         DO ji = 1, jpi 
     123            imld(ji,jj) = mbkt(ji,jj) + 1        ! Initialization to the number of w ocean point 
     124         END DO 
     125      END DO 
    120126 
    121127      DO jk = jpkm1, nlb10, -1         ! from the bottom to nlb10  
Note: See TracChangeset for help on using the changeset viewer.