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/LIM_SRC_2/limthd_2.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/LIM_SRC_2/limthd_2.F90

    r6140 r7037  
    121121       
    122122!!gm needed?  yes at least for some of these arrays  
     123!$OMP PARALLEL 
     124!$OMP WORKSHARE 
    123125      zdvosif(:,:) = 0.e0   ! variation of ice volume at surface 
    124126      zdvobif(:,:) = 0.e0   ! variation of ice volume at bottom 
     
    137139      rdq_ice(:,:) = 0.e0   ! heat content associated with rdm_ice 
    138140      zmsk (:,:,:) = 0.e0 
     141!$OMP END WORKSHARE NOWAIT 
    139142 
    140143      ! set to zero snow thickness smaller than epsi04 
     144!$OMP DO schedule(static) private(jj,ji) 
    141145      DO jj = 1, jpj 
    142146         DO ji = 1, jpi 
     
    144148         END DO 
    145149      END DO 
     150!$OMP END DO NOWAIT 
     151!$OMP END PARALLEL 
    146152!!gm better coded (do not use SIGN...) 
    147153!     WHERE( hsnif(:,:) < epsi04 )   hsnif(:,:) = 0.e0 
     
    154160      !-----------------------------------! 
    155161       
     162!$OMP PARALLEL DO schedule(static) private(jj,ji,zindg,zthsnice,zindb,za,zh) 
    156163      DO jj = 1, jpj 
    157164         DO ji = 1, jpi 
     
    196203      !-------------------------------------------------------------------------- 
    197204 
     205!$OMP PARALLEL DO schedule(static) private(jj,ji,zthsnice,zindb,zinda,zfricp,zfric_u,zfntlat,zpareff) 
    198206      DO jj = 1, jpj 
    199207         DO ji = 1, jpi 
     
    337345      ! Up-date sea ice thickness 
    338346      !-------------------------- 
     347!$OMP PARALLEL DO schedule(static) private(jj,ji) 
    339348      DO jj = 1, jpj 
    340349         DO ji = 1, jpi 
     
    348357      !-------------------------------------------------------- 
    349358      IF( fcor(1,1) < 0.e0 ) THEN 
     359!$OMP PARALLEL DO schedule(static) private(jj,ji) 
    350360         DO jj = 1, njeqm1 
    351361            DO ji = 1, jpi 
     
    418428      ! Update daily thermodynamic ice production.     
    419429      !------------------------------------------------------------------------------ 
     430!$OMP PARALLEL DO schedule(static) private(jj,ji) 
    420431      DO jj = 1, jpj 
    421432         DO ji = 1, jpi 
     
    453464      IF( iom_use('u_imasstr') ) THEN 
    454465         zztmp = 0.25 * rhoic 
     466!$OMP PARALLEL DO schedule(static) private(jj,ji,zuice_m,zhice_u) 
    455467         DO jj = 1, jpjm1  
    456468            DO ji = 1, jpim1   ! NO vector opt. 
     
    466478      IF( iom_use('v_imasstr') ) THEN 
    467479         zztmp = 0.25 * rhoic 
     480!$OMP PARALLEL DO schedule(static) private(jj,ji,zvice_m,zhice_v) 
    468481         DO jj = 1, jpjm1  
    469482            DO ji = 1, jpim1   ! NO vector opt. 
     
    480493      !! Fram Strait sea-ice transport (sea-ice + snow)  (in ORCA2 = 5 points) 
    481494      IF( iom_use('fram_trans') .and. cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN    ! ORCA R2 configuration 
     495!$OMP PARALLEL DO schedule(static) private(jj,ji,zrhoij,zrhoijm1,ztr_fram) 
    482496         DO jj = mj0(137), mj1(137) ! B grid 
    483497            IF( mj0(jj-1) >= nldj ) THEN 
     
    497511!! ce     A big warning because the model crashes on IDRIS/IBM SP6 with xlf 13.1.0.3, see ticket #761 
    498512!! ce     We Unroll the loop and everything works fine       
     513!$OMP PARALLEL DO schedule(static) private(jj,ji) 
    499514         DO jj = 1, jpj 
    500515            DO ji = 1, jpi 
Note: See TracChangeset for help on using the changeset viewer.