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 7753 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90 – NEMO

Ignore:
Timestamp:
2017-03-03T12:46:59+01:00 (7 years ago)
Author:
mocavero
Message:

Reverting trunk to remove OpenMP

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

    r7698 r7753  
    7171      ! Constrain the thickness of the smallest category above himin 
    7272      !---------------------------------------------------------------------- 
    73 !$OMP PARALLEL 
    74 !$OMP DO schedule(static) private(jj,ji,rswitch) 
    7573      DO jj = 1, jpj  
    7674         DO ji = 1, jpi 
     
    8785      ! ice concentration should not exceed amax  
    8886      !----------------------------------------------------- 
    89 !$OMP DO schedule(static) private(jj, ji) 
    90       DO jj = 1, jpj 
    91          DO ji = 1, jpi 
    92             at_i(ji,jj) = 0._wp 
    93          END DO 
    94       END DO 
     87      at_i(:,:) = 0._wp 
    9588      DO jl = 1, jpl 
    96 !$OMP DO schedule(static) private(jj, ji) 
    97          DO jj = 1, jpj 
    98             DO ji = 1, jpi 
    99                at_i(ji,jj) = a_i(ji,jj,jl) + at_i(ji,jj) 
    100             END DO 
    101          END DO 
     89         at_i(:,:) = a_i(:,:,jl) + at_i(:,:) 
    10290      END DO 
    10391 
    10492      DO jl  = 1, jpl 
    105 !$OMP DO schedule(static) private(jj, ji) 
    10693         DO jj = 1, jpj 
    10794            DO ji = 1, jpi 
     
    113100         END DO 
    114101      END DO 
    115 !$OMP END PARALLEL 
    116102 
    117103      !--------------------- 
     
    120106      IF (  nn_icesal == 2  ) THEN  
    121107         DO jl = 1, jpl 
    122 !$OMP PARALLEL DO schedule(static) private(jj,ji,zsal,rswitch) 
    123108            DO jj = 1, jpj  
    124109               DO ji = 1, jpi 
     
    149134      ! Ice drift 
    150135      !------------ 
    151 !$OMP PARALLEL DO schedule(static) private(jj, ji) 
    152136      DO jj = 2, jpjm1 
    153137         DO ji = 2, jpim1 
     
    164148      CALL lbc_lnk( v_ice(:,:), 'V', -1. ) 
    165149      !mask velocities 
    166 !$OMP PARALLEL 
    167 !$OMP DO schedule(static) private(jj, ji) 
    168       DO jj = 1, jpj 
    169          DO ji = 1, jpi 
    170             u_ice(ji,jj) = u_ice(ji,jj) * umask(ji,jj,1) 
    171             v_ice(ji,jj) = v_ice(ji,jj) * vmask(ji,jj,1) 
    172          END DO 
    173       END DO 
     150      u_ice(:,:) = u_ice(:,:) * umask(:,:,1) 
     151      v_ice(:,:) = v_ice(:,:) * vmask(:,:,1) 
    174152  
    175153      ! ------------------------------------------------- 
     
    177155      ! ------------------------------------------------- 
    178156      DO jl  = 1, jpl 
    179 !$OMP DO schedule(static) private(jj, ji) 
    180          DO jj = 1, jpj 
    181             DO ji = 1, jpi 
    182                oa_i(ji,jj,jl) = oa_i(ji,jj,jl) + a_i(ji,jj,jl) * rdt_ice / rday   ! ice natural aging 
    183                afx_thd(ji,jj) = afx_thd(ji,jj) + ( a_i(ji,jj,jl) - a_i_b(ji,jj,jl) ) * r1_rdtice 
    184             END DO 
    185          END DO 
     157         oa_i(:,:,jl) = oa_i(:,:,jl) + a_i(:,:,jl) * rdt_ice / rday   ! ice natural aging 
     158         afx_thd(:,:) = afx_thd(:,:) + ( a_i(:,:,jl) - a_i_b(:,:,jl) ) * r1_rdtice 
    186159      END DO 
    187160      afx_tot = afx_thd + afx_dyn 
    188161 
    189 !$OMP DO schedule(static) private(jj, ji) 
    190162      DO jj = 1, jpj 
    191163         DO ji = 1, jpi             
     
    201173         END DO 
    202174      END DO 
    203 !$OMP END PARALLEL 
    204175 
    205176      ! conservation test 
Note: See TracChangeset for help on using the changeset viewer.