Changeset 12537 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zbio.F90

Timestamp:

2020-03-11T16:02:54+01:00 (4 years ago)

Author:

aumont

Message:

Comments in routines have been revised and significantly augmented

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zbio.F90 (modified) (6 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/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
@@ -50 +52 @@
       !!
       !! ** 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  : - ???
@@ -62 +64 @@
       !
       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
       xdiss(:,:,:) = 1.
-!!gm the use of nmld should be better here?
       DO jk = 2, jpkm1
          DO jj = 1, jpj
             DO ji = 1, jpi
-!!gm  :  use nmln  and test on jk ...  less memory acces
                IF( gdepw_n(ji,jj,jk+1) > hmld(ji,jj) )   xdiss(ji,jj,jk) = 0.01
             END DO 
@@ -77 +76 @@
       END DO
 
-      CALL p4z_opt     ( kt, knt )     ! Optic: PAR in the water column
+      CALL p4z_opt     ( kt, knt )     ! Optics: PAR in the water column
       CALL p4z_sink    ( kt, knt )     ! vertical flux of particulate organic matter
       CALL p4z_fechem  ( kt, knt )     ! Iron chemistry/scavenging
@@ -87 +86 @@
          CALL p4z_mort ( kt      )     ! phytoplankton mortality
          !                             ! zooplankton sources/sinks routines 
-         CALL p4z_micro( kt, knt )           ! microzooplankton
-         CALL p4z_meso ( kt, knt )           ! mesozooplankton
+         CALL p4z_micro( kt, knt )     ! microzooplankton
+         CALL p4z_meso ( kt, knt )     ! mesozooplankton
       ELSE
          CALL p5z_lim  ( kt, knt )     ! co-limitations by the various nutrients
          CALL p5z_prod ( kt, knt )     ! phytoplankton growth rate over the global ocean. 
-         !                             ! (for each element : C, Si, Fe, Chl )
+         !                             ! (for each element : C, N, P, Si, Fe, Chl )
          CALL p5z_mort ( kt      )     ! phytoplankton mortality
          !                             ! zooplankton sources/sinks routines 
-         CALL p5z_micro( kt, knt )           ! microzooplankton
-         CALL p5z_meso ( kt, knt )           ! mesozooplankton
+         CALL p5z_micro( kt, knt )     ! microzooplankton
+         CALL p5z_meso ( kt, knt )     ! mesozooplankton
       ENDIF
       !
@@ -103 +102 @@
       CALL p4z_poc     ( kt, knt )     ! Remineralization of organic particles
       !
+      ! Ligand production. ln_ligand should be set .true. to activate
       IF( ln_ligand )  &
       & CALL p4z_ligand( kt, knt )
-      ! Update of the size of the organisms
+
+      ! Update of the size of the different phytoplankton groups
       IF (ln_p5z) THEN
          sized(:,:,:) = sizeda(:,:,:)