Location: FCU_ExcitationContraction_Coupling @ 0a35f86cc268 / FCU_EC_coupling.cellml

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2022-07-11 10:30:52+12:00
Desc:
With updated XB module and no more rootfinding
Permanent Source URI:
https://models.cellml.org/workspace/7a8/rawfile/0a35f86cc26838a052c3248e94224d5b384d89a3/FCU_EC_coupling.cellml

<?xml version='1.0' encoding='UTF-8'?>
<model name="FCU_EC_coupling" 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"/>
        <units name="per_sec" units_ref="per_sec"/>
        <units name="J_per_K_per_mol" units_ref="J_per_K_per_mol"/>
        <units name="fmol_per_L" units_ref="fmol_per_L"/>
        <units name="fmol_per_L_per_sec" units_ref="fmol_per_L_per_sec"/>
        <units name="per_sec_per_fmol_per_L" units_ref="per_sec_per_fmol_per_L"/>
        <units name="uM" units_ref="uM"/>
        <units name="mM_per_sec" units_ref="mM_per_sec"/>
        <units name="uM_per_sec" units_ref="uM_per_sec"/>
        <units name="pL" units_ref="pL"/>
        <units name="m_to_u" units_ref="m_to_u"/>
    </import>
    <import xlink:href="units_and_constants/constants_BG.cellml">
        <component component_ref="constants" name="constants"/>
    </import>
    <import xlink:href="units_and_constants/ion_valency.cellml">
        <component component_ref="ion_valences" name="ion_valences"/>
    </import>
    <import xlink:href="Ca_leak/BG_Ca_leak.cellml">
        <component component_ref="Ca_leak" name="Ca_leak"/>
    </import>
    <import xlink:href="crossbridge_TRPN/BG_crossbridge_TRPN.cellml">
        <component component_ref="crossbridge_TRPN" name="crossbridge_TRPN"/>
    </import>
    <import xlink:href="diadCaDiffusion/BG_diadCaDiffusion.cellml">
        <component component_ref="diadCaDiffusion" name="diadCaDiffusion"/>
    </import>
    <import xlink:href="LCC_diadic/BG_LCC_diadic.cellml">
        <component component_ref="LCC" name="LCC_diadic"/>
    </import>
    <import xlink:href="NCX/BG_NCX.cellml">
        <component component_ref="NCX" name="NCX"/>
    </import>
    <import xlink:href="pCa/BG_pCa.cellml">
        <component component_ref="pCa" name="pCa"/>
    </import>
    <import xlink:href="PLB/BG_PLB.cellml">
        <component component_ref="PLB" name="PLB"/>
    </import>
    <import xlink:href="RyR/BG_RyR.cellml">
        <component component_ref="RyR" name="RyR"/>
    </import>
    <import xlink:href="SERCA/BG_SERCA.cellml">
        <component component_ref="SERCA" name="SERCA"/>
    </import>
    <units name="m_per_s">
        <unit units="metre"/>
        <unit exponent="-1" units="second"/>
    </units>
    <component name="BG_parameters">
        <variable initial_value="1.19278e-07" name="kappa_leak" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="61.975" name="kappa_R_TRPNCa" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.70563" name="kappa_R_BU" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.295103" name="kappa_R_UW" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00694361" name="kappa_R_WS" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="124985" name="kappa_R_SU" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0319966" name="kappa_diadDf" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.45746" name="kappa_LCC_Ca1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.91895" name="kappa_LCC_Ca2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="306.532" name="kappa_d000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="15.4798" name="kappa_d010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="20.3814" name="kappa_d020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0998404" name="kappa_d001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00504193" name="kappa_d011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00663841" name="kappa_d021" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="29.5266" name="kappa_f1_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="146.621" name="kappa_f1_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00961708" name="kappa_f1_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.047756" name="kappa_f1_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.91912" name="kappa_f2_000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="9.52987" name="kappa_f2_100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.000625075" name="kappa_f2_001" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.00310397" name="kappa_f2_101" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="4187.14" name="kappa_f3_010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="20792.3" name="kappa_f3_110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.36379" name="kappa_f3_011" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.77224" name="kappa_f3_111" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="34185.2" name="kappa_fCa000" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="169755" name="kappa_fCa100" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1726.35" name="kappa_fCa010" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8572.61" name="kappa_fCa110" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="2272.98" name="kappa_fCa020" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="11287.1" name="kappa_fCa120" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="429216" name="kappa_1_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.06923e+08" name="kappa_2_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3622.96" name="kappa_3_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.06923e+08" name="kappa_4_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="429216" name="kappa_5_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.380974" name="kappa_6_NCX" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="10726.4" name="kappa_pCa_R1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.40741e-05" name="kappa_pCa_R2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="12.9276" name="kappa_PLBph1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.86249" name="kappa_PLBph2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.0640148" name="kappa_PLBd1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.200763" name="kappa_PLBd2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1269.27" name="kappa_Inh" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="1.13157e+06" name="kappa_RyR" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="8.50173" name="kappa_OC" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.164843" name="kappa_CCI" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="626.393" name="kappa_CII" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="52.1994" name="kappa_IO" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="3.59662e-07" name="kappa_SERCA_R1_2" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.73317" name="kappa_SERCA_R2_4" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="32063.4" name="kappa_SERCA_R2_2a" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6.73317" name="kappa_SERCA_R4_5" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="0.408566" name="kappa_SERCA_R5_6" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="313.54" name="kappa_SERCA_R6_8" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6675" name="kappa_SERCA_R8_9" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="6675" name="kappa_SERCA_R9_10" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="5.05704e-06" name="kappa_SERCA_R10_1" public_interface="out" units="fmol_per_sec"/>
        <variable initial_value="4.69975" name="K_Ca_SR" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.11404" name="K_Ca_i" public_interface="out" units="per_fmol"/>
        <variable initial_value="12.0566" name="K_TRPN" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.062444" name="K_Ca_TRPN" public_interface="out" units="per_fmol"/>
        <variable initial_value="9.09233" name="K_B_XB" public_interface="out" units="per_fmol"/>
        <variable initial_value="3.40963" name="K_U_XB" public_interface="out" units="per_fmol"/>
        <variable initial_value="22.2937" name="K_W_XB" public_interface="out" units="per_fmol"/>
        <variable initial_value="5.57343e-06" name="K_S_XB" public_interface="out" units="per_fmol"/>
        <variable initial_value="66.2656" name="K_Ca_D" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0265011" name="K_Ca_o" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.58796" name="K_000_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="31.4449" name="K_010_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="23.8827" name="K_020_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.319783" name="K_100_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="6.33236" name="K_110_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="4.80948" name="K_120_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="4875.4" name="K_001_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="96542.9" name="K_011_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="73325.1" name="K_021_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="981.805" name="K_101_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="19441.8" name="K_111_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="14766.2" name="K_121_LCC" public_interface="out" units="per_fmol"/>
        <variable initial_value="2.85254e+08" name="K_1_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="16.9249" name="K_2_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="34379.9" name="K_3_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="313010" name="K_4_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="6.01862e+07" name="K_5_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="929.686" name="K_6_NCX" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.012842" name="K_Na_i" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0107799" name="K_Na_o" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0126813" name="K_pCa" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0178289" name="K_pCa_Ca" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.482706" name="K_PLB" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0120458" name="K_PKACI" public_interface="out" units="per_fmol"/>
        <variable initial_value="2.99358" name="K_PLB_PKACI" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0189254" name="K_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.63848" name="K_PLBp_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="66.755" name="K_PLBp" public_interface="out" units="per_fmol"/>
        <variable initial_value="62.347" name="K_Ip" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0304897" name="K_Ip_PP1" public_interface="out" units="per_fmol"/>
        <variable initial_value="2234.81" name="K_C_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="30.332" name="K_CI_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.0957866" name="K_I_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="7.05739" name="K_O_RyR" public_interface="out" units="per_fmol"/>
        <variable initial_value="62.115" name="K_P1_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.88391" name="K_P2_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.271481" name="K_P2a_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="1292.79" name="K_P4_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00208939" name="K_P5_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="139.3" name="K_P6_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.710561" name="K_P8_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.30406" name="K_P9_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.3624" name="K_P10_SERCA" public_interface="out" units="per_fmol"/>
        <variable initial_value="486.59" name="K_H" public_interface="out" units="per_fmol"/>
        <variable initial_value="15.0787" name="K_MgATP" public_interface="out" units="per_fmol"/>
        <variable initial_value="1.53572e-05" name="K_MgADP" public_interface="out" units="per_fmol"/>
        <variable initial_value="0.00616948" name="K_Pi" public_interface="out" units="per_fmol"/>
        <variable initial_value="2" name="z_rCa" public_interface="out" units="dimensionless"/>
        <variable initial_value="2.1404" name="z_fd" public_interface="out" units="dimensionless"/>
        <variable initial_value="-1.1495" name="z_ff1" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.72162" name="z_ff2" public_interface="out" units="dimensionless"/>
        <variable initial_value="4.2933" name="z_ff3" public_interface="out" units="dimensionless"/>
        <variable initial_value="-2.1404" name="z_rd" public_interface="out" units="dimensionless"/>
        <variable initial_value="1.8993" name="z_rf1" public_interface="out" units="dimensionless"/>
        <variable initial_value="-0.52288" name="z_rf2" public_interface="out" units="dimensionless"/>
        <variable initial_value="0" name="z_rf3" public_interface="out" units="dimensionless"/>
        <variable initial_value="3" name="nNa_o_NCX" public_interface="out" units="dimensionless"/>
        <variable initial_value="3" name="nNa_i_NCX" public_interface="out" units="dimensionless"/>
        <variable initial_value="-0.55679" name="zf_NCX" public_interface="out" units="dimensionless"/>
        <variable initial_value="0.44321" name="zr_NCX" public_interface="out" units="dimensionless"/>
        <variable initial_value="1.1" name="n_Tm" public_interface="out" units="dimensionless"/>
        <variable initial_value="2" name="n_CaTRPN" public_interface="out" units="dimensionless"/>
    </component>
    <component name="environment">
        <variable name="time" public_interface="out" units="second"/>
        <variable initial_value="25.84" name="vol_myo" public_interface="out" units="pL"/>
        <variable initial_value="153400" name="C_m" public_interface="out" units="fF"/>
        <!--        var F: C_per_mol {pub: in};-->
        <variable name="n_CaTRPN" public_interface="in" units="dimensionless"/>
        <variable initial_value="-13039" name="q_membrane" public_interface="out" units="fC"/>
        <variable initial_value="3.68876E-07" name="q_Ca_D" public_interface="out" units="fmol"/>
        <variable initial_value="2.57" name="q_TRPN" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_Ca_TRPN" public_interface="out" units="fmol"/>
        <variable initial_value="0.005663912" name="q_Ca_i" public_interface="out" units="fmol"/>
        <variable initial_value="481.46" name="q_Na_i" public_interface="out" units="fmol"/>
        <variable initial_value="725.48" name="q_Na_o" public_interface="out" units="fmol"/>
        <variable initial_value="1" name="q_B_XB" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_U_XB" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_W_XB" public_interface="out" units="fmol"/>
        <variable initial_value="0" name="q_S_XB" public_interface="out" units="fmol"/>
        <variable initial_value="6.84" name="q_Ca_o" public_interface="out" units="fmol"/>
        <variable initial_value="4.028E+00" name="q_PLB" public_interface="out" units="fmol"/>
        <variable initial_value="2.234E-03" name="q_PKACI" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLB_PKACI" public_interface="out" units="fmol"/>
        <variable initial_value="3.382E-02" name="q_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLBp_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLBp" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_PLBp_baseline" public_interface="out" units="fmol"/>
        <variable initial_value="1.999E-03" name="q_Ip" public_interface="out" units="fmol"/>
        <variable initial_value="1e-18" name="q_Ip_PP1" public_interface="out" units="fmol"/>
        <variable initial_value="0.641" name="q_Ca_SR" public_interface="out" units="fmol"/>
        <variable initial_value="0.0032" name="q_pCa" public_interface="out" units="fmol"/>
        <variable initial_value="1e-9" name="q_pCa_Ca" public_interface="out" units="fmol"/>
        <variable initial_value="2.584e-6" name="q_H" public_interface="out" units="fmol"/>
        <variable initial_value="570" name="q_Pi" public_interface="out" units="fmol"/>
        <variable initial_value="1.3794" name="q_MgADP" public_interface="out" units="fmol"/>
        <variable initial_value="3.8" name="q_MgATP" public_interface="out" units="fmol"/>
        <!--     gate species-->
        <variable initial_value="1.3798e-05" name="q_000_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_010_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_020_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_100_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_110_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_120_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_001_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_011_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_021_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_101_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_111_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="1.3798e-05" name="q_121_LCC" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P1_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P2_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P3_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P4_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P5_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.24554e-07" name="q_P6_NCX" public_interface="out" units="fmol"/>
        <variable initial_value="4.13941e-05" name="q_C_RyR" public_interface="out" units="fmol"/>
        <variable initial_value="4.13941e-05" name="q_CI_RyR" public_interface="out" units="fmol"/>
        <variable initial_value="4.13941e-05" name="q_I_RyR" public_interface="out" units="fmol"/>
        <variable initial_value="4.13941e-05" name="q_O_RyR" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P1_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P2_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P2a_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P4_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P5_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P6_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P8_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P9_SERCA" public_interface="out" units="fmol"/>
        <variable initial_value="2.83036e-08" name="q_P10_SERCA" public_interface="out" units="fmol"/>
        <!-- Ca_leak imports-->
        <variable name="v_leak" public_interface="in" units="fmol_per_sec"/>
        <!-- crossbridge_TRPN imports-->
        <variable name="v_R_TRPNCa" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_BU" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_UW" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_WS" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_R_SU" public_interface="in" units="fmol_per_sec"/>
        <variable name="Tension" public_interface="in" units="kPa"/>
        <variable name="SL" public_interface="in" units="metre"/>
        <!-- diadCaDiffusion imports-->
        <variable name="v_diadDf" public_interface="in" units="fmol_per_sec"/>
        <!-- LCC_diadic imports-->
        <variable name="v_Ca_D_LCC" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Ca_o_LCC" public_interface="in" units="fmol_per_sec"/>
        <variable name="I_mem_LCC" public_interface="in" units="fA"/>
        <!-- NCX imports-->
        <variable name="v_Ca_o_NCX" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Ca_i_NCX" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Na_o_NCX" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Na_i_NCX" public_interface="in" units="fmol_per_sec"/>
        <variable name="I_mem_NCX" public_interface="in" units="fA"/>
        <!-- pCa imports-->
        <variable name="v_pCa_R1" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_pCa_R2" public_interface="in" units="fmol_per_sec"/>
        <variable name="I_mem_pCa" public_interface="in" units="fA"/>
        <!-- PLB imports-->
        <variable name="v_PLBph1" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBph2" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBd1" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_PLBd2" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Inh" public_interface="in" units="fmol_per_sec"/>
        <!-- RyR imports-->
        <variable name="v_RyR" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_RyRgate_Ca_D" public_interface="in" units="fmol_per_sec"/>
        <!-- SERCA imports-->
        <variable name="v_Ca_i_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Ca_SR_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_H_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_MgADP_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_MgATP_SERCA" public_interface="in" units="fmol_per_sec"/>
        <variable name="v_Pi_SERCA" public_interface="in" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_SR</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <minus/>
                            <ci>v_leak</ci>
                        </apply>
                        <ci>v_RyR</ci>
                    </apply>
                    <ci>v_Ca_SR_SERCA</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_i</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <apply>
                            <plus/>
                            <apply>
                                <times/>
                                <apply>
                                    <minus/>
                                    <ci>v_R_TRPNCa</ci>
                                </apply>
                                <ci>n_CaTRPN</ci>
                            </apply>
                            <ci>v_Ca_i_SERCA</ci>
                            <ci>v_leak</ci>
                            <ci>v_diadDf</ci>
                            <ci>v_Ca_i_NCX</ci>
                        </apply>
                        <ci>v_pCa_R1</ci>
                    </apply>
                    <ci>v_Ca_stim</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_TRPN</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_TRPNCa</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_TRPN</ci>
                </apply>
                <ci>v_R_TRPNCa</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_B_XB</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_BU</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_U_XB</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_R_BU</ci>
                        <ci>v_R_UW</ci>
                    </apply>
                    <ci>v_R_SU</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_W_XB</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_UW</ci>
                    <ci>v_R_WS</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_S_XB</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_R_WS</ci>
                    <ci>v_R_SU</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_D</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_diadDf</ci>
                    </apply>
                    <ci>v_Ca_D_LCC</ci>
                    <ci>v_RyR</ci>
                    <ci>v_RyRgate_Ca_D</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ca_o</ci>
                </apply>
                <apply>
                    <plus/>
                    <ci>v_Ca_o_LCC</ci>
                    <ci>v_Ca_o_NCX</ci>
                    <ci>v_pCa_R2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Na_i</ci>
                </apply>
                <ci>v_Na_i_NCX</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Na_o</ci>
                </apply>
                <ci>v_Na_o_NCX</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_pCa</ci>
                </apply>
                <apply>
                    <plus/>
                    <apply>
                        <minus/>
                        <ci>v_pCa_R1</ci>
                    </apply>
                    <ci>v_pCa_R2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_pCa_Ca</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_pCa_R1</ci>
                    <ci>v_pCa_R2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLB</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBd2</ci>
                    <ci>v_PLBph1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PKACI</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph2</ci>
                    <ci>v_PLBph1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLB_PKACI</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph1</ci>
                    <ci>v_PLBph2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PP1</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph1</ci>
                    <ci>v_PLBph2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLBp_PP1</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBd1</ci>
                    <ci>v_PLBd2</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_PLBp</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_PLBph2</ci>
                    <ci>v_PLBd1</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ip</ci>
                </apply>
                <apply>
                    <minus/>
                    <ci>v_Inh</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Ip_PP1</ci>
                </apply>
                <ci>v_Inh</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_H</ci>
                </apply>
                <ci>v_H_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_MgADP</ci>
                </apply>
                <ci>v_MgADP_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_MgATP</ci>
                </apply>
                <ci>v_MgATP_SERCA</ci>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_Pi</ci>
                </apply>
                <ci>v_Pi_SERCA</ci>
            </apply>
        </math>
        <!-- membrane potential change-->
        <variable name="V_m" public_interface="out" units="volt"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>V_m</ci>
                <apply>
                    <divide/>
                    <ci>q_membrane</ci>
                    <ci>C_m</ci>
                </apply>
            </apply>
        </math>
        <variable name="I_stim" units="fA"/>
        <variable initial_value="1" name="stimPeriod" units="second"/>
        <variable initial_value="0.001" name="stimDuration" units="second"/>
        <variable name="tPeriod" units="second"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>tPeriod</ci>
                <apply>
                    <minus/>
                    <ci>time</ci>
                    <apply>
                        <times/>
                        <apply>
                            <floor/>
                            <apply>
                                <divide/>
                                <ci>time</ci>
                                <ci>stimPeriod</ci>
                            </apply>
                        </apply>
                        <ci>stimPeriod</ci>
                    </apply>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>I_stim</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <divide/>
                            <apply>
                                <times/>
                                <cn cellml:units="dimensionless">0</cn>
                                <cn cellml:units="volt">0.03</cn>
                                <ci>C_m</ci>
                            </apply>
                            <ci>stimDuration</ci>
                        </apply>
                        <apply>
                            <and/>
                            <apply>
                                <geq/>
                                <ci>tPeriod</ci>
                                <cn cellml:units="second">0.3</cn>
                            </apply>
                            <apply>
                                <leq/>
                                <ci>tPeriod</ci>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="second">0.3</cn>
                                    <ci>stimDuration</ci>
                                </apply>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <cn cellml:units="fA">0.0</cn>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
        <variable name="v_Ca_stim" units="fmol_per_sec"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>v_Ca_stim</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0</cn>
                            <cn cellml:units="fmol_per_sec" type="e-notation">8<sep/>3</cn>
                        </apply>
                        <apply>
                            <and/>
                            <apply>
                                <geq/>
                                <ci>tPeriod</ci>
                                <cn cellml:units="second">0.3</cn>
                            </apply>
                            <apply>
                                <leq/>
                                <ci>tPeriod</ci>
                                <apply>
                                    <plus/>
                                    <cn cellml:units="second">0.3</cn>
                                    <ci>stimDuration</ci>
                                </apply>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <cn cellml:units="fmol_per_sec">0.0</cn>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
        <variable name="sum_I" units="fA"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>sum_I</ci>
                <apply>
                    <plus/>
                    <ci>I_mem_LCC</ci>
                    <ci>I_mem_NCX</ci>
                    <ci>I_mem_pCa</ci>
                    <ci>I_stim</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <apply>
                    <diff/>
                    <bvar>
                        <ci>time</ci>
                    </bvar>
                    <ci>q_membrane</ci>
                </apply>
                <ci>sum_I</ci>
            </apply>
        </math>
        <!-- changing sarcomere length-->
        <variable name="dSL" public_interface="out" units="m_per_s"/>
        <variable initial_value="19e-6" name="V_max" units="m_per_s"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>dSL</ci>
                <piecewise>
                    <piece>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">1</cn>
                            <ci>V_max</ci>
                        </apply>
                        <apply>
                            <lt/>
                            <ci>time</ci>
                            <cn cellml:units="second">0.06</cn>
                        </apply>
                    </piece>
                    <piece>
                        <apply>
                            <times/>
                            <apply>
                                <minus/>
                                <cn cellml:units="dimensionless">0.1</cn>
                            </apply>
                            <ci>V_max</ci>
                        </apply>
                        <apply>
                            <and/>
                            <apply>
                                <gt/>
                                <ci>time</ci>
                                <cn cellml:units="second">0.5</cn>
                            </apply>
                            <apply>
                                <lt/>
                                <ci>time</ci>
                                <cn cellml:units="second">1.1</cn>
                            </apply>
                        </apply>
                    </piece>
                    <otherwise>
                        <apply>
                            <times/>
                            <cn cellml:units="dimensionless">0</cn>
                            <ci>V_max</ci>
                        </apply>
                    </otherwise>
                </piecewise>
            </apply>
        </math>
        <!-- concentrations-->
        <variable initial_value="38" name="vol_cell" units="pL"/>
        <variable name="vol_SR" units="pL"/>
        <variable initial_value="5.182" name="vol_o" units="pL"/>
        <variable initial_value="0.0025375" name="vol_D" units="pL"/>
        <variable name="cCai" units="mM"/>
        <variable name="cCao" units="mM"/>
        <variable name="cCaSR" units="mM"/>
        <variable name="cCaD" units="mM"/>
        <variable name="cNai" units="mM"/>
        <variable name="pHi" units="dimensionless"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>vol_SR</ci>
                <apply>
                    <times/>
                    <ci>vol_cell</ci>
                    <cn cellml:units="dimensionless">0.035</cn>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>cCai</ci>
                <apply>
                    <divide/>
                    <ci>q_Ca_i</ci>
                    <ci>vol_myo</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>cCao</ci>
                <apply>
                    <divide/>
                    <ci>q_Ca_o</ci>
                    <ci>vol_o</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>cCaD</ci>
                <apply>
                    <divide/>
                    <ci>q_Ca_D</ci>
                    <ci>vol_D</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>cCaSR</ci>
                <apply>
                    <divide/>
                    <ci>q_Ca_SR</ci>
                    <ci>vol_SR</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>cNai</ci>
                <apply>
                    <divide/>
                    <ci>q_Na_i</ci>
                    <ci>vol_myo</ci>
                </apply>
            </apply>
            <apply>
                <eq/>
                <ci>pHi</ci>
                <apply>
                    <minus/>
                    <apply>
                        <log/>
                        <apply>
                            <divide/>
                            <ci>q_H</ci>
                            <ci>vol_myo</ci>
                        </apply>
                    </apply>
                </apply>
            </apply>
        </math>
        <!-- conservation of Ca-->
        <variable name="Ca_tot" units="fmol"/>
        <math xmlns="http://www.w3.org/1998/Math/MathML">
            <apply>
                <eq/>
                <ci>Ca_tot</ci>
                <apply>
                    <plus/>
                    <ci>q_Ca_i</ci>
                    <ci>q_Ca_o</ci>
                    <ci>q_Ca_TRPN</ci>
                    <ci>q_Ca_SR</ci>
                    <ci>q_Ca_D</ci>
                    <ci>q_pCa_Ca</ci>
                    <apply>
                        <times/>
                        <cn cellml:units="dimensionless">2</cn>
                        <apply>
                            <plus/>
                            <ci>q_P4_SERCA</ci>
                            <ci>q_P5_SERCA</ci>
                            <ci>q_P6_SERCA</ci>
                        </apply>
                    </apply>
                    <ci>q_P3_NCX</ci>
                    <ci>q_P4_NCX</ci>
                </apply>
            </apply>
        </math>
    </component>
    <connection>
        <map_components component_1="environment" component_2="Ca_leak"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="v_leak" variable_2="v_leak"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_TRPN" variable_2="q_TRPN"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="q_Ca_TRPN" variable_2="q_Ca_TRPN"/>
        <map_variables variable_1="q_B_XB" variable_2="q_B_XB"/>
        <map_variables variable_1="q_U_XB" variable_2="q_U_XB"/>
        <map_variables variable_1="q_W_XB" variable_2="q_W_XB"/>
        <map_variables variable_1="q_S_XB" variable_2="q_S_XB"/>
        <map_variables variable_1="v_R_TRPNCa" variable_2="v_R_TRPNCa"/>
        <map_variables variable_1="v_R_BU" variable_2="v_R_BU"/>
        <map_variables variable_1="v_R_UW" variable_2="v_R_UW"/>
        <map_variables variable_1="v_R_WS" variable_2="v_R_WS"/>
        <map_variables variable_1="v_R_SU" variable_2="v_R_SU"/>
        <map_variables variable_1="Tension" variable_2="T_total"/>
        <map_variables variable_1="SL" variable_2="SL"/>
        <map_variables variable_1="dSL" variable_2="dSL"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="diadCaDiffusion"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_Ca_D" variable_2="q_Ca_D"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="v_diadDf" variable_2="v_diadDf"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="LCC_diadic"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="q_Ca_D" variable_2="q_Ca_D"/>
        <map_variables variable_1="q_Ca_o" variable_2="q_Ca_o"/>
        <map_variables variable_1="q_membrane" variable_2="q_mem"/>
        <map_variables variable_1="C_m" variable_2="C_m"/>
        <map_variables variable_1="v_Ca_D_LCC" variable_2="v_Ca_D_LCC"/>
        <map_variables variable_1="v_Ca_o_LCC" variable_2="v_Ca_o_LCC"/>
        <map_variables variable_1="I_mem_LCC" variable_2="I_mem_LCC"/>
        <map_variables variable_1="q_000_LCC" variable_2="q_S000_LCC"/>
        <map_variables variable_1="q_010_LCC" variable_2="q_S010_LCC"/>
        <map_variables variable_1="q_020_LCC" variable_2="q_S020_LCC"/>
        <map_variables variable_1="q_100_LCC" variable_2="q_S100_LCC"/>
        <map_variables variable_1="q_110_LCC" variable_2="q_S110_LCC"/>
        <map_variables variable_1="q_120_LCC" variable_2="q_S120_LCC"/>
        <map_variables variable_1="q_001_LCC" variable_2="q_S001_LCC"/>
        <map_variables variable_1="q_011_LCC" variable_2="q_S011_LCC"/>
        <map_variables variable_1="q_021_LCC" variable_2="q_S021_LCC"/>
        <map_variables variable_1="q_101_LCC" variable_2="q_S101_LCC"/>
        <map_variables variable_1="q_111_LCC" variable_2="q_S111_LCC"/>
        <map_variables variable_1="q_121_LCC" variable_2="q_S121_LCC"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="NCX"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="C_m" variable_2="C_m"/>
        <map_variables variable_1="q_Na_i" variable_2="q_Na_i"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="q_Na_o" variable_2="q_Na_o"/>
        <map_variables variable_1="q_Ca_o" variable_2="q_Ca_o"/>
        <map_variables variable_1="v_Ca_o_NCX" variable_2="v_Ca_o_NCX"/>
        <map_variables variable_1="v_Ca_i_NCX" variable_2="v_Ca_i_NCX"/>
        <map_variables variable_1="v_Na_o_NCX" variable_2="v_Na_o_NCX"/>
        <map_variables variable_1="v_Na_i_NCX" variable_2="v_Na_i_NCX"/>
        <map_variables variable_1="q_P1_NCX" variable_2="q_P1_NCX"/>
        <map_variables variable_1="q_P2_NCX" variable_2="q_P2_NCX"/>
        <map_variables variable_1="q_P3_NCX" variable_2="q_P3_NCX"/>
        <map_variables variable_1="q_P4_NCX" variable_2="q_P4_NCX"/>
        <map_variables variable_1="q_P5_NCX" variable_2="q_P5_NCX"/>
        <map_variables variable_1="q_P6_NCX" variable_2="q_P6_NCX"/>
        <map_variables variable_1="q_membrane" variable_2="q_mem"/>
        <map_variables variable_1="I_mem_NCX" variable_2="I_mem_NCX"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="pCa"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="q_Ca_o" variable_2="q_Ca_o"/>
        <map_variables variable_1="q_pCa" variable_2="q_pCa"/>
        <map_variables variable_1="q_pCa_Ca" variable_2="q_pCa_Ca"/>
        <map_variables variable_1="v_pCa_R1" variable_2="v_pCa_R1"/>
        <map_variables variable_1="v_pCa_R2" variable_2="v_pCa_R2"/>
        <map_variables variable_1="I_mem_pCa" variable_2="I_mem_pCa"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="PLB"/>
        <map_variables variable_1="time" variable_2="time"/>
        <map_variables variable_1="q_PLB" variable_2="q_PLB"/>
        <map_variables variable_1="q_PKACI" variable_2="q_PKACI"/>
        <map_variables variable_1="q_PLB_PKACI" variable_2="q_PLB_PKACI"/>
        <map_variables variable_1="q_PP1" variable_2="q_PP1"/>
        <map_variables variable_1="q_PLBp_PP1" variable_2="q_PLBp_PP1"/>
        <map_variables variable_1="q_PLBp" variable_2="q_PLBp"/>
        <map_variables variable_1="q_Ip" variable_2="q_Ip"/>
        <map_variables variable_1="q_Ip_PP1" variable_2="q_Ip_PP1"/>
        <map_variables variable_1="v_PLBph1" variable_2="v_PLBph1"/>
        <map_variables variable_1="v_PLBph2" variable_2="v_PLBph2"/>
        <map_variables variable_1="v_PLBd1" variable_2="v_PLBd1"/>
        <map_variables variable_1="v_PLBd2" variable_2="v_PLBd2"/>
        <map_variables variable_1="v_Inh" variable_2="v_Inh"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="RyR"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
        <map_variables variable_1="q_Ca_D" variable_2="q_Ca_D"/>
        <map_variables variable_1="v_RyR" variable_2="v_RyR"/>
        <map_variables variable_1="v_RyRgate_Ca_D" variable_2="v_RyRgate_Ca_D"/>
        <map_variables variable_1="q_C_RyR" variable_2="q_C_RyR"/>
        <map_variables variable_1="q_O_RyR" variable_2="q_O_RyR"/>
        <map_variables variable_1="q_CI_RyR" variable_2="q_CI_RyR"/>
        <map_variables variable_1="q_I_RyR" variable_2="q_I_RyR"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="SERCA"/>
        <map_variables variable_1="time" variable_2="t"/>
        <map_variables variable_1="q_H" variable_2="q_H"/>
        <map_variables variable_1="q_Ca_i" variable_2="q_Ca_i"/>
        <map_variables variable_1="q_Ca_SR" variable_2="q_Ca_SR"/>
        <map_variables variable_1="q_MgATP" variable_2="q_MgATP"/>
        <map_variables variable_1="q_MgADP" variable_2="q_MgADP"/>
        <map_variables variable_1="q_Pi" variable_2="q_Pi"/>
        <map_variables variable_1="q_PLBp" variable_2="q_PLBp"/>
        <map_variables variable_1="q_PLBp_baseline" variable_2="q_PLBp_baseline"/>
        <map_variables variable_1="v_Ca_i_SERCA" variable_2="v_Ca_i_SERCA"/>
        <map_variables variable_1="v_Ca_SR_SERCA" variable_2="v_Ca_SR_SERCA"/>
        <map_variables variable_1="v_H_SERCA" variable_2="v_H_SERCA"/>
        <map_variables variable_1="v_MgADP_SERCA" variable_2="v_MgADP_SERCA"/>
        <map_variables variable_1="v_MgATP_SERCA" variable_2="v_MgATP_SERCA"/>
        <map_variables variable_1="v_Pi_SERCA" variable_2="v_Pi_SERCA"/>
        <map_variables variable_1="q_P1_SERCA" variable_2="q_P1_SERCA"/>
        <map_variables variable_1="q_P2_SERCA" variable_2="q_P2_SERCA"/>
        <map_variables variable_1="q_P2a_SERCA" variable_2="q_P2a_SERCA"/>
        <map_variables variable_1="q_P4_SERCA" variable_2="q_P4_SERCA"/>
        <map_variables variable_1="q_P5_SERCA" variable_2="q_P5_SERCA"/>
        <map_variables variable_1="q_P6_SERCA" variable_2="q_P6_SERCA"/>
        <map_variables variable_1="q_P8_SERCA" variable_2="q_P8_SERCA"/>
        <map_variables variable_1="q_P9_SERCA" variable_2="q_P9_SERCA"/>
        <map_variables variable_1="q_P10_SERCA" variable_2="q_P10_SERCA"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="Ca_leak"/>
        <map_variables variable_1="kappa_leak" variable_2="kappa_leak"/>
        <map_variables variable_1="K_Ca_SR" variable_2="K_Ca_SR"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="environment" component_2="BG_parameters"/>
        <map_variables variable_1="n_CaTRPN" variable_2="n_CaTRPN"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="kappa_R_TRPNCa" variable_2="kappa_R_TRPNCa"/>
        <map_variables variable_1="kappa_R_BU" variable_2="kappa_R_BU"/>
        <map_variables variable_1="kappa_R_UW" variable_2="kappa_R_UW"/>
        <map_variables variable_1="kappa_R_WS" variable_2="kappa_R_WS"/>
        <map_variables variable_1="kappa_R_SU" variable_2="kappa_R_SU"/>
        <map_variables variable_1="K_TRPN" variable_2="K_TRPN"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
        <map_variables variable_1="K_Ca_TRPN" variable_2="K_Ca_TRPN"/>
        <map_variables variable_1="K_B_XB" variable_2="K_B_XB"/>
        <map_variables variable_1="K_U_XB" variable_2="K_U_XB"/>
        <map_variables variable_1="K_W_XB" variable_2="K_W_XB"/>
        <map_variables variable_1="K_S_XB" variable_2="K_S_XB"/>
        <map_variables variable_1="n_Tm" variable_2="n_Tm"/>
        <map_variables variable_1="n_CaTRPN" variable_2="n_CaTRPN"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="diadCaDiffusion"/>
        <map_variables variable_1="kappa_diadDf" variable_2="kappa_diadDf"/>
        <map_variables variable_1="K_Ca_D" variable_2="K_Ca_D"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="LCC_diadic"/>
        <map_variables variable_1="kappa_LCC_Ca1" variable_2="kappa_LCC_Ca1"/>
        <map_variables variable_1="kappa_LCC_Ca2" variable_2="kappa_LCC_Ca2"/>
        <map_variables variable_1="kappa_d000" variable_2="kappa_d000"/>
        <map_variables variable_1="kappa_d010" variable_2="kappa_d010"/>
        <map_variables variable_1="kappa_d020" variable_2="kappa_d020"/>
        <map_variables variable_1="kappa_d001" variable_2="kappa_d001"/>
        <map_variables variable_1="kappa_d011" variable_2="kappa_d011"/>
        <map_variables variable_1="kappa_d021" variable_2="kappa_d021"/>
        <map_variables variable_1="kappa_f1_000" variable_2="kappa_f1_000"/>
        <map_variables variable_1="kappa_f1_100" variable_2="kappa_f1_100"/>
        <map_variables variable_1="kappa_f1_001" variable_2="kappa_f1_001"/>
        <map_variables variable_1="kappa_f1_101" variable_2="kappa_f1_101"/>
        <map_variables variable_1="kappa_f2_000" variable_2="kappa_f2_000"/>
        <map_variables variable_1="kappa_f2_100" variable_2="kappa_f2_100"/>
        <map_variables variable_1="kappa_f2_001" variable_2="kappa_f2_001"/>
        <map_variables variable_1="kappa_f2_101" variable_2="kappa_f2_101"/>
        <map_variables variable_1="kappa_f3_010" variable_2="kappa_f3_010"/>
        <map_variables variable_1="kappa_f3_110" variable_2="kappa_f3_110"/>
        <map_variables variable_1="kappa_f3_011" variable_2="kappa_f3_011"/>
        <map_variables variable_1="kappa_f3_111" variable_2="kappa_f3_111"/>
        <map_variables variable_1="kappa_fCa000" variable_2="kappa_fCa000"/>
        <map_variables variable_1="kappa_fCa100" variable_2="kappa_fCa100"/>
        <map_variables variable_1="kappa_fCa010" variable_2="kappa_fCa010"/>
        <map_variables variable_1="kappa_fCa110" variable_2="kappa_fCa110"/>
        <map_variables variable_1="kappa_fCa020" variable_2="kappa_fCa020"/>
        <map_variables variable_1="kappa_fCa120" variable_2="kappa_fCa120"/>
        <map_variables variable_1="K_Ca_D" variable_2="K_Ca_D"/>
        <map_variables variable_1="K_Ca_o" variable_2="K_Ca_o"/>
        <map_variables variable_1="K_000_LCC" variable_2="K_000_LCC"/>
        <map_variables variable_1="K_010_LCC" variable_2="K_010_LCC"/>
        <map_variables variable_1="K_020_LCC" variable_2="K_020_LCC"/>
        <map_variables variable_1="K_100_LCC" variable_2="K_100_LCC"/>
        <map_variables variable_1="K_110_LCC" variable_2="K_110_LCC"/>
        <map_variables variable_1="K_120_LCC" variable_2="K_120_LCC"/>
        <map_variables variable_1="K_001_LCC" variable_2="K_001_LCC"/>
        <map_variables variable_1="K_011_LCC" variable_2="K_011_LCC"/>
        <map_variables variable_1="K_021_LCC" variable_2="K_021_LCC"/>
        <map_variables variable_1="K_101_LCC" variable_2="K_101_LCC"/>
        <map_variables variable_1="K_111_LCC" variable_2="K_111_LCC"/>
        <map_variables variable_1="K_121_LCC" variable_2="K_121_LCC"/>
        <map_variables variable_1="z_rCa" variable_2="z_rCa"/>
        <map_variables variable_1="z_fd" variable_2="z_fd"/>
        <map_variables variable_1="z_ff1" variable_2="z_ff1"/>
        <map_variables variable_1="z_ff2" variable_2="z_ff2"/>
        <map_variables variable_1="z_ff3" variable_2="z_ff3"/>
        <map_variables variable_1="z_rd" variable_2="z_rd"/>
        <map_variables variable_1="z_rf1" variable_2="z_rf1"/>
        <map_variables variable_1="z_rf2" variable_2="z_rf2"/>
        <map_variables variable_1="z_rf3" variable_2="z_rf3"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="NCX"/>
        <map_variables variable_1="kappa_1_NCX" variable_2="kappa_1_NCX"/>
        <map_variables variable_1="kappa_2_NCX" variable_2="kappa_2_NCX"/>
        <map_variables variable_1="kappa_3_NCX" variable_2="kappa_3_NCX"/>
        <map_variables variable_1="kappa_4_NCX" variable_2="kappa_4_NCX"/>
        <map_variables variable_1="kappa_5_NCX" variable_2="kappa_5_NCX"/>
        <map_variables variable_1="kappa_6_NCX" variable_2="kappa_6_NCX"/>
        <map_variables variable_1="K_1_NCX" variable_2="K_1_NCX"/>
        <map_variables variable_1="K_2_NCX" variable_2="K_2_NCX"/>
        <map_variables variable_1="K_3_NCX" variable_2="K_3_NCX"/>
        <map_variables variable_1="K_4_NCX" variable_2="K_4_NCX"/>
        <map_variables variable_1="K_5_NCX" variable_2="K_5_NCX"/>
        <map_variables variable_1="K_6_NCX" variable_2="K_6_NCX"/>
        <map_variables variable_1="K_Na_i" variable_2="K_Na_i"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
        <map_variables variable_1="K_Na_o" variable_2="K_Na_o"/>
        <map_variables variable_1="K_Ca_o" variable_2="K_Ca_o"/>
        <map_variables variable_1="nNa_o_NCX" variable_2="nNa_o_NCX"/>
        <map_variables variable_1="nNa_i_NCX" variable_2="nNa_i_NCX"/>
        <map_variables variable_1="zf_NCX" variable_2="zf_NCX"/>
        <map_variables variable_1="zr_NCX" variable_2="zr_NCX"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="pCa"/>
        <map_variables variable_1="kappa_pCa_R1" variable_2="kappa_pCa_R1"/>
        <map_variables variable_1="kappa_pCa_R2" variable_2="kappa_pCa_R2"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
        <map_variables variable_1="K_Ca_o" variable_2="K_Ca_o"/>
        <map_variables variable_1="K_pCa" variable_2="K_pCa"/>
        <map_variables variable_1="K_pCa_Ca" variable_2="K_pCa_Ca"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="PLB"/>
        <map_variables variable_1="kappa_PLBph1" variable_2="kappa_PLBph1"/>
        <map_variables variable_1="kappa_PLBph2" variable_2="kappa_PLBph2"/>
        <map_variables variable_1="kappa_PLBd1" variable_2="kappa_PLBd1"/>
        <map_variables variable_1="kappa_PLBd2" variable_2="kappa_PLBd2"/>
        <map_variables variable_1="kappa_Inh" variable_2="kappa_Inh"/>
        <map_variables variable_1="K_PLB" variable_2="K_PLB"/>
        <map_variables variable_1="K_PKACI" variable_2="K_PKACI"/>
        <map_variables variable_1="K_PLB_PKACI" variable_2="K_PLB_PKACI"/>
        <map_variables variable_1="K_PP1" variable_2="K_PP1"/>
        <map_variables variable_1="K_PLBp_PP1" variable_2="K_PLBp_PP1"/>
        <map_variables variable_1="K_PLBp" variable_2="K_PLBp"/>
        <map_variables variable_1="K_Ip" variable_2="K_Ip"/>
        <map_variables variable_1="K_Ip_PP1" variable_2="K_Ip_PP1"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="RyR"/>
        <map_variables variable_1="kappa_RyR" variable_2="kappa_RyR"/>
        <map_variables variable_1="kappa_OC" variable_2="kappa_OC"/>
        <map_variables variable_1="kappa_CCI" variable_2="kappa_CCI"/>
        <map_variables variable_1="kappa_CII" variable_2="kappa_CII"/>
        <map_variables variable_1="kappa_IO" variable_2="kappa_IO"/>
        <map_variables variable_1="K_Ca_SR" variable_2="K_Ca_SR"/>
        <map_variables variable_1="K_Ca_D" variable_2="K_Ca_D"/>
        <map_variables variable_1="K_C_RyR" variable_2="K_C_RyR"/>
        <map_variables variable_1="K_CI_RyR" variable_2="K_CI_RyR"/>
        <map_variables variable_1="K_I_RyR" variable_2="K_I_RyR"/>
        <map_variables variable_1="K_O_RyR" variable_2="K_O_RyR"/>
    </connection>
    <connection>
        <map_components component_1="BG_parameters" component_2="SERCA"/>
        <map_variables variable_1="kappa_SERCA_R1_2" variable_2="kappa_SERCA_R1_2"/>
        <map_variables variable_1="kappa_SERCA_R2_4" variable_2="kappa_SERCA_R2_4"/>
        <map_variables variable_1="kappa_SERCA_R2_2a" variable_2="kappa_SERCA_R2_2a"/>
        <map_variables variable_1="kappa_SERCA_R4_5" variable_2="kappa_SERCA_R4_5"/>
        <map_variables variable_1="kappa_SERCA_R5_6" variable_2="kappa_SERCA_R5_6"/>
        <map_variables variable_1="kappa_SERCA_R6_8" variable_2="kappa_SERCA_R6_8"/>
        <map_variables variable_1="kappa_SERCA_R8_9" variable_2="kappa_SERCA_R8_9"/>
        <map_variables variable_1="kappa_SERCA_R9_10" variable_2="kappa_SERCA_R9_10"/>
        <map_variables variable_1="kappa_SERCA_R10_1" variable_2="kappa_SERCA_R10_1"/>
        <map_variables variable_1="K_P1_SERCA" variable_2="K_P1_SERCA"/>
        <map_variables variable_1="K_P2_SERCA" variable_2="K_P2_SERCA"/>
        <map_variables variable_1="K_P2a_SERCA" variable_2="K_P2a_SERCA"/>
        <map_variables variable_1="K_P4_SERCA" variable_2="K_P4_SERCA"/>
        <map_variables variable_1="K_P5_SERCA" variable_2="K_P5_SERCA"/>
        <map_variables variable_1="K_P6_SERCA" variable_2="K_P6_SERCA"/>
        <map_variables variable_1="K_P8_SERCA" variable_2="K_P8_SERCA"/>
        <map_variables variable_1="K_P9_SERCA" variable_2="K_P9_SERCA"/>
        <map_variables variable_1="K_P10_SERCA" variable_2="K_P10_SERCA"/>
        <map_variables variable_1="K_H" variable_2="K_H"/>
        <map_variables variable_1="K_Ca_i" variable_2="K_Ca_i"/>
        <map_variables variable_1="K_Ca_SR" variable_2="K_Ca_SR"/>
        <map_variables variable_1="K_MgATP" variable_2="K_MgATP"/>
        <map_variables variable_1="K_MgADP" variable_2="K_MgADP"/>
        <map_variables variable_1="K_Pi" variable_2="K_Pi"/>
        <map_variables variable_1="K_PLBp" variable_2="K_PLBp"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="Ca_leak"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="crossbridge_TRPN"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="diadCaDiffusion"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="LCC_diadic"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="NCX"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="pCa"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="PLB"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="RyR"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <connection>
        <map_components component_1="constants" component_2="SERCA"/>
        <map_variables variable_1="R" variable_2="R"/>
        <map_variables variable_1="T" variable_2="T"/>
        <map_variables variable_1="F" variable_2="F"/>
    </connection>
    <!--    def map between constants and environment for
        vars F and F;
    enddef;-->
    <connection>
        <map_components component_1="ion_valences" component_2="LCC_diadic"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
    </connection>
    <connection>
        <map_components component_1="ion_valences" component_2="pCa"/>
        <map_variables variable_1="zCa" variable_2="zCa"/>
    </connection>
</model>