Commit b7e38dff0f7e4c984bb667c59730bf4dc406969d

Authored by dumoda01
1 parent 4aaa729a
Exists in master and in 1 other branch snow

Ajout d'un terme de broutage du p1 par le z1 avec sa preference r14 correspondan…

…te. Ajout de amratio et hmratio comme parametres de la namelist bio_ismer.nml.
Showing 2 changed files with 41 additions and 30 deletions   Show diff stats
nml/bio_ismer.nml
@@ -19,39 +19,42 @@ @@ -19,39 +19,42 @@
19 ! ca2 = temp-dependence coeff for p2 19 ! ca2 = temp-dependence coeff for p2
20 ! ch1 = temp-dependence coeff for z1 20 ! ch1 = temp-dependence coeff for z1
21 ! ch2 = temp-dependence coeff for z2 21 ! ch2 = temp-dependence coeff for z2
  22 +! amratio = Mass ratio between p2 and p1
  23 +! hmratio = Mass ratio between z2 and z1
22 ! vp1 = maximum flagellate uptake rate by flagellates [1/day] 24 ! vp1 = maximum flagellate uptake rate by flagellates [1/day]
23 ! vp2 = maximum diatom uptake rate by diatoms [1/day] 25 ! vp2 = maximum diatom uptake rate by diatoms [1/day]
24 ! alpha1 = slope of the flagellate PI-curve [m2/(W day)] 26 ! alpha1 = slope of the flagellate PI-curve [m2/(W day)]
25 ! alpha2 = slope of the diatom PI-curve [m2/(W day)] 27 ! alpha2 = slope of the diatom PI-curve [m2/(W day)]
26 ! inib1 = inhibition slope of the flagellate PI-curve (pos.) [m2/(W day)] 28 ! inib1 = inhibition slope of the flagellate PI-curve (pos.) [m2/(W day)]
27 ! inib2 = inhibition slope of the PI-curve (pos.) [m2/(W day)] 29 ! inib2 = inhibition slope of the PI-curve (pos.) [m2/(W day)]
28 -! kn1 = half sat. constant nitrate uptake by fla [mmol n/m3]  
29 -! ka1 = half sat. constant ammonium uptake by fla [mmol n/m3]  
30 -! kn2 = half sat. constant nitrate uptake by diatoms [mmol n/m3]  
31 -! ka2 = half sat. constant ammonium uptake by diatoms [mmol n/m3]  
32 -! mu11 = pico-phytoplankton mortality rate [1/day]  
33 -! mu12 = nano-phytoplankton mortality rate [1/day] 30 +! kn1 = half sat. constant nitrate uptake by p1 [mmol n/m3]
  31 +! ka1 = half sat. constant ammonium uptake by p1 [mmol n/m3]
  32 +! kn2 = half sat. constant nitrate uptake by p2 [mmol n/m3]
  33 +! ka2 = half sat. constant ammonium uptake by p2 [mmol n/m3]
  34 +! mu11 = mortality rate for p1 [1/day]
  35 +! mu12 = mortality rate for p2 [1/day]
34 ! k5 = half sat. constant phy. mortality [mmol n/m3] 36 ! k5 = half sat. constant phy. mortality [mmol n/m3]
35 ! gamma = exudation fraction [-] 37 ! gamma = exudation fraction [-]
36 -! w_p1 = flagellate settling velocity [m/day]  
37 -! w_p2 = diatom settling velocity [m/day] 38 +! w_p1 = settling velocity for p1 [m/day]
  39 +! w_p2 = settling velocity for p2 [m/day]
38 ! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)] 40 ! theta = phytoplancton buoyancy parameter [m3 day/(mmol N)]
39 -! g1max = maximum microzooplankton ingestion rate [1/day]  
40 -! g2max = maximum mesozooplankton ingestion rate [1/day] 41 +! g1max = maximum ingestion rate for z1 [1/day]
  42 +! g2max = maximum ingestion rate for z2 [1/day]
