Rendering of the source text

<?xml version='1.0'?>
<model cmeta:id="Tong_Choi_Kharche_Holden_Zhang_Taggart_2011" name="Tong_Choi_Kharche_Holden_Zhang_Taggart_2011" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#" xmlns:cmeta="http://www.cellml.org/metadata/1.0#">
    <!-- units-->
    <units name="msec">
        <unit prefix="milli" units="second"/>
    </units>
    <units name="joule_per_kelvin_per_kilomole">
        <unit units="joule"/>
        <unit exponent="-1" units="kelvin"/>
        <unit exponent="-1" prefix="kilo" units="mole"/>
    </units>
    <units name="mM">
        <unit prefix="milli" units="mole"/>
        <unit exponent="-1" units="liter"/>
    </units>
    <units name="mM4">
        <unit exponent="4" units="mM"/>
    </units>
    <units name="mM_per_msec">
        <unit units="mM"/>
        <unit exponent="-1" prefix="milli" units="second"/>
    </units>
    <units name="M_per_msec">
        <unit units="mole"/>
        <unit exponent="-1" prefix="milli" units="second"/>
    </units>
    <units name="coulomb_per_mole">
        <unit units="coulomb"/>
        <unit exponent="-1" units="mole"/>
    </units>
    <units name="uA_per_uF">
        <unit prefix="micro" units="ampere"/>
        <unit exponent="-1" prefix="micro" units="farad"/>
    </units>
    <units name="per_msec">
        <unit exponent="-1" prefix="milli" units="second"/>
    </units>
    <units name="mV">
        <unit prefix="milli" units="volt"/>
    </units>
    <units name="per_mV">
        <unit exponent="-1" units="mV"/>
    </units>
    <units name="uF_per_cm2">
        <unit prefix="micro" units="farad"/>
        <unit exponent="-2" prefix="centi" units="meter"/>
    </units>
    <units name="cm2_per_uL">
        <unit exponent="2" prefix="centi" units="meter"/>
        <unit exponent="-1" prefix="micro" units="liter"/>
    </units>
    <units name="nS_per_pF">
        <unit prefix="nano" units="siemens"/>
        <unit exponent="-1" prefix="pico" units="farad"/>
    </units>
    <units name="pA_per_pF">
        <unit prefix="pico" units="ampere"/>
        <unit exponent="-1" prefix="pico" units="farad"/>
    </units>
    <units name="nM">
        <unit prefix="nano" units="mole"/>
        <unit exponent="-1" units="liter"/>
    </units>
    <units name="uN">
        <unit prefix="micro" units="newton"/>
    </units>
    <units name="nM_to_mM">
        <unit units="mM"/>
        <unit exponent="-1" units="nM"/>
    </units>
    <units name="mM_to_M">
        <unit units="mole"/>
        <unit exponent="-1" units="mM"/>
    </units>
    <!-- Model - stimulation interface-->
    <component name="interface">
        <variable initial_value="0" name="I_hold" public_interface="none" units="pA_per_pF"/>
        <!--below default-->
        <variable initial_value="-0.5" name="I_test" public_interface="none" units="pA_per_pF"/>
        <!--below for test comparison Tong, 2014 Fig 8
var I_test: pA_per_pF {init: -0.2, pub: none};
below for test comparison Tong, 2014 Fig 11-->
        <!--var I_test: pA_per_pF {init: -0.9, pub: none};-->
        <variable initial_value="1000" name="test_start" public_interface="none" units="msec"/>
        <!--below default-->
        <variable initial_value="3000" name="test_end" public_interface="none" units="msec"/>
        <!--below for test comparison Tong, 2014 Fig 8
var test_end: msec {init: 50000, pub: none};
below for test comparison Tong, 2014 Fig 11
var test_end: msec {init: 10000, pub: none};-->
        <variable name="Ist" private_interface="out" units="pA_per_pF"/>
        <variable name="time" private_interface="in" units="msec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Ist</ci>
                <piecewise>
                    <piece>
                        <ci>I_test</ci>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <ci>test_start</ci>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <ci>test_end</ci>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <ci>I_hold</ci>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="interface" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="interface" component_2="membrane_potential"/>
        <map_variables variable_1="Ist" variable_2="Ist"/>
    </connection>
    <!-- Model - Time variable-->
    <component name="environment">
        <variable initial_value="0" name="time" public_interface="out" units="msec"/>
    </component>
    <!-- Model - V variable-->
    <component name="membrane_potential">
        <variable name="I_tot" public_interface="none" units="pA_per_pF"/>
        <variable name="I_Ca_tot" public_interface="out" units="pA_per_pF"/>
        <variable initial_value="-53.90915441282156" name="v" public_interface="out" units="mV"/>
        <variable name="Ist" public_interface="in" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="Cm" public_interface="in" units="uF_per_cm2"/>
        <variable name="ina" public_interface="in" units="pA_per_pF"/>
        <variable name="ical" public_interface="in" units="pA_per_pF"/>
        <variable name="icat" public_interface="in" units="pA_per_pF"/>
        <variable name="ib" public_interface="in" units="pA_per_pF"/>
        <variable name="ik1" public_interface="in" units="pA_per_pF"/>
        <!--below without ik2
var ik2: pA_per_pF {pub: in};
below without ika
var ika: pA_per_pF {pub: in};-->
        <variable name="iBKa" public_interface="in" units="pA_per_pF"/>
        <variable name="iBKab" public_interface="in" units="pA_per_pF"/>
        <!--var ih: pA_per_pF {pub: in};-->
        <variable name="icl" public_interface="in" units="pA_per_pF"/>
        <variable name="insna" public_interface="in" units="pA_per_pF"/>
        <variable name="insca" public_interface="in" units="pA_per_pF"/>
        <variable name="insk" public_interface="in" units="pA_per_pF"/>
        <!--below without inak
var inak: pA_per_pF {pub: in};
below without ncx current
var inaca: pA_per_pF {pub: in};-->
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>I_Ca_tot</ci>
                <apply>
                    <plus/>
                    <ci>ical</ci>
                    <ci>icat</ci>
                    <ci>insca</ci>
                </apply>
            </apply>
            <!--below original total with all currents
