Ajout de la variable sigma = gross growth rate

src/extras/bio/bio_fasham.F90

@@ -116,11 +116,13 @@
    REALTYPE                  ::  vp       = 1.5
    REALTYPE                  ::  alpha    = 0.065
 !CHG1 Intensité optimale (I_opt))
-   REALTYPE                  ::  I_opt     = 10.0
+   REALTYPE                  ::  I_opt    = 10.0
 !CHG1 Paramere d'inibition (inib)
    REALTYPE                  ::  inib     = 0.05
 !CHG2 Paramere de flottaison du phytoplancton
-   REALTYPE                  ::  theta     = 1.e7
+   REALTYPE                  ::  theta    = 0.0
+   REALTYPE                  ::  w_pmin   = -0.06
+   REALTYPE                  ::  w_pmax   = -0.38
    REALTYPE                  ::  k1       = 0.2
    REALTYPE                  ::  k2       = 0.8
    REALTYPE                  ::  mu1      = 0.05
@@ -143,9 +145,9 @@
    REALTYPE                  ::  eta      = 0.0
    REALTYPE                  ::  mu4      = 0.02
    REALTYPE                  ::  w_d      = -2.0
-   REALTYPE, public          ::  kc=0.03
-   REALTYPE                  ::  aa=0.62
-   REALTYPE                  ::  g2=20.0
+   REALTYPE, public          ::  kc       = 0.03
+   REALTYPE                  ::  aa       = 0.70
+   REALTYPE                  ::  g2       = 14.0
    integer                   ::  out_unit
    integer, parameter        ::  p=1,z=2,b=3,d=4,n=5,a=6,l=7
 !EOP
@@ -182,9 +184,9 @@
                         p_initial,z_initial,b_initial,d_initial,n_initial, &
                         a_initial,l_initial,p0,z0,b0,vp,alpha,k1,k2,mu1,k5,&
                         gamma,w_p,gmax,k3,beta,mu2,k6,delta,epsi,r1,r2,r3, &
-                        vb,k4,mu3,eta,mu4,w_d,kc,aa,g2,                    &
+                        vb,k4,mu3,eta,mu4,w_d,kc,                          &
                         I_opt,inib,                                        &    !CHG1
-                        theta                                                   !CHG2
+                        theta,w_pmax,w_pmin                                     !CHG2
 !EOP
 !-----------------------------------------------------------------------
 !BOC
@@ -512,7 +514,7 @@
 !  Original author(s): Hans Burchard, Karsten Bolding
 !
 ! !LOCAL VARIABLES:
-   REALTYPE                   :: ff,fac,fac2,min67
+   REALTYPE                   :: ff,fac,fac2,min67,sigma
    REALTYPE                   :: Ps                           !CHG1
    integer                    :: i,j,ci
 !EOP
@@ -539,13 +541,16 @@
 !      ff= vp*((par(ci)/I_opt)*exp(1-(par(ci)/I_opt)))**2
 ! Platt et al. (1980) - inhibition
       Ps= vp/((alpha/(alpha+inib))*(alpha/(alpha+inib))**(inib/alpha))
-      ff= Ps*(1-exp(-1*alpha*par(ci)/Ps))*exp(-1*inib*par(ci)/Ps)
+      ff= Ps*(1.-exp(-1.*alpha*par(ci)/Ps))*exp(-1.*inib*par(ci)/Ps)
 ! --------------------------------------------------------------------
 
       fac=(cc(z,ci)+z0)/(k3*(r1*cc(p,ci)+r2*cc(b,ci)+r3*cc(d,ci))+  &
                       r1*cc(p,ci)**2+r2*cc(b,ci)**2+r3*cc(d,ci)**2)
       min67=min(cc(a,ci),eta*cc(l,ci))
 
+      ! Gross growth rate / Taux de croissance brut (Kunh and Radach 1997)
+      sigma     =ff*(cc(n,ci)/k1+cc(a,ci)/k2)/(1.+cc(n,ci)/k1+cc(a,ci)/k2)
+
       dd(p,d,ci)=mu1*(cc(p,ci)+p0)/(k5+cc(p,ci)+p0)*cc(p,ci)  &
                  +(1.-beta)*gmax*r1*cc(p,ci)**2*fac
       dd(p,l,ci)=gamma*ff*(cc(n,ci)/k1+cc(a,ci)/k2)/    &
@@ -564,9 +569,9 @@
       dd(a,b,ci)=vb*min67/(k4+min67+cc(l,ci))*(cc(b,ci)+b0)
       dd(l,b,ci)=vb*cc(l,ci)/(k4+min67+cc(l,ci))*(cc(b,ci)+b0)
 
-! Update de la vitesse de chute du phytoplancton en fonction de dd(n,p,ci) et dd(a,p,ci)
-      ws(p,ci)=w_p*exp(-theta*(dd(n,p,ci)+dd(a,p,ci))/(cc(n,ci)+p0))    !CHG2
-!      ws(p,ci)=w_p
+! Taux de chute du phytoplancton en fonction de la croissance
+!      ws(p,ci)=w_pmin+(w_pmax-w_pmin)*exp(-theta*sigma)
+      ws(p,ci)=w_p
 
       do i=1,numc
          do j=1,numc