Changeset 2251 for CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/trc_oce.F90

Timestamp:

05/26/14 12:41:58 (10 years ago)

Author:

aclsce

Message:

Modifed IPSLCM6 NEMO sources to fit with NEMO revision 4624

File:

• CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/trc_oce.F90 (modified) (4 diffs)

CONFIG/UNIFORM/v6/IPSLCM6/SOURCES/NEMO/trc_oce.F90

@@ -67 +67 @@
    !!----------------------------------------------------------------------
    !! NEMO/OPA 3.3 , NEMO Consortium (2010)
-   !! $Id: trc_oce.F90 3770 2013-02-04 10:51:28Z cetlod $ 
+   !! $Id: trc_oce.F90 4623 2014-04-22 15:21:41Z acc $ 
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
@@ -178 +178 @@
       zrgb(1,61) = 10.000   ;   zrgb(2,61) = 0.47804   ;   zrgb(3,61) = 0.27178   ;   zrgb(4,61) = 0.56870
       !
+      !
       prgb(:,:) = zrgb(2:4,:)
       !
@@ -270 +271 @@
       !!----------------------------------------------------------------------
       !
-      ! 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) )
       !
@@ -276 +277 @@
       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
       !