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

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

@@ -45 +45 @@
    REAL(wp), PUBLIC ::  epsher2      !: growth efficiency
    REAL(wp), PUBLIC ::  epsher2min   !: minimum growth efficiency at high food for grazing 2
+   REAL(wp), PUBLIC ::  xsigma2      !: Width of the predation window
+   REAL(wp), PUBLIC ::  xsigma2del   !: Maximum width of the predation window at low food density
    REAL(wp), PUBLIC ::  grazflux     !: mesozoo flux feeding rate
    REAL(wp), PUBLIC ::  xfracmig     !: Fractional biomass of meso that performs DVM
@@ -169 +171 @@
                ! ----------------------------------------------------------
                zsigma = 1.0 - zdenom**2/(0.05**2+zdenom**2)
-               zsigma = 0.5 + 1.0 * zsigma
+               zsigma = xsigma2 + xsigma2del * zsigma
                ! Nanophytoplankton and diatoms are the only preys considered
                ! to be close enough to have potential interference
@@ -186 +188 @@
                !   Mesozooplankton regular grazing on the different preys
                !   ------------------------------------------------------
-               zgrazdc   = zgraze2  * ztmp3 * zdenom
-               zgraznc   = zgraze2  * ztmp1 * zdenom
-               zgrazpoc  = zgraze2  * ztmp2 * zdenom
-               zgrazz    = zgraze2  * ztmp4 * zdenom
-
+               zgrazdc   = zgraze2  * ztmp3 * zdenom  ! diatoms
+               zgraznc   = zgraze2  * ztmp1 * zdenom  ! nanophytoplankton
+               zgrazpoc  = zgraze2  * ztmp2 * zdenom  ! small POC
+               zgrazz    = zgraze2  * ztmp4 * zdenom  ! microzooplankton
+
+               ! Ingestion rates of the Fe content of the different preys
                zgraznf   = zgraznc  * trb(ji,jj,jk,jpnfe) / ( trb(ji,jj,jk,jpphy) + rtrn)
                zgrazdf   = zgrazdc  * trb(ji,jj,jk,jpdfe) / ( trb(ji,jj,jk,jpdia) + rtrn)
@@ -228 +231 @@
                zgrazfffp = zproport * zgrazfffp
                zgrazfffg = zproport * zgrazfffg
+
+               ! Total ingestion rates in C, N, Fe
                zgraztotc = zgrazdc + zgrazz + zgraznc + zgrazpoc + zgrazffep + zgrazffeg
                zgraztotn = zgrazdc * quotad(ji,jj,jk) + zgrazz + zgraznc * quotan(ji,jj,jk)   &
@@ -248 +253 @@
                zepshert  = MIN( 1., zgrasratn, zgrasratf / ferat3)
                zbeta     = MAX(0., (epsher2 - epsher2min) )
+               ! Food quantity deprivation of GGE
                zepsherf  = epsher2min + 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
                ! 
@@ -258 +266 @@
                ! according to a infinite chain of predators (ANderson et al., 2013)
                zmortzgoc = unass2 / ( 1. - epsher2 ) * ztortz + zrespz
+
+               ! Update of the trends
                tra(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) - zmortz + zepsherv * zgraztotc
                tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazdc
@@ -502 +512 @@
          WRITE(numout,*) '      Fraction of mesozoo excretion as DOM           sigma2       =', sigma2
          WRITE(numout,*) '      half sturation constant for grazing 2          xkgraz2      =', xkgraz2
+         WRITE(numout,*) '      Width of the grazing window                     xsigma2     =', xsigma2
+         WRITE(numout,*) '      Maximum additional width of the grazing window  xsigma2del  =', xsigma2del
          WRITE(numout,*) '      Diurnal vertical migration of mesozoo.         ln_dvm_meso  =', ln_dvm_meso
          WRITE(numout,*) '      Fractional biomass of meso  that performs DVM  xfracmig     =', xfracmig