Commit ac88286bd5135039ca1d2a1fdefa97c214a74b53

Authored by dumoda01
1 parent 8ea84af1
Exists in master and in 1 other branch snow

Ajout de broutage du mesozoo sur le microzoo (preference r24).

Showing 2 changed files with 114 additions and 108 deletions   Show diff stats
nml/bio_ismer.nml
... ... @@ -4,109 +4,111 @@
4 4 !
5 5 ! numc= number of compartments for geobiochemical model
6 6 !
7   -! p1_initial= initial flagellate concentration [mmol n/m3]
8   -! p2_initial= initial diatom concentration [mmol n/m3]
9   -! z1_initial= initial micro-zooplakton concentration [mmol n/m3]
10   -! z2_initial= initial meso-zooplakton concentration [mmol n/m3]
11   -! b_initial= initial bacteria concentration [mmol n/m3]
12   -! d_initial= initial detritus concentration [mmol n/m3]
13   -! n_initial= *** see obs.nml *** [mmol n/m3]
14   -! a_initial= *** see obs.nml *** [mmol n/m3]
15   -! l_initial= initial LDON concentration [mmol n/m3]
16   -! p0 = minimum phytoplankton concentration [mmol n/m3]
17   -! z0 = minimum zooplakton concentration [mmol n/m3]
18   -! b0 = minimum bacteria concentration [mmol n/m3]
19   -! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] !CHG2
20   -! vp1 = maximum flagellate uptake rate by flagellates [1/day]
21   -! vp2 = maximum diatom uptake rate by diatoms [1/day]
22   -! alpha1 = slope of the flagellate PI-curve [m2/(W day)]
23   -! alpha2 = slope of the diatom PI-curve [m2/(W day)]
24   -! inib1 = inhibition slope of the flagellate PI-curve (pos.) [m2/(W day)] !CHG1
25   -! inib2 = inhibition slope of the PI-curve (pos.) [m2/(W day)] !CHG1
26   -! kn1 = half sat. constant nitrate uptake by fla [mmol n/m3]
27   -! ka1 = half sat. constant ammonium uptake by fla [mmol n/m3]
28   -! kn2 = half sat. constant nitrate uptake by diatoms [mmol n/m3]
29   -! ka2 = half sat. constant ammonium uptake by diatoms [mmol n/m3]
30   -! mu1 = phytoplankton (fla & dia) mortality rate [1/day]
31   -! k5 = half sat. constant phy. mortality (fla & dia) [mmol n/m3]
32   -! gamma = exudation fraction [-]
33   -! w_p1 = flagellate settling velocity [m/day]
34   -! w_p2 = diatom settling velocity [m/day]
35   -! g1max = maximum microzooplankton ingestion rate [1/day]
36   -! g2max = maximum mesozooplankton ingestion rate [1/day]
37   -! k3 = half saturation constant ingestion [mmol n/m3]
38   -! beta = grazing efficiency [-]
39   -! mu2 = maximum zooplankton loss rate (mcz & msz) [1/day]
40   -! k6 = half saturation zooplankton loss (mcz & msz) [mmol n/m3]
41   -! delta = fractional zooplankton loss to LDON (mcz & msz) [-]
42   -! epsi = fractional zooplankton loss to ammonium (mcz & msz) [-]
43   -! r11 = mcz grazing preference on flagellates [-]
44   -! r12 = mcz grazing preference on bacteria [-]
45   -! r13 = mcz grazing preference on detritus [-]
46   -! r21 = msz grazing preference on flagellates [-]
47   -! r22 = msz grazing preference on diatoms [-]
48   -! r23 = msz grazing preference on detritus [-]
49   -! vb = maximum bacterial uptake rate [1/day]
50   -! k4 = half saturation bacterial uptake [mmol n/m3]
51   -! mu3 = bacteria excretion rate [1/day]
52   -! eta = uptake ratio ammonium:LDON [-]
53   -! mu4 = detritus breakdown rate [1/day]
54   -! mu5 = nitrification rate [1/day]
55   -! w_d = detritus settling velocity [m/day]
56   -! kc = attenuation constant for the self shading effect [m**2/mmol N]
  7 +! p1_init = initial flagellate concentration [mmol n/m3]
  8 +! p2_init = initial diatom concentration [mmol n/m3]
  9 +! z1_init = initial micro-zooplakton concentration [mmol n/m3]
  10 +! z2_init = initial meso-zooplakton concentration [mmol n/m3]
  11 +! b_init = initial bacteria concentration [mmol n/m3]
  12 +! d_init = initial detritus concentration [mmol n/m3]
  13 +! n_init = *** see obs.nml *** [mmol n/m3]
  14 +! a_init = *** see obs.nml *** [mmol n/m3]
  15 +! l_init = initial LDON concentration [mmol n/m3]
  16 +! p0 = minimum phytoplankton concentration [mmol n/m3]
  17 +! z0 = minimum zooplakton concentration [mmol n/m3]
  18 +! b0 = minimum bacteria concentration [mmol n/m3]
  19 +! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] !CHG2
  20 +! vp1 = maximum flagellate uptake rate by flagellates [1/day]
  21 +! vp2 = maximum diatom uptake rate by diatoms [1/day]
  22 +! alpha1 = slope of the flagellate PI-curve [m2/(W day)]
  23 +! alpha2 = slope of the diatom PI-curve [m2/(W day)]
  24 +! inib1 = inhibition slope of the flagellate PI-curve (pos.) [m2/(W day)] !CHG1
  25 +! inib2 = inhibition slope of the PI-curve (pos.) [m2/(W day)] !CHG1
  26 +! kn1 = half sat. constant nitrate uptake by fla [mmol n/m3]
  27 +! ka1 = half sat. constant ammonium uptake by fla [mmol n/m3]
  28 +! kn2 = half sat. constant nitrate uptake by diatoms [mmol n/m3]
  29 +! ka2 = half sat. constant ammonium uptake by diatoms [mmol n/m3]
  30 +! mu1 = phytoplankton (fla & dia) mortality rate [1/day]
  31 +! k5 = half sat. constant phy. mortality (fla & dia) [mmol n/m3]
  32 +! gamma = exudation fraction [-]
  33 +! w_p1 = flagellate settling velocity [m/day]
  34 +! w_p2 = diatom settling velocity [m/day]
  35 +! g1max = maximum microzooplankton ingestion rate [1/day]
  36 +! g2max = maximum mesozooplankton ingestion rate [1/day]
  37 +! k3 = half saturation constant ingestion [mmol n/m3]
  38 +! beta = grazing efficiency [-]
  39 +! mu2 = maximum zooplankton loss rate (mcz & msz) [1/day]
  40 +! k6 = half saturation zooplankton loss (mcz & msz) [mmol n/m3]
  41 +! delta = fractional zooplankton loss to LDON (mcz & msz) [-]
  42 +! epsi = fractional zooplankton loss to ammonium (mcz & msz) [-]
  43 +! r11 = mcz grazing preference on flagellates [-]
  44 +! r12 = mcz grazing preference on bacteria [-]
  45 +! r13 = mcz grazing preference on detritus [-]
  46 +! r21 = msz grazing preference on flagellates [-]
  47 +! r22 = msz grazing preference on diatoms [-]
  48 +! r23 = msz grazing preference on detritus [-]
  49 +! r24 = msz grazing preference on mcz [-]
  50 +! vb = maximum bacterial uptake rate [1/day]
  51 +! k4 = half saturation bacterial uptake [mmol n/m3]
  52 +! mu3 = bacteria excretion rate [1/day]
  53 +! eta = uptake ratio ammonium:LDON [-]
  54 +! mu4 = detritus breakdown rate [1/day]
  55 +! mu5 = nitrification rate [1/day]
  56 +! w_d = detritus settling velocity [m/day]
  57 +! kc = attenuation constant for the self shading effect [m**2/mmol N]