41 ! k3 = half saturation constant ingestion [mmol n/m3] 43 ! k3 = half saturation constant ingestion [mmol n/m3]
42 ! beta = grazing efficiency [-] 44 ! beta = grazing efficiency [-]
43 -! k6 = half saturation zooplankton loss (mcz & msz) [mmol n/m3]  
44 -! mu21 = maximum micro-zooplankton loss rate [1/day]  
45 -! mu22 = maximum meso-zooplankton loss rate [1/day]  
46 -! delta = fractional zooplankton loss to LDON (mcz & msz) [-]  
47 -! epsi = fractional zooplankton loss to ammonium (mcz & msz) [-]  
48 -! r11 = mcz grazing preference on flagellates [-]  
49 -! r12 = mcz grazing preference on bacteria [-]  
50 -! r13 = mcz grazing preference on detritus [-]  
51 -! r21 = msz grazing preference on flagellates [-]  
52 -! r22 = msz grazing preference on diatoms [-]  
53 -! r23 = msz grazing preference on detritus [-]  
54 -! r24 = msz grazing preference on mcz [-] 45 +! k6 = half saturation zooplankton loss (z1 & z2) [mmol n/m3]
  46 +! mu21 = maximum loss rate for z1 [1/day]
  47 +! mu22 = maximum loss rate for z2 [1/day]
  48 +! delta = fractional zooplankton loss to LDON (z1 & z2) [-]
  49 +! epsi = fractional zooplankton loss to ammonium (z1 & z2) [-]
  50 +! r11 = z1 grazing preference on p1 [-]
  51 +! r12 = z1 grazing preference on p2 [-]
  52 +! r13 = z1 grazing preference on bacteria [-]
  53 +! r14 = z1 grazing preference on detritus [-]
  54 +! r21 = z2 grazing preference on p1 [-]
  55 +! r22 = z2 grazing preference on p2 [-]
  56 +! r23 = z2 grazing preference on detritus [-]
  57 +! r24 = z2 grazing preference on z1 [-]
55 ! vb = maximum bacterial uptake rate [1/day] 58 ! vb = maximum bacterial uptake rate [1/day]
56 ! k4 = half saturation bacterial uptake [mmol n/m3] 59 ! k4 = half saturation bacterial uptake [mmol n/m3]
57 ! mu3 = bacteria excretion rate [1/day] 60 ! mu3 = bacteria excretion rate [1/day]
@@ -78,6 +81,8 @@ @@ -78,6 +81,8 @@
78 ca2 = 14.58 81 ca2 = 14.58
79 ch1 = 3.265 82 ch1 = 3.265
80 ch2 = 24.923 83 ch2 = 24.923
  84 + amratio = 200
  85 + hmratio = 1000
81 vp1 = 0.02 86 vp1 = 0.02
82 vp2 = 0.8 87 vp2 = 0.8
83 alpha1 = 0.02 88 alpha1 = 0.02
@@ -109,8 +114,9 @@ @@ -109,8 +114,9 @@
109 delta = 0.1 114 delta = 0.1
110 epsi = 0.70 115 epsi = 0.70
111 r11 = 0.55 116 r11 = 0.55
112 - r12 = 0.40 117 + r12 = 0.30
113 r13 = 0.05 118 r13 = 0.05
  119 + r14 = 0.10
114 r21 = 0.50 120 r21 = 0.50
115 r22 = 0.30 121 r22 = 0.30
116 r23 = 0.05 122 r23 = 0.05
src/extras/bio/bio_ismer.F90
@@ -69,6 +69,8 @@ @@ -69,6 +69,8 @@
69 REALTYPE :: ca2 = 14.58 69 REALTYPE :: ca2 = 14.58
70 REALTYPE :: ch1 = 3.265 70 REALTYPE :: ch1 = 3.265
71 REALTYPE :: ch2 = 24.923 71 REALTYPE :: ch2 = 24.923
  72 + REALTYPE :: amratio = 200
  73 + REALTYPE :: hmratio = 1000
72 REALTYPE :: vp1 = 1.5 74 REALTYPE :: vp1 = 1.5
73 REALTYPE :: alpha1 = 0.065 75 REALTYPE :: alpha1 = 0.065
74 REALTYPE :: inib1 = 0.05 76 REALTYPE :: inib1 = 0.05
@@ -100,8 +102,9 @@ @@ -100,8 +102,9 @@
100 REALTYPE :: delta = 0.1 102 REALTYPE :: delta = 0.1
101 REALTYPE :: epsi = 0.70 103 REALTYPE :: epsi = 0.70
102 REALTYPE :: r11 = 0.55 104 REALTYPE :: r11 = 0.55
103 - REALTYPE :: r12 = 0.4 105 + REALTYPE :: r12 = 0.30
104 REALTYPE :: r13 = 0.05 106 REALTYPE :: r13 = 0.05
  107 + REALTYPE :: r14 = 0.10
