Changeset 4623 for trunk/NEMOGCM/NEMO

Timestamp:

2014-04-22T17:21:41+02:00 (10 years ago)

Author:

acc

Message:

#1305. Fix trc_oce_ext_lev function in trc_oce.F90 so that the correct maximum light penetration level is determined even if all a processing elements wet points are above the maximum penetration depth

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90 (modified) (2 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

@@ -267 +267 @@
       !!----------------------------------------------------------------------
       !
-      ! It is not necessary to compute anything bellow the following depth
+      ! It is not necessary to compute anything below the following depth
       zhext = prldex * ( LOG(10._wp) * zprec + LOG(pqsr_frc) )
       !
@@ -273 +273 @@
       pjl = jpkm1
       DO jk = jpkm1, 1, -1
-         zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) )
-         IF( zem >= zhext )   pjl = jk                       ! last T-level reached by Qsr
+         IF(SUM(tmask(:,:,jk)) > 0 ) THEN
+            zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) )
+            IF( zem >= zhext )   pjl = jk                       ! last T-level reached by Qsr
+         ELSE
+            pjl = jk                                            ! or regional sea-bed depth 
+         ENDIF
       END DO
       !