def model FCU_calciumToolbox 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 "IPR/BG_IPR.cellml" for
comp IPR using comp IPR;
enddef;
def import using "pCa/BG_pCa.cellml" for
comp pCa using comp pCa;
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: 3.07502e-09, pub: out};
var kappa_R_IPR_main: fmol_per_sec {init: 2.9078e+08, pub: out};
var kappa_R1_a_IPR: fmol_per_sec {init: 0.0125066, pub: out};
var kappa_R2_a_IPR: fmol_per_sec {init: 16.6342, pub: out};
var kappa_R3_a_IPR: fmol_per_sec {init: 0.291916, pub: out};
var kappa_R4_a_IPR: fmol_per_sec {init: 1.25066, pub: out};
var kappa_R5_a_IPR: fmol_per_sec {init: 5.31398e-05, pub: out};
var kappa_R6_a_IPR: fmol_per_sec {init: 0.00973583, pub: out};
var kappa_R7_a_IPR: fmol_per_sec {init: 1.18239, pub: out};
var kappa_R8_a_IPR: fmol_per_sec {init: 4829.13, pub: out};
var kappa_R9_a_IPR: fmol_per_sec {init: 0.0955807, pub: out};
var kappa_R10_a_IPR: fmol_per_sec {init: 0.00718894, pub: out};
var kappa_R11_a_IPR: fmol_per_sec {init: 2.23094, pub: out};
var kappa_R12_a_IPR: fmol_per_sec {init: 152486, pub: out};
var kappa_R1_b_IPR: fmol_per_sec {init: 0.0125066, pub: out};
var kappa_R2_b_IPR: fmol_per_sec {init: 16.6342, pub: out};
var kappa_R3_b_IPR: fmol_per_sec {init: 0.291916, pub: out};
var kappa_R4_b_IPR: fmol_per_sec {init: 1.25066, pub: out};
var kappa_R5_b_IPR: fmol_per_sec {init: 5.31398e-05, pub: out};
var kappa_R6_b_IPR: fmol_per_sec {init: 0.00973583, pub: out};
var kappa_R7_b_IPR: fmol_per_sec {init: 1.18239, pub: out};
var kappa_R8_b_IPR: fmol_per_sec {init: 4829.13, pub: out};
var kappa_R9_b_IPR: fmol_per_sec {init: 0.0955807, pub: out};
var kappa_R10_b_IPR: fmol_per_sec {init: 0.00718894, pub: out};
var kappa_R11_b_IPR: fmol_per_sec {init: 2.23094, pub: out};
var kappa_R12_b_IPR: fmol_per_sec {init: 152486, pub: out};
var kappa_R1_c_IPR: fmol_per_sec {init: 0.0125066, pub: out};
var kappa_R2_c_IPR: fmol_per_sec {init: 16.6342, pub: out};
var kappa_R3_c_IPR: fmol_per_sec {init: 0.291916, pub: out};
var kappa_R4_c_IPR: fmol_per_sec {init: 1.25066, pub: out};
var kappa_R5_c_IPR: fmol_per_sec {init: 5.31398e-05, pub: out};
var kappa_R6_c_IPR: fmol_per_sec {init: 0.00973583, pub: out};
var kappa_R7_c_IPR: fmol_per_sec {init: 1.18239, pub: out};
var kappa_R8_c_IPR: fmol_per_sec {init: 4829.13, pub: out};
var kappa_R9_c_IPR: fmol_per_sec {init: 0.0955807, pub: out};
var kappa_R10_c_IPR: fmol_per_sec {init: 0.00718894, pub: out};
var kappa_R11_c_IPR: fmol_per_sec {init: 2.23094, pub: out};
var kappa_R12_c_IPR: fmol_per_sec {init: 152486, pub: out};
var kappa_pCa_R1: fmol_per_sec {init: 1654.75, pub: out};
var kappa_pCa_R2: fmol_per_sec {init: 0.00021208, pub: out};
var kappa_SERCA_R1_2: fmol_per_sec {init: 2.90729e-06, pub: out};
var kappa_SERCA_R2_4: fmol_per_sec {init: 0.166593, pub: out};
var kappa_SERCA_R2_2a: fmol_per_sec {init: 40567.4, pub: out};
var kappa_SERCA_R4_5: fmol_per_sec {init: 0.166593, pub: out};
var kappa_SERCA_R5_6: fmol_per_sec {init: 0.41262, pub: out};
var kappa_SERCA_R6_8: fmol_per_sec {init: 3.47304, pub: out};
var kappa_SERCA_R8_9: fmol_per_sec {init: 12551.8, pub: out};
var kappa_SERCA_R9_10: fmol_per_sec {init: 12551.8, pub: out};
var kappa_SERCA_R10_1: fmol_per_sec {init: 6.07543e-05, pub: out};
var K_Ca_SR: per_fmol {init: 391.219, pub: out};
var K_Ca_i: per_fmol {init: 20.1363, pub: out};
var K_IP3: per_fmol {init: 70087.3, pub: out};
var K_Sa_000_IPR: per_fmol {init: 30.734, pub: out};
var K_Sa_001_IPR: per_fmol {init: 2.31077, pub: out};
var K_Sa_010_IPR: per_fmol {init: 1.31674, pub: out};
var K_Sa_011_IPR: per_fmol {init: 0.0990011, pub: out};
var K_Sa_100_IPR: per_fmol {init: 7233.34, pub: out};
var K_Sa_101_IPR: per_fmol {init: 3948.08, pub: out};
var K_Sa_110_IPR: per_fmol {init: 0.000341014, pub: out};
var K_Sa_111_IPR: per_fmol {init: 169.148, pub: out};
var K_Sb_000_IPR: per_fmol {init: 30.734, pub: out};
var K_Sb_001_IPR: per_fmol {init: 2.31077, pub: out};
var K_Sb_010_IPR: per_fmol {init: 1.31674, pub: out};
var K_Sb_011_IPR: per_fmol {init: 0.0990011, pub: out};
var K_Sb_100_IPR: per_fmol {init: 7233.34, pub: out};
var K_Sb_101_IPR: per_fmol {init: 3948.08, pub: out};
var K_Sb_110_IPR: per_fmol {init: 0.000341014, pub: out};
var K_Sb_111_IPR: per_fmol {init: 169.148, pub: out};
var K_Sc_000_IPR: per_fmol {init: 30.734, pub: out};
var K_Sc_001_IPR: per_fmol {init: 2.31077, pub: out};
var K_Sc_010_IPR: per_fmol {init: 1.31674, pub: out};
var K_Sc_011_IPR: per_fmol {init: 0.0990011, pub: out};
var K_Sc_100_IPR: per_fmol {init: 7233.34, pub: out};
var K_Sc_101_IPR: per_fmol {init: 3948.08, pub: out};
var K_Sc_110_IPR: per_fmol {init: 0.000341014, pub: out};
var K_Sc_111_IPR: per_fmol {init: 169.148, pub: out};
var K_Ca_o: per_fmol {init: 5.02046e-05, pub: out};
var K_pCa: per_fmol {init: 0.0449472, pub: out};
var K_pCa_Ca: per_fmol {init: 0.0116935, pub: out};
var K_P1_SERCA: per_fmol {init: 10.9783, pub: out};
var K_P2_SERCA: per_fmol {init: 0.170651, pub: out};
var K_P2a_SERCA: per_fmol {init: 0.157113, pub: out};
var K_P4_SERCA: per_fmol {init: 38259, pub: out};
var K_P5_SERCA: per_fmol {init: 0.00151486, pub: out};
var K_P6_SERCA: per_fmol {init: 9208.24, pub: out};
var K_P8_SERCA: per_fmol {init: 0.00677857, pub: out};
var K_P9_SERCA: per_fmol {init: 0.50779, pub: out};
var K_P10_SERCA: per_fmol {init: 0.0830361, pub: out};
var K_H: per_fmol {init: 3108.77, pub: out};
var K_MgATP: per_fmol {init: 7.72815, pub: out};
var K_MgADP: per_fmol {init: 1.68438e-07, pub: out};
var K_Pi: per_fmol {init: 0.00212752, pub: out};
var K_PLBp: per_fmol {init: 78.0122, 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_IP3: fmol {init: 1e-888, pub: out};
var q_Sa_000_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_001_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_010_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_011_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_100_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_101_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_110_IPR: fmol {init: 1e-888, pub: out};
var q_Sa_111_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_000_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_001_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_010_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_011_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_100_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_101_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_110_IPR: fmol {init: 1e-888, pub: out};
var q_Sb_111_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_000_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_001_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_010_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_011_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_100_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_101_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_110_IPR: fmol {init: 1e-888, pub: out};
var q_Sc_111_IPR: fmol {init: 1e-888, pub: out};
var q_Ca_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_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};
var q_PLBp: fmol {init: 1e-888, pub: out};

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


