Location: ICC_Lees_Green_2014 @ 904cc945e815 / Simulations / Ano1_Lees_2021.cellml

Author:
Leyla <noroozbabaee@gmail.com>
Date:
2022-03-22 01:06:54+13:00
Desc:
Updatinf The Lees-Green-2014.
Permanent Source URI:
https://models.cellml.org/workspace/83c/rawfile/904cc945e8156453ad1a9330e8bd626a48c8b42d/Simulations/Ano1_Lees_2021.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="ICC" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
    <units name="conductance_units">
        <unit prefix="nano" units="siemens"/>
    </units>
    <units name="millifarads">
        <unit prefix="milli" units="farad"/>
    </units>
    <units name="voltage_units">
        <unit prefix="milli" units="volt"/>
    </units>
    <units name="volume_units">
        <unit units="litre"/>
    </units>
    <units name="Inverse_Voltage">
        <unit exponent="-1" units="voltage_units"/>
    </units>
    <units name="Inverse_Voltage_per_time_units">
        <unit units="Inverse_Voltage"/>
        <unit exponent="-1" units="time_units"/>
    </units>
    <units name="time_units">
        <unit units="second"/>
    </units>
    <units name="current_units">
        <unit units="voltage_units"/>
        <unit units="conductance_units"/>
    </units>
    <units name="rate_constants_units">
        <unit exponent="-1" units="time_units"/>
    </units>
    <units name="capacitance_units">
        <unit prefix="nano" units="farad"/>
    </units>
    <units name="Temperature_units">
        <unit units="kelvin"/>
    </units>
    <units name="Inverse_Temperature_units_times_conductance">
        <unit exponent="-1" units="Temperature_units"/>
        <unit units="conductance_units"/>
    </units>
    <units name="F_units">
        <unit exponent="-1" prefix="nano" units="mole"/>
        <unit prefix="micro" units="coulomb"/>
    </units>
    <units name="R_units">
        <unit prefix="nano" units="joule"/>
        <unit exponent="-1" prefix="nano" units="mole"/>
        <unit exponent="-1" units="kelvin"/>
    </units>
    <units name="millimolar">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="micromolar">
        <unit prefix="micro" units="mole"/>
        <unit exponent="-1" units="litre"/>
    </units>
    <units name="per_micromolar">
        <unit units="litre"/>
        <unit exponent="-1" prefix="micro" units="mole"/>
    </units>
    <units name="per_millimolar">
        <unit units="litre"/>
        <unit exponent="-1" prefix="milli" units="mole"/>
    </units>
    <units name="millimolar_per_second">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="litre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="per_micromolarsecond">
        <unit exponent="-1" prefix="micro" units="mole"/>
        <unit exponent="1" units="litre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="micromolar_per_second">
        <unit prefix="micro" units="mole"/>
        <unit exponent="-1" units="litre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <units name="millimolar_per_second_per_millivolt">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="litre"/>
        <unit exponent="-1" units="second"/>
        <unit exponent="-1" prefix="milli" units="volt"/>
    </units>
    <units name="millimole_to_nanomole">
        <unit prefix="nano" units="mole"/>
        <unit exponent="-1" prefix="milli" units="mole"/>
    </units>
    <units name="microcoulomb_to_picocoulomb">
        <unit prefix="pico" units="coulomb"/>
        <unit exponent="-1" prefix="micro" units="coulomb"/>
    </units>
    <units name="micrometer">
        <unit prefix="micro" units="meter"/>
    </units>
    <units name="micrometer2_per_time">
        <unit exponent="2" prefix="micro" units="meter"/>
        <unit exponent="-1" units="second"/>
    </units>
    <component name="Time">
        <variable cmeta:id="time" name="time" private_interface="out" public_interface="out" units="time_units"/>
    </component>
    <component name="Parameters">
        <variable initial_value="1e-9" name="p2m" public_interface="out" units="dimensionless"/>
        <variable initial_value="1e15" name="umc2L" public_interface="out" units="dimensionless"/>
        <variable initial_value="310" name="T" public_interface="out" units="Temperature_units"/>
        <variable initial_value="96.4846" name="F" public_interface="out" units="F_units"/>
        <variable initial_value="8.3144" name="R" public_interface="out" units="R_units"/>
        <variable initial_value="2" name="z_ca" 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="-1" name="z_cl" public_interface="out" units="dimensionless"/>
        <variable initial_value="2000" name="Ca_o" public_interface="out" units="micromolar"/>
        <variable initial_value="78000" name="Cl_i" public_interface="out" units="micromolar"/>
        <variable initial_value="166e3" name="Cl_o" public_interface="out" units="micromolar"/>
        <variable initial_value="140e3" name="K_i" public_interface="out" units="micromolar"/>
        <variable initial_value="5e3" name="K_o" public_interface="out" units="micromolar"/>
        <variable initial_value="30e3" name="Na_i" public_interface="out" units="micromolar"/>
        <variable initial_value="140e3" name="Na_o" public_interface="out" units="micromolar"/>
    </component>
    <component name="Volume">
        <variable initial_value="1.0e-12" name="cellVol" units="litre"/>
        <variable initial_value="0.1" name="cellPropER" units="dimensionless"/>
        <variable initial_value="0.7" name="cellPropCyto" units="dimensionless"/>
        <variable name="volER" public_interface="out" units="litre"/>
        <variable name="volCyto" public_interface="out" units="litre"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>volER</ci>
                <apply>
                    <times/>
                    <ci>cellVol</ci>
                    <ci>cellPropER</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>volCyto</ci>
                <apply>
                    <times/>
                    <ci>cellVol</ci>
                    <ci>cellPropCyto</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Nernst">
        <variable name="T" public_interface="in" units="Temperature_units"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="R" public_interface="in" units="R_units"/>
        <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="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="Cl_i" public_interface="in" units="micromolar"/>
        <variable name="Cl_o" public_interface="in" units="micromolar"/>
        <variable name="K_i" public_interface="in" units="micromolar"/>
        <variable name="K_o" public_interface="in" units="micromolar"/>
        <variable name="Na_i" public_interface="in" units="micromolar"/>
        <variable name="Na_o" public_interface="in" units="micromolar"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <variable name="Ca_o" public_interface="in" units="micromolar"/>
        <variable initial_value="0" name="E_NSCC" public_interface="out" units="voltage_units"/>
        <variable initial_value="0" name="E_NSV" public_interface="out" units="voltage_units"/>
        <variable initial_value="0" name="E_SOC" public_interface="out" units="voltage_units"/>
        <variable name="E_Ca" public_interface="out" units="voltage_units"/>
        <variable name="E_Cl" public_interface="out" units="voltage_units"/>
        <variable name="E_Na" public_interface="out" units="voltage_units"/>
        <variable name="E_K" public_interface="out" units="voltage_units"/>
        <!--E_Cl = R*T/(F*z_cl)*ln(Cl_o/Cl_i);-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>E_Ca</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>F</ci>
                            <ci>z_ca</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Ca_o</ci>
                            <ci>Ca_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>E_Na</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>F</ci>
                            <ci>z_na</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Na_o</ci>
                            <ci>Na_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>E_K</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>T</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>F</ci>
                            <ci>z_k</ci>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>K_o</ci>
                            <ci>K_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>E_Cl</ci>
                <cn cellml:units="dimensionless">0</cn>
            </apply>
        </math>
    </component>
    <component name="d_Ano1">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="J_SOC" public_interface="in" units="micromolar_per_second"/>
        <variable initial_value="250" name="Dc" units="micrometer2_per_time"/>
        <variable initial_value="75" name="Dm" units="micrometer2_per_time"/>
        <variable initial_value="50" name="Bm" units="micromolar"/>
        <variable initial_value="1" name="Km" units="micromolar"/>
        <variable initial_value="0.01248" name="kc" public_interface="out" units="Inverse_Voltage"/>
        <variable initial_value="50" name="nSOC" public_interface="out" units="dimensionless"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <!--var Ca: micromolar {pub: out};-->
        <variable name="rad" units="micrometer"/>
        <variable name="mouthCa" units="micromolar_per_second"/>
        <variable name="localCa" units="micromolar"/>
        <variable name="d_Ano1_tau" units="time_units"/>
        <variable initial_value="10" name="d_Ano1_tscale" units="dimensionless"/>
        <variable initial_value="0.05" name="colocalRadius" units="micrometer"/>
        <variable name="d_Ano1_inf" units="dimensionless"/>
        <variable initial_value="0.01" name="d_Ano1" public_interface="out" units="dimensionless"/>
        <variable initial_value="2" name="Ano1_n" units="dimensionless"/>
        <variable initial_value="-100" name="Ano1_Vh" units="voltage_units"/>
        <variable initial_value="0.0156" name="Ano1_s" units="voltage_units"/>
        <variable initial_value="1" name="localCa_clamp" units="micromolar"/>
        <variable initial_value="0" name="localCa_bl" units="micromolar"/>
        <variable initial_value="0" name="localCa_init" units="micromolar"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>rad</ci>
                <ci>colocalRadius</ci>
            </apply>
            <apply>
                <eq/>
                <ci>mouthCa</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <cn cellml:units="dimensionless" type="e-notation">1<sep/>15</cn>
                        <ci>J_SOC</ci>
                        <ci>volCyto</ci>
                    </apply>
                    <ci>nSOC</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>localCa</ci>
                <piecewise>
                    <piece>
                        <ci>localCa_init</ci>
                        <apply>
                            <leq/>
                            <ci>time</ci>
                            <cn cellml:units="time_units">1</cn>
                        </apply>
                    </piece>
                    <piece>
                        <ci>localCa_clamp</ci>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">1</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">9</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <ci>localCa_bl</ci>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>d_Ano1_tau</ci>
                <apply>
                    <times/>
                    <ci>d_Ano1_tscale</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <cn cellml:units="time_units">81.63</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="per_micromolar">0.57</cn>
                                    </apply>
                                    <ci>localCa</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="time_units">76.17</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <cn cellml:units="time_units">0.05374</cn>
                                    </apply>
                                    <ci>localCa</ci>
                                </apply>
                            </apply>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <ci>Vm</ci>
                                    <apply>
                                        <times/>
                                        <cn cellml:units="Inverse_Voltage">70.3</cn>
                                        <apply>
                                            <exp/>
                                            <apply>
                                                <times/>
                                                <cn cellml:units="per_micromolar">0.153</cn>
                                                <ci>localCa</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <cn cellml:units="dimensionless" type="e-notation">1<sep/>-3</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>d_Ano1_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <ci>Ano1_Vh</ci>
                                        <ci>Vm</ci>
                                    </apply>
                                    <ci>Ano1_s</ci>
                                </apply>
                            </apply>
                        </apply>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <ci>localCa</ci>
                                    <apply>
                                        <times/>
                                        <cn cellml:units="micromolar">1.3</cn>
                                        <apply>
                                            <exp/>
                                            <apply>
                                                <times/>
                                                <apply>
                                                    <minus/>
                                                    <ci>kc</ci>
                                                </apply>
                                                <ci>Vm</ci>
                                            </apply>
                                        </apply>
                                    </apply>
                                </apply>
                                <apply>
                                    <minus/>
                                    <ci>Ano1_n</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_Ano1</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_Ano1_inf</ci>
                        <ci>d_Ano1</ci>
                    </apply>
                    <ci>d_Ano1_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_Ano1">
        <variable initial_value="0.5" name="g_Ano1" units="conductance_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="J_SOC" private_interface="out" public_interface="in" units="micromolar_per_second"/>
        <variable name="Ca_i" private_interface="out" public_interface="in" units="micromolar"/>
        <variable name="volCyto" private_interface="out" public_interface="in" units="litre"/>
        <variable name="E_Cl" public_interface="in" units="voltage_units"/>
        <variable name="I_Ano1" public_interface="out" units="current_units"/>
        <variable name="d_Ano1" private_interface="in" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_Ano1</ci>
                <apply>
                    <times/>
                    <ci>g_Ano1</ci>
                    <ci>d_Ano1</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Cl</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_CaT">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-26" name="d_CaT_Vh" units="voltage_units"/>
        <variable initial_value="-6" name="d_CaT_S" units="voltage_units"/>
        <variable initial_value="0.006" name="d_CaT_tau" units="time_units"/>
        <variable name="d_CaT_inf" units="dimensionless"/>
        <variable initial_value="0.001271016263081" name="d_CaT" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_CaT_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_CaT_Vh</ci>
                                </apply>
                                <ci>d_CaT_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_CaT</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_CaT_inf</ci>
                        <ci>d_CaT</ci>
                    </apply>
                    <ci>d_CaT_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="f_CaT">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-66" name="f_CaT_Vh" units="voltage_units"/>
        <variable initial_value="6" name="f_CaT_S" units="voltage_units"/>
        <variable initial_value="0.04" name="f_CaT_tau" units="time_units"/>
        <variable name="f_CaT_inf" units="dimensionless"/>
        <variable initial_value="0.5" name="f_CaT" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f_CaT_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>f_CaT_Vh</ci>
                                </apply>
                                <ci>f_CaT_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f_CaT</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>f_CaT_inf</ci>
                        <ci>f_CaT</ci>
                    </apply>
                    <ci>f_CaT_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <!--Calculate ion currents CaT for H(NaV,NSV,CaV)-->
    <component name="I_CaT">
        <variable initial_value="0" name="g_CaT" units="conductance_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="p2m" public_interface="in" units="dimensionless"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="E_Ca" public_interface="in" units="voltage_units"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="f_CaT" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="d_CaT" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_CaT" public_interface="out" units="current_units"/>
        <variable name="J_CaT" public_interface="out" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_CaT</ci>
                <apply>
                    <times/>
                    <ci>g_CaT</ci>
                    <ci>d_CaT</ci>
                    <ci>f_CaT</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Ca</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_CaT</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <minus/>
                            <ci>p2m</ci>
                        </apply>
                        <ci>I_CaT</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>z_ca</ci>
                        <ci>F</ci>
                        <ci>volCyto</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_CaV">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-26" name="d_CaV_Vh" units="voltage_units"/>
        <variable initial_value="-6" name="d_CaV_S" units="voltage_units"/>
        <variable initial_value="0.006" name="d_CaV_tau" units="time_units"/>
        <variable name="d_CaV_inf" units="dimensionless"/>
        <variable initial_value="0.001271016263081" name="d_CaV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_CaV_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_CaV_Vh</ci>
                                </apply>
                                <ci>d_CaV_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_CaV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_CaV_inf</ci>
                        <ci>d_CaV</ci>
                    </apply>
                    <ci>d_CaV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_CaV">
        <variable initial_value="0" name="g_CaV" units="conductance_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="p2m" public_interface="in" units="dimensionless"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="E_Ca" public_interface="in" units="voltage_units"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="d_CaV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_CaV" public_interface="out" units="current_units"/>
        <variable name="J_CaV" public_interface="out" units="micromolar_per_second"/>
        <variable name="volER" public_interface="in" units="litre"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_CaV</ci>
                <apply>
                    <times/>
                    <ci>g_CaV</ci>
                    <ci>d_CaV</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Ca</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_CaV</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <minus/>
                            <ci>p2m</ci>
                        </apply>
                        <ci>I_CaV</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>z_ca</ci>
                        <ci>F</ci>
                        <ci>volER</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_NaV">
        <variable name="time" public_interface="in" units="time_units"/>
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable initial_value="-39" name="d_NaV_Vh" units="voltage_units"/>
        <variable initial_value="-5.6" name="d_NaV_S" units="voltage_units"/>
        <variable initial_value="50" name="d_NaV_t1" units="rate_constants_units"/>
        <variable initial_value="1050" name="d_NaV_t2" units="rate_constants_units"/>
        <variable initial_value="-0.06" name="d_NaV_t3" units="Inverse_Voltage"/>
        <variable initial_value="-20" name="d_NaV_t4" units="voltage_units"/>
        <variable name="d_NaV_tau" units="time_units"/>
        <variable name="d_NaV_inf" units="dimensionless"/>
        <variable initial_value="0.0006607800216060858" name="d_NaV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_NaV_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_NaV_Vh</ci>
                                </apply>
                                <ci>d_NaV_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>d_NaV_tau</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <ci>d_NaV_t1</ci>
                        <apply>
                            <divide/>
                            <apply>
                                <minus/>
                                <ci>d_NaV_t2</ci>
                                <ci>d_NaV_t1</ci>
                            </apply>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <times/>
                                        <ci>d_NaV_t3</ci>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>d_NaV_t4</ci>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_NaV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_NaV_inf</ci>
                        <ci>d_NaV</ci>
                    </apply>
                    <ci>d_NaV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="f_NaV">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-55" name="f_NaV_Vh" units="voltage_units"/>
        <variable initial_value="-4.5" name="f_NaV_S" units="voltage_units"/>
        <variable initial_value="1" name="f_NaV_tau" units="time_units"/>
        <variable name="f_NaV_inf" units="dimensionless"/>
        <variable initial_value="1" name="f_NaV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f_NaV_inf</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">0.53</cn>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>f_NaV_Vh</ci>
                                        </apply>
                                        <ci>f_NaV_S</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                    <cn cellml:units="dimensionless">0.53</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f_NaV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>f_NaV_inf</ci>
                        <ci>f_NaV</ci>
                    </apply>
                    <ci>f_NaV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_NaV">
        <variable initial_value="0" name="g_NaV" units="conductance_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="f_NaV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="d_NaV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_NaV" public_interface="out" units="current_units"/>
        <variable name="E_Na" public_interface="in" units="voltage_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_NaV</ci>
                <apply>
                    <times/>
                    <ci>g_NaV</ci>
                    <ci>d_NaV</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Na</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--Ion current specific to H(NaV)-->
        <!--I_NaV = g_NaV*d_NaV*f_NaV*(Vm-E_Na);-->
        <!--Ion current specific to L(NaV)-->
    </component>
    <component name="P_SOC">
        <variable initial_value="200" name="SOC_h" public_interface="out" units="micromolar"/>
        <variable name="Ca_er" public_interface="in" units="micromolar"/>
        <variable initial_value="8" name="SOC_n" public_interface="out" units="dimensionless"/>
        <variable name="P_SOC" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>P_SOC</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>Ca_er</ci>
                                <ci>SOC_h</ci>
                            </apply>
                            <ci>SOC_n</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <!--H(NaV,CaV)-->
    <component name="I_SOC">
        <variable initial_value="0.0" name="g_SOC" units="conductance_units"/>
        <variable name="Ca_er" private_interface="out" public_interface="in" units="micromolar"/>
        <variable name="p2m" public_interface="in" units="dimensionless"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="E_Ca" public_interface="in" units="voltage_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable initial_value="1.0" name="SOCPropCa" public_interface="out" units="dimensionless"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="P_SOC" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_SOC" public_interface="out" units="current_units"/>
        <variable name="J_SOC" public_interface="out" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_SOC</ci>
                <apply>
                    <times/>
                    <ci>g_SOC</ci>
                    <ci>P_SOC</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Ca</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_SOC</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <minus/>
                            <ci>p2m</ci>
                        </apply>
                        <ci>SOCPropCa</ci>
                        <ci>I_SOC</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>z_ca</ci>
                        <ci>F</ci>
                        <ci>volCyto</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_KV">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-20" name="d_KV_Vh" units="voltage_units"/>
        <variable initial_value="-4.6" name="d_KV_S" units="voltage_units"/>
        <variable initial_value="0.003" name="d_KV_tau" units="time_units"/>
        <variable initial_value="0.01" name="d_KV" public_interface="out" units="dimensionless"/>
        <variable name="d_KV_inf" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_KV_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_KV_Vh</ci>
                                </apply>
                                <ci>d_KV_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_KV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_KV_inf</ci>
                        <ci>d_KV</ci>
                    </apply>
                    <ci>d_KV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="f_KV">
        <variable name="time" public_interface="in" units="time_units"/>
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable initial_value="-45" name="f_KV_Vh" units="voltage_units"/>
        <variable initial_value="4.4" name="f_KV_S" units="voltage_units"/>
        <variable initial_value="1" name="f_KV_a" units="dimensionless"/>
        <variable initial_value="0.045" name="f_KV_tau" units="time_units"/>
        <variable name="f_KV_inf" units="dimensionless"/>
        <variable initial_value="1" name="f_KV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f_KV_inf</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <divide/>
                            <ci>f_KV_a</ci>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>f_KV_Vh</ci>
                                        </apply>
                                        <ci>f_KV_S</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                    <ci>f_KV_a</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f_KV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>f_KV_inf</ci>
                        <ci>f_KV</ci>
                    </apply>
                    <ci>f_KV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <!--KV (kt), L(NaV)-->
    <component name="I_KV">
        <variable initial_value="0" name="g_KV" units="conductance_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="E_K" public_interface="in" units="voltage_units"/>
        <variable name="f_KV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="d_KV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_KV" public_interface="out" units="current_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_KV</ci>
                <apply>
                    <times/>
                    <ci>g_KV</ci>
                    <ci>d_KV</ci>
                    <ci>f_KV</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_K</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_NSV">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-30" name="d_NSV_Vh" units="voltage_units"/>
        <variable initial_value="-5" name="d_NSV_S" units="voltage_units"/>
        <variable initial_value="0.005" name="d_NSV_tau" units="time_units"/>
        <variable name="d_NSV_inf" units="dimensionless"/>
        <variable initial_value="0.001271016263081" name="d_NSV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_NSV_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_NSV_Vh</ci>
                                </apply>
                                <ci>d_NSV_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_NSV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_NSV_inf</ci>
                        <ci>d_NSV</ci>
                    </apply>
                    <ci>d_NSV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="f_NSV">
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="time" public_interface="in" units="time_units"/>
        <variable initial_value="-50" name="f_NSV_Vh" units="voltage_units"/>
        <variable initial_value="5" name="f_NSV_S" units="voltage_units"/>
        <variable initial_value="0.8" name="f_NSV_tau" units="time_units"/>
        <variable name="f_NSV_inf" units="dimensionless"/>
        <variable initial_value="0.5" name="f_NSV" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f_NSV_inf</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">0.3</cn>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>f_NSV_Vh</ci>
                                        </apply>
                                        <ci>f_NSV_S</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                    <cn cellml:units="dimensionless">0.3</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f_NSV</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>f_NSV_inf</ci>
                        <ci>f_NSV</ci>
                    </apply>
                    <ci>f_NSV_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <!--H(NSV)-->
    <component name="I_NSV">
        <variable initial_value="0" name="g_NSV" units="conductance_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="E_NSV" public_interface="in" units="voltage_units"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="p2m" public_interface="in" units="dimensionless"/>
        <variable name="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="f_NSV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="d_NSV" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="I_NSV" public_interface="out" units="current_units"/>
        <variable name="J_NSV" public_interface="out" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_NSV</ci>
                <apply>
                    <times/>
                    <ci>g_NSV</ci>
                    <ci>d_NSV</ci>
                    <ci>f_NSV</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_NSV</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_NSV</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <minus/>
                            <ci>p2m</ci>
                        </apply>
                        <ci>I_NSV</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>z_ca</ci>
                        <ci>F</ci>
                        <ci>volCyto</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_BK">
        <variable initial_value="0" name="g_BK" units="conductance_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="E_K" public_interface="in" units="voltage_units"/>
        <variable name="I_BK" public_interface="out" units="current_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_BK</ci>
                <apply>
                    <times/>
                    <ci>g_BK</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_K</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="d_KERG">
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable initial_value="-35" name="d_KERG_Vh" units="voltage_units"/>
        <variable initial_value="-7.1" name="d_KERG_S" units="voltage_units"/>
        <variable initial_value="0.15" name="d_KERG_ta" units="rate_constants_units"/>
        <variable initial_value="-60" name="d_KERG_tV" units="voltage_units"/>
        <variable initial_value="20" name="d_KERG_tS" units="voltage_units"/>
        <variable name="d_KERG_tau" units="time_units"/>
        <variable name="d_KERG_inf" units="dimensionless"/>
        <variable initial_value="0.001" name="d_KERG" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>d_KERG_tau</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <times/>
                                        <ci>d_KERG_ta</ci>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>d_KERG_tV</ci>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>d_KERG_inf</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <ci>Vm</ci>
                                    <ci>d_KERG_Vh</ci>
                                </apply>
                                <ci>d_KERG_S</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>d_KERG</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>d_KERG_inf</ci>
                        <ci>d_KERG</ci>
                    </apply>
                    <ci>d_KERG_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="f_KERG">
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable initial_value="-42" name="f_KERG_Vh" units="voltage_units"/>
        <variable initial_value="4" name="f_KERG_S" units="voltage_units"/>
        <variable initial_value="1.0" name="f_KERG_a" units="dimensionless"/>
        <variable initial_value="0.003" name="f_KERG_tau" units="time_units"/>
        <variable name="f_KERG_inf" units="dimensionless"/>
        <variable initial_value="0.001" name="f_KERG" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f_KERG_inf</ci>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <apply>
                            <divide/>
                            <ci>f_KERG_a</ci>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <apply>
                                    <exp/>
                                    <apply>
                                        <divide/>
                                        <apply>
                                            <minus/>
                                            <ci>Vm</ci>
                                            <ci>f_KERG_Vh</ci>
                                        </apply>
                                        <ci>f_KERG_S</ci>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                    <ci>f_KERG_a</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f_KERG</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>f_KERG_inf</ci>
                        <ci>f_KERG</ci>
                    </apply>
                    <ci>f_KERG_tau</ci>
                </apply>
            </apply>
        </math>
    </component>
    <!--L(NaV)-->
    <component name="I_KERG">
        <variable initial_value="0.0" name="g_KERG" units="conductance_units"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="f_KERG" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="d_KERG" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="E_K" public_interface="in" units="voltage_units"/>
        <variable name="I_KERG" public_interface="out" units="current_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_KERG</ci>
                <apply>
                    <times/>
                    <ci>g_KERG</ci>
                    <ci>d_KERG</ci>
                    <ci>f_KERG</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_K</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_stim">
        <variable name="I_stim" public_interface="out" units="current_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_stim</ci>
                <piecewise>
                    <piece>
                        <cn cellml:units="current_units">0</cn>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">100</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">900</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <cn cellml:units="current_units">0</cn>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
    </component>
    <component name="P_NSCC">
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <variable name="P_NSCC" public_interface="out" units="dimensionless"/>
        <variable initial_value="4" name="n_NSCC" units="dimensionless"/>
        <variable initial_value="1.8" name="Ca_NSCC" units="micromolar"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>P_NSCC</ci>
                <apply>
                    <divide/>
                    <apply>
                        <power/>
                        <ci>Ca_i</ci>
                        <ci>n_NSCC</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <power/>
                            <ci>Ca_NSCC</ci>
                            <ci>n_NSCC</ci>
                        </apply>
                        <apply>
                            <power/>
                            <ci>Ca_i</ci>
                            <ci>n_NSCC</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_NSCC">
        <variable initial_value="0.0" name="g_NSCC" units="rate_constants_units"/>
        <variable name="Ca_i" private_interface="out" public_interface="in" units="micromolar"/>
        <variable name="I_NSCC" public_interface="out" units="current_units"/>
        <variable name="P_NSCC" private_interface="in" units="dimensionless"/>
        <variable name="Vm" private_interface="out" public_interface="in" units="voltage_units"/>
        <variable name="E_NSCC" public_interface="in" units="voltage_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_NSCC</ci>
                <apply>
                    <times/>
                    <ci>g_NSCC</ci>
                    <ci>P_NSCC</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_NSCC</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="J_SERCA">
        <variable initial_value="0.0" name="Ke" units="micromolar"/>
        <variable name="J_SERCA" public_interface="out" units="micromolar_per_second"/>
        <variable initial_value="80" name="Ve" units="micromolar_per_second"/>
        <variable initial_value="0.01" name="fc" public_interface="out" units="dimensionless"/>
        <variable initial_value="1" name="fe" public_interface="out" units="dimensionless"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_SERCA</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>Ve</ci>
                        <apply>
                            <power/>
                            <ci>Ca_i</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <apply>
                            <power/>
                            <ci>Ke</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                        <apply>
                            <power/>
                            <ci>Ca_i</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="Ca_er">
        <variable name="time" public_interface="in" units="time_units"/>
        <variable name="fe" public_interface="in" units="dimensionless"/>
        <variable name="volER" public_interface="in" units="litre"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable initial_value="290" name="Ca_er" public_interface="out" units="micromolar"/>
        <variable name="J_SERCA" public_interface="in" units="micromolar_per_second"/>
        <variable name="J_IPR" public_interface="in" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Ca_er</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>fe</ci>
                        <apply>
                            <minus/>
                            <ci>J_SERCA</ci>
                            <ci>J_IPR</ci>
                        </apply>
                        <ci>volCyto</ci>
                    </apply>
                    <ci>volER</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="y_bind">
        <variable name="time" public_interface="in" units="time_units"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <variable initial_value="4000" name="k1" public_interface="out" units="per_micromolarsecond"/>
        <variable initial_value="2.0" name="k2" units="per_micromolarsecond"/>
        <variable initial_value="4000" name="k3" units="per_micromolarsecond"/>
        <variable initial_value="2.0" name="k4" units="per_micromolarsecond"/>
        <variable initial_value="200" name="k5" public_interface="out" units="per_micromolarsecond"/>
        <variable initial_value="520" name="k_1" units="rate_constants_units"/>
        <variable initial_value="2.10" name="k_2" units="rate_constants_units"/>
        <variable initial_value="3772" name="k_3" units="rate_constants_units"/>
        <variable initial_value="0.29" name="k_4" units="rate_constants_units"/>
        <variable initial_value="16.4" name="k_5" units="rate_constants_units"/>
        <variable name="IP3" public_interface="in" units="micromolar"/>
        <variable name="K1" public_interface="out" units="micromolar"/>
        <variable name="K2" public_interface="out" units="micromolar"/>
        <variable name="K3" public_interface="out" units="micromolar"/>
        <variable name="K4" public_interface="out" units="micromolar"/>
        <variable name="K5" public_interface="out" units="micromolar"/>
        <variable initial_value="0.4" name="y" public_interface="out" units="dimensionless"/>
        <variable name="phi1" units="rate_constants_units"/>
        <variable name="phi2" units="rate_constants_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>K1</ci>
                <apply>
                    <divide/>
                    <ci>k_1</ci>
                    <ci>k1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K2</ci>
                <apply>
                    <divide/>
                    <ci>k_2</ci>
                    <ci>k2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K3</ci>
                <apply>
                    <divide/>
                    <ci>k_3</ci>
                    <ci>k3</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K4</ci>
                <apply>
                    <divide/>
                    <ci>k_4</ci>
                    <ci>k4</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K5</ci>
                <apply>
                    <divide/>
                    <ci>k_5</ci>
                    <ci>k5</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>phi1</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <ci>k_4</ci>
                                <ci>K2</ci>
                                <ci>K1</ci>
                            </apply>
                            <apply>
                                <times/>
                                <ci>k_2</ci>
                                <ci>K4</ci>
                                <ci>IP3</ci>
                            </apply>
                        </apply>
                        <ci>Ca_i</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>K4</ci>
                        <ci>K2</ci>
                        <apply>
                            <plus/>
                            <ci>K1</ci>
                            <ci>IP3</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>phi2</ci>
                <apply>
                    <divide/>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>k_2</ci>
                            <ci>IP3</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>k_4</ci>
                            <ci>K3</ci>
                        </apply>
                    </apply>
                    <apply>
                        <plus/>
                        <ci>K3</ci>
                        <ci>IP3</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>y</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <ci>phi1</ci>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>y</ci>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <ci>phi2</ci>
                        <ci>y</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="J_IPR">
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable name="Ca_i" private_interface="out" public_interface="in" units="micromolar"/>
        <variable name="Ca_er" public_interface="in" units="micromolar"/>
        <variable initial_value="7" name="kipr" units="rate_constants_units"/>
        <variable name="Pipr" units="dimensionless"/>
        <variable initial_value="0.01" name="Jer" units="rate_constants_units"/>
        <variable initial_value="0.5" name="IP3" private_interface="out" units="micromolar"/>
        <variable name="K1" private_interface="in" units="micromolar"/>
        <variable name="K5" private_interface="in" units="micromolar"/>
        <variable name="y" private_interface="in" public_interface="out" units="dimensionless"/>
        <variable name="J_SERCA" public_interface="in" units="micromolar_per_second"/>
        <variable name="J_IPR" public_interface="out" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Pipr</ci>
                <apply>
                    <power/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>IP3</ci>
                            <ci>Ca_i</ci>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless">1</cn>
                                <ci>y</ci>
                            </apply>
                        </apply>
                        <apply>
                            <times/>
                            <apply>
                                <plus/>
                                <ci>IP3</ci>
                                <ci>K1</ci>
                            </apply>
                            <apply>
                                <plus/>
                                <ci>Ca_i</ci>
                                <ci>K5</ci>
                            </apply>
                        </apply>
                    </apply>
                    <cn cellml:units="dimensionless">3</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_IPR</ci>
                <apply>
                    <times/>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <ci>kipr</ci>
                            <ci>Pipr</ci>
                        </apply>
                        <ci>Jer</ci>
                    </apply>
                    <apply>
                        <minus/>
                        <ci>Ca_er</ci>
                        <ci>Ca_i</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_NCX">
        <variable initial_value="0.0" name="switch_NCX" units="dimensionless"/>
        <variable name="z_ca" public_interface="in" units="dimensionless"/>
        <variable name="F" public_interface="in" units="F_units"/>
        <variable name="volCyto" public_interface="in" units="litre"/>
        <variable name="p2m" public_interface="in" units="dimensionless"/>
        <variable name="SOCPropCa" public_interface="in" units="dimensionless"/>
        <variable initial_value="225" name="v_NCX" units="micromolar_per_second"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <variable name="J_NCX" public_interface="out" units="micromolar_per_second"/>
        <variable name="I_NCX" public_interface="out" units="current_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_NCX</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>switch_NCX</ci>
                        <ci>v_NCX</ci>
                        <ci>Ca_i</ci>
                    </apply>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">3.6</cn>
                        <ci>Ca_i</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_NCX</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>switch_NCX</ci>
                            <ci>z_ca</ci>
                            <ci>F</ci>
                            <ci>volCyto</ci>
                        </apply>
                        <apply>
                            <minus/>
                            <ci>p2m</ci>
                        </apply>
                    </apply>
                    <ci>SOCPropCa</ci>
                    <ci>J_NCX</ci>
                </apply>
            </apply>
        </math>
    </component>
    <component name="J_PMCA">
        <variable initial_value="0.0" name="g_PMCA" units="dimensionless"/>
        <variable name="Ca_i" public_interface="in" units="micromolar"/>
        <variable initial_value="50" name="J_PMCA_max" units="micromolar_per_second"/>
        <variable initial_value="0.1" name="K_PMCA" units="micromolar"/>
        <variable initial_value="2" name="n_PMCA" units="dimensionless"/>
        <variable name="J_PMCA" public_interface="out" units="micromolar_per_second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_PMCA</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>g_PMCA</ci>
                        <ci>J_PMCA_max</ci>
                        <cn cellml:units="dimensionless">1</cn>
                    </apply>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>K_PMCA</ci>
                                <ci>Ca_i</ci>
                            </apply>
                            <ci>n_PMCA</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="I_BNa">
        <variable initial_value="0" name="g_BNa" units="conductance_units"/>
        <variable name="Vm" public_interface="in" units="voltage_units"/>
        <variable name="E_Na" public_interface="in" units="voltage_units"/>
        <variable name="I_BNa" public_interface="out" units="current_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_BNa</ci>
                <apply>
                    <times/>
                    <ci>g_BNa</ci>
                    <apply>
                        <minus/>
                        <ci>Vm</ci>
                        <ci>E_Na</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="ICC_Membrane">
        <variable name="Vm" private_interface="out" public_interface="out" units="voltage_units"/>
        <variable name="time" private_interface="out" public_interface="in" units="time_units"/>
        <variable initial_value="0.025" name="Cm" units="capacitance_units"/>
        <variable name="fc" private_interface="in" units="dimensionless"/>
        <variable name="f_CaT" private_interface="in" units="dimensionless"/>
        <variable name="d_CaT" private_interface="in" units="dimensionless"/>
        <variable name="I_CaT" private_interface="in" units="current_units"/>
        <variable name="J_CaT" private_interface="in" units="micromolar_per_second"/>
        <variable name="d_CaV" private_interface="in" units="dimensionless"/>
        <variable name="I_CaV" private_interface="in" units="current_units"/>
        <variable name="J_CaV" private_interface="in" units="micromolar_per_second"/>
        <variable name="f_NaV" private_interface="in" units="dimensionless"/>
        <variable name="d_NaV" private_interface="in" units="dimensionless"/>
        <variable name="I_NaV" private_interface="in" units="current_units"/>
        <variable name="f_NSV" private_interface="in" units="dimensionless"/>
        <variable name="d_NSV" private_interface="in" units="dimensionless"/>
        <variable name="I_NSV" private_interface="in" units="current_units"/>
        <variable name="J_NSV" private_interface="in" units="micromolar_per_second"/>
        <variable name="I_NCX" private_interface="in" units="current_units"/>
        <variable name="J_NCX" private_interface="in" units="micromolar_per_second"/>
        <variable name="f_KERG" private_interface="in" units="dimensionless"/>
        <variable name="d_KERG" private_interface="in" units="dimensionless"/>
        <variable name="I_KERG" private_interface="in" units="current_units"/>
        <variable name="f_KV" private_interface="in" units="dimensionless"/>
        <variable name="d_KV" private_interface="in" units="dimensionless"/>
        <variable name="I_NSCC" private_interface="in" units="current_units"/>
        <variable name="I_stim" private_interface="in" units="current_units"/>
        <variable name="I_KV" private_interface="in" units="current_units"/>
        <variable name="I_BNa" private_interface="in" units="current_units"/>
        <variable name="I_BK" private_interface="in" units="current_units"/>
        <variable name="d_Ano1" private_interface="in" units="dimensionless"/>
        <variable name="I_Ano1" private_interface="in" units="current_units"/>
        <variable name="I_SOC" private_interface="in" units="current_units"/>
        <variable name="J_SOC" private_interface="in" units="micromolar_per_second"/>
        <variable name="J_IPR" private_interface="in" units="micromolar_per_second"/>
        <variable name="J_PMCA" private_interface="in" units="micromolar_per_second"/>
        <variable name="J_SERCA" private_interface="in" units="micromolar_per_second"/>
        <variable name="Ca_er" private_interface="in" units="micromolar"/>
        <variable initial_value="0.14" name="Ca_i" private_interface="out" public_interface="out" units="micromolar"/>
        <variable name="Iion" public_interface="out" units="current_units"/>
        <variable initial_value="-100" name="V_clamp" units="voltage_units"/>
        <variable initial_value="-100" name="V_bl" units="voltage_units"/>
        <variable initial_value="0" name="V_init" units="voltage_units"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>Ca_i</ci>
                </apply>
                <cn cellml:units="dimensionless">0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Vm</ci>
                <piecewise>
                    <piece>
                        <ci>V_init</ci>
                        <apply>
                            <leq/>
                            <ci>time</ci>
                            <cn cellml:units="time_units">1</cn>
                        </apply>
                    </piece>
                    <piece>
                        <ci>V_clamp</ci>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">1</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <cn cellml:units="time_units">9</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <ci>V_bl</ci>
                    </otherwise>
                </piecewise>
            </apply>
            <apply>
                <eq/>
                <ci>Iion</ci>
                <cn cellml:units="dimensionless">0</cn>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="Parameters" component_2="Nernst"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <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="Cl_i" variable_2="Cl_i"/>
        <map_variables variable_1="Cl_o" variable_2="Cl_o"/>
        <map_variables variable_1="K_i" variable_2="K_i"/>
        <map_variables variable_1="K_o" variable_2="K_o"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
        <map_variables variable_1="Na_o" variable_2="Na_o"/>
        <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="R" variable_2="R"/>
    </connection>
    <connection>
        <map_components component_1="d_CaT" component_2="I_CaT"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_CaT" variable_2="d_CaT"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="f_CaT" component_2="I_CaT"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="f_CaT" variable_2="f_CaT"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_CaT"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_CaT"/>
        <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
    </connection>
    <connection>
        <map_components component_1="Parameters" component_2="I_CaT"/>
        <map_variables variable_1="p2m" variable_2="p2m"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="d_CaV" component_2="I_CaV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_CaV" variable_2="d_CaV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_CaV"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
        <map_variables variable_1="volER" variable_2="volER"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_CaV"/>
        <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
    </connection>
    <connection>
        <map_components component_1="Parameters" component_2="I_CaV"/>
        <map_variables variable_1="p2m" variable_2="p2m"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="Parameters" component_2="I_NSV"/>
        <map_variables variable_1="p2m" variable_2="p2m"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="d_NSV" component_2="I_NSV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_NSV" variable_2="d_NSV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="f_NSV" component_2="I_NSV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="f_NSV" variable_2="f_NSV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_NSV"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_NSV"/>
        <map_variables variable_1="E_NSV" variable_2="E_NSV"/>
    </connection>
    <connection>
        <map_components component_1="d_NaV" component_2="I_NaV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_NaV" variable_2="d_NaV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="f_NaV" component_2="I_NaV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="f_NaV" variable_2="f_NaV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_NaV"/>
        <map_variables variable_1="E_Na" variable_2="E_Na"/>
    </connection>
    <connection>
        <map_components component_1="d_KV" component_2="I_KV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_KV" variable_2="d_KV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="f_KV" component_2="I_KV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="f_KV" variable_2="f_KV"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_KV"/>
        <map_variables variable_1="E_K" variable_2="E_K"/>
    </connection>
    <connection>
        <map_components component_1="d_KERG" component_2="I_KERG"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_KERG" variable_2="d_KERG"/>
    </connection>
    <connection>
        <map_components component_1="f_KERG" component_2="I_KERG"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="f_KERG" variable_2="f_KERG"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_KERG"/>
        <map_variables variable_1="E_K" variable_2="E_K"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_SOC"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
    </connection>
    <connection>
        <map_components component_1="Parameters" component_2="I_SOC"/>
        <map_variables variable_1="p2m" variable_2="p2m"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="P_SOC" component_2="I_SOC"/>
        <map_variables variable_1="P_SOC" variable_2="P_SOC"/>
        <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_SOC"/>
        <map_variables variable_1="E_Ca" variable_2="E_Ca"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_BK"/>
        <map_variables variable_1="E_K" variable_2="E_K"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_NCX"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
    </connection>
    <connection>
        <map_components component_1="I_SOC" component_2="I_NCX"/>
        <map_variables variable_1="SOCPropCa" variable_2="SOCPropCa"/>
    </connection>
    <connection>
        <map_components component_1="Parameters" component_2="I_NCX"/>
        <map_variables variable_1="p2m" variable_2="p2m"/>
        <map_variables variable_1="z_ca" variable_2="z_ca"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_BNa"/>
        <map_variables variable_1="E_Na" variable_2="E_Na"/>
    </connection>
    <connection>
        <map_components component_1="ICC_Membrane" component_2="I_BNa"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="I_BNa" variable_2="I_BNa"/>
    </connection>
    <connection>
        <map_components component_1="Ca_er" component_2="J_IPR"/>
        <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
        <map_variables variable_1="J_IPR" variable_2="J_IPR"/>
    </connection>
    <connection>
        <map_components component_1="J_SERCA" component_2="J_IPR"/>
        <map_variables variable_1="J_SERCA" variable_2="J_SERCA"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="Ca_er"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
        <map_variables variable_1="volER" variable_2="volER"/>
    </connection>
    <connection>
        <map_components component_1="I_SOC" component_2="Ca_er"/>
        <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
    </connection>
    <connection>
        <map_components component_1="Ca_er" component_2="J_SERCA"/>
        <map_variables variable_1="fe" variable_2="fe"/>
        <map_variables variable_1="J_SERCA" variable_2="J_SERCA"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_NSCC"/>
        <map_variables variable_1="E_NSCC" variable_2="E_NSCC"/>
    </connection>
    <connection>
        <map_components component_1="I_NSCC" component_2="P_NSCC"/>
        <map_variables variable_1="P_NSCC" variable_2="P_NSCC"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="Volume" component_2="I_Ano1"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
    </connection>
    <connection>
        <map_components component_1="I_SOC" component_2="I_Ano1"/>
        <map_variables variable_1="J_SOC" variable_2="J_SOC"/>
    </connection>
    <connection>
        <map_components component_1="d_Ano1" component_2="I_Ano1"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="d_Ano1" variable_2="d_Ano1"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="volCyto" variable_2="volCyto"/>
        <map_variables variable_1="J_SOC" variable_2="J_SOC"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="I_Ano1"/>
        <map_variables variable_1="E_Cl" variable_2="E_Cl"/>
    </connection>
    <connection>
        <map_components component_1="y_bind" component_2="J_IPR"/>
        <map_variables variable_1="K1" variable_2="K1"/>
        <map_variables variable_1="K5" variable_2="K5"/>
        <map_variables variable_1="IP3" variable_2="IP3"/>
        <map_variables variable_1="y" variable_2="y"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_Ano1" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_Ano1" variable_2="d_Ano1"/>
        <map_variables variable_1="I_Ano1" variable_2="I_Ano1"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="I_NSCC" component_2="ICC_Membrane"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="I_NSCC" variable_2="I_NSCC"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="Ca_er" component_2="ICC_Membrane"/>
        <map_variables variable_1="Ca_er" variable_2="Ca_er"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="J_SERCA" component_2="ICC_Membrane"/>
        <map_variables variable_1="fc" variable_2="fc"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="J_SERCA" variable_2="J_SERCA"/>
    </connection>
    <connection>
        <map_components component_1="J_IPR" component_2="ICC_Membrane"/>
        <map_variables variable_1="J_IPR" variable_2="J_IPR"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_stim" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="I_stim" variable_2="I_stim"/>
    </connection>
    <connection>
        <map_components component_1="J_PMCA" component_2="ICC_Membrane"/>
        <map_variables variable_1="J_PMCA" variable_2="J_PMCA"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="I_NCX" component_2="ICC_Membrane"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
        <map_variables variable_1="I_NCX" variable_2="I_NCX"/>
        <map_variables variable_1="J_NCX" variable_2="J_NCX"/>
    </connection>
    <connection>
        <map_components component_1="Nernst" component_2="ICC_Membrane"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="I_CaT" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_CaT" variable_2="d_CaT"/>
        <map_variables variable_1="f_CaT" variable_2="f_CaT"/>
        <map_variables variable_1="I_CaT" variable_2="I_CaT"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="J_CaT" variable_2="J_CaT"/>
    </connection>
    <connection>
        <map_components component_1="I_CaV" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_CaV" variable_2="d_CaV"/>
        <map_variables variable_1="I_CaV" variable_2="I_CaV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="J_CaV" variable_2="J_CaV"/>
    </connection>
    <connection>
        <map_components component_1="I_NaV" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_NaV" variable_2="d_NaV"/>
        <map_variables variable_1="f_NaV" variable_2="f_NaV"/>
        <map_variables variable_1="I_NaV" variable_2="I_NaV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
    </connection>
    <connection>
        <map_components component_1="I_KV" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_KV" variable_2="d_KV"/>
        <map_variables variable_1="f_KV" variable_2="f_KV"/>
        <map_variables variable_1="I_KV" variable_2="I_KV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
    </connection>
    <connection>
        <map_components component_1="I_NSV" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="d_NSV" variable_2="d_NSV"/>
        <map_variables variable_1="f_NSV" variable_2="f_NSV"/>
        <map_variables variable_1="I_NSV" variable_2="I_NSV"/>
        <map_variables variable_1="J_NSV" variable_2="J_NSV"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
    </connection>
    <connection>
        <map_components component_1="I_KERG" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="I_KERG" variable_2="I_KERG"/>
        <map_variables variable_1="f_KERG" variable_2="f_KERG"/>
        <map_variables variable_1="d_KERG" variable_2="d_KERG"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
    </connection>
    <connection>
        <map_components component_1="I_SOC" component_2="ICC_Membrane"/>
        <map_variables variable_1="I_SOC" variable_2="I_SOC"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
        <map_variables variable_1="J_SOC" variable_2="J_SOC"/>
    </connection>
    <connection>
        <map_components component_1="I_BK" component_2="ICC_Membrane"/>
        <map_variables variable_1="I_BK" variable_2="I_BK"/>
        <map_variables variable_1="Vm" variable_2="Vm"/>
    </connection>
    <connection>
        <map_components component_1="Time" component_2="ICC_Membrane"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <group>
        <relationship_ref relationship="encapsulation"/>
        <component_ref component="ICC_Membrane">
            <component_ref component="I_CaT">
                <component_ref component="d_CaT"/>
                <component_ref component="f_CaT"/>
            </component_ref>
            <component_ref component="Parameters"/>
            <component_ref component="Volume"/>
            <component_ref component="Nernst"/>
            <component_ref component="I_BK"/>
            <component_ref component="I_BNa"/>
            <component_ref component="J_SERCA"/>
            <component_ref component="J_PMCA"/>
            <component_ref component="I_NCX"/>
            <component_ref component="I_stim"/>
            <component_ref component="J_IPR">
                <component_ref component="y_bind"/>
            </component_ref>
            <component_ref component="Ca_er"/>
            <component_ref component="I_SOC">
                <component_ref component="P_SOC"/>
            </component_ref>
            <component_ref component="I_KERG">
                <component_ref component="d_KERG"/>
                <component_ref component="f_KERG"/>
            </component_ref>
            <component_ref component="I_NaV">
                <component_ref component="d_NaV"/>
                <component_ref component="f_NaV"/>
            </component_ref>
            <component_ref component="I_CaV">
                <component_ref component="d_CaV"/>
            </component_ref>
            <component_ref component="I_NSV">
                <component_ref component="d_NSV"/>
                <component_ref component="f_NSV"/>
            </component_ref>
            <component_ref component="I_KV">
                <component_ref component="d_KV"/>
                <component_ref component="f_KV"/>
            </component_ref>
            <component_ref component="I_NSCC">
                <component_ref component="P_NSCC"/>
            </component_ref>
            <component_ref component="I_Ano1">
                <component_ref component="d_Ano1"/>
            </component_ref>
        </component_ref>
    </group>
</model>