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

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

@@ -15 +15 @@
    USE trc             !  passive tracers common variables 
    USE sms_pisces      !  PISCES Source Minus Sink variables
-   USE p4zlim
+   USE p4zlim          !  PISCES nutrient limitation term of PISCES std
    USE p5zlim          !  Phytoplankton limitation terms
    USE iom             !  I/O manager
@@ -50 +50 @@
    REAL(wp), PUBLIC ::  srespir     !: half sturation constant for grazing 1 
    REAL(wp), PUBLIC ::  ssigma      !: Fraction excreted as semi-labile DOM
+   REAL(wp), PUBLIC ::  xsigma      !: Width of the grazing window
+   REAL(wp), PUBLIC ::  xsigmadel   !: Maximum additional width of the grazing window at low food density
    LOGICAL,  PUBLIC ::  bmetexc     !: Use of excess carbon for respiration
 
@@ -161 +163 @@
                ! ----------------------------------------------------------
                zsigma = 1.0 - zdenom**2/(0.05**2+zdenom**2)
-               zsigma = 0.5 + 1.0*zsigma
+               zsigma = xsigma + xsigmadel * zsigma
                zdiffpn = exp( -ABS(log(0.5 * sizep(ji,jj,jk) / (3.0 * sizen(ji,jj,jk) + rtrn )) )**2 / zsigma**2 )
                zdiffdn = exp( -ABS(log(3.0 * sizen(ji,jj,jk) / (5.0 * sized(ji,jj,jk) + rtrn )) )**2 / zsigma**2)
@@ -179 +181 @@
                !   Microzooplankton regular grazing on the different preys
                !   -------------------------------------------------------
-               zgraznc   = zgraze  * ztmp1  * zdenom
+               !   Nanophytoplankton
+               zgraznc   = zgraze  * ztmp1  * zdenom  
                zgraznn   = zgraznc * trb(ji,jj,jk,jpnph) / (trb(ji,jj,jk,jpphy) + rtrn)
                zgraznp   = zgraznc * trb(ji,jj,jk,jppph) / (trb(ji,jj,jk,jpphy) + rtrn)
                zgraznf   = zgraznc * trb(ji,jj,jk,jpnfe) / (trb(ji,jj,jk,jpphy) + rtrn)
+
+               ! Picophytoplankton
                zgrazpc   = zgraze  * ztmp2  * zdenom
                zgrazpn   = zgrazpc * trb(ji,jj,jk,jpnpi) / (trb(ji,jj,jk,jppic) + rtrn)
                zgrazpp   = zgrazpc * trb(ji,jj,jk,jpppi) / (trb(ji,jj,jk,jppic) + rtrn)
                zgrazpf   = zgrazpc * trb(ji,jj,jk,jppfe) / (trb(ji,jj,jk,jppic) + rtrn)
+               ! Microzooplankton
                zgrazz    = zgraze  * ztmp5   * zdenom
+
+               ! small POC
                zgrazpoc  = zgraze  * ztmp3   * zdenom
                zgrazpon  = zgrazpoc * trb(ji,jj,jk,jppon) / ( trb(ji,jj,jk,jppoc) + rtrn )
                zgrazpop  = zgrazpoc * trb(ji,jj,jk,jppop) / ( trb(ji,jj,jk,jppoc) + rtrn )
                zgrazpof  = zgrazpoc* trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
+
+               ! Diatoms
                zgrazdc   = zgraze  * ztmp4  * zdenom
                zgrazdn   = zgrazdc * trb(ji,jj,jk,jpndi) / (trb(ji,jj,jk,jpdia) + rtrn)
@@ -197 +207 @@
                zgrazdf   = zgrazdc * trb(ji,jj,jk,jpdfe) / (trb(ji,jj,jk,jpdia) + rtrn)
                !
+               ! Total ingestion rates in C, P, Fe, N
                zgraztotc = zgraznc + zgrazpoc + zgrazdc + zgrazz + zgrazpc
                zgraztotn = zgraznn + zgrazpn + zgrazpon + zgrazdn + zgrazz * no3rat3
@@ -221 +232 @@
                zepshert  = MIN( 1., zgrasratn/ no3rat3, zgrasratp/ po4rat3, zgrasratf / ferat3)
                zbeta     = MAX( 0., (epsher - epshermin) )
+               ! Food density deprivation of GGE
                zepsherf  = epshermin + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta )
+               ! Food quality deprivation of GGE
                zepsherq  = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 )
+               ! Actual GGE
                zepsherv  = zepsherf * zepshert * zepsherq
 
                ! Respiration of microzooplankton
                ! Excess carbon in the food is used preferentially
+               ! when activated by zmetexcess
                ! ------------------------------------------------
                zexcess  = zgraztotc * zepsherf * (1.0 - zepshert) * zmetexcess
@@ -380 +395 @@
          &                xprefp, xprefd, xprefz, xthreshdia, xthreshphy, &
          &                xthreshpic, xthreshpoc, xthreshzoo, xthresh, xkgraz, &
-         &                epsher, epshermin, ssigma, srespir, unassc, unassn, unassp
+         &                epsher, epshermin, ssigma, srespir, unassc, unassn, unassp,   &
+         &                xsigma, xsigmadel   
       !!----------------------------------------------------------------------
       !
@@ -420 +436 @@
          WRITE(numout,*) '    half sturation constant for grazing 1           xkgraz      =', xkgraz
          WRITE(numout,*) '    Use of excess carbon for respiration            bmetexc     =', bmetexc
+         WRITE(numout,*) '      Width of the grazing window                     xsigma      =', xsigma
+         WRITE(numout,*) '      Maximum additional width of the grazing window  xsigmadel   =', xsigmadel
       ENDIF
       !