I_tot = ina+ih+inaca+inak+ical+icat+icl+ik1+ik2+ika+iBKa+iBKab+insna+insk+insca+ib;
below without 'hyperpol' current ih
I_tot = ina+inaca+inak+ical+icat+icl+ik1+ik2+ika+iBKa+iBKab+insna+insk+insca+ib;
below without K2 current ik2
I_tot = ina+inaca+inak+ical+icat+icl+ik1+ika+iBKa+iBKab+insna+insk+insca+ib;
below without Ka current ika
I_tot = ina+inaca+inak+ical+icat+icl+ik1+iBKa+iBKab+insna+insk+insca+ib;
below without N/K current inak
I_tot = ina+inaca+ical+icat+icl+ik1+iBKa+iBKab+insna+insk+insca+ib;
below without NCX current inaca; note jnaca still computed...-->
            <apply>
                <eq/>
                <ci>I_tot</ci>
                <apply>
                    <plus/>
                    <ci>ina</ci>
                    <ci>ical</ci>
                    <ci>icat</ci>
                    <ci>icl</ci>
                    <ci>ik1</ci>
                    <ci>iBKa</ci>
                    <ci>iBKab</ci>
                    <ci>insna</ci>
                    <ci>insk</ci>
                    <ci>insca</ci>
                    <ci>ib</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>v</ci>
                </apply>
                <apply>
                    <minus/>
                    <apply>
                        <plus/>
                        <ci>I_tot</ci>
                        <ci>Ist</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="membrane_potential" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="membrane_potential" component_2="parameters"/>
        <map_variables variable_1="Cm" variable_2="Cm"/>
        <map_variables variable_1="v" variable_2="v"/>
    </connection>
    <!-- Model - Cai variable-->
    <component name="Ca_Concentrations">
        <variable name="J_tot" public_interface="none" units="mM_per_msec"/>
        <variable name="J_Ca_mem" public_interface="none" units="mM_per_msec"/>
        <variable initial_value="0.0001161881607214449" name="cai" public_interface="out" units="mM"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="jnaca" public_interface="in" units="mM_per_msec"/>
        <variable name="jpmca" public_interface="in" units="mM_per_msec"/>
        <variable name="I_Ca_tot" public_interface="in" units="pA_per_pF"/>
        <variable name="buff" public_interface="in" units="dimensionless"/>
        <variable name="AV" public_interface="in" units="cm2_per_uL"/>
        <variable name="zca" public_interface="in" units="dimensionless"/>
        <variable name="frdy" public_interface="in" units="coulomb_per_mole"/>
        <variable name="Cm" public_interface="in" units="uF_per_cm2"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>J_Ca_mem</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>AV</ci>
                            <ci>Cm</ci>
                            <ci>buff</ci>
                        </apply>
                        <apply>
                            <times/>
                            <ci>zca</ci>
                            <ci>frdy</ci>
                        </apply>
                    </apply>
                    <ci>I_Ca_tot</ci>
                </apply>
            </apply>
        </math>
        <variable name="jcamem_plot" public_interface="none" units="M_per_msec"/>
        <variable name="jpmca_plot" public_interface="none" units="M_per_msec"/>
        <variable name="jnaca_plot" public_interface="none" units="M_per_msec"/>
        <variable initial_value="1000" name="conversion" public_interface="none" units="mM_to_M"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>jcamem_plot</ci>
                <apply>
                    <times/>
                    <ci>J_Ca_mem</ci>
                    <ci>conversion</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>jpmca_plot</ci>
                <apply>
                    <times/>
                    <ci>jpmca</ci>
                    <ci>conversion</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>jnaca_plot</ci>
                <apply>
                    <times/>
                    <ci>jnaca</ci>
                    <ci>conversion</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>J_tot</ci>
                <apply>
                    <plus/>
                    <ci>J_Ca_mem</ci>
                    <ci>jnaca</ci>
                    <ci>jpmca</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>cai</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>J_tot</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="Ca_Concentrations" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Ca_Concentrations" component_2="membrane_potential"/>
        <map_variables variable_1="I_Ca_tot" variable_2="I_Ca_tot"/>
    </connection>
    <connection>
        <map_components component_1="Ca_Concentrations" component_2="parameters"/>
        <map_variables variable_1="cai" variable_2="cai"/>
        <map_variables variable_1="buff" variable_2="buff"/>
        <map_variables variable_1="AV" variable_2="AV"/>
        <map_variables variable_1="zca" variable_2="zca"/>
        <map_variables variable_1="frdy" variable_2="frdy"/>
        <map_variables variable_1="Cm" variable_2="Cm"/>
    </connection>
    <!-- Model - Parameters-->
    <component name="parameters">
        <variable name="ki" public_interface="out" units="mM"/>
        <variable name="nai" public_interface="out" units="mM"/>
        <variable name="cli" public_interface="out" units="mM"/>
        <variable name="ko" public_interface="out" units="mM"/>
        <variable name="cao" public_interface="out" units="mM"/>
        <variable name="nao" public_interface="out" units="mM"/>
        <variable name="clo" public_interface="out" units="mM"/>
        <variable name="mgo" public_interface="out" units="mM"/>
        <variable name="buff" public_interface="out" units="dimensionless"/>
        <variable name="AV" public_interface="out" units="cm2_per_uL"/>
        <variable name="zca" public_interface="out" units="dimensionless"/>
        <variable name="zna" public_interface="out" units="dimensionless"/>
        <variable name="zk" public_interface="out" units="dimensionless"/>
        <variable name="R" public_interface="out" units="joule_per_kelvin_per_kilomole"/>
        <variable name="frdy" public_interface="out" units="coulomb_per_mole"/>
        <variable name="temp" public_interface="out" units="kelvin"/>
        <variable name="Cm" public_interface="out" units="uF_per_cm2"/>
        <variable name="gna" public_interface="out" units="nS_per_pF"/>
        <variable name="gcal" public_interface="out" units="nS_per_pF"/>
        <variable name="ecal" public_interface="out" units="mV"/>
        <variable name="kmca" public_interface="out" units="mM"/>
        <variable name="gcat" public_interface="out" units="nS_per_pF"/>
        <variable name="ecat" public_interface="out" units="mV"/>
        <variable name="gkca" public_interface="out" units="nS_per_pF"/>
        <variable name="gb" public_interface="out" units="nS_per_pF"/>
        <variable name="gk1" public_interface="out" units="nS_per_pF"/>
        <!--below removed for ik2
var gk2: nS_per_pF {pub: out};-->
        <variable name="gbka" public_interface="out" units="dimensionless"/>
        <variable name="gbkab" public_interface="out" units="dimensionless"/>
        <!--below removed for ika
var gka: nS_per_pF {pub: out};-->
        <variable name="gcl" public_interface="out" units="nS_per_pF"/>
        <!--below removed for ih
var gh: nS_per_pF {pub: out};-->
        <variable name="gns" public_interface="out" units="nS_per_pF"/>
        <variable name="PnsK" public_interface="out" units="dimensionless"/>
        <variable name="PnsNa" public_interface="out" units="dimensionless"/>
        <variable name="PnsCa" public_interface="out" units="dimensionless"/>
        <variable name="PnsCs" public_interface="out" units="dimensionless"/>
        <variable name="gnsCa" public_interface="out" units="dimensionless"/>
        <variable name="gnsNa" public_interface="out" units="dimensionless"/>
        <variable name="gnsK" public_interface="out" units="dimensionless"/>
        <variable name="gnsCs" public_interface="out" units="dimensionless"/>
        <!--below removed for inak-->
        <variable name="PK" public_interface="out" units="dimensionless"/>
        <variable name="PNa" public_interface="out" units="dimensionless"/>
        <variable name="Jpmca" public_interface="out" units="mM_per_msec"/>
        <variable name="Kmpmca" public_interface="out" units="mM"/>
        <variable name="npmca" public_interface="out" units="dimensionless"/>
        <variable name="Jnaca" public_interface="out" units="mM_per_msec"/>
        <variable name="Kmallo" public_interface="out" units="mM"/>
        <variable name="nallo" public_interface="out" units="dimensionless"/>
        <variable name="ksat" public_interface="out" units="dimensionless"/>
        <variable name="xgamma" public_interface="out" units="dimensionless"/>
        <variable name="Kmnai" public_interface="out" units="mM"/>
        <variable name="Kmcai" public_interface="out" units="mM"/>
        <variable name="Kmnao" public_interface="out" units="mM"/>
        <variable name="Kmcao" public_interface="out" units="mM"/>
        <variable name="Fmax" public_interface="out" units="uN"/>
        <variable name="FKm" public_interface="out" units="nM"/>
        <variable name="Fn" public_interface="out" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>ki</ci>
                <cn cellml:units="mM">140.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>nai</ci>
                <cn cellml:units="mM">4.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>cli</ci>
                <cn cellml:units="mM">46.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>ko</ci>
                <cn cellml:units="mM">6.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>cao</ci>
                <cn cellml:units="mM">2.5</cn>
            </apply>
            <apply>
                <eq/>
                <ci>nao</ci>
                <cn cellml:units="mM">130.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>clo</ci>
                <cn cellml:units="mM">130.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>mgo</ci>
                <cn cellml:units="mM">0.5</cn>
            </apply>
            <!--below base value from cellml
buff = 0.015{dimensionless};-->
            <!--below value from spreadsheet test
'Final parms gK 0.7 gNa 1.4x gCaL 1.0x with buff scaling 1.05x beta'
buff = 0.0158{dimensionless};-->
            <!--above 1.05x not enough Ca cyt
below try 1.1x beta
buff = 0.0165{dimensionless};-->
            <!--even higher to 1.125x beta-->
            <apply>
                <eq/>
                <ci>buff</ci>
                <cn cellml:units="dimensionless">0.0169</cn>
            </apply>
            <apply>
                <eq/>
                <ci>AV</ci>
                <cn cellml:units="cm2_per_uL">4.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>zca</ci>
                <cn cellml:units="dimensionless">2.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>zna</ci>
                <cn cellml:units="dimensionless">1.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>zk</ci>
                <cn cellml:units="dimensionless">1.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>R</ci>
                <cn cellml:units="joule_per_kelvin_per_kilomole">8314.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>frdy</ci>
                <cn cellml:units="coulomb_per_mole">96485.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>temp</ci>
                <cn cellml:units="kelvin">308.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Cm</ci>
                <cn cellml:units="uF_per_cm2">1.0</cn>
            </apply>
            <!--below commented out; base value given in cellml
gna = 0{nS_per_pF};-->
            <!--below from spreadsheet test
'Final parms gK 0.7 gNa 1.4x gCaL 1.0x with buff scaling 1.05x beta'
here, gNa used from usmc_model_v3p9 , g_na = mean([0.0028 0.125]) = 0.0639
=> 1.4x0.0639 = 0.0895-->
            <apply>
                <eq/>
                <ci>gna</ci>
                <cn cellml:units="nS_per_pF">0.0895</cn>
            </apply>
            <!--below for test comparison with Tong, 2014 Fig 11
