Changeset 15083

Timestamp:

2021-07-05T18:29:04+02:00 (3 years ago)

Author:

techene

Message:

#2605 RGBc debug : ORCA2 tests are now reproducible

File:

• NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90 (modified) (1 diff)

NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/TRA/traqsr.F90

@@ -789 +789 @@
       DO jk = jpkm1, 1, -1
          IF( SUM( tmask(:,:,jk) ) > 0 ) THEN    ! ocean point at that level
-            zdw  = MAXVAL( gdepw_0(:,:,jk+1) * wmask(:,:,jk) )    ! max w-depth at jk+1 level
-            ze3t = MINVAL(   e3t_0(:,:,jk  )                 )    ! minimum e3t at jk   level
-            zhext =  - pL * LOG( zcoef * ze3t )                   ! extinction depth
-            IF( zdw >= zhext )   klev = jk                        ! last T-level reached by Qsr
+            zdw  = MAXVAL( gdepw_0(:,:,jk+1) *       wmask(:,:,jk)         )    ! max w-depth at jk+1 level
+            ze3t = MINVAL(   e3t_0(:,:,jk  ) , mask=(wmask(:,:,jk+1)==1)   )    ! minimum e3t at jk   level
+            zhext =  - pL * LOG( zcoef * ze3t )                                 ! extinction depth
+            IF( zdw >= zhext )   klev = jk                                      ! last T-level reached by Qsr
          ELSE                                   ! only land point at level jk
-            klev = jk                                             ! local domain sea-bed level 
+            klev = jk                                                           ! local domain sea-bed level 
          ENDIF
       END DO