Location: Glucose Uptake in Enterocyte @ c052b0c46028 / SEDML_files / Model.cellml

Author:
nima <nafs080@aucklanduni.ac.nz>
Date:
2020-09-30 10:48:46+13:00
Desc:
Renamed the files
Permanent Source URI:
https://models.fieldml.org/workspace/572/rawfile/c052b0c460280139dad150937fbee4fa6a026505/SEDML_files/Model.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="Nima_2018" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#">
    <!-- Unit and value for K_s (GLUT2) is fixed
 NHE3 and NBC were fixed also sign of J_NHE3 changed from + to -  
 decrease D_P_Na from 1e-8 t0 1e-9
JENaC was fixed (-to+) and the sign for J_Enac and I_Enac changed from - to + 
CFTR was fixed (-to+) nad the sign for J_Enac and I_Enac changed from - to + -->
    <units name="per_second">
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_second2">
        <unit exponent="-1" units="second"/>
    </units>
    <!--    def unit per_second5 as
        unit second {expo: -5};
    enddef;-->
    <units name="volt_per_second">
        <unit units="volt"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_meter">
        <unit exponent="-1" units="meter"/>
    </units>
    <units name="per_volt_per_second">
        <unit exponent="-1" units="volt"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_volt4_per_second">
        <unit exponent="-4" units="volt"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_volt3_per_second">
        <unit exponent="-3" units="volt"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_volt2_per_second">
        <unit exponent="-2" units="volt"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_mol">
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="umol">
        <unit prefix="micro" units="mole"/>
    </units>
    <units name="per_umol">
        <unit exponent="-1" units="umol"/>
    </units>
    <units name="M">
        <unit units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="per_M">
        <unit exponent="-1" units="M"/>
    </units>
    <units name="m2">
        <unit exponent="2" units="metre"/>
    </units>
    <units name="cm">
        <unit prefix="centi" units="metre"/>
    </units>
    <units name="um2">
        <unit exponent="2" prefix="micro" units="metre"/>
    </units>
    <units name="um">
        <unit prefix="micro" units="meter"/>
    </units>
    <units name="um3">
        <unit exponent="3" prefix="micro" units="metre"/>
    </units>
    <units name="m3">
        <unit exponent="3" units="metre"/>
    </units>
    <!--    def unit litre as
        unit m3 {pref: milli};
    enddef;-->
    <units name="mol_per_um2">
        <unit units="mole"/>
        <unit exponent="-1" units="um2"/>
    </units>
    <units name="mM">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="mV">
        <unit prefix="milli" units="volt"/>
    </units>
    <units name="J_per_K_per_mol">
        <unit exponent="-1" units="mole"/>
        <unit units="joule"/>
        <unit exponent="-1" units="kelvin"/>
    </units>
    <units name="pA_per_cm2">
        <unit prefix="pico" units="ampere"/>
        <unit exponent="-2" prefix="centi" units="metre"/>
    </units>
    <units name="uA">
        <unit prefix="micro" units="ampere"/>
    </units>
    <units name="uA_per_um2">
        <unit units="uA"/>
        <unit exponent="-1" units="um2"/>
    </units>
    <units name="uA_per_A">
        <unit units="uA"/>
        <unit exponent="-1" units="ampere"/>
    </units>
    <units name="umol_per_mol">
        <unit units="umol"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="C_per_mol">
        <unit exponent="-1" units="mole"/>
        <unit units="coulomb"/>
    </units>
    <units name="per_M_per_second">
        <unit exponent="-1" units="M"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_M2_per_second">
        <unit exponent="-2" units="M"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_M4_per_second2">
        <unit exponent="-4" units="M"/>
        <unit exponent="-2" units="second"/>
    </units>
    <units name="per_M4_per_second">
        <unit exponent="-4" units="M"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="uM">
        <unit prefix="micro" units="M"/>
    </units>
    <units name="M2">
        <unit exponent="2" units="M"/>
    </units>
    <units name="M3">
        <unit exponent="3" units="M"/>
    </units>
    <units name="M3_per_second">
        <unit exponent="3" units="M"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_M3_per_second5">
        <unit exponent="-3" units="M"/>
        <unit exponent="-5" units="second"/>
    </units>
    <units name="M_per_second">
        <unit units="M"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="cm2">
        <unit exponent="2" prefix="centi" units="metre"/>
    </units>
    <units name="cm2_per_s">
        <unit units="cm2"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="mole_per_m2">
        <unit units="mole"/>
        <unit exponent="-1" units="m2"/>
    </units>
    <!--    def unit uM_per_cm2_per_s as
        unit uM;
        unit cm2 {expo: -1};
        unit second {expo: -1};
    enddef;-->
    <units name="Farad">
        <unit exponent="-1" units="volt"/>
        <unit units="coulomb"/>
    </units>
    <units name="uF">
        <unit prefix="micro" units="Farad"/>
    </units>
    <units name="uF_per_cm2">
        <unit prefix="micro" units="Farad"/>
        <unit exponent="-1" units="cm2"/>
    </units>
    <units name="umol_per_cm2_per_s">
        <unit prefix="micro" units="mole"/>
        <unit exponent="-1" units="cm2"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="umol_per_s">
        <unit prefix="micro" units="mole"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="uA_per_cm2">
        <unit units="uA"/>
        <unit exponent="-2" prefix="centi" units="metre"/>
    </units>
    <units name="m_per_s">
        <unit units="meter"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="m3_mole_per_litre_per_umol">
        <unit units="m3"/>
        <unit units="mole"/>
        <unit exponent="-1" units="litre"/>
        <unit exponent="-1" units="umol"/>
    </units>
    <units name="litre_umol_per_mol_per_m3">
        <unit units="litre"/>
        <unit units="umol"/>
        <unit exponent="-1" units="mole"/>
        <unit exponent="-1" units="m3"/>
    </units>
    <units name="S">
        <unit units="siemens"/>
    </units>
    <units name="S_per_m2">
        <unit units="siemens"/>
        <unit exponent="-1" units="m2"/>
    </units>
    <units name="per_m2">
        <unit exponent="-1" units="m2"/>
    </units>
    <units name="cm_per_s">
        <unit units="cm"/>
        <unit exponent="-1" units="second"/>
    </units>
    <component name="parameters">
        <variable name="time" public_interface="out" units="second"/>
        <variable initial_value="-1" name="z_Cl" public_interface="out" units="dimensionless"/>
        <variable initial_value="1" name="z_Na" public_interface="out" units="dimensionless"/>
        <variable initial_value="1" name="z_k" public_interface="out" units="dimensionless"/>
        <variable initial_value="96485" name="F" public_interface="out" units="C_per_mol"/>
        <variable initial_value="8.314" name="R" public_interface="out" units="J_per_K_per_mol"/>
        <variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
        <variable initial_value="1e-5" name="capacitance" public_interface="out" units="uF"/>
        <variable initial_value="8.4e-10" name="A_Apical" public_interface="out" units="m2"/>
        <variable initial_value="5e-10" name="A_Basol" public_interface="out" units="m2"/>
        <variable initial_value="2e-15" name="v_cell" public_interface="out" units="m3"/>
        <variable initial_value="0" name="P_Na" public_interface="out" units="m_per_s"/>
        <variable initial_value="0" name="P_k" public_interface="out" units="m_per_s"/>
        <variable initial_value="0" name="P_Cl" public_interface="out" units="m_per_s"/>
        <variable initial_value="1e-8" name="P_P_Na" public_interface="out" units="m_per_s"/>
        <variable initial_value="5e-8" name="P_P_k" public_interface="out" units="m_per_s"/>
        <variable initial_value="2.2e-8" name="P_P_Cl" public_interface="out" units="m_per_s"/>
        <variable initial_value="0" name="P_Bl_Na" public_interface="out" units="m_per_s"/>
        <variable initial_value="0" name="P_Bl_k" public_interface="out" units="m_per_s"/>
        <variable initial_value="0" name="P_Bl_Cl" public_interface="out" units="m_per_s"/>
        <variable initial_value="14000" name="k0_12" public_interface="out" units="per_M2_per_second"/>
        <variable initial_value="300" name="k0_21" public_interface="out" units="per_second"/>
        <variable initial_value="15000" name="k0_23" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="20" name="k0_32" public_interface="out" units="per_second"/>
        <variable initial_value="0.1" name="k0_25" public_interface="out" units="per_second"/>
        <variable initial_value="0.01" name="k0_52" public_interface="out" units="per_second"/>
        <variable initial_value="50" name="k0_34" public_interface="out" units="per_second"/>
        <variable initial_value="50" name="k0_43" public_interface="out" units="per_second"/>
        <variable initial_value="6" name="k0_45" public_interface="out" units="per_second"/>
        <variable initial_value="450" name="k0_54" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="10" name="k0_56" public_interface="out" units="per_second"/>
        <variable initial_value="200" name="k0_65" public_interface="out" units="per_M2_per_second"/>
        <variable initial_value="25" name="k0_61" public_interface="out" units="per_second"/>
        <variable initial_value="600" name="k0_16" public_interface="out" units="per_second"/>
        <variable initial_value="0.7" name="delta" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.3" name="alpha_p" public_interface="out" units="dimensionless"/>
        <variable initial_value="0" name="alpha_pp" public_interface="out" units="dimensionless"/>
        <variable initial_value="2" name="n" public_interface="out" units="dimensionless"/>
        <variable initial_value="-2" name="z_c" public_interface="out" units="dimensionless"/>
        <variable name="v_mc" public_interface="in" units="volt"/>
        <variable name="mu" public_interface="out" units="dimensionless"/>
        <!--alpha_p = 1{dimensionless}-delta-alpha_pp;-->
        <!--C_T = 1e6{umol_per_mol}*N_C/N_Avo;-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>mu</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>F</ci>
                        <ci>v_mc</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Apical_concentrations">
        <variable initial_value="0.105" name="Na_m" public_interface="out" units="M"/>
        <variable initial_value="0.004" name="K_m" public_interface="out" units="M"/>
        <variable initial_value="0.05" name="glucose_m" public_interface="out" units="M"/>
        <variable initial_value="0.106" name="Cl_m" public_interface="out" units="M"/>
        <!--        var k: per_second {init: 0.2, pub: out};
        var glucose_0: M {init: 20e-3, pub: out};
        var glucose_f: M {init: 0.1e-3, pub: out};
        var time: second {pub: in};-->
        <!--        glucose_m = sel
            case (time >= 0{second}) and (time < 60{second}):
                0.1e-3{M};
            case time >= 60{second}:
                (glucose_0-glucose_f)*exp(-k*(time-60{second}))+glucose_f;
        endsel;-->
    </component>
    <component name="Apical_voltage">
        <variable initial_value="-0.046" name="v_mc" public_interface="out" units="volt"/>
        <variable name="time" public_interface="in" units="second"/>
        <variable name="capacitance" public_interface="in" units="uF"/>
        <variable name="I_D_P_Na" public_interface="in" units="uA"/>
        <variable name="I_D_P_K" public_interface="in" units="uA"/>
        <variable name="I_D_P_Cl" public_interface="in" units="uA"/>
        <variable name="I_NHE3_H" public_interface="in" units="uA"/>
        <variable name="I_NHE3_Na" public_interface="in" units="uA"/>
        <variable name="I_BK" public_interface="in" units="uA"/>
        <!--        var I_AE1_HCO3: uA {pub: in};
        var I_AE1_Cl: uA {pub: in};-->
        <variable name="I_ENaC" public_interface="in" units="uA"/>
        <!--var I_NaCl: uA {pub: in};-->
        <variable name="I_Na_SGLT" public_interface="in" units="uA"/>
        <variable name="I_CFTR" public_interface="in" units="uA"/>
        <!--var I_CFTR_HCO3: uA {pub: in};-->
        <!--var I_Gl_SGLT: uA {pub: in};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>v_mc</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>I_Na_SGLT</ci>
                            </apply>
                            <ci>I_BK</ci>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>I_CFTR</ci>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>I_ENaC</ci>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <plus/>
                                    <ci>I_D_P_Na</ci>
                                    <ci>I_D_P_K</ci>
                                    <ci>I_D_P_Cl</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>capacitance</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Cell_concentration">
        <variable initial_value="0.047" name="Na_i" public_interface="out" units="M"/>
        <variable initial_value="0.003" name="glucose_i" public_interface="out" units="M"/>
        <variable initial_value="0.129" name="K_i" public_interface="out" units="M"/>
        <variable initial_value="6e-5" name="Ca_i" public_interface="out" units="M"/>
        <variable initial_value="0.057" name="Cl_i" public_interface="out" units="M"/>
        <variable name="HCO3_i" public_interface="out" units="M"/>
        <variable name="H_int" public_interface="out" units="M"/>
        <variable initial_value="0.003" name="HCO3_m" public_interface="out" units="M"/>
        <variable initial_value="0.041" name="HCO3_s" public_interface="out" units="M"/>
        <variable name="H_Ext" public_interface="out" units="M"/>
        <variable name="H_s" public_interface="out" units="M"/>
        <variable name="v_cell" public_interface="in" units="m3"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="J_Na_SGLT" public_interface="in" units="umol_per_s"/>
        <variable name="J_Gl_SGLT" public_interface="in" units="umol_per_s"/>
        <variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
        <variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
        <variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_NBC_HCO3" public_interface="in" units="umol_per_s"/>
        <!--var J_NaCl: umol_per_s {pub: in};-->
        <variable name="J_NHE3_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_NHE3_H" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_Bl_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
        <variable name="J_A_GLUT" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_K" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_Bl_K" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_Cl" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_Bl_Cl" public_interface="in" units="umol_per_s"/>
        <variable name="J_AE1_Cl" public_interface="in" units="umol_per_s"/>
        <variable name="J_AE1_HCO3" public_interface="in" units="umol_per_s"/>
        <variable name="J_CFTR" public_interface="in" units="umol_per_s"/>
        <variable name="J_CLC2" public_interface="in" units="umol_per_s"/>
        <variable name="J_CFTR_HCO3" public_interface="in" units="umol_per_s"/>
        <variable name="J_ENaC" public_interface="in" units="umol_per_s"/>
        <variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
        <variable name="J_IK" public_interface="in" units="umol_per_s"/>
        <variable name="J_BK" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
        <!--        var J_buf_C: M_per_second {pub: in};
        var J_buf_A: M_per_second {pub: in};-->
        <variable name="time" public_interface="in" units="second"/>
        <variable initial_value="7.2" name="pH_int" public_interface="out" units="per_second"/>
        <variable initial_value="7.2" name="pH_Ext" public_interface="out" units="per_second"/>
        <variable initial_value="7.2" name="pH_s" public_interface="out" units="per_second"/>
        <variable initial_value="1680000" name="rho" public_interface="out" units="per_meter"/>
        <variable name="beta" public_interface="out" units="M"/>
        <variable initial_value="6.3" name="p_K" public_interface="out" units="dimensionless"/>
        <variable initial_value="40" name="P_CO2" public_interface="out" units="dimensionless"/>
        <variable initial_value="1e-5" name="s" public_interface="out" units="dimensionless"/>
        <variable initial_value="1" name="theta_1" units="dimensionless"/>
        <variable initial_value="1" name="theta_2" units="dimensionless"/>
        <variable initial_value="1" name="theta_3" units="dimensionless"/>
        <variable initial_value="1" name="theta_4" units="dimensionless"/>
        <variable initial_value="1" name="theta_5" units="dimensionless"/>
        <variable initial_value="1" name="theta_6" units="dimensionless"/>
        <variable initial_value="1" name="theta_7" units="dimensionless"/>
        <variable initial_value="1" name="theta_8" units="dimensionless"/>
        <variable initial_value="1" name="theta_9" units="dimensionless"/>
        <variable initial_value="1" name="theta_10" units="dimensionless"/>
        <variable initial_value="1" name="theta_11" units="dimensionless"/>
        <variable initial_value="1" name="theta_12" units="dimensionless"/>
        <variable initial_value="1" name="theta_13" units="dimensionless"/>
        <variable initial_value="1" name="theta_14" units="dimensionless"/>
        <variable initial_value="1" name="theta_15" units="dimensionless"/>
        <variable initial_value="1" name="theta_16" units="dimensionless"/>
        <variable initial_value="1" name="theta_17" units="dimensionless"/>
        <variable initial_value="1" name="theta_18" units="dimensionless"/>
        <variable initial_value="1" name="theta_19" units="dimensionless"/>
        <variable initial_value="1" name="theta_20" units="dimensionless"/>
        <variable initial_value="1" name="theta_21" units="dimensionless"/>
        <variable initial_value="1" name="theta_22" units="dimensionless"/>
        <variable initial_value="1" name="theta_23" units="dimensionless"/>
        <variable initial_value="1" name="theta_24" units="dimensionless"/>
        <variable initial_value="1" name="theta_25" units="dimensionless"/>
        <variable initial_value="1" name="theta_26" units="dimensionless"/>
        <variable initial_value="1" name="theta_27" units="dimensionless"/>
        <variable name="K_Concentration" public_interface="out" units="M"/>
        <variable name="Na_Concentration" public_interface="out" units="M"/>
        <variable name="Cl_Concentration" public_interface="out" units="M"/>
        <variable name="Total_Gl_Flux_Apical" public_interface="out" units="umol_per_s"/>
        <!--var Concentration_Na: M {pub: out};-->
        <!-- Concentration_Na = Na_m+Na_i+Na_s;-->
        <!--var K_g: per_second {init: 0.0015, pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Total_Gl_Flux_Apical</ci>
                <apply>
                    <plus/>
                    <ci>J_A_GLUT</ci>
                    <ci>J_Gl_SGLT</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Na_Concentration</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <apply>
                                <plus/>
                                <ci>theta_15</ci>
                            </apply>
                            <ci>J_ENaC</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_3</ci>
                            <ci>J_NHE3_Na</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K_Concentration</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">1</cn>
                                    <ci>J_K_NaK</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">1</cn>
                                    <ci>J_NaKCC</ci>
                                </apply>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>J_D_K</ci>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>J_D_Bl_K</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0</cn>
                            <ci>J_IK</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Cl_Concentration</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    <apply>
                        <minus/>
                        <apply>
                            <times/>
                            <apply>
                                <plus/>
                                <ci>theta_11</ci>
                            </apply>
                            <ci>J_CFTR</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_14</ci>
                            <ci>J_AE1_Cl</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--var K_g: per_second {init: 0.0015, pub: out};-->
            <apply>
                <eq/>
                <ci>beta</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2.3</cn>
                    <ci>HCO3_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Na_i</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
                        <ci>v_cell</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <minus/>
                            <apply>
                                <times/>
                                <ci>theta_1</ci>
                                <ci>J_Na_SGLT</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>theta_2</ci>
                                <ci>J_Na_NaK</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_15</ci>
                            <ci>J_ENaC</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_16</ci>
                            <ci>J_NaKCC</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_19</ci>
                            <ci>J_NBC_Na</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_4</ci>
                            <ci>J_D_Na</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_5</ci>
                            <ci>J_D_Bl_Na</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_3</ci>
                            <ci>J_NHE3_Na</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>glucose_i</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
                        <ci>v_cell</ci>
                    </apply>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>theta_26</ci>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>J_A_GLUT</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>theta_6</ci>
                                <ci>J_Gl_SGLT</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_7</ci>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>J_GLUT</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>K_i</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
                        <ci>v_cell</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>theta_8</ci>
                            <ci>J_K_NaK</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_17</ci>
                            <ci>J_NaKCC</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_9</ci>
                            <ci>J_D_K</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_10</ci>
                            <ci>J_D_Bl_K</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_27</ci>
                            <ci>J_IK</ci>
                        </apply>
                        <ci>J_BK</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Cl_i</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
                        <ci>v_cell</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <minus/>
                            <apply>
                                <times/>
                                <apply>
                                    <minus/>
                                    <ci>theta_11</ci>
                                </apply>
                                <ci>J_CFTR</ci>
                            </apply>
                            <ci>J_CLC2</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_18</ci>
                            <cn cellml:units="dimensionless">2</cn>
                            <ci>J_NaKCC</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_12</ci>
                            <ci>J_D_Cl</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_13</ci>
                            <ci>J_D_Bl_Cl</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_14</ci>
                            <ci>J_AE1_Cl</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--ode(HCO3_i, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*theta_20*J_NBC_HCO3+theta_25*J_AE1_HCO3+theta_21*J_buf_C;-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>pH_int</ci>
                </apply>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
                            <ci>rho</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>beta</ci>
                            <ci>A_Apical</ci>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <apply>
                                <times/>
                                <ci>theta_20</ci>
                                <ci>J_NBC_HCO3</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>theta_22</ci>
                                <ci>J_NHE3_H</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>theta_25</ci>
                            <ci>J_AE1_HCO3</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--ode(pH_Ext, time) = 1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*(theta_22*J_NHE3_H-theta_25*J_AE1_HCO3-theta_26*J_CFTR_HCO3);-->
            <!--ode(pH_s, time) = -1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*theta_20*J_NBC_HCO3;-->
            <!--ode(HCO3_m, time) = theta_24*-J_AE1_HCO3+theta_23*J_buf_A;-->
            <apply>
                <eq/>
                <ci>H_int</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <apply>
                        <minus/>
                        <ci>pH_int</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>HCO3_i</ci>
                <apply>
                    <times/>
                    <ci>s</ci>
                    <ci>P_CO2</ci>
                    <apply>
                        <power/>
                        <cn cellml:units="dimensionless">10</cn>
                        <apply>
                            <minus/>
                            <ci>pH_int</ci>
                            <ci>p_K</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>H_Ext</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <apply>
                        <minus/>
                        <ci>pH_Ext</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>H_s</ci>
                <apply>
                    <power/>
                    <cn cellml:units="dimensionless">10</cn>
                    <apply>
                        <minus/>
                        <ci>pH_s</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--        HCO3_m = s*P_CO2*pow(10{dimensionless}, pH_Ext-p_K);
        HCO3_s = s*P_CO2*pow(10{dimensionless}, pH_s-p_K);-->
    </component>
    <component name="Basolateral_concentrations">
        <variable initial_value="0.155" name="Na_s" public_interface="out" units="M"/>
        <variable initial_value="0.005" name="glucose_s" public_interface="out" units="M"/>
        <variable initial_value="0.119" name="Cl_s" public_interface="out" units="M"/>
        <variable initial_value="0.0054" name="K_s" public_interface="out" units="M"/>
        <variable name="v_cell" public_interface="in" units="m3"/>
        <variable name="time" public_interface="in" units="second"/>
        <variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
        <variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
        <variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
        <variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_P_K" public_interface="in" units="umol_per_s"/>
        <variable name="J_IK" public_interface="in" units="umol_per_s"/>
        <variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
        <variable name="J_D_P_Cl" public_interface="in" units="umol_per_s"/>
        <variable name="J_CLC2" public_interface="in" units="umol_per_s"/>
        <variable name="v_b" public_interface="out" units="m3"/>
        <variable initial_value="1" name="m" public_interface="out" units="dimensionless"/>
        <variable initial_value="0" name="theta" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_b</ci>
                <apply>
                    <times/>
                    <ci>m</ci>
                    <ci>v_cell</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Cl_s</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>theta</ci>
                        <cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
                        <apply>
                            <plus/>
                            <apply>
                                <minus/>
                                <ci>J_CLC2</ci>
                                <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">2</cn>
                                    <ci>J_NaKCC</ci>
                                </apply>
                            </apply>
                            <ci>J_D_P_Cl</ci>
                        </apply>
                    </apply>
                    <ci>v_b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>K_s</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>theta</ci>
                        <cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
                        <apply>
                            <minus/>
                            <apply>
                                <minus/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <minus/>
                                        <ci>J_K_NaK</ci>
                                    </apply>
                                    <ci>J_NaKCC</ci>
                                </apply>
                                <ci>J_IK</ci>
                            </apply>
                            <ci>J_D_P_K</ci>
                        </apply>
                    </apply>
                    <ci>v_b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Na_s</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>theta</ci>
                        <cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
                        <apply>
                            <minus/>
                            <apply>
                                <minus/>
                                <apply>
                                    <minus/>
                                    <ci>J_Na_NaK</ci>
                                    <ci>J_NaKCC</ci>
                                </apply>
                                <ci>J_NBC_Na</ci>
                            </apply>
                            <ci>J_D_P_Na</ci>
                        </apply>
                    </apply>
                    <ci>v_b</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>glucose_s</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>theta</ci>
                        <cn cellml:units="mole_per_umol" type="e-notation">1<sep/>-9</cn>
                        <ci>J_GLUT</ci>
                    </apply>
                    <ci>v_b</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Basol_voltage">
        <variable initial_value="-0.05" name="v_sc" public_interface="out" units="volt"/>
        <variable name="time" public_interface="in" units="second"/>
        <variable name="capacitance" public_interface="in" units="uF"/>
        <variable name="I_Na_NaK" public_interface="in" units="uA"/>
        <variable name="I_K_NaK" public_interface="in" units="uA"/>
        <variable name="I_D_P_Na" public_interface="in" units="uA"/>
        <variable name="I_D_P_K" public_interface="in" units="uA"/>
        <variable name="I_D_P_Cl" public_interface="in" units="uA"/>
        <variable name="I_NBC_Na" public_interface="in" units="uA"/>
        <variable name="I_NBC_HCO3" public_interface="in" units="uA"/>
        <variable name="I_IK" public_interface="in" units="uA"/>
        <variable name="I_CLC2" public_interface="in" units="uA"/>
        <!--var I_GLUT: uA {pub: in};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>v_sc</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <plus/>
                            <apply>
                                <minus/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <minus/>
                                        <apply>
                                            <minus/>
                                            <apply>
                                                <plus/>
                                                <ci>I_CLC2</ci>
                                                <apply>
                                                    <minus/>
                                                    <ci>I_Na_NaK</ci>
                                                </apply>
                                                <ci>I_K_NaK</ci>
                                                <apply>
                                                    <times/>
                                                    <cn cellml:units="dimensionless">1</cn>
                                                    <ci>I_NBC_Na</ci>
                                                </apply>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <cn cellml:units="dimensionless">1</cn>
                                                <ci>I_NBC_HCO3</ci>
                                            </apply>
                                        </apply>
                                        <apply>
                                            <times/>
                                            <cn cellml:units="dimensionless">1</cn>
                                            <ci>I_D_P_Na</ci>
                                        </apply>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <cn cellml:units="dimensionless">1</cn>
                                        <ci>I_D_P_K</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <times/>
                                    <cn cellml:units="dimensionless">1</cn>
                                    <ci>I_D_P_Cl</ci>
                                </apply>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>I_IK</ci>
                            </apply>
                        </apply>
                    </apply>
                    <ci>capacitance</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="rate_constants">
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="glucose_i" public_interface="in" units="M"/>
        <variable name="glucose_m" public_interface="in" units="M"/>
        <variable name="k_12" public_interface="out" units="per_second"/>
        <variable name="ks_12" public_interface="out" units="per_M2_per_second"/>
        <variable name="k_21" public_interface="out" units="per_second"/>
        <variable name="k_32" public_interface="out" units="per_second"/>
        <variable name="k_25" public_interface="out" units="per_second"/>
        <variable name="k_52" public_interface="out" units="per_second"/>
        <variable name="k_34" public_interface="out" units="per_second"/>
        <variable name="k_43" public_interface="out" units="per_second"/>
        <variable name="k_45" public_interface="out" units="per_second"/>
        <variable name="k_54" public_interface="out" units="per_second"/>
        <variable name="k_56" public_interface="out" units="per_second"/>
        <variable name="k_65" public_interface="out" units="per_second"/>
        <variable name="k_61" public_interface="out" units="per_second"/>
        <variable name="k_16" public_interface="out" units="per_second"/>
        <variable name="k0_12" public_interface="in" units="per_M2_per_second"/>
        <variable name="k0_21" public_interface="in" units="per_second"/>
        <variable name="k0_32" public_interface="in" units="per_second"/>
        <variable name="k0_25" public_interface="in" units="per_second"/>
        <variable name="k0_52" public_interface="in" units="per_second"/>
        <variable name="k0_34" public_interface="in" units="per_second"/>
        <variable name="k0_43" public_interface="in" units="per_second"/>
        <variable name="k0_45" public_interface="in" units="per_second"/>
        <variable name="k0_54" public_interface="in" units="per_M_per_second"/>
        <variable name="k0_56" public_interface="in" units="per_second"/>
        <variable name="k0_65" public_interface="in" units="per_M2_per_second"/>
        <variable name="k0_61" public_interface="in" units="per_second"/>
        <variable name="k0_16" public_interface="in" units="per_second"/>
        <variable name="delta" public_interface="in" units="dimensionless"/>
        <variable name="n" public_interface="in" units="dimensionless"/>
        <variable name="z_c" public_interface="in" units="dimensionless"/>
        <variable name="z_Na" public_interface="in" units="dimensionless"/>
        <variable name="alpha_p" public_interface="in" units="dimensionless"/>
        <variable name="mu" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>ks_12</ci>
                <apply>
                    <times/>
                    <ci>k0_12</ci>
                    <apply>
                        <exp/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <minus/>
                                    <ci>n</ci>
                                </apply>
                                <ci>alpha_p</ci>
                                <ci>mu</ci>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_12</ci>
                <apply>
                    <times/>
                    <ci>k0_12</ci>
                    <apply>
                        <exp/>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <ci>alpha_p</ci>
                            </apply>
                            <ci>mu</ci>
                        </apply>
                    </apply>
                    <apply>
                        <power/>
                        <ci>Na_m</ci>
                        <cn cellml:units="dimensionless">2</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_21</ci>
                <apply>
                    <times/>
                    <ci>k0_21</ci>
                    <apply>
                        <exp/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>n</ci>
                                <ci>z_Na</ci>
                                <ci>alpha_p</ci>
                                <ci>mu</ci>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_32</ci>
                <ci>k0_32</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_34</ci>
                <ci>k0_34</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_43</ci>
                <ci>k0_43</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_45</ci>
                <ci>k0_45</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_54</ci>
                <apply>
                    <times/>
                    <ci>k0_54</ci>
                    <ci>glucose_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_25</ci>
                <ci>k0_25</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_56</ci>
                <ci>k0_56</ci>
            </apply>
            <apply>
                <eq/>
                <ci>k_65</ci>
                <apply>
                    <times/>
                    <ci>k0_65</ci>
                    <apply>
                        <power/>
                        <ci>Na_i</ci>
                        <cn cellml:units="dimensionless">2</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_61</ci>
                <apply>
                    <times/>
                    <ci>k0_61</ci>
                    <apply>
                        <exp/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>z_c</ci>
                                <ci>delta</ci>
                                <ci>mu</ci>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_16</ci>
                <apply>
                    <times/>
                    <ci>k0_16</ci>
                    <apply>
                        <exp/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <minus/>
                                    <ci>z_c</ci>
                                </apply>
                                <ci>delta</ci>
                                <ci>mu</ci>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>k_52</ci>
                <ci>k0_52</ci>
            </apply>
        </math>
        <!--k0_54_temp = k0_23*k_34*k_45/(k_43*k_32);-->
    </component>
    <component name="phenomonological_constants">
        <variable name="epsilon" public_interface="out" units="per_second"/>
        <variable name="lambda" public_interface="out" units="per_M3_per_second5"/>
        <variable name="chi" public_interface="out" units="M"/>
        <variable name="alpha" public_interface="out" units="M3"/>
        <variable name="beta" public_interface="out" units="M2"/>
        <variable name="gamma" public_interface="out" units="M3_per_second"/>
        <variable name="phi" public_interface="out" units="M_per_second"/>
        <!--var k_12: per_second {pub: in};-->
        <variable name="k_21" public_interface="in" units="per_second"/>
        <!--var k_23: per_second {pub: in};-->
        <variable name="k_32" public_interface="in" units="per_second"/>
        <variable name="k_34" public_interface="in" units="per_second"/>
        <variable name="k_43" public_interface="in" units="per_second"/>
        <variable name="k_45" public_interface="in" units="per_second"/>
        <variable name="k_54" public_interface="in" units="per_second"/>
        <variable name="k_25" public_interface="in" units="per_second"/>
        <variable name="k_52" public_interface="in" units="per_second"/>
        <variable name="k_56" public_interface="in" units="per_second"/>
        <variable name="k_65" public_interface="in" units="per_second"/>
        <variable name="k_61" public_interface="in" units="per_second"/>
        <variable name="k_16" public_interface="in" units="per_second"/>
        <variable name="ks_12" public_interface="in" units="per_M2_per_second"/>
        <variable name="k0_23" public_interface="in" units="per_M_per_second"/>
        <!--var C_T: umol {pub: in};-->
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="glucose_m" public_interface="in" units="M"/>
        <variable initial_value="40000000" name="n_SGLT" public_interface="out" units="dimensionless"/>
        <variable name="J_Na_SGLT" public_interface="out" units="umol_per_s"/>
        <variable name="I_Na_SGLT" public_interface="out" units="uA"/>
        <variable name="J_Gl_SGLT" public_interface="out" units="umol_per_s"/>
        <variable name="J_GL_SGLT1" public_interface="out" units="umol_per_s"/>
        <!--var I_Gl_SGLT: uA {pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>lambda</ci>
                <apply>
                    <times/>
                    <ci>ks_12</ci>
                    <ci>k0_23</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>chi</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>lambda</ci>
                    </apply>
                    <ci>ks_12</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>beta</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>lambda</ci>
                    </apply>
                    <ci>k0_23</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_34</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>alpha</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>lambda</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_54</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_65</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                            <ci>k_61</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_56</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_16</ci>
                            <ci>k_21</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gamma</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>lambda</ci>
                    </apply>
                    <ci>k_16</ci>
                    <ci>k_21</ci>
                    <ci>k_65</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_54</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                            <ci>k_52</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                            <ci>k_52</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>phi</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>lambda</ci>
                    </apply>
                    <ci>ks_12</ci>
                    <ci>k_56</ci>
                    <ci>k_61</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_34</ci>
                            <ci>k_45</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_45</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_25</ci>
                            <ci>k_32</ci>
                            <ci>k_43</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>epsilon</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                        </apply>
                        <ci>lambda</ci>
                    </apply>
                    <ci>ks_12</ci>
                    <ci>k0_23</ci>
                    <ci>k_34</ci>
                    <ci>k_45</ci>
                    <ci>k_56</ci>
                    <ci>k_61</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_Na_SGLT</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <minus/>
                                    <cn cellml:units="dimensionless">2</cn>
                                </apply>
                                <ci>n_SGLT</ci>
                            </apply>
                            <cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>epsilon</ci>
                                <apply>
                                    <power/>
                                    <ci>Na_m</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                </apply>
                                <ci>glucose_m</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>phi</ci>
                                <apply>
                                    <power/>
                                    <ci>Na_m</ci>
                                    <cn cellml:units="dimensionless">2</cn>
                                </apply>
                            </apply>
                            <ci>gamma</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>alpha</ci>
                        <apply>
                            <times/>
                            <ci>beta</ci>
                            <ci>glucose_m</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>chi</ci>
                            <apply>
                                <power/>
                                <ci>Na_m</ci>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <apply>
                                <power/>
                                <ci>Na_m</ci>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                            <ci>glucose_m</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_Na_SGLT</ci>
                <apply>
                    <times/>
                    <ci>J_Na_SGLT</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_Gl_SGLT</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <cn cellml:units="dimensionless">2</cn>
                    </apply>
                    <ci>J_Na_SGLT</ci>
                </apply>
            </apply>
            <!--to change the unit to picomol/hour-->
            <apply>
                <eq/>
                <ci>J_GL_SGLT1</ci>
                <apply>
                    <times/>
                    <ci>J_Gl_SGLT</ci>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                </apply>
            </apply>
        </math>
        <!--I_Gl_SGLT = -J_Gl_SGLT*F;-->
    </component>
    <component name="CFTR">
        <variable name="z_Cl" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="Cl_m" public_interface="in" units="M"/>
        <variable name="V_A_Cl" public_interface="out" units="volt"/>
        <variable name="V_A_HCO3" public_interface="out" units="volt"/>
        <variable name="v_mc" public_interface="in" units="volt"/>
        <variable name="J_CFTR" public_interface="out" units="umol_per_s"/>
        <!--        var J_CFTR1: umol_per_s {pub: out};-->
        <variable name="J_CFTR_HCO3" public_interface="out" units="umol_per_s"/>
        <variable name="I_CFTR" public_interface="out" units="uA"/>
        <variable name="I_CFTR_HCO3" public_interface="out" units="uA"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="1" name="G_CFTR" public_interface="out" units="S_per_m2"/>
        <variable name="HCO3_i" public_interface="in" units="M"/>
        <variable name="HCO3_m" public_interface="in" units="M"/>
        <variable initial_value="-1" name="z_HCO3" public_interface="out" units="dimensionless"/>
        <!-- alpha is multipled by area in order to have a dimensionless unit-->
        <!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};
        J_CFTR1 = J_CFTR*36e8{dimensionless};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_CFTR</ci>
                <apply>
                    <divide/>
                    <ci>I_CFTR</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_A_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_Cl</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Cl_m</ci>
                            <ci>Cl_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_CFTR</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <ci>G_CFTR</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_mc</ci>
                        <ci>V_A_Cl</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_A_HCO3</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_HCO3</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>HCO3_m</ci>
                            <ci>HCO3_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_CFTR_HCO3</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless" type="e-notation">0.25<sep/>6</cn>
                    <ci>G_CFTR</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_mc</ci>
                        <ci>V_A_HCO3</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_CFTR_HCO3</ci>
                <apply>
                    <divide/>
                    <ci>I_CFTR_HCO3</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="CLC2">
        <variable name="z_Cl" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="Cl_s" public_interface="in" units="M"/>
        <variable name="V_A_Cl" public_interface="out" units="volt"/>
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="J_CLC2" public_interface="out" units="umol_per_s"/>
        <variable name="I_CLC2" public_interface="out" units="uA"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="1" name="G_CLC2" public_interface="out" units="S_per_m2"/>
        <!-- alpha is multipled by area in order to have a dimensionless unit-->
        <!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_CLC2</ci>
                <apply>
                    <divide/>
                    <ci>I_CLC2</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_A_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_Cl</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Cl_s</ci>
                            <ci>Cl_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_CLC2</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <ci>G_CLC2</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_sc</ci>
                        <ci>V_A_Cl</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="ENaC">
        <variable name="z_Na" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="V_A_Na" public_interface="out" units="volt"/>
        <variable name="v_mc" public_interface="in" units="volt"/>
        <variable name="J_ENaC" public_interface="out" units="umol_per_s"/>
        <variable name="J_ENaC1" public_interface="out" units="umol_per_s"/>
        <variable name="I_ENaC" public_interface="out" units="uA"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="2" name="G_ENaC" public_interface="out" units="S_per_m2"/>
        <!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
        <!--var alpha_ENaC: per_m2 {init: 0.1, pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_ENaC</ci>
                <apply>
                    <divide/>
                    <ci>I_ENaC</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_ENaC1</ci>
                <apply>
                    <times/>
                    <ci>J_ENaC</ci>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_A_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_Na</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Na_m</ci>
                            <ci>Na_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_ENaC</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
                    </apply>
                    <ci>G_ENaC</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_mc</ci>
                        <ci>V_A_Na</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="NaKCC">
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="S_C" public_interface="out" units="M"/>
        <variable name="J_NaKCC" public_interface="out" units="umol_per_s"/>
        <variable name="I_NaKCC" public_interface="out" units="uA"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="K_i" public_interface="in" units="M"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="u_NaKCC" public_interface="out" units="per_second"/>
        <variable initial_value="157.55" name="gamma_1" public_interface="out" units="per_second2"/>
        <variable initial_value="20000000" name="gamma_2" public_interface="out" units="per_M4_per_second2"/>
        <variable initial_value="1.03" name="gamma_3" public_interface="out" units="per_second"/>
        <variable initial_value="1380000" name="gamma_4" public_interface="out" units="per_M4_per_second"/>
        <variable initial_value="1e-9" name="alpha_NaKCC" public_interface="out" units="mole_per_m2"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>S_C</ci>
                <apply>
                    <times/>
                    <ci>Na_i</ci>
                    <ci>Cl_i</ci>
                    <ci>Cl_i</ci>
                    <ci>K_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>u_NaKCC</ci>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>gamma_1</ci>
                        <apply>
                            <times/>
                            <ci>gamma_2</ci>
                            <ci>S_C</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>gamma_3</ci>
                        <apply>
                            <times/>
                            <ci>gamma_4</ci>
                            <ci>S_C</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NaKCC</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    <ci>alpha_NaKCC</ci>
                    <ci>u_NaKCC</ci>
                    <ci>A_Apical</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NaKCC</ci>
                <apply>
                    <times/>
                    <ci>J_NaKCC</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="A_GLUT2">
        <variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
        <variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
        <variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
        <variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
        <variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
        <variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
        <!--        var C_1:  M {init: 1e-3, pub: out};
        var C_2:  M {init: 0, pub: out};-->
        <variable name="glucose_i" public_interface="in" units="M"/>
        <variable name="glucose_m" public_interface="in" units="M"/>
        <!--        var CS_1:  M {pub: out};
        var CS_2:  M {pub: out};-->
        <variable name="K_S" public_interface="out" units="M"/>
        <variable name="J_A_GLUT" public_interface="out" units="umol_per_s"/>
        <variable name="C_tot" public_interface="out" units="umol"/>
        <variable initial_value="420000000" name="n_GLUT" public_interface="out" units="dimensionless"/>
        <variable name="D" public_interface="out" units="per_second"/>
        <variable name="J_A_GLUT1" public_interface="out" units="umol_per_s"/>
        <!--var time: second {init: 0, pub: out};-->
        <variable name="v_cell" public_interface="in" units="m3"/>
        <!--var I_GLUT: uA {pub: out};-->
        <!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
        <!--        D = sel
            case (glucose_m > 0.02{M}):
                (1{dimensionless}+glucose_i/K_S)*(K_32*glucose_m/K_S+K_41)+(1{dimensionless}+glucose_m/K_S)*(K_23*glucose_i/K_S+K_14);
            case glucose_m < 0.02{M}:
                1e100{per_second};
        endsel;-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>C_tot</ci>
                <apply>
                    <divide/>
                    <ci>n_GLUT</ci>
                    <cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K_S</ci>
                <apply>
                    <divide/>
                    <ci>K_21</ci>
                    <ci>K_12</ci>
                </apply>
            </apply>
            <!--        CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
        CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
            <apply>
                <eq/>
                <ci>J_A_GLUT</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <ci>C_tot</ci>
                        <ci>D</ci>
                    </apply>
                    <apply>
                        <minus/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>K_32</ci>
                                <ci>K_14</ci>
                                <ci>glucose_m</ci>
                            </apply>
                            <ci>K_S</ci>
                        </apply>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>K_23</ci>
                                <ci>K_41</ci>
                                <ci>glucose_i</ci>
                            </apply>
                            <ci>K_S</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_A_GLUT1</ci>
                <apply>
                    <times/>
                    <ci>J_A_GLUT</ci>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>D</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <divide/>
                                <ci>glucose_m</ci>
                                <ci>K_S</ci>
                            </apply>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>K_23</ci>
                                    <ci>glucose_i</ci>
                                </apply>
                                <ci>K_S</ci>
                            </apply>
                            <ci>K_14</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <divide/>
                                <ci>glucose_i</ci>
                                <ci>K_S</ci>
                            </apply>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>K_32</ci>
                                    <ci>glucose_m</ci>
                                </apply>
                                <ci>K_S</ci>
                            </apply>
                            <ci>K_41</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="GLUT2">
        <variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
        <variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
        <variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
        <variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
        <variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
        <variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
        <variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
        <!--        var C_1:  M {init: 1e-3, pub: out};
        var C_2:  M {init: 0, pub: out};-->
        <variable name="glucose_i" public_interface="in" units="M"/>
        <variable name="glucose_s" public_interface="in" units="M"/>
        <!--        var CS_1:  M {pub: out};
        var CS_2:  M {pub: out};-->
        <variable name="K_S" public_interface="out" units="M"/>
        <variable name="J_GLUT" public_interface="out" units="umol_per_s"/>
        <variable name="C_tot" public_interface="out" units="umol"/>
        <variable initial_value="140000000" name="n_GLUT" public_interface="out" units="dimensionless"/>
        <variable name="D" public_interface="out" units="per_second"/>
        <variable name="J_GLUT1" public_interface="out" units="umol_per_s"/>
        <!--var time: second {init: 0, pub: out};-->
        <variable name="v_cell" public_interface="in" units="m3"/>
        <!--var I_GLUT: uA {pub: out};-->
        <!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>C_tot</ci>
                <apply>
                    <divide/>
                    <ci>n_GLUT</ci>
                    <cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K_S</ci>
                <apply>
                    <divide/>
                    <ci>K_21</ci>
                    <ci>K_12</ci>
                </apply>
            </apply>
            <!--        CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
        CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
            <apply>
                <eq/>
                <ci>J_GLUT</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <ci>C_tot</ci>
                        <ci>D</ci>
                    </apply>
                    <apply>
                        <minus/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>K_23</ci>
                                <ci>K_41</ci>
                                <ci>glucose_i</ci>
                            </apply>
                            <ci>K_S</ci>
                        </apply>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <ci>K_32</ci>
                                <ci>K_14</ci>
                                <ci>glucose_s</ci>
                            </apply>
                            <ci>K_S</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_GLUT1</ci>
                <apply>
                    <times/>
                    <ci>J_GLUT</ci>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>D</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <divide/>
                                <ci>glucose_i</ci>
                                <ci>K_S</ci>
                            </apply>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>K_32</ci>
                                    <ci>glucose_s</ci>
                                </apply>
                                <ci>K_S</ci>
                            </apply>
                            <ci>K_41</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <divide/>
                                <ci>glucose_s</ci>
                                <ci>K_S</ci>
                            </apply>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>K_23</ci>
                                    <ci>glucose_i</ci>
                                </apply>
                                <ci>K_S</ci>
                            </apply>
                            <ci>K_14</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="NaK_channel">
        <variable name="time" public_interface="in" units="second"/>
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="K_s" public_interface="in" units="M"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="TO_NaK" public_interface="out" units="per_second"/>
        <variable initial_value="5000000" name="n_NaK" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.01" name="K_Na_i_half" public_interface="out" units="M"/>
        <variable initial_value="0.0015" name="K_K_s_half" public_interface="out" units="M"/>
        <variable name="I_Na_NaK" public_interface="out" units="uA"/>
        <variable name="I_K_NaK" public_interface="out" units="uA"/>
        <variable name="J_NaK" public_interface="out" units="umol_per_s"/>
        <variable name="J_Na_NaK" public_interface="out" units="umol_per_s"/>
        <variable name="J_K_NaK" public_interface="out" units="umol_per_s"/>
        <!--var J_K_NaK: umol_per_s {pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_K_NaK</ci>
                <apply>
                    <times/>
                    <ci>J_K_NaK</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_Na_NaK</ci>
                <apply>
                    <times/>
                    <ci>J_Na_NaK</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NaK</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>TO_NaK</ci>
                                        <ci>n_NaK</ci>
                                    </apply>
                                    <cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
                                </apply>
                                <apply>
                                    <power/>
                                    <ci>Na_i</ci>
                                    <cn cellml:units="dimensionless">1.36</cn>
                                </apply>
                                <cn cellml:units="dimensionless">1.13</cn>
                            </apply>
                            <apply>
                                <plus/>
                                <apply>
                                    <power/>
                                    <ci>K_Na_i_half</ci>
                                    <cn cellml:units="dimensionless">1.36</cn>
                                </apply>
                                <apply>
                                    <power/>
                                    <ci>Na_i</ci>
                                    <cn cellml:units="dimensionless">1.36</cn>
                                </apply>
                            </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1.3</cn>
                        <ci>K_s</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>K_K_s_half</ci>
                        <ci>K_s</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>TO_NaK</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <apply>
                                <times/>
                                <cn cellml:units="per_volt4_per_second" type="e-notation">5.46<sep/>+4</cn>
                                <apply>
                                    <power/>
                                    <ci>v_sc</ci>
                                    <cn cellml:units="dimensionless">4</cn>
                                </apply>
                            </apply>
                            <apply>
                                <times/>
                                <cn cellml:units="per_volt3_per_second" type="e-notation">2.43<sep/>+2</cn>
                                <apply>
                                    <power/>
                                    <ci>v_sc</ci>
                                    <cn cellml:units="dimensionless">3</cn>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="per_volt2_per_second" type="e-notation">2<sep/>+3</cn>
                            <apply>
                                <power/>
                                <ci>v_sc</ci>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <cn cellml:units="per_volt_per_second">160</cn>
                        <ci>v_sc</ci>
                    </apply>
                    <cn cellml:units="per_second">51</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_Na_NaK</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">3</cn>
                    <ci>J_NaK</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_K_NaK</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">2</cn>
                    <ci>J_NaK</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="NBC">
        <variable name="time" public_interface="in" units="second"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="v_cell" public_interface="in" units="m3"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="Na_s" public_interface="in" units="M"/>
        <variable name="HCO3_i" public_interface="in" units="M"/>
        <variable name="HCO3_s" public_interface="in" units="M"/>
        <variable initial_value="0.08" name="g_NBC" public_interface="out" units="S_per_m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <!--kelvin {init: 300};-->
        <variable initial_value="2" name="n" units="dimensionless"/>
        <!--var eps: M {init: 1e-12};-->
        <variable name="E_NBC" public_interface="out" units="volt"/>
        <variable name="J_NBC_Na" public_interface="out" units="umol_per_s"/>
        <variable name="J_NBC_HCO3" public_interface="out" units="umol_per_s"/>
        <variable name="I_NBC_Na" public_interface="out" units="uA"/>
        <variable name="I_NBC_HCO3" public_interface="out" units="uA"/>
        <variable initial_value="1e-18" name="eps" units="M"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>E_NBC</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>n</ci>
                            </apply>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <ci>Na_s</ci>
                                <ci>Na_i</ci>
                            </apply>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <ci>HCO3_s</ci>
                                    <ci>HCO3_i</ci>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NBC_Na</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                        <ci>g_NBC</ci>
                        <ci>A_Apical</ci>
                        <apply>
                            <minus/>
                            <ci>v_sc</ci>
                            <ci>E_NBC</ci>
                        </apply>
                    </apply>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NBC_HCO3</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">3.0</cn>
                    <ci>J_NBC_Na</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NBC_Na</ci>
                <apply>
                    <times/>
                    <ci>J_NBC_Na</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NBC_HCO3</ci>
                <apply>
                    <times/>
                    <ci>J_NBC_HCO3</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="NHE3">
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="v_cell" public_interface="in" units="m3"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="H_Ext" public_interface="in" units="M"/>
        <variable name="H_int" public_interface="in" units="M"/>
        <variable initial_value="0" name="NH4_ext" public_interface="out" units="M"/>
        <variable initial_value="0" name="NH4_int" public_interface="out" units="M"/>
        <variable name="J_NHE3_Na" public_interface="out" units="umol_per_s"/>
        <variable name="J_NHE3_H1" public_interface="out" units="umol_per_s"/>
        <variable name="J_NHE3_H" public_interface="out" units="umol_per_s"/>
        <variable name="J_NHE3_NH4" public_interface="out" units="umol_per_s"/>
        <variable name="J_NHE3_Na_Max" public_interface="out" units="umol_per_s"/>
        <variable name="I_NHE3_H" public_interface="out" units="uA"/>
        <variable name="I_NHE3_Na" public_interface="out" units="uA"/>
        <variable initial_value="0.001" name="x_T" units="M"/>
        <variable name="sigma" units="per_second"/>
        <variable name="P_Na" units="per_second"/>
        <variable name="P_H" units="per_second"/>
        <variable name="P_NH4" units="per_second"/>
        <variable initial_value="1.6" name="P0_Na" units="per_second"/>
        <variable initial_value="0.48" name="P0_H" units="per_second"/>
        <variable initial_value="1.6" name="P0_NH4" units="per_second"/>
        <variable initial_value="0.03" name="K_Na" units="M"/>
        <variable initial_value="7.2e-8" name="K_H" units="M"/>
        <variable initial_value="0.027" name="K_NH4" units="M"/>
        <variable initial_value="1e-9" name="K_I" units="M"/>
        <variable initial_value="0" name="f_m" units="dimensionless"/>
        <variable initial_value="2" name="f_M" units="dimensionless"/>
        <variable name="alpha_ext_Na" units="dimensionless"/>
        <variable name="alpha_int_Na" units="dimensionless"/>
        <variable name="beta_ext_H" units="dimensionless"/>
        <variable name="beta_int_H" units="dimensionless"/>
        <variable name="gamma_ext_NH4" units="dimensionless"/>
        <variable name="gamma_int_NH4" units="dimensionless"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <!--        var pH_Ext: dimensionless;
        var pH_int: dimensionless;
        var p_K: dimensionless {init: 6.1, pub: out};
        var P_CO2: dimensionless {init: 40, pub: out};
        var s: dimensionless {init: 0.03e-3, pub: out};
        var HCO3_i: M {pub: in};
        var HCO3_m: M {pub: in};-->
        <!--        pH_Ext = p_K+log(HCO3_m/(s*P_CO2));
        pH_int = p_K+log(HCO3_i/(s*P_CO2));-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>alpha_ext_Na</ci>
                <apply>
                    <divide/>
                    <ci>Na_m</ci>
                    <ci>K_Na</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>alpha_int_Na</ci>
                <apply>
                    <divide/>
                    <ci>Na_i</ci>
                    <ci>K_Na</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>beta_ext_H</ci>
                <apply>
                    <divide/>
                    <ci>H_Ext</ci>
                    <ci>K_H</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>beta_int_H</ci>
                <apply>
                    <divide/>
                    <ci>H_int</ci>
                    <ci>K_H</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gamma_ext_NH4</ci>
                <apply>
                    <divide/>
                    <ci>NH4_ext</ci>
                    <ci>K_NH4</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gamma_int_NH4</ci>
                <apply>
                    <divide/>
                    <ci>NH4_int</ci>
                    <ci>K_NH4</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>P_Na</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>P0_Na</ci>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>f_M</ci>
                                <ci>H_int</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>f_m</ci>
                                <ci>K_I</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>H_int</ci>
                        <ci>K_I</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>P_H</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>P0_H</ci>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>f_M</ci>
                                <ci>H_int</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>f_m</ci>
                                <ci>K_I</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>H_int</ci>
                        <ci>K_I</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>P_NH4</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>P0_NH4</ci>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>f_M</ci>
                                <ci>H_int</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>f_m</ci>
                                <ci>K_I</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>H_int</ci>
                        <ci>K_I</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sigma</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>alpha_ext_Na</ci>
                            <ci>beta_ext_H</ci>
                            <ci>gamma_ext_NH4</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_Na</ci>
                                <ci>alpha_int_Na</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_H</ci>
                                <ci>beta_int_H</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_NH4</ci>
                                <ci>gamma_int_NH4</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>alpha_int_Na</ci>
                            <ci>beta_int_H</ci>
                            <ci>gamma_int_NH4</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_Na</ci>
                                <ci>alpha_ext_Na</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_H</ci>
                                <ci>beta_ext_H</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_NH4</ci>
                                <ci>gamma_ext_NH4</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NHE3_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
                            </apply>
                            <ci>v_cell</ci>
                            <ci>x_T</ci>
                        </apply>
                        <ci>sigma</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>P_Na</ci>
                            <ci>P_H</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>alpha_int_Na</ci>
                                    <ci>beta_ext_H</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>alpha_ext_Na</ci>
                                    <ci>beta_int_H</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_Na</ci>
                            <ci>P_NH4</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>alpha_int_Na</ci>
                                    <ci>gamma_ext_NH4</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>alpha_ext_Na</ci>
                                    <ci>gamma_int_NH4</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NHE3_H</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
                            </apply>
                            <ci>v_cell</ci>
                            <ci>x_T</ci>
                        </apply>
                        <ci>sigma</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>P_Na</ci>
                            <ci>P_H</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>alpha_ext_Na</ci>
                                    <ci>beta_int_H</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>alpha_int_Na</ci>
                                    <ci>beta_ext_H</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_H</ci>
                            <ci>P_NH4</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>beta_int_H</ci>
                                    <ci>gamma_ext_NH4</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>beta_ext_H</ci>
                                    <ci>gamma_int_NH4</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NHE3_H1</ci>
                <apply>
                    <times/>
                    <ci>J_NHE3_H</ci>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NHE3_NH4</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <ci>x_T</ci>
                        <ci>sigma</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>P_Na</ci>
                            <ci>P_NH4</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>alpha_ext_Na</ci>
                                    <ci>gamma_int_NH4</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>alpha_int_Na</ci>
                                    <ci>gamma_ext_NH4</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_H</ci>
                            <ci>P_NH4</ci>
                            <apply>
                                <minus/>
                                <apply>
                                    <times/>
                                    <ci>beta_ext_H</ci>
                                    <ci>gamma_int_NH4</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>gamma_ext_NH4</ci>
                                    <ci>beta_int_H</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NHE3_Na_Max</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>x_T</ci>
                        <ci>P_Na</ci>
                        <ci>P_H</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>P_Na</ci>
                        <ci>P_H</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NHE3_H</ci>
                <apply>
                    <times/>
                    <ci>J_NHE3_H</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NHE3_Na</ci>
                <apply>
                    <times/>
                    <ci>J_NHE3_Na</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="AE1">
        <!--var A_Apical: m2 {pub: in};-->
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="v_cell" public_interface="in" units="m3"/>
        <variable name="HCO3_i" public_interface="in" units="M"/>
        <variable name="HCO3_m" public_interface="in" units="M"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="Cl_m" public_interface="in" units="M"/>
        <variable name="J_AE1_HCO3" public_interface="out" units="umol_per_s"/>
        <!--var J_AE1_HCO31: umol_per_s {pub: out};-->
        <variable name="J_AE1_Cl" public_interface="out" units="umol_per_s"/>
        <!--var J_AE1_Cl1: umol_per_s {pub: out};-->
        <variable name="I_AE1_HCO3" public_interface="out" units="uA"/>
        <variable name="I_AE1_Cl" public_interface="out" units="uA"/>
        <variable name="J_HCO3_influx" public_interface="out" units="umol_per_s"/>
        <variable name="J_Cl_influx" public_interface="out" units="umol_per_s"/>
        <variable initial_value="0.198" name="K_HCO3_m" units="M"/>
        <variable initial_value="0.198" name="K_HCO3_i" units="M"/>
        <variable initial_value="0.05" name="K_Cl_s" units="M"/>
        <variable initial_value="0.05" name="K_Cl_i" units="M"/>
        <variable initial_value="1247" name="P_HCO3_m" units="per_second"/>
        <variable initial_value="135" name="P_HCO3_i" units="per_second"/>
        <variable initial_value="562" name="P_Cl_s" units="per_second"/>
        <variable initial_value="61" name="P_Cl_i" units="per_second"/>
        <variable name="beta_ext" units="dimensionless"/>
        <variable name="beta_int" units="dimensionless"/>
        <variable name="gamma_ext" units="dimensionless"/>
        <variable name="gamma_int" units="dimensionless"/>
        <variable name="sigma" units="per_second"/>
        <variable initial_value="0.001" name="x_Tmax" units="M"/>
        <variable initial_value="0.172" name="K_I" units="M"/>
        <variable name="x_T" units="M"/>
        <variable name="x_ext" units="M"/>
        <variable name="x_int" units="M"/>
        <variable name="Jo_bm" units="umol_per_s"/>
        <variable name="Ji_bm" units="umol_per_s"/>
        <variable name="Js_bm" units="umol_per_s"/>
        <variable name="Jo_cm" units="umol_per_s"/>
        <variable name="Ji_cm" units="umol_per_s"/>
        <variable name="Js_cm" units="umol_per_s"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>x_T</ci>
                <apply>
                    <divide/>
                    <ci>x_Tmax</ci>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <divide/>
                            <ci>HCO3_i</ci>
                            <ci>K_I</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>beta_ext</ci>
                <apply>
                    <divide/>
                    <ci>HCO3_m</ci>
                    <ci>K_HCO3_m</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>beta_int</ci>
                <apply>
                    <divide/>
                    <ci>HCO3_i</ci>
                    <ci>K_HCO3_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gamma_ext</ci>
                <apply>
                    <divide/>
                    <ci>Cl_m</ci>
                    <ci>K_Cl_s</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gamma_int</ci>
                <apply>
                    <divide/>
                    <ci>Cl_i</ci>
                    <ci>K_Cl_i</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>sigma</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>beta_ext</ci>
                            <ci>gamma_ext</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_HCO3_i</ci>
                                <ci>beta_int</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_Cl_i</ci>
                                <ci>gamma_int</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>beta_int</ci>
                            <ci>gamma_int</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_HCO3_m</ci>
                                <ci>beta_ext</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_Cl_s</ci>
                                <ci>gamma_ext</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_AE1_HCO3</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
                            <ci>v_cell</ci>
                            <ci>x_Tmax</ci>
                        </apply>
                        <ci>sigma</ci>
                    </apply>
                    <apply>
                        <minus/>
                        <apply>
                            <times/>
                            <ci>P_HCO3_i</ci>
                            <ci>beta_int</ci>
                            <ci>P_Cl_s</ci>
                            <ci>gamma_ext</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_HCO3_m</ci>
                            <ci>beta_ext</ci>
                            <ci>P_Cl_i</ci>
                            <ci>gamma_int</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--J_AE1_HCO31 = J_AE1_HCO3*36e8{dimensionless};
J_AE1_Cl1 = -J_AE1_HCO31;-->
            <apply>
                <eq/>
                <ci>J_AE1_Cl</ci>
                <ci>J_AE1_HCO3</ci>
            </apply>
            <apply>
                <eq/>
                <ci>I_AE1_HCO3</ci>
                <apply>
                    <times/>
                    <ci>J_AE1_HCO3</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_AE1_Cl</ci>
                <apply>
                    <times/>
                    <ci>J_AE1_Cl</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>x_ext</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>x_T</ci>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_HCO3_i</ci>
                                <ci>beta_int</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_Cl_i</ci>
                                <ci>gamma_int</ci>
                            </apply>
                        </apply>
                    </apply>
                    <ci>sigma</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>x_int</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>x_T</ci>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>P_HCO3_m</ci>
                                <ci>beta_ext</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>P_Cl_s</ci>
                                <ci>gamma_ext</ci>
                            </apply>
                        </apply>
                    </apply>
                    <ci>sigma</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_HCO3_influx</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <ci>x_T</ci>
                        <ci>sigma</ci>
                    </apply>
                    <ci>P_HCO3_m</ci>
                    <ci>beta_ext</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>P_HCO3_i</ci>
                            <ci>beta_int</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_Cl_i</ci>
                            <ci>gamma_int</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_Cl_influx</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <ci>x_T</ci>
                        <ci>sigma</ci>
                    </apply>
                    <ci>P_Cl_s</ci>
                    <ci>gamma_ext</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>P_HCO3_i</ci>
                            <ci>beta_int</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>P_Cl_i</ci>
                            <ci>gamma_int</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Jo_bm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_m</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_i</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <ci>K_HCO3_i</ci>
                                <apply>
                                    <times/>
                                    <ci>P_HCO3_i</ci>
                                    <ci>HCO3_i</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ji_bm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_m</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_i</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <ci>K_HCO3_m</ci>
                                <apply>
                                    <times/>
                                    <ci>P_HCO3_m</ci>
                                    <ci>HCO3_m</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Js_bm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_m</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_HCO3_i</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Jo_cm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_s</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_i</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <ci>K_Cl_i</ci>
                                <apply>
                                    <times/>
                                    <ci>P_Cl_i</ci>
                                    <ci>Cl_i</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Ji_cm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_s</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_i</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <ci>K_Cl_s</ci>
                                <apply>
                                    <times/>
                                    <ci>P_Cl_s</ci>
                                    <ci>Cl_m</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Js_cm</ci>
                <apply>
                    <power/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>x_T</ci>
                        </apply>
                        <apply>
                            <plus/>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_s</ci>
                            </apply>
                            <apply>
                                <divide/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>P_Cl_i</ci>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <!--    def comp Hydration as
        var time: second {pub: in};
        var v_cell: m3 {pub: in};
        var CO2_i: M {init: 1e-3, pub: out};
        var CO2_s: M {init: 1.2e-3, pub: out};
        var CO2_m: M {init: 1.2e-3, pub: out};
        var H_int: M {pub: in};
        var H_Ext: M {pub: in};
        var HCO3_i: M {pub: in};
        var HCO3_m: M {pub: in};
        var p_CO2: per_second {init: 5e-6, pub: out};
        var k_f: per_second {init: 0.03, pub: out};
        var k_r: per_M_per_second {init: 20e-3, pub: out};
        var J_CDF_A: umol_per_s {pub: out};
        var J_CDF_B: umol_per_s {pub: out};
        var J_buf_C: umol_per_s {pub: out};
        var J_buf_C1: umol_per_s {pub: out};
        var J_buf_A: umol_per_s {pub: out};
        var J_D_P_CO2: umol_per_s {pub: in};-->
    <!--        J_CDF_A = -1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_m);
        J_CDF_B = 1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_s);
        J_buf_C = 1e9{dimensionless}*v_cell*(k_f*CO2_i-k_r*HCO3_i*H_int);
        J_buf_C1 = J_buf_C*36e8{dimensionless};
        J_buf_A = 1e9{dimensionless}*v_cell*(k_f*CO2_m-k_r*HCO3_m*H_Ext);
        ode(CO2_i, time) = -J_buf_C+J_CDF_A+J_CDF_B;
        ode(CO2_m, time) = -J_buf_A-J_CDF_A-0{dimensionless}*J_D_P_CO2;
    enddef;-->
    <component name="IK">
        <variable name="z_k" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="K_i" public_interface="in" units="M"/>
        <variable name="K_s" public_interface="in" units="M"/>
        <variable name="Ca_i" public_interface="in" units="M"/>
        <variable name="V_A_IK" public_interface="out" units="volt"/>
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="J_IK" public_interface="out" units="umol_per_s"/>
        <variable name="I_IK" public_interface="out" units="uA"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="1.7" name="G_IK" public_interface="out" units="S_per_m2"/>
        <!--var P_IK: dimensionless {pub: out};-->
        <variable initial_value="0.0003" name="K_s_IK" public_interface="out" units="M"/>
        <variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
        <!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
        <!--var alpha_IK: per_m2 {init: 1, pub: out};-->
        <!--P_IK = 1{dimensionless}/(1{dimensionless}+pow(K_s_IK/Ca_i, eta));-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>V_A_IK</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_k</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>K_s</ci>
                            <ci>K_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_IK</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
                    </apply>
                    <ci>G_IK</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_sc</ci>
                        <ci>V_A_IK</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_IK</ci>
                <apply>
                    <divide/>
                    <ci>I_IK</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="BK">
        <variable name="z_k" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="K_i" public_interface="in" units="M"/>
        <variable name="K_m" public_interface="in" units="M"/>
        <variable name="Ca_i" public_interface="in" units="M"/>
        <variable name="V_A_BK" public_interface="out" units="volt"/>
        <variable name="v_mc" public_interface="in" units="volt"/>
        <variable name="J_BK" public_interface="out" units="umol_per_s"/>
        <variable name="I_BK" public_interface="out" units="uA"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable initial_value="1" name="G_BK" public_interface="out" units="S_per_m2"/>
        <variable name="P_BK" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.0003" name="K_s_BK" public_interface="out" units="M"/>
        <variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
        <!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
        <!--var alpha_BK: per_m2 {init: 0.1, pub: out};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>P_BK</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>K_s_BK</ci>
                                <ci>Ca_i</ci>
                            </apply>
                            <ci>eta</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_A_BK</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_k</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>K_m</ci>
                            <ci>K_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_BK</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
                    </apply>
                    <ci>G_BK</ci>
                    <ci>P_BK</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_mc</ci>
                        <ci>V_A_BK</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_BK</ci>
                <apply>
                    <divide/>
                    <ci>I_BK</ci>
                    <ci>F</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Apical_Diffusion">
        <variable name="v_mc" public_interface="in" units="volt"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="K_m" public_interface="in" units="M"/>
        <variable name="Cl_m" public_interface="in" units="M"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="K_i" public_interface="in" units="M"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="z_Na" public_interface="in" units="dimensionless"/>
        <variable name="z_k" public_interface="in" units="dimensionless"/>
        <variable name="z_Cl" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="P_Na" public_interface="in" units="m_per_s"/>
        <variable name="P_k" public_interface="in" units="m_per_s"/>
        <variable name="P_Cl" public_interface="in" units="m_per_s"/>
        <variable name="J_D_Na" public_interface="out" units="umol_per_s"/>
        <variable name="J_D_NA" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_Na" public_interface="out" units="uA"/>
        <variable name="J_D_K" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_K" public_interface="out" units="uA"/>
        <variable name="J_D_Cl" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_Cl" public_interface="out" units="uA"/>
        <variable name="J_D_Apical" public_interface="out" units="umol_per_s"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_D_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_Na</ci>
                                    <ci>F</ci>
                                    <ci>z_Na</ci>
                                    <ci>v_mc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>Na_m</ci>
                                <apply>
                                    <times/>
                                    <ci>Na_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_Na</ci>
                                                <ci>F</ci>
                                                <ci>v_mc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_Na</ci>
                                        <ci>F</ci>
                                        <ci>v_mc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Apical</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_NA</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
                    <ci>J_D_Na</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_Na</ci>
                <apply>
                    <times/>
                    <ci>J_D_Na</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_K</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_k</ci>
                                    <ci>F</ci>
                                    <ci>z_k</ci>
                                    <ci>v_mc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>K_m</ci>
                                <apply>
                                    <times/>
                                    <ci>K_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_k</ci>
                                                <ci>F</ci>
                                                <ci>v_mc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_k</ci>
                                        <ci>F</ci>
                                        <ci>v_mc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Apical</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_K</ci>
                <apply>
                    <times/>
                    <ci>J_D_K</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_Cl</ci>
                                    <ci>F</ci>
                                    <ci>z_Cl</ci>
                                    <ci>v_mc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>Cl_m</ci>
                                <apply>
                                    <times/>
                                    <ci>Cl_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_Cl</ci>
                                                <ci>F</ci>
                                                <ci>v_mc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_Cl</ci>
                                        <ci>F</ci>
                                        <ci>v_mc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Apical</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_Cl</ci>
                <apply>
                    <times/>
                    <ci>J_D_Cl</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Apical</ci>
                <apply>
                    <plus/>
                    <ci>J_D_Na</ci>
                    <ci>J_D_K</ci>
                    <ci>J_D_Cl</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Basolateral_Diffusion">
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="Na_s" public_interface="in" units="M"/>
        <variable name="K_s" public_interface="in" units="M"/>
        <variable name="Cl_s" public_interface="in" units="M"/>
        <variable name="Na_i" public_interface="in" units="M"/>
        <variable name="K_i" public_interface="in" units="M"/>
        <variable name="Cl_i" public_interface="in" units="M"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="z_Na" public_interface="in" units="dimensionless"/>
        <variable name="z_k" public_interface="in" units="dimensionless"/>
        <variable name="z_Cl" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="A_Basol" public_interface="in" units="m2"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="P_Bl_Na" public_interface="in" units="m_per_s"/>
        <variable name="P_Bl_k" public_interface="in" units="m_per_s"/>
        <variable name="P_Bl_Cl" public_interface="in" units="m_per_s"/>
        <variable name="J_D_Bl_Na" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_Bl_Na" public_interface="out" units="uA"/>
        <variable name="J_D_Bl_K" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_Bl_K" public_interface="out" units="uA"/>
        <variable name="J_D_Bl_Cl" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_Bl_Cl" public_interface="out" units="uA"/>
        <variable name="J_D_Basolateral" public_interface="out" units="umol_per_s"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_D_Bl_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_Bl_Na</ci>
                                    <ci>F</ci>
                                    <ci>z_Na</ci>
                                    <ci>v_sc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>Na_s</ci>
                                <apply>
                                    <times/>
                                    <ci>Na_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_Na</ci>
                                                <ci>F</ci>
                                                <ci>v_sc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_Na</ci>
                                        <ci>F</ci>
                                        <ci>v_sc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Basol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_Bl_Na</ci>
                <apply>
                    <times/>
                    <ci>J_D_Bl_Na</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Bl_K</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_Bl_k</ci>
                                    <ci>F</ci>
                                    <ci>z_k</ci>
                                    <ci>v_sc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>K_s</ci>
                                <apply>
                                    <times/>
                                    <ci>K_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_k</ci>
                                                <ci>F</ci>
                                                <ci>v_sc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_k</ci>
                                        <ci>F</ci>
                                        <ci>v_sc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Basol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_Bl_K</ci>
                <apply>
                    <times/>
                    <ci>J_D_Bl_K</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Bl_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
                                    </apply>
                                    <ci>P_Bl_Cl</ci>
                                    <ci>F</ci>
                                    <ci>z_Cl</ci>
                                    <ci>v_sc</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>T</ci>
                                </apply>
                            </apply>
                            <apply>
                                <minus/>
                                <ci>Cl_s</ci>
                                <apply>
                                    <times/>
                                    <ci>Cl_i</ci>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <times/>
                                                <ci>z_Cl</ci>
                                                <ci>F</ci>
                                                <ci>v_sc</ci>
                                            </apply>
                                            <apply>
                                                <times/>
                                                <ci>R</ci>
                                                <ci>T</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <ci>z_Cl</ci>
                                        <ci>F</ci>
                                        <ci>v_sc</ci>
                                    </apply>
                                    <apply>
                                        <times/>
                                        <ci>R</ci>
                                        <ci>T</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <ci>A_Basol</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_Bl_Cl</ci>
                <apply>
                    <times/>
                    <ci>J_D_Bl_Cl</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Basolateral</ci>
                <apply>
                    <plus/>
                    <ci>J_D_Bl_Na</ci>
                    <ci>J_D_Bl_K</ci>
                    <ci>J_D_Bl_Cl</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Paracellular_voltage">
        <variable name="v_ms" public_interface="out" units="volt"/>
        <variable name="v_sc" public_interface="in" units="volt"/>
        <variable name="v_mc" public_interface="in" units="volt"/>
        <!--        var time: second {pub: in};
        var capacitance: uF {pub: in};
        var I_D_P_Na: uA {pub: in};
        var I_D_P_K: uA {pub: in};
        var I_D_P_Cl: uA {pub: in};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_ms</ci>
                <apply>
                    <minus/>
                    <ci>v_mc</ci>
                    <ci>v_sc</ci>
                </apply>
            </apply>
        </math>
        <!--ode(v_ms, time) = (I_D_P_Na+I_D_P_K+I_D_P_Cl)/capacitance;-->
    </component>
    <component name="Paracellular_Diffusion">
        <variable name="v_ms" public_interface="in" units="volt"/>
        <variable name="Na_m" public_interface="in" units="M"/>
        <variable name="K_m" public_interface="in" units="M"/>
        <variable name="Cl_m" public_interface="in" units="M"/>
        <variable name="Na_s" public_interface="in" units="M"/>
        <variable name="K_s" public_interface="in" units="M"/>
        <variable name="Cl_s" public_interface="in" units="M"/>
        <variable name="F" public_interface="in" units="C_per_mol"/>
        <variable name="z_Na" public_interface="in" units="dimensionless"/>
        <variable name="z_k" public_interface="in" units="dimensionless"/>
        <variable name="z_Cl" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="J_per_K_per_mol"/>
        <variable name="A_Apical" public_interface="in" units="m2"/>
        <variable name="T" public_interface="in" units="kelvin"/>
        <variable name="P_P_Na" public_interface="in" units="m_per_s"/>
        <variable name="P_P_k" public_interface="in" units="m_per_s"/>
        <variable name="P_P_Cl" public_interface="in" units="m_per_s"/>
        <variable name="J_D_P_Na" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_P_Na" public_interface="out" units="uA"/>
        <variable name="J_D_P_K" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_P_K" public_interface="out" units="uA"/>
        <variable name="J_D_P_Cl" public_interface="out" units="umol_per_s"/>
        <variable name="I_D_P_Cl" public_interface="out" units="uA"/>
        <variable name="J_D_Paracellular" public_interface="out" units="umol_per_s"/>
        <variable name="V_P_Na" public_interface="out" units="volt"/>
        <variable name="V_P_K" public_interface="out" units="volt"/>
        <variable name="V_P_Cl" public_interface="out" units="volt"/>
        <variable initial_value="5" name="G_P_Na" public_interface="out" units="S_per_m2"/>
        <variable initial_value="3" name="G_P_K" public_interface="out" units="S_per_m2"/>
        <variable initial_value="5" name="G_P_Cl" public_interface="out" units="S_per_m2"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>V_P_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_Na</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Na_m</ci>
                            <ci>Na_s</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_P_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <ci>G_P_Na</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_ms</ci>
                        <ci>V_P_Na</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_P_Na</ci>
                <apply>
                    <divide/>
                    <ci>I_D_P_Na</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_P_K</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_k</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>K_m</ci>
                            <ci>K_s</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_P_K</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <ci>G_P_K</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_ms</ci>
                        <ci>V_P_K</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_P_K</ci>
                <apply>
                    <divide/>
                    <ci>I_D_P_K</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>V_P_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>z_Cl</ci>
                            <ci>F</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Cl_m</ci>
                            <ci>Cl_s</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_D_P_Cl</ci>
                <apply>
                    <times/>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
                    </apply>
                    <ci>G_P_Cl</ci>
                    <ci>A_Apical</ci>
                    <apply>
                        <minus/>
                        <ci>v_ms</ci>
                        <ci>V_P_Cl</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_P_Cl</ci>
                <apply>
                    <divide/>
                    <ci>I_D_P_Cl</ci>
                    <ci>F</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_D_Paracellular</ci>
                <apply>
                    <plus/>
                    <ci>J_D_P_Na</ci>
                    <ci>J_D_P_K</ci>
                    <ci>J_D_P_Cl</ci>
                </apply>
            </apply>
        </math>
    </component>
    <!--    def comp interstitium_concentration as
        var glucose_s: M {init: 0.005, pub: out};-->
    <!--        //var J_GLUT: umol_per_s {pub: in};
        //var J_G_Blood: umol_per_s {pub: in};
        //var v_cell: m3 {pub: in};
        //var time: second {pub: in};-->
    <!--        //ode(glucose_s, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*(-J_GLUT-J_G_Blood);
    enddef;-->
    <!--    def comp Blood_concentration as
        var glucose_s: M {pub: in};
        var glucose_B: M {init: 0.012, pub: out};
        var J_G_Blood: umol_per_s {pub: out};
        var K: dimensionless {init: 1e-6, pub: out};-->
    <!--        //var v_cell: m3 {pub: in};-->
    <!--        var time: second {pub: in};-->
    <!--        J_G_Blood = K*(glucose_s-glucose_B);
        ode(glucose_B, time) = 0{dimensionless};
    enddef;-->
    <connection>
        <map_components component_1="rate_constants" component_2="phenomonological_constants"/>
        <map_variables variable_1="ks_12" variable_2="ks_12"/>
        <!--vars k_12 and k_12;-->
        <map_variables variable_1="k_21" variable_2="k_21"/>
        <!--vars k_23 and k_23;-->
        <map_variables variable_1="k_32" variable_2="k_32"/>
        <map_variables variable_1="k_25" variable_2="k_25"/>
        <map_variables variable_1="k_52" variable_2="k_52"/>
        <map_variables variable_1="k_34" variable_2="k_34"/>
        <map_variables variable_1="k_43" variable_2="k_43"/>
        <map_variables variable_1="k_45" variable_2="k_45"/>
        <map_variables variable_1="k_54" variable_2="k_54"/>
        <map_variables variable_1="k_56" variable_2="k_56"/>
        <map_variables variable_1="k_65" variable_2="k_65"/>
        <map_variables variable_1="k_61" variable_2="k_61"/>
        <map_variables variable_1="k_16" variable_2="k_16"/>
    </connection>
    <connection>
        <map_components component_1="phenomonological_constants" component_2="Apical_concentrations"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
        <map_variables variable_1="glucose_m" variable_2="glucose_m"/>
    </connection>
    <connection>
        <map_components component_1="phenomonological_constants" component_2="Cell_concentration"/>
        <map_variables variable_1="J_Na_SGLT" variable_2="J_Na_SGLT"/>
        <map_variables variable_1="J_Gl_SGLT" variable_2="J_Gl_SGLT"/>
    </connection>
    <connection>
        <map_components component_1="phenomonological_constants" component_2="parameters"/>
        <map_variables variable_1="k0_23" variable_2="k0_23"/>
        <!--vars C_T and C_T;-->
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="phenomonological_constants" component_2="Apical_voltage"/>
        <map_variables variable_1="I_Na_SGLT" variable_2="I_Na_SGLT"/>
        <!--vars I_Gl_SGLT and I_Gl_SGLT;-->
    </connection>
    <connection>
        <map_components component_1="rate_constants" component_2="Apical_concentrations"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
        <map_variables variable_1="glucose_m" variable_2="glucose_m"/>
    </connection>
    <connection>
        <map_components component_1="rate_constants" component_2="Cell_concentration"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="glucose_i" variable_2="glucose_i"/>
    </connection>
    <connection>
        <map_components component_1="rate_constants" component_2="parameters"/>
        <map_variables variable_1="k0_12" variable_2="k0_12"/>
        <map_variables variable_1="k0_21" variable_2="k0_21"/>
        <!--vars k0_23 and k0_23;-->
        <map_variables variable_1="k0_32" variable_2="k0_32"/>
        <map_variables variable_1="k0_25" variable_2="k0_25"/>
        <map_variables variable_1="k0_52" variable_2="k0_52"/>
        <map_variables variable_1="k0_34" variable_2="k0_34"/>
        <map_variables variable_1="k0_43" variable_2="k0_43"/>
        <map_variables variable_1="k0_45" variable_2="k0_45"/>
        <map_variables variable_1="k0_54" variable_2="k0_54"/>
        <map_variables variable_1="k0_56" variable_2="k0_56"/>
        <map_variables variable_1="k0_65" variable_2="k0_65"/>
        <map_variables variable_1="k0_61" variable_2="k0_61"/>
        <map_variables variable_1="k0_16" variable_2="k0_16"/>
        <map_variables variable_1="delta" variable_2="delta"/>
        <map_variables variable_1="alpha_p" variable_2="alpha_p"/>
        <!--vars alpha_pp and alpha_pp;-->
        <!--vars C_T and C_T;-->
        <map_variables variable_1="n" variable_2="n"/>
        <map_variables variable_1="z_c" variable_2="z_c"/>
        <map_variables variable_1="z_Na" variable_2="z_Na"/>
        <map_variables variable_1="mu" variable_2="mu"/>
    </connection>
    <connection>
        <map_components component_1="CFTR" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Cl" variable_2="z_Cl"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="CLC2" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Cl" variable_2="z_Cl"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="CFTR" component_2="Cell_concentration"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
        <map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
        <map_variables variable_1="J_CFTR" variable_2="J_CFTR"/>
        <map_variables variable_1="J_CFTR_HCO3" variable_2="J_CFTR_HCO3"/>
    </connection>
    <connection>
        <map_components component_1="CLC2" component_2="Cell_concentration"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <!--        vars HCO3_i and HCO3_i;
        vars HCO3_s and HCO3_s;-->
        <map_variables variable_1="J_CLC2" variable_2="J_CLC2"/>
        <!--vars J_CFTR_HCO3 and J_CFTR_HCO3;-->
    </connection>
    <connection>
        <map_components component_1="CFTR" component_2="Apical_concentrations"/>
        <map_variables variable_1="Cl_m" variable_2="Cl_m"/>
    </connection>
    <connection>
        <map_components component_1="CLC2" component_2="Basolateral_concentrations"/>
        <map_variables variable_1="Cl_s" variable_2="Cl_s"/>
        <map_variables variable_1="J_CLC2" variable_2="J_CLC2"/>
    </connection>
    <connection>
        <map_components component_1="CFTR" component_2="Apical_voltage"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
        <map_variables variable_1="I_CFTR" variable_2="I_CFTR"/>
        <!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
    </connection>
    <connection>
        <map_components component_1="CLC2" component_2="Basol_voltage"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
        <map_variables variable_1="I_CLC2" variable_2="I_CLC2"/>
        <!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
    </connection>
    <connection>
        <map_components component_1="ENaC" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Na" variable_2="z_Na"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="ENaC" component_2="Cell_concentration"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="J_ENaC" variable_2="J_ENaC"/>
    </connection>
    <connection>
        <map_components component_1="ENaC" component_2="Apical_concentrations"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
    </connection>
    <connection>
        <map_components component_1="ENaC" component_2="Apical_voltage"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
        <map_variables variable_1="I_ENaC" variable_2="I_ENaC"/>
    </connection>
    <connection>
        <map_components component_1="IK" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_k" variable_2="z_k"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="IK" component_2="Cell_concentration"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="J_IK" variable_2="J_IK"/>
    </connection>
    <connection>
        <map_components component_1="IK" component_2="Basolateral_concentrations"/>
        <map_variables variable_1="K_s" variable_2="K_s"/>
        <map_variables variable_1="J_IK" variable_2="J_IK"/>
    </connection>
    <connection>
        <map_components component_1="IK" component_2="Basol_voltage"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
        <map_variables variable_1="I_IK" variable_2="I_IK"/>
    </connection>
    <connection>
        <map_components component_1="BK" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_k" variable_2="z_k"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="BK" component_2="Cell_concentration"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="J_BK" variable_2="J_BK"/>
    </connection>
    <connection>
        <map_components component_1="BK" component_2="Apical_concentrations"/>
        <map_variables variable_1="K_m" variable_2="K_m"/>
    </connection>
    <connection>
        <map_components component_1="BK" component_2="Apical_voltage"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
        <map_variables variable_1="I_BK" variable_2="I_BK"/>
    </connection>
    <connection>
        <map_components component_1="NaKCC" component_2="Cell_concentration"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
    </connection>
    <connection>
        <map_components component_1="NaKCC" component_2="parameters"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="GLUT2"/>
        <map_variables variable_1="glucose_i" variable_2="glucose_i"/>
        <map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="GLUT2"/>
        <map_variables variable_1="glucose_s" variable_2="glucose_s"/>
        <map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
        <!--vars J_GLUT and J_GLUT;-->
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="parameters"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
    </connection>
    <!--    def map between interstitium_concentration and parameters for
        vars time and time;
        vars v_cell and v_cell;
    enddef;-->
    <connection>
        <map_components component_1="parameters" component_2="GLUT2"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
    </connection>
    <!--    def map between GLUT2 and Basol_voltage for
        vars I_GLUT and I_GLUT;
    enddef;-->
    <!--    def map between Cell_concentration and Na_Cl for
        vars Na_i and Na_i;
        vars Cl_i and Cl_i;
        vars J_NaCl and J_NaCl;
    enddef;-->
    <connection>
        <map_components component_1="Cell_concentration" component_2="NHE3"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="H_int" variable_2="H_int"/>
        <map_variables variable_1="H_Ext" variable_2="H_Ext"/>
        <!--        vars HCO3_i and HCO3_i;
        vars HCO3_m and HCO3_m;-->
        <!--        vars Cl_i and Cl_i;-->
        <map_variables variable_1="J_NHE3_Na" variable_2="J_NHE3_Na"/>
        <map_variables variable_1="J_NHE3_H" variable_2="J_NHE3_H"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="NaK_channel"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="NaK_channel"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
        <map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="NaK_channel"/>
        <map_variables variable_1="K_s" variable_2="K_s"/>
        <map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
        <map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="NaKCC"/>
        <map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
    </connection>
    <connection>
        <map_components component_1="Basol_voltage" component_2="NaK_channel"/>
        <map_variables variable_1="I_Na_NaK" variable_2="I_Na_NaK"/>
        <map_variables variable_1="I_K_NaK" variable_2="I_K_NaK"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="NBC"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <!--        vars Cl_i and Cl_i;-->
        <map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
        <map_variables variable_1="J_NBC_HCO3" variable_2="J_NBC_HCO3"/>
        <map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
        <map_variables variable_1="HCO3_s" variable_2="HCO3_s"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="NBC"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
    </connection>
    <connection>
        <map_components component_1="Basol_voltage" component_2="NBC"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
        <map_variables variable_1="I_NBC_HCO3" variable_2="I_NBC_HCO3"/>
        <map_variables variable_1="I_NBC_Na" variable_2="I_NBC_Na"/>
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="NBC"/>
        <map_variables variable_1="Na_s" variable_2="Na_s"/>
        <map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
    </connection>
    <!--    def map between Cell_concentration and NBC for
        vars Na_i and Na_i;
    enddef;-->
    <connection>
        <map_components component_1="Cell_concentration" component_2="Apical_Diffusion"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="J_D_Na" variable_2="J_D_Na"/>
        <map_variables variable_1="J_D_K" variable_2="J_D_K"/>
        <map_variables variable_1="J_D_Cl" variable_2="J_D_Cl"/>
    </connection>
    <!--        def map between Apical_concentrations and Cell_concentration for
            vars time and time;
        enddef;-->
    <connection>
        <map_components component_1="Apical_concentrations" component_2="Apical_Diffusion"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
        <map_variables variable_1="K_m" variable_2="K_m"/>
        <map_variables variable_1="Cl_m" variable_2="Cl_m"/>
        <!--vars F and F;-->
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="Apical_Diffusion"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Na" variable_2="z_Na"/>
        <map_variables variable_1="z_k" variable_2="z_k"/>
        <map_variables variable_1="z_Cl" variable_2="z_Cl"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
        <map_variables variable_1="P_Na" variable_2="P_Na"/>
        <map_variables variable_1="P_k" variable_2="P_k"/>
        <map_variables variable_1="P_Cl" variable_2="P_Cl"/>
    </connection>
    <!--    def map between parameters and Cell_concentration for
        vars A_Apical and A_Apical;
    enddef;-->
    <!--    def map between parameters and Apical_voltage for
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Apical_Diffusion" component_2="Apical_voltage"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="Apical_voltage"/>
        <map_variables variable_1="capacitance" variable_2="capacitance"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="Basolateral_Diffusion"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="J_D_Bl_Na" variable_2="J_D_Bl_Na"/>
        <map_variables variable_1="J_D_Bl_K" variable_2="J_D_Bl_K"/>
        <map_variables variable_1="J_D_Bl_Cl" variable_2="J_D_Bl_Cl"/>
    </connection>
    <!--    def map between Basolateral_concentrations and Cell_concentration for
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="Basolateral_Diffusion"/>
        <map_variables variable_1="Na_s" variable_2="Na_s"/>
        <map_variables variable_1="K_s" variable_2="K_s"/>
        <map_variables variable_1="Cl_s" variable_2="Cl_s"/>
        <!--vars F and F;-->
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="Basolateral_Diffusion"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Na" variable_2="z_Na"/>
        <map_variables variable_1="z_k" variable_2="z_k"/>
        <map_variables variable_1="z_Cl" variable_2="z_Cl"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Basol" variable_2="A_Basol"/>
        <map_variables variable_1="P_Bl_Na" variable_2="P_Bl_Na"/>
        <map_variables variable_1="P_Bl_k" variable_2="P_Bl_k"/>
        <map_variables variable_1="P_Bl_Cl" variable_2="P_Bl_Cl"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="Cell_concentration"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="NHE3"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
        <!--        vars time and time;
        vars R and R;
        vars T and T;-->
    </connection>
    <connection>
        <map_components component_1="Apical_voltage" component_2="NHE3"/>
        <!--vars v_sc and v_sc;-->
        <map_variables variable_1="I_NHE3_H" variable_2="I_NHE3_H"/>
        <map_variables variable_1="I_NHE3_Na" variable_2="I_NHE3_Na"/>
    </connection>
    <connection>
        <map_components component_1="Apical_concentrations" component_2="NHE3"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
    </connection>
    <!--    def map between Cell_concentration and NBC for
        vars Na_i and Na_i;
    enddef;-->
    <!--    def map between AE1 and NBC for
        vars HCO3_s and HCO3_s;
        vars HCO3_i and HCO3_i;
    enddef; -->
    <connection>
        <map_components component_1="parameters" component_2="AE1"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
        <!--vars A_Apical and A_Apical;-->
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="AE1"/>
        <map_variables variable_1="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
        <map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
        <map_variables variable_1="J_AE1_Cl" variable_2="J_AE1_Cl"/>
        <map_variables variable_1="J_AE1_HCO3" variable_2="J_AE1_HCO3"/>
    </connection>
    <connection>
        <map_components component_1="Apical_concentrations" component_2="AE1"/>
        <map_variables variable_1="Cl_m" variable_2="Cl_m"/>
    </connection>
    <!--    def map between Apical_voltage and AE1 for
        //vars v_sc and v_sc;-->
    <!--        vars I_AE1_Cl and I_AE1_Cl;
        vars I_AE1_HCO3 and I_AE1_HCO3;
    enddef;-->
    <!--    def map between Basolateral_concentrations and Basol_voltage for
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Basolateral_Diffusion" component_2="Basol_voltage"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="Basol_voltage"/>
        <map_variables variable_1="capacitance" variable_2="capacitance"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Apical_concentrations" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="Na_m" variable_2="Na_m"/>
        <map_variables variable_1="K_m" variable_2="K_m"/>
        <map_variables variable_1="Cl_m" variable_2="Cl_m"/>
        <!--vars F and F;-->
    </connection>
    <!--    def map between Apical_concentrations and parameters for
        vars time and time;-->
    <!--        //vars F and F;
    enddef;-->
    <connection>
        <map_components component_1="parameters" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="z_Na" variable_2="z_Na"/>
        <map_variables variable_1="z_k" variable_2="z_k"/>
        <map_variables variable_1="z_Cl" variable_2="z_Cl"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="A_Apical" variable_2="A_Apical"/>
        <map_variables variable_1="P_P_Na" variable_2="P_P_Na"/>
        <map_variables variable_1="P_P_k" variable_2="P_P_k"/>
        <map_variables variable_1="P_P_Cl" variable_2="P_P_Cl"/>
    </connection>
    <connection>
        <map_components component_1="Cell_concentration" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
    </connection>
    <!--    def map between Apical_concentrations and Paracellular_voltage for
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Paracellular_Diffusion" component_2="Paracellular_voltage"/>
        <!--        vars I_D_P_Na and I_D_P_Na;
        vars I_D_P_K and I_D_P_K;
        vars I_D_P_Cl and I_D_P_Cl;-->
        <map_variables variable_1="v_ms" variable_2="v_ms"/>
    </connection>
    <!--    def map between parameters and Paracellular_voltage for
        vars capacitance and capacitance;
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Basol_voltage" component_2="Paracellular_voltage"/>
        <map_variables variable_1="v_sc" variable_2="v_sc"/>
    </connection>
    <connection>
        <map_components component_1="Apical_voltage" component_2="Paracellular_voltage"/>
        <map_variables variable_1="v_mc" variable_2="v_mc"/>
    </connection>
    <connection>
        <map_components component_1="Apical_voltage" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
        <map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
        <map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
    </connection>
    <connection>
        <map_components component_1="Basol_voltage" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
        <map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
        <map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
    </connection>
    <connection>
        <map_components component_1="Basolateral_concentrations" component_2="Paracellular_Diffusion"/>
        <map_variables variable_1="Na_s" variable_2="Na_s"/>
        <map_variables variable_1="K_s" variable_2="K_s"/>
        <map_variables variable_1="Cl_s" variable_2="Cl_s"/>
        <map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
        <map_variables variable_1="J_D_P_K" variable_2="J_D_P_K"/>
        <map_variables variable_1="J_D_P_Cl" variable_2="J_D_P_Cl"/>
        <!--vars F and F;-->
    </connection>
    <!--    def map between interstitium_concentration and Blood_concentration for
        //vars J_G_Blood and J_G_Blood;-->
    <!--        vars glucose_s and glucose_s;
    enddef;-->
    <!--    def map between parameters and Blood_concentration for
        vars time and time;
    enddef;-->
    <connection>
        <map_components component_1="Cell_concentration" component_2="A_GLUT2"/>
        <map_variables variable_1="glucose_i" variable_2="glucose_i"/>
        <map_variables variable_1="J_A_GLUT" variable_2="J_A_GLUT"/>
    </connection>
    <connection>
        <map_components component_1="Apical_concentrations" component_2="A_GLUT2"/>
        <map_variables variable_1="glucose_m" variable_2="glucose_m"/>
        <!--vars J_GLUT and J_GLUT;-->
    </connection>
    <connection>
        <map_components component_1="parameters" component_2="A_GLUT2"/>
        <map_variables variable_1="v_cell" variable_2="v_cell"/>
    </connection>
    <!--    def map between Cell_concentration and Hydration for
        vars H_int and H_int;
        vars HCO3_i and HCO3_i;
        vars J_buf_C and J_buf_C;
        vars H_Ext and H_Ext;
        vars HCO3_m and HCO3_m;
        vars J_buf_A and J_buf_A;
    enddef;-->
    <!--    def map between parameters and Hydration for
        vars time and time;
        vars v_cell and v_cell;
    enddef;-->
    <!--    def map between Paracellular_Diffusion and Hydration for
        vars CO2_i and CO2_i;
        vars CO2_m and CO2_m;
        vars J_D_P_CO2 and J_D_P_CO2;
    enddef;-->
</model>