gna = 0.125{nS_per_pF};-->
            <!--below base value given in cellml-->
            <apply>
                <eq/>
                <ci>gcal</ci>
                <cn cellml:units="nS_per_pF">0.6</cn>
            </apply>
            <!--below for test comparison with Tong, 2014 Fig 11
gcal = 1.5{nS_per_pF};-->
            <apply>
                <eq/>
                <ci>ecal</ci>
                <cn cellml:units="mV">45.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>kmca</ci>
                <cn cellml:units="mM">0.001</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gcat</ci>
                <cn cellml:units="nS_per_pF">0.058</cn>
            </apply>
            <apply>
                <eq/>
                <ci>ecat</ci>
                <cn cellml:units="mV">42.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gkca</ci>
                <cn cellml:units="nS_per_pF">0.8</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gb</ci>
                <cn cellml:units="nS_per_pF">0.004</cn>
            </apply>
            <!--below default opencor val
gk1 = 0.52{nS_per_pF};
below from usmc_model_v3p9 
'comparison with full tong reproduce their fig 12'-->
            <!--gk1 = 0.6315{nS_per_pF};-->
            <!--below from spreadsheet test
'Final parms gK 0.7 gNa 1.4x gCaL 1.0x with buff scaling 1.05x beta'-->
            <!--gk1 = 0.7{nS_per_pF};-->
            <!--below for test comparison with Tong, 2014 Fig 8-->
            <!--gk1 = 0.75{nS_per_pF};-->
            <!--below for test comparison with Tong, 2014 Fig 11
note overall maximum K current is about 0.8...-->
            <apply>
                <eq/>
                <ci>gk1</ci>
                <cn cellml:units="nS_per_pF">0.7</cn>
            </apply>
            <!--gk1 = 0.85{nS_per_pF};-->
            <!--below removed for ik2
gk2 = 0.032{nS_per_pF};-->
            <apply>
                <eq/>
                <ci>gbka</ci>
                <cn cellml:units="dimensionless">0.2</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gbkab</ci>
                <cn cellml:units="dimensionless">0.1</cn>
            </apply>
            <!--below removed for ika
gka = 0.16{nS_per_pF};-->
            <apply>
                <eq/>
                <ci>gcl</ci>
                <cn cellml:units="nS_per_pF">0.1875</cn>
            </apply>
            <!--below removed for ih
gh = 0.0542{nS_per_pF};-->
            <apply>
                <eq/>
                <ci>gns</ci>
                <cn cellml:units="nS_per_pF">0.0123</cn>
            </apply>
            <apply>
                <eq/>
                <ci>PnsK</ci>
                <cn cellml:units="dimensionless">1.3</cn>
            </apply>
            <apply>
                <eq/>
                <ci>PnsNa</ci>
                <cn cellml:units="dimensionless">0.9</cn>
            </apply>
            <apply>
                <eq/>
                <ci>PnsCa</ci>
                <cn cellml:units="dimensionless">0.89</cn>
            </apply>
            <apply>
                <eq/>
                <ci>PnsCs</ci>
                <cn cellml:units="dimensionless">1.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gnsCa</ci>
                <cn cellml:units="dimensionless">0.5</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gnsNa</ci>
                <cn cellml:units="dimensionless">1.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gnsK</ci>
                <cn cellml:units="dimensionless">1.19</cn>
            </apply>
            <apply>
                <eq/>
                <ci>gnsCs</ci>
                <cn cellml:units="dimensionless">1.6</cn>
            </apply>
            <!--below without inak-->
            <apply>
                <eq/>
                <ci>PK</ci>
                <cn cellml:units="dimensionless">1.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>PNa</ci>
                <cn cellml:units="dimensionless">0.35</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Jpmca</ci>
                <cn cellml:units="mM_per_msec">0.00000035</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmpmca</ci>
                <cn cellml:units="mM">0.0005</cn>
            </apply>
            <apply>
                <eq/>
                <ci>npmca</ci>
                <cn cellml:units="dimensionless">2</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Jnaca</ci>
                <cn cellml:units="mM_per_msec">0.0000035</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmallo</ci>
                <cn cellml:units="mM">0.0003</cn>
            </apply>
            <apply>
                <eq/>
                <ci>nallo</ci>
                <cn cellml:units="dimensionless">4</cn>
            </apply>
            <apply>
                <eq/>
                <ci>ksat</ci>
                <cn cellml:units="dimensionless">0.27</cn>
            </apply>
            <apply>
                <eq/>
                <ci>xgamma</ci>
                <cn cellml:units="dimensionless">0.35</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmnai</ci>
                <cn cellml:units="mM">30.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmcai</ci>
                <cn cellml:units="mM">0.007</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmnao</ci>
                <cn cellml:units="mM">87.5</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Kmcao</ci>
                <cn cellml:units="mM">1.3</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Fmax</ci>
                <cn cellml:units="uN">3.0</cn>
            </apply>
            <apply>
                <eq/>
                <ci>FKm</ci>
                <cn cellml:units="nM">161.301</cn>
            </apply>
            <apply>
                <eq/>
                <ci>Fn</ci>
                <cn cellml:units="dimensionless">3.60205</cn>
            </apply>
        </math>
        <variable name="vFRT" public_interface="out" units="dimensionless"/>
        <variable name="ena" public_interface="out" units="mV"/>
        <variable name="ek" public_interface="out" units="mV"/>
        <!--below removed for ih 
var eh: mV {pub: out};-->
        <variable name="ecl" public_interface="out" units="mV"/>
        <variable name="enscc" public_interface="out" units="mV"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>vFRT</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>v</ci>
                        <ci>frdy</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>R</ci>
                        <ci>temp</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ena</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>temp</ci>
                        </apply>
                        <ci>frdy</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>nao</ci>
                            <ci>nai</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ek</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>temp</ci>
                        </apply>
                        <ci>frdy</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>ko</ci>
                            <ci>ki</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below removed for ih
eh = R*temp/frdy*ln((ko+PNa/PK*nao)/(ki+PNa/PK*nai));-->
            <apply>
                <eq/>
                <ci>ecl</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>temp</ci>
                        </apply>
                        <ci>frdy</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <ci>cli</ci>
                            <ci>clo</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>enscc</ci>
                <apply>
                    <times/>
                    <apply>
                        <divide/>
                        <apply>
                            <times/>
                            <ci>R</ci>
                            <ci>temp</ci>
                        </apply>
                        <ci>frdy</ci>
                    </apply>
                    <apply>
                        <ln/>
                        <apply>
                            <divide/>
                            <apply>
                                <plus/>
                                <apply>
                                    <times/>
                                    <ci>PnsK</ci>
                                    <ci>ko</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>PnsNa</ci>
                                    <ci>nao</ci>
                                </apply>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <cn cellml:units="dimensionless">4</cn>
                                        <ci>PnsCa</ci>
                                        <ci>cao</ci>
                                    </apply>
                                    <apply>
                                        <plus/>
                                        <cn cellml:units="dimensionless">1.0</cn>
                                        <apply>
                                            <exp/>
                                            <ci>vFRT</ci>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                            <apply>
                                <plus/>
                                <apply>
                                    <times/>
                                    <ci>PnsK</ci>
                                    <ci>ki</ci>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>PnsNa</ci>
                                    <ci>nai</ci>
                                </apply>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <times/>
                                        <cn cellml:units="dimensionless">4</cn>
                                        <ci>PnsCa</ci>
                                        <ci>cai</ci>
                                    </apply>
                                    <apply>
                                        <plus/>
                                        <cn cellml:units="dimensionless">1.0</cn>
                                        <apply>
                                            <exp/>
                                            <ci>vFRT</ci>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <!-- Model - Force and w variable-->
    <component name="Ca_dependent_Force">
        <variable name="wss" public_interface="none" units="dimensionless"/>
        <variable name="wtc" public_interface="none" units="msec"/>
        <variable initial_value="1e-6" name="conversion" public_interface="none" units="nM_to_mM"/>
        <variable name="Force" public_interface="none" units="uN"/>
        <variable initial_value="0.2345238135343783" name="w" public_interface="none" units="dimensionless"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="Fmax" public_interface="in" units="uN"/>
        <variable name="FKm" public_interface="in" units="nM"/>
        <variable name="Fn" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>wss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <ci>FKm</ci>
                                    <ci>conversion</ci>
                                </apply>
                                <ci>cai</ci>
                            </apply>
                            <ci>Fn</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>wtc</ci>
                <apply>
                    <times/>
                    <cn cellml:units="msec">4000.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">0.234845</cn>
                        <apply>
                            <divide/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless">1.0</cn>
                                <cn cellml:units="dimensionless">0.234845</cn>
                            </apply>
                            <apply>
                                <plus/>
                                <cn cellml:units="dimensionless">1.0</cn>
                                <apply>
                                    <power/>
                                    <apply>
                                        <divide/>
                                        <ci>cai</ci>
                                        <apply>
                                            <times/>
                                            <ci>FKm</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                    </apply>
                                    <ci>Fn</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>Force</ci>
                <apply>
                    <times/>
                    <ci>Fmax</ci>
                    <apply>
                        <minus/>
                        <ci>w</ci>
                        <cn cellml:units="dimensionless">0.2345</cn>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>w</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>wss</ci>
                        <ci>w</ci>
                    </apply>
                    <ci>wtc</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="Ca_dependent_Force" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="Ca_dependent_Force" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
    </connection>
    <connection>
        <map_components component_1="Ca_dependent_Force" component_2="parameters"/>
        <map_variables variable_1="Fmax" variable_2="Fmax"/>
        <map_variables variable_1="FKm" variable_2="FKm"/>
        <map_variables variable_1="Fn" variable_2="Fn"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_Na">
        <variable name="mss" public_interface="none" units="dimensionless"/>
        <variable name="hss" public_interface="none" units="dimensionless"/>
        <!--below commented out for ss-approximation