57 58 !-------------------------------------------------------------------------------
58 59 &bio_ismer_nml
59   - numc= 9
60   - p1_initial= 0.012
61   - p2_initial= 0.012
62   - z1_initial= 0.012
63   - z2_initial= 0.012
64   - b_initial= 0.001
65   - d_initial= 0.01
66   - l_initial= 0.1
67   - p0= 0.0001
68   - z0= 0.0001
69   - b0= 0.0001
70   - vp1= 0.3
71   - vp2= 0.3
72   - alpha1= 0.04
73   - alpha2= 0.04
74   - inib1= 0.06
75   - inib2= 0.06
76   - kn1= 1.0
77   - ka1= 0.8
78   - kn2= 1.0
79   - ka2= 0.8
80   - mu1= 0.05
81   - k5= 0.2
82   - gamma= 0.05
83   - w_p1= -0.38
84   - w_p2= -0.10
85   - theta= 0.0
86   - w_p1min= -0.06
87   - w_p1max= -0.38
88   - w_p2min= -0.06
89   - w_p2max= -0.38
90   - g1max= 1.0
91   - g2max= 1.0
92   - k3= 1.0
93   - beta= 0.625
94   - mu2= 0.3
95   - k6= 0.2
96   - delta= 0.1
97   - epsi= 0.70
98   - r11= 0.55
99   - r12= 0.4
100   - r13= 0.05
101   - r21= 0.55
102   - r22= 0.4
103   - r23= 0.05
104   - vb= 0.24
105   - k4= 0.5
106   - mu3= 0.03
107   - eta= 0.0
108   - mu4= 0.02
109   - mu5= 0.00
110   - w_d= -5.0
111   - kc= 0.03
  60 + numc = 9
  61 + p1_init = 0.012
  62 + p2_init = 0.012
  63 + z1_init = 0.012
  64 + z2_init = 0.012
  65 + b_init = 0.001
  66 + d_init = 0.01
  67 + l_init = 0.1
  68 + p0 = 0.0001
  69 + z0 = 0.0001
  70 + b0 = 0.0001
  71 + vp1 = 1.0
  72 + vp2 = 0.8
  73 + alpha1 = 0.04
  74 + alpha2 = 0.04
  75 + inib1 = 0.006
  76 + inib2 = 0.006
  77 + kn1 = 1.0
  78 + ka1 = 0.8
  79 + kn2 = 1.0
  80 + ka2 = 0.8
  81 + mu1 = 0.05
  82 + k5 = 0.2
  83 + gamma = 0.05
  84 + w_p1 =-0.10
  85 + w_p2 =-0.38
  86 + theta = 0.0
  87 + w_p1min =-0.01
  88 + w_p1max =-0.10
  89 + w_p2min =-0.05
  90 + w_p2max =-0.38
  91 + g1max = 1.0
  92 + g2max = 1.0
  93 + k3 = 1.0
  94 + beta = 0.625
  95 + mu2 = 0.3
  96 + k6 = 0.2
  97 + delta = 0.1
  98 + epsi = 0.70
  99 + r11 = 0.55
  100 + r12 = 0.40
  101 + r13 = 0.05
  102 + r21 = 0.50
  103 + r22 = 0.30
  104 + r23 = 0.05
  105 + r24 = 0.15
  106 + vb = 0.24
  107 + k4 = 0.5
  108 + mu3 = 0.03
  109 + eta = 0.0
  110 + mu4 = 0.02
  111 + mu5 = 0.00
  112 + w_d = -5.0
  113 + kc = 0.03
