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 "PLB/BG_PLB.cellml" for
comp PLB using comp PLB;
enddef;
def import using "SERCA/BG_SERCA.cellml" for
comp SERCA using comp SERCA;
enddef;
def import using "RyR/BG_RyR.cellml" for
comp RyR using comp RyR;
enddef;
def import using "LCC/BG_LCC.cellml" for
comp LCC using comp LCC;
enddef;
def import using "TRPN/BG_TRPN.cellml" for
comp TRPN using comp TRPN;
enddef;

def comp BG_parameters as
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_SERCA_R1_2: fmol_per_sec {init: 2.24963e-05, pub: out};
var kappa_SERCA_R2_4: fmol_per_sec {init: 108329, pub: out};
var kappa_SERCA_R2_2a: fmol_per_sec {init: 2.77146e+06, pub: out};
var kappa_SERCA_R4_5: fmol_per_sec {init: 108329, pub: out};
var kappa_SERCA_R5_6: fmol_per_sec {init: 3442.07, pub: out};
var kappa_SERCA_R6_8: fmol_per_sec {init: 3.07427e+06, pub: out};
var kappa_SERCA_R8_9: fmol_per_sec {init: 1.0306e+07, pub: out};
var kappa_SERCA_R9_10: fmol_per_sec {init: 1.0306e+07, pub: out};
var kappa_SERCA_R10_1: fmol_per_sec {init: 0.00565005, pub: out};
var kappa_RyR: fmol_per_sec {init: 3.58172e+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_Ca1: fmol_per_sec {init: 8.37314, pub: out};
var kappa_Ca2: fmol_per_sec {init: 15.6575, pub: out};
var kappa_K1: fmol_per_sec {init: 0.0929479, pub: out};
var kappa_K2: fmol_per_sec {init: 0.17381, pub: out};
var kappa_d000: fmol_per_sec {init: 942.411, pub: out};
var kappa_d010: fmol_per_sec {init: 23.2508, pub: out};
var kappa_d020: fmol_per_sec {init: 43.4783, pub: out};
var kappa_d001: fmol_per_sec {init: 0.175516, pub: out};
var kappa_d011: fmol_per_sec {init: 0.00433026, pub: out};
var kappa_d021: fmol_per_sec {init: 0.00809745, pub: out};
var kappa_f1_000: fmol_per_sec {init: 1.73616, pub: out};
var kappa_f1_100: fmol_per_sec {init: 8.62132, pub: out};
var kappa_f1_001: fmol_per_sec {init: 0.000323344, pub: out};
var kappa_f1_101: fmol_per_sec {init: 0.00160565, pub: out};
var kappa_f2_000: fmol_per_sec {init: 3.17968, pub: out};
var kappa_f2_100: fmol_per_sec {init: 15.7895, pub: out};
var kappa_f2_001: fmol_per_sec {init: 0.000592187, pub: out};
var kappa_f2_101: fmol_per_sec {init: 0.00294065, pub: out};
var kappa_f3_010: fmol_per_sec {init: 8932.14, pub: out};
var kappa_f3_110: fmol_per_sec {init: 44354.8, pub: out};
var kappa_f3_011: fmol_per_sec {init: 1.66353, pub: out};
var kappa_f3_111: fmol_per_sec {init: 8.2607, pub: out};
var kappa_fCa000: fmol_per_sec {init: 60096.5, pub: out};
var kappa_fCa100: fmol_per_sec {init: 298424, pub: out};
var kappa_fCa010: fmol_per_sec {init: 1482.67, pub: out};
var kappa_fCa110: fmol_per_sec {init: 7362.59, pub: out};
var kappa_fCa020: fmol_per_sec {init: 2772.56, pub: out};
var kappa_fCa120: fmol_per_sec {init: 13767.8, pub: out};
var kappa_R_TRPNCa: fmol_per_sec {init: 223.356, pub: out};
var K_PLB: per_fmol {init: 0.00235741, pub: out};
var K_PKACI: per_fmol {init: 0.395191, pub: out};
var K_PLB_PKACI: per_fmol {init: 0.638527, pub: out};
var K_PP1: per_fmol {init: 0.361989, pub: out};
var K_PLBp_PP1: per_fmol {init: 0.203757, pub: out};
var K_PLBp: per_fmol {init: 0.326014, pub: out};
var K_Ip: per_fmol {init: 5.41093, pub: out};
var K_Ip_PP1: per_fmol {init: 0.0673793, pub: out};
var K_P1_SERCA: per_fmol {init: 906.879, pub: out};
var K_P2_SERCA: per_fmol {init: 2605.18, pub: out};
var K_P2a_SERCA: per_fmol {init: 271.664, pub: out};
var K_P4_SERCA: per_fmol {init: 6950.17, pub: out};
var K_P5_SERCA: per_fmol {init: 0.0214514, pub: out};
var K_P6_SERCA: per_fmol {init: 1228.84, pub: out};
var K_P8_SERCA: per_fmol {init: 76.0193, pub: out};
var K_P9_SERCA: per_fmol {init: 73.055, pub: out};
var K_P10_SERCA: per_fmol {init: 105.473, pub: out};
var K_H: per_fmol {init: 278.106, pub: out};
var K_Cai: per_fmol {init: 0.052177, pub: out};
var K_Ca_SR: per_fmol {init: 0.052177, pub: out};
var K_MgATP: per_fmol {init: 1072.81, pub: out};
var K_MgADP: per_fmol {init: 1.34257e-05, pub: out};
var K_P: per_fmol {init: 0.0245735, pub: out};
var K_C_RyR: per_fmol {init: 115.191, pub: out};
var K_CI_RyR: per_fmol {init: 1.15191, pub: out};
var K_I_RyR: per_fmol {init: 0.00197471, pub: out};
var K_O_RyR: per_fmol {init: 0.197471, pub: out};
var K_Cao: per_fmol {init: 0.052177, pub: out};
var K_Ki: per_fmol {init: 0.00369486, pub: out};
var K_Ko: per_fmol {init: 0.00369486, pub: out};
var K_000_LCC: per_fmol {init: 0.0150147, pub: out};
var K_010_LCC: per_fmol {init: 0.608585, pub: out};
var K_020_LCC: per_fmol {init: 0.325451, pub: out};
var K_100_LCC: per_fmol {init: 0.00302366, pub: out};
var K_110_LCC: per_fmol {init: 0.122556, pub: out};
var K_120_LCC: per_fmol {init: 0.0655392, pub: out};
var K_001_LCC: per_fmol {init: 80.6198, pub: out};
var K_011_LCC: per_fmol {init: 3267.72, pub: out};
var K_021_LCC: per_fmol {init: 1747.47, pub: out};
var K_101_LCC: per_fmol {init: 16.2352, pub: out};
var K_111_LCC: per_fmol {init: 658.051, pub: out};
var K_121_LCC: per_fmol {init: 351.905, pub: out};
var K_TRPN: per_fmol {init: 7.25115, pub: out};
var K_Ca_TRPN: per_fmol {init: 0.02603, 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};
                // stimulus
                // ramp UP and ramp DOWN
                var stimSt: second {init: 3.5e-4};
                var stimDur: second {init: 0.25e-4};
                var tRamp: second {init: 1.8e-4};
                var stimMag: fmol {init: 1e1};
                var stimHolding: fmol {init: 1e-5};  
                var m: fmol_per_sec;  
                m = stimMag/tRamp;   
                q_L_B1_init = sel          
                    case (time < stimSt) and (time > stimSt-tRamp):   
                        stimHolding+m*(time-stimSt+tRamp);
                    case (time >= stimSt) and (time < stimSt+stimDur): 
                        stimMag+stimHolding;   
                    case (time < stimSt+tRamp+stimDur) and (time >= stimSt+stimDur):  
                        stimHolding+-m*(time-stimSt-tRamp-stimDur);   
                    otherwise:             
                        stimHolding;       
                endsel;
var q_PLB_init: fmol {init: 1e-888};
var q_PKACI_init: fmol {init: 1e-888};
var q_PLB_PKACI_init: fmol {init: 1e-888};
var q_PP1_init: fmol {init: 1e-888};
var q_PLBp_PP1_init: fmol {init: 1e-888};
var q_PLBp_init: fmol {init: 1e-888};
var q_Ip_init: fmol {init: 1e-888};
var q_Ip_PP1_init: fmol {init: 1e-888};
var q_P1_SERCA_init: fmol {init: 1e-888};
var q_P2_SERCA_init: fmol {init: 1e-888};
var q_P2a_SERCA_init: fmol {init: 1e-888};
var q_P4_SERCA_init: fmol {init: 1e-888};
var q_P5_SERCA_init: fmol {init: 1e-888};
var q_P6_SERCA_init: fmol {init: 1e-888};
var q_P8_SERCA_init: fmol {init: 1e-888};
var q_P9_SERCA_init: fmol {init: 1e-888};
var q_P10_SERCA_init: fmol {init: 1e-888};
var q_H_init: fmol {init: 1e-888};
var q_Cai_init: fmol {init: 1e-888};
var q_Ca_SR_init: fmol {init: 1e-888};
var q_MgATP_init: fmol {init: 1e-888};
var q_MgADP_init: fmol {init: 1e-888};
var q_P_init: fmol {init: 1e-888};
var q_C_RyR_init: fmol {init: 1e-888};
var q_CI_RyR_init: fmol {init: 1e-888};
var q_I_RyR_init: fmol {init: 1e-888};
var q_O_RyR_init: fmol {init: 1e-888};
var q_Cao_init: fmol {init: 1e-888};
var q_Ki_init: fmol {init: 1e-888};
var q_Ko_init: fmol {init: 1e-888};
var q_000_LCC_init: fmol {init: 1e-888};
var q_010_LCC_init: fmol {init: 1e-888};
var q_020_LCC_init: fmol {init: 1e-888};
var q_100_LCC_init: fmol {init: 1e-888};
var q_110_LCC_init: fmol {init: 1e-888};
var q_120_LCC_init: fmol {init: 1e-888};
var q_001_LCC_init: fmol {init: 1e-888};
var q_011_LCC_init: fmol {init: 1e-888};
var q_021_LCC_init: fmol {init: 1e-888};
var q_101_LCC_init: fmol {init: 1e-888};
var q_111_LCC_init: fmol {init: 1e-888};
var q_121_LCC_init: fmol {init: 1e-888};
var q_TRPN_init: fmol {init: 1e-888};
var q_Ca_TRPN_init: fmol {init: 1e-888};

// Global value
var q_PLB: fmol {pub: out};
var q_PKACI: fmol {pub: out};
var q_PLB_PKACI: fmol {pub: out};
var q_PP1: fmol {pub: out};
var q_PLBp_PP1: fmol {pub: out};
var q_PLBp: fmol {pub: out};
var q_Ip: fmol {pub: out};
var q_Ip_PP1: fmol {pub: out};
var q_P1_SERCA: fmol {pub: out};
var q_P2_SERCA: fmol {pub: out};
var q_P2a_SERCA: fmol {pub: out};
var q_P4_SERCA: fmol {pub: out};
var q_P5_SERCA: fmol {pub: out};
var q_P6_SERCA: fmol {pub: out};
var q_P8_SERCA: fmol {pub: out};
var q_P9_SERCA: fmol {pub: out};
var q_P10_SERCA: fmol {pub: out};
var q_H: fmol {pub: out};
var q_Cai: fmol {pub: out};
var q_Ca_SR: fmol {pub: out};
var q_MgATP: fmol {pub: out};
var q_MgADP: fmol {pub: out};
var q_P: fmol {pub: out};
var q_C_RyR: fmol {pub: out};
var q_CI_RyR: fmol {pub: out};
var q_I_RyR: fmol {pub: out};
var q_O_RyR: fmol {pub: out};
var q_Cao: fmol {pub: out};
var q_Ki: fmol {pub: out};
var q_Ko: fmol {pub: out};
var q_000_LCC: fmol {pub: out};
var q_010_LCC: fmol {pub: out};
var q_020_LCC: fmol {pub: out};
var q_100_LCC: fmol {pub: out};
var q_110_LCC: fmol {pub: out};
var q_120_LCC: fmol {pub: out};
var q_001_LCC: fmol {pub: out};
var q_011_LCC: fmol {pub: out};
var q_021_LCC: fmol {pub: out};
var q_101_LCC: fmol {pub: out};
var q_111_LCC: fmol {pub: out};
var q_121_LCC: fmol {pub: out};
var q_TRPN: fmol {pub: out};
var q_Ca_TRPN: fmol {pub: out};

// PLB imports
var q_PLB_mPLB: fmol {pub: in};
var q_PKACI_mPLB: fmol {pub: in};
var q_PLB_PKACI_mPLB: fmol {pub: in};
var q_PP1_mPLB: fmol {pub: in};
var q_PLBp_PP1_mPLB: fmol {pub: in};
var q_PLBp_mPLB: fmol {pub: in};
var q_Ip_mPLB: fmol {pub: in};
var q_Ip_PP1_mPLB: fmol {pub: in};


// SERCA imports
var q_P1_SERCA_mSERCA: fmol {pub: in};
var q_P2_SERCA_mSERCA: fmol {pub: in};
var q_P2a_SERCA_mSERCA: fmol {pub: in};
var q_P4_SERCA_mSERCA: fmol {pub: in};
var q_P5_SERCA_mSERCA: fmol {pub: in};
var q_P6_SERCA_mSERCA: fmol {pub: in};
var q_P8_SERCA_mSERCA: fmol {pub: in};
var q_P9_SERCA_mSERCA: fmol {pub: in};
var q_P10_SERCA_mSERCA: fmol {pub: in};
var q_H_mSERCA: fmol {pub: in};
var q_Cai_mSERCA: fmol {pub: in};
var q_Ca_SR_mSERCA: fmol {pub: in};
var q_MgATP_mSERCA: fmol {pub: in};
var q_MgADP_mSERCA: fmol {pub: in};
var q_P_mSERCA: fmol {pub: in};


// RyR imports
var q_Ca_SR_mRyR: fmol {pub: in};
var q_Cai_mRyR: fmol {pub: in};
var q_C_RyR_mRyR: fmol {pub: in};
var q_CI_RyR_mRyR: fmol {pub: in};
var q_I_RyR_mRyR: fmol {pub: in};
var q_O_RyR_mRyR: fmol {pub: in};


// LCC imports
var q_Cai_mLCC: fmol {pub: in};
var q_Cao_mLCC: fmol {pub: in};
var q_Ki_mLCC: fmol {pub: in};
var q_Ko_mLCC: fmol {pub: in};
var q_000_LCC_mLCC: fmol {pub: in};
var q_010_LCC_mLCC: fmol {pub: in};
var q_020_LCC_mLCC: fmol {pub: in};
var q_100_LCC_mLCC: fmol {pub: in};
var q_110_LCC_mLCC: fmol {pub: in};
var q_120_LCC_mLCC: fmol {pub: in};
var q_001_LCC_mLCC: fmol {pub: in};
var q_011_LCC_mLCC: fmol {pub: in};
var q_021_LCC_mLCC: fmol {pub: in};
var q_101_LCC_mLCC: fmol {pub: in};
var q_111_LCC_mLCC: fmol {pub: in};
var q_121_LCC_mLCC: fmol {pub: in};


// TRPN imports
var q_Cai_mTRPN: fmol {pub: in};
var q_TRPN_mTRPN: fmol {pub: in};
var q_Ca_TRPN_mTRPN: fmol {pub: in};


q_PLB = q_PLB_init + q_PLB_mPLB ;
q_PKACI = q_PKACI_init + q_PKACI_mPLB ;
q_PLB_PKACI = q_PLB_PKACI_init + q_PLB_PKACI_mPLB ;
q_PP1 = q_PP1_init + q_PP1_mPLB ;
q_PLBp_PP1 = q_PLBp_PP1_init + q_PLBp_PP1_mPLB ;
q_PLBp = q_PLBp_init + q_PLBp_mPLB ;
q_Ip = q_Ip_init + q_Ip_mPLB ;
q_Ip_PP1 = q_Ip_PP1_init + q_Ip_PP1_mPLB ;
q_P1_SERCA = q_P1_SERCA_init + q_P1_SERCA_mSERCA ;
q_P2_SERCA = q_P2_SERCA_init + q_P2_SERCA_mSERCA ;
q_P2a_SERCA = q_P2a_SERCA_init + q_P2a_SERCA_mSERCA ;
q_P4_SERCA = q_P4_SERCA_init + q_P4_SERCA_mSERCA ;
q_P5_SERCA = q_P5_SERCA_init + q_P5_SERCA_mSERCA ;
q_P6_SERCA = q_P6_SERCA_init + q_P6_SERCA_mSERCA ;
q_P8_SERCA = q_P8_SERCA_init + q_P8_SERCA_mSERCA ;
q_P9_SERCA = q_P9_SERCA_init + q_P9_SERCA_mSERCA ;
q_P10_SERCA = q_P10_SERCA_init + q_P10_SERCA_mSERCA ;
q_H = q_H_init + q_H_mSERCA ;
q_Cai = q_Cai_init + q_Cai_mSERCA  + q_Cai_mRyR  + q_Cai_mLCC  + q_Cai_mTRPN ;
q_Ca_SR = q_Ca_SR_init + q_Ca_SR_mSERCA  + q_Ca_SR_mRyR ;
q_MgATP = q_MgATP_init + q_MgATP_mSERCA ;
q_MgADP = q_MgADP_init + q_MgADP_mSERCA ;
q_P = q_P_init + q_P_mSERCA ;
q_C_RyR = q_C_RyR_init + q_C_RyR_mRyR ;
q_CI_RyR = q_CI_RyR_init + q_CI_RyR_mRyR ;
q_I_RyR = q_I_RyR_init + q_I_RyR_mRyR ;
q_O_RyR = q_O_RyR_init + q_O_RyR_mRyR ;
q_Cao = q_Cao_init + q_Cao_mLCC ;
q_Ki = q_Ki_init + q_Ki_mLCC ;
q_Ko = q_Ko_init + q_Ko_mLCC ;
q_000_LCC = q_000_LCC_init + q_000_LCC_mLCC ;
q_010_LCC = q_010_LCC_init + q_010_LCC_mLCC ;
q_020_LCC = q_020_LCC_init + q_020_LCC_mLCC ;
q_100_LCC = q_100_LCC_init + q_100_LCC_mLCC ;
q_110_LCC = q_110_LCC_init + q_110_LCC_mLCC ;
q_120_LCC = q_120_LCC_init + q_120_LCC_mLCC ;
q_001_LCC = q_001_LCC_init + q_001_LCC_mLCC ;
q_011_LCC = q_011_LCC_init + q_011_LCC_mLCC ;
q_021_LCC = q_021_LCC_init + q_021_LCC_mLCC ;
q_101_LCC = q_101_LCC_init + q_101_LCC_mLCC ;
q_111_LCC = q_111_LCC_init + q_111_LCC_mLCC ;
q_121_LCC = q_121_LCC_init + q_121_LCC_mLCC ;
q_TRPN = q_TRPN_init + q_TRPN_mTRPN ;
q_Ca_TRPN = q_Ca_TRPN_init + q_Ca_TRPN_mTRPN ;
enddef;

def map between environment and PLB for
vars time and time;
vars q_PLB_mPLB and q_PLB;
vars q_PLB and q_PLB_global;
vars q_PKACI_mPLB and q_PKACI;
vars q_PKACI and q_PKACI_global;
vars q_PLB_PKACI_mPLB and q_PLB_PKACI;
vars q_PLB_PKACI and q_PLB_PKACI_global;
vars q_PP1_mPLB and q_PP1;
vars q_PP1 and q_PP1_global;
vars q_PLBp_PP1_mPLB and q_PLBp_PP1;
vars q_PLBp_PP1 and q_PLBp_PP1_global;
vars q_PLBp_mPLB and q_PLBp;
vars q_PLBp and q_PLBp_global;
vars q_Ip_mPLB and q_Ip;
vars q_Ip and q_Ip_global;
vars q_Ip_PP1_mPLB and q_Ip_PP1;
vars q_Ip_PP1 and q_Ip_PP1_global;
enddef;

def map between environment and SERCA for
vars time and time;
vars q_P1_SERCA_mSERCA and q_P1_SERCA;
vars q_P1_SERCA and q_P1_SERCA_global;
vars q_P2_SERCA_mSERCA and q_P2_SERCA;
vars q_P2_SERCA and q_P2_SERCA_global;
vars q_P2a_SERCA_mSERCA and q_P2a_SERCA;
vars q_P2a_SERCA and q_P2a_SERCA_global;
vars q_P4_SERCA_mSERCA and q_P4_SERCA;
vars q_P4_SERCA and q_P4_SERCA_global;
vars q_P5_SERCA_mSERCA and q_P5_SERCA;
vars q_P5_SERCA and q_P5_SERCA_global;
vars q_P6_SERCA_mSERCA and q_P6_SERCA;
vars q_P6_SERCA and q_P6_SERCA_global;
vars q_P8_SERCA_mSERCA and q_P8_SERCA;
vars q_P8_SERCA and q_P8_SERCA_global;
vars q_P9_SERCA_mSERCA and q_P9_SERCA;
vars q_P9_SERCA and q_P9_SERCA_global;
vars q_P10_SERCA_mSERCA and q_P10_SERCA;
vars q_P10_SERCA and q_P10_SERCA_global;
vars q_H_mSERCA and q_H;
vars q_H and q_H_global;
vars q_Cai_mSERCA and q_Cai;
vars q_Cai and q_Cai_global;
vars q_Ca_SR_mSERCA and q_Ca_SR;
vars q_Ca_SR and q_Ca_SR_global;
vars q_MgATP_mSERCA and q_MgATP;
vars q_MgATP and q_MgATP_global;
vars q_MgADP_mSERCA and q_MgADP;
vars q_MgADP and q_MgADP_global;
vars q_P_mSERCA and q_P;
vars q_P and q_P_global;
enddef;

def map between environment and RyR for
vars time and time;
vars q_Ca_SR_mRyR and q_Ca_SR;
vars q_Ca_SR and q_Ca_SR_global;
vars q_Cai_mRyR and q_Cai;
vars q_Cai and q_Cai_global;
vars q_C_RyR_mRyR and q_C_RyR;
vars q_C_RyR and q_C_RyR_global;
vars q_CI_RyR_mRyR and q_CI_RyR;
vars q_CI_RyR and q_CI_RyR_global;
vars q_I_RyR_mRyR and q_I_RyR;
vars q_I_RyR and q_I_RyR_global;
vars q_O_RyR_mRyR and q_O_RyR;
vars q_O_RyR and q_O_RyR_global;
enddef;

def map between environment and LCC for
vars time and time;
vars q_Cai_mLCC and q_Cai;
vars q_Cai and q_Cai_global;
vars q_Cao_mLCC and q_Cao;
vars q_Cao and q_Cao_global;
vars q_Ki_mLCC and q_Ki;
vars q_Ki and q_Ki_global;
vars q_Ko_mLCC and q_Ko;
vars q_Ko and q_Ko_global;
vars q_000_LCC_mLCC and q_000_LCC;
vars q_000_LCC and q_000_LCC_global;
vars q_010_LCC_mLCC and q_010_LCC;
vars q_010_LCC and q_010_LCC_global;
vars q_020_LCC_mLCC and q_020_LCC;
vars q_020_LCC and q_020_LCC_global;
vars q_100_LCC_mLCC and q_100_LCC;
vars q_100_LCC and q_100_LCC_global;
vars q_110_LCC_mLCC and q_110_LCC;
vars q_110_LCC and q_110_LCC_global;
vars q_120_LCC_mLCC and q_120_LCC;
vars q_120_LCC and q_120_LCC_global;
vars q_001_LCC_mLCC and q_001_LCC;
vars q_001_LCC and q_001_LCC_global;
vars q_011_LCC_mLCC and q_011_LCC;
vars q_011_LCC and q_011_LCC_global;
vars q_021_LCC_mLCC and q_021_LCC;
vars q_021_LCC and q_021_LCC_global;
vars q_101_LCC_mLCC and q_101_LCC;
vars q_101_LCC and q_101_LCC_global;
vars q_111_LCC_mLCC and q_111_LCC;
vars q_111_LCC and q_111_LCC_global;
vars q_121_LCC_mLCC and q_121_LCC;
vars q_121_LCC and q_121_LCC_global;
enddef;

def map between environment and TRPN for
vars time and time;
vars q_Cai_mTRPN and q_Cai;
vars q_Cai and q_Cai_global;
vars q_TRPN_mTRPN and q_TRPN;
vars q_TRPN and q_TRPN_global;
vars q_Ca_TRPN_mTRPN and q_Ca_TRPN;
vars q_Ca_TRPN and q_Ca_TRPN_global;
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 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_Cai and K_Cai;
vars K_Ca_SR and K_Ca_SR;
vars K_MgATP and K_MgATP;
vars K_MgADP and K_MgADP;
vars K_P and K_P;
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_Cai and K_Cai;
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 LCC for
vars kappa_Ca1 and kappa_Ca1;
vars kappa_Ca2 and kappa_Ca2;
vars kappa_K1 and kappa_K1;
vars kappa_K2 and kappa_K2;
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_Cai and K_Cai;
vars K_Cao and K_Cao;
vars K_Ki and K_Ki;
vars K_Ko and K_Ko;
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 TRPN for
vars kappa_R_TRPNCa and kappa_R_TRPNCa;
vars K_Cai and K_Cai;
vars K_TRPN and K_TRPN;
vars K_Ca_TRPN and K_Ca_TRPN;
enddef;

def map between constants and PLB 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;
def map between constants and RyR for
	vars R and R;
	vars T and T;
enddef;
def map between constants and LCC for
	vars R and R;
	vars T and T;
enddef;
def map between constants and TRPN for
	vars R and R;
	vars T and T;
enddef;

enddef;