New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
fabm.yaml in NEMO/branches/UKMO/AMM15_v3_6_STABLE_new_ersem/NEMOGCM/CONFIG/AMM7_FABM_ERSEM/EXP2020 – NEMO

source: NEMO/branches/UKMO/AMM15_v3_6_STABLE_new_ersem/NEMOGCM/CONFIG/AMM7_FABM_ERSEM/EXP2020/fabm.yaml @ 15480

Last change on this file since 15480 was 15480, checked in by jcastill, 3 years ago

Changes as in the git branch NEMO-FABMv1-ERSEM

File size: 102.7 KB
Line 
1check_conservation: false
2require_initialization: true
3instances:
4  zenithAngle:
5    model: ersem/zenith_angle
6  light:
7    model: ersem/light_iop_ady
8    parameters:
9      a0w: 0.03                                            # absorption coefficient of clear water (1/m), default = 0.036
10      b0w: 0.0015                                          # backscatter coefficient of clear water (1/m), default = 0.0016
11      pEIR_eow: 0.5                                        # photosynthetically active fraction of shortwave radiation (-), default = 0.5
12      relax: 0.143                                         # relaxation towards satellite gelbstoff absorption (1/d), default=0.033
13    initialization:
14      ADY: 0.05
15  N1:
16    long_name: phosphate
17    model: ersem/pelagic_base
18    parameters:
19      composition: p                                       # elemental composition
20      c0: 0.001                                            # background carbon concentration (mg C/m^3), default = 0.0
21    initialization:
22      p: 0.4                                               # phosphorus (mmol P/m^3)
23  N3:
24    long_name: nitrate
25    model: ersem/pelagic_base
26    parameters:
27      composition: n                                       # elemental composition
28      c0: 0.001                                            # background carbon concentration (mg C/m^3), default = 0.0
29    initialization:
30      n: 8.0                                               # nitrogen (mmol N/m^3)
31  N4:
32    long_name: ammonium
33    model: ersem/pelagic_base
34    parameters:
35      composition: n                                       # elemental composition
36      c0: 0.001                                            # background carbon concentration (mg C/m^3), default = 0.0
37    initialization:
38      n: 0.1                                               # nitrogen (mmol N/m^3)
39  N5:
40    long_name: silicate
41    model: ersem/pelagic_base
42    parameters:
43      composition: s                                       # elemental composition
44      s0: 0.0003                                           # background silicon concentration (mmol Si/m^3), default = 0.0
45    initialization:
46      s: 4.5                                               # silicate (mmol Si/m^3)
47  O2:
48    long_name: oxygen
49    model: ersem/oxygen
50    parameters:
51      ISWO2: 2                                             # saturation formulation (1: legacy ERSEM, 2: Weiss 1970)
52    initialization:
53      o: 300.0                                             # oxygen (mmol O_2/m^3)
54  O3:
55    long_name: carbonate
56    model: ersem/carbonate
57    parameters:
58      iswCO2: 1                                            # carbonate system diagnostics (0: off, 1: on), default = 1
59      iswASFLUX: 1                                         # air-sea CO2 exchange (0: none, 1: Nightingale and Liss, 2: Wanninkhof 1992 without chemical enhancement, 3: Wanninkhof 1992 with chemical enhancement, 4: Wanninkhof 1999), default = 1
60      iswtalk: 2                                           # alkalinity formulation (1-4: from salinity and temperature, 5: dynamic alkalinity), default = 5
61      iswbioalk: 1
62    initialization:
63      c: 2130.0                                            # total dissolved inorganic carbon (mmol C/m^3)
64  R1:
65    long_name: labile dissolved organic matter
66    model: ersem/pelagic_base
67    parameters:
68      composition: cnp                                     # elemental composition
69      c0: 0.0034                                           # background carbon concentration (mg C/m^3), default = 0.0
70    initialization:
71      c: 10.0                                              # carbon (mg C/m^3)
72      n: 0.14                                              # nitrogen (mmol N/m^3)
73      p: 0.01                                              # phosphorus (mmol P/m^3)
74  R2:
75    long_name: semi-labile dissolved organic matter
76    model: ersem/pelagic_base
77    parameters:
78      composition: c                                       # elemental composition
79      c0: 0.0033                                           # background carbon concentration (mg C/m^3), default = 0.0
80    initialization:
81      c: 12.0                                              # carbon (mg C/m^3)
82  R3:
83    long_name: semi-refractory dissolved organic matter
84    model: ersem/pelagic_base
85    parameters:
86      composition: c                                       # elemental composition
87      c0: 0.0033                                           # background carbon concentration (mg C/m^3), default = 0.0
88    initialization:
89      c: 12.0                                              # carbon (mg C/m^3)
90  R4:
91    long_name: small-size pom
92    model: ersem/pelagic_base
93    parameters:
94      composition: cnpf                                    # elemental composition
95      rm: 1.0                                              # sinking velocity (m/d), default = 0.0
96      ndeposition: 2                                           # number of target pools where the sedimentation flux goes (-)
97      qxc2: 0.1
98      qxn2: 0.1
99      qxp2: 0.06
100      c0: 0.0033                                           # background carbon concentration (mg C/m^3), default = 0.0
101    initialization:
102      c: 7.2                                               # carbon (mg C/m^3)
103      n: 0.1                                               # nitrogen (mmol N/m^3)
104      p: 0.007                                             # phosphorus (mmol P/m^3)
105    coupling:
106      deposition_target1: Q6/surface
107      deposition_target2: Q7/surface
108  R6:
109    long_name: medium-size pom
110    model: ersem/pelagic_base
111    parameters:
112      composition: cnpsf                                   # elemental composition
113      rm: 5.0                                              # sinking velocity (m/d), default = 0.0
114      ndeposition: 2                                           # number of target pools where the sedimentation flux goes (-)
115      qxc2: 0.1
116      qxn2: 0.1
117      qxp2: 0.06
118      qxs2: 0.0
119      c0: 0.0033                                           # background carbon concentration (mg C/m^3), default = 0.0
120    initialization:
121      c: 17.0                                              # carbon (mg C/m^3)
122      n: 0.24                                              # nitrogen (mmol N/m^3)
123      p: 0.02                                              # phosphorus (mmol P/m^3)
124      s: 0.1                                               # silicate (mmol Si/m^3)
125    coupling:
126      deposition_target1: Q6/surface
127      deposition_target2: Q7/surface
128  R8:
129    long_name: large-size pom
130    model: ersem/pelagic_base
131    parameters:
132      composition: cnps                                    # elemental composition
133      rm: 10.0                                             # sinking velocity (m/d), default = 0.0
134      ndeposition: 2                                       # number of target pools where the sedimentation flux goes (-)
135      qxc2: 0.1
136      qxn2: 0.1
137      qxp2: 0.06
138      qxs2: 0.0
139      c0: 0.0033                                           # background carbon concentration (mg C/m^3), default = 0.0
140    initialization:
141      c: 0.17                                              # carbon (mg C/m^3)
142      n: 0.0024                                            # nitrogen (mmol N/m^3)
143      p: 0.0002                                            # phosphorus (mmol P/m^3)
144      s: 0.001                                             # silicate (mmol Si/m^3)
145    coupling:
146      deposition_target1: Q6/surface
147      deposition_target2: Q7/surface
148  Q1:
149    long_name: benthic dissolved organic matter
150    model: ersem/benthic_base
151    parameters:
152      composition: cnp                                     # elemental composition
153      c0: 0.001                                            # background carbon concentration (mg C/m^2)
154    initialization:
155      c: 18.9                                              # carbon (mg C/m^2)
156      n: 0.6                                               # nitrogen (mmol N/m^2)
157      p: 0.0035                                            # phosphorus (mmol P/m^2)
158  Q6:
159    long_name: benthic particulate organic matter
160    model: ersem/benthic_column_particulate_matter
161    parameters:
162      composition: cnpsf                                   # elemental composition
163      resuspension: true                                   # enable resuspension, default = false
164      c0: 0.001                                            # background carbon concentration (mg C/m^2)
165    initialization:
166      c: 3052.0                                            # carbon (mg C/m^2)
167      n: 3.5                                               # nitrogen (mmol N/m^2)
168      p: 0.5                                               # phosphorus (mmol P/m^2)
169      s: 149.0                                             # silicate (mmol Si/m^2)
170      pen_depth_c: 0.0173                                  # penetration depth of carbon (m)
171      pen_depth_n: 0.02                                    # penetration depth of nitrogen (m)
172      pen_depth_p: 0.025                                   # penetration depth of phosphorus (m)
173      pen_depth_s: 0.0313                                  # penetration depth of silicate (m)
174    coupling:
175      RP: R6
176  Q7:
177    long_name: benthic refractory matter
178    model: ersem/benthic_column_particulate_matter
179    parameters:
180      composition: cnp                                     # elemental composition
181      burial: true                                         # enable burial, default = false
182      c0: 0.001                                            # background carbon concentration (mg C/m^2)
183    initialization:
184      c: 30520.0                                           # carbon (mg C/m^2)
185      n: 80.0                                              # nitrogen (mmol N/m^2)
186      p: 8.0                                               # phosphorus (mmol P/m^2)
187      pen_depth_c: 0.2                                     # penetration depth of carbon (m)
188      pen_depth_n: 0.25                                    # penetration depth of nitrogen (m)
189      pen_depth_p: 0.313                                   # penetration depth of phosphorus (m)
190    coupling:
191      burial_target: Q17
192  Q17:
193    long_name: benthic buried matter
194    model: ersem/benthic_base
195    parameters:
196      composition: cnp                                     # elemental composition
197    initialization:
198      c: 0.0                                               # carbon (mg C/m^2)
199      n: 0.0                                               # nitrogen (mmol N/m^2)
200      p: 0.0                                               # phosphorus (mmol P/m^2)
201  B1:
202    long_name: bacteria
203    model: ersem/bacteria_docdyn
204    parameters:
205      iswBlim: 2                                           # nutrient limitation (1: minimum of inorganic and organic availability, 2: additive availability)
206      q10: 2.0                                             # Q_10 temperature coefficient (-)
207      chdo: 0.31                                           # Michaelis-Menten constant for oxygen limitation (-)
208      chn: 0.5                                             # Michaelis-Menten constant for nitrate limitation (mmol N/m^3)
209      chp: 0.1                                             # Michaelis-Menten constant for phosphate limitation (mmol P/m^3)
210      sd: 0.05                                             # specific mortality at reference temperature (1/d)
211      sum: 2.2                                             # maximum specific uptake at reference temperature (1/d)
212      pu: 0.6                                              # efficiency at high oxygen levels (-)
213      puo: 0.2                                             # efficiency at low oxygen levels (-)
214      srs: 0.1                                             # specific rest respiration at reference temperature (1/d)
215      sR1: 1.0                                             # maximum turn-over rate of DOM (1/d), default = 1.0
216      qpc: 0.0019                                          # maximum phosphorus to carbon ratio (mmol P/mg C)
217      qnc: 0.0167                                          # maximum nitrogen to carbon ratio (mmol N/mg C)
218      ur_O2: 0.1                                           # oxygen consumed per carbon respired (mmol O_2/mg C)
219      sR1N1: 0.0                                           # mineralisation rate of labile dissolved organic phosphorus (1/d)
220      sR1N4: 0.0                                           # mineralisation rate of labile dissolved organic nitrogen (1/d)
221      fsink: 7e-05                                         # scavenging rate for iron (1/d)
222      c0: 0.01                                             # background carbon concentration (mg C/m^3)
223      nRP: 3                                               # number of substrates, default = 0
224      sRP1R1: 0.01                                         # remineralisation of substrate 1 to DOM (1/d)
225      sRP2R1: 0.0025                                       # remineralisation of substrate 2 to DOM (1/d)
226      sRP3R1: 0.001                                        # remineralisation of substrate 3 to DOM (1/d)
227      rR2: 0.0075                                          # fraction of semi-labile DOC available to bacteria (-)
228      rR3: 0.0025                                          # fraction of semi-refractory DOC available to bacteria (-)
229      frR3: 0.3                                            # fraction of activity respiration converted to semi-refractory DOC (-)
230    initialization:
231      c: 15.7                                              # carbon (mg C/m^3)
232      n: 0.26                                              # nitrogen (mmol N/m^3)
233      p: 0.029                                             # phosphorus (mmol P/m^3)
234    coupling:
235      RP1: R4
236      RP2: R6
237      RP3: R8
238      N1p: N1/p                                            # phosphate (mmol P/m^3)
239      N4n: N4/n                                            # ammonium (mmol N/m^3)
240      R1c: R1/c                                            # labile dissolved organic carbon (mg C/m^3)
241      R1p: R1/p                                            # labile dissolved organic phosphorus (mmol P/m^3)
242      R1n: R1/n                                            # labile dissolved organic nitrogen (mmol N/m^3)
243      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
244      R3c: R3/c                                            # semi-refractory DOC (mg C/m^3)
245      O3c: O3/c                                            # carbon dioxide (mmol C/m^3)
246      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
247  pel_nit:
248    long_name: pelagic nitrification
249    model: ersem/nitrification
250    parameters:
251      q10: 2.0                                             # Q_10 temperature coefficient (-)
252      ISWph: 1                                             # pH impact on nitrification (0: off, 1: on)
253      sN4N3: 0.5                                           # specific nitrification rate (1/d)
254      chN3o: 2700.0                                        # Michaelis-Menten constant for cubic oxygen dependence of nitrification ((mmol O_2/m^3)^3)
255      chN4n: 0.5                                           # Michaelis-Menten constant for cubic ammonium dependence of nitrification ((mmol N/m^3)^3), default = 0.0
256    coupling:
257      N3n: N3/n                                            # nitrate (mmol N/m^3)
258      N4n: N4/n                                            # ammonium (mmol N/m^3)
259      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
260  P1:
261    long_name: diatoms
262    model: ersem/primary_producer
263    parameters:
264      sum: 1.375                                           # maximum specific productivity at reference temperature (1/d)
265      q10: 2.0                                             # Q_10 temperature coefficient (-)
266      srs: 0.04                                            # specific rest respiration at reference temperature (1/d)
267      pu_ea: 0.2                                           # excreted fraction of primary production (-)
268      pu_ra: 0.2                                           # respired fraction of primary production (-)
269      qnlc: 0.0042                                         # minimum nitrogen to carbon ratio (mmol N/mg C)
270      qplc: 0.0001                                         # minimum phosphorus to carbon ratio (mmol P/mg C)
271      xqcp: 1.0                                            # threshold for phosphorus limitation (relative to Redfield ratio) (-)
272      xqcn: 1.0                                            # threshold for nitrogen limitation (relative to Redfield ratio) (-)
273      xqp: 2.0                                             # maximum phosphorus to carbon ratio (relative to Redfield ratio) (-)
274      xqn: 1.075                                           # maximum nitrogen to carbon ratio (relative to Redfield ratio) (-)
275      qun3: 0.0025                                         # nitrate affinity (m^3/mg C/d)
276      qun4: 0.0025                                         # ammonium affinity (m^3/mg C/d)
277      qurp: 0.003                                          # phosphate affinity (m^3/mg C/d)
278      snplux: 1.0                                          # specific tendency of luxury uptake of nutrients towards maximum quota (1/d), default = 1.0
279      use_Si: true                                         # use silicate, default = false
280      qsc: 0.0118                                          # maximum silicate to carbon ratio (mmol Si/mg C)
281      chs: 0.2                                             # Michaelis-Menten constant for silicate limitation (mmol/m^3)
282      sdo: 0.05                                            # 1.1 of minimal specific lysis rate (1/d)
283      alpha: 4.0                                           # initial slope of PI-curve (mg C m^2/mg Chl/W/d)
284      beta: 0.07                                           # photoinhibition parameter (mg C m^2/mg Chl/W/d)
285      phim: 0.06                                           # maximum effective chlorophyll to carbon photosynthesis ratio (mg Chl/mg C)
286      Limnut: 1                                            # nitrogen-phosphorus colimitation formulation (0: geometric mean, 1: minimum, 2: harmonic mean)
287      docdyn: true                                         # use dynamic ratio of labile to semi-labile DOM production, default = false
288      uB1c_O2: 0.11                                        # oxygen produced per unit of carbon fixed (mmol O_2/mg C)
289      urB1_O2: 0.1                                         # oxygen consumed per unit of carbon respired (mmol O_2/mg C)
290      iopABS: 0.007                                        # specific shortwave absorption (m^2/mg Chl), default = 0.008
291      iopBBS: 0.00048                                      # specific shortwave backscatter (m^2/mg Chl), default = 0.003
292      c0: 0.0001                                           # background carbon concentration (mg C/m^3), default = 0.0
293      resm: 5.0                                            # maximum nutrient-limitation-induced sinking velocity (m/d), default = 0.0
294      esni: 0.7                                            # level of nutrient limitation below which sinking commences (-)
295      cenh: false                                          # enable atmospheric CO2 influence on photosynthesis, default = false
296      ndeposition: 3                                       # number of target pools for sedimentation, default = 1
297      qxc2: 0.2                                            # fraction of carbon sinking into deposition target 2 (-)
298      qxc3: 0.005                                          # fraction of carbon sinking into deposition target 3 (-)
299      qxn2: 0.2                                            # fraction of nitrogen sinking into deposition target 2 (-)
300      qxn3: 0.005                                          # fraction of nitrogen sinking into deposition target 3 (-)
301      qxp2: 0.24                                           # fraction of phosphorus sinking into deposition target 2 (-)
302      qxp3: 0.003                                          # fraction of phosphorus sinking into deposition target 3 (-)
303      qxs2: 0.0                                            # fraction of silicate sinking into deposition target 2 (-)
304      qxs3: 0.0                                            # fraction of silicate sinking into deposition target 3 (-)
305    initialization:
306      c: 8.0                                               # carbon (mg C/m^3)
307      n: 0.1114                                            # nitrogen (mmol N/m^3)
308      p: 0.009                                             # phosphorus (mmol P/m^3)
309      s: 0.128                                             # silicate (mmol Si/m^3)
310      Chl: 0.4                                             # chlorophyll a (mg/m^3)
311    coupling:
312      RP: R6
313      deposition_target1: Q6/surface
314      deposition_target2: Q1
315      deposition_target3: Q7/surface
316      N1p: N1/p                                            # phosphate (mmol P/m^3)
317      N3n: N3/n                                            # nitrate (mmol N/m^3)
318      N4n: N4/n                                            # ammonium (mmol N/m^3)
319      N5s: N5/s                                            # silicate (mmol Si/m^3)
320      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
321      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
322      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
323      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
324      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
325      O3c: O3/c                                            # carbon dioxide (mmol C/m^3)
326  P2:
327    long_name: nanophytoplankton
328    model: ersem/primary_producer
329    parameters:
330      sum: 1.625                                           # maximum specific productivity at reference temperature (1/d)
331      q10: 2.0                                             # Q_10 temperature coefficient (-)
332      srs: 0.04                                            # specific rest respiration at reference temperature (1/d)
333      pu_ea: 0.2                                           # excreted fraction of primary production (-)
334      pu_ra: 0.2                                           # respired fraction of primary production (-)
335      qnlc: 0.005                                          # minimum nitrogen to carbon ratio (mmol N/mg C)
336      qplc: 0.000225                                       # minimum phosphorus to carbon ratio (mmol P/mg C)
337      xqcp: 1.0                                            # threshold for phosphorus limitation (relative to Redfield ratio) (-)
338      xqcn: 1.0                                            # threshold for nitrogen limitation (relative to Redfield ratio) (-)
339      xqp: 2.0                                             # maximum phosphorus to carbon ratio (relative to Redfield ratio) (-)
340      xqn: 1.075                                           # maximum nitrogen to carbon ratio (relative to Redfield ratio) (-)
341      qun3: 0.004                                          # nitrate affinity (m^3/mg C/d)
342      qun4: 0.004                                          # ammonium affinity (m^3/mg C/d)
343      qurp: 0.004                                          # phosphate affinity (m^3/mg C/d)
344      snplux: 1.0                                          # specific tendency of luxury uptake of nutrients towards maximum quota (1/d), default = 1.0
345      use_Si: false                                        # use silicate, default = false
346      sdo: 0.05                                            # 1.1 of minimal specific lysis rate (1/d)
347      alpha: 5.0                                           # initial slope of PI-curve (mg C m^2/mg Chl/W/d)
348      beta: 0.1                                            # photoinhibition parameter (mg C m^2/mg Chl/W/d)
349      phim: 0.025                                          # maximum effective chlorophyll to carbon photosynthesis ratio (mg Chl/mg C)
350      Limnut: 1                                            # nitrogen-phosphorus colimitation formulation (0: geometric mean, 1: minimum, 2: harmonic mean)
351      docdyn: true                                         # use dynamic ratio of labile to semi-labile DOM production, default = false
352      uB1c_O2: 0.11                                        # oxygen produced per unit of carbon fixed (mmol O_2/mg C)
353      urB1_O2: 0.1                                         # oxygen consumed per unit of carbon respired (mmol O_2/mg C)
354      iopABS: 0.0041                                       # specific shortwave absorption (m^2/mg Chl), default = 0.008
355      iopBBS: 0.003                                        # specific shortwave backscatter (m^2/mg Chl), default = 0.003
356      c0: 0.001                                            # background carbon concentration (mg C/m^3), default = 0.0
357      calcify: true                                        # calcify, default = false
358      resm: 0.0                                            # maximum nutrient-limitation-induced sinking velocity (m/d), default = 0.0
359      esni: 0.7                                            # level of nutrient limitation below which sinking commences (-)
360      cenh: false                                          # enable atmospheric CO2 influence on photosynthesis, default = false
361      ndeposition: 3                                       # number of target pools for sedimentation, default = 1
362      qxc2: 0.5                                            # fraction of carbon sinking into deposition target 2 (-)
363      qxc3: 0.05                                           # fraction of carbon sinking into deposition target 3 (-)
364      qxn2: 0.5                                            # fraction of nitrogen sinking into deposition target 2 (-)
365      qxn3: 0.05                                           # fraction of nitrogen sinking into deposition target 3 (-)
366      qxp2: 0.6                                            # fraction of phosphorus sinking into deposition target 2 (-)
367      qxp3: 0.03                                           # fraction of phosphorus sinking into deposition target 3 (-)
368    initialization:
369      c: 5.9                                               # carbon (mg C/m^3)
370      n: 0.0926                                            # nitrogen (mmol N/m^3)
371      p: 0.0036                                            # phosphorus (mmol P/m^3)
372      Chl: 0.3                                             # chlorophyll a (mg/m^3)
373    coupling:
374      RP: R4
375      deposition_target1: Q6/surface
376      deposition_target2: Q1
377      deposition_target3: Q7/surface
378      N1p: N1/p                                            # phosphate (mmol P/m^3)
379      N3n: N3/n                                            # nitrate (mmol N/m^3)
380      N4n: N4/n                                            # ammonium (mmol N/m^3)
381      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
382      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
383      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
384      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
385      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
386      O3c: O3/c                                            # carbon dioxide (mmol C/m^3)
387      RainR: L2/RainR                                      # rain ratio (PIC : POC) (-)
388      L2c: L2/c                                            # free calcite (mg C/m^3)
389  P3:
390    long_name: picophytoplankton
391    model: ersem/primary_producer
392    parameters:
393      sum: 2.0                                             # maximum specific productivity at reference temperature (1/d)
394      q10: 2.0                                             # Q_10 temperature coefficient (-)
395      srs: 0.045                                           # specific rest respiration at reference temperature (1/d)
396      pu_ea: 0.2                                           # excreted fraction of primary production (-)
397      pu_ra: 0.2                                           # respired fraction of primary production (-)
398      qnlc: 0.006                                          # minimum nitrogen to carbon ratio (mmol N/mg C)
399      qplc: 0.00035                                        # minimum phosphorus to carbon ratio (mmol P/mg C)
400      xqcp: 1.0                                            # threshold for phosphorus limitation (relative to Redfield ratio) (-)
401      xqcn: 1.0                                            # threshold for nitrogen limitation (relative to Redfield ratio) (-)
402      xqp: 1.5                                             # maximum phosphorus to carbon ratio (relative to Redfield ratio) (-)
403      xqn: 1.05                                            # maximum nitrogen to carbon ratio (relative to Redfield ratio) (-)
404      qun3: 0.006                                          # nitrate affinity (m^3/mg C/d)
405      qun4: 0.007                                          # ammonium affinity (m^3/mg C/d)
406      qurp: 0.006                                          # phosphate affinity (m^3/mg C/d)
407      snplux: 1.0                                          # specific tendency of luxuary uptake of nutrients towards maximum quoata (1/d)
408      use_Si: false                                        # use silicate, default = false
409      sdo: 0.055                                           # 1.1 of minimal specific lysis rate (1/d)
410      alpha: 6.0                                           # initial slope of PI-curve (mg C m^2/mg Chl/W/d)
411      beta: 0.12                                           # photoinhibition parameter (mg C m^2/mg Chl/W/d)
412      phim: 0.015                                          # maximum effective chlorophyll to carbon photosynthesis ratio (mg Chl/mg C)
413      Limnut: 1                                            # nitrogen-phosphorus colimitation formulation (0: geometric mean, 1: minimum, 2: harmonic mean)
414      docdyn: true                                         # use dynamic ratio of labile to semi-labile DOM production, default = false
415      uB1c_O2: 0.11                                        # oxygen produced per unit of carbon fixed (mmol O_2/mg C)
416      urB1_O2: 0.1                                         # oxygen consumed per unit of carbon respired (mmol O_2/mg C)
417      iopABS: 0.023                                        # specific shortwave absorption (m^2/mg Chl), default = 0.008
418      iopBBS: 0.003                                        # specific shortwave backscatter (m^2/mg Chl), default = 0.003
419      c0: 0.0088                                           # background carbon concentration (mg C/m^3), default = 0.0
420      resm: 0.0                                            # maximum nutrient-limitation-induced sinking velocity (m/d), default = 0.0
421      esni: 0.7                                            # level of nutrient limitation below which sinking commences (-)
422      cenh: false                                          # enable atmospheric CO2 influence on photosynthesis, default = false
423      ndeposition: 3                                       # number of target pools for sedimentation, default = 1
424      qxc2: 0.5                                            # fraction of carbon sinking into deposition target 2 (-)
425      qxc3: 0.05                                           # fraction of carbon sinking into deposition target 3 (-)
426      qxn2: 0.5                                            # fraction of nitrogen sinking into deposition target 2 (-)
427      qxn3: 0.05                                           # fraction of nitrogen sinking into deposition target 3 (-)
428      qxp2: 0.6                                            # fraction of phosphorus sinking into deposition target 2 (-)
429      qxp3: 0.03                                           # fraction of phosphorus sinking into deposition target 3 (-)
430    initialization:
431      c: 5.9                                               # carbon (mg C/m^3)
432      n: 0.0926                                            # nitrogen (mmol N/m^3)
433      p: 0.0036                                            # phosphorus (mmol P/m^3)
434      Chl: 0.3                                             # chlorophyll a (mg/m^3)
435    coupling:
436      RP: R4
437      deposition_target1: Q6/surface
438      deposition_target2: Q1
439      deposition_target3: Q7/surface
440      N1p: N1/p                                            # phosphate (mmol P/m^3)
441      N3n: N3/n                                            # nitrate (mmol N/m^3)
442      N4n: N4/n                                            # ammonium (mmol N/m^3)
443      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
444      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
445      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
446      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
447      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
448      O3c: O3/c                                            # carbon dioxide (mmol C/m^3)
449  P4:
450    long_name: microphytoplankton
451    model: ersem/primary_producer
452    parameters:
453      sum: 1.125                                           # maximum specific productivity at reference temperature (1/d)
454      q10: 2.0                                             # Q_10 temperature coefficient (-)
455      srs: 0.035                                           # specific rest respiration at reference temperature (1/d)
456      pu_ea: 0.2                                           # excreted fraction of primary production (-)
457      pu_ra: 0.2                                           # respired fraction of primary production (-)
458      qnlc: 0.0042                                         # minimum nitrogen to carbon ratio (mmol N/mg C)
459      qplc: 0.0001                                         # minimum phosphorus to carbon ratio (mmol P/mg C)
460      xqcp: 1.0                                            # threshold for phosphorus limitation (relative to Redfield ratio) (-)
461      xqcn: 1.0                                            # threshold for nitrogen limitation (relative to Redfield ratio) (-)
462      xqp: 2.7                                             # maximum phosphorus to carbon ratio (relative to Redfield ratio) (-)
463      xqn: 1.1                                             # maximum nitrogen to carbon ratio (relative to Redfield ratio) (-)
464      qun3: 0.002                                          # nitrate affinity (m^3/mg C/d)
465      qun4: 0.002                                          # ammonium affinity (m^3/mg C/d)
466      qurp: 0.002                                          # phosphate affinity (m^3/mg C/d)
467      snplux: 1.0                                          # specific tendency of luxuary uptake of nutrients towards maximum quoata (1/d)
468      use_Si: false                                        # use silicate, default = false
469      sdo: 0.045                                           # 1.1 of minimal specific lysis rate (1/d)
470      alpha: 3.0                                           # initial slope of PI-curve (mg C m^2/mg Chl/W/d)
471      beta: 0.06                                           # photoinhibition parameter (mg C m^2/mg Chl/W/d)
472      phim: 0.045                                          # maximum effective chlorophyll to carbon photosynthesis ratio (mg Chl/mg C)
473      Limnut: 1                                            # nitrogen-phosphorus colimitation formulation (0: geometric mean, 1: minimum, 2: harmonic mean)
474      docdyn: true                                         # use dynamic ratio of labile to semi-labile DOM production, default = false
475      uB1c_O2: 0.11                                        # oxygen produced per unit of carbon fixed (mmol O_2/mg C)
476      urB1_O2: 0.1                                         # oxygen consumed per unit of carbon respired (mmol O_2/mg C)
477      iopABS: 0.008                                        # specific shortwave absorption (m^2/mg Chl), default = 0.008
478      iopBBS: 0.00048                                      # specific shortwave backscatter (m^2/mg Chl), default = 0.003
479      c0: 0.0001                                           # background carbon concentration (mg C/m^3), default = 0.0
480      resm: 5.0                                            # maximum nutrient-limitation-induced sinking velocity (m/d), default = 0.0
481      esni: 0.7                                            # level of nutrient limitation below which sinking commences (-)
482      cenh: false                                          # enable atmospheric CO2 influence on photosynthesis, default = false
483      ndeposition: 3                                       # number of target pools for sedimentation, default = 1
484      qxc2: 0.5                                            # fraction of carbon sinking into deposition target 2 (-)
485      qxc3: 0.05                                           # fraction of carbon sinking into deposition target 3 (-)
486      qxn2: 0.5                                            # fraction of nitrogen sinking into deposition target 2 (-)
487      qxn3: 0.05                                           # fraction of nitrogen sinking into deposition target 3 (-)
488      qxp2: 0.6                                            # fraction of phosphorus sinking into deposition target 2 (-)
489      qxp3: 0.03                                           # fraction of phosphorus sinking into deposition target 3 (-)
490    initialization:
491      c: 5.9                                               # carbon (mg C/m^3)
492      n: 0.0926                                            # nitrogen (mmol N/m^3)
493      p: 0.0036                                            # phosphorus (mmol P/m^3)
494      Chl: 0.3                                             # chlorophyll a (mg/m^3)
495    coupling:
496      RP: R6
497      deposition_target1: Q6/surface
498      deposition_target2: Q1
499      deposition_target3: Q7/surface
500      N1p: N1/p                                            # phosphate (mmol P/m^3)
501      N3n: N3/n                                            # nitrate (mmol N/m^3)
502      N4n: N4/n                                            # ammonium (mmol N/m^3)
503      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
504      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
505      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
506      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
507      O2o: O2/o                                            # oxygen (mmol O_2/m^3)
508      O3c: O3/c                                            # carbon dioxide (mmol C/m^3)
509  Z4:
510    long_name: mesozooplankton
511    model: ersem/mesozooplankton
512    parameters:
513      q10: 2.0                                             # Q_10 temperature coefficient (-)
514      minfood: 12.0                                        # Michaelis-Menten constant to perceive food (mg C/m^3)
515      chuc: 36.0                                           # Michaelis-Menten constant for food uptake (mg C/m^3)
516      sum: 1.0                                             # maximum specific uptake at reference temperature (1/d)
517      pu: 0.6                                              # assimilation efficiency (-)
518      pu_ea: 0.5                                           # fraction of unassimilated prey that is excreted (not respired) (-)
519      pu_eaR: 0.9                                          # fraction of unassimilated detritus that is excreted (not respired) (-)
520      pe_R1: 0.5                                           # dissolved fraction of excreted/dying matter (-)
521      srs: 0.015                                           # specific rest respiration at reference temperature (1/d)
522      sd: 0.05                                             # basal mortality (1/d)
523      sdo: 0.2                                             # maximum mortality due to oxygen limitation (1/d)
524      chro: 7.81                                           # Michaelis-Menten constant for oxygen limitation (-)
525      qpc: 0.000786                                        # phosphorus to carbon ratio (mmol P/mg C)
526      qnc: 0.0126                                          # nitrogen to carbon ratio (mmol N/mg C)
527      Minprey: 300.0                                       # food threshold for overwintering state (mg C/m^2)
528      repw: 0.0025                                         # specific overwintering respiration (1/d)
529      mort: 0.0025                                         # specific overwintering mortality (1/d)
530      R1R2: 1.0                                            # labile fraction of produced dissolved organic carbon (-)
531      xR1p: 1.2                                            # transfer of phosphorus to DOM, relative to POM (-), default = 0.0
532      xR1n: 1.0                                            # transfer of nitrogen to DOM, relative to POM (-), default = 0.0
533      urB1_O2: 0.1                                         # oxygen consumed per carbon respired (mmol O_2/mg C)
534      gutdiss: 0.5                                         # fraction of prey calcite that dissolves after ingestion (-)
535      c0: 0.0033                                           # background concentration (mg C/m^3)
536      nprey: 9                                             # number of prey types, default = 0
537      suprey1: 0.15                                        # relative affinity for prey type 1 (-)
538      suprey2: 0.05                                        # relative affinity for prey type 2 (-)
539      suprey3: 0.0                                         # relative affinity for prey type 3 (-)
540      suprey4: 0.15                                        # relative affinity for prey type 4 (-)
541      suprey5: 0.25                                        # relative affinity for prey type 5 (-)
542      suprey6: 0.25                                        # relative affinity for prey type 6 (-)
543      suprey7: 0.05                                        # relative affinity for prey type 7 (-)
544      suprey8: 0.0                                         # relative affinity for prey type 8 (-)
545      suprey9: 0.1                                         # relative affinity for prey type 9 (-)
546      prey9ispom: true                                     # prey type 9 is detritus, default = false
547    initialization:
548      c: 1.2                                               # carbon (mg C/m^3)
549    coupling:
550      prey1: P1
551      prey2: P2
552      prey3: P3
553      prey4: P4
554      prey5: Z4
555      prey6: Z5
556      prey7: Z6
557      prey8: B1
558      prey9: R6
559      RP: R8
560      N1p: N1/p                                            # phosphate (mmol P/m^3)
561      N4n: N4/n                                            # ammonium (mmol N/m^3)
562      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
563      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
564      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
565      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
566      O2o: O2/o                                            # oxygen source (mmol O_2/m^3)
567      O3c: O3/c                                            # carbon dioxide sink (mmol C/m^3)
568      L2c: L2/c                                            # calcite (mg C/m^3)
569  Z5:
570    long_name: microzooplankton
571    model: ersem/microzooplankton
572    parameters:
573      q10: 2.0                                             # Q_10 temperature coefficient (-)
574      minfood: 12.0                                        # Michaelis-Menten constant to perceive food (mg C/m^3)
575      chuc: 32.0                                           # Michaelis-Menten constant for food uptake (mg C/m^3)
576      sum: 1.25                                            # maximum specific uptake at reference temperature (1/d)
577      pu: 0.5                                              # assimilation efficiency (-)
578      pu_ea: 0.5                                           # fraction of unassimilated prey that is excreted (not respired) (-)
579      pe_R1: 0.5                                           # dissolved fraction of excreted/dying matter (-)
580      srs: 0.02                                            # specific rest respiration at reference temperature (1/d)
581      sd: 0.05                                             # basal mortality (1/d)
582      sdo: 0.25                                            # maximum mortality due to oxygen limitation (1/d)
583      chro: 7.81                                           # Michaelis-Menten constant for oxygen limitation (-)
584      qpc: 0.001                                           # maximum phosphorus to carbon ratio (mmol P/mg C)
585      qnc: 0.0167                                          # maximum nitrogen to carbon ratio (mmol N/mg C)
586      stempp: 0.5                                          # specific excretion rate of excess phosphorus (1/d)
587      stempn: 0.5                                          # specific excretion rate of excess nitrogen (1/d)
588      R1R2: 1.0                                            # labile fraction of produced dissolved organic carbon (-)
589      xR1p: 1.2                                            # transfer of phosphorus to DOM, relative to POM (-), default = 1.0
590      xR1n: 1.0                                            # transfer of nitrogen to DOM, relative to POM (-), default = 1.0
591      urB1_O2: 0.1                                         # oxygen consumed per carbon respired (mmol O_2/mg C)
592      gutdiss: 0.5                                         # fraction of prey calcite that dissolves after ingestion (-)
593      c0: 0.0033                                           # background concentration (mg C/m^3)
594      nprey: 7                                             # number of prey types, default = 0
595      suprey1: 0.1                                         # relative affinity for prey type 1 (-)
596      suprey2: 0.15                                        # relative affinity for prey type 2 (-)
597      suprey3: 0.15                                        # relative affinity for prey type 3 (-)
598      suprey4: 0.15                                        # relative affinity for prey type 4 (-)
599      suprey5: 0.1                                         # relative affinity for prey type 5 (-)
600      suprey6: 0.15                                        # relative affinity for prey type 6 (-)
601      suprey7: 0.2                                         # relative affinity for prey type 7 (-)
602    initialization:
603      c: 7.2                                               # carbon (mg C/m^3)
604      n: 0.12                                              # nitrogen (mmol N/m^3)
605      p: 0.0113                                            # phosphorus (mmol P/m^3)
606    coupling:
607      prey1: B1
608      prey2: P1
609      prey3: P2
610      prey4: P3
611      prey5: P4
612      prey6: Z5
613      prey7: Z6
614      RP: R6
615      N1p: N1/p                                            # phosphate (mmol P/m^3)
616      N4n: N4/n                                            # ammonium (mmol N/m^3)
617      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
618      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
619      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
620      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
621      O2o: O2/o                                            # oxygen source (mmol O_2/m^3)
622      O3c: O3/c                                            # carbon dioxide sink (mmol C/m^3)
623      L2c: L2/c                                            # calcite (mg C/m^3)
624  Z6:
625    long_name: nanoflagellates
626    model: ersem/microzooplankton
627    parameters:
628      q10: 2.0                                             # Q_10 temperature coefficient (-)
629      minfood: 12.0                                        # Michaelis-Menten constant to perceive food (mg C/m^3)
630      chuc: 28.0                                           # Michaelis-Menten constant for food uptake (mg C/m^3)
631      sum: 1.5                                             # maximum specific uptake at reference temperature (1/d)
632      pu: 0.4                                              # assimilation efficiency (-)
633      pu_ea: 0.5                                           # fraction of unassimilated prey that is excreted (not respired) (-)
634      pe_R1: 0.5                                           # dissolved fraction of excreted/dying matter (-)
635      srs: 0.025                                           # specific rest respiration at reference temperature (1/d)
636      sd: 0.05                                             # basal mortality (1/d)
637      sdo: 0.3                                             # maximum mortality due to oxygen limitation (1/d)
638      chro: 7.81                                           # Michaelis-Menten constant for oxygen limitation (-)
639      qpc: 0.001                                           # maximum phosphorus to carbon ratio (mmol P/mg C)
640      qnc: 0.0167                                          # maximum nitrogen to carbon ratio (mmol N/mg C)
641      stempp: 0.5                                          # specific excretion rate of excess phosphorus (1/d)
642      stempn: 0.5                                          # specific excretion rate of excess nitrogen (1/d)
643      R1R2: 1.0                                            # labile fraction of produced dissolved organic carbon (-)
644      xR1p: 1.2                                            # transfer of phosphorus to DOM, relative to POM (-), default = 1.0
645      xR1n: 1.0                                            # transfer of nitrogen to DOM, relative to POM (-), default = 1.0
646      urB1_O2: 0.1                                         # oxygen consumed per carbon respired (mmol O_2/mg C)
647      gutdiss: 0.5                                         # fraction of prey calcite that dissolves after ingestion (-)
648      c0: 0.0033                                           # background concentration (mg C/m^3)
649      nprey: 4                                             # number of prey types, default = 0
650      suprey1: 0.45                                        # relative affinity for prey type 1 (-)
651      suprey2: 0.15                                        # relative affinity for prey type 2 (-)
652      suprey3: 0.25                                        # relative affinity for prey type 3 (-)
653      suprey4: 0.15                                        # relative affinity for prey type 4 (-)
654    initialization:
655      c: 2.421                                             # carbon (mg C/m^3)
656      n: 0.0505                                            # nitrogen (mmol N/m^3)
657      p: 0.047                                             # phosphorus (mmol P/m^3)
658    coupling:
659      prey1: B1
660      prey2: P2
661      prey3: P3
662      prey4: Z6
663      RP: R4
664      N1p: N1/p                                            # phosphate (mmol P/m^3)
665      N4n: N4/n                                            # ammonium (mmol N/m^3)
666      R1c: R1/c                                            # dissolved organic carbon (mg C/m^3)
667      R1p: R1/p                                            # dissolved organic phosphorus (mmol P/m^3)
668      R1n: R1/n                                            # dissolved organic nitrogen (mmol N/m^3)
669      R2c: R2/c                                            # semi-labile dissolved organic carbon (mg C/m^3)
670      RPs: R6/s                                            # particulate organic silicate (mmol Si/m^3)
671      O2o: O2/o                                            # oxygen source (mmol O_2/m^3)
672      O3c: O3/c                                            # carbon dioxide sink (mmol C/m^3)
673      L2c: L2/c                                            # calcite (mg C/m^3)
674  L2:
675    long_name: calcite
676    model: ersem/calcification
677    parameters:
678      iswcal: 1                                            # calcification/dissolution dependence on calcite saturation (1: power law, 2: hyperbolic)
679      ncalc: 0.81                                          # power of the calcification law (Ridgwell et al. 2007, mineral calcite) (-)
680      ndiss: 2.22                                          # power of the dissolution law (Keir 1980) (-)
681      Rain0: 0.6                                           # maximum rain ratio from PISCES (-)
682      sedL2: 10.0                                          # sinking velocity (m/d)
683      sL2O3: 0.03                                          # maximum specific dissolution rate (1/d), default = 1.0
684      c0: 1e-05                                            # background concentration (mg C/m^3), default = 0.0
685      ndeposition: 1                                       # number of target pools for sedimentation, default = 1
686    initialization:
687      c: 0.05                                              # carbon (mg C/m^3)
688    coupling:
689      deposition_target1: bL2
690      om_cal: O3/Om_cal                                    # calcite saturation (-)
691      O3c: O3/c                                            # total dissolved inorganic carbon (mmol C/m^3)
692  bL2:
693    long_name: benthic calcite
694    model: ersem/benthic_calcite
695    parameters:
696      iswcal: 1                                            # dissolution dependence on calcite saturation (0: none, 1: power law, 2: hyperbolic)
697      ndiss: 2.22                                          # power of the dissolution law (Keir 1980) (-)
698      fdissmax: 30.0                                       # maximum specific dissolution rate (1/d), default = 0.0
699      fdissmin: 0.05                                       # minimum specific dissolution rate (1/d), default = 0.03
700      c0: 1e-05                                            # background calcite concentration (mg C/m^2), default = 0.0
701    initialization:
702      c: 0.05                                              # carbon (mg C/m^2)
703    coupling:
704      Om_Cal: O3/Om_cal                                    # calcite saturation (-)
705      O3c: O3/c                                            # dissolved inorganic carbon (mmol/m^3)
706  ben_col:
707    long_name: benthic column
708    model: ersem/benthic_column
709    parameters:
710      qPW: 0.4                                             # sediment porosity (-)
711      EDZ_mix: 20.0                                        # equilibrium diffusive speed between sediment surface water (d/m)
712      d_tot: 0.3                                           # depth of sediment column (m)
713      Etur: 2e-06                                          # basal bioturbation rate (m^2/d)
714      mtur: 10.0                                           # maximum relative turbation enhancement (-)
715      htur: 116.0                                          # Michaelis-Menten constant for bioturbation (mg C/m^2/d)
716      dtur: 0.02                                           # bioturbation depth (m)
717      EDZ_1: 5e-05                                         # diffusivity in oxygenated layer (m^2/d)
718      EDZ_2: 5e-05                                         # diffusivity in oxidized layer (m^2/d)
719      EDZ_3: 5e-05                                         # diffusivity in anoxic layer (m^2/d)
720      irr_min: 2.0                                         # minimum diffusion enhancement through bioirrigation (-)
721      mirr: 10.0                                           # maximum relative diffusion enhancement due to bioirrigation (-)
722      hirr: 101.0                                          # Michaelis-Menten constant for bioirrigation (mg C/m^2/d)
723    initialization:
724      D1m: 0.009                                           # depth of bottom interface of oxygenated layer (m)
725      D2m: 0.05                                            # depth of bottom interface of oxidized layer (m)
726  K1:
727    long_name: benthic phosphate
728    model: ersem/benthic_column_dissolved_matter
729    parameters:
730      composition: p                                       # composition (any combination of c,n,p,s,o,a)
731      ads1: 100.0                                          # adsorption in layer 1 (total:dissolved) (-), default = 1.0
732      ads2: 100.0                                          # adsorption in layer 2 (total:dissolved) (-), default = 1.0
733      ads3: 2.0                                            # adsorption in layer 3 (total:dissolved) (-), default = 1.0
734      c0: 0.0003                                           # background concentration (mg C/m^2)
735    initialization:
736      p: 6.8                                               # phosphorus (mmol/m^2)
737    coupling:
738      p_pel: N1/p                                          # pelagic phosphorus (mmol/m^3)
739  K3:
740    long_name: benthic nitrate
741    model: ersem/benthic_column_dissolved_matter
742    parameters:
743      composition: n                                       # composition (any combination of c,n,p,s,o,a)
744      last_layer: 2                                        # sediment layer where concentration drops to zero, default = 3
745      relax: 5.0                                           # rate of relaxation towards equilibrium concentration profile (1/d)
746      minD: 0.0001                                         # minimum depth of zero-concentration isocline (m)
747      c0: 0.0003                                           # background concentration (mg C/m^2)
748    initialization:
749      n: 0.584                                             # nitrogen (mmol/m^2)
750    coupling:
751      n_pel: N3/n                                          # pelagic nitrogen (mmol/m^3)
752  K4:
753    long_name: benthic ammonium
754    model: ersem/benthic_column_dissolved_matter
755    parameters:
756      composition: n                                       # composition (any combination of c,n,p,s,o,a)
757      ads1: 3.0                                            # adsorption in layer 1 (total:dissolved) (-), default = 1.0
758      ads2: 3.0                                            # adsorption in layer 2 (total:dissolved) (-), default = 1.0
759      ads3: 3.0                                            # adsorption in layer 3 (total:dissolved) (-), default = 1.0
760      correction: true                                     # move losses in oxygenic layer to deeper layers if pelagic concentration is limiting, default = false
761      c0: 0.0003                                           # background concentration (mg C/m^2)
762    initialization:
763      n: 25.0                                              # nitrogen (mmol/m^2)
764    coupling:
765      n_pel: N4/n                                          # pelagic nitrogen (mmol/m^3)
766  K5:
767    long_name: benthic silicate
768    model: ersem/benthic_column_dissolved_matter
769    parameters:
770      composition: s                                       # composition (any combination of c,n,p,s,o,a)
771      c0: 9.e-05                                           # background concentration (mg C/m^2)
772    initialization:
773      s: 20.1                                              # silicate (mmol/m^2)
774    coupling:
775      s_pel: N5/s                                          # pelagic silicate (mmol/m^3)
776  G2:
777    long_name: benthic oxygen
778    model: ersem/benthic_column_dissolved_matter
779    parameters:
780      composition: o                                       # composition (any combination of c,n,p,s,o,a)
781      last_layer: 1                                        # sediment layer where concentration drops to zero, default = 3
782      relax: 5.0                                           # rate of relaxation towards equilibrium concentration profile (1/d)
783      minD: 0.0001                                         # minimum depth of zero-concentration isocline (m)
784    initialization:
785      o: 0.0225                                            # oxygen (mmol/m^2)
786    coupling:
787      o_pel: O2/o                                          # pelagic oxygen (mmol/m^3)
788  G3:
789    long_name: benthic dissolved inorganic carbon
790    model: ersem/benthic_column_dissolved_matter
791    parameters:
792      composition: c                                       # composition (any combination of c,n,p,s,o,a)
793    initialization:
794      c: 375.0                                             # carbon (mmol/m^2)
795    coupling:
796      c_pel: O3/c                                          # pelagic carbon (mmol/m^3)
797  Q6s_aerobic:
798    long_name: aerobic particulate silicate
799    model: ersem/benthic_column_particulate_matter_layer
800    parameters:
801      composition: s                                       # elemental composition
802      remin: 0.007                                         # remineralization rate (1/d), default = 0.0
803      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
804      variable_maximum_depth: true                         # maximum depth is variable, default = false
805    coupling:
806      Q: Q6
807      maximum_depth: ben_col/D1m                           # maximum depth (m)
808      s_remin_target: K5/per_layer/s1                      # sink for remineralized silicate (mmol Si/m^2)
809  Q6s_anaerobic:
810    long_name: anaerobic particulate silicate
811    model: ersem/benthic_column_particulate_matter_layer
812    parameters:
813      composition: s                                       # elemental composition
814      remin: 0.007                                         # remineralization rate (1/d), default = 0.0
815      variable_minimum_depth: true                         # minimum depth is variable, default = false
816      maximum_depth: 0.3                                   # maximum depth (m), default = 0.0
817      source_depth_distribution: 2                         # vertical distribution of changes (1: constant absolute rate, 2: constant relative rate, 3: constant carbon-based relative rate), default = 1
818    coupling:
819      Q: Q6
820      minimum_depth: ben_col/D1m                           # minimum depth (m)
821      s_remin_target: K5/per_layer/s2                      # sink for remineralized silicate (mmol Si/m^2)
822  ben_nit:
823    long_name: benthic nitrogen cycle
824    model: ersem/benthic_nitrogen_cycle
825    parameters:
826      q10nit: 2.0                                          # Q_10 temperature coefficient (-)
827      hM4M3: 10.0                                          # Michaelis-Menten constant for nitrate limitation (mmol/m^3)
828      ISWph: 1                                             # pH impact on nitrification (0: off, 1: on), default = 0
829      sM4M3: 4.0                                           # maximum nitrification rate at 10 degrees Celsius (1/d)
830      xno3: 2.0                                            # oxygen consumed per nitrate produced (mol O_2/mol N)
831      pammon: 0.5                                          # fraction of oxygen demand fulfilled by denitrification under anaerobic conditions (-)
832      pdenit: 0.05                                         # fraction of denitrification producing dinitrogen gas (remainder produces ammonium) (-)
833      xn2: 1.25                                            # oxygen produced per N2 produced (-)
834      hM3G4: 1.0                                           # Michaelis-Menten constant for nitrate limitation of denitrification (mmol N/m^3)
835    coupling:
836      K3n: K3/per_layer/n1                                 # nitrate (mmol N/m^2)
837      K4n: K4/per_layer/n1                                 # ammonium (mmol N/m^2)
838      G2o: G2/per_layer/o1                                 # oxygen (mmol O_2/m^2)
839      N4n: N4/n                                            # pelagic ammonium (mmol N/m^3)
840      benTA: zero_hz                                       # benthic alkalinity in aerobic layer (mEq/m^2)
841      benTA2: zero_hz                                      # benthic alkalinity in anaerobic layer (mEq/m^2)
842      K3n2: K3/per_layer/n2                                # benthic nitrate in 2nd layer (mmol N/m^2)
843      K4n2: K4/per_layer/n2                                # benthic ammonium in 2nd layer (mmol N/m^2)
844      G2o2: G2/per_layer/o2                                # oxygen in 2nd layer (mmol O_2/m^2)
845      layer2_thickness: ben_col/layer2_thickness           # thickness of 2nd layer (m)
846  H1:
847    long_name: benthic aerobic bacteria
848    model: ersem/benthic_bacteria
849    parameters:
850      qnc: 0.0167                                          # nitrogen to carbon ratio (mmol N/mg C)
851      qpc: 0.00125                                         # phosphorus to carbon ratio (mmol P/mg C)
852      c0: 1.e-06                                           # background carbon concentration (mg C/m^2)
853      q10: 2.0                                             # Q_10 temperature coefficient (-)
854      dd: 0.001                                            # Michaelis-Menten constant for oxygen limitation through layer thickness (m)
855      nfood: 3                                             # number of food sources
856      su1: 0.0005                                          # base affinity for food source 1 (m^2/mg C/d)
857      su2: 2.e-05                                          # base affinity for food source 2 (m^2/mg C/d)
858      suf2: 0.0002                                         # additional affinity for food source 2 subject to nutrient limitation (m^2/mg C/d), default = 0.0
859      puinc2: 2.0                                          # affinity for nutrients relative to carbon in food source 2 (-), default = 1.0
860      pue2: 0.1                                            # fraction of consumed food source 2 that is excreted (-), default = 0.0
861      su3: 2.e-06                                          # base affinity for food source 3 (m^2/mg C/d)
862      pue3: 0.1                                            # fraction of consumed food source 3 that is excreted (-), default = 0.0
863      pur: 0.3                                             # fraction of consumed carbon that is respired (-)
864      sr: 0.02                                             # specific rest respiration (1/d)
865      pdQ1: 0.1                                            # fraction of dying matter that is dissolved (-)
866      sd: 0.05                                             # specific maximum mortality related to oxygen limitation (1/d)
867    initialization:
868      c: 25.0                                              # carbon (mg C/m^2)
869    coupling:
870      Q1: Q1
871      Q6: H1_Q6
872      food1: Q1
873      food2: H1_Q6
874      food3: H1_Q7
875      Dm: ben_col/D1m                                      # depth interval available to bacteria (m)
876      K4n: K4/per_layer/n1_pw                              # ammonium (mmol N/m^2)
877      K1p: K1/per_layer/p1_pw                              # phosphate (mmol N/m^2)
878      G2o: G2/per_layer/o1                                 # oxygen (mmol O_2/m^2)
879      G3c: G3/per_layer/c1                                 # dissolved inorganic carbon (mmol C/m^2)
880      benTA: zero_hz                                       # benthic alkalinity (mEq/m^2)
881  H1_Q6:
882    long_name: benthic aerobic bacteria POM
883    model: ersem/benthic_column_particulate_matter_layer
884    parameters:
885      composition: cnps                                    # elemental composition
886      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
887      variable_maximum_depth: true                         # maximum depth is variable, default = false
888    coupling:
889      Q: Q6
890      maximum_depth: ben_col/D1m                           # maximum depth (m)
891  H1_Q7:
892    long_name: benthic aerobic bacteria refractory POM
893    model: ersem/benthic_column_particulate_matter_layer
894    parameters:
895      composition: cnp                                     # elemental composition
896      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
897      variable_maximum_depth: true                         # maximum depth is variable, default = false
898    coupling:
899      Q: Q7
900      maximum_depth: ben_col/D1m                           # maximum depth (m)
901  H2:
902    long_name: benthic anaerobic bacteria
903    model: ersem/benthic_bacteria
904    parameters:
905      qnc: 0.0167                                          # nitrogen to carbon ratio (mmol N/mg C)
906      qpc: 0.00125                                         # phosphorus to carbon ratio (mmol P/mg C)
907      c0: 1.e-06                                           # background carbon concentration (mg C/m^2)
908      q10: 2.0                                             # Q_10 temperature coefficient (-)
909      dd: 0.01                                             # Michaelis-Menten constant for oxygen limitation through layer thickness (m)
910      nfood: 2                                             # number of food sources
911      su1: 2e-05                                           # base affinity for food source 1 (m^2/mg C/d)
912      suf1: 0.0002                                         # additional affinity for food source 1 subject to nutrient limitation (m^2/mg C/d), default = 0.0
913      puinc1: 2.0                                          # affinity for nutrients relative to carbon in food source 1 (-), default = 1.0
914      pue1: 0.0                                            # fraction of consumed food source 1 that is excreted (-), default = 0.0
915      su2: 2e-06                                           # base affinity for food source 2 (m^2/mg C/d)
916      pue2: 0.0                                            # fraction of consumed food source 2 that is excreted (-), default = 0.0
917      pur: 0.3                                             # fraction of consumed carbon that is respired (-)
918      sr: 0.02                                             # specific rest respiration (1/d)
919      pdQ1: 0.0                                            # fraction of dying matter that is dissolved (-)
920      sd: 0.05                                             # specific maximum mortality related to oxygen limitation (1/d)
921    initialization:
922      c: 400.0                                             # carbon (mg C/m^2)
923    coupling:
924      Q1: Q1
925      Q6: H2_Q6
926      food1: H2_Q6
927      food2: H2_Q7
928      Dm: ben_col/layer2_thickness                         # depth interval available to bacteria (m)
929      K4n: K4/per_layer/n2_pw                              # ammonium (mmol N/m^2)
930      K1p: K1/per_layer/p2_pw                              # phosphate (mmol N/m^2)
931      G2o: ben_nit/K6_calculator/K6                        # oxygen (mmol O_2/m^2)
932      G3c: G3/per_layer/c2                                 # dissolved inorganic carbon (mmol C/m^2)
933      benTA: zero_hz                                       # benthic alkalinity (mEq/m^2)
934  H2_Q6:
935    long_name: benthic anaerobic bacteria POM
936    model: ersem/benthic_column_particulate_matter_layer
937    parameters:
938      composition: cnps                                    # elemental composition
939      variable_minimum_depth: true                         # minimum depth is variable, default = false
940      maximum_depth: 0.3                                   # maximum depth (m), default = 0.0
941      source_depth_distribution: 3                         # vertical distribution of changes (1: constant absolute rate, 2: constant relative rate, 3: constant carbon-based relative rate), default = 1
942    coupling:
943      Q: Q6
944      minimum_depth: ben_col/D1m                           # minimum depth (m)
945  H2_Q7:
946    long_name: benthic anaerobic bacteria refractory POM
947    model: ersem/benthic_column_particulate_matter_layer
948    parameters:
949      composition: cnp                                     # elemental composition
950      variable_minimum_depth: true                         # minimum depth is variable, default = false
951      maximum_depth: 0.3                                   # maximum depth (m), default = 0.0
952      source_depth_distribution: 3                         # vertical distribution of changes (1: constant absolute rate, 2: constant relative rate, 3: constant carbon-based relative rate), default = 1
953    coupling:
954      Q: Q7
955      minimum_depth: ben_col/D1m                           # minimum depth (m)
956  Y2:
957    long_name: benthic deposit feeder
958    model: ersem/benthic_fauna
959    parameters:
960      qnc: 0.0119                                          # nitrogen to carbon ratio (mmol N/mg C)
961      qpc: 0.000792                                        # phosphorus to carbon ratio (mmol P/mg C)
962      c0: 0.001                                            # background carbon concentration (mg C/m^2)
963      q10: 2.0                                             # Q_10 temperature coefficient (-)
964      rlO2: 0.0                                            # minimum pelagic oxygen concentration (mmol O2/m^3)
965      hO2: 0.0                                             # Michaelis-Menten constant for oxygen limitation (mmol O2/m^3)
966      xcl: 2500.0                                          # abundance above which crowding reduces food uptake (mg C/m^2)
967      xcs: 3000.0                                          # Michaelis-Menten constant for the impact of crowding (mg C/m^2)
968      xch: 5000.0                                          # abundance determining asymptotic threshold of crowding limitation (-> xch/(Yc+xch) for Yc-> inf) (mg C/m^2)
969      su: 0.11                                             # specific maximum uptake at reference temperature (1/d)
970      lu: 250.0                                            # Michaelis-Menten constant for food preference as function of food concentration (mg C/m^2)
971      hu: 3000.0                                           # Michaelis-Menten constant for gross carbon uptake (mg C/m^2)
972      pue: 0.35                                            # fraction of carbon in consumed live food that goes to faeces (-)
973      pueQ: 0.8                                            # fraction of carbon in consumed detritus that goes to faeces (-)
974      pudil: 0.8                                           # relative nutrient content of faeces (-)
975      sd: 0.001                                            # specific mortality at reference temperature (1/d)
976      sdmO2: 1.0                                           # specific maximum additional mortality due to oxygen stress (1/d)
977      sdc: 0.022                                           # specific maximum additional mortality induced by cold temperature (1/d)
978      xdc: -0.7                                            # e-folding temperature factor of cold mortality response (1/degree_C)
979      sr: 0.0027                                           # specific rest respiration at reference temperature (1/d)
980      pur: 0.35                                            # fraction of assimilated food that is respired (-)
981      nfood: 4                                             # number of food sources, default = 0
982      pufood1: 1.0                                         # preference for food source 1 (-)
983      pufood2: 1.0                                         # preference for food source 2 (-)
984      pufood3: 0.1                                         # preference for food source 3 (-)
985      pufood4: 1.0                                         # preference for food source 4 (-)
986      food3ispom: true                                     # food source 3 is detritus, default = false
987      ptur: 1.0                                            # relative contribution to bioturbation (-), default = 0.0
988      pirr: 1.0                                            # relative contribution to bioirrigation (-), default = 0.0
989    initialization:
990      c: 3000.0                                            # carbon (mg C/m^2)
991    coupling:
992      food1: H1
993      food2: H2
994      food3: Y2_Q6_source
995      food4: Y4
996      Q: Y2_Q6_sink
997      Dm: ben_col/D1m                                      # depth of limiting layer for uptake (m)
998      O2o: O2/o                                            # pelagic oxygen (mmol O_2/m^3)
999      G3c: G3/per_layer/c1                                 # carbon dioxide (mmol C/m^2)
1000      G2o: G2/per_layer/o1                                 # oxygen (mmol O_2/m^2)
1001      K1p: K1/per_layer/p1                                 # benthic phosphate in aerobic layer (mmol P/m^2)
1002      K4n: K4/per_layer/n1                                 # benthic ammonium in aerobic layer (mmol N/m^2)
1003      K1p2: K1/per_layer/p2                                # benthic phosphate in anaerobic layer (mmol P/m^2)
1004      K4n2: K4/per_layer/n2                                # benthic ammonium in anaerobic layer (mmol N/m^2)
1005      benTA: zero_hz                                       # benthic alkalinity in aerobic layer (mEq/m^2)
1006      benTA2: zero_hz                                      # benthic alkalinity in anaerobic layer (mEq/m^2)
1007      food2c_an: H2/c                                      # food source 2 carbon in anaerobic layer (mg C/m^2)
1008      food3c_an: Y2_Q6_source_anaerobic/c                  # food source 3 carbon in anaerobic layer (mg C/m^2)
1009  Y2_Q6_source:
1010    long_name: benthic deposit feeder total POM
1011    model: ersem/benthic_column_particulate_matter_layer
1012    parameters:
1013      composition: cnps                                    # elemental composition
1014      minimum_depth: 0.0025                                # minimum depth (m), default = 0.0
1015      maximum_depth: 0.3                                   # maximum depth (m), default = 0.0
1016    coupling:
1017      Q: Q6
1018  Y2_Q6_source_anaerobic:
1019    long_name: benthic deposit feeder anaerobic POM
1020    model: ersem/benthic_column_particulate_matter_layer
1021    parameters:
1022      composition: cnps                                    # elemental composition
1023      variable_minimum_depth: true                         # minimum depth is variable, default = false
1024      maximum_depth: 0.3                                   # maximum depth (m), default = 0.0
1025    coupling:
1026      Q: Q6
1027      minimum_depth: ben_col/D1m                           # minimum depth (m)
1028  Y2_Q6_sink:
1029    long_name: benthic deposit feeder POM sink
1030    model: ersem/benthic_column_particulate_matter_layer
1031    parameters:
1032      composition: cnps                                    # elemental composition
1033      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
1034      variable_maximum_depth: true                         # maximum depth is variable, default = false
1035    coupling:
1036      Q: Q6
1037      maximum_depth: Q6/pen_depth_c                        # maximum depth (m)
1038  Y3:
1039    long_name: benthic suspension feeder
1040    model: ersem/benthic_fauna
1041    parameters:
1042      qnc: 0.0119                                          # nitrogen to carbon ratio (mmol N/mg C)
1043      qpc: 0.000792                                        # phosphorus to carbon ratio (mmol P/mg C)
1044      c0: 0.001                                            # background carbon concentration (mg C/m^2)
1045      q10: 2.0                                             # Q_10 temperature coefficient (-)
1046      rlO2: 0.0                                            # minimum pelagic oxygen concentration (mmol O2/m^3)
1047      hO2: 0.0                                             # Michaelis-Menten constant for oxygen limitation (mmol O2/m^3)
1048      xcl: 2500.0                                          # abundance above which crowding reduces food uptake (mg C/m^2)
1049      xcs: 3000.0                                          # Michaelis-Menten constant for the impact of crowding (mg C/m^2)
1050      xch: 5000.0                                          # abundance determining asymptotic threshold of crowding limitation (-> xch/(Yc+xch) for Yc-> inf) (mg C/m^2)
1051      su: 0.09                                             # specific maximum uptake at reference temperature (1/d)
1052      lu: 10.0                                             # Michaelis-Menten constant for food preference as function of food concentration (mg C/m^2)
1053      hu: 300.0                                            # Michaelis-Menten constant for gross carbon uptake (mg C/m^2)
1054      pue: 0.35                                            # fraction of carbon in consumed live food that goes to faeces (-)
1055      pueQ: 0.85                                           # fraction of carbon in consumed detritus that goes to faeces (-)
1056      pudil: 0.8                                           # relative nutrient content of faeces (-)
1057      sd: 0.001                                            # specific mortality at reference temperature (1/d)
1058      sdmO2: 1.0                                           # specific maximum additional mortality due to oxygen stress (1/d)
1059      sdc: 0.022                                           # specific maximum additional mortality induced by cold temperature (1/d)
1060      xdc: -0.7                                            # e-folding temperature factor of cold mortality response (1/degree_C)
1061      sr: 0.0027                                           # specific rest respiration at reference temperature (1/d)
1062      pur: 0.4                                             # fraction of assimilated food that is respired (-)
1063      nfood: 6                                             # number of food sources, default = 0
1064      food1ispel: true                                     # food source 1 is pelagic, default = false
1065      food2ispel: true                                     # food source 2 is pelagic, default = false
1066      food3ispel: true                                     # food source 3 is pelagic, default = false
1067      food4ispel: true                                     # food source 4 is pelagic, default = false
1068      food5_ll: true                                       # availability of food source 5 is limited by aerobic layer height, default = false
1069      pufood1: 1.0                                         # preference for food source 1 (-)
1070      pufood2: 1.0                                         # preference for food source 2 (-)
1071      pufood3: 1.0                                         # preference for food source 3 (-)
1072      pufood4: 1.0                                         # preference for food source 4 (-)
1073      pufood5: 1.0                                         # preference for food source 5 (-)
1074      pufood6: 0.1                                         # preference for food source 6 (-)
1075      food4ispom: true                                     # food source 4 is detritus, default = false
1076      food6ispom: true                                     # food source 6 is detritus, default = false
1077      dwat: 1.0                                            # water layer accessible for pelagic food uptake (m), default = 1.0
1078      dQ6: 0.0025                                          # depth of available sediment layer (m), default = 0.0
1079    initialization:
1080      c: 2000.0                                            # carbon (mg C/m^2)
1081    coupling:
1082      food1: P1
1083      food2: P2
1084      food3: P3
1085      food4: R6
1086      food5: H1
1087      food6: Y3_Q6
1088      Q: Y3_Q6
1089      Dm: ben_col/D1m                                      # depth of limiting layer for uptake (m)
1090      O2o: O2/o                                            # pelagic oxygen (mmol O_2/m^3)
1091      G3c: G3/per_layer/c1                                 # carbon dioxide (mmol C/m^2)
1092      G2o: G2/per_layer/o1                                 # oxygen (mmol O_2/m^2)
1093      K1p: K1/per_layer/p1                                 # benthic phosphate in aerobic layer (mmol P/m^2)
1094      K4n: K4/per_layer/n1                                 # benthic ammonium in aerobic layer (mmol N/m^2)
1095      K1p2: K1/per_layer/p2                                # benthic phosphate in anaerobic layer (mmol P/m^2)
1096      K4n2: K4/per_layer/n2                                # benthic ammonium in anaerobic layer (mmol N/m^2)
1097      benTA: zero_hz                                       # benthic alkalinity in aerobic layer (mEq/m^2)
1098      benTA2: zero_hz                                      # benthic alkalinity in anaerobic layer (mEq/m^2)
1099  Y3_Q6:
1100    long_name: benthic suspension feeder POM
1101    model: ersem/benthic_column_particulate_matter_layer
1102    parameters:
1103      composition: cnps                                    # elemental composition
1104      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
1105      maximum_depth: 0.0025                                # maximum depth (m), default = 0.0
1106    coupling:
1107      Q: Q6
1108  Y4:
1109    long_name: benthic meiofauna
1110    model: ersem/benthic_fauna
1111    parameters:
1112      qnc: 0.0119                                          # nitrogen to carbon ratio (mmol N/mg C)
1113      qpc: 0.000792                                        # phosphorus to carbon ratio (mmol P/mg C)
1114      c0: 0.001                                            # background carbon concentration (mg C/m^2)
1115      q10: 2.0                                             # Q_10 temperature coefficient (-)
1116      rlO2: 0.0                                            # minimum pelagic oxygen concentration (mmol O2/m^3)
1117      hO2: 0.0                                             # Michaelis-Menten constant for oxygen limitation (mmol O2/m^3)
1118      xcl: 25000.0                                         # abundance above which crowding reduces food uptake (mg C/m^2)
1119      xcs: 3000.0                                          # Michaelis-Menten constant for the impact of crowding (mg C/m^2)
1120      xch: 5000.0                                          # abundance determining asymptotic threshold of crowding limitation (-> xch/(Yc+xch) for Yc-> inf) (mg C/m^2)
1121      su: 0.4                                              # specific maximum uptake at reference temperature (1/d)
1122      lu: 50.0                                             # Michaelis-Menten constant for food preference as function of food concentration (mg C/m^2)
1123      hu: 1000.0                                           # Michaelis-Menten constant for gross carbon uptake (mg C/m^2)
1124      pue: 0.25                                            # fraction of carbon in consumed live food that goes to faeces (-)
1125      pueQ: 0.4                                            # fraction of carbon in consumed detritus that goes to faeces (-)
1126      pudil: 0.8                                           # relative nutrient content of faeces (-)
1127      sd: 0.01                                             # specific mortality at reference temperature (1/d)
1128      sdmO2: 1.0                                           # specific maximum additional mortality due to oxygen stress (1/d)
1129      sdc: 0.02                                            # specific maximum additional mortality induced by cold temperature (1/d)
1130      xdc: -0.7                                            # e-folding temperature factor of cold mortality response (1/degree_C)
1131      sr: 0.01                                             # specific rest respiration at reference temperature (1/d)
1132      pur: 0.45                                            # fraction of assimilated food that is respired (-)
1133      nfood: 4                                             # number of food sources, default = 0
1134      pufood1: 1.0                                         # preference for food source 1 (-)
1135      pufood2: 1.0                                         # preference for food source 2 (-)
1136      pufood3: 1.0                                         # preference for food source 3 (-)
1137      pufood4: 0.3                                         # preference for food source 4 (-)
1138      food4ispom: true                                     # food source 4 is detritus, default = false
1139      pirr: 0.2                                            # relative contribution to bioirrigation (-), default = 0.0
1140    initialization:
1141      c: 200.0                                             # carbon (mg C/m^2)
1142    coupling:
1143      food1: H1
1144      food2: H2
1145      food3: Y4
1146      food4: Y4_Q6
1147      Q: Y4_Q6
1148      Dm: ben_col/D1m                                      # depth of limiting layer for uptake (m)
1149      O2o: O2/o                                            # pelagic oxygen (mmol O_2/m^3)
1150      G3c: G3/per_layer/c1                                 # carbon dioxide (mmol C/m^2)
1151      G2o: G2/per_layer/o1                                 # oxygen (mmol O_2/m^2)
1152      K1p: K1/per_layer/p1                                 # benthic phosphate in aerobic layer (mmol P/m^2)
1153      K4n: K4/per_layer/n1                                 # benthic ammonium in aerobic layer (mmol N/m^2)
1154      K1p2: K1/per_layer/p2                                # benthic phosphate in anaerobic layer (mmol P/m^2)
1155      K4n2: K4/per_layer/n2                                # benthic ammonium in anaerobic layer (mmol N/m^2)
1156      benTA: zero_hz                                       # benthic alkalinity in aerobic layer (mEq/m^2)
1157      benTA2: zero_hz                                      # benthic alkalinity in anaerobic layer (mEq/m^2)
1158  Y4_Q6:
1159    long_name: benthic meiofauna POM
1160    model: ersem/benthic_column_particulate_matter_layer
1161    parameters:
1162      composition: cnps                                    # elemental composition
1163      minimum_depth: 0.0                                   # minimum depth (m), default = 0.0
1164      maximum_depth: 0.03                                  # maximum depth (m), default = 0.0
1165    coupling:
1166      Q: Q6
1167  N3_flux:
1168    model: external_surface_flux
1169    coupling:
1170      target: N3/n
1171  N4_flux:
1172    model: external_surface_flux
1173    coupling:
1174      target: N4/n
1175  erosion:
1176    model: ersem/benthic_erosion
1177  #Diagnostics:
1178  PON:
1179   long_name: particulate organic nitrogen
1180   model: weighted_sum
1181   parameters:
1182      n: 11
1183      weight5: 0.0126
1184   coupling:
1185      term1: P1/n
1186      term2: P2/n
1187      term3: P3/n
1188      term4: P4/n
1189      term5: Z4/c
1190      term6: Z5/n
1191      term7: Z6/n
1192      term8: B1/n
1193      term9: R4/n
1194      term10: R6/n
1195      term11: R8/n
1196  POP:
1197   long_name: particulate organic phosphorus
1198   model: weighted_sum
1199   parameters:
1200      n: 11
1201      weight5: 0.000786
1202   coupling:
1203      term1: P1/p
1204      term2: P2/p
1205      term3: P3/p
1206      term4: P4/p
1207      term5: Z4/c
1208      term6: Z5/p
1209      term7: Z6/p
1210      term8: B1/p
1211      term9: R4/p
1212      term10: R6/p
1213      term11: R8/p
1214  POSi:
1215   model: weighted_sum
1216   parameters:
1217      n: 3
1218   coupling:
1219      term1: P1/s
1220      term2: R6/s
1221      term3: R8/s
1222  P1_fPIRIc:
1223    long_name: diatom loss to organic carbon
1224    model: weighted_sum
1225    parameters:
1226      n: 3
1227    coupling:
1228      term1: P1/fPIR1c
1229      term2: P1/fPIR2c
1230      term3: P1/fPIRPc
1231  P2_fPIRIc:
1232    long_name: nanophytoplankton loss to organic carbon
1233    model: weighted_sum
1234    parameters:
1235      n: 3
1236    coupling:
1237      term1: P2/fPIR1c
1238      term2: P2/fPIR2c
1239      term3: P2/fPIRPc
1240  P3_fPIRIc:
1241    long_name: picophytoplankton loss to organic carbon
1242    model: weighted_sum
1243    parameters:
1244      n: 3
1245    coupling:
1246      term1: P3/fPIR1c
1247      term2: P3/fPIR2c
1248      term3: P3/fPIRPc
1249  P4_fPIRIc:
1250    long_name: microphytoplankton loss to organic carbon
1251    model: weighted_sum
1252    parameters:
1253      n: 3
1254    coupling:
1255      term1: P4/fPIR1c
1256      term2: P4/fPIR2c
1257      term3: P4/fPIRPc
1258  phytoplankton_loss_to_labile_DOC:
1259   long_name: phytoplankton loss to labile dissolved organic carbon
1260   model: weighted_sum
1261   parameters:
1262      n: 4
1263   coupling:
1264      term1: P1/fPIR1c
1265      term2: P2/fPIR1c
1266      term3: P3/fPIR1c
1267      term4: P4/fPIR1c
1268  phytoplankton_loss_to_labile_DON:
1269   long_name: phytoplankton loss to labile dissolved organic nitrogen
1270   model: weighted_sum
1271   parameters:
1272      n: 4
1273   coupling:
1274      term1: P1/fPIR1n
1275      term2: P2/fPIR1n
1276      term3: P3/fPIR1n
1277      term4: P4/fPIR1n
1278  phytoplankton_loss_to_labile_DOP:
1279   long_name: phytoplankton loss to labile dissolved organic phosphorus
1280   model: weighted_sum
1281   parameters:
1282      n: 4
1283   coupling:
1284      term1: P1/fPIR1p
1285      term2: P2/fPIR1p
1286      term3: P3/fPIR1p
1287      term4: P4/fPIR1p
1288  phytoplankton_loss_to_non_labile_DOC:
1289   long_name: phytoplankton loss to non-labile dissolved organic carbon
1290   model: weighted_sum
1291   parameters:
1292      n: 4
1293   coupling:
1294      term1: P1/fPIR2c
1295      term2: P2/fPIR2c
1296      term3: P3/fPIR2c
1297      term4: P4/fPIR2c
1298  phytoplankton_loss_to_POC:
1299   long_name: phytoplankton loss to particulate organic carbon
1300   model: weighted_sum
1301   parameters:
1302      n: 4
1303   coupling:
1304      term1: P1/fPIRPc
1305      term2: P2/fPIRPc
1306      term3: P3/fPIRPc
1307      term4: P4/fPIRPc
1308  phytoplankton_loss_to_PON:
1309   long_name: phytoplankton loss to particulate organic nitrogen
1310   model: weighted_sum
1311   parameters:
1312      n: 4
1313   coupling:
1314      term1: P1/fPIRPn
1315      term2: P2/fPIRPn
1316      term3: P3/fPIRPn
1317      term4: P4/fPIRPn
1318  phytoplankton_loss_to_POP:
1319   long_name: phytoplankton loss to particulate organic phosphorus
1320   model: weighted_sum
1321   parameters:
1322      n: 4
1323   coupling:
1324      term1: P1/fPIRPp
1325      term2: P2/fPIRPp
1326      term3: P3/fPIRPp
1327      term4: P4/fPIRPp
1328  phytoplankton_uptake_oxN:
1329   long_name: phytoplankton uptake of oxidized nitrogen
1330   model: weighted_sum
1331   parameters:
1332      n: 4
1333   coupling:
1334      term1: P1/fN3PIn
1335      term2: P2/fN3PIn
1336      term3: P3/fN3PIn
1337      term4: P4/fN3PIn
1338  phytoplankton_uptake_redN:
1339   long_name: phytoplankton uptake of reduced nitrogen
1340   model: weighted_sum
1341   parameters:
1342      n: 4
1343   coupling:
1344      term1: P1/fN4PIn
1345      term2: P2/fN4PIn
1346      term3: P3/fN4PIn
1347      term4: P4/fN4PIn
1348  phytoplankton_uptake_DIP:
1349   long_name: phytoplankton uptake of dissolved inorganic phosphorus
1350   model: weighted_sum
1351   parameters:
1352      n: 4
1353   coupling:
1354      term1: P1/fN1PIp
1355      term2: P2/fN1PIp
1356      term3: P3/fN1PIp
1357      term4: P4/fN1PIp
1358  Z4_fphytoc:
1359   long_name: mesozooplankton grazing of phytoplankton carbon
1360   model: weighted_sum
1361   parameters:
1362      n: 4
1363   coupling:
1364      term1: Z4/fprey1c
1365      term2: Z4/fprey2c
1366      term3: Z4/fprey3c
1367      term4: Z4/fprey4c
1368  Z4_fphyton:
1369   long_name: mesozooplankton grazing of phytoplankton nitrogen
1370   model: weighted_sum
1371   parameters:
1372      n: 4
1373   coupling:
1374      term1: Z4/fprey1n
1375      term2: Z4/fprey2n
1376      term3: Z4/fprey3n
1377      term4: Z4/fprey4n
1378  Z4_fphytop:
1379   long_name: mesozooplankton grazing of phytoplankton phosphorus
1380   model: weighted_sum
1381   parameters:
1382      n: 4
1383   coupling:
1384      term1: Z4/fprey1p
1385      term2: Z4/fprey2p
1386      term3: Z4/fprey3p
1387      term4: Z4/fprey4p
1388  Z4_fzooc:
1389   long_name: mesozooplankton predation of phytoplankton carbon
1390   model: weighted_sum
1391   parameters:
1392      n: 3
1393   coupling:
1394      term1: Z4/fprey5c
1395      term2: Z4/fprey6c
1396      term3: Z4/fprey7c
1397  Z4_fzoon:
1398   long_name: mesozooplankton predation of phytoplankton nitrogen
1399   model: weighted_sum
1400   parameters:
1401      n: 3
1402   coupling:
1403      term1: Z4/fprey5n
1404      term2: Z4/fprey6n
1405      term3: Z4/fprey7n
1406  Z4_fzoop:
1407   long_name: mesozooplankton predation of zooplankton phosphorus
1408   model: weighted_sum
1409   parameters:
1410      n: 3
1411   coupling:
1412      term1: Z4/fprey5p
1413      term2: Z4/fprey6p
1414      term3: Z4/fprey7p
1415  Z5_fphytoc:
1416   long_name: microzooplankton grazing of phytoplankton carbon
1417   model: weighted_sum
1418   parameters:
1419      n: 4
1420   coupling:
1421      term1: Z5/fprey2c
1422      term2: Z5/fprey3c
1423      term3: Z5/fprey4c
1424      term4: Z5/fprey5c
1425  Z5_fphyton:
1426   long_name: microzooplankton grazing of phytoplankton nitrogen
1427   model: weighted_sum
1428   parameters:
1429      n: 4
1430   coupling:
1431      term1: Z5/fprey2n
1432      term2: Z5/fprey3n
1433      term3: Z5/fprey4n
1434      term4: Z5/fprey5n
1435  Z5_fphytop:
1436   long_name: microzooplankton grazing of phytoplankton phosphorus
1437   model: weighted_sum
1438   parameters:
1439      n: 4
1440   coupling:
1441      term1: Z5/fprey2p
1442      term2: Z5/fprey3p
1443      term3: Z5/fprey4p
1444      term4: Z5/fprey5p
1445  Z5_fzooc:
1446   long_name: microzooplankton predation of zooplankton carbon
1447   model: weighted_sum
1448   parameters:
1449      n: 2
1450   coupling:
1451      term1: Z5/fprey6c
1452      term2: Z5/fprey7c
1453  Z5_fzoon:
1454   long_name: microzooplankton predation of zooplankton nitrogen
1455   model: weighted_sum
1456   parameters:
1457      n: 2
1458   coupling:
1459      term1: Z5/fprey6n
1460      term2: Z5/fprey7n
1461  Z5_fzoop:
1462   long_name: microzooplankton predation of zooplankton phosphorus
1463   model: weighted_sum
1464   parameters:
1465      n: 2
1466   coupling:
1467      term1: Z5/fprey6p
1468      term2: Z5/fprey7p
1469  Z6_fphytoc:
1470   long_name: heterotrophic flagellate grazing of phytoplankton carbon
1471   model: weighted_sum
1472   parameters:
1473      n: 2
1474   coupling:
1475      term1: Z6/fprey2c
1476      term2: Z6/fprey3c
1477  Z6_fphyton:
1478   long_name: heterotrophic flagellate grazing of phytoplankton nitrogen
1479   model: weighted_sum
1480   parameters:
1481      n: 2
1482   coupling:
1483      term1: Z6/fprey2n
1484      term2: Z6/fprey3n
1485  Z6_fphytop:
1486   long_name: heterotrophic flagellate grazing of phytoplankton phosphorus
1487   model: weighted_sum
1488   parameters:
1489      n: 2
1490   coupling:
1491      term1: Z6/fprey2p
1492      term2: Z6/fprey3p
1493  total_microzooplankton_loss_to_POC:
1494   long_name: total microzooplankton loss to particulate organic carbon
1495   model: weighted_sum
1496   parameters:
1497      n: 2
1498   coupling:
1499      term1: Z5/fZIRPc
1500      term2: Z6/fZIRPc
1501  total_microzooplankton_loss_to_PON:
1502   long_name: total microzooplankton loss to particulate organic nitrogen
1503   model: weighted_sum
1504   parameters:
1505      n: 2
1506   coupling:
1507      term1: Z5/fZIRPn
1508      term2: Z6/fZIRPn
1509  total_microzooplankton_loss_to_POP:
1510   long_name: total microzooplankton loss to particulate organic phosphorus
1511   model: weighted_sum
1512   parameters:
1513      n: 2
1514   coupling:
1515      term1: Z5/fZIRPp
1516      term2: Z6/fZIRPp
1517  total_microzooplankton_loss_to_DOC:
1518   long_name: total microzooplankton loss to dissolved organic carbon
1519   model: weighted_sum
1520   parameters:
1521      n: 2
1522   coupling:
1523      term1: Z5/fZIRDc
1524      term2: Z6/fZIRDc
1525  total_microzooplankton_loss_to_DON:
1526   long_name: total microzooplankton loss to dissolved organic nitrogen
1527   model: weighted_sum
1528   parameters:
1529      n: 2
1530   coupling:
1531      term1: Z5/fZIRDn
1532      term2: Z6/fZIRDn
1533  total_microzooplankton_loss_to_DOP:
1534   long_name: total microzooplankton loss to dissolved organic phosphorus
1535   model: weighted_sum
1536   parameters:
1537      n: 2
1538   coupling:
1539      term1: Z5/fZIRDp
1540      term2: Z6/fZIRDp
1541  total_microzooplankton_respiration:
1542   long_name: total microzooplankton respiration
1543   model: weighted_sum
1544   parameters:
1545      n: 2
1546   coupling:
1547      term1: Z5/fZIO3c
1548      term2: Z6/fZIO3c
1549  total_microzooplankton_loss_to_DIN:
1550   long_name: total microzooplankton loss to dissolved inorganic nitrogen
1551   model: weighted_sum
1552   parameters:
1553      n: 2
1554   coupling:
1555      term1: Z5/fZINIn
1556      term2: Z6/fZINIn
1557  total_microzooplankton_loss_to_DIP:
1558   long_name: total microzooplankton loss to dissolved inorganic phosphorus
1559   model: weighted_sum
1560   parameters:
1561      n: 2
1562   coupling:
1563      term1: Z5/fZINIp
1564      term2: Z6/fZINIp
1565  gross_calcification_rate:
1566    long_name: gross calcification rate
1567    model: weighted_sum
1568    parameters:
1569      n: 4
1570    coupling:
1571      term1: P2/calcification
1572      term2: Z4/calcification
1573      term3: Z5/calcification
1574      term4: Z6/calcification
1575  ben_PON:
1576   long_name: benthic particulate organic nitrogen
1577   model: horizontal_weighted_sum
1578   parameters:
1579      n: 2
1580   coupling:
1581      term1: Q6/n
1582      term2: Q7/n
1583  ben_POP:
1584   long_name: benthic particulate organic phosphate
1585   model: horizontal_weighted_sum
1586   parameters:
1587      n: 2
1588   coupling:
1589      term1: Q6/p
1590      term2: Q7/p
1591  benthic_denitrification_ammonification:
1592    long_name: benthic denitrification and ammonification
1593    model: horizontal_weighted_sum
1594    parameters:
1595      n: 2
1596    coupling:
1597      term1: ben_nit/jM3G4n
1598      term2: ben_nit/jM3M4n
1599  Y3_fphytoc:
1600   long_name: benthic suspension feeder grazing of phytoplankton carbon
1601   model: horizontal_weighted_sum
1602   parameters:
1603      n: 3
1604   coupling:
1605      term1: Y3/fprey1c
1606      term2: Y3/fprey2c
1607      term3: Y3/fprey3c
1608  Y3_fphyton:
1609   long_name: benthic suspension feeder grazing of phytoplankton nitrogen
1610   model: horizontal_weighted_sum
1611   parameters:
1612      n: 3
1613   coupling:
1614      term1: Y3/fprey1n
1615      term2: Y3/fprey2n
1616      term3: Y3/fprey3n
1617  Y3_fphytop:
1618   long_name: benthic suspension feeder grazing of phytoplankton phosphorus
1619   model: horizontal_weighted_sum
1620   parameters:
1621      n: 3
1622   coupling:
1623      term1: Y3/fprey1p
1624      term2: Y3/fprey2p
1625      term3: Y3/fprey3p
1626  Ymacro_fYG3c:
1627   long_name: benthic macrofauna respiration
1628   model: horizontal_weighted_sum
1629   parameters:
1630      n: 2
1631   coupling:
1632      term1: Y2/fYG3c
1633      term2: Y3/fYG3c
1634  Ymacro_fYQPc:
1635   long_name: benthic macrofauna loass to particulate organic carbon
1636   model: horizontal_weighted_sum
1637   parameters:
1638      n: 2
1639   coupling:
1640      term1: Y2/fYQPc
1641      term2: Y3/fYQPc
1642  Ymacro_fYQPn:
1643   long_name: benthic macrofauna loass to particulate organic nitrogen
1644   model: horizontal_weighted_sum
1645   parameters:
1646      n: 2
1647   coupling:
1648      term1: Y2/fYQPn
1649      term2: Y3/fYQPn
1650  Ymacro_fYQPp:
1651   long_name: benthic macrofauna loass to particulate organic phosphorus
1652   model: horizontal_weighted_sum
1653   parameters:
1654      n: 2
1655   coupling:
1656      term1: Y2/fYQPp
1657      term2: Y3/fYQPp
1658  Ymacro_fYKIn:
1659   long_name: benthic macrofauna loass to dissolved inorganic nitrogen
1660   model: horizontal_weighted_sum
1661   parameters:
1662      n: 2
1663   coupling:
1664      term1: Y2/fYKIn
1665      term2: Y3/fYKIn
1666  Ymacro_fYK1p:
1667   long_name: benthic macrofauna loass to dissolved inorganic phosphorus
1668   model: horizontal_weighted_sum
1669   parameters:
1670      n: 2
1671   coupling:
1672      term1: Y2/fYK1p
1673      term2: Y3/fYK1p
1674  Q6_remins:
1675    long_name: benthic remineralisation of silicate
1676    model: horizontal_weighted_sum
1677    parameters:
1678      n: 1
1679      weight1: 0.007
1680    coupling:
1681      term1: Q6/s
1682  net_PelBen_POC:
1683   long_name: flux of particulate organic carbon from pelagic to benthos
1684   model: horizontal_weighted_sum
1685   parameters:
1686      n: 17
1687      weight17: -1.
1688   coupling:
1689      term1: R4/dep1c
1690      term2: R4/dep2c
1691      term3: R6/dep1c
1692      term4: R6/dep2c
1693      term5: R8/dep1c
1694      term6: R8/dep2c
1695      term7: P1/dep1c
1696      term8: P1/dep2c
1697      term9: P1/dep3c
1698      term10: P4/dep1c
1699      term11: P4/dep2c
1700      term12: P4/dep3c
1701      term13: Y3/fprey1c
1702      term14: Y3/fprey2c
1703      term15: Y3/fprey3c
1704      term16: Y3/fprey4c
1705      term17: Q6/resuspension_flux_c
1706  net_PelBen_PON:
1707   long_name: flux of particulate organic nitrogen from pelagic to benthos
1708   model: horizontal_weighted_sum
1709   parameters:
1710      n: 17
1711      weight17: -1.
1712   coupling:
1713      term1: R4/dep1n
1714      term2: R4/dep2n
1715      term3: R6/dep1n
1716      term4: R6/dep2n
1717      term5: R8/dep1n
1718      term6: R8/dep2n
1719      term7: P1/dep1n
1720      term8: P1/dep2n
1721      term9: P1/dep3n
1722      term10: P4/dep1n
1723      term11: P4/dep2n
1724      term12: P4/dep3n
1725      term13: Y3/fprey1n
1726      term14: Y3/fprey2n
1727      term15: Y3/fprey3n
1728      term16: Y3/fprey4n
1729      term17: Q6/resuspension_flux_n
1730  net_PelBen_POP:
1731   long_name: flux of particulate organic phosphorus from pelagic to benthos
1732   model: horizontal_weighted_sum
1733   parameters:
1734      n: 17
1735      weight17: -1.
1736   coupling:
1737      term1: R4/dep1p
1738      term2: R4/dep2p
1739      term3: R6/dep1p
1740      term4: R6/dep2p
1741      term5: R8/dep1p
1742      term6: R8/dep2p
1743      term7: P1/dep1p
1744      term8: P1/dep2p
1745      term9: P1/dep3p
1746      term10: P4/dep1p
1747      term11: P4/dep2p
1748      term12: P4/dep3p
1749      term13: Y3/fprey1p
1750      term14: Y3/fprey2p
1751      term15: Y3/fprey3p
1752      term16: Y3/fprey4p
1753      term17: Q6/resuspension_flux_p
1754  net_PelBen_POSi:
1755   long_name: flux of particulate organic silicate from pelagic to benthos
1756   model: horizontal_weighted_sum
1757   parameters:
1758      n: 6
1759      weight6: -1.
1760   coupling:
1761      term1: R6/dep1s
1762      term2: R8/dep1s
1763      term3: P1/dep1s
1764      term4: Y3/fprey1s
1765      term5: Y3/fprey4s
1766      term6: Q6/resuspension_flux_s
Note: See TracBrowser for help on using the repository browser.