Location: BG_LCC @ 79e47069637e / BG_LCC.cellml
- Author:
- Shelley Fong <s.fong@auckland.ac.nz>
- Date:
- 2021-08-06 13:28:44+12:00
- Desc:
- Tidy code
- Permanent Source URI:
- https://models.physiomeproject.org/workspace/6d7/rawfile/79e47069637ef50896a1b52e8c60ff2152df2176/BG_LCC.cellml
<?xml version='1.0' encoding='UTF-8'?>
<!-- L_type Ca channel with gates
same as Pan and LRd
flux of charge (dq_membrane/dt) = sum of all currents-->
<!-- 5 Aug repurposing this code to use as general, modular LCC channel
6 Aug updating parameters with values for an isolated LCC system (no other channels present when finding BG parameters)-->
<model name="LCC" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#" xmlns:xlink="http://www.w3.org/1999/xlink">
<import xlink:href="units_and_constants/units_BG.cellml">
<units name="mM" units_ref="mM"/>
<units name="fmol" units_ref="fmol"/>
<units name="per_fmol" units_ref="per_fmol"/>
<units name="J_per_mol" units_ref="J_per_mol"/>
<units name="fmol_per_sec" units_ref="fmol_per_sec"/>
<units name="C_per_mol" units_ref="C_per_mol"/>
<units name="J_per_C" units_ref="J_per_C"/>
<units name="microm3" units_ref="microm3"/>
<units name="fF" units_ref="fF"/>
<units name="fC" units_ref="fC"/>
<units name="fA" units_ref="fA"/>
<units name="per_second" units_ref="per_second"/>
<units name="millivolt" units_ref="millivolt"/>
</import>
<component name="membrane">
<variable name="t" public_interface="out" units="second"/>
<variable name="Vmem" units="volt"/>
<variable initial_value="-13039" name="q_membrane" public_interface="out" units="fC"/>
<variable initial_value="153400" name="C_m" public_interface="out" units="fF"/>
<variable initial_value="2494.34" name="RT" public_interface="out" units="J_per_mol"/>
<variable initial_value="96485" name="F" public_interface="out" units="C_per_mol"/>
<variable initial_value="1e4" name="vol_o" public_interface="out" units="microm3"/>
<variable initial_value="9.3276" name="q_Cao" public_interface="out" units="fmol"/>
<variable initial_value="0.00456" name="q_Cai" public_interface="out" units="fmol"/>
<variable initial_value="27.9828" name="q_Ko" public_interface="out" units="fmol"/>
<variable initial_value="5510" name="q_Ki" public_interface="out" units="fmol"/>
<variable name="v_CaGHK_i" public_interface="in" units="fmol_per_sec"/>
<variable name="v_CaGHK_o" public_interface="in" units="fmol_per_sec"/>
<variable name="v_KGHK_i" public_interface="in" units="fmol_per_sec"/>
<variable name="v_KGHK_o" public_interface="in" units="fmol_per_sec"/>
<variable name="v_E" units="fA"/>
<variable initial_value="2" name="zCa" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="zK" public_interface="out" units="dimensionless"/>
<!--{pub: in};
v_E = zCa*F*(v_LCC_Ca1+v_LCC_Ca2+v_LCC_K1+v_LCC_K2+v_LCC_Na1+v_LCC_Na2);-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_E</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>F</ci>
<apply>
<plus/>
<ci>v_CaGHK_o</ci>
<ci>v_KGHK_o</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Cai</ci>
</apply>
<ci>v_CaGHK_i</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Cao</ci>
</apply>
<ci>v_CaGHK_o</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Ki</ci>
</apply>
<ci>v_KGHK_i</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Ko</ci>
</apply>
<ci>v_KGHK_o</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_membrane</ci>
</apply>
<apply>
<minus/>
<ci>v_E</ci>
</apply>
</apply>
<!--q_membrane = -10e4{fC}+floor(t/1{dimensionless})*2e3{dimensionless};-->
<apply>
<eq/>
<ci>Vmem</ci>
<apply>
<divide/>
<ci>q_membrane</ci>
<ci>C_m</ci>
</apply>
</apply>
</math>
<!-- I_LR = P_Ca * z*F*V_norm * (gamma_Cai * cCai_st * exp(V_norm) - gamma_Cao* cCao_st) / (exp(V_norm) - 1{dimensionless});
I_GHK = P_GHK * z*F*V_norm * (cCai_st * exp(V_norm) - cCao_st) / (exp(V_norm) - 1{dimensionless});
I_LR = P_K * z*F*V_norm * (gamma_Ki * cKi_st * exp(V_norm) - gamma_Ko* cKo_st) / (exp(V_norm) - 1{dimensionless}) ;
I_GHK = P_GHK * z*F*V_norm * (cKi_st * exp(V_norm) - cKo_st) / (exp(V_norm) - 1{dimensionless});-->
</component>
<component name="BG_LCC">
<variable name="t" private_interface="out" public_interface="in" units="second"/>
<variable name="RT" private_interface="out" public_interface="in" units="J_per_mol"/>
<variable name="F" private_interface="out" public_interface="in" units="C_per_mol"/>
<variable name="C_m" public_interface="in" units="fF"/>
<!-- BG constants (using GHK for ion permeability)-->
<variable initial_value="0.0102432" name="K_Cai" units="per_fmol"/>
<variable initial_value="0.0751143" name="K_Cao" units="per_fmol"/>
<variable initial_value="0.00578984" name="K_Ki" units="per_fmol"/>
<variable initial_value="0.0424574" name="K_Ko" units="per_fmol"/>
<variable initial_value="115.689" name="kappa_LCC_Ca1" units="fmol_per_sec"/>
<variable initial_value="216.336" name="kappa_LCC_Ca2" units="fmol_per_sec"/>
<variable initial_value="0.160891" name="kappa_LCC_K1" units="fmol_per_sec"/>
<variable initial_value="0.300862" name="kappa_LCC_K2" units="fmol_per_sec"/>
<variable name="Cao" units="mM"/>
<variable name="vol_o" public_interface="in" units="microm3"/>
<variable name="zCa" public_interface="in" units="dimensionless"/>
<variable name="zK" public_interface="in" units="dimensionless"/>
<!-- state variables-->
<variable name="q_Cai" public_interface="in" units="fmol"/>
<variable name="q_Cao" public_interface="in" units="fmol"/>
<variable name="q_Ki" public_interface="in" units="fmol"/>
<variable name="q_Ko" public_interface="in" units="fmol"/>
<variable name="q_E" public_interface="in" units="fC"/>
<variable name="mu_E" private_interface="out" units="J_per_C"/>
<variable name="mu_mem_Ca" units="J_per_mol"/>
<variable name="mu_mem_K" units="J_per_mol"/>
<variable name="v_CaGHK_i" public_interface="out" units="fmol_per_sec"/>
<variable name="v_CaGHK_o" public_interface="out" units="fmol_per_sec"/>
<variable name="v_KGHK_i" public_interface="out" units="fmol_per_sec"/>
<variable name="v_KGHK_o" public_interface="out" units="fmol_per_sec"/>
<!-- BG potential energy-->
<variable name="mu_Cai" private_interface="out" units="J_per_mol"/>
<variable name="mu_Cao" units="J_per_mol"/>
<variable name="mu_Ki" units="J_per_mol"/>
<variable name="mu_Ko" units="J_per_mol"/>
<variable name="mu_S111" private_interface="in" units="J_per_mol"/>
<variable name="mu_S121" private_interface="in" units="J_per_mol"/>
<!-- BG reaction parameters-->
<variable name="Af_LCC_Ca1" units="J_per_mol"/>
<variable name="Ar_LCC_Ca1" units="J_per_mol"/>
<variable name="Am_LCC_Ca1" units="J_per_mol"/>
<variable name="v_LCC_Ca1" public_interface="out" units="fmol_per_sec"/>
<variable name="Af_LCC_Ca2" units="J_per_mol"/>
<variable name="Ar_LCC_Ca2" units="J_per_mol"/>
<variable name="Am_LCC_Ca2" units="J_per_mol"/>
<variable name="v_LCC_Ca2" public_interface="out" units="fmol_per_sec"/>
<variable name="Af_LCC_K1" units="J_per_mol"/>
<variable name="Ar_LCC_K1" units="J_per_mol"/>
<variable name="Am_LCC_K1" units="J_per_mol"/>
<variable name="v_LCC_K1" public_interface="out" units="fmol_per_sec"/>
<variable name="Af_LCC_K2" units="J_per_mol"/>
<variable name="Ar_LCC_K2" units="J_per_mol"/>
<variable name="Am_LCC_K2" units="J_per_mol"/>
<variable name="v_LCC_K2" public_interface="out" units="fmol_per_sec"/>
<!-- fCa gate variable velocities-->
<variable name="v_fCa00" private_interface="in" units="fmol_per_sec"/>
<variable name="v_fCa01" private_interface="in" units="fmol_per_sec"/>
<variable name="v_fCa02" private_interface="in" units="fmol_per_sec"/>
<variable name="v_fCa10" private_interface="in" units="fmol_per_sec"/>
<variable name="v_fCa11" private_interface="in" units="fmol_per_sec"/>
<variable name="v_fCa12" private_interface="in" units="fmol_per_sec"/>
<!-- potential energies-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu_Cai</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cai</ci>
<ci>q_Cai</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Cao</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Cao</ci>
<ci>q_Cao</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Ki</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Ki</ci>
<ci>q_Ki</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_Ko</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_Ko</ci>
<ci>q_Ko</ci>
</apply>
</apply>
</apply>
</apply>
<!-- Constitutive relations - Convention: positive flow from inside to outside-->
<apply>
<eq/>
<ci>Cao</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>-3</cn>
<ci>q_Cao</ci>
</apply>
<ci>vol_o</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_mem_Ca</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>F</ci>
<ci>mu_E</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>mu_mem_K</ci>
<apply>
<times/>
<ci>zK</ci>
<ci>F</ci>
<ci>mu_E</ci>
</apply>
</apply>
<!-- Re_GHK:LCC_Ca1-->
<apply>
<eq/>
<ci>Af_LCC_Ca1</ci>
<apply>
<plus/>
<ci>mu_Cai</ci>
<ci>mu_mem_Ca</ci>
<ci>mu_S111</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_LCC_Ca1</ci>
<apply>
<plus/>
<ci>mu_Cao</ci>
<ci>mu_S111</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_Ca1</ci>
<ci>mu_mem_Ca</ci>
</apply>
<apply>
<eq/>
<ci>v_LCC_Ca1</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>kappa_LCC_Ca1</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_Ca1</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_Ca1</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_Ca1</ci>
<cn cellml:units="J_per_mol">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>kappa_LCC_Ca1</ci>
<ci>Am_LCC_Ca1</ci>
</apply>
<ci>RT</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Am_LCC_Ca1</ci>
<ci>RT</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_Ca1</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_Ca1</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<!-- Re_GHK:LCC_Ca2-->
<apply>
<eq/>
<ci>Af_LCC_Ca2</ci>
<apply>
<plus/>
<ci>mu_Cai</ci>
<ci>mu_mem_Ca</ci>
<ci>mu_S121</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_LCC_Ca2</ci>
<apply>
<plus/>
<ci>mu_Cao</ci>
<ci>mu_S121</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_Ca2</ci>
<ci>mu_mem_Ca</ci>
</apply>
<apply>
<eq/>
<ci>v_LCC_Ca2</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>kappa_LCC_Ca2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_Ca2</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_Ca2</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_Ca2</ci>
<cn cellml:units="J_per_mol">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>kappa_LCC_Ca2</ci>
<ci>Am_LCC_Ca2</ci>
</apply>
<ci>RT</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Am_LCC_Ca2</ci>
<ci>RT</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_Ca2</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_Ca2</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<!-- Re_GHK:LCC_K1-->
<apply>
<eq/>
<ci>Af_LCC_K1</ci>
<apply>
<plus/>
<ci>mu_Ki</ci>
<ci>mu_mem_K</ci>
<ci>mu_S111</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_LCC_K1</ci>
<apply>
<plus/>
<ci>mu_Ko</ci>
<ci>mu_S111</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_K1</ci>
<ci>mu_mem_K</ci>
</apply>
<apply>
<eq/>
<ci>v_LCC_K1</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>kappa_LCC_K1</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_K1</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_K1</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_K1</ci>
<cn cellml:units="J_per_mol">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>kappa_LCC_K1</ci>
<ci>Am_LCC_K1</ci>
</apply>
<ci>RT</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Am_LCC_K1</ci>
<ci>RT</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_K1</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_K1</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<!-- Re_GHK:LCC_K2-->
<apply>
<eq/>
<ci>Af_LCC_K2</ci>
<apply>
<plus/>
<ci>mu_Ki</ci>
<ci>mu_mem_K</ci>
<ci>mu_S121</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_LCC_K2</ci>
<apply>
<plus/>
<ci>mu_Ko</ci>
<ci>mu_S121</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_K2</ci>
<ci>mu_mem_K</ci>
</apply>
<apply>
<eq/>
<ci>v_LCC_K2</ci>
<piecewise>
<piece>
<apply>
<times/>
<ci>kappa_LCC_K2</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_K2</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_K2</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Am_LCC_K2</ci>
<cn cellml:units="J_per_mol">0</cn>
</apply>
</piece>
<otherwise>
<apply>
<times/>
<apply>
<divide/>
<apply>
<divide/>
<apply>
<times/>
<ci>kappa_LCC_K2</ci>
<ci>Am_LCC_K2</ci>
</apply>
<ci>RT</ci>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Am_LCC_K2</ci>
<ci>RT</ci>
</apply>
</apply>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_LCC_K2</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_LCC_K2</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</otherwise>
</piecewise>
</apply>
<apply>
<eq/>
<ci>mu_E</ci>
<apply>
<divide/>
<ci>q_E</ci>
<ci>C_m</ci>
</apply>
</apply>
<!-- ***********************************CHECK v_E *********************************// -->
<!--v_E = zCa*F*(v_LCC_Ca1+v_LCC_Ca2+v_LCC_K1+v_LCC_K2+v_LCC_Na1+v_LCC_Na2);-->
<!-- output total fluxes of ions-->
<apply>
<eq/>
<ci>v_CaGHK_i</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>zCa</ci>
<apply>
<plus/>
<ci>v_fCa10</ci>
<ci>v_fCa11</ci>
<ci>v_fCa02</ci>
<ci>v_fCa01</ci>
<ci>v_fCa00</ci>
<ci>v_fCa12</ci>
</apply>
</apply>
<ci>v_LCC_Ca1</ci>
</apply>
<ci>v_LCC_Ca2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_CaGHK_o</ci>
<apply>
<plus/>
<ci>v_LCC_Ca2</ci>
<ci>v_LCC_Ca1</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_KGHK_i</ci>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_LCC_K1</ci>
</apply>
<ci>v_LCC_K2</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>v_KGHK_o</ci>
<apply>
<plus/>
<ci>v_LCC_K2</ci>
<ci>v_LCC_K1</ci>
</apply>
</apply>
</math>
</component>
<component name="LCC_states_BG">
<variable name="t" public_interface="in" units="second"/>
<variable name="Vmem" public_interface="in" units="volt"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="RT" public_interface="in" units="J_per_mol"/>
<!-- gating parameters-->
<variable initial_value="132.798" name="kappa_d000" units="fmol_per_sec"/>
<variable initial_value="3.27634" name="kappa_d010" units="fmol_per_sec"/>
<variable initial_value="6.12666" name="kappa_d020" units="fmol_per_sec"/>
<variable initial_value="0.525901" name="kappa_d001" units="fmol_per_sec"/>
<variable initial_value="0.0129748" name="kappa_d011" units="fmol_per_sec"/>
<variable initial_value="0.0242625" name="kappa_d021" units="fmol_per_sec"/>
<variable initial_value="0.244648" name="kappa_f1_000" units="fmol_per_sec"/>
<variable initial_value="1.21486" name="kappa_f1_100" units="fmol_per_sec"/>
<variable initial_value="0.000968841" name="kappa_f1_001" units="fmol_per_sec"/>
<variable initial_value="0.00481102" name="kappa_f1_101" units="fmol_per_sec"/>
<variable initial_value="0.448058" name="kappa_f2_000" units="fmol_per_sec"/>
<variable initial_value="2.22495" name="kappa_f2_100" units="fmol_per_sec"/>
<variable initial_value="0.00177438" name="kappa_f2_001" units="fmol_per_sec"/>
<variable initial_value="0.00881111" name="kappa_f2_101" units="fmol_per_sec"/>
<variable initial_value="1258.66" name="kappa_f3_010" units="fmol_per_sec"/>
<variable initial_value="6250.17" name="kappa_f3_110" units="fmol_per_sec"/>
<variable initial_value="4.98447" name="kappa_f3_011" units="fmol_per_sec"/>
<variable initial_value="24.7516" name="kappa_f3_111" units="fmol_per_sec"/>
<variable initial_value="180068" name="kappa_fCa000" units="fmol_per_sec"/>
<variable initial_value="894171" name="kappa_fCa100" units="fmol_per_sec"/>
<variable initial_value="4442.55" name="kappa_fCa010" units="fmol_per_sec"/>
<variable initial_value="22060.6" name="kappa_fCa110" units="fmol_per_sec"/>
<variable initial_value="8307.45" name="kappa_fCa020" units="fmol_per_sec"/>
<variable initial_value="41252.7" name="kappa_fCa120" units="fmol_per_sec"/>
<variable initial_value="3.66543" name="K_000_LCC" units="per_fmol"/>
<variable initial_value="148.569" name="K_010_LCC" units="per_fmol"/>
<variable initial_value="79.4498" name="K_020_LCC" units="per_fmol"/>
<variable initial_value="0.738142" name="K_100_LCC" units="per_fmol"/>
<variable initial_value="29.9187" name="K_110_LCC" units="per_fmol"/>
<variable initial_value="15.9996" name="K_120_LCC" units="per_fmol"/>
<variable initial_value="925.578" name="K_001_LCC" units="per_fmol"/>
<variable initial_value="37516" name="K_011_LCC" units="per_fmol"/>
<variable initial_value="20062.3" name="K_021_LCC" units="per_fmol"/>
<variable initial_value="186.392" name="K_101_LCC" units="per_fmol"/>
<variable initial_value="7554.95" name="K_111_LCC" units="per_fmol"/>
<variable initial_value="4040.14" name="K_121_LCC" units="per_fmol"/>
<!-- BG mol quantities -->
<variable initial_value="3.6989e-07" name="q_S000_LCC" units="fmol"/>
<variable initial_value="7.3239e-05" name="q_S010_LCC" units="fmol"/>
<variable initial_value="3.6989e-07" name="q_S020_LCC" units="fmol"/>
<variable initial_value="3.7363e-09" name="q_S100_LCC" units="fmol"/>
<variable initial_value="7.3979e-07" name="q_S110_LCC" units="fmol"/>
<variable initial_value="3.7363e-09" name="q_S120_LCC" units="fmol"/>
<variable initial_value="4.1099e-08" name="q_S001_LCC" units="fmol"/>
<variable initial_value="8.1377e-06" name="q_S011_LCC" units="fmol"/>
<variable initial_value="4.1099e-08" name="q_S021_LCC" units="fmol"/>
<variable initial_value="4.1514e-10" name="q_S101_LCC" units="fmol"/>
<variable initial_value="8.2199e-08" name="q_S111_LCC" units="fmol"/>
<variable initial_value="4.1514e-10" name="q_S121_LCC" units="fmol"/>
<!-- gating variable potential energy-->
<variable name="mu_Cai" public_interface="in" units="J_per_mol"/>
<variable name="mu_S000_LCC" units="J_per_mol"/>
<variable name="mu_S010_LCC" units="J_per_mol"/>
<variable name="mu_S020_LCC" units="J_per_mol"/>
<variable name="mu_S100_LCC" units="J_per_mol"/>
<variable name="mu_S110_LCC" units="J_per_mol"/>
<variable name="mu_S120_LCC" units="J_per_mol"/>
<variable name="mu_S001_LCC" units="J_per_mol"/>
<variable name="mu_S011_LCC" units="J_per_mol"/>
<variable name="mu_S021_LCC" units="J_per_mol"/>
<variable name="mu_S101_LCC" units="J_per_mol"/>
<variable name="mu_S111_LCC" public_interface="out" units="J_per_mol"/>
<variable name="mu_S121_LCC" public_interface="out" units="J_per_mol"/>
<variable name="Af_fCa00" units="J_per_mol"/>
<variable name="Ar_fCa00" units="J_per_mol"/>
<variable name="Af_fCa01" units="J_per_mol"/>
<variable name="Ar_fCa01" units="J_per_mol"/>
<variable name="Af_fCa02" units="J_per_mol"/>
<variable name="Ar_fCa02" units="J_per_mol"/>
<variable name="Af_fCa10" units="J_per_mol"/>
<variable name="Ar_fCa10" units="J_per_mol"/>
<variable name="Af_fCa11" units="J_per_mol"/>
<variable name="Ar_fCa11" units="J_per_mol"/>
<variable name="Af_fCa12" units="J_per_mol"/>
<variable name="Ar_fCa12" units="J_per_mol"/>
<variable name="Af_d000" units="J_per_mol"/>
<variable name="Ar_d000" units="J_per_mol"/>
<variable name="Ar_d010" units="J_per_mol"/>
<variable name="Af_d010" units="J_per_mol"/>
<variable name="Af_d020" units="J_per_mol"/>
<variable name="Ar_d020" units="J_per_mol"/>
<variable name="Af_f1_000" units="J_per_mol"/>
<variable name="Ar_f1_000" units="J_per_mol"/>
<variable name="Af_f2_000" units="J_per_mol"/>
<variable name="Ar_f2_000" units="J_per_mol"/>
<variable name="Af_f3_010" units="J_per_mol"/>
<variable name="Ar_f3_010" units="J_per_mol"/>
<variable name="Af_f1_100" units="J_per_mol"/>
<variable name="Ar_f1_100" units="J_per_mol"/>
<variable name="Af_f2_100" units="J_per_mol"/>
<variable name="Ar_f2_100" units="J_per_mol"/>
<variable name="Af_f3_110" units="J_per_mol"/>
<variable name="Ar_f3_110" units="J_per_mol"/>
<variable name="Af_d001" units="J_per_mol"/>
<variable name="Ar_d001" units="J_per_mol"/>
<variable name="Af_d011" units="J_per_mol"/>
<variable name="Ar_d011" units="J_per_mol"/>
<variable name="Af_d021" units="J_per_mol"/>
<variable name="Ar_d021" units="J_per_mol"/>
<variable name="Af_f1_001" units="J_per_mol"/>
<variable name="Ar_f1_001" units="J_per_mol"/>
<variable name="Af_f2_001" units="J_per_mol"/>
<variable name="Ar_f2_001" units="J_per_mol"/>
<variable name="Af_f3_011" units="J_per_mol"/>
<variable name="Ar_f3_011" units="J_per_mol"/>
<variable name="Af_f1_101" units="J_per_mol"/>
<variable name="Ar_f1_101" units="J_per_mol"/>
<variable name="Af_f2_101" units="J_per_mol"/>
<variable name="Ar_f2_101" units="J_per_mol"/>
<variable name="Af_f3_111" units="J_per_mol"/>
<variable name="Ar_f3_111" units="J_per_mol"/>
<!-- velocities-->
<variable name="v_fCa00" public_interface="out" units="fmol_per_sec"/>
<variable name="v_fCa01" public_interface="out" units="fmol_per_sec"/>
<variable name="v_fCa02" public_interface="out" units="fmol_per_sec"/>
<variable name="v_fCa10" public_interface="out" units="fmol_per_sec"/>
<variable name="v_fCa11" public_interface="out" units="fmol_per_sec"/>
<variable name="v_fCa12" public_interface="out" units="fmol_per_sec"/>
<variable name="v_f1_000" units="fmol_per_sec"/>
<variable name="v_f2_000" units="fmol_per_sec"/>
<variable name="v_f3_010" units="fmol_per_sec"/>
<variable name="v_f1_100" units="fmol_per_sec"/>
<variable name="v_f2_100" units="fmol_per_sec"/>
<variable name="v_f3_110" units="fmol_per_sec"/>
<variable name="v_d000" units="fmol_per_sec"/>
<variable name="v_d010" units="fmol_per_sec"/>
<variable name="v_d020" units="fmol_per_sec"/>
<variable name="v_d001" units="fmol_per_sec"/>
<variable name="v_d011" units="fmol_per_sec"/>
<variable name="v_d021" units="fmol_per_sec"/>
<variable name="v_f1_001" units="fmol_per_sec"/>
<variable name="v_f2_001" units="fmol_per_sec"/>
<variable name="v_f3_011" units="fmol_per_sec"/>
<variable name="v_f1_101" units="fmol_per_sec"/>
<variable name="v_f2_101" units="fmol_per_sec"/>
<variable name="v_f3_111" units="fmol_per_sec"/>
<variable name="v_S000_LCC" units="fmol_per_sec"/>
<variable name="v_S010_LCC" units="fmol_per_sec"/>
<variable name="v_S020_LCC" units="fmol_per_sec"/>
<variable name="v_S100_LCC" units="fmol_per_sec"/>
<variable name="v_S110_LCC" units="fmol_per_sec"/>
<variable name="v_S120_LCC" units="fmol_per_sec"/>
<variable name="v_S001_LCC" units="fmol_per_sec"/>
<variable name="v_S011_LCC" units="fmol_per_sec"/>
<variable name="v_S021_LCC" units="fmol_per_sec"/>
<variable name="v_S101_LCC" units="fmol_per_sec"/>
<variable name="v_S111_LCC" units="fmol_per_sec"/>
<variable name="v_S121_LCC" units="fmol_per_sec"/>
<!-- zf and zr-->
<variable initial_value="2" name="zCa" units="dimensionless"/>
<variable initial_value="2.1404" name="z_fd" units="dimensionless"/>
<variable initial_value="-1.1495" name="z_ff1" units="dimensionless"/>
<variable initial_value="0.72162" name="z_ff2" units="dimensionless"/>
<variable initial_value="4.2933" name="z_ff3" units="dimensionless"/>
<variable initial_value="-2.1404" name="z_rd" units="dimensionless"/>
<variable initial_value="1.8993" name="z_rf1" units="dimensionless"/>
<variable initial_value="-0.52288" name="z_rf2" units="dimensionless"/>
<variable initial_value="0" name="z_rf3" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Af_fCa00</ci>
<ci>mu_S000_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa00</ci>
<apply>
<plus/>
<ci>mu_S001_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa00</ci>
<apply>
<times/>
<ci>kappa_fCa000</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa00</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa00</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_fCa01</ci>
<ci>mu_S010_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa01</ci>
<apply>
<plus/>
<ci>mu_S011_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa01</ci>
<apply>
<times/>
<ci>kappa_fCa010</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa01</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa01</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_fCa02</ci>
<ci>mu_S020_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa02</ci>
<apply>
<plus/>
<ci>mu_S021_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa02</ci>
<apply>
<times/>
<ci>kappa_fCa020</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa02</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa02</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_fCa10</ci>
<ci>mu_S100_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa10</ci>
<apply>
<plus/>
<ci>mu_S101_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa10</ci>
<apply>
<times/>
<ci>kappa_fCa100</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa10</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa10</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_fCa11</ci>
<ci>mu_S110_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa11</ci>
<apply>
<plus/>
<ci>mu_S111_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa11</ci>
<apply>
<times/>
<ci>kappa_fCa110</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa11</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa11</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_fCa12</ci>
<ci>mu_S120_LCC</ci>
</apply>
<apply>
<eq/>
<ci>Ar_fCa12</ci>
<apply>
<plus/>
<ci>mu_S121_LCC</ci>
<apply>
<times/>
<ci>zCa</ci>
<ci>mu_Cai</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_fCa12</ci>
<apply>
<times/>
<ci>kappa_fCa120</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_fCa12</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_fCa12</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- -->
<apply>
<eq/>
<ci>Af_d000</ci>
<apply>
<plus/>
<ci>mu_S000_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d000</ci>
<apply>
<plus/>
<ci>mu_S100_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d000</ci>
<apply>
<times/>
<ci>kappa_d000</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d000</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d000</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_d010</ci>
<apply>
<plus/>
<ci>mu_S010_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d010</ci>
<apply>
<plus/>
<ci>mu_S110_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d010</ci>
<apply>
<times/>
<ci>kappa_d010</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d010</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d010</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_d020</ci>
<apply>
<plus/>
<ci>mu_S020_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d020</ci>
<apply>
<plus/>
<ci>mu_S120_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d020</ci>
<apply>
<times/>
<ci>kappa_d020</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d020</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d020</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f1_000</ci>
<apply>
<plus/>
<ci>mu_S000_LCC</ci>
<apply>
<times/>
<ci>z_ff1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f1_000</ci>
<apply>
<plus/>
<ci>mu_S010_LCC</ci>
<apply>
<times/>
<ci>z_rf1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f1_000</ci>
<apply>
<times/>
<ci>kappa_f1_000</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f1_000</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f1_000</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f2_000</ci>
<apply>
<plus/>
<ci>mu_S000_LCC</ci>
<apply>
<times/>
<ci>z_ff2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f2_000</ci>
<apply>
<plus/>
<ci>mu_S020_LCC</ci>
<apply>
<times/>
<ci>z_rf2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f2_000</ci>
<apply>
<times/>
<ci>kappa_f2_000</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f2_000</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f2_000</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f3_010</ci>
<apply>
<plus/>
<ci>mu_S010_LCC</ci>
<apply>
<times/>
<ci>z_ff3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f3_010</ci>
<apply>
<plus/>
<ci>mu_S020_LCC</ci>
<apply>
<times/>
<ci>z_rf3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f3_010</ci>
<apply>
<times/>
<ci>kappa_f3_010</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f3_010</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f3_010</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f1_100</ci>
<apply>
<plus/>
<ci>mu_S100_LCC</ci>
<apply>
<times/>
<ci>z_ff1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f1_100</ci>
<apply>
<plus/>
<ci>mu_S110_LCC</ci>
<apply>
<times/>
<ci>z_rf1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f1_100</ci>
<apply>
<times/>
<ci>kappa_f1_100</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f1_100</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f1_100</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f2_100</ci>
<apply>
<plus/>
<ci>mu_S100_LCC</ci>
<apply>
<times/>
<ci>z_ff2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f2_100</ci>
<apply>
<plus/>
<ci>mu_S120_LCC</ci>
<apply>
<times/>
<ci>z_rf2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f2_100</ci>
<apply>
<times/>
<ci>kappa_f2_100</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f2_100</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f2_100</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f3_110</ci>
<apply>
<plus/>
<ci>mu_S110_LCC</ci>
<apply>
<times/>
<ci>z_ff3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f3_110</ci>
<apply>
<plus/>
<ci>mu_S120_LCC</ci>
<apply>
<times/>
<ci>z_rf3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f3_110</ci>
<apply>
<times/>
<ci>kappa_f3_110</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f3_110</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f3_110</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_d001</ci>
<apply>
<plus/>
<ci>mu_S001_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d001</ci>
<apply>
<plus/>
<ci>mu_S101_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d001</ci>
<apply>
<times/>
<ci>kappa_d001</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d001</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d001</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_d011</ci>
<apply>
<plus/>
<ci>mu_S011_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d011</ci>
<apply>
<plus/>
<ci>mu_S111_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d011</ci>
<apply>
<times/>
<ci>kappa_d011</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d011</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d011</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_d021</ci>
<apply>
<plus/>
<ci>mu_S021_LCC</ci>
<apply>
<times/>
<ci>z_fd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_d021</ci>
<apply>
<plus/>
<ci>mu_S121_LCC</ci>
<apply>
<times/>
<ci>z_rd</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_d021</ci>
<apply>
<times/>
<ci>kappa_d021</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_d021</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_d021</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f1_001</ci>
<apply>
<plus/>
<ci>mu_S001_LCC</ci>
<apply>
<times/>
<ci>z_ff1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f1_001</ci>
<apply>
<plus/>
<ci>mu_S011_LCC</ci>
<apply>
<times/>
<ci>z_rf1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f1_001</ci>
<apply>
<times/>
<ci>kappa_f1_001</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f1_001</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f1_001</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f2_001</ci>
<apply>
<plus/>
<ci>mu_S001_LCC</ci>
<apply>
<times/>
<ci>z_ff2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f2_001</ci>
<apply>
<plus/>
<ci>mu_S021_LCC</ci>
<apply>
<times/>
<ci>z_rf2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f2_001</ci>
<apply>
<times/>
<ci>kappa_f2_001</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f2_001</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f2_001</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f3_011</ci>
<apply>
<plus/>
<ci>mu_S011_LCC</ci>
<apply>
<times/>
<ci>z_ff3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f3_011</ci>
<apply>
<plus/>
<ci>mu_S021_LCC</ci>
<apply>
<times/>
<ci>z_rf3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f3_011</ci>
<apply>
<times/>
<ci>kappa_f3_011</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f3_011</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f3_011</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f1_101</ci>
<apply>
<plus/>
<ci>mu_S101_LCC</ci>
<apply>
<times/>
<ci>z_ff1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f1_101</ci>
<apply>
<plus/>
<ci>mu_S111_LCC</ci>
<apply>
<times/>
<ci>z_rf1</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f1_101</ci>
<apply>
<times/>
<ci>kappa_f1_101</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f1_101</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f1_101</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f2_101</ci>
<apply>
<plus/>
<ci>mu_S101_LCC</ci>
<apply>
<times/>
<ci>z_ff2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f2_101</ci>
<apply>
<plus/>
<ci>mu_S121_LCC</ci>
<apply>
<times/>
<ci>z_rf2</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f2_101</ci>
<apply>
<times/>
<ci>kappa_f2_101</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f2_101</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f2_101</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Af_f3_111</ci>
<apply>
<plus/>
<ci>mu_S111_LCC</ci>
<apply>
<times/>
<ci>z_ff3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ar_f3_111</ci>
<apply>
<plus/>
<ci>mu_S121_LCC</ci>
<apply>
<times/>
<ci>z_rf3</ci>
<ci>F</ci>
<ci>Vmem</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>v_f3_111</ci>
<apply>
<times/>
<ci>kappa_f3_111</ci>
<apply>
<minus/>
<apply>
<exp/>
<apply>
<divide/>
<ci>Af_f3_111</ci>
<ci>RT</ci>
</apply>
</apply>
<apply>
<exp/>
<apply>
<divide/>
<ci>Ar_f3_111</ci>
<ci>RT</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<!-- C:S000_LCC-->
<apply>
<eq/>
<ci>v_S000_LCC</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_d000</ci>
</apply>
<ci>v_f1_000</ci>
</apply>
<ci>v_f2_000</ci>
</apply>
<ci>v_fCa00</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S000_LCC</ci>
</apply>
<ci>v_S000_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S000_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_000_LCC</ci>
<ci>q_S000_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S010_LCC-->
<apply>
<eq/>
<ci>v_S010_LCC</ci>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_d010</ci>
</apply>
<ci>v_f3_010</ci>
</apply>
<ci>v_f1_100</ci>
</apply>
<ci>v_fCa01</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S010_LCC</ci>
</apply>
<ci>v_S010_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S010_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_010_LCC</ci>
<ci>q_S010_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S020_LCC-->
<apply>
<eq/>
<ci>v_S020_LCC</ci>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_d020</ci>
</apply>
<ci>v_f3_010</ci>
<ci>v_f2_000</ci>
</apply>
<ci>v_fCa02</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S020_LCC</ci>
</apply>
<ci>v_S020_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S020_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_020_LCC</ci>
<ci>q_S020_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S100_LCC-->
<apply>
<eq/>
<ci>v_S100_LCC</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_f2_100</ci>
</apply>
<ci>v_d000</ci>
</apply>
<ci>v_f1_100</ci>
</apply>
<ci>v_fCa10</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S100_LCC</ci>
</apply>
<ci>v_S100_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S100_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_100_LCC</ci>
<ci>q_S100_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S110_LCC-->
<apply>
<eq/>
<ci>v_S110_LCC</ci>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_f1_100</ci>
<ci>v_d010</ci>
</apply>
<ci>v_fCa11</ci>
</apply>
<ci>v_f3_110</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S110_LCC</ci>
</apply>
<ci>v_S110_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S110_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_110_LCC</ci>
<ci>q_S110_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S120_LCC-->
<apply>
<eq/>
<ci>v_S120_LCC</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_f2_100</ci>
<ci>v_d020</ci>
</apply>
<ci>v_fCa12</ci>
</apply>
<ci>v_f3_110</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S120_LCC</ci>
</apply>
<ci>v_S120_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S120_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_120_LCC</ci>
<ci>q_S120_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S001_LCC-->
<apply>
<eq/>
<ci>v_S001_LCC</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_d001</ci>
</apply>
<ci>v_f1_001</ci>
</apply>
<ci>v_f2_001</ci>
</apply>
<ci>v_fCa00</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S001_LCC</ci>
</apply>
<ci>v_S001_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S001_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_001_LCC</ci>
<ci>q_S001_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S011_LCC-->
<apply>
<eq/>
<ci>v_S011_LCC</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_d011</ci>
</apply>
<ci>v_f3_011</ci>
</apply>
<ci>v_f1_101</ci>
<ci>v_fCa01</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S011_LCC</ci>
</apply>
<ci>v_S011_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S011_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_011_LCC</ci>
<ci>q_S011_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S021_LCC-->
<apply>
<eq/>
<ci>v_S021_LCC</ci>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_d021</ci>
</apply>
<ci>v_f3_011</ci>
<ci>v_f2_001</ci>
<ci>v_fCa02</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S021_LCC</ci>
</apply>
<ci>v_S021_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S021_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_021_LCC</ci>
<ci>q_S021_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S101_LCC-->
<apply>
<eq/>
<ci>v_S101_LCC</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_f2_101</ci>
</apply>
<ci>v_d001</ci>
</apply>
<ci>v_f1_101</ci>
</apply>
<ci>v_fCa10</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S101_LCC</ci>
</apply>
<ci>v_S101_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S101_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_101_LCC</ci>
<ci>q_S101_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S111_LCC-->
<apply>
<eq/>
<ci>v_S111_LCC</ci>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_f1_101</ci>
<ci>v_d011</ci>
<ci>v_fCa11</ci>
</apply>
<ci>v_f3_111</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S111_LCC</ci>
</apply>
<ci>v_S111_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S111_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_111_LCC</ci>
<ci>q_S111_LCC</ci>
</apply>
</apply>
</apply>
</apply>
<!-- C:S121_LCC-->
<apply>
<eq/>
<ci>v_S121_LCC</ci>
<apply>
<plus/>
<ci>v_f2_101</ci>
<ci>v_d021</ci>
<ci>v_fCa12</ci>
<ci>v_f3_111</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_S121_LCC</ci>
</apply>
<ci>v_S121_LCC</ci>
</apply>
<apply>
<eq/>
<ci>mu_S121_LCC</ci>
<apply>
<times/>
<ci>RT</ci>
<apply>
<ln/>
<apply>
<times/>
<ci>K_121_LCC</ci>
<ci>q_S121_LCC</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<connection>
<map_components component_1="membrane" component_2="BG_LCC"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="RT" variable_2="RT"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="q_Cai" variable_2="q_Cai"/>
<map_variables variable_1="q_Cao" variable_2="q_Cao"/>
<map_variables variable_1="q_Ki" variable_2="q_Ki"/>
<map_variables variable_1="q_Ko" variable_2="q_Ko"/>
<map_variables variable_1="q_membrane" variable_2="q_E"/>
<map_variables variable_1="C_m" variable_2="C_m"/>
<map_variables variable_1="vol_o" variable_2="vol_o"/>
<map_variables variable_1="v_CaGHK_i" variable_2="v_CaGHK_i"/>
<map_variables variable_1="v_CaGHK_o" variable_2="v_CaGHK_o"/>
<map_variables variable_1="v_KGHK_i" variable_2="v_KGHK_i"/>
<map_variables variable_1="v_KGHK_o" variable_2="v_KGHK_o"/>
<map_variables variable_1="zCa" variable_2="zCa"/>
<map_variables variable_1="zK" variable_2="zK"/>
</connection>
<connection>
<map_components component_1="BG_LCC" component_2="LCC_states_BG"/>
<map_variables variable_1="t" variable_2="t"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="RT" variable_2="RT"/>
<map_variables variable_1="mu_E" variable_2="Vmem"/>
<map_variables variable_1="mu_Cai" variable_2="mu_Cai"/>
<map_variables variable_1="mu_S111" variable_2="mu_S111_LCC"/>
<map_variables variable_1="mu_S121" variable_2="mu_S121_LCC"/>
<map_variables variable_1="v_fCa00" variable_2="v_fCa00"/>
<map_variables variable_1="v_fCa01" variable_2="v_fCa01"/>
<map_variables variable_1="v_fCa02" variable_2="v_fCa02"/>
<map_variables variable_1="v_fCa10" variable_2="v_fCa10"/>
<map_variables variable_1="v_fCa11" variable_2="v_fCa11"/>
<map_variables variable_1="v_fCa12" variable_2="v_fCa12"/>
</connection>
<group>
<relationship_ref relationship="encapsulation"/>
<component_ref component="BG_LCC">
<component_ref component="LCC_states_BG"/>
</component_ref>
</group>
</model>