<?xml version='1.0' encoding='UTF-8'?>
<model name="mahajan_shieferaw" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
    <!-- Not sure how to do this in CellML 2.0-->
    <import xlink:href="../units.cellml">
        <units name="u_per_m" units_ref="u_per_m"/>
        <units name="ms" units_ref="ms"/>
        <units name="per_ms" units_ref="per_ms"/>
        <units name="mV" units_ref="mV"/>
        <units name="per_mV" units_ref="per_mV"/>
        <units name="uS_per_nF" units_ref="uS_per_nF"/>
        <units name="nA_per_nF" units_ref="nA_per_nF"/>
        <units name="molar" units_ref="molar"/>
        <units name="mM" units_ref="mM"/>
        <units name="uM" units_ref="uM"/>
        <units name="mM4" units_ref="mM4"/>
        <units name="mV_per_uM" units_ref="mV_per_uM"/>
        <units name="uM_per_ms" units_ref="uM_per_ms"/>
        <units name="per_uM_per_ms" units_ref="per_uM_per_ms"/>
        <units name="J_per_moleK" units_ref="J_per_moleK"/>
        <units name="coulomb_per_mmole" units_ref="coulomb_per_mmole"/>
        <units name="cm_per_s" units_ref="cm_per_s"/>
        <units name="mA_per_cm2" units_ref="mA_per_cm2"/>
        <units name="mmole_per_coulomb_cm" units_ref="mmole_per_coulomb_cm"/>
    </import>
    <import xlink:href="channels/mahajan_2008_LType_calcium_channel.cellml">
        <component component_ref="ICaL" name="ICaL"/>
    </import>
    <import xlink:href="channels/na_channel.cellml">
        <component component_ref="INa" name="INa"/>
    </import>
    <import xlink:href="channels/k1_channel.cellml">
        <component component_ref="IK1" name="IK1"/>
    </import>
    <import xlink:href="channels/kr_channel.cellml">
        <component component_ref="IKr" name="IKr"/>
    </import>
    <import xlink:href="channels/ks_channel.cellml">
        <component component_ref="IKs" name="IKs"/>
    </import>
    <import xlink:href="channels/to_channel.cellml">
        <component component_ref="Ito" name="Ito"/>
    </import>
    <!-- Dummy to use the same concentration component-->
    <component name="IPz1">
        <variable initial_value="0.0" name="xiCaPz1" public_interface="out" units="nA_per_nF"/>
        <variable initial_value="0.0" name="xiNaPz1" public_interface="out" units="nA_per_nF"/>
        <variable initial_value="0.0" name="xiKPz1" public_interface="out" units="nA_per_nF"/>
        <variable initial_value="0.0" name="jCaPz1" public_interface="out" units="uM_per_ms"/>
        <variable initial_value="0.0" name="jNaPz1" public_interface="out" units="uM_per_ms"/>
        <variable initial_value="0.0" name="jKPz1" public_interface="out" units="uM_per_ms"/>
    </component>
    <import xlink:href="pumps/nak_pump.cellml">
        <component component_ref="INaK" name="INaK"/>
    </import>
    <import xlink:href="pumps/naca_pump.cellml">
        <component component_ref="INaCa" name="INaCa"/>
    </import>
    <import xlink:href="concentrations/natrium.cellml">
        <component component_ref="Na" name="Na"/>
    </import>
    <import xlink:href="concentrations/calcium.cellml">
        <component component_ref="Ca" name="Ca"/>
        <component component_ref="Ileak_Iup_Ixfer" name="Ileak_Iup_Ixfer"/>
        <component component_ref="Irel" name="Irel"/>
    </import>
    <component name="interface">
        <variable name="time" private_interface="out" public_interface="in" units="ms"/>
        <variable name="i_stim" private_interface="out" public_interface="in" units="nA_per_nF"/>
        <variable name="pressure" private_interface="out" public_interface="in" units="newton"/>
    </component>
    <group>
        <relationship_ref relationship="encapsulation"/>
        <component_ref component="interface">
            <component_ref component="cell_membrane"/>
            <component_ref component="INa"/>
            <component_ref component="ICaL"/>
            <component_ref component="IK1"/>
            <component_ref component="IKr"/>
            <component_ref component="IKs"/>
            <component_ref component="Ito"/>
            <component_ref component="INaK"/>
            <component_ref component="INaCa"/>
            <component_ref component="Irel"/>
            <component_ref component="Ca"/>
            <component_ref component="Ileak_Iup_Ixfer"/>
            <component_ref component="Na"/>
            <component_ref component="reversal_potentials"/>
            <component_ref component="IPz1"/>
            <!-- How to get rid of this?-->
            <component_ref component="environment"/>
        </component_ref>
    </group>
    <component name="environment">
        <variable initial_value="8.314472" name="R" units="J_per_moleK"/>
        <variable initial_value="308" name="T" units="kelvin"/>
        <variable initial_value="96.4853415" name="F" public_interface="out" units="coulomb_per_mmole"/>
        <variable initial_value="5.4" name="K_o" public_interface="out" units="mM"/>
        <variable initial_value="1.8" name="Ca_o" public_interface="out" units="mM"/>
        <variable initial_value="136" name="Na_o" public_interface="out" units="mM"/>
        <variable name="FonRT" public_interface="out" units="per_mV"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>FonRT</ci>
                <apply>
                    <divide/>
                    <ci>F</ci>
                    <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>T</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="cell_membrane">
        <variable name="time" public_interface="in" units="ms"/>
        <variable initial_value="-87.169816169406" name="V" public_interface="out" units="mV"/>
        <variable name="xik1" public_interface="in" units="nA_per_nF"/>
        <variable name="xito" public_interface="in" units="nA_per_nF"/>
        <variable name="xiNaK" public_interface="in" units="nA_per_nF"/>
        <variable initial_value="8" name="wca" public_interface="out" units="mV_per_uM"/>
        <variable name="xiNaCa" public_interface="in" units="nA_per_nF"/>
        <variable name="xica" public_interface="in" units="nA_per_nF"/>
        <variable name="xina" public_interface="in" units="nA_per_nF"/>
        <variable name="xikr" public_interface="in" units="nA_per_nF"/>
        <variable name="xiks" public_interface="in" units="nA_per_nF"/>
        <variable name="i_stim" public_interface="in" units="nA_per_nF"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>V</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <ci>xina</ci>
                        <ci>xik1</ci>
                        <ci>xikr</ci>
                        <ci>xiks</ci>
                        <ci>xito</ci>
                        <ci>xiNaCa</ci>
                        <ci>xica</ci>
                        <ci>xiNaK</ci>
                        <ci>i_stim</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <component name="reversal_potentials">
        <variable name="ek" public_interface="out" units="mV"/>
        <variable name="eks" public_interface="out" units="mV"/>
        <variable name="ecas" public_interface="out" units="mV"/>
        <variable name="ena" public_interface="out" units="mV"/>
        <variable initial_value="140" name="K_i" units="mM"/>
        <variable name="K_o" public_interface="in" units="mM"/>
        <variable name="Na_i" public_interface="in" units="mM"/>
        <variable name="Na_o" public_interface="in" units="mM"/>
        <variable name="Ca_submem" public_interface="in" units="mM"/>
        <variable name="Ca_o" public_interface="in" units="mM"/>
        <variable name="FonRT" public_interface="in" units="per_mV"/>
        <variable initial_value="0.01833" name="prNaK" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>ek</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>FonRT</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>K_o</ci>
                            <ci>K_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>eks</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>FonRT</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <apply>
                                <plus/>
                                <ci>K_o</ci>
                                <apply>
                                    <times/>
                                    <ci>prNaK</ci>
                                    <ci>Na_o</ci>
                                </apply>
                            </apply>
                            <apply>
                                <plus/>
                                <ci>K_i</ci>
                                <apply>
                                    <times/>
                                    <ci>prNaK</ci>
                                    <ci>Na_i</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ena</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <ci>FonRT</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Na_o</ci>
                            <ci>Na_i</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ecas</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1</cn>
                        <apply>
                            <times/>
                            <ci>FonRT</ci>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>Ca_o</ci>
                            <ci>Ca_submem</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="cell_membrane" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="i_stim" variable_2="i_stim"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="environment"/>
        <map_variables variable_1="K_o" variable_2="K_o"/>
        <map_variables variable_1="Na_o" variable_2="Na_o"/>
        <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
        <map_variables variable_1="FonRT" variable_2="FonRT"/>
    </connection>
    <connection>
        <map_components component_1="INa" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="ICaL" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="ICaL" component_2="environment"/>
        <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
        <map_variables variable_1="F" variable_2="F"/>
        <map_variables variable_1="FonRT" variable_2="FonRT"/>
    </connection>
    <connection>
        <map_components component_1="IKr" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="IKr" component_2="environment"/>
        <map_variables variable_1="K_o" variable_2="K_o"/>
    </connection>
    <connection>
        <map_components component_1="IKs" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="IK1" component_2="environment"/>
        <map_variables variable_1="K_o" variable_2="K_o"/>
    </connection>
    <connection>
        <map_components component_1="Ito" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="INaK" component_2="environment"/>
        <map_variables variable_1="K_o" variable_2="K_o"/>
        <map_variables variable_1="Na_o" variable_2="Na_o"/>
        <map_variables variable_1="FonRT" variable_2="FonRT"/>
    </connection>
    <connection>
        <map_components component_1="INaCa" component_2="environment"/>
        <map_variables variable_1="Ca_o" variable_2="Ca_o"/>
        <map_variables variable_1="Na_o" variable_2="Na_o"/>
        <map_variables variable_1="FonRT" variable_2="FonRT"/>
    </connection>
    <connection>
        <map_components component_1="Irel" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Ca" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Na" component_2="interface"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="INa"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xina" variable_2="xina"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="ICaL"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xica" variable_2="xica"/>
        <map_variables variable_1="wca" variable_2="wca"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="IKr"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xikr" variable_2="xikr"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="IKs"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xiks" variable_2="xiks"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="Ito"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xito" variable_2="xito"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="IK1"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xik1" variable_2="xik1"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="INaK"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xiNaK" variable_2="xiNaK"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="INaCa"/>
        <map_variables variable_1="V" variable_2="V"/>
        <map_variables variable_1="xiNaCa" variable_2="xiNaCa"/>
        <map_variables variable_1="wca" variable_2="wca"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="Irel"/>
        <map_variables variable_1="V" variable_2="V"/>
    </connection>
    <connection>
        <map_components component_1="cell_membrane" component_2="Na"/>
        <map_variables variable_1="wca" variable_2="wca"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="IK1"/>
        <map_variables variable_1="ek" variable_2="ek"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="IKr"/>
        <map_variables variable_1="ek" variable_2="ek"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="Ito"/>
        <map_variables variable_1="ek" variable_2="ek"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="IKs"/>
        <map_variables variable_1="eks" variable_2="eks"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="INa"/>
        <map_variables variable_1="ena" variable_2="ena"/>
    </connection>
    <connection>
        <map_components component_1="reversal_potentials" component_2="Na"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
    </connection>
    <connection>
        <map_components component_1="INa" component_2="Na"/>
        <map_variables variable_1="xina" variable_2="xina"/>
    </connection>
    <connection>
        <map_components component_1="INaK" component_2="Na"/>
        <map_variables variable_1="xiNaK" variable_2="xiNaK"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
    </connection>
    <connection>
        <map_components component_1="IPz1" component_2="Na"/>
        <map_variables variable_1="xiNaPz1" variable_2="xiNaPz1"/>
    </connection>
    <connection>
        <map_components component_1="INaCa" component_2="Na"/>
        <map_variables variable_1="xiNaCa" variable_2="xiNaCa"/>
        <map_variables variable_1="Na_i" variable_2="Na_i"/>
    </connection>
    <connection>
        <map_components component_1="INaCa" component_2="Ca"/>
        <map_variables variable_1="jNaCa" variable_2="jNaCa"/>
        <map_variables variable_1="Ca_submem" variable_2="Ca_submem"/>
        <map_variables variable_1="csm" variable_2="csm"/>
    </connection>
    <connection>
        <map_components component_1="Ca" component_2="reversal_potentials"/>
        <map_variables variable_1="Ca_submem" variable_2="Ca_submem"/>
    </connection>
    <connection>
        <map_components component_1="ICaL" component_2="Ca"/>
        <map_variables variable_1="Ca_dyad" variable_2="Ca_dyad"/>
        <map_variables variable_1="csm" variable_2="csm"/>
        <map_variables variable_1="jca" variable_2="jca"/>
    </connection>
    <connection>
        <map_components component_1="IKs" component_2="Ca"/>
        <map_variables variable_1="Ca_i" variable_2="Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="ICaL" component_2="Irel"/>
        <map_variables variable_1="po" variable_2="po"/>
        <map_variables variable_1="rxa" variable_2="rxa"/>
    </connection>
    <connection>
        <map_components component_1="IPz1" component_2="Ca"/>
        <map_variables variable_1="jCaPz1" variable_2="jCaPz1"/>
    </connection>
</model>
