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 6748 for branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90 – NEMO

Ignore:
Timestamp:
2016-06-28T11:53:56+02:00 (8 years ago)
Author:
mocavero
Message:

GYRE hybrid parallelization

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90

    r6140 r6748  
    5656      IF ( PRESENT(kkk) ) ik=kkk 
    5757      SELECT CASE( cdgrid ) 
    58       CASE( 'U' )  ; zglam(:,:) = glamu(:,:) ; zgphi(:,:) = gphiu(:,:) ; zmask(nldi:nlei,nldj:nlej) = umask(nldi:nlei,nldj:nlej,ik) 
    59       CASE( 'V' )  ; zglam(:,:) = glamv(:,:) ; zgphi(:,:) = gphiv(:,:) ; zmask(nldi:nlei,nldj:nlej) = vmask(nldi:nlei,nldj:nlej,ik) 
    60       CASE( 'F' )  ; zglam(:,:) = glamf(:,:) ; zgphi(:,:) = gphif(:,:) ; zmask(nldi:nlei,nldj:nlej) = fmask(nldi:nlei,nldj:nlej,ik) 
    61       CASE DEFAULT ; zglam(:,:) = glamt(:,:) ; zgphi(:,:) = gphit(:,:) ; zmask(nldi:nlei,nldj:nlej) = tmask(nldi:nlei,nldj:nlej,ik) 
     58      CASE( 'U' ) 
     59!$OMP PARALLEL WORKSHARE 
     60      zglam(:,:) = glamu(:,:) ; zgphi(:,:) = gphiu(:,:) ; zmask(nldi:nlei,nldj:nlej) = umask(nldi:nlei,nldj:nlej,ik) 
     61!$OMP END PARALLEL WORKSHARE 
     62      CASE( 'V' )   
     63!$OMP PARALLEL WORKSHARE 
     64      zglam(:,:) = glamv(:,:) ; zgphi(:,:) = gphiv(:,:) ; zmask(nldi:nlei,nldj:nlej) = vmask(nldi:nlei,nldj:nlej,ik) 
     65!$OMP END PARALLEL WORKSHARE 
     66      CASE( 'F' )   
     67!$OMP PARALLEL WORKSHARE 
     68      zglam(:,:) = glamf(:,:) ; zgphi(:,:) = gphif(:,:) ; zmask(nldi:nlei,nldj:nlej) = fmask(nldi:nlei,nldj:nlej,ik) 
     69!$OMP END PARALLEL WORKSHARE 
     70      CASE DEFAULT  
     71!$OMP PARALLEL WORKSHARE 
     72      zglam(:,:) = glamt(:,:) ; zgphi(:,:) = gphit(:,:) ; zmask(nldi:nlei,nldj:nlej) = tmask(nldi:nlei,nldj:nlej,ik) 
     73!$OMP END PARALLEL WORKSHARE 
    6274      END SELECT 
    6375 
    … …  
    7183         zglam(:,:) = zglam(:,:) - plon 
    7284      END IF 
    73  
     85!$OMP PARALLEL WORKSHARE 
    7486      zgphi(:,:) = zgphi(:,:) - plat 
    7587      zdist(:,:) = zglam(:,:) * zglam(:,:) + zgphi(:,:) * zgphi(:,:) 
    76        
     88!$OMP END PARALLEL WORKSHARE 
    7789      IF( lk_mpp ) THEN   
    7890         CALL mpp_minloc( zdist(:,:), zmask, zmini, kii, kjj) 
Note: See TracChangeset for help on using the changeset viewer.