 |
OGS
|
|
OGS
|
<OpenGeoSysProject>
<meshes>
<mesh>Wetland_1d.vtu</mesh>
<mesh>Wetland_1d_inlet.vtu</mesh>
<mesh>Wetland_1d_outlet.vtu</mesh>
<mesh>Wetland_1d_ReactiveDomain.vtu</mesh>
</meshes>
<processes>
<process>
<name>hc</name>
<type>ComponentTransport</type>
<integration_order>2</integration_order>
<coupling_scheme>staggered</coupling_scheme>
<process_variables>
<pressure>pressure</pressure>
<concentration>H</concentration>
<concentration>Do</concentration>
<concentration>Sf</concentration>
<concentration>Sa</concentration>
<concentration>Sin</concentration>
<concentration>Snh</concentration>
<concentration>Sno</concentration>
<concentration>Sulphide</concentration>
<concentration>Sso</concentration>
<concentration>Xs_m</concentration>
<concentration>Xi_m</concentration>
</process_variables>
<specific_body_force>0</specific_body_force>
<secondary_variables>
<secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity" />
</secondary_variables>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>H</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Do</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sf</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sa</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sin</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Snh</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sno</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sulphide</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sso</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xs_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xi_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
</components>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<components>
<component>
<name>Aeration</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xi</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xs</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_m</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xa_on_Snh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
<component>
<name>Xi_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
</components>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Parameter</type>
<parameter_name>kappa0</parameter_name>
</property>
<property>
<name>porosity</name>
<type>Parameter</type>
<parameter_name>porosity</parameter_name>
</property>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.45</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.00</value>
</property>
</properties>
</medium>
<medium id="1">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>H</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Do</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sf</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sa</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sin</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Snh</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sno</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sulphide</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sso</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xs_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xi_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
</components>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<components>
<component>
<name>Aeration</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xi</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xs</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_m</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xa_on_Snh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
<component>
<name>Xi_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
</components>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Parameter</type>
<parameter_name>kappa1</parameter_name>
</property>
<property>
<name>porosity</name>
<type>Parameter</type>
<parameter_name>porosity</parameter_name>
</property>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.4</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.00</value>
</property>
</properties>
</medium>
<medium id="2">
<phases>
<phase>
<type>AqueousLiquid</type>
<components>
<component>
<name>H</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Do</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sf</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sa</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sin</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Snh</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sno</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sulphide</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Sso</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xs_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
<component>
<name>Xi_m</name>
<properties>
<property>
<name>pore_diffusion</name>
<type>Constant</type>
<value>0.0e-9</value>
</property>
<property>
<name>retardation_factor</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>decay_rate</name>
<type>Parameter</type>
<parameter_name>decay</parameter_name>
</property>
</properties>
</component>
</components>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>1e3</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<components>
<component>
<name>Aeration</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xi</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>AttDetachment_Xs</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_m</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Hydrolysis_Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sf</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xh_on_Sa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xa_on_Snh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Growth_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Aerobic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Anoxic_Growth_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Lysis_Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0</value>
</property>
</properties>
</component>
<component>
<name>Xh</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xa</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xfb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xamb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xasrb</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xsob</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-3</value>
</property>
</properties>
</component>
<component>
<name>Xs_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
<component>
<name>Xi_im</name>
<properties>
<property>
<name>molality</name>
<type>Constant</type>
<value>1.0e-4</value>
</property>
</properties>
</component>
</components>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Parameter</type>
<parameter_name>kappa0</parameter_name>
</property>
<property>
<name>porosity</name>
<type>Parameter</type>
<parameter_name>porosity</parameter_name>
</property>
<property>
<name>longitudinal_dispersivity</name>
<type>Constant</type>
<value>0.45</value>
</property>
<property>
<name>transversal_dispersivity</name>
<type>Constant</type>
<value>0.00</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<global_process_coupling>
<max_iter>6</max_iter>
<convergence_criteria>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
</convergence_criteria>
</global_process_coupling>
<processes>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
<process ref="hc">
<nonlinear_solver>basic_picard</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<reltol>1e-14</reltol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0.0</t_initial>
<t_end>28800</t_end>
<timesteps>
<pair>
<repeat>4</repeat>
<delta_t>7200</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>Wetland_1d</prefix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>4</each_steps>
</pair>
</timesteps>
<variables>
<variable>pressure</variable>
<variable>darcy_velocity</variable>
<variable>H</variable>
<variable>Do</variable>
<variable>Sf</variable>
<variable>Sa</variable>
<variable>Sin</variable>
<variable>Snh</variable>
<variable>Sno</variable>
<variable>Sulphide</variable>
<variable>Sso</variable>
<variable>Xs_m</variable>
<variable>Xi_m</variable>
</variables>
<suffix>ts{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<chemical_system chemical_solver="Phreeqc">
<mesh>Wetland_1d_ReactiveDomain</mesh>
<linear_solver>general_linear_solver</linear_solver>
<database>cwm1.dat</database>
<solution>
<temperature>10</temperature>
<pressure>1</pressure>
<pe>4</pe>
<charge_balance>pH</charge_balance>
<components>
<component>Do</component>
<component>Sf</component>
<component>Sa</component>
<component>Sin</component>
<component>Snh</component>
<component>Sno</component>
<component>Sulphide</component>
<component>Sso</component>
<component>Xs_m</component>
<component>Xi_m</component>
</components>
</solution>
<kinetic_reactants>
<kinetic_reactant>
<name>Aeration</name>
<chemical_formula>Do 1</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>AttDetachment_Xi</name>
<chemical_formula>Xi_m -1</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>AttDetachment_Xs</name>
<chemical_formula>Xs_m -1</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Hydrolysis_Xs_m</name>
<chemical_formula>Sf 1 Sin 0 Snh 0.01 Xs_m -1</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Hydrolysis_Xs_im</name>
<chemical_formula>Sf 1 Sin 0 Snh 0.01</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Aerobic_Growth_Xh_on_Sf</name>
<chemical_formula>Do -0.587301587301587 Sf -1.58730158730159 Snh -0.022380952</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Anoxic_Growth_Xh_on_Sf</name>
<chemical_formula>Sf -1.58730158730159 Sno -0.205350205 Snh -0.02238095</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Aerobic_Growth_Xh_on_Sa</name>
<chemical_formula>Do -0.587301587301587 Sa -1.58730158730159 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Anoxic_Growth_Xh_on_Sa</name>
<chemical_formula>Sa -1.58730158730159 Sno -0.205350205 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xh</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Aerobic_Growth_Xa_on_Snh</name>
<chemical_formula>Do -18.04166667 Sno 4.166666667 Snh -4.236666667</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xa</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Growth_Xfb</name>
<chemical_formula>Sf -18.86792453 Sa 17.86792453 Snh 0.496037736</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xfb</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Growth_Xamb</name>
<chemical_formula>Sa -31.25 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xamb</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Growth_Xasrb</name>
<chemical_formula>Sa -20 Sso -9.5 Sulphide 9.5 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xasrb</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Aerobic_Growth_Xsob</name>
<chemical_formula>Do -15.66666667 Sso 8.333333333 Sulphide -8.333333333 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Anoxic_Growth_Xsob</name>
<chemical_formula>Sno -8.380952381 Sso 8.333333333 Sulphide -8.333333333 Snh -0.07</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Lysis_Xsob</name>
<chemical_formula>Sf 0.05 Snh 0.0315</chemical_formula>
<fix_amount>true</fix_amount>
</kinetic_reactant>
<kinetic_reactant>
<name>Xh</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xa</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xfb</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xamb</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xasrb</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xsob</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xs_im</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
<kinetic_reactant>
<name>Xi_im</name>
<chemical_formula>Do 0</chemical_formula>
</kinetic_reactant>
</kinetic_reactants>
<rates>
<rate>
<kinetic_reactant>Aeration</kinetic_reactant>
<expression>
<statement>So2 = MOL("Do")</statement>
<statement>So2_sat = CALC_VALUE("So_sat")</statement>
<statement>if So2 > So2_sat then goto 8</statement>
<statement>kLa = 2 / 3600</statement>
<statement>r_So2 = So2_sat - So2</statement>
<statement>rates = CALC_VALUE("kLa_theta") * kLa * r_So2</statement>
<statement>dGrowth= rates*time</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>AttDetachment_Xi</kinetic_reactant>
<expression>
<statement>Xi_mob = MOL("Xi_m")</statement>
<statement>Xi_immob = KIN("Xi_im")</statement>
<statement>IF Xi_immob >= 30.0 THEN GOTO 10</statement>
<statement>k_att = 0.1 /3600</statement>
<statement>k_det = 0</statement>
<statement>r_Att = k_att * Xi_mob</statement>
<statement>r_Det = k_det * Xi_immob</statement>
<statement>rates = r_Att - r_Det</statement>
<statement>dGrowth= rates*time</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>AttDetachment_Xs</kinetic_reactant>
<expression>
<statement>Xs_mob = MOL("Xs_m")</statement>
<statement>Xs_immob = KIN("Xs_im")</statement>
<statement>k_att = 0.1 /3600</statement>
<statement>k_det = 0</statement>
<statement>r_Att = k_att * Xs_mob</statement>
<statement>r_Det = k_det * Xs_immob</statement>
<statement>rates = r_Att - r_Det</statement>
<statement>dGrowth= rates*time</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Hydrolysis_Xs_m</kinetic_reactant>
<expression>
<statement>Xs_m = MOL("Xs_m")</statement>
<statement>Xfb = KIN("Xfb")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Xs = (Xs_m/(Xh+Xfb))</statement>
<statement>r_Xs = r_Xs / (CALC_VALUE("Kx") + r_Xs)</statement>
<statement>rates = CALC_VALUE("Kh") * r_Xs * (Xh + CALC_VALUE("eta_h") * Xfb)</statement>
<statement>dHydrolysis = rates * time</statement>
<statement>dHydrolysis = dHydrolysis</statement>
<statement>SAVE dHydrolysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Hydrolysis_Xs_im</kinetic_reactant>
<expression>
<statement>Xs_im = KIN("Xs_im")</statement>
<statement>Xfb = KIN("Xfb")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Xs = (Xs_im/(Xh+Xfb))</statement>
<statement>r_Xs = r_Xs / (CALC_VALUE("Kx") + r_Xs)</statement>
<statement>rates = CALC_VALUE("Kh") * r_Xs * (Xh + CALC_VALUE("eta_h") * Xfb)</statement>
<statement>dHydrolysis = rates * time</statement>
<statement>dHydrolysis = dHydrolysis</statement>
<statement>SAVE dHydrolysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Aerobic_Growth_Xh_on_Sf</kinetic_reactant>
<expression>
<statement>Sf = MOL("Sf")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Sa = MOL("Sa")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Sf = Sf/(CALC_VALUE("Ksf")+ Sf)</statement>
<statement>r_So2 = So2/(CALC_VALUE("Koh")+ So2)</statement>
<statement>r_Sa = Sf/(Sf + Sa)</statement>
<statement>r_Snh = Snh /(CALC_VALUE("Knhh")+ Snh)</statement>
<statement>r_Sh2s = CALC_VALUE("Kh2s_h") / (CALC_VALUE("Kh2s_h") + Sh2S )</statement>
<statement>rates = CALC_VALUE("muh") * r_Sf * r_So2 * r_Sa * r_Snh * r_Sh2s * Xh * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Anoxic_Growth_Xh_on_Sf</kinetic_reactant>
<expression>
<statement>Sf = MOL("Sf")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Sa = MOL("Sa")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Sf = Sf/(CALC_VALUE("Ksf")+ Sf)</statement>
<statement>r_Koh = CALC_VALUE("Koh")/(CALC_VALUE("Koh")+ So2)</statement>
<statement>r_Sa = Sf/(Sf + Sa)</statement>
<statement>r_Snh = Snh /(CALC_VALUE("Knhh")+ Snh)</statement>
<statement>r_Sno = Sno /(CALC_VALUE("Knoh")+ Sno)</statement>
<statement>r_Sh2s = CALC_VALUE("Kh2s_h") / (CALC_VALUE("Kh2s_h") + Sh2S)</statement>
<statement>rates = CALC_VALUE("eta_g") * CALC_VALUE("muh") * r_Sf * r_Sa * r_Koh * r_Snh * r_Sno * r_Sh2s * Xh * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Aerobic_Growth_Xh_on_Sa</kinetic_reactant>
<expression>
<statement>Sf = MOL("Sf")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Sa = MOL("Sa")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Sa1 = Sa /(CALC_VALUE("Ksa")+ Sa)</statement>
<statement>r_Sa2 = Sa /(Sf + Sa)</statement>
<statement>r_So2 = So2/(CALC_VALUE("Koh")+ So2)</statement>
<statement>r_Snh = Snh /(CALC_VALUE("Knhh")+ Snh)</statement>
<statement>r_Sh2s = CALC_VALUE("Kh2s_h") / (CALC_VALUE("Kh2s_h") + Sh2S)</statement>
<statement>rates = CALC_VALUE("muh") * r_Sa1 * r_Sa2 * r_So2 * r_Snh * r_Sh2s * Xh * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates * time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Anoxic_Growth_Xh_on_Sa</kinetic_reactant>
<expression>
<statement>Sf = MOL("Sf")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Sa = MOL("Sa")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xh = KIN("Xh")</statement>
<statement>r_Sa1 = Sa/(CALC_VALUE("Ksa")+ Sa)</statement>
<statement>r_Sa2 = Sa/(Sf + Sa)</statement>
<statement>r_Koh = CALC_VALUE("Koh")/(CALC_VALUE("Koh")+ So2)</statement>
<statement>r_Snh = Snh /(CALC_VALUE("Knhh")+ Snh)</statement>
<statement>r_Sno = Sno /(CALC_VALUE("Knoh")+ Sno)</statement>
<statement>r_Sh2s = CALC_VALUE("Kh2s_h") / (CALC_VALUE("Kh2s_h") + Sh2S)</statement>
<statement>rates = CALC_VALUE("eta_g") * CALC_VALUE("muh") * r_Sa1 * r_Sa2 * r_Koh * r_Sno * r_Snh * r_Sh2s * Xh * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xh</kinetic_reactant>
<expression>
<statement>Xh = KIN("Xh")</statement>
<statement>rates = CALC_VALUE("bh") * Xh</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Aerobic_Growth_Xa_on_Snh</kinetic_reactant>
<expression>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xa = KIN("Xa")</statement>
<statement>r_So2_a = So2 /(CALC_VALUE("Koa")+ So2)</statement>
<statement>r_Snh_a = Snh /(CALC_VALUE("Knha")+ Snh)</statement>
<statement>r_Sh2s_a = CALC_VALUE("Kh2s_a") / (CALC_VALUE("Kh2s_a") + Sh2S)</statement>
<statement>rates = CALC_VALUE("mu_a") * r_So2_a * r_Snh_a * r_Sh2s_a * Xa * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xa</kinetic_reactant>
<expression>
<statement>Xa = KIN("Xa")</statement>
<statement>rates = CALC_VALUE("ba") * Xa</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Growth_Xfb</kinetic_reactant>
<expression>
<statement>Sf = MOL("Sf")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xfb = KIN("Xfb")</statement>
<statement>r_Sf_fb = Sf/(CALC_VALUE("Ksfb")+ Sf)</statement>
<statement>r_Ko_fb = CALC_VALUE("Kofb")/(CALC_VALUE("Kofb")+ So2)</statement>
<statement>r_Snh_fb = Snh /(CALC_VALUE("Knhfb")+ Snh)</statement>
<statement>r_Kno_fb = CALC_VALUE("Knofb")/(CALC_VALUE("Knofb")+ Sno)</statement>
<statement>r_Sh2s_fb = CALC_VALUE("Kh2s_fb") / (CALC_VALUE("Kh2s_fb") + Sh2S)</statement>
<statement>rates = CALC_VALUE("mu_fb") * r_Sf_fb * r_Ko_fb * r_Snh_fb * r_Kno_fb * r_Sh2s_fb * Xfb * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xfb</kinetic_reactant>
<expression>
<statement>Xfb = KIN("Xfb")</statement>
<statement>rates = CALC_VALUE("b_fb") * Xfb</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Growth_Xamb</kinetic_reactant>
<expression>
<statement>Sa = MOL("Sa")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xamb = KIN("Xamb")</statement>
<statement>r_Sa_amb = Sa/(CALC_VALUE("Ksamb")+ Sa)</statement>
<statement>r_Ko_amb = CALC_VALUE("Koamb")/(CALC_VALUE("Koamb")+ So2)</statement>
<statement>r_Snh_amb = Snh /(CALC_VALUE("Knhamb")+ Snh)</statement>
<statement>r_Kno_amb = CALC_VALUE("Knoamb") /(CALC_VALUE("Knoamb")+ Sno)</statement>
<statement>r_Sh2s_amb = CALC_VALUE("Kh2s_amb") / (CALC_VALUE("Kh2s_amb") + Sh2S)</statement>
<statement>rates = CALC_VALUE("mu_amb") * r_Sa_amb * r_Ko_amb * r_Snh_amb * r_Kno_amb * r_Sh2s_amb * Xamb * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xamb</kinetic_reactant>
<expression>
<statement>Xamb = KIN("Xamb")</statement>
<statement>rates = CALC_VALUE("b_amb") * Xamb</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Growth_Xasrb</kinetic_reactant>
<expression>
<statement>Sa = MOL("Sa")</statement>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Sso4 = MOL("Sso")</statement>
<statement>Xasrb = KIN("Xasrb")</statement>
<statement>r_Sa_asrb = Sa/(CALC_VALUE("Ks_asrb")+ Sa)</statement>
<statement>r_Ko_asrb = CALC_VALUE("Ko_asrb")/(CALC_VALUE("Ko_asrb")+ So2)</statement>
<statement>r_Snh_asrb = Snh /(CALC_VALUE("Knh_asrb")+ Snh)</statement>
<statement>r_Kno_asrb = CALC_VALUE("Kno_asrb") /(CALC_VALUE("Kno_asrb")+ Sno)</statement>
<statement>r_Sh2s_asrb = CALC_VALUE("Kh2s_asrb") / (CALC_VALUE("Kh2s_asrb") + Sh2S)</statement>
<statement>r_Sso4_asrb = Sso4 / (CALC_VALUE("Kso_asrb") + Sso4)</statement>
<statement>rates = CALC_VALUE("mu_asrb") * r_Sa_asrb * r_Ko_asrb * r_Snh_asrb * r_Kno_asrb * r_Sh2s_asrb * r_Sso4_asrb * Xasrb * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xasrb</kinetic_reactant>
<expression>
<statement>Xasrb = KIN("Xasrb")</statement>
<statement>rates = CALC_VALUE("b_asrb") * Xasrb</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Aerobic_Growth_Xsob</kinetic_reactant>
<expression>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xsob = KIN("Xsob")</statement>
<statement>r_So2_sob = So2 /(CALC_VALUE("Ko_sob")+ So2)</statement>
<statement>r_Snh_sob = Snh /(CALC_VALUE("Knh_sob")+ Snh)</statement>
<statement>r_Sh2s_sob = Sh2s / (CALC_VALUE("Kh2s_sob") + Sh2S)</statement>
<statement>rates = CALC_VALUE("mu_sob") * r_So2_sob * r_Sh2s_sob * r_Snh_sob * Xsob * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Anoxic_Growth_Xsob</kinetic_reactant>
<expression>
<statement>So2 = MOL("Do")</statement>
<statement>Snh = MOL("Snh")</statement>
<statement>Sno = MOL("Sno")</statement>
<statement>Sh2S = MOL("Sulphide")</statement>
<statement>Xsob = KIN("Xsob")</statement>
<statement>r_Ko2_sob = CALC_VALUE("Ko_sob") /(CALC_VALUE("Ko_sob")+ So2)</statement>
<statement>r_Snh_sob = Snh /(CALC_VALUE("Knh_sob")+ Snh)</statement>
<statement>r_Sh2s_sob = Sh2s / (CALC_VALUE("Kh2s_sob") + Sh2S)</statement>
<statement>r_Sno_sob = Sno / (CALC_VALUE("Kno_sob") + Sno)</statement>
<statement>rates = CALC_VALUE("mu_sob") * CALC_VALUE("eta_sob") * r_Ko2_sob * r_Sh2s_sob * r_Snh_sob * r_Sno_sob * Xsob * CALC_VALUE("bac_growth_limit")</statement>
<statement>dGrowth= rates*time</statement>
<statement>dGrowth = dGrowth</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Lysis_Xsob</kinetic_reactant>
<expression>
<statement>Xsob = KIN("Xsob")</statement>
<statement>rates = CALC_VALUE("b_sob") * Xsob</statement>
<statement>dLysis = rates * time</statement>
<statement>dLysis = dLysis</statement>
<statement>SAVE dLysis</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xs_im</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("AttDetachment_Xs")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Hydrolysis_Xs_im")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xh")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xa")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xfb")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xamb")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xasrb")</statement>
<statement>dGrowth = dGrowth + 0.85*KIN_DELTA("Lysis_Xsob")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xi_im</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("AttDetachment_Xi")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xh")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xa")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xfb")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xamb")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xasrb")</statement>
<statement>dGrowth = dGrowth + 0.1*KIN_DELTA("Lysis_Xsob")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xh</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Aerobic_Growth_Xh_on_Sf")</statement>
<statement>dGrowth = dGrowth + KIN_DELTA("Anoxic_Growth_Xh_on_Sf")</statement>
<statement>dGrowth = dGrowth + KIN_DELTA("Aerobic_Growth_Xh_on_Sa")</statement>
<statement>dGrowth = dGrowth + KIN_DELTA("Anoxic_Growth_Xh_on_Sa")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xh")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xa</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Aerobic_Growth_Xa_on_Snh")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xa")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xfb</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Growth_Xfb")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xfb")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xamb</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Growth_Xamb")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xamb")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xasrb</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Growth_Xasrb")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xasrb")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
<rate>
<kinetic_reactant>Xsob</kinetic_reactant>
<expression>
<statement>dGrowth = KIN_DELTA("Aerobic_Growth_Xsob")</statement>
<statement>dGrowth = dGrowth + KIN_DELTA("Anoxic_Growth_Xsob")</statement>
<statement>dGrowth = dGrowth - KIN_DELTA("Lysis_Xsob")</statement>
<statement>SAVE dGrowth</statement>
</expression>
</rate>
</rates>
<knobs>
<max_iter>100</max_iter>
<relative_convergence_tolerance>1e-12</relative_convergence_tolerance>
<tolerance>1e-15</tolerance>
<step_size>100</step_size>
<scaling>0</scaling>
</knobs>
</chemical_system>
<parameters>
<parameter>
<name>kappa0</name>
<type>Constant</type>
<values>1.0e-7</values>
</parameter>
<parameter>
<name>kappa1</name>
<type>Constant</type>
<values>1.0e-8</values>
</parameter>
<parameter>
<name>porosity</name>
<type>Constant</type>
<value>0.38</value>
</parameter>
<parameter>
<name>decay</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>p0</name>
<type>Constant</type>
<value>8829</value>
</parameter>
<parameter>
<name>p_Neumann_inlet</name>
<type>Constant</type>
<value>5.555e-3</value>
</parameter>
<parameter>
<name>p_Dirichlet_right</name>
<type>Constant</type>
<value>8829</value>
</parameter>
<parameter>
<name>c_CurveScaled</name>
<type>Constant</type>
<value>1.0</value>
</parameter>
<parameter>
<name>c_H</name>
<type>Constant</type>
<value>1e-7</value>
</parameter>
<parameter>
<name>c0_CWM1</name>
<type>Constant</type>
<value>1e-4</value>
</parameter>
<parameter>
<name>c_Do_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Do</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Sf_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Sf</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Sa_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Sa</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Sin_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Sin</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Xsm_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Xsm</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Xim_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Xim</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Snh_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Snh</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Sno_inlet</name>
<type>CurveScaled</type>
<curve>c_t_Sno</curve>
<parameter>c_CurveScaled</parameter>
</parameter>
<parameter>
<name>c_Sso_inlet</name>
<type>Constant</type>
<value>5.66e-2</value>
</parameter>
<parameter>
<name>c_Sulphide_inlet</name>
<type>Constant</type>
<value>8.6e-3</value>
</parameter>
</parameters>
<curves>
<curve>
<name>c_t_Do</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>9.00E-04 2.00E-04 0.00018 0.00021 0.00055 0.00074 0.00054 0.00093 0.00085 0.00055 0.00038 0.00034 0.00037 0.00049 0.00062 4.00E-04 0.00072 0.00031 0.00078 0.00061 0.00072 0.00089 0.00054 0.00125 0.00083 0.00052 0.00064 0.00071 0.00087 0.00059 0.00101 0.00097 0.00063 0.00043 0.00068</values>
</curve>
<curve>
<name>c_t_Sf</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>2.29E-01 1.92E-01 0.187574 0.2529343 0.2244159 0.2772126 0.2811686 0.2861437 0.209862 0.1226828 0.3143407 0.3353002 0.3534289 0.103283 0.0971518 5.82E-02 0.1027765 0.1537516 0.1811263 0.0806371 0.1088622 0.1340445 0.1565852 0.1631798 0.1713685 0.1830446 0.1945754 0.1854146 0.1791338 0.1705848 0.1620868 0.154141 0.1461953 0.1365357 0.0911085</values>
</curve>
<curve>
<name>c_t_Sa</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>1.53E-01 1.28E-01 0.1250493 0.1686229 0.1496106 0.1848084 0.1874457 0.1907625 0.139908 0.0817885 0.2095605 0.2235335 0.2356192 0.0688553 0.0647679 3.88E-02 0.0685177 0.1025011 0.1207509 0.053758 0.0725748 0.089363 0.1043901 0.1087865 0.1142457 0.1220298 0.129717 0.1236097 0.1194225 0.1137232 0.1080579 0.1027607 0.0974635 0.0910238 0.060739</values>
</curve>
<curve>
<name>c_t_Sin</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>1.48E-02 1.64E-02 0.0179034 0.0192179 0.0208604 0.0176998 0.0205697 0.0217411 0.0229125 0.0242513 0.0198148 0.0220722 0.023953 0.0227709 0.022143 2.19E-02 0.0230327 0.0214195 0.022098 0.0226606 0.0236999 0.0245914 0.0245276 0.0241703 0.0229284 0.024961 0.0272358 0.0285153 0.0249947 0.0252544 0.0255125 0.0247668 0.0240212 0.0261319 0.0236982</values>
</curve>
<curve>
<name>c_t_Xsm</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>8.98E-02 8.83E-02 0.1103771 0.0021094 7.69E-05 0.0305699 0.0307991 0.0920976 0.1937629 0.3099518 0.0143083 0.001983 0.0161316 0.0884845 0.0836935 9.82E-02 0.1125255 0.0328301 0.0345647 0.1685278 0.2101795 0.1536543 0.1033509 0.0664341 0.078712 0.0885889 0.0987081 0.105789 0.0698635 0.0227084 0.0682557 0.0432921 0.1179384 0.1872539 0.3011724</values>
</curve>
<curve>
<name>c_t_Xim</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>6.93E-04 2.89E-02 0 0.0021094 2.07E-02 0 0.0117404 0 0.0220868 0.0576956 0.0143083 0.0415899 0.0161316 0.0783689 0.0690938 3.78E-02 0.063681 0.1475792 0.0345647 0.1729577 0.0701732 0.0845583 0.0979891 0.1074293 0.1642352 0.1330902 0.1520012 0.1076714 0.1255855 0.1067292 0.0690872 0.0740393 0.0953817 0.1854669 0.0673989</values>
</curve>
<curve>
<name>c_t_Snh</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>7.96E-02 7.96E-02 0.0699766 0.0750936 7.94E-02 0.0759865 0.0828972 0.0836815 0.0661563 0.0523349 0.0631047 0.0659777 0.072438 0.0523659 0.0526066 4.36E-02 0.0624447 0.0511158 0.0494541 0.0462939 0.053686 0.050184 0.0545944 0.0602861 0.0600376 0.0635007 0.0505101 0.0696194 0.0513876 0.0511003 0.0441663 0.0596649 0.0537869 0.0509217 0.0561862</values>
</curve>
<curve>
<name>c_t_Sno</name>
<coords>0 604800 1209600 1814400 2419200 3024000 3628800 4233600 4838400 5529600 6048000 6652800 7257600 8470800 9075600 10198800 10890000 11667600 12099600 13914000 14518800 15210000 15814800 16333200 16938000 17542800 18147600 18752400 19357200 19962000 20563200 21168000 21772800 22377600 22896000</coords>
<values>7.39E-05 7.39E-05 9.64E-05 9.64E-05 9.64E-05 7.39E-05 7.39E-05 7.39E-05 8.30E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 9.64E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05 7.39E-05</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>p0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_outlet</mesh>
<type>Dirichlet</type>
<parameter>p_Dirichlet_right</parameter>
</boundary_condition>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Neumann</type>
<parameter>p_Neumann_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>H</name>
<components>1</components>
<order>1</order>
<initial_condition>c_H</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_H</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Do</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Do_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sf</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sf_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sa</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sa_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sin</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sin_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Xs_m</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Xsm_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Xi_m</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Xim_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Snh</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Snh_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sno</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sno_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sso</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sso_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>Sulphide</name>
<components>1</components>
<order>1</order>
<initial_condition>c0_CWM1</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>Wetland_1d_inlet</mesh>
<type>Dirichlet</type>
<parameter>c_Sulphide_inlet</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_picard</name>
<type>Picard</type>
<max_iter>10</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>ILUT</precon_type>
<max_iteration_step>1000</max_iteration_step>
<error_tolerance>1e-14</error_tolerance>
</eigen>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
1.11.0