def model FCU_EC_coupling as
 def import using "units_and_constants/units_BG.cellml" for
            unit mM using unit mM;
unit fmol using unit fmol;
unit per_fmol using unit per_fmol;
            unit J_per_mol using unit J_per_mol;
unit fmol_per_sec using unit fmol_per_sec;
            unit C_per_mol using unit C_per_mol;
  unit J_per_C using unit J_per_C;
            unit microm3 using unit microm3;
  unit fF using unit fF;
            unit fC using unit fC;
  unit fA using unit fA;
            unit per_second using unit per_second;
  unit millivolt using unit millivolt;
            unit per_sec using unit per_sec;
  unit J_per_K_per_mol using unit J_per_K_per_mol;
            unit fmol_per_L using unit fmol_per_L;
  unit fmol_per_L_per_sec using unit fmol_per_L_per_sec;
            unit per_sec_per_fmol_per_L using unit per_sec_per_fmol_per_L;
  unit uM using unit uM;
            unit mM_per_sec using unit mM_per_sec;
  unit uM_per_sec using unit uM_per_sec;
            unit pL using unit pL;
  unit m_to_u using unit m_to_u;
 enddef;
def import using "units_and_constants/constants_BG.cellml" for
                comp constants using comp constants;
enddef;

def import using "Ca_leak/BG_Ca_leak.cellml" for
comp Ca_leak using comp Ca_leak;
enddef;
def import using "crossbridge_TRPN/BG_crossbridge_TRPN.cellml" for
comp crossbridge_TRPN using comp crossbridge_TRPN;
enddef;
def import using "diadCaDiffusion/BG_diadCaDiffusion.cellml" for
comp diadCaDiffusion using comp diadCaDiffusion;
enddef;
def import using "LCC_diadic/BG_LCC_diadic.cellml" for
comp LCC_diadic using comp LCC_diadic;
enddef;
def import using "NCX/BG_NCX.cellml" for
comp NCX using comp NCX;
enddef;
def import using "pCa/BG_pCa.cellml" for
comp pCa using comp pCa;
enddef;
def import using "PLB/BG_PLB.cellml" for
comp PLB using comp PLB;
enddef;
def import using "RyR/BG_RyR.cellml" for
comp RyR using comp RyR;
enddef;
def import using "SERCA/BG_SERCA.cellml" for
comp SERCA using comp SERCA;
enddef;

