Changeset 1176 for trunk/NEMO/TOP_SRC/LOBSTER/trcopt.F90

Timestamp:

2008-09-11T18:26:59+02:00 (16 years ago)

Author:

cetlod

Message:

update LOBSTER modules to take into account new trends organization, see ticket:248

File:

• trunk/NEMO/TOP_SRC/LOBSTER/trcopt.F90 (modified) (2 diffs)

trunk/NEMO/TOP_SRC/LOBSTER/trcopt.F90

@@ -48 +48 @@
       INTEGER, INTENT( in ) ::   kt   ! index of the time stepping
       INTEGER  ::   ji, jj, jk
-      INTEGER , DIMENSION(jpi,jpj)     ::   itabe           ! euphotic layer last k index
-      INTEGER , DIMENSION(jpi,jpj,jpk) ::   imaske          ! euphotic layer mask
       REAL(wp) ::   zpig                                    ! total pigment
       REAL(wp) ::   zkr                                     ! total absorption coefficient in red
@@ -103 +101 @@
       END DO
 
-      ! determination of euphotic layer depth (xze)
-      ! -------------------------------------
+      ! 3. Determination of euphotic layer depth
+      ! ----------------------------------------
 
-      DO jk = 1, jpk                   ! imaske equal 1 in the euphotic layer, and 0 without
+      ! Euphotic layer bottom level
+      neln(:,:) = 1                                           ! initialisation of EL level
+      DO jk = 1, jpk
          DO jj = 1, jpj
-            DO ji = 1,jpi
-               IF( xpar(ji,jj,jk) >= zpar100(ji,jj) ) THEN
-                  imaske(ji,jj,jk) = 1
-               ELSE
-                  imaske(ji,jj,jk) = 0
-               ENDIF
-            END DO
+           DO ji = 1, jpi
+              IF( xpar(ji,jj,jk) >= zpar100(ji,jj) )   neln(ji,jj) = jk+1 ! 1rst T-level strictly below EL bottom
+              !                                                  ! nb. this is to ensure compatibility with
+              !                                                  ! nmld_trc definition in trd_mld_trc_zint
+           END DO
          END DO
-      END DO
-      !                                ! sum of imaske Cover the vertical with a minimim value of 1 
-      itabe(:,:) = 1                   !   surface value setto 1 <=> set a ninimum value to 1
-      DO jk = 2, jpk
-         DO jj = 1, jpj
-            DO ji = 1,jpi
-               itabe(ji,jj) = itabe(ji,jj) + imaske(ji,jj,jk)
-            END DO
+      ENDDO
+
+      ! Euphotic layer depth
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            xze(ji,jj) = fsdepw( ji, jj, neln(ji,jj) )            ! exact EL depth
          END DO
-      END DO
-      DO jj = 1, jpj                  ! converte the number of level into depth
-         DO ji = 1,jpi
-            xze(ji,jj) = fsdepw(ji,jj,itabe(ji,jj)+1)
-         END DO
-      END DO
-      !
+      ENDDO 
+
 
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)