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

Ignore:
Timestamp:
2016-10-18T15:32:04+02:00 (7 years ago)
Author:
mocavero
Message:

ORCA2_LIM_PISCES hybrid version update

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_r6519_HPC_4/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

    r6748 r7037  
    325325! need to be like this to compute the pressure gradient with ISF. If not, level beneath the ISF are not aligned (sum(e3t) /= depth) 
    326326! define e3t_0 and e3w_0 as the differences between gdept and gdepw respectively 
    327 !$OMP PARALLEL DO schedule(static) private(jk) 
    328327         DO jk = 1, jpkm1 
    329328            e3t_1d(jk) = gdepw_1d(jk+1)-gdepw_1d(jk)  
    330329         END DO 
    331330         e3t_1d(jpk) = e3t_1d(jpk-1)   ! we don't care because this level is masked in NEMO 
    332 !$OMP PARALLEL DO schedule(static) private(jk) 
    333331         DO jk = 2, jpk 
    334332            e3w_1d(jk) = gdept_1d(jk) - gdept_1d(jk-1)  
     
    429427!$OMP END PARALLEL WORKSHARE 
    430428               ELSE                                                ! z-coordinate (zco or zps): step-like topography 
    431 !$OMP PARALLEL WORKSHARE 
     429!$OMP PARALLEL 
     430!$OMP WORKSHARE 
    432431                  idta(:,:) = jpkm1 
    433 !$OMP END PARALLEL WORKSHARE 
     432!$OMP END WORKSHARE 
    434433                  DO jk = 1, jpkm1 
     434!$OMP WORKSHARE 
    435435                     WHERE( gdept_1d(jk) < zdta(:,:) .AND. zdta(:,:) <= gdept_1d(jk+1) )   idta(:,:) = jk 
     436!$OMP END WORKSHARE 
    436437                  END DO 
     438!$OMP END PARALLEL 
    437439               ENDIF 
    438440            ELSE 
     
    506508 
    507509         !                                            ! local domain level and meter bathymetries (mbathy,bathy) 
    508 !$OMP PARALLEL WORKSHARE 
     510!$OMP PARALLEL 
     511!$OMP WORKSHARE 
    509512         mbathy(:,:) = 0                                   ! set to zero extra halo points 
    510513         bathy (:,:) = 0._wp                               ! (require for mpp case) 
    511 !$OMP END PARALLEL WORKSHARE 
    512 !$OMP PARALLEL DO schedule(static) private(jj, ji) 
     514!$OMP END WORKSHARE 
     515!$OMP DO schedule(static) private(jj, ji) 
    513516         DO jj = 1, nlcj                                   ! interior values 
    514517            DO ji = 1, nlci 
     
    517520            END DO 
    518521         END DO 
    519 !$OMP PARALLEL WORKSHARE 
     522!$OMP END DO NOWAIT 
     523!$OMP WORKSHARE 
    520524         risfdep(:,:)=0.e0 
    521525         misfdep(:,:)=1 
    522 !$OMP END PARALLEL WORKSHARE 
     526!$OMP END WORKSHARE NOWAIT 
     527!$OMP END PARALLEL 
    523528         ! 
    524529         DEALLOCATE( idta, zdta ) 
     
    532537            CALL iom_get  ( inum, jpdom_data, 'Bathy_level', bathy ) 
    533538            CALL iom_close( inum ) 
    534 !$OMP PARALLEL WORKSHARE 
    535539            mbathy(:,:) = INT( bathy(:,:) ) 
    536540            ! initialisation isf variables 
     541!$OMP PARALLEL WORKSHARE 
    537542            risfdep(:,:)=0._wp ; misfdep(:,:)=1              
    538543!$OMP END PARALLEL WORKSHARE 
     
    11961201         END DO 
    11971202      ELSE ! no cavity 
     1203!$OMP PARALLEL  
     1204!$OMP WORKSHARE  
    11981205         gde3w_0(:,:,1) = 0.5_wp * e3w_0(:,:,1) 
     1206!$OMP END WORKSHARE 
    11991207         DO jk = 2, jpk 
    1200 !$OMP PARALLEL DO schedule(static) private(jj, ji)  
     1208!$OMP DO schedule(static) private(jj, ji)  
    12011209             DO jj =1, jpj  
    12021210                DO ji=1, jpi 
     
    12051213             END DO 
    12061214         END DO 
     1215!$OMP END PARALLEL  
    12071216      END IF 
    12081217      ! 
Note: See TracChangeset for help on using the changeset viewer.