def comp BG_parameters as
var kappa_leak: fmol_per_sec {init: 1.59936e-07, pub: out};
var kappa_R_TRPNCa: fmol_per_sec {init: 62.0526, pub: out};
var kappa_R_BU: fmol_per_sec {init: 3.44525, pub: out};
var kappa_R_UW: fmol_per_sec {init: 0.324557, pub: out};
var kappa_R_WS: fmol_per_sec {init: 0.00763664, pub: out};
var kappa_R_SU: fmol_per_sec {init: 137459, pub: out};
var kappa_diadDf: fmol_per_sec {init: 0.0429033, pub: out};
var kappa_LCC_Ca1: fmol_per_sec {init: 2.42576, pub: out};
var kappa_LCC_Ca2: fmol_per_sec {init: 3.19386, pub: out};
var kappa_d000: fmol_per_sec {init: 211.625, pub: out};
var kappa_d010: fmol_per_sec {init: 10.687, pub: out};
var kappa_d020: fmol_per_sec {init: 14.071, pub: out};
var kappa_d001: fmol_per_sec {init: 0.123928, pub: out};
var kappa_d011: fmol_per_sec {init: 0.00625836, pub: out};
var kappa_d021: fmol_per_sec {init: 0.00824002, pub: out};
var kappa_f1_000: fmol_per_sec {init: 20.3846, pub: out};
var kappa_f1_100: fmol_per_sec {init: 101.225, pub: out};
var kappa_f1_001: fmol_per_sec {init: 0.0119373, pub: out};
var kappa_f1_101: fmol_per_sec {init: 0.0592778, pub: out};
var kappa_f2_000: fmol_per_sec {init: 1.32493, pub: out};
var kappa_f2_100: fmol_per_sec {init: 6.57925, pub: out};
var kappa_f2_001: fmol_per_sec {init: 0.000775884, pub: out};
var kappa_f2_101: fmol_per_sec {init: 0.00385284, pub: out};
var kappa_f3_010: fmol_per_sec {init: 2890.73, pub: out};
var kappa_f3_110: fmol_per_sec {init: 14354.6, pub: out};
var kappa_f3_011: fmol_per_sec {init: 1.69282, pub: out};
var kappa_f3_111: fmol_per_sec {init: 8.40614, pub: out};
var kappa_fCa000: fmol_per_sec {init: 42432.9, pub: out};
var kappa_fCa100: fmol_per_sec {init: 210711, pub: out};
var kappa_fCa010: fmol_per_sec {init: 2142.86, pub: out};
var kappa_fCa110: fmol_per_sec {init: 10640.9, pub: out};
var kappa_fCa020: fmol_per_sec {init: 2821.37, pub: out};
var kappa_fCa120: fmol_per_sec {init: 14010.2, pub: out};
var kappa_1_NCX: fmol_per_sec {init: 135753, pub: out};
var kappa_2_NCX: fmol_per_sec {init: 5.3773e+07, pub: out};
var kappa_3_NCX: fmol_per_sec {init: 19913.7, pub: out};
var kappa_4_NCX: fmol_per_sec {init: 5.3773e+07, pub: out};
var kappa_5_NCX: fmol_per_sec {init: 135753, pub: out};
var kappa_6_NCX: fmol_per_sec {init: 2.56034, pub: out};
var kappa_pCa_R1: fmol_per_sec {init: 11542.5, pub: out};
var kappa_pCa_R2: fmol_per_sec {init: 1.51449e-05, pub: out};
var kappa_PLBph1: fmol_per_sec {init: 45.5263, pub: out};
var kappa_PLBph2: fmol_per_sec {init: 6.55904, pub: out};
var kappa_PLBd1: fmol_per_sec {init: 0.386674, pub: out};
var kappa_PLBd2: fmol_per_sec {init: 1.21269, pub: out};
var kappa_Inh: fmol_per_sec {init: 431.435, pub: out};
var kappa_RyR: fmol_per_sec {init: 1.51729e+06, pub: out};
var kappa_OC: fmol_per_sec {init: 8.83262, pub: out};
var kappa_CCI: fmol_per_sec {init: 0.12618, pub: out};
var kappa_CII: fmol_per_sec {init: 883.262, pub: out};
var kappa_IO: fmol_per_sec {init: 73.6052, pub: out};
var kappa_SERCA_R1_2: fmol_per_sec {init: 0.000238983, pub: out};
var kappa_SERCA_R2_4: fmol_per_sec {init: 8043.9, pub: out};
var kappa_SERCA_R2_2a: fmol_per_sec {init: 9.60791e+06, pub: out};
var kappa_SERCA_R4_5: fmol_per_sec {init: 8043.9, pub: out};
var kappa_SERCA_R5_6: fmol_per_sec {init: 2400.02, pub: out};
var kappa_SERCA_R6_8: fmol_per_sec {init: 2.35087e+06, pub: out};
var kappa_SERCA_R8_9: fmol_per_sec {init: 5.66119e+06, pub: out};
var kappa_SERCA_R9_10: fmol_per_sec {init: 5.66119e+06, pub: out};
var kappa_SERCA_R10_1: fmol_per_sec {init: 0.00951056, pub: out};
var K_Ca_SR: per_fmol {init: 3.505, pub: out};
var K_Ca_i: per_fmol {init: 0.830832, pub: out};
var K_TRPN: per_fmol {init: 2.90497, pub: out};
var K_Ca_TRPN: per_fmol {init: 0.124732, pub: out};
var K_B_XB: per_fmol {init: 8.2672, pub: out};
var K_U_XB: per_fmol {init: 3.1002, pub: out};
var K_W_XB: per_fmol {init: 20.2705, pub: out};
var K_S_XB: per_fmol {init: 5.06764e-06, pub: out};
var K_Ca_D: per_fmol {init: 49.4198, pub: out};
var K_Ca_o: per_fmol {init: 0.0197641, pub: out};
var K_000_LCC: per_fmol {init: 2.30012, pub: out};
var K_010_LCC: per_fmol {init: 45.5471, pub: out};
var K_020_LCC: per_fmol {init: 34.5934, pub: out};
var K_100_LCC: per_fmol {init: 0.463197, pub: out};
var K_110_LCC: per_fmol {init: 9.17225, pub: out};
var K_120_LCC: per_fmol {init: 6.9664, pub: out};
var K_001_LCC: per_fmol {init: 3927.77, pub: out};
var K_011_LCC: per_fmol {init: 77777.9, pub: out};
var K_021_LCC: per_fmol {init: 59072.9, pub: out};
var K_101_LCC: per_fmol {init: 790.972, pub: out};
var K_111_LCC: per_fmol {init: 15662.9, pub: out};
var K_121_LCC: per_fmol {init: 11896.1, pub: out};
var K_1_NCX: per_fmol {init: 5.34633e+08, pub: out};
var K_2_NCX: per_fmol {init: 38.284, pub: out};
var K_3_NCX: per_fmol {init: 23987, pub: out};
var K_4_NCX: per_fmol {init: 43815.8, pub: out};
var K_5_NCX: per_fmol {init: 9.13027e+08, pub: out};
var K_6_NCX: per_fmol {init: 138.33, pub: out};
var K_Na_i: per_fmol {init: 0.0480675, pub: out};
var K_Na_o: per_fmol {init: 0.0403493, pub: out};
var K_pCa: per_fmol {init: 0.0158017, pub: out};
var K_pCa_Ca: per_fmol {init: 0.0165683, pub: out};
var K_PLB: per_fmol {init: 0.00313835, pub: out};
var K_PKACI: per_fmol {init: 0.526105, pub: out};
var K_PLB_PKACI: per_fmol {init: 0.850051, pub: out};
var K_PP1: per_fmol {init: 0.481905, pub: out};
var K_PLBp_PP1: per_fmol {init: 0.271255, pub: out};
var K_PLBp: per_fmol {init: 0.434012, pub: out};
var K_Ip: per_fmol {init: 7.20341, pub: out};
var K_Ip_PP1: per_fmol {init: 0.0896999, pub: out};
var K_C_RyR: per_fmol {init: 3962.58, pub: out};
var K_CI_RyR: per_fmol {init: 39.6258, pub: out};
var K_I_RyR: per_fmol {init: 0.06793, pub: out};
var K_O_RyR: per_fmol {init: 6.793, pub: out};
var K_P1_SERCA: per_fmol {init: 11929.9, pub: out};
var K_P2_SERCA: per_fmol {init: 8436.04, pub: out};
var K_P2a_SERCA: per_fmol {init: 2695.7, pub: out};
var K_P4_SERCA: per_fmol {init: 3.21983e+06, pub: out};
var K_P5_SERCA: per_fmol {init: 1.05832, pub: out};
var K_P6_SERCA: per_fmol {init: 55279.8, pub: out};
var K_P8_SERCA: per_fmol {init: 506.98, pub: out};
var K_P9_SERCA: per_fmol {init: 4575.01, pub: out};
var K_P10_SERCA: per_fmol {init: 2155.5, pub: out};
var K_H: per_fmol {init: 1078.99, pub: out};
var K_MgATP: per_fmol {init: 351.562, pub: out};
var K_MgADP: per_fmol {init: 1.96018e-05, pub: out};
var K_Pi: per_fmol {init: 0.0508218, pub: out};
enddef;
    def comp environment as
                var time: second {pub: out};
                var vol_myo: pL {init: 34.4, pub: out};
                var freq: dimensionless {init: 500};
