Commit 5652bb2b authored by Dany Dumont's avatar Dany Dumont

Switch logique pour rendre les hydrocarures toxique pour le zooplancton dans bio_gsj

parent 74bf91b3
......@@ -111,6 +111,7 @@
g2max = 1.0
k3 = 1.0
beta = 0.625
txc4zoo = .true.
mu21 = 0.3
mu21max = 1.5
bhc = 11.1
......
......@@ -64,7 +64,8 @@
REALTYPE :: p0 = 0.0
REALTYPE :: z0 = 0.0
REALTYPE :: b0 = 0.0
LOGICAL :: mte = .true.
LOGICAL :: mte = .false.
LOGICAL :: txc4zoo = .false.
REALTYPE :: ca1 = 3.61
REALTYPE :: ca2 = 14.58
REALTYPE :: ch1 = 3.265
......@@ -177,6 +178,7 @@
! mu22 : Maximum large zooplankton loss rate
! k6 : Half saturation constant of zooplankton loss (small and large)
!! k7 : Half saturation constant of zooplankton loss by hydrocarbon toxicity (small and large)
! txc4zoo : Logical to make hydrocarbon toxic for zooplankton
! bhc : Slope factor of zooplankton mortality for sigmoïd exponential type curve
! mhc : Factor of zooplankton mortality for Michaelis-Menten exponential type curve
! hc0 : Hydrocarbon sill concentration for zooplankton mortality by hydrocarbons
......@@ -241,7 +243,7 @@
r11,r12,r13,r14,r21,r22,r23,r24, &
vb1,vb2,k4,k10,w_h,mu3,etaa,etah,mu4,w_d,kc,mu5, &
theta,w_p1max,w_p1min,w_p2min,w_p2max, &
mte,ca1,ca2,ch1,ch2,amratio,hmratio
txc4zoo,mte,ca1,ca2,ch1,ch2,amratio,hmratio
!EOP
!-----------------------------------------------------------------------
......@@ -766,10 +768,15 @@
nitlim2(ci) =qn2
ammlim2(ci) =qa2
! Hydrocarbon inhibition factors
!mu21hc = mu21+(mu21max-mu21)/(1d0 + exp(-1d0*(cc(hc,ci)-k6hc)/bhc))
mu21hc = mu21+((mu21max-mu21)*(cc(hc,ci))**mhc)/((cc(hc,ci))**mhc + k6hc**mhc)
mu22hc = mu21hc
! Hydrocarbon inhibition factors
!mu21hc = mu21+(mu21max-mu21)/(1d0 + exp(-1d0*(cc(hc,ci)-k6hc)/bhc))
if (txc4zoo) then
mu21hc = mu21+((mu21max-mu21)*(cc(hc,ci))**mhc)/((cc(hc,ci))**mhc + k6hc**mhc)
mu22hc = mu21hc
else
mu21hc = mu21
mu22hc = mu22
endif
! Nutrient uptake by pico- and nano-phytoplankton
dd(n,p1,ci) =amr1*ff1*qn1*(cc(p1,ci)+p0)
......
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