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 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: 0.0108383, pub: out};
var kappa_SERCA_R2_4: fmol_per_sec {init: 192.314, pub: out};
var kappa_SERCA_R2_2a: fmol_per_sec {init: 1.00791e+08, pub: out};
var kappa_SERCA_R4_5: fmol_per_sec {init: 192.314, pub: out};
var kappa_SERCA_R5_6: fmol_per_sec {init: 1072.41, pub: out};
var kappa_SERCA_R6_8: fmol_per_sec {init: 8.43318e+06, pub: out};
var kappa_SERCA_R8_9: fmol_per_sec {init: 1.22619e+06, pub: out};
var kappa_SERCA_R9_10: fmol_per_sec {init: 1.22619e+06, pub: out};
var kappa_SERCA_R10_1: fmol_per_sec {init: 0.00890543, pub: out};
var kappa_RyR: fmol_per_sec {init: 34574.2, 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: 115.689, pub: out};
var kappa_Ca2: fmol_per_sec {init: 216.336, pub: out};
var kappa_K1: fmol_per_sec {init: 0.160891, pub: out};
var kappa_K2: fmol_per_sec {init: 0.300862, pub: out};
var kappa_d000: fmol_per_sec {init: 132.798, pub: out};
var kappa_d010: fmol_per_sec {init: 3.27634, pub: out};
var kappa_d020: fmol_per_sec {init: 6.12666, pub: out};
var kappa_d001: fmol_per_sec {init: 0.525901, pub: out};
var kappa_d011: fmol_per_sec {init: 0.0129748, pub: out};
var kappa_d021: fmol_per_sec {init: 0.0242625, pub: out};
var kappa_f1_000: fmol_per_sec {init: 0.244648, pub: out};
var kappa_f1_100: fmol_per_sec {init: 1.21486, pub: out};
var kappa_f1_001: fmol_per_sec {init: 0.000968841, pub: out};
var kappa_f1_101: fmol_per_sec {init: 0.00481102, pub: out};
var kappa_f2_000: fmol_per_sec {init: 0.448058, pub: out};
var kappa_f2_100: fmol_per_sec {init: 2.22495, pub: out};
var kappa_f2_001: fmol_per_sec {init: 0.00177438, pub: out};
var kappa_f2_101: fmol_per_sec {init: 0.00881111, pub: out};
var kappa_f3_010: fmol_per_sec {init: 1258.66, pub: out};
var kappa_f3_110: fmol_per_sec {init: 6250.17, pub: out};
var kappa_f3_011: fmol_per_sec {init: 4.98447, pub: out};
var kappa_f3_111: fmol_per_sec {init: 24.7516, pub: out};
var kappa_fCa000: fmol_per_sec {init: 180068, pub: out};
var kappa_fCa100: fmol_per_sec {init: 894171, pub: out};
var kappa_fCa010: fmol_per_sec {init: 4442.55, pub: out};
var kappa_fCa110: fmol_per_sec {init: 22060.6, pub: out};
var kappa_fCa020: fmol_per_sec {init: 8307.45, pub: out};
var kappa_fCa120: fmol_per_sec {init: 41252.7, 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: 99.4026, pub: out};
var K_P2_SERCA: per_fmol {init: 5.40736, pub: out};
var K_P2a_SERCA: per_fmol {init: 7.46995, pub: out};
var K_P4_SERCA: per_fmol {init: 3.91498e+06, pub: out};
var K_P5_SERCA: per_fmol {init: 0.0688518, pub: out};
var K_P6_SERCA: per_fmol {init: 447.967, pub: out};
var K_P8_SERCA: per_fmol {init: 3.64069, pub: out};
var K_P9_SERCA: per_fmol {init: 614.023, pub: out};
var K_P10_SERCA: per_fmol {init: 66.9176, pub: out};
var K_H: per_fmol {init: 3684.23, pub: out};
var K_Ca_i: per_fmol {init: 27.1814, pub: out};
var K_Ca_sr: per_fmol {init: 0.143954, pub: out};
var K_MgATP: per_fmol {init: 20.3153, pub: out};
var K_MgADP: per_fmol {init: 0.000118208, pub: out};
var K_Pi: per_fmol {init: 0.142238, pub: out};
var K_Ca_SR: per_fmol {init: 5.40528, pub: out};
var K_Ca_di: per_fmol {init: 5.40528, 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_Cai: per_fmol {init: 0.0113152, pub: out};
var K_Cao: per_fmol {init: 0.0113152, pub: out};
var K_Ki: per_fmol {init: 0.00639576, pub: out};
var K_Ko: per_fmol {init: 0.00639576, pub: out};
var K_000_LCC: per_fmol {init: 0.106553, pub: out};
var K_010_LCC: per_fmol {init: 4.31886, pub: out};
var K_020_LCC: per_fmol {init: 2.30959, pub: out};
var K_100_LCC: per_fmol {init: 0.0214576, pub: out};
var K_110_LCC: per_fmol {init: 0.86973, pub: out};
var K_120_LCC: per_fmol {init: 0.465103, pub: out};
var K_001_LCC: per_fmol {init: 26.9063, pub: out};
var K_011_LCC: per_fmol {init: 1090.58, pub: out};
var K_021_LCC: per_fmol {init: 583.207, pub: out};
var K_101_LCC: per_fmol {init: 5.41838, pub: out};
var K_111_LCC: per_fmol {init: 219.621, pub: out};
var K_121_LCC: per_fmol {init: 117.446, 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_Ca_i_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_Pi_init: fmol {init: 1e-888};
var q_Ca_SR_init: fmol {init: 1e-888};
var q_Ca_di_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_Cai_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};

// mass conservation checks
 var L_B1_T: fmol;
            var R_B1_T: fmol;
 var Gs_T: fmol;
 var adenosine_T: fmol;
        L_B1_T = q_L_B1+q_LR_B1Gs+q_LR_B1+q_LR_B1_aby+q_LR_B1_aby_T;
            R_B1_T = q_R_B1+q_R_B1Gs+q_LR_B1+q_LR_B1Gs+q_R_B1_aby+q_R_B1_aby_T+q_LR_B1_aby+q_LR_B1_aby_T;
            Gs_T = q_Gs+q_R_B1Gs+q_LR_B1Gs+q_a_Gs_GTP+q_a_Gs_GDP+q_R_B1_aby+q_R_B1_aby_T+q_LR_B1_aby+q_LR_B1_aby_T;
            adenosine_T = q_cAMP+q_PDE_cAMP+q_five_AMP+q_ATP+q_AC_ATP+q_a_Gs_GTP_AC_ATP+q_FSK_AC_ATP;

// 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_Ca_i: fmol {pub: out};
var q_Ca_sr: fmol {pub: out};
var q_MgATP: fmol {pub: out};
var q_MgADP: fmol {pub: out};
var q_Pi: fmol {pub: out};
var q_Ca_SR: fmol {pub: out};
var q_Ca_di: 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_Cai: 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};

// 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_Ca_i_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_Pi_mSERCA: fmol {pub: in};


// RyR imports
var q_Ca_SR_mRyR: fmol {pub: in};
var q_Ca_di_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};


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_Ca_i = q_Ca_i_init + q_Ca_i_mSERCA ;
q_Ca_sr = q_Ca_sr_init + q_Ca_sr_mSERCA ;
q_MgATP = q_MgATP_init + q_MgATP_mSERCA ;
q_MgADP = q_MgADP_init + q_MgADP_mSERCA ;
q_Pi = q_Pi_init + q_Pi_mSERCA ;
q_Ca_SR = q_Ca_SR_init + q_Ca_SR_mRyR ;
q_Ca_di = q_Ca_di_init + q_Ca_di_mRyR ;
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_Cai = q_Cai_init + q_Cai_mLCC ;
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 ;
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_Ca_i_mSERCA and q_Ca_i;
vars q_Ca_i and q_Ca_i_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_Pi_mSERCA and q_Pi;
vars q_Pi and q_Pi_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_Ca_di_mRyR and q_Ca_di;
vars q_Ca_di and q_Ca_di_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 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_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 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_di and K_Ca_di;
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 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;

enddef;