var q_Ca_SR: fmol {init: 1e-888, pub: out};
var q_Ca_i: fmol {init: 1e-888, pub: out};
var q_TRPN: fmol {init: 1e-888, pub: out};
var q_Ca_TRPN: fmol {init: 1e-888, pub: out};
var q_B_XB: fmol {init: 1e-888, pub: out};
var q_U_XB: fmol {init: 1e-888, pub: out};
var q_W_XB: fmol {init: 1e-888, pub: out};
var q_S_XB: fmol {init: 1e-888, pub: out};
var q_Ca_D: fmol {init: 1e-888, pub: out};
var q_Ca_o: fmol {init: 1e-888, pub: out};
var q_000_LCC: fmol {init: 1e-888, pub: out};
var q_010_LCC: fmol {init: 1e-888, pub: out};
var q_020_LCC: fmol {init: 1e-888, pub: out};
var q_100_LCC: fmol {init: 1e-888, pub: out};
var q_110_LCC: fmol {init: 1e-888, pub: out};
var q_120_LCC: fmol {init: 1e-888, pub: out};
var q_001_LCC: fmol {init: 1e-888, pub: out};
var q_011_LCC: fmol {init: 1e-888, pub: out};
var q_021_LCC: fmol {init: 1e-888, pub: out};
var q_101_LCC: fmol {init: 1e-888, pub: out};
var q_111_LCC: fmol {init: 1e-888, pub: out};
var q_121_LCC: fmol {init: 1e-888, pub: out};
var q_1_NCX: fmol {init: 1e-888, pub: out};
var q_2_NCX: fmol {init: 1e-888, pub: out};
var q_3_NCX: fmol {init: 1e-888, pub: out};
var q_4_NCX: fmol {init: 1e-888, pub: out};
var q_5_NCX: fmol {init: 1e-888, pub: out};
var q_6_NCX: fmol {init: 1e-888, pub: out};
var q_Na_i: fmol {init: 1e-888, pub: out};
var q_Na_o: fmol {init: 1e-888, pub: out};
var q_pCa: fmol {init: 1e-888, pub: out};
var q_pCa_Ca: fmol {init: 1e-888, pub: out};
var q_PLB: fmol {init: 1e-888, pub: out};
var q_PKACI: fmol {init: 1e-888, pub: out};
var q_PLB_PKACI: fmol {init: 1e-888, pub: out};
var q_PP1: fmol {init: 1e-888, pub: out};
var q_PLBp_PP1: fmol {init: 1e-888, pub: out};
var q_PLBp: fmol {init: 1e-888, pub: out};
var q_Ip: fmol {init: 1e-888, pub: out};
var q_Ip_PP1: fmol {init: 1e-888, pub: out};
var q_C_RyR: fmol {init: 1e-888, pub: out};
var q_CI_RyR: fmol {init: 1e-888, pub: out};
var q_I_RyR: fmol {init: 1e-888, pub: out};
var q_O_RyR: fmol {init: 1e-888, pub: out};
var q_P1_SERCA: fmol {init: 1e-888, pub: out};
var q_P2_SERCA: fmol {init: 1e-888, pub: out};
var q_P2a_SERCA: fmol {init: 1e-888, pub: out};
var q_P4_SERCA: fmol {init: 1e-888, pub: out};
var q_P5_SERCA: fmol {init: 1e-888, pub: out};
var q_P6_SERCA: fmol {init: 1e-888, pub: out};
var q_P8_SERCA: fmol {init: 1e-888, pub: out};
var q_P9_SERCA: fmol {init: 1e-888, pub: out};
var q_P10_SERCA: fmol {init: 1e-888, pub: out};
var q_H: fmol {init: 1e-888, pub: out};
var q_MgATP: fmol {init: 1e-888, pub: out};
var q_MgADP: fmol {init: 1e-888, pub: out};
var q_Pi: fmol {init: 1e-888, pub: out};

