Changeset 14385 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zbio.F90

Timestamp:

2021-02-03T16:03:34+01:00 (3 years ago)

Author:

cetlod

Message:

dev_r11708_aumont_PISCES_QUOTA : merge with the trunk

File:

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

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

@@ -2 +2 @@
    !!======================================================================
    !!                         ***  MODULE p4zbio  ***
-   !! TOP :   PISCES bio-model
+   !! TOP :   PISCES biogeochemical model
+   !!         This module is for both PISCES and PISCES-QUOTA
    !!======================================================================
    !! History :   1.0  !  2004     (O. Aumont) Original code
    !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
+   !!             3.6  ! 2015 (O. Aumont) PISCES-QUOTA
    !!----------------------------------------------------------------------
    !!   p4z_bio        :   computes the interactions between the different
@@ -53 +55 @@
       !!
       !! ** Purpose :   Ecosystem model in the whole ocean: computes the
-      !!              different interactions between the different compartments
-      !!              of PISCES
+      !!                different interactions between the different compartments
+      !!                of PISCES
       !!
       !! ** Method  : - ???
@@ -66 +68 @@
       !
       IF( ln_timing )   CALL timing_start('p4z_bio')
-      !
-      !     ASSIGN THE SHEAR RATE THAT IS USED FOR AGGREGATION
-      !     OF PHYTOPLANKTON AND DETRITUS
 
+      ! ASSIGN THE SHEAR RATE THAT IS USED FOR AGGREGATION
+      ! OF PHYTOPLANKTON AND DETRITUS. Shear rate is supposed to equal 1
+      ! in the mixed layer and 0.1 below the mixed layer.
       xdiss(:,:,:) = 1.
-!!gm the use of nmld should be better here?
       DO_3D( 1, 1, 1, 1, 2, jpkm1 )
-!!gm  :  use nmln  and test on jk ...  less memory acces
          IF( gdepw(ji,jj,jk+1,Kmm) > hmld(ji,jj) )   xdiss(ji,jj,jk) = 0.01
       END_3D
@@ -81 +81 @@
       CALL p4z_fechem  ( kt, knt, Kbb, Kmm, Krhs )     ! Iron chemistry/scavenging
       !
-      IF( ln_p4z ) THEN
+      IF( ln_p4z ) THEN  ! PISCES standard
+         ! Phytoplankton only sources/sinks terms
          CALL p4z_lim  ( kt, knt, Kbb, Kmm       )     ! co-limitations by the various nutrients
          CALL p4z_prod ( kt, knt, Kbb, Kmm, Krhs )     ! phytoplankton growth rate over the global ocean. 
          !                                          ! (for each element : C, Si, Fe, Chl )
          CALL p4z_mort ( kt,      Kbb,      Krhs )     ! phytoplankton mortality
-         !                                          ! zooplankton sources/sinks routines 
+         ! zooplankton sources/sinks routines 
          CALL p4z_micro( kt, knt, Kbb,      Krhs )     ! microzooplankton
-         CALL p4z_meso ( kt, knt, Kbb,      Krhs )     ! mesozooplankton
-      ELSE
+         CALL p4z_meso ( kt, knt, Kbb, Kmm, Krhs )     ! mesozooplankton
+      ELSE  ! PISCES-QUOTA
+         ! Phytoplankton only sources/sinks terms
          CALL p5z_lim  ( kt, knt, Kbb, Kmm       )     ! co-limitations by the various nutrients
          CALL p5z_prod ( kt, knt, Kbb, Kmm, Krhs )     ! phytoplankton growth rate over the global ocean. 
          !                                          ! (for each element : C, Si, Fe, Chl )
          CALL p5z_mort ( kt,      Kbb,      Krhs      )     ! phytoplankton mortality
-         !                                          ! zooplankton sources/sinks routines 
+         !  zooplankton sources/sinks routines 
          CALL p5z_micro( kt, knt, Kbb,      Krhs )           ! microzooplankton
-         CALL p5z_meso ( kt, knt, Kbb,      Krhs )           ! mesozooplankton
+         CALL p5z_meso ( kt, knt, Kbb, Kmm, Krhs )           ! mesozooplankton
       ENDIF
       !
@@ -103 +105 @@
       CALL p4z_poc     ( kt, knt, Kbb, Kmm, Krhs )     ! Remineralization of organic particles
       !
+      ! Ligand production. ln_ligand should be set .true. to activate
       IF( ln_ligand )  &
       & CALL p4z_ligand( kt, knt, Kbb,      Krhs )
+
+      ! Update of the size of the different phytoplankton groups
+      sized(:,:,:) = MAX(1.0, sizeda(:,:,:) )
+      sizen(:,:,:) = MAX(1.0, sizena(:,:,:) )
+      IF (ln_p5z) THEN
+         sizep(:,:,:) = MAX(1.0, sizepa(:,:,:) )
+      ENDIF
       !                                                             !
       IF(sn_cfctl%l_prttrc)   THEN  ! print mean trends (used for debugging)