// IPR imports
var v_R_IPR_main: fmol_per_sec {pub: in};
var v_R1_a_IPR: fmol_per_sec {pub: in};
var v_R2_a_IPR: fmol_per_sec {pub: in};
var v_R3_a_IPR: fmol_per_sec {pub: in};
var v_R4_a_IPR: fmol_per_sec {pub: in};
var v_R5_a_IPR: fmol_per_sec {pub: in};
var v_R6_a_IPR: fmol_per_sec {pub: in};
var v_R7_a_IPR: fmol_per_sec {pub: in};
var v_R8_a_IPR: fmol_per_sec {pub: in};
var v_R9_a_IPR: fmol_per_sec {pub: in};
var v_R10_a_IPR: fmol_per_sec {pub: in};
var v_R11_a_IPR: fmol_per_sec {pub: in};
var v_R12_a_IPR: fmol_per_sec {pub: in};
var v_R1_b_IPR: fmol_per_sec {pub: in};
var v_R2_b_IPR: fmol_per_sec {pub: in};
var v_R3_b_IPR: fmol_per_sec {pub: in};
var v_R4_b_IPR: fmol_per_sec {pub: in};
var v_R5_b_IPR: fmol_per_sec {pub: in};
var v_R6_b_IPR: fmol_per_sec {pub: in};
var v_R7_b_IPR: fmol_per_sec {pub: in};
var v_R8_b_IPR: fmol_per_sec {pub: in};
var v_R9_b_IPR: fmol_per_sec {pub: in};
var v_R10_b_IPR: fmol_per_sec {pub: in};
var v_R11_b_IPR: fmol_per_sec {pub: in};
var v_R12_b_IPR: fmol_per_sec {pub: in};
var v_R1_c_IPR: fmol_per_sec {pub: in};
var v_R2_c_IPR: fmol_per_sec {pub: in};
var v_R3_c_IPR: fmol_per_sec {pub: in};
var v_R4_c_IPR: fmol_per_sec {pub: in};
var v_R5_c_IPR: fmol_per_sec {pub: in};
var v_R6_c_IPR: fmol_per_sec {pub: in};
var v_R7_c_IPR: fmol_per_sec {pub: in};
var v_R8_c_IPR: fmol_per_sec {pub: in};
var v_R9_c_IPR: fmol_per_sec {pub: in};
var v_R10_c_IPR: fmol_per_sec {pub: in};
var v_R11_c_IPR: fmol_per_sec {pub: in};
var v_R12_c_IPR: 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};