// Ca_leak imports
var v_leak: fmol_per_sec {pub: in};


// crossbridge_TRPN imports
var v_R_TRPNCa: fmol_per_sec {pub: in};
var v_R_BU: fmol_per_sec {pub: in};
var v_R_UW: fmol_per_sec {pub: in};
var v_R_WS: fmol_per_sec {pub: in};
var v_R_SU: fmol_per_sec {pub: in};


// diadCaDiffusion imports
var v_diadDf: fmol_per_sec {pub: in};


// LCC_diadic imports
var v_LCC_Ca1: fmol_per_sec {pub: in};
var v_LCC_Ca2: fmol_per_sec {pub: in};
var v_d000: fmol_per_sec {pub: in};
var v_d010: fmol_per_sec {pub: in};
var v_d020: fmol_per_sec {pub: in};
var v_d001: fmol_per_sec {pub: in};
var v_d011: fmol_per_sec {pub: in};
var v_d021: fmol_per_sec {pub: in};
var v_f1_000: fmol_per_sec {pub: in};
var v_f1_100: fmol_per_sec {pub: in};
var v_f1_001: fmol_per_sec {pub: in};
var v_f1_101: fmol_per_sec {pub: in};
var v_f2_000: fmol_per_sec {pub: in};
var v_f2_100: fmol_per_sec {pub: in};
var v_f2_001: fmol_per_sec {pub: in};
var v_f2_101: fmol_per_sec {pub: in};
var v_f3_010: fmol_per_sec {pub: in};
var v_f3_110: fmol_per_sec {pub: in};
var v_f3_011: fmol_per_sec {pub: in};
var v_f3_111: fmol_per_sec {pub: in};
var v_fCa000: fmol_per_sec {pub: in};
var v_fCa100: fmol_per_sec {pub: in};
var v_fCa010: fmol_per_sec {pub: in};
var v_fCa110: fmol_per_sec {pub: in};
var v_fCa020: fmol_per_sec {pub: in};
var v_fCa120: fmol_per_sec {pub: in};


// NCX imports
var v_1_NCX: fmol_per_sec {pub: in};
var v_2_NCX: fmol_per_sec {pub: in};
var v_3_NCX: fmol_per_sec {pub: in};
var v_4_NCX: fmol_per_sec {pub: in};
var v_5_NCX: fmol_per_sec {pub: in};
var v_6_NCX: fmol_per_sec {pub: in};


// pCa imports
var v_pCa_R1: fmol_per_sec {pub: in};
var v_pCa_R2: fmol_per_sec {pub: in};


// PLB imports
var v_PLBph1: fmol_per_sec {pub: in};
var v_PLBph2: fmol_per_sec {pub: in};
var v_PLBd1: fmol_per_sec {pub: in};
var v_PLBd2: fmol_per_sec {pub: in};
var v_Inh: fmol_per_sec {pub: in};


// RyR imports
var v_RyR: fmol_per_sec {pub: in};
var v_OC: fmol_per_sec {pub: in};
var v_CCI: fmol_per_sec {pub: in};
var v_CII: fmol_per_sec {pub: in};
var v_IO: fmol_per_sec {pub: in};


