<?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_CO2_i" units="mole"/>
        <variable initial_value="1.1877" name="q_CO2_o" units="mole"/>
        <variable name="v_CO2_m" public_interface="out" units="mol_per_s"/>
        <variable initial_value="1" name="E_CO2_0" units="mole"/>
        <variable initial_value="1e6" name="kappa_CO2_3" units="mol_per_s"/>
        <variable initial_value="1e6" name="kappa_CO2_4" units="mol_per_s"/>
        <variable initial_value="96" name="K_CO2_m" units="per_mol"/>
        <variable initial_value="1e-8" name="K_CO2_s" units="per_mol"/>
        <variable name="t" units="second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>t</ci>
                    </bvar>
                    <ci>q_CO2_i</ci>
                </apply>
                <cn cellml:units="mol_per_s">0.02</cn>
            </apply>
            <apply>
                <eq/>
                <ci>v_CO2_m</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>E_CO2_0</ci>
                        <ci>kappa_CO2_3</ci>
                        <ci>kappa_CO2_4</ci>
                        <ci>K_CO2_m</ci>
                        <ci>K_CO2_s</ci>
                        <apply>
                            <minus/>
                            <ci>q_CO2_o</ci>
                            <ci>q_CO2_i</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <cn cellml:units="dim">2</cn>
                            <ci>kappa_CO2_4</ci>
                        </apply>
                        <apply>
                            <times/>
                            <apply>
                                <plus/>
                                <ci>kappa_CO2_3</ci>
                                <ci>kappa_CO2_4</ci>
                            </apply>
                            <ci>K_CO2_s</ci>
                            <apply>
                                <plus/>
                                <ci>q_CO2_o</ci>
                                <ci>q_CO2_i</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="dim">2</cn>
                            <ci>kappa_CO2_3</ci>
                            <ci>K_CO2_s</ci>
                            <ci>K_CO2_s</ci>
                            <ci>q_CO2_o</ci>
                            <ci>q_CO2_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
</model>