var mtc: msec {pub: none};-->
        <variable name="htc" public_interface="none" units="msec"/>
        <!--below commented out for ss-approximation
var m: dimensionless {init: 0.1253518889572223, pub: none};-->
        <variable initial_value="0.404599170710196" name="h" public_interface="none" units="dimensionless"/>
        <variable name="ina" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="ena" public_interface="in" units="mV"/>
        <variable name="gna" public_interface="in" units="nS_per_pF"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>mss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">35.9584</cn>
                                    </apply>
                                </apply>
                                <cn cellml:units="mV">9.24013</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>hss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>v</ci>
                                    <cn cellml:units="mV">57.0</cn>
                                </apply>
                                <cn cellml:units="mV">8.0</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below commented out for ss-approximation
mtc = 0.25{msec}+7.0{msec}/(1.0{dimensionless}+exp((v+38.0{mV})/10.0{mV}));-->
            <apply>
                <eq/>
                <ci>htc</ci>
                <apply>
                    <plus/>
                    <cn cellml:units="msec">0.9</cn>
                    <apply>
                        <divide/>
                        <cn cellml:units="msec">1002.85</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">47.5</cn>
                                    </apply>
                                    <cn cellml:units="mV">1.5</cn>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below for ss-approximation calc-->
            <!--ina = gna*mss*mss*mss*h*(v-ena);-->
            <!--below with both m and h ss-->
            <apply>
                <eq/>
                <ci>ina</ci>
                <apply>
                    <times/>
                    <ci>gna</ci>
                    <ci>mss</ci>
                    <ci>mss</ci>
                    <ci>mss</ci>
                    <ci>hss</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ena</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--below commented out for ss-approximation-->
        <!--ode(h, time) = (hss-h)/htc;-->
    </component>
    <connection>
        <map_components component_1="I_Na" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_Na" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="ina" variable_2="ina"/>
    </connection>
    <connection>
        <map_components component_1="I_Na" component_2="parameters"/>
        <map_variables variable_1="ena" variable_2="ena"/>
        <map_variables variable_1="gna" variable_2="gna"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_CaL">
        <variable name="dss" public_interface="none" units="dimensionless"/>
        <variable name="fss" public_interface="none" units="dimensionless"/>
        <variable name="fca" public_interface="none" units="dimensionless"/>
        <!--below commented out for ss approximation-->
        <variable name="f2tc" public_interface="none" units="msec"/>
        <!--below commented out for ss approximation-->
        <variable initial_value="0.9065967263076083" name="f2" public_interface="none" units="dimensionless"/>
        <variable name="ical" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="ecal" public_interface="in" units="mV"/>
        <variable name="gcal" public_interface="in" units="nS_per_pF"/>
        <variable name="kmca" public_interface="in" units="mM"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>fca</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>cai</ci>
                                <ci>kmca</ci>
                            </apply>
                            <cn cellml:units="dimensionless">4.0</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>dss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">22.0</cn>
                                    </apply>
                                </apply>
                                <cn cellml:units="mV">7.0</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>fss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>v</ci>
                                    <cn cellml:units="mV">38.0</cn>
                                </apply>
                                <cn cellml:units="mV">7.0</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below commented out for ss approximation-->
            <apply>
                <eq/>
                <ci>f2tc</ci>
                <apply>
                    <times/>
                    <cn cellml:units="msec">90.9699</cn>
                    <apply>
                        <minus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <times/>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="dimensionless">1.0</cn>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <plus/>
                                                <ci>v</ci>
                                                <cn cellml:units="mV">13.9629</cn>
                                            </apply>
                                            <cn cellml:units="mV">45.3782</cn>
                                        </apply>
                                    </apply>
                                </apply>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="dimensionless">1.0</cn>
                                    <apply>
                                        <exp/>
                                        <apply>
                                            <divide/>
                                            <apply>
                                                <minus/>
                                                <apply>
                                                    <plus/>
                                                    <ci>v</ci>
                                                    <cn cellml:units="mV">9.49866</cn>
                                                </apply>
                                            </apply>
                                            <cn cellml:units="mV">3.3945</cn>
                                        </apply>
                                    </apply>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below for ss-calculation-->
            <apply>
                <eq/>
                <ci>ical</ci>
                <apply>
                    <times/>
                    <ci>gcal</ci>
                    <ci>fca</ci>
                    <ci>dss</ci>
                    <ci>dss</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0.8</cn>
                            <ci>fss</ci>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0.2</cn>
                            <ci>f2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ecal</ci>
                    </apply>
                </apply>
            </apply>
            <!--below commented out for ss approximation-->
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>f2</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>fss</ci>
                        <ci>f2</ci>
                    </apply>
                    <ci>f2tc</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="I_CaL" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_CaL" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="ical" variable_2="ical"/>
    </connection>
    <connection>
        <map_components component_1="I_CaL" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
    </connection>
    <connection>
        <map_components component_1="I_CaL" component_2="parameters"/>
        <map_variables variable_1="ecal" variable_2="ecal"/>
        <map_variables variable_1="gcal" variable_2="gcal"/>
        <map_variables variable_1="kmca" variable_2="kmca"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_CaT">
        <variable name="bss" public_interface="none" units="dimensionless"/>
        <variable name="gss" public_interface="none" units="dimensionless"/>
        <!--below commented out for ss approximation-->
        <variable name="icat" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="ecat" public_interface="in" units="mV"/>
        <variable name="gcat" public_interface="in" units="nS_per_pF"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>bss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">54.23</cn>
                                    </apply>
                                </apply>
                                <cn cellml:units="mV">9.88</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gss</ci>
                <apply>
                    <plus/>
                    <cn cellml:units="dimensionless">0.02</cn>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">0.98</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">72.978</cn>
                                    </apply>
                                    <cn cellml:units="mV">4.64</cn>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below ss-approximation-->
            <apply>
                <eq/>
                <ci>icat</ci>
                <apply>
                    <times/>
                    <ci>gcat</ci>
                    <ci>bss</ci>
                    <ci>bss</ci>
                    <ci>gss</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ecat</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--below commented out for ss approximation-->
    </component>
    <connection>
        <map_components component_1="I_CaT" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_CaT" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="icat" variable_2="icat"/>
    </connection>
    <connection>
        <map_components component_1="I_CaT" component_2="parameters"/>
        <map_variables variable_1="ecat" variable_2="ecat"/>
        <map_variables variable_1="gcat" variable_2="gcat"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_b">
        <variable name="ib" public_interface="out" units="pA_per_pF"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="ek" public_interface="in" units="mV"/>
        <variable name="gb" public_interface="in" units="nS_per_pF"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>ib</ci>
                <apply>
                    <times/>
                    <ci>gb</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ek</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="I_b" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="ib" variable_2="ib"/>
    </connection>
    <connection>
        <map_components component_1="I_b" component_2="parameters"/>
        <map_variables variable_1="ek" variable_2="ek"/>
        <map_variables variable_1="gb" variable_2="gb"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_K1">
        <variable name="qss" public_interface="none" units="dimensionless"/>
        <variable name="rss" public_interface="none" units="dimensionless"/>
        <variable name="qtc" public_interface="none" units="msec"/>
        <variable name="r1tc" public_interface="none" units="msec"/>
        <variable name="r2tc" public_interface="none" units="msec"/>
        <variable initial_value="0.2060363247740295" name="q" public_interface="none" units="dimensionless"/>
        <variable initial_value="0.1922244113609531" name="r1" public_interface="none" units="dimensionless"/>
        <variable initial_value="0.1932803618375963" name="r2" public_interface="none" units="dimensionless"/>
        <variable name="ik1" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="ek" public_interface="in" units="mV"/>
        <variable name="gk1" public_interface="in" units="nS_per_pF"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>qss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">0.978613</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <minus/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">18.6736</cn>
                                    </apply>
                                </apply>
                                <cn cellml:units="mV">26.6606</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>rss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>v</ci>
                                    <cn cellml:units="mV">63.0</cn>
                                </apply>
                                <cn cellml:units="mV">6.3</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>qtc</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="msec">500.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>v</ci>
                                    <cn cellml:units="mV">60.71</cn>
                                </apply>
                                <cn cellml:units="mV">15.79</cn>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>r1tc</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="msec">5000.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <apply>
                                    <plus/>
                                    <ci>v</ci>
                                    <cn cellml:units="mV">62.7133</cn>
                                </apply>
                                <cn cellml:units="mV">35.8611</cn>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>r2tc</ci>
                <apply>
                    <plus/>
                    <cn cellml:units="msec">30000.0</cn>
                    <apply>
                        <divide/>
                        <cn cellml:units="msec">220000.0</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">22.0</cn>
                                    </apply>
                                    <cn cellml:units="mV">4.0</cn>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ik1</ci>
                <apply>
                    <times/>
                    <ci>gk1</ci>
                    <ci>q</ci>
                    <ci>q</ci>
                    <apply>
                        <plus/>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0.38</cn>
                            <ci>r1</ci>
                        </apply>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0.63</cn>
                            <ci>r2</ci>
                        </apply>
                    </apply>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ek</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>qss</ci>
                        <ci>q</ci>
                    </apply>
                    <ci>qtc</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>r1</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>rss</ci>
                        <ci>r1</ci>
                    </apply>
                    <ci>r1tc</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>r2</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>rss</ci>
                        <ci>r2</ci>
                    </apply>
                    <ci>r2tc</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="I_K1" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_K1" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="ik1" variable_2="ik1"/>
    </connection>
    <connection>
        <map_components component_1="I_K1" component_2="parameters"/>
        <map_variables variable_1="ek" variable_2="ek"/>
        <map_variables variable_1="gk1" variable_2="gk1"/>
    </connection>
    <!-- Model - currents-->
    <!-- Model - currents-->
    <!-- Model - currents-->
    <component name="I_BKa">
        <variable name="xass_z" public_interface="none" units="dimensionless"/>
        <variable name="xass_vh" public_interface="none" units="mV"/>
        <variable initial_value="1000.0" name="conversion" public_interface="none" units="mM_to_M"/>
        <variable name="xass" public_interface="none" units="dimensionless"/>
        <!--below commented out for ss approximation-->
        <variable name="iBKa" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="ek" public_interface="in" units="mV"/>
        <variable name="gkca" public_interface="in" units="nS_per_pF"/>
        <variable name="gbka" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="joule_per_kelvin_per_kilomole"/>
        <variable name="frdy" public_interface="in" units="coulomb_per_mole"/>
        <variable name="temp" public_interface="in" units="kelvin"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>xass_z</ci>
                <apply>
                    <plus/>
                    <apply>
                        <divide/>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.749234</cn>
                        </apply>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <minus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">0.0630535</cn>
                                    </apply>
                                    <cn cellml:units="mole">0.161942</cn>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">8.38384</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">1538.29</cn>
                                    </apply>
                                    <cn cellml:units="mole">739.057</cn>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>xass_vh</ci>
                <apply>
                    <minus/>
                    <apply>
                        <divide/>
                        <cn cellml:units="mV">5011.47</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">0.237503</cn>
                                    </apply>
                                    <cn cellml:units="mole">0.000239278</cn>
                                </apply>
                                <cn cellml:units="dimensionless">0.42291</cn>
                            </apply>
                        </apply>
                    </apply>
                    <cn cellml:units="mV">37.5137</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>xass</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <ci>xass_z</ci>
                                    </apply>
                                    <ci>frdy</ci>
                                    <apply>
                                        <minus/>
                                        <ci>v</ci>
                                        <ci>xass_vh</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>temp</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below for ss computation-->
            <apply>
                <eq/>
                <ci>iBKa</ci>
                <apply>
                    <times/>
                    <ci>gkca</ci>
                    <ci>gbka</ci>
                    <ci>xass</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ek</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--below commented out for ss approximation-->
    </component>
    <connection>
        <map_components component_1="I_BKa" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_BKa" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="iBKa" variable_2="iBKa"/>
    </connection>
    <connection>
        <map_components component_1="I_BKa" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
    </connection>
    <connection>
        <map_components component_1="I_BKa" component_2="parameters"/>
        <map_variables variable_1="ek" variable_2="ek"/>
        <map_variables variable_1="gkca" variable_2="gkca"/>
        <map_variables variable_1="gbka" variable_2="gbka"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="frdy" variable_2="frdy"/>
        <map_variables variable_1="temp" variable_2="temp"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_BKab">
        <variable name="xabss_z" public_interface="none" units="dimensionless"/>
        <variable name="xabss_vh" public_interface="none" units="mV"/>
        <variable initial_value="1000.0" name="conversion" public_interface="none" units="mM_to_M"/>
        <variable name="xabss" public_interface="none" units="dimensionless"/>
        <!--below commented out for ss approximation-->
        <variable name="iBKab" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="ek" public_interface="in" units="mV"/>
        <variable name="gkca" public_interface="in" units="nS_per_pF"/>
        <variable name="gbkab" public_interface="in" units="dimensionless"/>
        <variable name="R" public_interface="in" units="joule_per_kelvin_per_kilomole"/>
        <variable name="frdy" public_interface="in" units="coulomb_per_mole"/>
        <variable name="temp" public_interface="in" units="kelvin"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>xabss_z</ci>
                <apply>
                    <plus/>
                    <apply>
                        <divide/>
                        <apply>
                            <minus/>
                            <cn cellml:units="dimensionless">0.681249</cn>
                        </apply>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <minus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">0.218988</cn>
                                    </apply>
                                    <cn cellml:units="mole">0.428335</cn>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">1.40001</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">228.71</cn>
                                    </apply>
                                    <cn cellml:units="mole">684.946</cn>
                                </apply>
                                <cn cellml:units="dimensionless">2</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>xabss_vh</ci>
                <apply>
                    <minus/>
                    <apply>
                        <divide/>
                        <cn cellml:units="mV">8540.23</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <apply>
                                            <times/>
                                            <ci>cai</ci>
                                            <ci>conversion</ci>
                                        </apply>
                                        <cn cellml:units="mole">0.401189</cn>
                                    </apply>
                                    <cn cellml:units="mole">0.00399115</cn>
                                </apply>
                                <cn cellml:units="dimensionless">0.668054</cn>
                            </apply>
                        </apply>
                    </apply>
                    <cn cellml:units="mV">109.275</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>xabss</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <exp/>
                            <apply>
                                <divide/>
                                <apply>
                                    <times/>
                                    <apply>
                                        <minus/>
                                        <ci>xabss_z</ci>
                                    </apply>
                                    <ci>frdy</ci>
                                    <apply>
                                        <minus/>
                                        <ci>v</ci>
                                        <ci>xabss_vh</ci>
                                    </apply>
                                </apply>
                                <apply>
                                    <times/>
                                    <ci>R</ci>
                                    <ci>temp</ci>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <!--below for ss computation-->
            <apply>
                <eq/>
                <ci>iBKab</ci>
                <apply>
                    <times/>
                    <ci>gkca</ci>
                    <ci>gbkab</ci>
                    <ci>xabss</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ek</ci>
                    </apply>
                </apply>
            </apply>
        </math>
        <!--below commented out for ss approximation-->
    </component>
    <connection>
        <map_components component_1="I_BKab" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_BKab" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="iBKab" variable_2="iBKab"/>
    </connection>
    <connection>
        <map_components component_1="I_BKab" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
    </connection>
    <connection>
        <map_components component_1="I_BKab" component_2="parameters"/>
        <map_variables variable_1="ek" variable_2="ek"/>
        <map_variables variable_1="gkca" variable_2="gkca"/>
        <map_variables variable_1="gbkab" variable_2="gbkab"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="frdy" variable_2="frdy"/>
        <map_variables variable_1="temp" variable_2="temp"/>
    </connection>
    <!-- Model - currents-->
    <!-- Model - currents-->
    <component name="I_Cl">
        <variable name="css" public_interface="none" units="dimensionless"/>
        <variable name="ctc" public_interface="none" units="msec"/>
        <variable name="K1cl" public_interface="none" units="mM"/>
        <variable name="K2cl" public_interface="none" units="dimensionless"/>
        <variable initial_value="0.0003764413740731269" name="c" public_interface="none" units="dimensionless"/>
        <variable name="icl" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="ecl" public_interface="in" units="mV"/>
        <variable name="gcl" public_interface="in" units="nS_per_pF"/>
        <variable name="vFRT" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>K1cl</ci>
                <apply>
                    <times/>
                    <cn cellml:units="mM">0.0006</cn>
                    <apply>
                        <exp/>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">2.53</cn>
                            <ci>vFRT</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>K2cl</ci>
                <apply>
                    <times/>
                    <cn cellml:units="dimensionless">0.1</cn>
                    <apply>
                        <exp/>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless">5.0</cn>
                            </apply>
                            <ci>vFRT</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>css</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <times/>
                            <ci>K2cl</ci>
                            <apply>
                                <plus/>
                                <apply>
                                    <power/>
                                    <apply>
                                        <divide/>
                                        <ci>K1cl</ci>
                                        <ci>cai</ci>
                                    </apply>
                                    <cn cellml:units="dimensionless">2</cn>
                                </apply>
                                <apply>
                                    <divide/>
                                    <ci>K1cl</ci>
                                    <ci>cai</ci>
                                </apply>
                                <cn cellml:units="dimensionless">1.0</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>ctc</ci>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <cn cellml:units="msec">160.0</cn>
                    </apply>
                    <apply>
                        <divide/>
                        <cn cellml:units="msec">210.0</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <plus/>
                                        <ci>v</ci>
                                        <cn cellml:units="mV">4.56</cn>
                                    </apply>
                                    <cn cellml:units="mV">11.62</cn>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                    <apply>
                        <divide/>
                        <cn cellml:units="msec">170.0</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <exp/>
                                <apply>
                                    <divide/>
                                    <apply>
                                        <minus/>
                                        <apply>
                                            <plus/>
                                            <ci>v</ci>
                                            <cn cellml:units="mV">25.5</cn>
                                        </apply>
                                    </apply>
                                    <cn cellml:units="mV">11.62</cn>
                                </apply>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>icl</ci>
                <apply>
                    <times/>
                    <ci>gcl</ci>
                    <ci>c</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>ecl</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>c</ci>
                </apply>
                <apply>
                    <divide/>
                    <apply>
                        <minus/>
                        <ci>css</ci>
                        <ci>c</ci>
                    </apply>
                    <ci>ctc</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="I_Cl" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_Cl" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="icl" variable_2="icl"/>
    </connection>
    <connection>
        <map_components component_1="I_Cl" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
    </connection>
    <connection>
        <map_components component_1="I_Cl" component_2="parameters"/>
        <map_variables variable_1="ecl" variable_2="ecl"/>
        <map_variables variable_1="gcl" variable_2="gcl"/>
        <map_variables variable_1="vFRT" variable_2="vFRT"/>
    </connection>
    <!-- Model - currents-->
    <component name="I_ns">
        <variable name="fmg" public_interface="none" units="dimensionless"/>
        <variable name="gs_nao" public_interface="none" units="dimensionless"/>
        <variable name="gs_cao" public_interface="none" units="dimensionless"/>
        <variable name="gs_ko" public_interface="none" units="dimensionless"/>
        <variable initial_value="1e-8" name="tinyamount" public_interface="none" units="mM"/>
        <variable name="insna" public_interface="out" units="pA_per_pF"/>
        <variable name="insca" public_interface="out" units="pA_per_pF"/>
        <variable name="insk" public_interface="out" units="pA_per_pF"/>
        <variable name="time" public_interface="in" units="msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="ko" public_interface="in" units="mM"/>
        <variable name="cao" public_interface="in" units="mM"/>
        <variable name="nao" public_interface="in" units="mM"/>
        <variable name="mgo" public_interface="in" units="mM"/>
        <variable name="enscc" public_interface="in" units="mV"/>
        <variable name="gns" public_interface="in" units="nS_per_pF"/>
        <variable name="gnsCa" public_interface="in" units="dimensionless"/>
        <variable name="gnsNa" public_interface="in" units="dimensionless"/>
        <variable name="gnsK" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>fmg</ci>
                <apply>
                    <plus/>
                    <cn cellml:units="dimensionless">0.108043</cn>
                    <apply>
                        <divide/>
                        <cn cellml:units="dimensionless">0.903902</cn>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <ci>mgo</ci>
                                    <cn cellml:units="mM">0.281007</cn>
                                </apply>
                                <cn cellml:units="dimensionless">1.29834</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gs_nao</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <cn cellml:units="dimensionless">0.0123</cn>
                        </apply>
                        <cn cellml:units="dimensionless">0.03</cn>
                    </apply>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <cn cellml:units="mM">150.0</cn>
                                <apply>
                                    <plus/>
                                    <ci>nao</ci>
                                    <ci>tinyamount</ci>
                                </apply>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gs_ko</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <cn cellml:units="dimensionless">0.0123</cn>
                        </apply>
                        <cn cellml:units="dimensionless">0.03</cn>
                    </apply>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <cn cellml:units="mM">150.0</cn>
                                <apply>
                                    <plus/>
                                    <ci>ko</ci>
                                    <ci>tinyamount</ci>
                                </apply>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>gs_cao</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <apply>
                            <divide/>
                            <cn cellml:units="dimensionless">1</cn>
                            <cn cellml:units="dimensionless">0.000525</cn>
                        </apply>
                        <cn cellml:units="dimensionless">0.03</cn>
                    </apply>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <cn cellml:units="mM">150.0</cn>
                                <apply>
                                    <plus/>
                                    <ci>cao</ci>
                                    <ci>tinyamount</ci>
                                </apply>
                            </apply>
                            <cn cellml:units="dimensionless">2</cn>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>insca</ci>
                <apply>
                    <times/>
                    <ci>fmg</ci>
                    <ci>gs_cao</ci>
                    <ci>gnsCa</ci>
                    <ci>gns</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>enscc</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>insna</ci>
                <apply>
                    <times/>
                    <ci>fmg</ci>
                    <ci>gs_nao</ci>
                    <ci>gnsNa</ci>
                    <ci>gns</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>enscc</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>insk</ci>
                <apply>
                    <times/>
                    <ci>fmg</ci>
                    <ci>gs_ko</ci>
                    <ci>gnsK</ci>
                    <ci>gns</ci>
                    <apply>
                        <minus/>
                        <ci>v</ci>
                        <ci>enscc</ci>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="I_ns" component_2="environment"/>
        <map_variables variable_1="time" variable_2="time"/>
    </connection>
    <connection>
        <map_components component_1="I_ns" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
        <map_variables variable_1="insna" variable_2="insna"/>
        <map_variables variable_1="insca" variable_2="insca"/>
        <map_variables variable_1="insk" variable_2="insk"/>
    </connection>
    <connection>
        <map_components component_1="I_ns" component_2="parameters"/>
        <map_variables variable_1="enscc" variable_2="enscc"/>
        <map_variables variable_1="gns" variable_2="gns"/>
        <map_variables variable_1="gnsNa" variable_2="gnsNa"/>
        <map_variables variable_1="gnsCa" variable_2="gnsCa"/>
        <map_variables variable_1="gnsK" variable_2="gnsK"/>
        <map_variables variable_1="nao" variable_2="nao"/>
        <map_variables variable_1="cao" variable_2="cao"/>
        <map_variables variable_1="ko" variable_2="ko"/>
        <map_variables variable_1="mgo" variable_2="mgo"/>
    </connection>
    <!-- Model - currents-->
    <!-- Model - currents-->
    <!-- Model - flux-->
    <component name="J_NaCa">
        <variable name="f1naca" public_interface="none" units="dimensionless"/>
        <variable name="f2naca" public_interface="none" units="dimensionless"/>
        <variable name="fallo" public_interface="none" units="dimensionless"/>
        <variable name="naca_Eup" public_interface="none" units="mM4"/>
        <variable name="naca_Ed1" public_interface="none" units="dimensionless"/>
        <variable name="naca_Ed2" public_interface="none" units="mM4"/>
        <variable name="naca_Ed3" public_interface="none" units="mM4"/>
        <variable initial_value="-1" name="jnaca_sign" public_interface="none" units="dimensionless"/>
        <variable name="jnaca" public_interface="out" units="mM_per_msec"/>
        <variable name="v" public_interface="in" units="mV"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="cao" public_interface="in" units="mM"/>
        <variable name="nao" public_interface="in" units="mM"/>
        <variable name="nai" public_interface="in" units="mM"/>
        <variable name="Jnaca" public_interface="in" units="mM_per_msec"/>
        <variable name="Kmallo" public_interface="in" units="mM"/>
        <variable name="nallo" public_interface="in" units="dimensionless"/>
        <variable name="ksat" public_interface="in" units="dimensionless"/>
        <variable name="xgamma" public_interface="in" units="dimensionless"/>
        <variable name="Kmnai" public_interface="in" units="mM"/>
        <variable name="Kmcai" public_interface="in" units="mM"/>
        <variable name="Kmnao" public_interface="in" units="mM"/>
        <variable name="Kmcao" public_interface="in" units="mM"/>
        <variable name="vFRT" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>f1naca</ci>
                <apply>
                    <exp/>
                    <apply>
                        <times/>
                        <apply>
                            <minus/>
                            <ci>xgamma</ci>
                            <cn cellml:units="dimensionless">1.0</cn>
                        </apply>
                        <ci>vFRT</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>f2naca</ci>
                <apply>
                    <exp/>
                    <apply>
                        <times/>
                        <ci>xgamma</ci>
                        <ci>vFRT</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>fallo</ci>
                <apply>
                    <divide/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>Kmallo</ci>
                                <ci>cai</ci>
                            </apply>
                            <ci>nallo</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>naca_Eup</ci>
                <apply>
                    <minus/>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>nai</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>cao</ci>
                        <ci>f2naca</ci>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>nao</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>cai</ci>
                        <ci>f1naca</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>naca_Ed1</ci>
                <apply>
                    <plus/>
                    <cn cellml:units="dimensionless">1.0</cn>
                    <apply>
                        <times/>
                        <ci>ksat</ci>
                        <ci>f1naca</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>naca_Ed2</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <ci>Kmcao</ci>
                        <apply>
                            <power/>
                            <ci>nai</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>Kmnao</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>cai</ci>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>Kmnai</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>cao</ci>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <divide/>
                                <ci>cai</ci>
                                <ci>Kmcai</ci>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>naca_Ed3</ci>
                <apply>
                    <plus/>
                    <apply>
                        <times/>
                        <ci>cao</ci>
                        <apply>
                            <power/>
                            <ci>nai</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>nao</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>cai</ci>
                    </apply>
                    <apply>
                        <times/>
                        <apply>
                            <power/>
                            <ci>nao</ci>
                            <cn cellml:units="dimensionless">3</cn>
                        </apply>
                        <ci>Kmcai</ci>
                        <apply>
                            <plus/>
                            <cn cellml:units="dimensionless">1.0</cn>
                            <apply>
                                <power/>
                                <apply>
                                    <divide/>
                                    <ci>nai</ci>
                                    <ci>Kmnai</ci>
                                </apply>
                                <cn cellml:units="dimensionless">3</cn>
                            </apply>
                        </apply>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>jnaca</ci>
                <apply>
                    <divide/>
                    <apply>
                        <times/>
                        <ci>jnaca_sign</ci>
                        <ci>Jnaca</ci>
                        <ci>fallo</ci>
                        <ci>naca_Eup</ci>
                    </apply>
                    <apply>
                        <times/>
                        <ci>naca_Ed1</ci>
                        <apply>
                            <plus/>
                            <ci>naca_Ed2</ci>
                            <ci>naca_Ed3</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="J_NaCa" component_2="membrane_potential"/>
        <map_variables variable_1="v" variable_2="v"/>
    </connection>
    <connection>
        <map_components component_1="J_NaCa" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
        <map_variables variable_1="jnaca" variable_2="jnaca"/>
    </connection>
    <connection>
        <map_components component_1="J_NaCa" component_2="parameters"/>
        <map_variables variable_1="nao" variable_2="nao"/>
        <map_variables variable_1="nai" variable_2="nai"/>
        <map_variables variable_1="cao" variable_2="cao"/>
        <map_variables variable_1="Jnaca" variable_2="Jnaca"/>
        <map_variables variable_1="Kmallo" variable_2="Kmallo"/>
        <map_variables variable_1="nallo" variable_2="nallo"/>
        <map_variables variable_1="ksat" variable_2="ksat"/>
        <map_variables variable_1="xgamma" variable_2="xgamma"/>
        <map_variables variable_1="Kmnai" variable_2="Kmnai"/>
        <map_variables variable_1="Kmcai" variable_2="Kmcai"/>
        <map_variables variable_1="Kmnao" variable_2="Kmnao"/>
        <map_variables variable_1="Kmcao" variable_2="Kmcao"/>
        <map_variables variable_1="vFRT" variable_2="vFRT"/>
    </connection>
    <!-- Model - flux-->
    <component name="J_PMCA">
        <variable name="jpmca" public_interface="out" units="mM_per_msec"/>
        <variable name="cai" public_interface="in" units="mM"/>
        <variable name="Jpmca" public_interface="in" units="mM_per_msec"/>
        <variable name="Kmpmca" public_interface="in" units="mM"/>
        <variable name="npmca" public_interface="in" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>jpmca</ci>
                <apply>
                    <divide/>
                    <ci>Jpmca</ci>
                    <apply>
                        <plus/>
                        <cn cellml:units="dimensionless">1.0</cn>
                        <apply>
                            <power/>
                            <apply>
                                <divide/>
                                <ci>Kmpmca</ci>
                                <ci>cai</ci>
                            </apply>
                            <ci>npmca</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="J_PMCA" component_2="Ca_Concentrations"/>
        <map_variables variable_1="cai" variable_2="cai"/>
        <map_variables variable_1="jpmca" variable_2="jpmca"/>
    </connection>
    <connection>
        <map_components component_1="J_PMCA" component_2="parameters"/>
        <map_variables variable_1="Jpmca" variable_2="Jpmca"/>
        <map_variables variable_1="Kmpmca" variable_2="Kmpmca"/>
        <map_variables variable_1="npmca" variable_2="npmca"/>
    </connection>
    <!-- Model - Group-->
    <group>
        <relationship_ref relationship="encapsulation"/>
        <component_ref component="interface">
            <component_ref component="environment"/>
            <component_ref component="parameters"/>
            <component_ref component="membrane_potential"/>
            <component_ref component="Ca_Concentrations"/>
            <component_ref component="Ca_dependent_Force"/>
            <component_ref component="I_Na"/>
            <component_ref component="I_CaL"/>
            <component_ref component="I_CaT"/>
            <component_ref component="I_b"/>
            <component_ref component="I_K1"/>
            <!--below removed for reduced
