Changeset 4634 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

Timestamp:

2014-05-12T22:46:18+02:00 (10 years ago)

Author:

clem

Message:

major changes in heat budget

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90 (modified) (1 diff)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/OPA_SRC/trc_oce.F90

@@ -271 +271 @@
       !
       ! Level of light extinction
+!cr      pjl = jpkm1
+!cr      DO jk = jpkm1, 1, -1
+!cr         zem = MAXVAL( fsdepw(:,:,jk+1) * tmask(:,:,jk) )
+!cr         IF( zem >= zhext )   pjl = jk                       ! last T-level reached by Qsr
+!cr      END DO
       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
       !