- Author:
- Leyla <noroozbabaee@gmail.com>
- Date:
- 2022-05-10 14:01:08+12:00
- Desc:
- Adding Tong_2011 to PMR
- Permanent Source URI:
- https://models.cellml.org/workspace/85c/rawfile/a03f680a69226c515cd789723c8036028406852b/Experiments/Figure_1/Fig1_Tong_2011.cellml
<?xml version='1.0' encoding='UTF-8'?>
<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.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" 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-->
<!-- def comp interface as
var I_hold: pA_per_pF {init: 0, pub: none};
var I_test: pA_per_pF {init: 0.0, pub: none};//origin -0.5
var test_start: msec {init: 1000, pub: none};
var test_end: msec {init: 3000, pub: none};
var Ist: pA_per_pF {priv: out};
var time: msec {priv: in};-->
<!-- Ist = sel
case (time > test_start) and (time < test_end):
I_test;
otherwise:
I_hold;
endsel;
enddef;-->
<!-- def map between interface and environment for
vars time and time;
enddef;-->
<!-- def map between interface and membrane_potential for
vars Ist and Ist;
enddef;-->
<!-- Model - Time variable-->
<component name="environment">
<variable initial_value="-130" name="vv" public_interface="out" units="mV"/>
</component>
<!-- Model - V variable-->
<!-- def comp membrane_potential as
var I_tot: pA_per_pF {pub: none};
var I_Ca_tot: pA_per_pF {pub: out};
var v: mV {init: -153.99, pub: out}; //origin:-53.90915441282156, table:-53.999
var Ist: pA_per_pF {pub: in};
var time: msec {pub: in};
var Cm: uF_per_cm2 {pub: in};
var ina: pA_per_pF {pub: in};
var ical: pA_per_pF {pub: in};
var icat: pA_per_pF {pub: in};
var ib: pA_per_pF {pub: in};
var ik1: pA_per_pF {pub: in};
var ik2: pA_per_pF {pub: in};
var ika: pA_per_pF {pub: in};
var iBKa: pA_per_pF {pub: in};
var iBKab: pA_per_pF {pub: in};
var ih: pA_per_pF {pub: in};
var icl: pA_per_pF {pub: in};
var insna: pA_per_pF {pub: in};
var insca: pA_per_pF {pub: in};
var insk: pA_per_pF {pub: in};
var inak: pA_per_pF {pub: in};
var inaca: pA_per_pF {pub: in};-->
<!-- I_Ca_tot = ical+icat+insca;
I_tot = ina+ih+inaca+inak+ical+icat+icl+ik1+ik2+ika+iBKa+iBKab+insna+insk+insca+ib;
ode(v, time) = -(I_tot+Ist);
enddef;-->
<!-- def map between membrane_potential and environment for
vars time and time;
enddef;-->
<!-- def map between membrane_potential and parameters for
vars Cm and Cm;
vars v and v;
enddef;-->
<!-- // Model - Cai variable-->
<!-- def comp Ca_Concentrations as
var J_tot: mM_per_msec {pub: none};
var J_Ca_mem: mM_per_msec {pub: none};
var cai: mM {init: 0.0001161881607214449, pub: out};
var time: msec {pub: in};
var jnaca: mM_per_msec {pub: in};
var jpmca: mM_per_msec {pub: in};
var I_Ca_tot: pA_per_pF {pub: in};
var buff: dimensionless {pub: in};
var AV: cm2_per_uL {pub: in};
var zca: dimensionless {pub: in};
var frdy: coulomb_per_mole {pub: in};
var Cm: uF_per_cm2 {pub: in};-->
<!-- J_Ca_mem = AV*Cm*buff/(zca*frdy)*I_Ca_tot;-->
<!-- var jcamem_plot: M_per_msec {pub: none};
var jpmca_plot: M_per_msec {pub: none};
var jnaca_plot: M_per_msec {pub: none};
var conversion: mM_to_M {init: 1000, pub: none};-->
<!-- jcamem_plot = J_Ca_mem*conversion;
jpmca_plot = jpmca*conversion;
jnaca_plot = jnaca*conversion;
J_tot = J_Ca_mem+jnaca+jpmca;
ode(cai, time) = -J_tot;
enddef;-->
<!-- def map between Ca_Concentrations and environment for
vars time and time;
enddef;-->
<!-- def map between Ca_Concentrations and membrane_potential for
vars I_Ca_tot and I_Ca_tot;
enddef;-->
<!-- def map between Ca_Concentrations and parameters for
vars cai and cai;
vars buff and buff;
vars AV and AV;
vars zca and zca;
vars frdy and frdy;
vars Cm and Cm;
enddef;-->
<!-- // 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"/>
<variable name="gk2" public_interface="out" units="nS_per_pF"/>
<variable name="gbka" public_interface="out" units="dimensionless"/>
<variable name="gbkab" public_interface="out" units="dimensionless"/>
<variable name="gka" public_interface="out" units="nS_per_pF"/>
<variable name="gcl" public_interface="out" units="nS_per_pF"/>
<variable name="gh" public_interface="out" units="nS_per_pF"/>
<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"/>
<variable name="ginak" public_interface="out" units="pA_per_pF"/>
<variable name="nakKmko" public_interface="out" units="mM"/>
<variable name="nakKmnai" public_interface="out" units="mM"/>
<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>
<apply>
<eq/>
<ci>buff</ci>
<cn cellml:units="dimensionless">0.015</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>
<apply>
<eq/>
<ci>gna</ci>
<cn cellml:units="nS_per_pF">0</cn>
</apply>
<apply>
<eq/>
<ci>gcal</ci>
<cn cellml:units="nS_per_pF">0.6</cn>
</apply>
<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>
<apply>
<eq/>
<ci>gk1</ci>
<cn cellml:units="nS_per_pF">0.52</cn>
</apply>
<apply>
<eq/>
<ci>gk2</ci>
<cn cellml:units="nS_per_pF">0.032</cn>
</apply>
<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>
<apply>
<eq/>
<ci>gka</ci>
<cn cellml:units="nS_per_pF">0.16</cn>
</apply>
<apply>
<eq/>
<ci>gcl</ci>
<cn cellml:units="nS_per_pF">0.1875</cn>
</apply>
<apply>
<eq/>
<ci>gh</ci>
<cn cellml:units="nS_per_pF">0.0542</cn>
</apply>
<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>
<apply>
<eq/>
<ci>ginak</ci>
<cn cellml:units="pA_per_pF">1.7</cn>
</apply>
<apply>
<eq/>
<ci>nakKmko</ci>
<cn cellml:units="mM">2.0</cn>
</apply>
<apply>
<eq/>
<ci>nakKmnai</ci>
<cn cellml:units="mM">22.0</cn>
</apply>
<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>
<!-- var v: mV {pub: in};
var cai: mM {pub: in};-->
</component>
<!-- // Model - Force and w variable-->
<!-- def comp Ca_dependent_Force as
var wss: dimensionless {pub: none};
var wtc: msec {pub: none};
var conversion: nM_to_mM {init: 1e-6, pub: none};
var Force: uN {pub: none};
var w: dimensionless {init: 0.2345238135343783, pub: none};
var time: msec {pub: in};
var cai: mM {pub: in};
var Fmax: uN {pub: in};
var FKm: nM {pub: in};
var Fn: dimensionless {pub: in};-->
<!-- wss = 1.0{dimensionless}/(1.0{dimensionless}+pow(FKm*conversion/cai, Fn));
wtc = 4000.0{msec}*(0.234845{dimensionless}+(1.0{dimensionless}-0.234845{dimensionless})/(1.0{dimensionless}+pow(cai/(FKm*conversion), Fn)));
Force = Fmax*(w-0.2345{dimensionless});
ode(w, time) = (wss-w)/wtc;
enddef;-->
<!-- def map between Ca_dependent_Force and environment for
vars time and time;
enddef;-->
<!-- def map between Ca_dependent_Force and Ca_Concentrations for
vars cai and cai;
enddef;-->
<!-- def map between Ca_dependent_Force and parameters for
vars Fmax and Fmax;
vars FKm and FKm;
vars Fn and Fn;
enddef;-->
<!-- Model - currents-->
<!-- def comp I_Na as
var mss: dimensionless {pub: none};
var hss: dimensionless {pub: none};
var mtc: msec {pub: none};
var htc: msec {pub: none};
var m: dimensionless {init: 0.1253518889572223, pub: none};
var h: dimensionless {init: 0.404599170710196, pub: none};
var ina: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ena: mV {pub: in};
var gna: nS_per_pF {pub: in};-->
<!-- mss = 1.0{dimensionless}/(1.0{dimensionless}+exp(-(v+35.9584{mV})/9.24013{mV}));
hss = 1.0{dimensionless}/(1.0{dimensionless}+exp((v+57.0{mV})/8.0{mV}));
mtc = 0.25{msec}+7.0{msec}/(1.0{dimensionless}+exp((v+38.0{mV})/10.0{mV}));
htc = 0.9{msec}+1002.85{msec}/(1.0{dimensionless}+sqr((v+47.5{mV})/1.5{mV}));
ina = gna*m*m*m*h*(v-ena);
ode(m, time) = (mss-m)/mtc;
ode(h, time) = (hss-h)/htc;
enddef;-->
<!-- def map between I_Na and environment for
vars time and time;
enddef;-->
<!-- def map between I_Na and membrane_potential for
vars v and v;
vars ina and ina;
enddef;-->
<!-- def map between I_Na and parameters for
vars ena and ena;
vars gna and gna;
enddef;-->
<!-- Model - currents-->
<component name="I_CaL">
<!-- var dss: dimensionless {pub: none};
var fss: dimensionless {pub: none};
var fca: dimensionless {pub: none};
var dtc: msec {pub: none};
var f1tc: msec {pub: none};
var f2tc: msec {pub: none};-->
<!-- var dss_clamp: dimensionless{pub: out};
var fss_clamp: dimensionless {pub: out};
var dtc_clamp: msec{pub: out};
var f1tc_clamp: msec {pub: out};
var f2tc_clamp: msec {pub: out};-->
<variable name="dss_vv" public_interface="out" units="dimensionless"/>
<variable name="fss_vv" public_interface="out" units="dimensionless"/>
<variable name="dtc_vv" public_interface="out" units="msec"/>
<variable name="f1tc_vv" public_interface="out" units="msec"/>
<variable name="f2tc_vv" public_interface="out" units="msec"/>
<variable initial_value="-60" name="V_hold" public_interface="none" units="mV"/>
<variable initial_value="0.0" name="V_test" public_interface="none" units="mV"/>
<!--origin -0.5-->
<variable initial_value="1000" name="test_s" public_interface="none" units="msec"/>
<variable initial_value="6000" name="test_e" public_interface="none" units="msec"/>
<!-- var v_clamp: mV {pub: out};-->
<!-- v_clamp = sel
case (time > test_s) and (time < test_e):
V_test;
otherwise:
V_hold;
endsel;
var d: dimensionless {init: 0.01036961357784695, pub: none};
var f1: dimensionless {init: 0.9065941499695301, pub: none};
var f2: dimensionless {init: 0.9065967263076083, pub: none};-->
<variable initial_value="0.0102" name="d" public_interface="none" units="dimensionless"/>
<variable initial_value="0.9077" name="f1" public_interface="none" units="dimensionless"/>
<variable initial_value="0.9077" name="f2" public_interface="none" units="dimensionless"/>
<!--var ical: pA_per_pF {pub: out};-->
<!-- var v: mV {pub: in};-->
<variable name="vv" public_interface="in" units="mV"/>
<!--var cai: mM {pub: in};
var ecal: mV {pub: in};-->
<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>dss_vv</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>vv</ci>
<cn cellml:units="mV">22.0</cn>
</apply>
</apply>
<cn cellml:units="mV">7.0</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- checked-->
<apply>
<eq/>
<ci>fss_vv</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>vv</ci>
<cn cellml:units="mV">38.0</cn>
</apply>
<cn cellml:units="mV">7.0</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- checked-->
<apply>
<eq/>
<ci>dtc_vv</ci>
<apply>
<plus/>
<cn cellml:units="msec">2.29</cn>
<apply>
<divide/>
<cn cellml:units="msec">5.7</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1.0</cn>
<apply>
<power/>
<apply>
<divide/>
<apply>
<plus/>
<ci>vv</ci>
<cn cellml:units="mV">29.97</cn>
</apply>
<cn cellml:units="mV">9.0</cn>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- checked-->
<apply>
<eq/>
<ci>f1tc_vv</ci>
<cn cellml:units="msec">12.0</cn>
</apply>
<!-- checked-->
<apply>
<eq/>
<ci>f2tc_vv</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>vv</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>vv</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>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>vv</ci>
</bvar>
<ci>d</ci>
</apply>
<cn cellml:units="dimensionless">0</cn>
</apply>
</math>
<!-- ode(d, time) = (dss-d)/dtc;
ode(f1, time) = (fss-f1)/f1tc;
ode(f2, time) = (fss-f2)/f2tc;
var f22: dimensionless {init: 0.9077, pub: none};
ode(f22, vv) = 0{dimensionless};-->
</component>
<connection>
<map_components component_1="I_CaL" component_2="environment"/>
<map_variables variable_1="vv" variable_2="vv"/>
</connection>
<!-- def map between I_CaL and membrane_potential for
vars v and v;
vars ical and ical;
enddef;-->
<!-- def map between I_CaL and Ca_Concentrations for
vars cai and cai;
enddef;-->
<connection>
<map_components component_1="I_CaL" component_2="parameters"/>
<!-- vars ecal and ecal;-->
<map_variables variable_1="gcal" variable_2="gcal"/>
<map_variables variable_1="kmca" variable_2="kmca"/>
</connection>
<!-- Model - currents-->
<!-- def comp I_CaT as
var bss: dimensionless {pub: none};
var gss: dimensionless {pub: none};
var btc: msec {pub: none};
var gtc: msec {pub: none};
var b: dimensionless {init: 0.508117603077852, pub: none};
var g: dimensionless {init: 0.03582573962705717, pub: none};
var icat: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ecat: mV {pub: in};
var gcat: nS_per_pF {pub: in};-->
<!-- bss = 1.0{dimensionless}/(1.0{dimensionless}+exp(-(v+54.23{mV})/9.88{mV}));
gss = 0.02{dimensionless}+0.98{dimensionless}/(1.0{dimensionless}+exp((v+72.978{mV})/4.64{mV}));
btc = 0.45{msec}+3.9{msec}/(1.0{dimensionless}+sqr((v+66.0{mV})/26.0{mV}));
gtc = 150.0{msec}*(1.0{dimensionless}-1.0{dimensionless}/((1.0{dimensionless}+exp((v-417.43{mV})/203.18{mV}))*(1.0{dimensionless}+exp(-(v+61.11{mV})/8.07{mV}))));
icat = gcat*b*b*g*(v-ecat);
ode(b, time) = (bss-b)/btc;
ode(g, time) = (gss-g)/gtc;
enddef;-->
<!-- def map between I_CaT and environment for
vars time and time;
enddef;-->
<!-- def map between I_CaT and membrane_potential for
vars v and v;
vars icat and icat;
enddef;-->
<!-- def map between I_CaT and parameters for
vars ecat and ecat;
vars gcat and gcat;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_b as
var ib: pA_per_pF {pub: out};
var v: mV {pub: in};
var ek: mV {pub: in};
var gb: nS_per_pF {pub: in};-->
<!-- ib = gb*(v-ek);
enddef;-->
<!-- def map between I_b and membrane_potential for
vars v and v;
vars ib and ib;
enddef;-->
<!-- def map between I_b and parameters for
vars ek and ek;
vars gb and gb;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_K1 as
var qss: dimensionless {pub: none};
var rss: dimensionless {pub: none};
var qtc: msec {pub: none};
var r1tc: msec {pub: none};
var r2tc: msec {pub: none};
var q: dimensionless {init: 0.2060363247740295, pub: none};
var r1: dimensionless {init: 0.1922244113609531, pub: none};
var r2: dimensionless {init: 0.1932803618375963, pub: none};
var ik1: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ek: mV {pub: in};
var gk1: nS_per_pF {pub: in};-->
<!-- qss = 0.978613{dimensionless}/(1.0{dimensionless}+exp(-(v+18.6736{mV})/26.6606{mV}));
rss = 1.0{dimensionless}/(1.0{dimensionless}+exp((v+63.0{mV})/6.3{mV}));
qtc = 500.0{msec}/(1.0{dimensionless}+sqr((v+60.71{mV})/15.79{mV}));
r1tc = 5000.0{msec}/(1.0{dimensionless}+sqr((v+62.7133{mV})/35.8611{mV}));
r2tc = 30000.0{msec}+220000.0{msec}/(1.0{dimensionless}+exp((v+22.0{mV})/4.0{mV}));
ik1 = gk1*q*q*(0.38{dimensionless}*r1+0.63{dimensionless}*r2)*(v-ek);
ode(q, time) = (qss-q)/qtc;
ode(r1, time) = (rss-r1)/r1tc;
ode(r2, time) = (rss-r2)/r2tc;
enddef;-->
<!-- def map between I_K1 and environment for
vars time and time;
enddef;-->
<!-- def map between I_K1 and membrane_potential for
vars v and v;
vars ik1 and ik1;
enddef;-->
<!-- def map between I_K1 and parameters for
vars ek and ek;
vars gk1 and gk1;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_K2 as
var pss: dimensionless {pub: none};
var kss: dimensionless {pub: none};
var ptc: msec {pub: none};
var k1tc: msec {pub: none};
var k2tc: msec {pub: none};
var p: dimensionless {init: 0.1174074734567931, pub: none};
var k1: dimensionless {init: 0.9968385770271651, pub: none};
var k2: dimensionless {init: 0.9968408069904307, pub: none};
var ik2: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ek: mV {pub: in};
var gk2: nS_per_pF {pub: in};-->
<!-- pss = 0.948{dimensionless}/(1.0{dimensionless}+exp(-(v+17.91{mV})/18.4{mV}));
kss = 1.0{dimensionless}/(1.0{dimensionless}+exp((v+21.2{mV})/5.7{mV}));
ptc = 100.0{msec}/(1.0{dimensionless}+sqr((v+64.1{mV})/28.67{mV}));
k1tc = 1000000.0{msec}*(1.0{dimensionless}-1.0{dimensionless}/((1.0{dimensionless}+exp((v-315.0{mV})/50.0{mV}))*(1.0{dimensionless}+exp(-(v+74.9{mV})/8.0{mV}))));
k2tc = 2500000.0{msec}*(1.0{dimensionless}-1.0{dimensionless}/((1.0{dimensionless}+exp((v-132.868{mV})/25.3992{mV}))*(1.0{dimensionless}+exp(-(v+24.9203{mV})/2.67915{mV}))));
ik2 = gk2*p*p*(0.75{dimensionless}*k1+0.25{dimensionless}*k2)*(v-ek);
ode(p, time) = (pss-p)/ptc;
ode(k1, time) = (kss-k1)/k1tc;
ode(k2, time) = (kss-k2)/k2tc;
enddef;-->
<!-- def map between I_K2 and environment for
vars time and time;
enddef;-->
<!-- def map between I_K2 and membrane_potential for
vars v and v;
vars ik2 and ik2;
enddef;-->
<!-- def map between I_K2 and parameters for
vars ek and ek;
vars gk2 and gk2;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_Ka as
var sss: dimensionless {pub: none};
var xss: dimensionless {pub: none};
var stc: msec {pub: none};
var xtc: msec {pub: none};
var s: dimensionless {init: 0.0307583106982354, pub: none};
var x: dimensionless {init: 0.08785242843398365, pub: none};
var ika: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ek: mV {pub: in};
var gka: nS_per_pF {pub: in};-->
<!-- sss = 1.0{dimensionless}/(1.0{dimensionless}+exp(-(v+27.79{mV})/7.57{mV}));
xss = 0.02{dimensionless}+0.98{dimensionless}/(1.0{dimensionless}+exp((v+69.5{mV})/6.0{mV}));
stc = 17.0{msec}/(1.0{dimensionless}+sqr((v+20.5232{mV})/35.0{mV}));
xtc = 7.5{msec}+10.0{msec}/(1.0{dimensionless}+sqr((v+34.1765{mV})/120.0{mV}));
ika = gka*s*x*(v-ek);
ode(s, time) = (sss-s)/stc;
ode(x, time) = (xss-x)/xtc;
enddef;-->
<!-- def map between I_Ka and environment for
vars time and time;
enddef;-->
<!-- def map between I_Ka and membrane_potential for
vars v and v;
vars ika and ika;
enddef;-->
<!-- def map between I_Ka and parameters for
vars ek and ek;
vars gka and gka;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_BKa as
var xass_z: dimensionless {pub: none};
var xass_vh: mV {pub: none};
var conversion: mM_to_M {init: 1000.0, pub: none};
var xass: dimensionless {pub: none};
var xatc: msec {pub: none};
var xa: dimensionless {init: 0.0003569126518797985, pub: none};
var iBKa: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var cai: mM {pub: in};
var ek: mV {pub: in};
var gkca: nS_per_pF {pub: in};
var gbka: dimensionless {pub: in};
var R: joule_per_kelvin_per_kilomole {pub: in};
var frdy: coulomb_per_mole {pub: in};
var temp: kelvin {pub: in};-->
<!-- xass_z = -0.749234{dimensionless}/(1.0{dimensionless}+sqr((cai*conversion-0.0630535{mole})/0.161942{mole}))+8.38384{dimensionless}/(1.0{dimensionless}+sqr((cai*conversion+1538.29{mole})/739.057{mole}));
xass_vh = 5011.47{mV}/(1.0{dimensionless}+pow((cai*conversion+0.237503{mole})/0.000239278{mole}, 0.42291{dimensionless}))-37.5137{mV};
xass = 1.0{dimensionless}/(1.0{dimensionless}+exp(-xass_z*frdy*(v-xass_vh)/(R*temp)));
xatc = 2.40914{msec}/(1.0{dimensionless}+sqr((v-158.779{mV})/-52.1497{mV}));
iBKa = gkca*gbka*xa*(v-ek);
ode(xa, time) = (xass-xa)/xatc;
enddef;-->
<!-- def map between I_BKa and environment for
vars time and time;
enddef;-->
<!-- def map between I_BKa and membrane_potential for
vars v and v;
vars iBKa and iBKa;
enddef;-->
<!-- def map between I_BKa and Ca_Concentrations for
vars cai and cai;
enddef;-->
<!-- def map between I_BKa and parameters for
vars ek and ek;
vars gkca and gkca;
vars gbka and gbka;
vars R and R;
vars frdy and frdy;
vars temp and temp;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_BKab as
var xabss_z: dimensionless {pub: none};
var xabss_vh: mV {pub: none};
var conversion: mM_to_M {init: 1000.0, pub: none};
var xabss: dimensionless {pub: none};
var xabtc: msec {pub: none};
var xab: dimensionless {init: 0.002220456569762898, pub: none};
var iBKab: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var cai: mM {pub: in};
var ek: mV {pub: in};
var gkca: nS_per_pF {pub: in};
var gbkab: dimensionless {pub: in};
var R: joule_per_kelvin_per_kilomole {pub: in};
var frdy: coulomb_per_mole {pub: in};
var temp: kelvin {pub: in};-->
<!-- xabss_z = -0.681249{dimensionless}/(1.0{dimensionless}+sqr((cai*conversion-0.218988{mole})/0.428335{mole}))+1.40001{dimensionless}/(1.0{dimensionless}+sqr((cai*conversion+228.71{mole})/684.946{mole}));
xabss_vh = 8540.23{mV}/(1.0{dimensionless}+pow((cai*conversion+0.401189{mole})/0.00399115{mole}, 0.668054{dimensionless}))-109.275{mV};
xabss = 1.0{dimensionless}/(1.0{dimensionless}+exp(-xabss_z*frdy*(v-xabss_vh)/(R*temp)));
xabtc = 13.8049{msec}/(1.0{dimensionless}+sqr((v-153.019{mV})/66.4952{mV}));
iBKab = gkca*gbkab*xab*(v-ek);
ode(xab, time) = (xabss-xab)/xabtc;
enddef;-->
<!-- def map between I_BKab and environment for
vars time and time;
enddef;-->
<!-- def map between I_BKab and membrane_potential for
vars v and v;
vars iBKab and iBKab;
enddef;-->
<!-- def map between I_BKab and Ca_Concentrations for
vars cai and cai;
enddef;-->
<!-- def map between I_BKab and parameters for
vars ek and ek;
vars gkca and gkca;
vars gbkab and gbkab;
vars R and R;
vars frdy and frdy;
vars temp and temp;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_h as
var yss: dimensionless {pub: none};
var ytc: msec {pub: none};
var ya: per_msec {pub: none};
var yb: per_msec {pub: none};
var y: dimensionless {init: 0.002604864867063448, pub: none};
var ih: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var eh: mV {pub: in};
var gh: nS_per_pF {pub: in};-->
<!-- yss = 1.0{dimensionless}/(1.0{dimensionless}+exp((v+105.39{mV})/8.6553{mV}));
ya = 0.0000035{per_msec}*exp(-0.0497{per_mV}*v);
yb = 0.04003{per_msec}*exp(0.05211{per_mV}*v);
ytc = 1.0{dimensionless}/(ya+yb);
ih = gh*y*(v-eh);
ode(y, time) = (yss-y)/ytc;
enddef;-->
<!-- def map between I_h and environment for
vars time and time;
enddef;-->
<!-- def map between I_h and membrane_potential for
vars v and v;
vars ih and ih;
enddef;-->
<!-- def map between I_h and parameters for
vars eh and eh;
vars gh and gh;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_Cl as
var css: dimensionless {pub: none};
var ctc: msec {pub: none};
var K1cl: mM {pub: none};
var K2cl: dimensionless {pub: none};
var c: dimensionless {init: 0.0003764413740731269, pub: none};
var icl: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var cai: mM {pub: in};
var ecl: mV {pub: in};
var gcl: nS_per_pF {pub: in};
var vFRT: dimensionless {pub: in};-->
<!-- K1cl = 0.0006{mM}*exp(2.53{dimensionless}*vFRT);
K2cl = 0.1{dimensionless}*exp(-5.0{dimensionless}*vFRT);
css = 1.0{dimensionless}/(1.0{dimensionless}+K2cl*(sqr(K1cl/cai)+K1cl/cai+1.0{dimensionless}));
ctc = -160.0{msec}+210.0{msec}/(1.0{dimensionless}+exp((v+4.56{mV})/11.62{mV}))+170.0{msec}/(1.0{dimensionless}+exp(-(v+25.5{mV})/11.62{mV}));
icl = gcl*c*(v-ecl);
ode(c, time) = (css-c)/ctc;
enddef;-->
<!-- def map between I_Cl and environment for
vars time and time;
enddef;-->
<!-- def map between I_Cl and membrane_potential for
vars v and v;
vars icl and icl;
enddef;-->
<!-- def map between I_Cl and Ca_Concentrations for
vars cai and cai;
enddef;-->
<!-- def map between I_Cl and parameters for
vars ecl and ecl;
vars gcl and gcl;
vars vFRT and vFRT;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_ns as
var fmg: dimensionless {pub: none};
var gs_nao: dimensionless {pub: none};
var gs_cao: dimensionless {pub: none};
var gs_ko: dimensionless {pub: none};
var tinyamount: mM {init: 1e-8, pub: none};
var insna: pA_per_pF {pub: out};
var insca: pA_per_pF {pub: out};
var insk: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ko: mM {pub: in};
var cao: mM {pub: in};
var nao: mM {pub: in};
var mgo: mM {pub: in};
var enscc: mV {pub: in};
var gns: nS_per_pF {pub: in};
var gnsCa: dimensionless {pub: in};
var gnsNa: dimensionless {pub: in};
var gnsK: dimensionless {pub: in};-->
<!-- fmg = 0.108043{dimensionless}+0.903902{dimensionless}/(1.0{dimensionless}+pow(mgo/0.281007{mM}, 1.29834{dimensionless}));
gs_nao = 1{dimensionless}/0.0123{dimensionless}*0.03{dimensionless}/(1.0{dimensionless}+sqr(150.0{mM}/(nao+tinyamount)));
gs_ko = 1{dimensionless}/0.0123{dimensionless}*0.03{dimensionless}/(1.0{dimensionless}+sqr(150.0{mM}/(ko+tinyamount)));
gs_cao = 1{dimensionless}/0.000525{dimensionless}*0.03{dimensionless}/(1.0{dimensionless}+sqr(150.0{mM}/(cao+tinyamount)));
insca = fmg*gs_cao*gnsCa*gns*(v-enscc);
insna = fmg*gs_nao*gnsNa*gns*(v-enscc);
insk = fmg*gs_ko*gnsK*gns*(v-enscc);
enddef;-->
<!-- def map between I_ns and environment for
vars time and time;
enddef;-->
<!-- def map between I_ns and membrane_potential for
vars v and v;
vars insna and insna;
vars insca and insca;
vars insk and insk;
enddef;-->
<!-- def map between I_ns and parameters for
vars enscc and enscc;
vars gns and gns;
vars gnsNa and gnsNa;
vars gnsCa and gnsCa;
vars gnsK and gnsK;
vars nao and nao;
vars cao and cao;
vars ko and ko;
vars mgo and mgo;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_NaK as
var fnak: dimensionless {pub: none};
var knak: dimensionless {pub: none};
var nnak: dimensionless {pub: none};
var inak: pA_per_pF {pub: out};
var time: msec {pub: in};
var v: mV {pub: in};
var ginak: pA_per_pF {pub: in};
var ko: mM {pub: in};
var nao: mM {pub: in};
var nai: mM {pub: in};
var nakKmko: mM {pub: in};
var nakKmnai: mM {pub: in};
var vFRT: dimensionless {pub: in};-->
<!-- fnak = 1.0{dimensionless}/(1.0{dimensionless}+0.1245{dimensionless}*exp(-0.1{dimensionless}*vFRT)+0.00219{dimensionless}*exp(nao/49.71{mM})*exp(-1.9{dimensionless}*vFRT));
knak = 1.0{dimensionless}/(1.0{dimensionless}+pow(nakKmko/ko, 1.5{dimensionless}));
nnak = 1.0{dimensionless}/(1.0{dimensionless}+sqr(nakKmnai/nai));
inak = ginak*knak*nnak*fnak;
enddef;-->
<!-- def map between I_NaK and environment for
vars time and time;
enddef;-->
<!-- def map between I_NaK and membrane_potential for
vars v and v;
vars inak and inak;
enddef;-->
<!-- def map between I_NaK and parameters for
vars nao and nao;
vars nai and nai;
vars ko and ko;
vars ginak and ginak;
vars nakKmko and nakKmko;
vars nakKmnai and nakKmnai;
vars vFRT and vFRT;
enddef;-->
<!-- // Model - currents-->
<!-- def comp I_NaCa as
var inaca_sign: dimensionless {init: -1, pub: none};
var inaca: pA_per_pF {pub: out};
var jnaca: mM_per_msec {pub: in};
var buff: dimensionless {pub: in};
var AV: cm2_per_uL {pub: in};
var zca: dimensionless {pub: in};
var frdy: coulomb_per_mole {pub: in};
var Cm: uF_per_cm2 {pub: in};-->
<!-- inaca = 0.5{dimensionless}*zca*frdy/(AV*Cm*buff)*inaca_sign*jnaca;
enddef;-->
<!-- def map between I_NaCa and membrane_potential for
vars inaca and inaca;
enddef;-->
<!-- def map between I_NaCa and J_NaCa for
vars jnaca and jnaca;
enddef;-->
<!-- def map between I_NaCa and parameters for
vars buff and buff;
vars AV and AV;
vars zca and zca;
vars frdy and frdy;
vars Cm and Cm;
enddef;-->
<!-- // Model - flux-->
<!-- def comp J_NaCa as
var f1naca: dimensionless {pub: none};
var f2naca: dimensionless {pub: none};
var fallo: dimensionless {pub: none};
var naca_Eup: mM4 {pub: none};
var naca_Ed1: dimensionless {pub: none};
var naca_Ed2: mM4 {pub: none};
var naca_Ed3: mM4 {pub: none};
var jnaca_sign: dimensionless {init: -1, pub: none};
var jnaca: mM_per_msec {pub: out};
var v: mV {pub: in};
var cai: mM {pub: in};
var cao: mM {pub: in};
var nao: mM {pub: in};
var nai: mM {pub: in};
var Jnaca: mM_per_msec {pub: in};
var Kmallo: mM {pub: in};
var nallo: dimensionless {pub: in};
var ksat: dimensionless {pub: in};
var xgamma: dimensionless {pub: in};
var Kmnai: mM {pub: in};
var Kmcai: mM {pub: in};
var Kmnao: mM {pub: in};
var Kmcao: mM {pub: in};
var vFRT: dimensionless {pub: in};-->
<!-- f1naca = exp((xgamma-1.0{dimensionless})*vFRT);
f2naca = exp(xgamma*vFRT);
fallo = 1.0{dimensionless}/(1.0{dimensionless}+pow(Kmallo/cai, nallo));
naca_Eup = pow(nai, 3{dimensionless})*cao*f2naca-pow(nao, 3{dimensionless})*cai*f1naca;
naca_Ed1 = 1.0{dimensionless}+ksat*f1naca;
naca_Ed2 = Kmcao*pow(nai, 3{dimensionless})+pow(Kmnao, 3{dimensionless})*cai+pow(Kmnai, 3{dimensionless})*cao*(1.0{dimensionless}+cai/Kmcai);
naca_Ed3 = cao*pow(nai, 3{dimensionless})+pow(nao, 3{dimensionless})*cai+pow(nao, 3{dimensionless})*Kmcai*(1.0{dimensionless}+pow(nai/Kmnai, 3{dimensionless}));
jnaca = jnaca_sign*Jnaca*fallo*naca_Eup/(naca_Ed1*(naca_Ed2+naca_Ed3));
enddef;-->
<!-- def map between J_NaCa and membrane_potential for
vars v and v;
enddef;-->
<!-- def map between J_NaCa and Ca_Concentrations for
vars cai and cai;
vars jnaca and jnaca;
enddef;-->
<!-- def map between J_NaCa and parameters for
vars nao and nao;
vars nai and nai;
vars cao and cao;
vars Jnaca and Jnaca;
vars Kmallo and Kmallo;
vars nallo and nallo;
vars ksat and ksat;
vars xgamma and xgamma;
vars Kmnai and Kmnai;
vars Kmcai and Kmcai;
vars Kmnao and Kmnao;
vars Kmcao and Kmcao;
vars vFRT and vFRT;
enddef;-->
<!-- // Model - flux-->
<!-- def comp J_PMCA as
var jpmca: mM_per_msec {pub: out};
var cai: mM {pub: in};
var Jpmca: mM_per_msec {pub: in};
var Kmpmca: mM {pub: in};
var npmca: dimensionless {pub: in};-->
<!-- jpmca = Jpmca/(1.0{dimensionless}+pow(Kmpmca/cai, npmca));
enddef;-->
<!-- def map between J_PMCA and Ca_Concentrations for
vars cai and cai;
vars jpmca and jpmca;
enddef;-->
<!-- def map between J_PMCA and parameters for
vars Jpmca and Jpmca;
vars Kmpmca and Kmpmca;
vars npmca and npmca;
enddef;-->
<!-- Model - Group-->
<!-- def group as encapsulation for
// comp interface incl
comp environment;
comp parameters;-->
<!-- comp I_CaL;-->
<!--// endcomp;
enddef;-->
<!-- Model - Mapping variables-->
</model>