Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zlim.F90

Timestamp:

2020-07-02T20:34:16+02:00 (4 years ago)

Author:

aumont

Message:

update of the PISCES comments

File:

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

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zlim.F90

@@ -84 +84 @@
       !!
       !! ** Method  : - Limitation follows the Liebieg law of the minimum
+      !!              - Monod approach for N, P and Si. Quota approach 
+      !!                for Iron
       !!---------------------------------------------------------------------
       INTEGER, INTENT(in)  :: kt, knt
@@ -113 +115 @@
                ! Computation of a variable Ks of diatoms taking into account
                ! that increasing biomass is made of generally bigger cells
+               ! The allometric relationship is classical.
                !------------------------------------------------------------
                z1_trbphy   = 1. / ( trb(ji,jj,jk,jpphy) + rtrn )
@@ -156 +159 @@
 
                ! Michaelis-Menten Limitation term by nutrients: Nanophyto
+               ! Optimal parameterization by Smith and Pahlow series of 
+               ! papers is used. Optimal allocation is supposed independant
+               ! for all nutrients. 
                ! --------------------------------------------------------
-               ! Limitation of Fe uptake
+
+               ! Limitation of Fe uptake (Quota formalism)
                zfalim = (1.-fananof) / fananof
                xnanofer(ji,jj,jk) = (1. - fananof) * biron(ji,jj,jk) / ( biron(ji,jj,jk) + zfalim * concnfe(ji,jj,jk) )
@@ -181 +188 @@
                !   Michaelis-Menten Limitation term by nutrients : Diatoms
                !   -------------------------------------------------------
-               ! Limitation of Fe uptake
+               ! Limitation of Fe uptake (Quota formalism)
                zfalim = (1.-fadiatf) / fadiatf
                xdiatfer(ji,jj,jk) = (1. - fadiatf) * biron(ji,jj,jk) / ( biron(ji,jj,jk) + zfalim * concdfe(ji,jj,jk) )
@@ -256 +263 @@
             DO ji = 1, jpi
                ! denitrification factor computed from O2 levels
+               ! This factor diagnoses below which level of O2 denitrification
+               ! is active
                nitrfac(ji,jj,jk) = MAX(  0.e0, 0.4 * ( 6.e-6  - trb(ji,jj,jk,jpoxy) )    &
                   &                                / ( oxymin + trb(ji,jj,jk,jpoxy) )  )
@@ -261 +270 @@
                !
                ! redox factor computed from NO3 levels
+               ! This factor diagnoses below which level of NO3 additional redox
+               ! reactions are taking place.
                nitrfac2(ji,jj,jk) = MAX( 0.e0,       ( 1.E-6 - trb(ji,jj,jk,jpno3) )  &
                   &                                / ( 1.E-6 + trb(ji,jj,jk,jpno3) ) )
@@ -296 +307 @@
       !!----------------------------------------------------------------------
       INTEGER ::   ios   ! Local integer
-      !
+
+      ! Namelist block
       NAMELIST/namp4zlim/ concnno3, concdno3, concnnh4, concdnh4, concnfer, concdfer, concbfe,   &
          &                concbno3, concbnh4, xsizedia, xsizephy, xsizern, xsizerd,          & 
@@ -348 +360 @@
       !!----------------------------------------------------------------------
       !!                     ***  ROUTINE p5z_lim_alloc  ***
+      !! 
+      !            Allocation of the arrays used in this module
       !!----------------------------------------------------------------------
       USE lib_mpp , ONLY: ctl_stop