105 REALTYPE :: r21 = 0.50 108 REALTYPE :: r21 = 0.50
106 REALTYPE :: r22 = 0.30 109 REALTYPE :: r22 = 0.30
107 REALTYPE :: r23 = 0.05 110 REALTYPE :: r23 = 0.05
@@ -153,10 +156,10 @@ @@ -153,10 +156,10 @@
153 p0,z0,b0,vp1,alpha1,inib1,vp2,alpha2,inib2, & 156 p0,z0,b0,vp1,alpha1,inib1,vp2,alpha2,inib2, &
154 kn1,ka1,kn2,ka2,mu11,mu12,k5,gamma,w_p1,w_p2, & 157 kn1,ka1,kn2,ka2,mu11,mu12,k5,gamma,w_p1,w_p2, &
155 g1max,g2max,k3,beta,mu21,mu22,k6,delta,epsi, & 158 g1max,g2max,k3,beta,mu21,mu22,k6,delta,epsi, &
156 - r11,r12,r13,r21,r22,r23,r24, & 159 + r11,r12,r13,r14,r21,r22,r23,r24, &
157 vb,k4,mu3,eta,mu4,w_d,kc,mu5, & 160 vb,k4,mu3,eta,mu4,w_d,kc,mu5, &
158 theta,w_p1max,w_p1min,w_p2min,w_p2max, & 161 theta,w_p1max,w_p1min,w_p2min,w_p2max, &
159 - mte,ca1,ca2,ch1,ch2 162 + mte,ca1,ca2,ch1,ch2,amratio,hmratio
160 163
161 !EOP 164 !EOP
162 !----------------------------------------------------------------------- 165 !-----------------------------------------------------------------------
@@ -202,6 +205,7 @@ @@ -202,6 +205,7 @@
202 write(10,901) r11 205 write(10,901) r11
203 write(10,901) r12 206 write(10,901) r12
204 write(10,901) r13 207 write(10,901) r13
  208 + write(10,901) r14
205 write(10,901) r21 209 write(10,901) r21
206 write(10,901) r22 210 write(10,901) r22
207 write(10,901) r23 211 write(10,901) r23
@@ -219,6 +223,8 @@ @@ -219,6 +223,8 @@
219 write(10,901) ca2 223 write(10,901) ca2
220 write(10,901) ch1 224 write(10,901) ch1
221 write(10,901) ch2 225 write(10,901) ch2
  226 + write(10,901) amratio
  227 + write(10,901) hmratio
222 endif 228 endif
223 229
224 230
@@ -593,8 +599,6 @@ @@ -593,8 +599,6 @@
593 ! their respective metabolic rates, according to the MTE. 599 ! their respective metabolic rates, according to the MTE.
594 ! The same applies to micro-zooplankton (z1) versus 600 ! The same applies to micro-zooplankton (z1) versus
595 ! meso-zooplankton (z2) through fac4=hmratio**0.25. 601 ! meso-zooplankton (z2) through fac4=hmratio**0.25.
596 - amratio = 200.0  
597 - hmratio = 1000.0  
598 602
599 fac3 = amratio**0.25 603 fac3 = amratio**0.25
600 fac4 = hmratio**0.25 604 fac4 = hmratio**0.25
@@ -677,8 +681,9 @@ @@ -677,8 +681,9 @@
677 dd(b,d,ci) =(1.-beta)*g1max*r12*cc(b,ci)**2*fac1 681 dd(b,d,ci) =(1.-beta)*g1max*r12*cc(b,ci)**2*fac1
678 682
679 dd(p1,z1,ci)=hmr1*beta*g1max*r11*cc(p1,ci)**2*fac1 683 dd(p1,z1,ci)=hmr1*beta*g1max*r11*cc(p1,ci)**2*fac1
680 - dd(b,z1,ci) =hmr1*beta*g1max*r12*cc(b,ci)**2*fac1  
681 - dd(d,z1,ci) =hmr1*beta*g1max*r13*cc(d,ci)**2*fac1 684 + dd(p2,z1,ci)=hmr1*beta*g1max*r12*cc(p2,ci)**2*fac1
  685 + dd(b,z1,ci) =hmr1*beta*g1max*r13*cc(b,ci)**2*fac1
  686 + dd(d,z1,ci) =hmr1*beta*g1max*r14*cc(d,ci)**2*fac1
682 687
683 dd(p1,z2,ci)=hmr2*beta*g2max*r21*cc(p1,ci)**2*fac2 688 dd(p1,z2,ci)=hmr2*beta*g2max*r21*cc(p1,ci)**2*fac2
684 dd(p2,z2,ci)=hmr2*beta*g2max*r22*cc(p2,ci)**2*fac2 689 dd(p2,z2,ci)=hmr2*beta*g2max*r22*cc(p2,ci)**2*fac2