def model FCU_adenylylCyclase 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 "cAMP/BG_cAMP.cellml" for
comp cAMP using comp cAMP;
enddef;
def import using "GPCR_B1AR_reduced/BG_GPCR_B1AR_reduced.cellml" for
comp GPCR_B1AR_reduced using comp GPCR_B1AR_reduced;
enddef;
def import using "GPCR_M2_reduced/BG_GPCR_M2_reduced.cellml" for
comp GPCR_M2_reduced using comp GPCR_M2_reduced;
enddef;

def comp BG_parameters as
var kappa_1a: fmol_per_sec {init: 6.19718e+06, pub: out};
var kappa_1b: fmol_per_sec {init: 0.00129391, pub: out};
var kappa_2a: fmol_per_sec {init: 176114, pub: out};
var kappa_2b: fmol_per_sec {init: 0.0475228, pub: out};
var kappa_3a: fmol_per_sec {init: 2.64321e+08, pub: out};
var kappa_3b: fmol_per_sec {init: 3.0735e-17, pub: out};
var kappa_4a: fmol_per_sec {init: 30960.9, pub: out};
var kappa_4b: fmol_per_sec {init: 0.119081, pub: out};
var kappa_5: fmol_per_sec {init: 823.893, pub: out};
var kappa_6: fmol_per_sec {init: 3191.74, pub: out};
var kappa_7: fmol_per_sec {init: 479033, pub: out};
var kappa_GiAC: fmol_per_sec {init: 2417.47, pub: out};
var kappa_Rswitch_B1: fmol_per_sec {init: 1.02184, pub: out};
var kappa_LRswitch_B1: fmol_per_sec {init: 1.49856e-07, pub: out};
var kappa_C_B1: fmol_per_sec {init: 4082.77, pub: out};
var kappa_R_B1: fmol_per_sec {init: 5.98748e+06, pub: out};
var kappa_L_B1: fmol_per_sec {init: 4.27089e-08, pub: out};
var kappa_Act1_Gs: fmol_per_sec {init: 0.000215226, pub: out};
var kappa_Act2_Gs: fmol_per_sec {init: 167.999, pub: out};
var kappa_Hyd_Gs: fmol_per_sec {init: 0.00211433, pub: out};
var kappa_Reassoc_Gs: fmol_per_sec {init: 3.99549e-05, pub: out};
var kappa_InternR_B1: fmol_per_sec {init: 0.892127, pub: out};
var kappa_InternLR_B1: fmol_per_sec {init: 696367, pub: out};
var kappa_Rswitch_M2: fmol_per_sec {init: 241052, pub: out};
var kappa_LRswitch_M2: fmol_per_sec {init: 5.10328, pub: out};
var kappa_C_M2: fmol_per_sec {init: 97149.6, pub: out};
var kappa_R_M2: fmol_per_sec {init: 2.05674e+06, pub: out};
var kappa_L_M2: fmol_per_sec {init: 56.136, pub: out};
var kappa_Act1_Gi: fmol_per_sec {init: 0.00220056, pub: out};
var kappa_Act2_Gi: fmol_per_sec {init: 0.609877, pub: out};
var kappa_Hyd_Gi: fmol_per_sec {init: 1.60142, pub: out};
var kappa_Reassoc_Gi: fmol_per_sec {init: 4.03024e-05, pub: out};
var kappa_InternR_M2: fmol_per_sec {init: 3.61714e-08, pub: out};
var kappa_InternLR_M2: fmol_per_sec {init: 1.00248e-05, pub: out};
var K_ATP: per_fmol {init: 5.26836e-05, pub: out};
var K_cAMP: per_fmol {init: 0.0157529, pub: out};
var K_AC: per_fmol {init: 2.40711, pub: out};
var K_AC_ATP: per_fmol {init: 4.49332, pub: out};
var K_Gsa_GTP_AC: per_fmol {init: 9.10781, pub: out};
var K_Gsa_GTP_AC_ATP: per_fmol {init: 5.19946, pub: out};
var K_FSK_AC: per_fmol {init: 0.0606843, pub: out};
var K_FSK_AC_ATP: per_fmol {init: 0.0945821, pub: out};
var K_PDE: per_fmol {init: 1.73264, pub: out};
var K_PDE_cAMP: per_fmol {init: 1.22059, pub: out};
var K_five_AMP: per_fmol {init: 0.0157529, pub: out};
var K_IBMX: per_fmol {init: 0.0197325, pub: out};
var K_PDEinh: per_fmol {init: 35.2834, pub: out};
var K_Gsa_GTP: per_fmol {init: 0.274979, pub: out};
var K_FSK: per_fmol {init: 1.66559e-05, pub: out};
var K_Gia_GTP: per_fmol {init: 0.014522, pub: out};
var K_ACinh: per_fmol {init: 12.0249, pub: out};
var K_PPi: per_fmol {init: 9.722e-05, pub: out};
var K_RB1_inactive: per_fmol {init: 28448.3, pub: out};
var K_L_RB1_inactive: per_fmol {init: 1.93985e+11, pub: out};
var K_LB1: per_fmol {init: 695517, pub: out};
var K_RB1: per_fmol {init: 284.483, pub: out};
var K_Gs: per_fmol {init: 0.000727564, pub: out};
var K_RB1_Gs: per_fmol {init: 234.964, pub: out};
var K_L_RB1: per_fmol {init: 0.193985, pub: out};
var K_L_RB1_Gs: per_fmol {init: 0.000301016, pub: out};
var K_Gsbetagamma: per_fmol {init: 490.088, pub: out};
var K_Gsa_GDP: per_fmol {init: 52.2185, pub: out};
var K_GTP: per_fmol {init: 10.6947, pub: out};
var K_GDP: per_fmol {init: 3.83436e-09, pub: out};
var K_Pi: per_fmol {init: 7.65396e-09, pub: out};
var K_RB1_tag: per_fmol {init: 0.00642105, pub: out};
var K_L_RB1_tag: per_fmol {init: 8.22611e-09, pub: out};
var K_RB1_GRKArr: per_fmol {init: 3.25848e-08, pub: out};
var K_L_RB1_GRKArr: per_fmol {init: 4.17449e-14, pub: out};
var K_GRKArr: per_fmol {init: 147520, pub: out};
var K_RM2_inactive: per_fmol {init: 0.00120595, pub: out};
var K_L_RM2_inactive: per_fmol {init: 56.963, pub: out};
var K_LM2: per_fmol {init: 124.828, pub: out};
var K_RM2: per_fmol {init: 1.20595e-05, pub: out};
var K_Gi: per_fmol {init: 7.21292, pub: out};
var K_RM2_Gi: per_fmol {init: 0.0897681, pub: out};
var K_L_RM2: per_fmol {init: 5.6963e-07, pub: out};
var K_L_RM2_Gi: per_fmol {init: 0.000323902, pub: out};
var K_Gibetagamma: per_fmol {init: 36.7999, pub: out};
var K_Gia_GDP: per_fmol {init: 689.431, pub: out};
var K_RM2_tag: per_fmol {init: 1583.68, pub: out};
var K_L_RM2_tag: per_fmol {init: 5.71425, pub: out};
var K_RM2_GRKArr: per_fmol {init: 80.3667, pub: out};
var K_L_RM2_GRKArr: per_fmol {init: 0.28998, 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 stimSt1: second {init: 3.5e-4};
            var stimDur1: second {init: 0.25e-4};
            var tRamp1: second {init: 1.8e-4};
            var stimMag1: fmol {init: 1e1};
            var stimHolding1: fmol {init: 1e-5};
            var m1: fmol_per_sec;
            m1 = stimMag1/tRamp1;
            q_LB1 = sel
            case (time < stimSt1) and (time > stimSt1-tRamp1): 
            stimHolding1+m1*(time-stimSt1+tRamp1);
            case (time >= stimSt1) and (time < stimSt1+stimDur1): 
            stimMag1+stimHolding1;
            case (time < stimSt1+tRamp1+stimDur1) and (time >= stimSt1+stimDur1): 
            stimHolding1+-m1*(time-stimSt1-tRamp1-stimDur1);
            otherwise: 
            stimHolding1;
            endsel;
            var stimSt2: second {init: 7e-4};
            var stimDur2: second {init: 0.25e-4};
            var tRamp2: second {init: 1.8e-4};
            var stimMag2: fmol {init: 1e1};
            var stimHolding2: fmol {init: 1e-5};
            var m2: fmol_per_sec;
            m2 = stimMag2/tRamp2;
            q_LM2 = sel
            case (time < stimSt2) and (time > stimSt2-tRamp2): 
            stimHolding2+m2*(time-stimSt2+tRamp2);
            case (time >= stimSt2) and (time < stimSt2+stimDur2): 
            stimMag2+stimHolding2;
            case (time < stimSt2+tRamp2+stimDur2) and (time >= stimSt2+stimDur2): 
            stimHolding2+-m2*(time-stimSt2-tRamp2-stimDur2);
            otherwise: 
            stimHolding2;
            endsel;
var q_ATP: fmol {init: 1e-888, pub: out};
var q_cAMP: fmol {init: 1e-888, pub: out};
var q_AC: fmol {init: 1e-888, pub: out};
var q_AC_ATP: fmol {init: 1e-888, pub: out};
var q_Gsa_GTP_AC: fmol {init: 1e-888, pub: out};
var q_Gsa_GTP_AC_ATP: fmol {init: 1e-888, pub: out};
var q_FSK_AC: fmol {init: 1e-888, pub: out};
var q_FSK_AC_ATP: fmol {init: 1e-888, pub: out};
var q_PDE: fmol {init: 1e-888, pub: out};
var q_PDE_cAMP: fmol {init: 1e-888, pub: out};
var q_five_AMP: fmol {init: 1e-888, pub: out};
var q_IBMX: fmol {init: 1e-888, pub: out};
var q_PDEinh: fmol {init: 1e-888, pub: out};
var q_Gsa_GTP: fmol {init: 1e-888, pub: out};
var q_FSK: fmol {init: 1e-888, pub: out};
var q_Gia_GTP: fmol {init: 1e-888, pub: out};
var q_ACinh: fmol {init: 1e-888, pub: out};
var q_PPi: fmol {init: 1e-888, pub: out};
var q_RB1_inactive: fmol {init: 1e-888, pub: out};
var q_L_RB1_inactive: fmol {init: 1e-888, pub: out};
var q_LB1: fmol {init: 1e-888, pub: out};
var q_RB1: fmol {init: 1e-888, pub: out};
var q_Gs: fmol {init: 1e-888, pub: out};
var q_RB1_Gs: fmol {init: 1e-888, pub: out};
var q_L_RB1: fmol {init: 1e-888, pub: out};
var q_L_RB1_Gs: fmol {init: 1e-888, pub: out};
var q_Gsbetagamma: fmol {init: 1e-888, pub: out};
var q_Gsa_GDP: fmol {init: 1e-888, pub: out};
var q_GTP: fmol {init: 1e-888, pub: out};
var q_GDP: fmol {init: 1e-888, pub: out};
var q_Pi: fmol {init: 1e-888, pub: out};
var q_RB1_tag: fmol {init: 1e-888, pub: out};
var q_L_RB1_tag: fmol {init: 1e-888, pub: out};
var q_RB1_GRKArr: fmol {init: 1e-888, pub: out};
var q_L_RB1_GRKArr: fmol {init: 1e-888, pub: out};
var q_GRKArr: fmol {init: 1e-888, pub: out};
var q_RM2_inactive: fmol {init: 1e-888, pub: out};
var q_L_RM2_inactive: fmol {init: 1e-888, pub: out};
var q_LM2: fmol {init: 1e-888, pub: out};
var q_RM2: fmol {init: 1e-888, pub: out};
var q_Gi: fmol {init: 1e-888, pub: out};
var q_RM2_Gi: fmol {init: 1e-888, pub: out};
var q_L_RM2: fmol {init: 1e-888, pub: out};
var q_L_RM2_Gi: fmol {init: 1e-888, pub: out};
var q_Gibetagamma: fmol {init: 1e-888, pub: out};
var q_Gia_GDP: fmol {init: 1e-888, pub: out};
var q_RM2_tag: fmol {init: 1e-888, pub: out};
var q_L_RM2_tag: fmol {init: 1e-888, pub: out};
var q_RM2_GRKArr: fmol {init: 1e-888, pub: out};
var q_L_RM2_GRKArr: fmol {init: 1e-888, pub: out};

// mass conservation checks
 var LB1_T: fmol;
            var RB1_T: fmol;
            var Gs_T: fmol;
            var adenosine_T: fmol;
        LB1_T = q_L_RB1_inactive+q_LB1+q_L_RB1_Gs+q_L_RB1+q_L_RB1_tag+q_L_RB1_GRKArr;
            RB1_T = q_RB1_inactive+q_L_RB1_inactive+q_RB1+q_RB1_Gs+q_L_RB1+q_L_RB1_Gs +q_RB1_tag+q_L_RB1_tag+q_RB1_GRKArr+q_L_RB1_GRKArr;
            Gs_T = q_Gs+q_RB1_Gs+q_L_RB1_Gs+q_Gsa_GTP+q_Gsa_GDP;
            adenosine_T = q_cAMP+q_PDE_cAMP+q_five_AMP+q_ATP+q_AC_ATP+q_Gsa_GTP_AC_ATP+q_FSK_AC_ATP;

// cAMP imports
var v_1a: fmol_per_sec {pub: in};
var v_1b: fmol_per_sec {pub: in};
var v_2a: fmol_per_sec {pub: in};
var v_2b: fmol_per_sec {pub: in};
var v_3a: fmol_per_sec {pub: in};
var v_3b: fmol_per_sec {pub: in};
var v_4a: fmol_per_sec {pub: in};
var v_4b: fmol_per_sec {pub: in};
var v_5: fmol_per_sec {pub: in};
var v_6: fmol_per_sec {pub: in};
var v_7: fmol_per_sec {pub: in};
var v_GiAC: fmol_per_sec {pub: in};


// GPCR_B1AR_reduced imports
var v_Rswitch_B1: fmol_per_sec {pub: in};
var v_LRswitch_B1: fmol_per_sec {pub: in};
var v_C_B1: fmol_per_sec {pub: in};
var v_R_B1: fmol_per_sec {pub: in};
var v_L_B1: fmol_per_sec {pub: in};
var v_Act1_Gs: fmol_per_sec {pub: in};
var v_Act2_Gs: fmol_per_sec {pub: in};
var v_Hyd_Gs: fmol_per_sec {pub: in};
var v_Reassoc_Gs: fmol_per_sec {pub: in};
var v_InternR_B1: fmol_per_sec {pub: in};
var v_InternLR_B1: fmol_per_sec {pub: in};


// GPCR_M2_reduced imports
var v_Rswitch_M2: fmol_per_sec {pub: in};
var v_LRswitch_M2: fmol_per_sec {pub: in};
var v_C_M2: fmol_per_sec {pub: in};
var v_R_M2: fmol_per_sec {pub: in};
var v_L_M2: fmol_per_sec {pub: in};
var v_Act1_Gi: fmol_per_sec {pub: in};
var v_Act2_Gi: fmol_per_sec {pub: in};
var v_Hyd_Gi: fmol_per_sec {pub: in};
var v_Reassoc_Gi: fmol_per_sec {pub: in};
var v_InternR_M2: fmol_per_sec {pub: in};
var v_InternLR_M2: fmol_per_sec {pub: in};


ode(q_ATP, time) = vvv;
ode(q_cAMP, time) = vvv;
ode(q_AC, time) = vvv;
ode(q_AC_ATP, time) = vvv;
ode(q_Gsa_GTP_AC, time) = vvv;
ode(q_Gsa_GTP_AC_ATP, time) = vvv;
ode(q_FSK_AC, time) = vvv;
ode(q_FSK_AC_ATP, time) = vvv;
ode(q_PDE, time) = vvv;
ode(q_PDE_cAMP, time) = vvv;
ode(q_five_AMP, time) = vvv;
ode(q_IBMX, time) = vvv;
ode(q_PDEinh, time) = vvv;
ode(q_Gsa_GTP, time) = vvv;
ode(q_FSK, time) = vvv;
ode(q_Gia_GTP, time) = vvv;
ode(q_ACinh, time) = vvv;
ode(q_PPi, time) = vvv;
ode(q_RB1_inactive, time) = vvv;
ode(q_L_RB1_inactive, time) = vvv;
ode(q_LB1, time) = vvv;
ode(q_RB1, time) = vvv;
ode(q_Gs, time) = vvv;
ode(q_RB1_Gs, time) = vvv;
ode(q_L_RB1, time) = vvv;
ode(q_L_RB1_Gs, time) = vvv;
ode(q_Gsbetagamma, time) = vvv;
ode(q_Gsa_GDP, time) = vvv;
ode(q_GTP, time) = vvv;
ode(q_GDP, time) = vvv;
ode(q_Pi, time) = vvv;
ode(q_RB1_tag, time) = vvv;
ode(q_L_RB1_tag, time) = vvv;
ode(q_RB1_GRKArr, time) = vvv;
ode(q_L_RB1_GRKArr, time) = vvv;
ode(q_GRKArr, time) = vvv;
ode(q_RM2_inactive, time) = vvv;
ode(q_L_RM2_inactive, time) = vvv;
ode(q_LM2, time) = vvv;
ode(q_RM2, time) = vvv;
ode(q_Gi, time) = vvv;
ode(q_RM2_Gi, time) = vvv;
ode(q_L_RM2, time) = vvv;
ode(q_L_RM2_Gi, time) = vvv;
ode(q_Gibetagamma, time) = vvv;
ode(q_Gia_GDP, time) = vvv;
ode(q_RM2_tag, time) = vvv;
ode(q_L_RM2_tag, time) = vvv;
ode(q_RM2_GRKArr, time) = vvv;
ode(q_L_RM2_GRKArr, time) = vvv;
enddef;

def map between environment and cAMP for
vars time and time;
vars q_ATP and q_ATP;
vars q_cAMP and q_cAMP;
vars q_AC and q_AC;
vars q_AC_ATP and q_AC_ATP;
vars q_Gsa_GTP_AC and q_Gsa_GTP_AC;
vars q_Gsa_GTP_AC_ATP and q_Gsa_GTP_AC_ATP;
vars q_FSK_AC and q_FSK_AC;
vars q_FSK_AC_ATP and q_FSK_AC_ATP;
vars q_PDE and q_PDE;
vars q_PDE_cAMP and q_PDE_cAMP;
vars q_five_AMP and q_five_AMP;
vars q_IBMX and q_IBMX;
vars q_PDEinh and q_PDEinh;
vars q_Gsa_GTP and q_Gsa_GTP;
vars q_FSK and q_FSK;
vars q_Gia_GTP and q_Gia_GTP;
vars q_ACinh and q_ACinh;
vars q_PPi and q_PPi;
vars v_1a and v_1a;
vars v_1b and v_1b;
vars v_2a and v_2a;
vars v_2b and v_2b;
vars v_3a and v_3a;
vars v_3b and v_3b;
vars v_4a and v_4a;
vars v_4b and v_4b;
vars v_5 and v_5;
vars v_6 and v_6;
vars v_7 and v_7;
vars v_GiAC and v_GiAC;
enddef;

def map between environment and GPCR_B1AR_reduced for
vars time and time;
vars q_RB1_inactive and q_RB1_inactive;
vars q_L_RB1_inactive and q_L_RB1_inactive;
vars q_LB1 and q_LB1;
vars q_RB1 and q_RB1;
vars q_Gs and q_Gs;
vars q_RB1_Gs and q_RB1_Gs;
vars q_L_RB1 and q_L_RB1;
vars q_L_RB1_Gs and q_L_RB1_Gs;
vars q_Gsa_GTP and q_Gsa_GTP;
vars q_Gsbetagamma and q_Gsbetagamma;
vars q_Gsa_GDP and q_Gsa_GDP;
vars q_GTP and q_GTP;
vars q_GDP and q_GDP;
vars q_Pi and q_Pi;
vars q_RB1_tag and q_RB1_tag;
vars q_L_RB1_tag and q_L_RB1_tag;
vars q_RB1_GRKArr and q_RB1_GRKArr;
vars q_L_RB1_GRKArr and q_L_RB1_GRKArr;
vars q_GRKArr and q_GRKArr;
vars v_Rswitch_B1 and v_Rswitch_B1;
vars v_LRswitch_B1 and v_LRswitch_B1;
vars v_C_B1 and v_C_B1;
vars v_R_B1 and v_R_B1;
vars v_L_B1 and v_L_B1;
vars v_Act1_Gs and v_Act1_Gs;
vars v_Act2_Gs and v_Act2_Gs;
vars v_Hyd_Gs and v_Hyd_Gs;
vars v_Reassoc_Gs and v_Reassoc_Gs;
vars v_InternR_B1 and v_InternR_B1;
vars v_InternLR_B1 and v_InternLR_B1;
enddef;

def map between environment and GPCR_M2_reduced for
vars time and time;
vars q_RM2_inactive and q_RM2_inactive;
vars q_L_RM2_inactive and q_L_RM2_inactive;
vars q_LM2 and q_LM2;
vars q_RM2 and q_RM2;
vars q_Gi and q_Gi;
vars q_RM2_Gi and q_RM2_Gi;
vars q_L_RM2 and q_L_RM2;
vars q_L_RM2_Gi and q_L_RM2_Gi;
vars q_Gia_GTP and q_Gia_GTP;
vars q_Gibetagamma and q_Gibetagamma;
vars q_Gia_GDP and q_Gia_GDP;
vars q_GTP and q_GTP;
vars q_GDP and q_GDP;
vars q_Pi and q_Pi;
vars q_RM2_tag and q_RM2_tag;
vars q_L_RM2_tag and q_L_RM2_tag;
vars q_RM2_GRKArr and q_RM2_GRKArr;
vars q_L_RM2_GRKArr and q_L_RM2_GRKArr;
vars q_GRKArr and q_GRKArr;
vars v_Rswitch_M2 and v_Rswitch_M2;
vars v_LRswitch_M2 and v_LRswitch_M2;
vars v_C_M2 and v_C_M2;
vars v_R_M2 and v_R_M2;
vars v_L_M2 and v_L_M2;
vars v_Act1_Gi and v_Act1_Gi;
vars v_Act2_Gi and v_Act2_Gi;
vars v_Hyd_Gi and v_Hyd_Gi;
vars v_Reassoc_Gi and v_Reassoc_Gi;
vars v_InternR_M2 and v_InternR_M2;
vars v_InternLR_M2 and v_InternLR_M2;
enddef;

def map between BG_parameters and cAMP for
vars kappa_1a and kappa_1a;
vars kappa_1b and kappa_1b;
vars kappa_2a and kappa_2a;
vars kappa_2b and kappa_2b;
vars kappa_3a and kappa_3a;
vars kappa_3b and kappa_3b;
vars kappa_4a and kappa_4a;
vars kappa_4b and kappa_4b;
vars kappa_5 and kappa_5;
vars kappa_6 and kappa_6;
vars kappa_7 and kappa_7;
vars kappa_GiAC and kappa_GiAC;
vars K_ATP and K_ATP;
vars K_cAMP and K_cAMP;
vars K_AC and K_AC;
vars K_AC_ATP and K_AC_ATP;
vars K_Gsa_GTP_AC and K_Gsa_GTP_AC;
vars K_Gsa_GTP_AC_ATP and K_Gsa_GTP_AC_ATP;
vars K_FSK_AC and K_FSK_AC;
vars K_FSK_AC_ATP and K_FSK_AC_ATP;
vars K_PDE and K_PDE;
vars K_PDE_cAMP and K_PDE_cAMP;
vars K_five_AMP and K_five_AMP;
vars K_IBMX and K_IBMX;
vars K_PDEinh and K_PDEinh;
vars K_Gsa_GTP and K_Gsa_GTP;
vars K_FSK and K_FSK;
vars K_Gia_GTP and K_Gia_GTP;
vars K_ACinh and K_ACinh;
vars K_PPi and K_PPi;
enddef;
def map between BG_parameters and GPCR_B1AR_reduced for
vars kappa_Rswitch_B1 and kappa_Rswitch_B1;
vars kappa_LRswitch_B1 and kappa_LRswitch_B1;
vars kappa_C_B1 and kappa_C_B1;
vars kappa_R_B1 and kappa_R_B1;
vars kappa_L_B1 and kappa_L_B1;
vars kappa_Act1_Gs and kappa_Act1_Gs;
vars kappa_Act2_Gs and kappa_Act2_Gs;
vars kappa_Hyd_Gs and kappa_Hyd_Gs;
vars kappa_Reassoc_Gs and kappa_Reassoc_Gs;
vars kappa_InternR_B1 and kappa_InternR_B1;
vars kappa_InternLR_B1 and kappa_InternLR_B1;
vars K_RB1_inactive and K_RB1_inactive;
vars K_L_RB1_inactive and K_L_RB1_inactive;
vars K_LB1 and K_LB1;
vars K_RB1 and K_RB1;
vars K_Gs and K_Gs;
vars K_RB1_Gs and K_RB1_Gs;
vars K_L_RB1 and K_L_RB1;
vars K_L_RB1_Gs and K_L_RB1_Gs;
vars K_Gsa_GTP and K_Gsa_GTP;
vars K_Gsbetagamma and K_Gsbetagamma;
vars K_Gsa_GDP and K_Gsa_GDP;
vars K_GTP and K_GTP;
vars K_GDP and K_GDP;
vars K_Pi and K_Pi;
vars K_RB1_tag and K_RB1_tag;
vars K_L_RB1_tag and K_L_RB1_tag;
vars K_RB1_GRKArr and K_RB1_GRKArr;
vars K_L_RB1_GRKArr and K_L_RB1_GRKArr;
vars K_GRKArr and K_GRKArr;
enddef;
def map between BG_parameters and GPCR_M2_reduced for
vars kappa_Rswitch_M2 and kappa_Rswitch_M2;
vars kappa_LRswitch_M2 and kappa_LRswitch_M2;
vars kappa_C_M2 and kappa_C_M2;
vars kappa_R_M2 and kappa_R_M2;
vars kappa_L_M2 and kappa_L_M2;
vars kappa_Act1_Gi and kappa_Act1_Gi;
vars kappa_Act2_Gi and kappa_Act2_Gi;
vars kappa_Hyd_Gi and kappa_Hyd_Gi;
vars kappa_Reassoc_Gi and kappa_Reassoc_Gi;
vars kappa_InternR_M2 and kappa_InternR_M2;
vars kappa_InternLR_M2 and kappa_InternLR_M2;
vars K_RM2_inactive and K_RM2_inactive;
vars K_L_RM2_inactive and K_L_RM2_inactive;
vars K_LM2 and K_LM2;
vars K_RM2 and K_RM2;
vars K_Gi and K_Gi;
vars K_RM2_Gi and K_RM2_Gi;
vars K_L_RM2 and K_L_RM2;
vars K_L_RM2_Gi and K_L_RM2_Gi;
vars K_Gia_GTP and K_Gia_GTP;
vars K_Gibetagamma and K_Gibetagamma;
vars K_Gia_GDP and K_Gia_GDP;
vars K_GTP and K_GTP;
vars K_GDP and K_GDP;
vars K_Pi and K_Pi;
vars K_RM2_tag and K_RM2_tag;
vars K_L_RM2_tag and K_L_RM2_tag;
vars K_RM2_GRKArr and K_RM2_GRKArr;
vars K_L_RM2_GRKArr and K_L_RM2_GRKArr;
vars K_GRKArr and K_GRKArr;
enddef;

def map between constants and cAMP for
	vars R and R;
	vars T and T;
enddef;
def map between constants and GPCR_B1AR_reduced for
	vars R and R;
	vars T and T;
enddef;
def map between constants and GPCR_M2_reduced for
	vars R and R;
	vars T and T;
enddef;

enddef;