comp I_K2;
below removed for reduced
comp I_Ka;-->
            <component_ref component="I_BKa"/>
            <component_ref component="I_BKab"/>
            <!--below removed for reduced
comp I_h;-->
            <component_ref component="I_Cl"/>
            <component_ref component="I_ns"/>
            <!--below removed for inak
comp I_NaK;
below removed for ncx current inaca
comp I_NaCa;-->
            <component_ref component="J_NaCa"/>
            <component_ref component="J_PMCA"/>
        </component_ref>
    </group>
    <!-- Model - Mapping variables-->
    <documentation xmlns="http://cellml.org/tmp-documentation">
        <article xmlns="http://cellml.org/tmp-documentation">
            <articleinfo xmlns="http://cellml.org/tmp-documentation">
                <title xmlns="http://cellml.org/tmp-documentation">A computational model of the ionic currents, Ca2+ dynamics and action potentials underlying contraction of isolated uterine smooth muscle</title>
                <author xmlns="http://cellml.org/tmp-documentation">
                    <firstname xmlns="http://cellml.org/tmp-documentation">Wing Chiu</firstname>
                    <surname xmlns="http://cellml.org/tmp-documentation">Tong</surname>
                    <affiliation xmlns="http://cellml.org/tmp-documentation">
                        <shortaffil xmlns="http://cellml.org/tmp-documentation">Newcastle University</shortaffil>
                    </affiliation>
                </author>
            </articleinfo>
            <section xmlns="http://cellml.org/tmp-documentation" id="sec_status">
                <title xmlns="http://cellml.org/tmp-documentation">Model Status</title>
                <para xmlns="http://cellml.org/tmp-documentation">
           This CellML model is based on a translation of the original C code. The equations are validated using OpenCell and all units are consistent. The model can be run in both OpenCell (version 0.8) and OpenCOR (version 0.4) environments using any of the integration methods. To recreate the published results, for example figure 12 from the original paper, apply a suitable stimulating current and integrate the model for a total 10 sec with a time step of 0.02 ms. 
          </para>
            </section>
            <sect1 xmlns="http://cellml.org/tmp-documentation" id="sec_structure">
                <title xmlns="http://cellml.org/tmp-documentation">Model Structure</title>
                <para xmlns="http://cellml.org/tmp-documentation">ABSTRACT:  Uterine contractions during labor are discretely regulated by rhythmic action potentials (AP) of varying duration and form that serve to determine Ca2+ -dependent force production. We have employed a computational biology approach to develop a fuller understanding of the complexity of excitation-contraction coupling of uterine smooth muscle cells (USMC). Our overall aim is to establish a mathematical platform of sufficient biophysical detail to quantitatively describe known uterine E-C coupling parameters and thereby inform future empirical investigations of physiological and pathophysiological mechanisms governing normal and dysfunctional labors. From published and unpublished data we construct mathematical models for thirteen ionic currents of USMCs: Ca2+ currents (L- and T-type), Na+ current, and hyperpolarization-activated current, three voltage-gated K+ currents, Ca2+ -activated K+ current,Ca2+ -activated Cl- current, non-specific cation current, Na+ -Ca2+ exchanger, Na+ -K+ pump and background current. The magnitudes and kinetics of each current system in a spindle shaped single cell with a specified surface area:volume ratio is described by differential equations, in terms of maximal conductances, electrochemical gradient, voltage-dependent activation/inactivation gating variables and temporal changes in intracellular Ca2+ computed from known Ca2+ fluxes. These quantifications are validated by the reconstruction of the individual experimental ionic currents obtained under voltage-clamp. Phasic contraction is modelled in relation to the time constant of changing [Ca2+ ]i . This integrated model is validated by its reconstruction of the different action potential configurations of USMCs (spikes,plateau and bursts of spikes), the change from bursting to plateau type AP produced by estradiol, and of simultaneous experimental recordings of spontaneous AP, [Ca2+ ]i and phasic force. In summary, our advanced mathematical model provides a powerful tool to investigate the physiological ionic mechanisms underlying the genesis of uterine electrical excitation-contraction coupling of labor and parturition. This will furnish the evolution of descriptive and predictive quantitative models of myometrial electrogenesis at the whole cell and tissue levels.
