Changeset 5288 for branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90

Timestamp:

2015-05-22T18:22:18+02:00 (9 years ago)

Author:

aumont

Message:

various bug fixes and updates of PISCES quota

File:

• branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90 (modified) (3 diffs)

branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90

@@ -8 +8 @@
    !!             3.2  !  2009-04  (C. Ethe, G. Madec)  optimisation
    !!             3.4  !  2011-06  (O. Aumont, C. Ethe) Improve light availability of nano & diat
+   !!             3.6  !  2015-05  (O. Aumont) PISCES quota
    !!----------------------------------------------------------------------
 #if defined  key_pisces || defined key_pisces_quota
@@ -218 +219 @@
       neln(:,:) = 1                            !  ------------------------
       heup(:,:) = 300.
+      heup_01(:,:) = 300.
 
       DO jk = 2, nksrp
          DO jj = 1, jpj
            DO ji = 1, jpi
-              IF( etot(ji,jj,jk) * tmask(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )  THEN
+              IF( etot(ji,jj,jk) * tmask(ji,jj,jk) >= 0.01 * parlux * qsr(ji,jj) )  THEN
                  neln(ji,jj) = jk+1                    ! Euphotic level : 1rst T-level strictly below Euphotic layer
                  !                                     ! nb: ensure the compatibility with nmld_trc definition in trd_mld_trc_zint
                  heup(ji,jj) = fsdepw(ji,jj,jk+1)      ! Euphotic layer depth
               ENDIF
+              IF( etot(ji,jj,jk) * tmask(ji,jj,jk) >= 0.1 )  THEN
+                 heup_01(ji,jj) = fsdepw(ji,jj,jk+1)      ! Absolute euphotic layer depth
+              ENDIF
            END DO
         END DO
       END DO
  
-      heup(:,:) = MIN( 300., heup(:,:) )
+      heup(:,:)    = MIN( 300., heup(:,:) )
+      heup_01(:,:) = MIN( 300., heup_01(:,:) )
 
       !                                        !* mean light over the mixed layer
@@ -284 +290 @@
         IF( lk_iomput ) THEN
            IF( jnt == nrdttrc ) THEN
-              CALL iom_put( "Heup", heup(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht
+              CALL iom_put( "Heup", heup_01(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht
               CALL iom_put( "PAR" , emoy(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
            ENDIF