// SERCA imports
var v_SERCA_R1_2: fmol_per_sec {pub: in};
var v_SERCA_R2_4: fmol_per_sec {pub: in};
var v_SERCA_R2_2a: fmol_per_sec {pub: in};
var v_SERCA_R4_5: fmol_per_sec {pub: in};
var v_SERCA_R5_6: fmol_per_sec {pub: in};
var v_SERCA_R6_8: fmol_per_sec {pub: in};
var v_SERCA_R8_9: fmol_per_sec {pub: in};
var v_SERCA_R9_10: fmol_per_sec {pub: in};
var v_SERCA_R10_1: fmol_per_sec {pub: in};


ode(q_Ca_SR, time) = + v_mCa_leak  + v_mRyR  + v_mSERCA ;
ode(q_Ca_i, time) = + v_mCa_leak  + v_mcrossbridge_TRPN  + v_mdiadCaDiffusion  + v_mNCX  + v_mpCa  + v_mSERCA ;
ode(q_TRPN, time) = + v_mcrossbridge_TRPN ;
ode(q_Ca_TRPN, time) = + v_mcrossbridge_TRPN ;
ode(q_B_XB, time) = + v_mcrossbridge_TRPN ;
ode(q_U_XB, time) = + v_mcrossbridge_TRPN ;
ode(q_W_XB, time) = + v_mcrossbridge_TRPN ;
ode(q_S_XB, time) = + v_mcrossbridge_TRPN ;
ode(q_Ca_D, time) = + v_mdiadCaDiffusion  + v_mLCC_diadic  + v_mRyR ;
ode(q_Ca_o, time) = + v_mLCC_diadic  + v_mNCX  + v_mpCa ;
ode(q_000_LCC, time) = + v_mLCC_diadic ;
ode(q_010_LCC, time) = + v_mLCC_diadic ;
ode(q_020_LCC, time) = + v_mLCC_diadic ;
ode(q_100_LCC, time) = + v_mLCC_diadic ;
ode(q_110_LCC, time) = + v_mLCC_diadic ;
ode(q_120_LCC, time) = + v_mLCC_diadic ;
ode(q_001_LCC, time) = + v_mLCC_diadic ;
ode(q_011_LCC, time) = + v_mLCC_diadic ;
ode(q_021_LCC, time) = + v_mLCC_diadic ;
ode(q_101_LCC, time) = + v_mLCC_diadic ;
ode(q_111_LCC, time) = + v_mLCC_diadic ;
ode(q_121_LCC, time) = + v_mLCC_diadic ;
ode(q_1_NCX, time) = + v_mNCX ;
ode(q_2_NCX, time) = + v_mNCX ;
ode(q_3_NCX, time) = + v_mNCX ;
ode(q_4_NCX, time) = + v_mNCX ;
ode(q_5_NCX, time) = + v_mNCX ;
ode(q_6_NCX, time) = + v_mNCX ;
ode(q_Na_i, time) = + v_mNCX ;
ode(q_Na_o, time) = + v_mNCX ;
ode(q_pCa, time) = + v_mpCa ;
ode(q_pCa_Ca, time) = + v_mpCa ;
ode(q_PLB, time) = + v_mPLB ;
ode(q_PKACI, time) = + v_mPLB ;
ode(q_PLB_PKACI, time) = + v_mPLB ;
ode(q_PP1, time) = + v_mPLB ;
ode(q_PLBp_PP1, time) = + v_mPLB ;
ode(q_PLBp, time) = + v_mPLB ;
ode(q_Ip, time) = + v_mPLB ;
ode(q_Ip_PP1, time) = + v_mPLB ;
ode(q_C_RyR, time) = + v_mRyR ;
ode(q_CI_RyR, time) = + v_mRyR ;
ode(q_I_RyR, time) = + v_mRyR ;
ode(q_O_RyR, time) = + v_mRyR ;
ode(q_P1_SERCA, time) = + v_mSERCA ;
ode(q_P2_SERCA, time) = + v_mSERCA ;
ode(q_P2a_SERCA, time) = + v_mSERCA ;
ode(q_P4_SERCA, time) = + v_mSERCA ;
ode(q_P5_SERCA, time) = + v_mSERCA ;
ode(q_P6_SERCA, time) = + v_mSERCA ;
ode(q_P8_SERCA, time) = + v_mSERCA ;
ode(q_P9_SERCA, time) = + v_mSERCA ;
ode(q_P10_SERCA, time) = + v_mSERCA ;
ode(q_H, time) = + v_mSERCA ;
ode(q_MgATP, time) = + v_mSERCA ;
ode(q_MgADP, time) = + v_mSERCA ;
ode(q_Pi, time) = + v_mSERCA ;
enddef;

def map between environment and Ca_leak for
vars time and time;
vars q_Ca_SR and q_Ca_SR;
vars q_Ca_i and q_Ca_i;
vars v_leak and v_leak;
enddef;