112 114 /
... ...
src/extras/bio/bio_ismer.F90
... ... @@ -95,9 +95,10 @@
95 95 REALTYPE :: r11 = 0.55
96 96 REALTYPE :: r12 = 0.4
97 97 REALTYPE :: r13 = 0.05
98   - REALTYPE :: r21 = 0.55
99   - REALTYPE :: r22 = 0.4
  98 + REALTYPE :: r21 = 0.50
  99 + REALTYPE :: r22 = 0.30
100 100 REALTYPE :: r23 = 0.05
  101 + REALTYPE :: r24 = 0.15
101 102 REALTYPE :: vb = 1.2
102 103 REALTYPE :: k4 = 0.5
103 104 REALTYPE :: mu3 = 0.15
... ... @@ -145,7 +146,7 @@
145 146 p0,z0,b0,vp1,alpha1,inib1,vp2,alpha2,inib2, &
146 147 kn1,ka1,kn2,ka2,mu1,k5,gamma,w_p1,w_p2, &
147 148 g1max,g2max,k3,beta,mu2,k6,delta,epsi, &
148   - r11,r12,r13,r21,r22,r23, &
  149 + r11,r12,r13,r21,r22,r23,r24, &
149 150 vb,k4,mu3,eta,mu4,w_d,kc,mu5, &
150 151 theta,w_p1max,w_p1min,w_p2min,w_p2max
151 152 !EOP
... ... @@ -547,8 +548,10 @@
547 548 fac1=(cc(z1,ci)+z0)/(k3*(r11*cc(p1,ci)+r12*cc(b,ci)+r13*cc(d,ci))+ &
548 549 r11*cc(p1,ci)**2+r12*cc(b,ci)**2+r13*cc(d,ci)**2)
549 550  
550   - fac2=(cc(z2,ci)+z0)/(k3*(r21*cc(p1,ci)+r22*cc(p2,ci)+r23*cc(d,ci))+ &
551   - r21*cc(p1,ci)**2+r22*cc(p2,ci)**2+r23*cc(d,ci)**2)
  551 + fac2=(cc(z2,ci)+z0)/(k3*(r21*cc(p1,ci)+r22*cc(p2,ci) &
  552 + +r23*cc(d,ci)+r24*cc(z1,ci)) &
  553 + +r21*cc(p1,ci)**2+r22*cc(p2,ci)**2 &
  554 + +r23*cc(d,ci)**2+r24*cc(z1,ci)**2)
552 555  
553 556 minal=min(cc(a,ci),eta*cc(l,ci))
554 557  
... ... @@ -585,6 +588,7 @@
585 588 dd(p1,z2,ci)=fac3*beta*g2max*r21*cc(p1,ci)**2*fac2
586 589 dd(p2,z2,ci)=fac3*beta*g2max*r22*cc(p2,ci)**2*fac2
587 590 dd(d,z2,ci) =fac3*beta*g2max*r23*cc(d,ci)**2*fac2
  591 + dd(z1,z2,ci)=fac3*beta*g2max*r24*cc(z1,ci)**2*fac2
588 592  
589 593 dd(b,a,ci) =mu3*cc(b,ci)
590 594 dd(d,l,ci) =mu4*cc(d,ci)
... ...