Changeset 1180 for trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90

Timestamp:

2008-09-11T18:41:47+02:00 (16 years ago)

Author:

cetlod

Message:

update PISCES modules to couple with the sediment model, see ticket:249

File:

• trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90 (modified) (2 diffs)

trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90

@@ -27 +27 @@
       etot, enano, ediat,       &  !: PAR for phyto, nano and diat 
       emoy                         !: averaged PAR in the mixed layer
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &   !:
-      heup                         !: Depth of the euphotic zone
 
    !! * Module variables
@@ -210 +208 @@
 !     Computation of the euphotic depth
 !     ---------------------------------
-    
-      heup(:,:) = 300.e0
+      ! Euphotic layer bottom level
+      neln(:,:) = 1                                           ! initialisation of EL level
+      heup(:,:) = 300.
 
       DO jk = 2, jpkm1
          DO jj = 1, jpj
-            DO ji = 1, jpi
-               IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )   heup(ji,jj) = fsdepw(ji,jj,jk+1)
-            END DO
-         END DO
-      END DO
-
-      heup(:,:) = MIN( 300., heup(:,:) )
+           DO ji = 1, jpi
+              IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )  THEN
+                 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
+                heup(ji,jj) = fsdepw(ji,jj,jk+1)                 ! Euphotic layer depth
+             ENDIF
+          END DO
+        END DO
+     ENDDO
+ 
+     heup(:,:) = MIN( 300., heup(:,:) )
 
 !    Computation of the mean light over the mixed layer depth