def map between environment and crossbridge_TRPN for
vars time and time;
vars q_TRPN and q_TRPN;
vars q_Ca_i and q_Ca_i;
vars q_Ca_TRPN and q_Ca_TRPN;
vars q_B_XB and q_B_XB;
vars q_U_XB and q_U_XB;
vars q_W_XB and q_W_XB;
vars q_S_XB and q_S_XB;
vars v_R_TRPNCa and v_R_TRPNCa;
vars v_R_BU and v_R_BU;
vars v_R_UW and v_R_UW;
vars v_R_WS and v_R_WS;
vars v_R_SU and v_R_SU;
enddef;

def map between environment and diadCaDiffusion for
vars time and time;
vars q_Ca_D and q_Ca_D;
vars q_Ca_i and q_Ca_i;
vars v_diadDf and v_diadDf;
enddef;

def map between environment and LCC_diadic for
vars time and time;
vars q_Ca_D and q_Ca_D;
vars q_Ca_o and q_Ca_o;
vars q_000_LCC and q_000_LCC;
vars q_010_LCC and q_010_LCC;
vars q_020_LCC and q_020_LCC;
vars q_100_LCC and q_100_LCC;
vars q_110_LCC and q_110_LCC;
vars q_120_LCC and q_120_LCC;
vars q_001_LCC and q_001_LCC;
vars q_011_LCC and q_011_LCC;
vars q_021_LCC and q_021_LCC;
vars q_101_LCC and q_101_LCC;
vars q_111_LCC and q_111_LCC;
vars q_121_LCC and q_121_LCC;
vars v_LCC_Ca1 and v_LCC_Ca1;
vars v_LCC_Ca2 and v_LCC_Ca2;
vars v_d000 and v_d000;
vars v_d010 and v_d010;
vars v_d020 and v_d020;
vars v_d001 and v_d001;
vars v_d011 and v_d011;
vars v_d021 and v_d021;
vars v_f1_000 and v_f1_000;
vars v_f1_100 and v_f1_100;
vars v_f1_001 and v_f1_001;
vars v_f1_101 and v_f1_101;
vars v_f2_000 and v_f2_000;
vars v_f2_100 and v_f2_100;
vars v_f2_001 and v_f2_001;
vars v_f2_101 and v_f2_101;
vars v_f3_010 and v_f3_010;
vars v_f3_110 and v_f3_110;
vars v_f3_011 and v_f3_011;
vars v_f3_111 and v_f3_111;
vars v_fCa000 and v_fCa000;
vars v_fCa100 and v_fCa100;
vars v_fCa010 and v_fCa010;
vars v_fCa110 and v_fCa110;
vars v_fCa020 and v_fCa020;
vars v_fCa120 and v_fCa120;
enddef;

def map between environment and NCX for
vars time and time;
vars q_1_NCX and q_1_NCX;
vars q_2_NCX and q_2_NCX;
vars q_3_NCX and q_3_NCX;
vars q_4_NCX and q_4_NCX;
vars q_5_NCX and q_5_NCX;
vars q_6_NCX and q_6_NCX;
vars q_Na_i and q_Na_i;
vars q_Ca_i and q_Ca_i;
vars q_Na_o and q_Na_o;
vars q_Ca_o and q_Ca_o;
vars v_1_NCX and v_1_NCX;
vars v_2_NCX and v_2_NCX;
vars v_3_NCX and v_3_NCX;
vars v_4_NCX and v_4_NCX;
vars v_5_NCX and v_5_NCX;
vars v_6_NCX and v_6_NCX;
enddef;

def map between environment and pCa for
vars time and time;
vars q_Ca_i and q_Ca_i;
vars q_Ca_o and q_Ca_o;
vars q_pCa and q_pCa;
vars q_pCa_Ca and q_pCa_Ca;
vars v_pCa_R1 and v_pCa_R1;
vars v_pCa_R2 and v_pCa_R2;
enddef;

def map between environment and PLB for
vars time and time;
vars q_PLB and q_PLB;
vars q_PKACI and q_PKACI;
vars q_PLB_PKACI and q_PLB_PKACI;
vars q_PP1 and q_PP1;
vars q_PLBp_PP1 and q_PLBp_PP1;
vars q_PLBp and q_PLBp;
vars q_Ip and q_Ip;
vars q_Ip_PP1 and q_Ip_PP1;
vars v_PLBph1 and v_PLBph1;
vars v_PLBph2 and v_PLBph2;
vars v_PLBd1 and v_PLBd1;
vars v_PLBd2 and v_PLBd2;
vars v_Inh and v_Inh;
enddef;

def map between environment and RyR for
vars time and time;
vars q_Ca_SR and q_Ca_SR;
vars q_Ca_D and q_Ca_D;
vars q_C_RyR and q_C_RyR;
vars q_CI_RyR and q_CI_RyR;
vars q_I_RyR and q_I_RyR;
vars q_O_RyR and q_O_RyR;
vars v_RyR and v_RyR;
vars v_OC and v_OC;
vars v_CCI and v_CCI;
vars v_CII and v_CII;
vars v_IO and v_IO;
enddef;

