Commit 4e3703cc authored by Dany Dumont's avatar Dany Dumont

Ajout de la remineralisation de ammonium vers nitrate a un taux lineaire proportionnel a mu5

parent 6d7adc84
......@@ -149,6 +149,7 @@
REALTYPE :: bertha = 0.05
REALTYPE :: parcrit = 0.02
REALTYPE :: pmin = 0.05
REALTYPE :: mu5 = 0.00
REALTYPE, public :: kc = 0.03
integer :: out_unit
integer, parameter :: p=1,z=2,b=3,d=4,n=5,a=6,l=7,f=8
......@@ -190,7 +191,7 @@
vb,k4,mu3,eta,mu4,w_d,kc, &
I_opt,inib, &
theta,w_pmax,w_pmin, &
bertha,parcrit,w_zmax,pmin
bertha,parcrit,w_zmax,pmin,mu5
!EOP
!-----------------------------------------------------------------------
!BOC
......@@ -230,6 +231,7 @@
mu2 = mu2 /secs_pr_day
mu3 = mu3 /secs_pr_day
mu4 = mu4 /secs_pr_day
mu5 = mu5 /secs_pr_day
gmax = gmax /secs_pr_day
w_p = w_p /secs_pr_day
w_pmin = w_pmin /secs_pr_day !DD
......@@ -595,6 +597,7 @@
dd(d,z,ci)=beta*gmax*r3*cc(d,ci)**2*fac
dd(b,a,ci)=mu3*cc(b,ci)
dd(d,l,ci)=mu4*cc(d,ci)
dd(a,n,ci)=mu5*cc(a,ci)
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)
......@@ -609,7 +612,7 @@
!(food) concentration
!ws(z,ci) = -1.0*w_zmax*tanh(bertha*(par(ci)-parcrit))
if ((par(ci).lt.parcrit).and.(cc(p,ci) .gt. pmin)) then
ws(z,ci) = 0.0
ws(z,ci) = 0.0
else
ws(z,ci) = -1.0*w_zmax*tanh(bertha*(par(ci)-parcrit))
end if
......
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