<?xml version='1.0' encoding='UTF-8'?>
<model name="BG51" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
    <!-- Intracellular pH regulation-->
    <units name="dim">
        <unit units="dimensionless"/>
    </units>
    <units name="per_s">
        <unit exponent="-1" units="second"/>
    </units>
    <units name="mol_per_s">
        <unit units="mole"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_mol">
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="per_mol2">
        <unit exponent="-2" units="mole"/>
    </units>
    <units name="C_per_mol">
        <unit units="coulomb"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="J_per_mol">
        <unit units="joule"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="J_per_C">
        <unit units="joule"/>
        <unit exponent="-1" units="coulomb"/>
    </units>
    <component name="main">
        <variable initial_value="0" name="q_NH3_i" units="mole"/>
        <variable initial_value="0.14" name="q_NH3_o" units="mole"/>
        <variable initial_value="1e-3" name="K_NH3_o" units="per_mol"/>
        <variable initial_value="1e-3" name="K_NH3_i" units="per_mol"/>
        <variable name="Q_NH3_o" units="dim"/>
        <variable name="Q_NH3_i" units="dim"/>
        <variable name="t" units="second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Q_NH3_o</ci>
                <apply>
                    <times/>
                    <ci>K_NH3_o</ci>
                    <ci>q_NH3_o</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Q_NH3_i</ci>
                <apply>
                    <times/>
                    <ci>K_NH3_i</ci>
                    <ci>q_NH3_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_NH3_i</ci>
                </apply>
                <cn cellml:units="mol_per_s">0.01</cn>
            </apply>
        </math>
        <variable name="v_NH3_m" units="mol_per_s"/>
        <variable initial_value="4.801e+02" name="v_NH3_max" units="mol_per_s"/>
        <variable initial_value="1" name="gamma_NH3" units="dim"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_NH3_m</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>v_NH3_max</ci>
                        <apply>
                            <minus/>
                            <ci>Q_NH3_o</ci>
                            <ci>Q_NH3_i</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <cn cellml:units="dim">2</cn>
                                <ci>gamma_NH3</ci>
                            </apply>
                            <apply>
                                <plus/>
                                <cn cellml:units="dim">1</cn>
                                <ci>gamma_NH3</ci>
                            </apply>
                        </apply>
                        <ci>Q_NH3_o</ci>
                        <ci>Q_NH3_i</ci>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dim">2</cn>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="dim">1</cn>
                                    <ci>gamma_NH3</ci>
                                </apply>
                            </apply>
                            <ci>Q_NH3_o</ci>
                            <ci>Q_NH3_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
</model>
