Changeset 12349 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmeso.F90

Timestamp:

2020-01-28T17:17:13+01:00 (4 years ago)

Author:

aumont

Message:

update of the quota version of PISCES + some corrections of the GGE of zooplankton

File:

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

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

@@ -72 +72 @@
       REAL(wp) :: zgraze2, zdenom, zfact, zfood, zfoodlim, zproport
       REAL(wp) :: zmortzgoc, zfracc, zfracn, zfracp, zfracfe, zratio, zratio2
-      REAL(wp) :: zepsherf, zepshert, zepsherv, zrespirc, zrespirn, zrespirp, zbasresb, zbasresi
+      REAL(wp) :: zepsherf, zepshert, zepsherq, zepsherv, zrespirc, zrespirn, zrespirp, zbasresb, zbasresi
       REAL(wp) :: zgraztotc, zgraztotn, zgraztotp, zgraztotf, zbasresn, zbasresp, zbasresf
       REAL(wp) :: zgradoc, zgradon, zgradop, zgratmp, zgradoct, zgradont, zgrareft, zgradopt
@@ -85 +85 @@
       REAL(wp) :: zgrazffnp, zgrazffng, zgrazffpp, zgrazffpg
       CHARACTER (len=25) :: charout
-      REAL(wp) :: zrfact2, zmetexcess
+      REAL(wp) :: zrfact2, zmetexcess, zsigma, zdiffdn
       REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing, zfezoo2
       REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d, zz2ligprod
@@ -140 +140 @@
                !   An active switching parameterization is used here.
                !   We don't use the KTW parameterization proposed by 
-               !   Vallina et al. because it tends to produce to steady biomass
+               !   Vallina et al. because it tends to produce too steady biomass
                !   composition and the variance of Chl is too low as it grazes
-               !   too strongly on winning organisms. Thus, instead of a square
-               !   a 1.5 power value is used which decreases the pressure on the
-               !   most abundant species
+               !   too strongly on winning organisms. We use a generalized
+               !   switching parameterization proposed by Morozov and 
+               !   Petrovskii (2013)
                !   ------------------------------------------------------------  
-               ztmp1 = xpref2n * zcompaph**1.5
-               ztmp2 = xpref2m * zcompames**1.5
-               ztmp3 = xpref2c * zcompapoc**1.5
-               ztmp4 = xpref2d * zcompadi**1.5
-               ztmp5 = xpref2z * zcompaz**1.5
+               zsigma = 1.0 - zdenom**2/(0.05**2+zdenom**2)
+               zsigma = 0.5 + 1.0 * zsigma
+               zdiffdn = exp( -ABS(log(3.0 * sizen(ji,jj,jk) / (5.0 * sized(ji,jj,jk) + rtrn )) )**2 / zsigma**2 )
+               ztmp1 = xpref2n * zcompaph * ( zcompaph + zdiffdn * zcompadi )
+               ztmp2 = xpref2m * zcompames**2
+               ztmp3 = xpref2c * zcompapoc**2
+               ztmp4 = xpref2d * zcompadi * ( zdiffdn * zcompadi + zcompaph )
+               ztmp5 = xpref2z * zcompaz**2
                ztmptot = ztmp1 + ztmp2 + ztmp3 + ztmp4 + ztmp5 + rtrn
                ztmp1 = ztmp1 / ztmptot
@@ -238 +241 @@
                zbeta     = MAX(0., (epsher2 - epsher2min) )
                zepsherf  = epsher2min + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta )
-               zepsherv  = zepsherf * zepshert
-
+               zepsherq  = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 )
+               zepsherv  = zepsherf * zepshert * zepsherq
                !   Respiration of mesozooplankton
                !   Excess carbon in the food is used preferentially