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 8279 for branches/2017/wrk_OMP_test_for_Silvia/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90 – NEMO

Ignore:
Timestamp:
2017-07-04T17:46:48+02:00 (7 years ago)
Author:
mocavero
Message:

Implementation of OMP coarse-grained parallelization on ZDF new package

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/wrk_OMP_test_for_Silvia/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfric.F90

    r8056 r8279  
    164164      !                       !==  avm and avt = F(Richardson number)  ==! 
    165165      DO jk = 2, jpkm1 
    166          DO jj = k_Jstr, k_Jend 
    167             DO ji = k_Jstr, k_Iend        ! coefficient = F(richardson number) (avm-weighted Ri) 
     166         DO jj = tnldj, tnlej 
     167            DO ji = tnldi, tnlei        ! coefficient = F(richardson number) (avm-weighted Ri) 
    168168               zcfRi = 1._wp / (  1._wp + rn_alp * MAX(  0._wp , p_avm(ji,jj,jk) * rn2(ji,jj,jk) / ( p_sh2(ji,jj,jk) + 1.e-20 ) )  ) 
    169169               zav   = rn_avmri * zcfRi**nn_ric 
     
    180180      IF( ln_mldw ) THEN      !==  set a minimum value in the Ekman layer  ==! 
    181181         ! 
    182          DO jj = k_Jstr, k_Jend     !* Ekman depth 
    183             DO ji = k_Jstr, k_Iend 
     182         DO jj = tnldj, tnlej     !* Ekman depth 
     183            DO ji = tnldi, tnlei 
    184184               zustar = SQRT( taum(ji,jj) * r1_rau0 ) 
    185185               zed    = rn_ekmfc * zustar / ( ABS( ff_t(ji,jj) ) + rsmall )     ! Ekman depth 
     
    188188         END DO 
    189189         DO jk = 2, jpkm1           !* minimum mixing coeff. within the Ekman layer 
    190             DO jj = k_Jstr, k_Jend 
    191                DO ji = k_Jstr, k_Iend 
     190            DO jj = tnldj, tnlej 
     191               DO ji = tnldi, tnlei 
    192192                  IF( pdept(ji,jj,jk) < zh_ekm(ji,jj) ) THEN 
    193193                     p_avm(ji,jj,jk) = MAX(  p_avm(ji,jj,jk), rn_wvmix  ) * wmask(ji,jj,jk) 
Note: See TracChangeset for help on using the changeset viewer.