Changeset 14786 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zprod.F90

Timestamp:

2021-05-05T09:09:10+02:00 (3 years ago)

Author:

aumont

Message:

Various bug fixes in PISCES

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zprod.F90 (modified) (5 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zprod.F90

@@ -100 +100 @@
       ! of the thermal dependency
       ! Parameters are taken from Bissinger et al. (2008)
-      zprmaxn(:,:,:) = 0.8_wp * r1_rday * tgfunc(:,:,:)
+      zprmaxn(:,:,:) = 0.65_wp * r1_rday * tgfunc(:,:,:)
       zprmaxd(:,:,:) = zprmaxn(:,:,:)
 
@@ -109 +109 @@
       ! absolute light level definition of the euphotic zone
       ! ------------------------------------------------------------------------- 
-      DO_3D( 1, 1, 1, 1, 1, jpkm1 )
-         IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN
-            zval = MAX( 1., strn(ji,jj) )
-            IF( gdepw(ji,jj,jk+1,Kmm) <= hmld(ji,jj) ) THEN
-               zval = zval * MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn ))
+      IF ( ln_p4z_dcyc ) THEN    ! Diurnal cycle in PISCES
+
+         DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+            IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN
+               IF( gdepw(ji,jj,jk+1,Kmm) <= hmld(ji,jj) ) THEN
+                  zval = MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn ))
+               ENDIF
+               zmxl_chl(ji,jj,jk) = zval / 24.
+               zmxl_fac(ji,jj,jk) = 1.0 - exp( -0.26 * zval )
             ENDIF
-            zmxl_chl(ji,jj,jk) = zval / 24.
-            zmxl_fac(ji,jj,jk) = 1.0 - exp( -0.26 * zval )
-         ENDIF
-      END_3D
+         END_3D
+ 
+      ELSE ! No diurnal cycle in PISCES
+
+         DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+            IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN
+               zval = MAX( 1., strn(ji,jj) )
+               IF( gdepw(ji,jj,jk+1,Kmm) <= hmld(ji,jj) ) THEN
+                  zval = zval * MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn ))
+               ENDIF
+               zmxl_chl(ji,jj,jk) = zval / 24.
+               zmxl_fac(ji,jj,jk) = 1.0 - exp( -0.26 * zval )
+            ENDIF
+         END_3D
+
+      ENDIF
 
       zprbio(:,:,:) = zprmaxn(:,:,:) * zmxl_fac(:,:,:)
@@ -179 +195 @@
           zval = MIN( xnanopo4(ji,jj,jk), ( xnanonh4(ji,jj,jk) + xnanono3(ji,jj,jk) ) )   &
           &      * zprmaxn(ji,jj,jk) / ( zprbio(ji,jj,jk) + rtrn )
-          quotan(ji,jj,jk) = MIN( 1., 0.2 + 0.8 * zval )
+          quotan(ji,jj,jk) = MIN( 1., 0.3 + 0.7 * zval )
           zval = MIN( xdiatpo4(ji,jj,jk), ( xdiatnh4(ji,jj,jk) + xdiatno3(ji,jj,jk) ) )   &
           &      * zprmaxd(ji,jj,jk) / ( zprdia(ji,jj,jk) + rtrn )
-          quotad(ji,jj,jk) = MIN( 1., 0.2 + 0.8 * zval )
+          quotad(ji,jj,jk) = MIN( 1., 0.3 + 0.7 * zval )
       END_3D
 
@@ -256 +272 @@
             &          + xnanonh4(ji,jj,jk) ) * (1. - xnanofer(ji,jj,jk) ) )   &
             &          * xnanofer(ji,jj,jk) * zmax * tr(ji,jj,jk,jpphy,Kbb) * rfact2
-
             ! production terms of diatoms (C)
             zprorcad(ji,jj,jk) = zprdia(ji,jj,jk) * xlimdia(ji,jj,jk) * tr(ji,jj,jk,jpdia,Kbb) * rfact2
@@ -337 +352 @@
            consfe3(ji,jj,jk)   = zprodfer * 75.0 / ( rtrn + ( plig(ji,jj,jk) + 75.0 * (1.0 - plig(ji,jj,jk) ) )   &
            &                   * tr(ji,jj,jk,jpfer,Kbb) ) / rfact2
-           
            !
            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zpptot