Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zmicro.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/p4zmicro.F90 (modified) (7 diffs)

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

@@ -42 +42 @@
    REAL(wp), PUBLIC ::   epsher      !: growth efficiency for grazing 1 
    REAL(wp), PUBLIC ::   epshermin   !: minimum growth efficiency for grazing 1
+   REAL(wp), PUBLIC ::   xsigma      !: Width of the grazing window
+   REAL(wp), PUBLIC ::   xsigmadel   !: Maximum additional width of the grazing window at low food density 
+
 
    !!----------------------------------------------------------------------
@@ -59 +62 @@
       !!                parameterization.
       !!                All living compartments smaller than microzooplankton
-      !!                are potential preys of mesozooplankton
+      !!                are potential preys of microzooplankton
       !!
       !! ** Method  : - ???
@@ -146 +149 @@
                ! ----------------------------------------------------------
                zsigma = 1.0 - zdenom**2/(0.05**2+zdenom**2)
-               zsigma = 0.5 + 1.0*zsigma
+               zsigma = xsigma + xsigmadel * zsigma
                zdiffdn = exp( -ABS(log(1.67 * sizen(ji,jj,jk) / (5.0 * sized(ji,jj,jk) + rtrn )) )**2 / zsigma**2)
                ztmp1 = xprefn * zcompaph * ( zcompaph + zdiffdn * zcompadi ) / ( 1.0 + zdiffdn )
@@ -156 +159 @@
                ztmp3 = ztmp3 / ztmptot
 
-               zgraznc   = zgraze   * ztmp1 * zdenom
-               zgrazdc   = zgraze   * ztmp2 * zdenom
-               zgrazpoc  = zgraze   * ztmp3 * zdenom
-
+               ! Ingestion terms on the different preys of microzooplankton
+               zgraznc   = zgraze   * ztmp1 * zdenom  ! Nanophytoplankton
+               zgrazdc   = zgraze   * ztmp2 * zdenom  ! Diatoms
+               zgrazpoc  = zgraze   * ztmp3 * zdenom  ! POC
+
+               ! Ingestion terms on the iron content of the different preys
                zgraznf   = zgraznc  * trb(ji,jj,jk,jpnfe) / (trb(ji,jj,jk,jpphy) + rtrn)
                zgrazpof  = zgrazpoc * trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
                zgrazdf   = zgrazdc  * trb(ji,jj,jk,jpdfe) / (trb(ji,jj,jk,jpdia) + rtrn)
                !
+               ! Total ingestion rate in C, Fe, N units
                zgraztotc = zgraznc + zgrazpoc + zgrazdc
                zgraztotf = zgraznf + zgrazdf  + zgrazpof 
@@ -183 +189 @@
                zepshert  =  MIN( 1., zgrasratn, zgrasratf / ferat3)
                zbeta     = MAX(0., (epsher - epshermin) )
+               ! Food quantity deprivation of the GGE
                zepsherf  = epshermin + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta )
+               ! Food quality deprivation of the GGE
                zepsherq  = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 )
+               ! Actual GGE of microzooplankton
                zepsherv  = zepsherf * zepshert * zepsherq
+               ! Excretion of Fe
                zgrafer   = zgraztotc * MAX( 0. , ( 1. - unass ) * zgrasratf - ferat3 * zepsherv ) 
+               ! Excretion of C, N, P
                zgrarem   = zgraztotc * ( 1. - zepsherv - unass )
+               ! Egestion of C, N, P
                zgrapoc   = zgraztotc * unass
 
                !  Update of the TRA arrays
                !  ------------------------
+               ! Fraction of excretion as inorganic nutrients and DIC
                zgrarsig  = zgrarem * sigma1
                tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) + zgrarsig
@@ -287 +300 @@
       NAMELIST/namp4zzoo/ part, grazrat, resrat, mzrat, xprefn, xprefc, &
          &                xprefd,  xthreshdia,  xthreshphy,  xthreshpoc, &
-         &                xthresh, xkgraz, epsher, epshermin, sigma1, unass
+         &                xthresh, xkgraz, epsher, epshermin, sigma1, unass  &
+         &                xsigma, xsigmadel
       !!----------------------------------------------------------------------
       !
@@ -321 +335 @@
          WRITE(numout,*) '      Minimum efficicency of microzoo growth          epshermin   =', epshermin
          WRITE(numout,*) '      Fraction of microzoo excretion as DOM           sigma1      =', sigma1
-         WRITE(numout,*) '      half sturation constant for grazing 1           xkgraz      =', xkgraz
+         WRITE(numout,*) '      half saturation constant for grazing 1          xkgraz      =', xkgraz
+         WRITE(numout,*) '      Width of the grazing window                     xsigma      =', xsigma
+         WRITE(numout,*) '      Maximum additional width of the grazing window  xsigmadel   =', xsigmadel
+
       ENDIF
       !