</para>
                <para xmlns="http://cellml.org/tmp-documentation">
The original paper reference is cited below:
</para>
                <para xmlns="http://cellml.org/tmp-documentation">
A computational model of the ionic currents, Ca2+ dynamics and action potentials underlying contraction of isolated uterine smooth muscle. Wing-Chiu Tong, Cecilia Y. Choi, Sanjay Kharche, Arun V. Holden, Henggui Zhang and Michael J. Taggart, 2011, PLoS One, 6(4), e18685.  <ulink xmlns="http://cellml.org/tmp-documentation" url="http://www.ncbi.nlm.nih.gov/pubmed/21559514">PubMed ID: 21559514</ulink>
                </para>
                <informalfigure xmlns="http://cellml.org/tmp-documentation" float="0" id="fig_cell_diagram_faber">
                    <mediaobject xmlns="http://cellml.org/tmp-documentation">
                        <imageobject xmlns="http://cellml.org/tmp-documentation">
                            <objectinfo xmlns="http://cellml.org/tmp-documentation">
                                <title xmlns="http://cellml.org/tmp-documentation">cell diagram of the USMC models showing the ionic currents, pumps and exchangers within the sarcolemma</title>
                            </objectinfo>
                            <imagedata xmlns="http://cellml.org/tmp-documentation" fileref="USMC_Tong_et_al_2011.png"/>
                        </imageobject>
                    </mediaobject>
                    <caption xmlns="http://cellml.org/tmp-documentation">A schematic diagram of the Tong et. al. 2011 uterine smooth muscle cell model. </caption>
                </informalfigure>
            </sect1>
        </article>
    </documentation>
    <rdf:RDF xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="">
            <title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">A computational model of the ionic currents, Ca2+ dynamics and action potentials underlying contraction of isolated uterine smooth muscle</title>
            <dc:creator xmlns:dc="http://purl.org/dc/elements/1.1/" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                    <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Tong</vCard:Family>
                    <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Wing Chiu</vCard:Given>
                    <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
                </vCard:N>
                <vCard:EMAIL xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                    <rdf:value xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">winnie.tong@ncl.ac.uk</rdf:value>
                    <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:resource="http://imc.org/vCard/3.0#internet"/>
                </vCard:EMAIL>
                <vCard:ORG xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                    <vCard:Orgname xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Newcastle University</vCard:Orgname>
                    <vCard:Orgunit xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Institute of Cellular Medicine</vCard:Orgunit>
                </vCard:ORG>
            </dc:creator>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Tong_Choi_Kharche_Holden_Zhang_Taggart_2011">
            <reference xmlns="http://www.cellml.org/bqs/1.0#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Subject_citation"/>
            </reference>
            <reference xmlns="http://www.cellml.org/bqs/1.0#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Subject_Model_keywords"/>
            </reference>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Subject_citation">
            <Pubmed_id xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">21559514</Pubmed_id>
            <JournalArticle xmlns="http://www.cellml.org/bqs/1.0#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#citation"/>
            </JournalArticle>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#citation">
            <title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">A computational model of the ionic currents, Ca2+ dynamics and action potentials underlying contraction of isolated uterine smooth muscle</title>
            <creator xmlns="http://purl.org/dc/elements/1.1/">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Authors"/>
            </creator>
            <volume xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">6</volume>
            <first_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">e18685</first_page>
            <Journal xmlns="http://www.cellml.org/bqs/1.0#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Journal"/>
            </Journal>
            <last_page xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace"/>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Journal">
            <title xmlns="http://purl.org/dc/elements/1.1/" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">PLoS One</title>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Authors">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Seq"/>
            </rdf:type>
            <rdf:_1 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_1"/>
            </rdf:_1>
            <rdf:_2 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_2"/>
            </rdf:_2>
            <rdf:_3 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_3"/>
            </rdf:_3>
            <rdf:_4 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_4"/>
            </rdf:_4>
            <rdf:_5 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_5"/>
            </rdf:_5>
            <rdf:_6 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_6"/>
            </rdf:_6>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_1">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Tong</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">WC</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_2">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Choi</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">C</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Y</vCard:Other>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_3">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Kharche</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">S</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_4">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Holden</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">A</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">V</vCard:Other>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_5">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Zhang</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">H</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#"/>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Author_6">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.cellml.org/bqs/1.0#Person"/>
            </rdf:type>
            <vCard:N xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#" rdf:parseType="Resource" xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <vCard:Family xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">Taggart</vCard:Family>
                <vCard:Given xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">M</vCard:Given>
                <vCard:Other xmlns:vCard="http://www.w3.org/2001/vcard-rdf/3.0#">J</vCard:Other>
            </vCard:N>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Subject_Model_keywords">
            <subject xmlns="http://purl.org/dc/elements/1.1/">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Model_keywords"/>
            </subject>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#Model_keywords">
            <rdf:value xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#keywords"/>
            </rdf:value>
            <subject_type xmlns="http://www.cellml.org/bqs/1.0#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">keyword</subject_type>
        </rdf:Description>
        <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="#keywords">
            <rdf:type xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#">
                <rdf:Description xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" rdf:about="http://www.w3.org/1999/02/22-rdf-syntax-ns#Bag"/>
            </rdf:type>
            <rdf:_1 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">electrophysiology</rdf:_1>
            <rdf:_2 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">excitation-contraction_coupling</rdf:_2>
            <rdf:_3 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">calcium</rdf:_3>
            <rdf:_4 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">force</rdf:_4>
            <rdf:_5 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">pregnancy</rdf:_5>
            <rdf:_6 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">uterus</rdf:_6>
            <rdf:_7 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">myometrium</rdf:_7>
            <rdf:_8 xmlns:rdf="http://www.w3.org/1999/02/22-rdf-syntax-ns#" xml:lang="en" xmlns:xml="http://www.w3.org/XML/1998/namespace">smooth muscle</rdf:_8>
        </rdf:Description>
    </rdf:RDF>
</model>