def map between environment and SERCA for
vars time and time;
vars q_P1_SERCA and q_P1_SERCA;
vars q_P2_SERCA and q_P2_SERCA;
vars q_P2a_SERCA and q_P2a_SERCA;
vars q_P4_SERCA and q_P4_SERCA;
vars q_P5_SERCA and q_P5_SERCA;
vars q_P6_SERCA and q_P6_SERCA;
vars q_P8_SERCA and q_P8_SERCA;
vars q_P9_SERCA and q_P9_SERCA;
vars q_P10_SERCA and q_P10_SERCA;
vars q_H and q_H;
vars q_Ca_i and q_Ca_i;
vars q_Ca_SR and q_Ca_SR;
vars q_MgATP and q_MgATP;
vars q_MgADP and q_MgADP;
vars q_Pi and q_Pi;
vars v_SERCA_R1_2 and v_SERCA_R1_2;
vars v_SERCA_R2_4 and v_SERCA_R2_4;
vars v_SERCA_R2_2a and v_SERCA_R2_2a;
vars v_SERCA_R4_5 and v_SERCA_R4_5;
vars v_SERCA_R5_6 and v_SERCA_R5_6;
vars v_SERCA_R6_8 and v_SERCA_R6_8;
vars v_SERCA_R8_9 and v_SERCA_R8_9;
vars v_SERCA_R9_10 and v_SERCA_R9_10;
vars v_SERCA_R10_1 and v_SERCA_R10_1;
enddef;