// 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_mIPR  + v_mSERCA ;
ode(q_Ca_i, time) = + v_mCa_leak  + v_mIPR  + v_mpCa  + v_mSERCA ;
ode(q_IP3, time) = + v_mIPR ;
ode(q_Sa_000_IPR, time) = + v_mIPR ;
ode(q_Sa_001_IPR, time) = + v_mIPR ;
ode(q_Sa_010_IPR, time) = + v_mIPR ;
ode(q_Sa_011_IPR, time) = + v_mIPR ;
ode(q_Sa_100_IPR, time) = + v_mIPR ;
ode(q_Sa_101_IPR, time) = + v_mIPR ;
ode(q_Sa_110_IPR, time) = + v_mIPR ;
ode(q_Sa_111_IPR, time) = + v_mIPR ;
ode(q_Sb_000_IPR, time) = + v_mIPR ;
ode(q_Sb_001_IPR, time) = + v_mIPR ;
ode(q_Sb_010_IPR, time) = + v_mIPR ;
ode(q_Sb_011_IPR, time) = + v_mIPR ;
ode(q_Sb_100_IPR, time) = + v_mIPR ;
ode(q_Sb_101_IPR, time) = + v_mIPR ;
ode(q_Sb_110_IPR, time) = + v_mIPR ;
ode(q_Sb_111_IPR, time) = + v_mIPR ;
ode(q_Sc_000_IPR, time) = + v_mIPR ;
ode(q_Sc_001_IPR, time) = + v_mIPR ;
ode(q_Sc_010_IPR, time) = + v_mIPR ;
ode(q_Sc_011_IPR, time) = + v_mIPR ;
ode(q_Sc_100_IPR, time) = + v_mIPR ;
ode(q_Sc_101_IPR, time) = + v_mIPR ;
ode(q_Sc_110_IPR, time) = + v_mIPR ;
ode(q_Sc_111_IPR, time) = + v_mIPR ;
ode(q_Ca_o, time) = + v_mpCa ;
ode(q_pCa, time) = + v_mpCa ;
ode(q_pCa_Ca, time) = + v_mpCa ;
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 ;
ode(q_PLBp, 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 IPR for
vars time and time;
vars q_Ca_i and q_Ca_i;
vars q_Ca_SR and q_Ca_SR;
vars q_IP3 and q_IP3;
vars q_Sa_000_IPR and q_Sa_000_IPR;
vars q_Sa_001_IPR and q_Sa_001_IPR;
vars q_Sa_010_IPR and q_Sa_010_IPR;
vars q_Sa_011_IPR and q_Sa_011_IPR;
vars q_Sa_100_IPR and q_Sa_100_IPR;
vars q_Sa_101_IPR and q_Sa_101_IPR;
vars q_Sa_110_IPR and q_Sa_110_IPR;
vars q_Sa_111_IPR and q_Sa_111_IPR;
vars q_Sb_000_IPR and q_Sb_000_IPR;
vars q_Sb_001_IPR and q_Sb_001_IPR;
vars q_Sb_010_IPR and q_Sb_010_IPR;
vars q_Sb_011_IPR and q_Sb_011_IPR;
vars q_Sb_100_IPR and q_Sb_100_IPR;
vars q_Sb_101_IPR and q_Sb_101_IPR;
vars q_Sb_110_IPR and q_Sb_110_IPR;
vars q_Sb_111_IPR and q_Sb_111_IPR;
vars q_Sc_000_IPR and q_Sc_000_IPR;
vars q_Sc_001_IPR and q_Sc_001_IPR;
vars q_Sc_010_IPR and q_Sc_010_IPR;
vars q_Sc_011_IPR and q_Sc_011_IPR;
vars q_Sc_100_IPR and q_Sc_100_IPR;
vars q_Sc_101_IPR and q_Sc_101_IPR;
vars q_Sc_110_IPR and q_Sc_110_IPR;
vars q_Sc_111_IPR and q_Sc_111_IPR;
vars v_R_IPR_main and v_R_IPR_main;
vars v_R1_a_IPR and v_R1_a_IPR;
vars v_R2_a_IPR and v_R2_a_IPR;
vars v_R3_a_IPR and v_R3_a_IPR;
vars v_R4_a_IPR and v_R4_a_IPR;
vars v_R5_a_IPR and v_R5_a_IPR;
vars v_R6_a_IPR and v_R6_a_IPR;
vars v_R7_a_IPR and v_R7_a_IPR;
vars v_R8_a_IPR and v_R8_a_IPR;
vars v_R9_a_IPR and v_R9_a_IPR;
vars v_R10_a_IPR and v_R10_a_IPR;
vars v_R11_a_IPR and v_R11_a_IPR;
vars v_R12_a_IPR and v_R12_a_IPR;
vars v_R1_b_IPR and v_R1_b_IPR;
vars v_R2_b_IPR and v_R2_b_IPR;
vars v_R3_b_IPR and v_R3_b_IPR;
vars v_R4_b_IPR and v_R4_b_IPR;
vars v_R5_b_IPR and v_R5_b_IPR;
vars v_R6_b_IPR and v_R6_b_IPR;
vars v_R7_b_IPR and v_R7_b_IPR;
vars v_R8_b_IPR and v_R8_b_IPR;
vars v_R9_b_IPR and v_R9_b_IPR;
vars v_R10_b_IPR and v_R10_b_IPR;
vars v_R11_b_IPR and v_R11_b_IPR;
vars v_R12_b_IPR and v_R12_b_IPR;
vars v_R1_c_IPR and v_R1_c_IPR;
vars v_R2_c_IPR and v_R2_c_IPR;
vars v_R3_c_IPR and v_R3_c_IPR;
vars v_R4_c_IPR and v_R4_c_IPR;
vars v_R5_c_IPR and v_R5_c_IPR;
vars v_R6_c_IPR and v_R6_c_IPR;
vars v_R7_c_IPR and v_R7_c_IPR;
vars v_R8_c_IPR and v_R8_c_IPR;
vars v_R9_c_IPR and v_R9_c_IPR;
vars v_R10_c_IPR and v_R10_c_IPR;
vars v_R11_c_IPR and v_R11_c_IPR;
vars v_R12_c_IPR and v_R12_c_IPR;
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 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 q_PLBp and q_PLBp;
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 IPR for
vars kappa_R_IPR_main and kappa_R_IPR_main;
vars kappa_R1_a_IPR and kappa_R1_a_IPR;
vars kappa_R2_a_IPR and kappa_R2_a_IPR;
vars kappa_R3_a_IPR and kappa_R3_a_IPR;
vars kappa_R4_a_IPR and kappa_R4_a_IPR;
vars kappa_R5_a_IPR and kappa_R5_a_IPR;
vars kappa_R6_a_IPR and kappa_R6_a_IPR;
vars kappa_R7_a_IPR and kappa_R7_a_IPR;
vars kappa_R8_a_IPR and kappa_R8_a_IPR;
vars kappa_R9_a_IPR and kappa_R9_a_IPR;
vars kappa_R10_a_IPR and kappa_R10_a_IPR;
vars kappa_R11_a_IPR and kappa_R11_a_IPR;
vars kappa_R12_a_IPR and kappa_R12_a_IPR;
vars kappa_R1_b_IPR and kappa_R1_b_IPR;
vars kappa_R2_b_IPR and kappa_R2_b_IPR;
vars kappa_R3_b_IPR and kappa_R3_b_IPR;
vars kappa_R4_b_IPR and kappa_R4_b_IPR;
vars kappa_R5_b_IPR and kappa_R5_b_IPR;
vars kappa_R6_b_IPR and kappa_R6_b_IPR;
vars kappa_R7_b_IPR and kappa_R7_b_IPR;
vars kappa_R8_b_IPR and kappa_R8_b_IPR;
vars kappa_R9_b_IPR and kappa_R9_b_IPR;
vars kappa_R10_b_IPR and kappa_R10_b_IPR;
vars kappa_R11_b_IPR and kappa_R11_b_IPR;
vars kappa_R12_b_IPR and kappa_R12_b_IPR;
vars kappa_R1_c_IPR and kappa_R1_c_IPR;
vars kappa_R2_c_IPR and kappa_R2_c_IPR;
vars kappa_R3_c_IPR and kappa_R3_c_IPR;
vars kappa_R4_c_IPR and kappa_R4_c_IPR;
vars kappa_R5_c_IPR and kappa_R5_c_IPR;
vars kappa_R6_c_IPR and kappa_R6_c_IPR;
vars kappa_R7_c_IPR and kappa_R7_c_IPR;
vars kappa_R8_c_IPR and kappa_R8_c_IPR;
vars kappa_R9_c_IPR and kappa_R9_c_IPR;
vars kappa_R10_c_IPR and kappa_R10_c_IPR;
vars kappa_R11_c_IPR and kappa_R11_c_IPR;
vars kappa_R12_c_IPR and kappa_R12_c_IPR;
vars K_Ca_i and K_Ca_i;
vars K_Ca_SR and K_Ca_SR;
vars K_IP3 and K_IP3;
vars K_Sa_000_IPR and K_Sa_000_IPR;
vars K_Sa_001_IPR and K_Sa_001_IPR;
vars K_Sa_010_IPR and K_Sa_010_IPR;
vars K_Sa_011_IPR and K_Sa_011_IPR;
vars K_Sa_100_IPR and K_Sa_100_IPR;
vars K_Sa_101_IPR and K_Sa_101_IPR;
vars K_Sa_110_IPR and K_Sa_110_IPR;
vars K_Sa_111_IPR and K_Sa_111_IPR;
vars K_Sb_000_IPR and K_Sb_000_IPR;
vars K_Sb_001_IPR and K_Sb_001_IPR;
vars K_Sb_010_IPR and K_Sb_010_IPR;
vars K_Sb_011_IPR and K_Sb_011_IPR;
vars K_Sb_100_IPR and K_Sb_100_IPR;
vars K_Sb_101_IPR and K_Sb_101_IPR;
vars K_Sb_110_IPR and K_Sb_110_IPR;
vars K_Sb_111_IPR and K_Sb_111_IPR;
vars K_Sc_000_IPR and K_Sc_000_IPR;
vars K_Sc_001_IPR and K_Sc_001_IPR;
vars K_Sc_010_IPR and K_Sc_010_IPR;
vars K_Sc_011_IPR and K_Sc_011_IPR;
vars K_Sc_100_IPR and K_Sc_100_IPR;
vars K_Sc_101_IPR and K_Sc_101_IPR;
vars K_Sc_110_IPR and K_Sc_110_IPR;
vars K_Sc_111_IPR and K_Sc_111_IPR;
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 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;
vars K_PLBp and K_PLBp;
enddef;

def map between constants and Ca_leak for
	vars R and R;
	vars T and T;
enddef;
def map between constants and IPR 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 SERCA for
	vars R and R;
	vars T and T;
enddef;

enddef;