Ignore:
Timestamp:
2018-05-22T11:09:09+02:00 (2 years ago)
Author:
andmirek
Message:

#2001 few additionale changes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package_OMP/NEMOGCM/NEMO/OPA_SRC/TRA/traqsr.F90

    r6498 r9616  
    153153      IF( lk_qsr_bio .AND. ln_qsr_bio ) THEN      !  bio-model fluxes  : all vertical coordinates  ! 
    154154         !                                        ! ============================================== ! 
     155!$OMP PARALLEL DO  
    155156         DO jk = 1, jpkm1 
    156157            qsr_hc(:,:,jk) = r1_rau0_rcp * ( etot3(:,:,jk) - etot3(:,:,jk+1) ) 
    157158         END DO 
    158159         !                                        Add to the general trend 
     160!$OMP PARALLEL DO PRIVATE (z1_e3t) 
    159161         DO jk = 1, jpkm1 
    160162            DO jj = 2, jpjm1  
     
    198200                  ! 
    199201                  CALL fld_read( kt, 1, sf_chl )            ! Read Chl data and provides it at the current time step 
    200 !CDIR NOVERRCHK   ! 
     202!$OMP PARALLEL DO PRIVATE(zchl, zCtot, zze, zlogc, zlogc2, zlogc3, zCb, zCmax, & 
     203!$OMP&                    zpsimax, zdelpsi, zCze, zpsi, jk) SHARED (nksr) 
    201204                  DO jj = 1, jpj 
    202 !CDIR NOVERRCHK 
    203205                     DO ji = 1, jpi 
    204206                        zchl    = sf_chl(1)%fnow(ji,jj,1) 
     
    264266               END DO 
    265267               ! 
     268!$OMP PARALLEL DO 
    266269               DO jk = 1, nksr                                        ! compute and add qsr trend to ta 
    267270                  qsr_hc(:,:,jk) = r1_rau0_rcp * ( zea(:,:,jk) - zea(:,:,jk+1) ) 
     
    271274               ! 
    272275               IF ( ln_qsr_ice ) THEN    ! store attenuation coefficient of the first ocean level 
    273 !CDIR NOVERRCHK 
    274276                  DO jj = 1, jpj                                         ! Separation in R-G-B depending of the surface Chl 
    275 !CDIR NOVERRCHK 
    276277                     DO ji = 1, jpi 
    277278                        zchl = MIN( 10. , MAX( 0.03, zchl3d(ji,jj,1) ) ) 
     
    283284                  END DO 
    284285                  !  
     286!$OMP PARALLEL DO PRIVATE(zc0, zc1, zc2, zc3) 
    285287                  DO jj = 1, jpj 
    286288                     DO ji = 1, jpi 
     
    313315               zz0   =        rn_abs   * r1_rau0_rcp 
    314316               zz1   = ( 1. - rn_abs ) * r1_rau0_rcp 
     317!$OMP PARALLEL DO PRIVATE (zc0, zc1) 
    315318               DO jk = 1, nksr                    ! solar heat absorbed at T-point in the top 400m  
    316319                  DO jj = 1, jpj 
     
    324327               ! clem: store attenuation coefficient of the first ocean level 
    325328               IF ( ln_qsr_ice ) THEN 
     329!$OMP PARALLEL DO PRIVATE (zc0, zc1) 
    326330                  DO jj = 1, jpj 
    327331                     DO ji = 1, jpi 
     
    333337               ENDIF 
    334338            ELSE                                               !* constant volume: coef. computed one for all 
     339!$OMP PARALLEL DO 
    335340               DO jk = 1, nksr 
    336341                  DO jj = 2, jpjm1 
     
    351356         ! 
    352357         !                                        Add to the general trend 
     358!$OMP PARALLEL DO PRIVATE (z1_e3t) 
    353359         DO jk = 1, nksr 
    354360            DO jj = 2, jpjm1  
     
    538544                  ze3(:,:,1) = zcoef 
    539545                  zea(:,:,1) = tmask(:,:,1)                   ! = ( ze0+ze1+z2+ze3 ) * tmask 
    540                 
     546!$OMP PARALLEL DO PRIVATE(zc0, zc1, zc2, zc3) 
    541547                  DO jk = 2, nksr+1 
    542 !CDIR NOVERRCHK 
    543548                     DO jj = 1, jpj 
    544 !CDIR NOVERRCHK    
    545549                        DO ji = 1, jpi 
    546550                           zc0 = ze0(ji,jj,jk-1) * EXP( - e3t_0(ji,jj,jk-1) * xsi0r     ) 
     
    557561                  END DO  
    558562                  ! 
     563!$OMP PARALLEL DO 
    559564                  DO jk = 1, nksr 
    560565                     ! (ISF) no light penetration below the ice shelves 
     
    582587               zz0 =        rn_abs   * r1_rau0_rcp 
    583588               zz1 = ( 1. - rn_abs ) * r1_rau0_rcp 
     589!$OMP PARALLEL DO PRIVATE(zc0, zc1) 
    584590               DO jk = 1, nksr                    !*  solar heat absorbed at T-point computed once for all 
    585591                  DO jj = 1, jpj                              ! top 400 meters 
Note: See TracChangeset for help on using the changeset viewer.