def map between BG_parameters and Ca_leak for
vars kappa_leak and kappa_leak;
vars K_Ca_SR and K_Ca_SR;
vars K_Ca_i and K_Ca_i;
enddef;
def map between BG_parameters and crossbridge_TRPN for
vars kappa_R_TRPNCa and kappa_R_TRPNCa;
vars kappa_R_BU and kappa_R_BU;
vars kappa_R_UW and kappa_R_UW;
vars kappa_R_WS and kappa_R_WS;
vars kappa_R_SU and kappa_R_SU;
vars K_TRPN and K_TRPN;
vars K_Ca_i and K_Ca_i;
vars K_Ca_TRPN and K_Ca_TRPN;
vars K_B_XB and K_B_XB;
vars K_U_XB and K_U_XB;
vars K_W_XB and K_W_XB;
vars K_S_XB and K_S_XB;
enddef;
def map between BG_parameters and diadCaDiffusion for
vars kappa_diadDf and kappa_diadDf;
vars K_Ca_D and K_Ca_D;
vars K_Ca_i and K_Ca_i;
enddef;
def map between BG_parameters and LCC_diadic for
vars kappa_LCC_Ca1 and kappa_LCC_Ca1;
vars kappa_LCC_Ca2 and kappa_LCC_Ca2;
vars kappa_d000 and kappa_d000;
vars kappa_d010 and kappa_d010;
vars kappa_d020 and kappa_d020;
vars kappa_d001 and kappa_d001;
vars kappa_d011 and kappa_d011;
vars kappa_d021 and kappa_d021;
vars kappa_f1_000 and kappa_f1_000;
vars kappa_f1_100 and kappa_f1_100;
vars kappa_f1_001 and kappa_f1_001;
vars kappa_f1_101 and kappa_f1_101;
vars kappa_f2_000 and kappa_f2_000;
vars kappa_f2_100 and kappa_f2_100;
vars kappa_f2_001 and kappa_f2_001;
vars kappa_f2_101 and kappa_f2_101;
vars kappa_f3_010 and kappa_f3_010;
vars kappa_f3_110 and kappa_f3_110;
vars kappa_f3_011 and kappa_f3_011;
vars kappa_f3_111 and kappa_f3_111;
vars kappa_fCa000 and kappa_fCa000;
vars kappa_fCa100 and kappa_fCa100;
vars kappa_fCa010 and kappa_fCa010;
vars kappa_fCa110 and kappa_fCa110;
vars kappa_fCa020 and kappa_fCa020;
vars kappa_fCa120 and kappa_fCa120;
vars K_Ca_D and K_Ca_D;
vars K_Ca_o and K_Ca_o;
vars K_000_LCC and K_000_LCC;
vars K_010_LCC and K_010_LCC;
vars K_020_LCC and K_020_LCC;
vars K_100_LCC and K_100_LCC;
vars K_110_LCC and K_110_LCC;
vars K_120_LCC and K_120_LCC;
vars K_001_LCC and K_001_LCC;
vars K_011_LCC and K_011_LCC;
vars K_021_LCC and K_021_LCC;
vars K_101_LCC and K_101_LCC;
vars K_111_LCC and K_111_LCC;
vars K_121_LCC and K_121_LCC;
enddef;
def map between BG_parameters and NCX for
vars kappa_1_NCX and kappa_1_NCX;
vars kappa_2_NCX and kappa_2_NCX;
vars kappa_3_NCX and kappa_3_NCX;
vars kappa_4_NCX and kappa_4_NCX;
vars kappa_5_NCX and kappa_5_NCX;
vars kappa_6_NCX and kappa_6_NCX;
vars K_1_NCX and K_1_NCX;
vars K_2_NCX and K_2_NCX;
vars K_3_NCX and K_3_NCX;
vars K_4_NCX and K_4_NCX;
vars K_5_NCX and K_5_NCX;
vars K_6_NCX and K_6_NCX;
vars K_Na_i and K_Na_i;
vars K_Ca_i and K_Ca_i;
vars K_Na_o and K_Na_o;
vars K_Ca_o and K_Ca_o;
enddef;
def map between BG_parameters and pCa for
vars kappa_pCa_R1 and kappa_pCa_R1;
vars kappa_pCa_R2 and kappa_pCa_R2;
vars K_Ca_i and K_Ca_i;
vars K_Ca_o and K_Ca_o;
vars K_pCa and K_pCa;
vars K_pCa_Ca and K_pCa_Ca;
enddef;
def map between BG_parameters and PLB for
vars kappa_PLBph1 and kappa_PLBph1;
vars kappa_PLBph2 and kappa_PLBph2;
vars kappa_PLBd1 and kappa_PLBd1;
vars kappa_PLBd2 and kappa_PLBd2;
vars kappa_Inh and kappa_Inh;
vars K_PLB and K_PLB;
vars K_PKACI and K_PKACI;
vars K_PLB_PKACI and K_PLB_PKACI;
vars K_PP1 and K_PP1;
vars K_PLBp_PP1 and K_PLBp_PP1;
vars K_PLBp and K_PLBp;
vars K_Ip and K_Ip;
vars K_Ip_PP1 and K_Ip_PP1;
enddef;
def map between BG_parameters and RyR for
vars kappa_RyR and kappa_RyR;
vars kappa_OC and kappa_OC;
vars kappa_CCI and kappa_CCI;
vars kappa_CII and kappa_CII;
vars kappa_IO and kappa_IO;
vars K_Ca_SR and K_Ca_SR;
vars K_Ca_D and K_Ca_D;
vars K_C_RyR and K_C_RyR;
vars K_CI_RyR and K_CI_RyR;
vars K_I_RyR and K_I_RyR;
vars K_O_RyR and K_O_RyR;
enddef;
def map between BG_parameters and SERCA for
vars kappa_SERCA_R1_2 and kappa_SERCA_R1_2;
vars kappa_SERCA_R2_4 and kappa_SERCA_R2_4;
vars kappa_SERCA_R2_2a and kappa_SERCA_R2_2a;
vars kappa_SERCA_R4_5 and kappa_SERCA_R4_5;
vars kappa_SERCA_R5_6 and kappa_SERCA_R5_6;
vars kappa_SERCA_R6_8 and kappa_SERCA_R6_8;
vars kappa_SERCA_R8_9 and kappa_SERCA_R8_9;
vars kappa_SERCA_R9_10 and kappa_SERCA_R9_10;
vars kappa_SERCA_R10_1 and kappa_SERCA_R10_1;
vars K_P1_SERCA and K_P1_SERCA;
vars K_P2_SERCA and K_P2_SERCA;
vars K_P2a_SERCA and K_P2a_SERCA;
vars K_P4_SERCA and K_P4_SERCA;
vars K_P5_SERCA and K_P5_SERCA;
vars K_P6_SERCA and K_P6_SERCA;
vars K_P8_SERCA and K_P8_SERCA;
vars K_P9_SERCA and K_P9_SERCA;
vars K_P10_SERCA and K_P10_SERCA;
vars K_H and K_H;
vars K_Ca_i and K_Ca_i;
vars K_Ca_SR and K_Ca_SR;
vars K_MgATP and K_MgATP;
vars K_MgADP and K_MgADP;
vars K_Pi and K_Pi;
enddef;

def map between constants and Ca_leak for
	vars R and R;
	vars T and T;
enddef;
def map between constants and crossbridge_TRPN for
	vars R and R;
	vars T and T;
enddef;
def map between constants and diadCaDiffusion for
	vars R and R;
	vars T and T;
enddef;
def map between constants and LCC_diadic for
	vars R and R;
	vars T and T;
enddef;
def map between constants and NCX for
	vars R and R;
	vars T and T;
enddef;
def map between constants and pCa for
	vars R and R;
	vars T and T;
enddef;
def map between constants and PLB for
	vars R and R;
	vars T and T;
enddef;
def map between constants and RyR for
	vars R and R;
	vars T and T;
enddef;
def map between constants and SERCA for
	vars R and R;
	vars T and T;
enddef;

enddef;