Commit 1025a7c4 authored by dumoda01's avatar dumoda01

Ajout de diagnostics utiles : lumlim, nitlim, ammlim

parent 62daa075
......@@ -115,9 +115,12 @@
REALTYPE :: b0 = 0.0
REALTYPE :: vp = 1.5
REALTYPE :: alpha = 0.065
REALTYPE :: I_opt = 10.0 !CHG1
REALTYPE :: inib = 0.05 !CHG1
REALTYPE :: theta = 0.0 !CHG2
!CHG1 Intensité optimale (I_opt))
REALTYPE :: I_opt = 10.0
!CHG1 Paramere d'inibition (inib)
REALTYPE :: inib = 0.05
!CHG2 Paramere de flottaison du phytoplancton
REALTYPE :: theta = 0.0
REALTYPE :: w_pmin = -0.06
REALTYPE :: w_pmax = -0.38
REALTYPE :: k1 = 0.2
......@@ -261,9 +264,10 @@
cc(b,i)=b_initial
cc(d,i)=d_initial
! cc(n,i)=n_initial
! cc(n,i)=nit(i) !CHG3
cc(n,i)=nprof(i) !CHG3
! cc(a,i)=a_initial
cc(n,i)=aprof(i) !CHG5
! cc(n,i)=aprof(i) !CHG5
cc(a,i)=a_initial !CHG5
cc(l,i)=l_initial
end do
......@@ -510,7 +514,7 @@
! Original author(s): Hans Burchard, Karsten Bolding
!
! !LOCAL VARIABLES:
REALTYPE :: ff,fac,fac2,min67,sigma
REALTYPE :: ff,fac,fac2,min67,q1,q2
REALTYPE :: Ps !CHG1
integer :: i,j,ci
!EOP
......@@ -540,17 +544,21 @@
ff= Ps*(1.-exp(-1.*alpha*par(ci)/Ps))*exp(-1.*inib*par(ci)/Ps)
! --------------------------------------------------------------------
q1=(cc(n,ci)/k1)/(1.+cc(n,ci)/k1+cc(a,ci)/k2)
q2=(cc(a,ci)/k2)/(1.+cc(n,ci)/k1+cc(a,ci)/k2)
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 (Kuhn and Radach 1997)
sigma =ff*(cc(n,ci)/k1+cc(a,ci)/k2)/(1.+cc(n,ci)/k1+cc(a,ci)/k2)
! Light and nutrient limitation factors
lumlim(ci)=ff
nitlim(ci)=q1
ammlim(ci)=q2
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)/ &
(1.+cc(n,ci)/k1+cc(a,ci)/k2)*cc(p,ci)
dd(p,l,ci)=gamma*ff*(q1+q2)*cc(p,ci)
dd(b,d,ci)=(1.-beta)*gmax*r2*cc(b,ci)**2*fac
dd(p,z,ci)=beta*gmax*r1*cc(p,ci)**2*fac
dd(b,z,ci)=beta*gmax*r2*cc(b,ci)**2*fac
......@@ -560,8 +568,8 @@
dd(z,d,ci)=(1.-epsi-delta)*mu2*(cc(z,ci)+z0)/(k6+cc(z,ci)+z0)*cc(z,ci)
dd(z,a,ci)=epsi*mu2*(cc(z,ci)+z0)/(k6+cc(z,ci)+z0)*cc(z,ci)
dd(z,l,ci)=delta*mu2*(cc(z,ci)+z0)/(k6+cc(z,ci)+z0)*cc(z,ci)
dd(n,p,ci)=ff*cc(n,ci)/k1/(1.+cc(n,ci)/k1+cc(a,ci)/k2)*(cc(p,ci)+p0)
dd(a,p,ci)=ff*cc(a,ci)/k2/(1.+cc(n,ci)/k1+cc(a,ci)/k2)*(cc(p,ci)+p0)
dd(n,p,ci)=ff*q1*(cc(p,ci)+p0)
dd(a,p,ci)=ff*q2*(cc(p,ci)+p0)
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)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment