Location: BG_CaL @ 58424fcf9c90 / pan_code_LCC.txt

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2022-02-22 14:14:19+13:00
Desc:
Direct from Pan code
Permanent Source URI:
https://models.cellml.org/workspace/6d7/rawfile/58424fcf9c90f67e489c1f1e8b1cc20ec938dd58/pan_code_LCC.txt

def model cardiac_AP as
    def unit J_per_K_per_mol as
        unit joule;
        unit kelvin {expo: -1};
        unit mole {expo: -1};
    enddef;

    def unit C_per_mol as
        unit coulomb;
        unit mole {expo: -1};
    enddef;

    def unit fF as
        unit farad {pref: femto};
    enddef;

    def unit fmol as
        unit mole {pref: femto};
    enddef;

    def unit per_fmol as
        unit fmol {expo: -1};
    enddef;

    def unit fmol_per_sec as
        unit fmol;
        unit second {expo: -1};
    enddef;

    def unit fC as
        unit coulomb {pref: femto};
    enddef;

    def unit fA as
        unit ampere {pref: femto};
    enddef;

    def unit J_per_mol as
        unit joule;
        unit mole {expo: -1};
    enddef;

    def unit volt_per_sec as
        unit volt;
        unit second {expo: -1};
    enddef;

    def unit mM as
        unit mole {pref: milli};
        unit litre {expo: -1};
    enddef;

    def comp environment as
        // Time

        var t: second {init: 0};

        // Physical parameters

        var R: J_per_K_per_mol {init: 8.314};
        var T: kelvin {init: 310};
        var F: C_per_mol {init: 96485};

        // Bond graph parameters

        var C_m: fF {init: 153400};
        var K_000_LCC: per_fmol {init: 1.3174};
        var K_000_Na: per_fmol {init: 0.0073376};
        var K_001_LCC: per_fmol {init: 1565.0581};
        var K_001_Na: per_fmol {init: 23718.0222};
        var K_00_K: per_fmol {init: 128.6868};
        var K_010_LCC: per_fmol {init: 53.3957};
        var K_010_Na: per_fmol {init: 17499.8767};
        var K_011_LCC: per_fmol {init: 63435.6844};
        var K_011_Na: per_fmol {init: 56566767727.4036};
        var K_01_K: per_fmol {init: 22.3033};
        var K_020_LCC: per_fmol {init: 28.5543};
        var K_021_LCC: per_fmol {init: 33923.3476};
        var K_100_LCC: per_fmol {init: 0.26529};
        var K_100_Na: per_fmol {init: 1.5633e-05};
        var K_101_LCC: per_fmol {init: 315.1706};
        var K_101_Na: per_fmol {init: 50.5317};
        var K_10_K: per_fmol {init: 16.5777};
        var K_10_NaK: per_fmol {init: 5076.1941};
        var K_110_LCC: per_fmol {init: 10.7528};
        var K_110_Na: per_fmol {init: 37.2839};
        var K_111_LCC: per_fmol {init: 12774.6463};
        var K_111_Na: per_fmol {init: 120516694.0958};
        var K_11_K: per_fmol {init: 2.8732};
        var K_11_NaK: per_fmol {init: 10.4383};
        var K_120_LCC: per_fmol {init: 5.7502};
        var K_121_LCC: per_fmol {init: 6831.4667};
        var K_12_NaK: per_fmol {init: 1.7361};
        var K_13_NaK: per_fmol {init: 1.155};
        var K_14_NaK: per_fmol {init: 319.1183};
        var K_15_NaK: per_fmol {init: 6250238.4012};
        var K_1_NCX: per_fmol {init: 160144852.2746};
        var K_1_NaK: per_fmol {init: 321178904.6865};
        var K_200_Na: per_fmol {init: 9.9918e-08};
        var K_201_Na: per_fmol {init: 0.32298};
        var K_20_K: per_fmol {init: 8.5423};
        var K_210_Na: per_fmol {init: 0.2383};
        var K_211_Na: per_fmol {init: 770290.0911};
        var K_21_K: per_fmol {init: 1.4805};
        var K_2_NCX: per_fmol {init: 407.3479};
        var K_2_NaK: per_fmol {init: 31093.1227};
        var K_300_Na: per_fmol {init: 1.9159e-09};
        var K_301_Na: per_fmol {init: 0.006193};
        var K_310_Na: per_fmol {init: 0.0045694};
        var K_311_Na: per_fmol {init: 14770.0739};
        var K_3_NCX: per_fmol {init: 22191.7214};
        var K_3_NaK: per_fmol {init: 12.0404};
        var K_4_NCX: per_fmol {init: 119580.5289};
        var K_4_NaK: per_fmol {init: 3.1105};
        var K_5_NCX: per_fmol {init: 10363474.1803};
        var K_5_NaK: per_fmol {init: 3.2141};
        var K_6_NCX: per_fmol {init: 199.44};
        var K_6_NaK: per_fmol {init: 170.6981};
        var K_7_NaK: per_fmol {init: 10036.4321};
        var K_8_NaK: per_fmol {init: 31550.6969};
        var K_9_NaK: per_fmol {init: 9874266.4618};
        var K_CMDN: per_fmol {init: 22.0647};
        var K_CMDNCa: per_fmol {init: 0.044337};
        var K_C_K1: per_fmol {init: 1.329};
        var K_C_Kp: per_fmol {init: 41.6265};
        var K_Cae: per_fmol {init: 0.16293};
        var K_Cai: per_fmol {init: 0.022218};
        var K_H: per_fmol {init: 0.060187};
        var K_Ke: per_fmol {init: 0.059344};
        var K_Ki: per_fmol {init: 0.0080926};
        var K_MgADP: per_fmol {init: 7.1044e-05};
        var K_MgATP: per_fmol {init: 3.6718};
        var K_Nae: per_fmol {init: 0.017678};
        var K_Nai: per_fmol {init: 0.0024108};
        var K_O_K1: per_fmol {init: 15967.3543};
        var K_O_Kp: per_fmol {init: 145.607};
        var K_P: per_fmol {init: 0.060187};
        var K_TRPN: per_fmol {init: 37.1165};
        var K_TRPNCa: per_fmol {init: 0.015668};
        var kappa_10_NaK: fmol_per_sec {init: 1045757.658};
        var kappa_11_NaK: fmol_per_sec {init: 311526.5915};
        var kappa_12_NaK: fmol_per_sec {init: 1873089.8475};
        var kappa_13_NaK: fmol_per_sec {init: 168.3622};
        var kappa_14_NaK: fmol_per_sec {init: 22.4585};
        var kappa_15_NaK: fmol_per_sec {init: 0.0049006};
        var kappa_1_NCX: fmol_per_sec {init: 256862.245};
        var kappa_1_NaK: fmol_per_sec {init: 104.5833};
        var kappa_2_NCX: fmol_per_sec {init: 21141816.1231};
        var kappa_2_NaK: fmol_per_sec {init: 270075.4148};
        var kappa_3_NCX: fmol_per_sec {init: 6834.2872};
        var kappa_3_NaK: fmol_per_sec {init: 906610.9954};
        var kappa_4_NCX: fmol_per_sec {init: 21141816.1231};
        var kappa_4_NaK: fmol_per_sec {init: 3509448.0698};
        var kappa_5_NCX: fmol_per_sec {init: 256862.245};
        var kappa_5_NaK: fmol_per_sec {init: 5858.2952};
        var kappa_6_NCX: fmol_per_sec {init: 1.7753};
        var kappa_6_NaK: fmol_per_sec {init: 8.3375};
        var kappa_7_NaK: fmol_per_sec {init: 1.1523};
        var kappa_8_NaK: fmol_per_sec {init: 1.1055};
        var kappa_9_NaK: fmol_per_sec {init: 2150.4256};
        var kappa_K: fmol_per_sec {init: 43947.9368};
        var kappa_K1: fmol_per_sec {init: 25.0154};
        var kappa_Kp: fmol_per_sec {init: 252.3419};
        var kappa_LCC_Ca1: fmol_per_sec {init: 31.5432};
        var kappa_LCC_Ca2: fmol_per_sec {init: 58.9849};
        var kappa_LCC_K1: fmol_per_sec {init: 0.068076};
        var kappa_LCC_K2: fmol_per_sec {init: 0.1273};
        var kappa_Na: fmol_per_sec {init: 32.858};
        var kappa_X00: fmol_per_sec {init: 0.034684};
        var kappa_X01: fmol_per_sec {init: 0.20012};
        var kappa_X10: fmol_per_sec {init: 0.13462};
        var kappa_X11: fmol_per_sec {init: 0.77675};
        var kappa_b1: fmol_per_sec {init: 839.7697};
        var kappa_b2: fmol_per_sec {init: 1412.6324};
        var kappa_d000: fmol_per_sec {init: 369.4994};
        var kappa_d001: fmol_per_sec {init: 0.31102};
        var kappa_d010: fmol_per_sec {init: 9.1161};
        var kappa_d011: fmol_per_sec {init: 0.0076733};
        var kappa_d020: fmol_per_sec {init: 17.0469};
        var kappa_d021: fmol_per_sec {init: 0.014349};
        var kappa_f1_000: fmol_per_sec {init: 0.68071};
        var kappa_f1_001: fmol_per_sec {init: 0.00057297};
        var kappa_f1_100: fmol_per_sec {init: 3.3802};
        var kappa_f1_101: fmol_per_sec {init: 0.0028452};
        var kappa_f2_000: fmol_per_sec {init: 1.2467};
        var kappa_f2_001: fmol_per_sec {init: 0.0010494};
        var kappa_f2_100: fmol_per_sec {init: 6.1907};
        var kappa_f2_101: fmol_per_sec {init: 0.0052109};
        var kappa_f3_010: fmol_per_sec {init: 3502.1028};
        var kappa_f3_011: fmol_per_sec {init: 2.9478};
        var kappa_f3_110: fmol_per_sec {init: 17390.563};
        var kappa_f3_111: fmol_per_sec {init: 14.6381};
        var kappa_fCa00: fmol_per_sec {init: 106492.3187};
        var kappa_fCa01: fmol_per_sec {init: 2627.333};
        var kappa_fCa02: fmol_per_sec {init: 4913.0371};
        var kappa_fCa10: fmol_per_sec {init: 528814.1044};
        var kappa_fCa11: fmol_per_sec {init: 13046.6757};
        var kappa_fCa12: fmol_per_sec {init: 24396.9083};
        var kappa_gK1: fmol_per_sec {init: 848.2812};
        var kappa_gKp: fmol_per_sec {init: 24.0195};
        var kappa_h000: fmol_per_sec {init: 0.045709};
        var kappa_h001: fmol_per_sec {init: 1.4141e-08};
        var kappa_h100: fmol_per_sec {init: 21.4544};
        var kappa_h101: fmol_per_sec {init: 6.6373e-06};
        var kappa_h200: fmol_per_sec {init: 3356.675};
        var kappa_h201: fmol_per_sec {init: 0.0010384};
        var kappa_h300: fmol_per_sec {init: 175057.5863};
        var kappa_h301: fmol_per_sec {init: 0.054157};
        var kappa_i00: fmol_per_sec {init: 7.7389};
        var kappa_i10: fmol_per_sec {init: 60.0742};
        var kappa_i20: fmol_per_sec {init: 116.5836};
        var kappa_j000: fmol_per_sec {init: 0.017824};
        var kappa_j010: fmol_per_sec {init: 7.4736e-09};
        var kappa_j100: fmol_per_sec {init: 8.3662};
        var kappa_j110: fmol_per_sec {init: 3.5079e-06};
        var kappa_j200: fmol_per_sec {init: 1308.9427};
        var kappa_j210: fmol_per_sec {init: 0.00054883};
        var kappa_j300: fmol_per_sec {init: 68264.0876};
        var kappa_j310: fmol_per_sec {init: 0.028623};
        var kappa_m000: fmol_per_sec {init: 5117405.6443};
        var kappa_m001: fmol_per_sec {init: 1.5832};
        var kappa_m010: fmol_per_sec {init: 2.1457};
        var kappa_m011: fmol_per_sec {init: 6.6381e-07};
        var kappa_m100: fmol_per_sec {init: 1601300127.6001};
        var kappa_m101: fmol_per_sec {init: 495.389};
        var kappa_m110: fmol_per_sec {init: 671.4132};
        var kappa_m111: fmol_per_sec {init: 0.00020771};
        var kappa_m200: fmol_per_sec {init: 125266701374.7363};
        var kappa_m201: fmol_per_sec {init: 38753.3515};
        var kappa_m210: fmol_per_sec {init: 52523.3901};
        var kappa_m211: fmol_per_sec {init: 0.016249};
        var nK: dimensionless {init: 0.5};
        var nK_stim: dimensionless;
        var nNae_NCX: dimensionless {init: 3};
        var nNai_NCX: dimensionless {init: 3};
        var zCa: dimensionless {init: 2};
        var zF_5: C_per_mol {init: -5301.979};
        var zF_8: C_per_mol {init: -91183.021};
        var zK: dimensionless {init: 1};
        var zNa: dimensionless {init: 1};
        var z_fK1: dimensionless {init: 0.033569};
        var z_fX: dimensionless {init: 0.51923};
        var z_fXi: dimensionless {init: 0};
        var z_f_Kp: dimensionless {init: 0};
        var z_fd: dimensionless {init: 2.1404};
        var z_ff1: dimensionless {init: -1.1495};
        var z_ff2: dimensionless {init: 0.72162};
        var z_ff3: dimensionless {init: 4.2933};
        var z_fh: dimensionless {init: -4.1892};
        var z_fj: dimensionless {init: -4.0381};
        var z_fm: dimensionless {init: 0.49541};
        var z_rCa: dimensionless {init: 2};
        var z_rK1: dimensionless {init: 3.1153};
        var z_rX: dimensionless {init: -0.73169};
        var z_rXi: dimensionless {init: 0.83216};
        var z_r_Kp: dimensionless {init: -4.4669};
        var z_rd: dimensionless {init: -2.1404};
        var z_rf1: dimensionless {init: 1.8993};
        var z_rf2: dimensionless {init: -0.52288};
        var z_rf3: dimensionless {init: 0};
        var z_rh: dimensionless {init: 1.2995};
        var z_rj: dimensionless {init: 1.4281};
        var z_rm: dimensionless {init: -2.4284};
        var zf_NCX: dimensionless {init: -0.55679};
        var zr_NCX: dimensionless {init: 0.44321};

        // Inputs

        var I_stim: fA;

        nK_stim = F;

        var stimPeriod: second {init: 1};
        var stimDuration: second {init: 0.0001};
        var tPeriod: second;

        tPeriod = t-floor(t/stimPeriod)*stimPeriod;

        I_stim = sel
            case (tPeriod >= 0.3{second}) and (tPeriod <= 0.3{second}+stimDuration):
                0.03{volt}*C_m/stimDuration;
            otherwise:
                0.0{fA};
        endsel;

        // States (and initial conditions)

        var x_Cae: fmol {init: 9.3276};
        var x_Cai: fmol {init: 0.00456};
        var x_Ke: fmol {init: 27.9828};
        var x_Ki: fmol {init: 5510};
        var x_Nae: fmol {init: 725.48};
        var x_Nai: fmol {init: 380};
        var q_mem: fC {init: -13039};
        var x_CMDN: fmol {init: 1.8088};
        var x_CMDNCa: fmol {init: 0.0912};
        var x_TRPN: fmol {init: 2.1452};
        var x_TRPNCa: fmol {init: 0.51484};
        var x_P1_NCX: fmol {init: 0.00072174};
        var x_P2_NCX: fmol {init: 0.00072174};
        var x_P3_NCX: fmol {init: 0.00072174};
        var x_P4_NCX: fmol {init: 0.00072174};
        var x_P5_NCX: fmol {init: 0.00072174};
        var x_P6_NCX: fmol {init: 0.00072174};
        var x_H: fmol {init: 0.0030534};
        var x_MgADP: fmol {init: 1.33};
        var x_MgATP: fmol {init: 264.1};
        var x_P7: fmol {init: 0.0078543};
        var x_P_i: fmol {init: 15.0898};
        var x_P14: fmol {init: 0.0078543};
        var x_P15: fmol {init: 0.0078543};
        var x_P1: fmol {init: 0.0078543};
        var x_P2: fmol {init: 0.0078543};
        var x_P3: fmol {init: 0.0078543};
        var x_P4: fmol {init: 0.0078543};
        var x_P5: fmol {init: 0.0078543};
        var x_P6: fmol {init: 0.0078543};
        var x_P10: fmol {init: 0.0078543};
        var x_P11: fmol {init: 0.0078543};
        var x_P12: fmol {init: 0.0078543};
        var x_P13: fmol {init: 0.0078543};
        var x_P8: fmol {init: 0.0078543};
        var x_P9: fmol {init: 0.0078543};
        var x_C_Kp: fmol {init: 9.6314e-07};
        var x_O_Kp: fmol {init: 2.4078e-07};
        var x_C_K1: fmol {init: 5.6606e-06};
        var x_O_K1: fmol {init: 1.4151e-06};
        var x_S00_K: fmol {init: 4.5648e-06};
        var x_S01_K: fmol {init: 1.1412e-06};
        var x_S10_K: fmol {init: 2.2824e-06};
        var x_S11_K: fmol {init: 5.706e-07};
        var x_S20_K: fmol {init: 2.853e-07};
        var x_S21_K: fmol {init: 7.1325e-08};
        var x_S000_LCC: fmol {init: 3.6989e-07};
        var x_S010_LCC: fmol {init: 7.3239e-05};
        var x_S020_LCC: fmol {init: 3.6989e-07};
        var x_S100_LCC: fmol {init: 3.7363e-09};
        var x_S110_LCC: fmol {init: 7.3979e-07};
        var x_S120_LCC: fmol {init: 3.7363e-09};
        var x_S001_LCC: fmol {init: 4.1099e-08};
        var x_S011_LCC: fmol {init: 8.1377e-06};
        var x_S021_LCC: fmol {init: 4.1099e-08};
        var x_S101_LCC: fmol {init: 4.1514e-10};
        var x_S111_LCC: fmol {init: 8.2199e-08};
        var x_S121_LCC: fmol {init: 4.1514e-10};
        var x_S000_Na: fmol {init: 0.00012033};
        var x_S010_Na: fmol {init: 1.337e-05};
        var x_S100_Na: fmol {init: 4.0111e-05};
        var x_S110_Na: fmol {init: 4.4568e-06};
        var x_S200_Na: fmol {init: 4.4568e-06};
        var x_S210_Na: fmol {init: 4.952e-07};
        var x_S300_Na: fmol {init: 1.6507e-07};
        var x_S310_Na: fmol {init: 1.8341e-08};
        var x_S001_Na: fmol {init: 1.337e-05};
        var x_S011_Na: fmol {init: 1.4856e-06};
        var x_S101_Na: fmol {init: 4.4568e-06};
        var x_S111_Na: fmol {init: 4.952e-07};
        var x_S201_Na: fmol {init: 4.952e-07};
        var x_S211_Na: fmol {init: 5.5022e-08};
        var x_S301_Na: fmol {init: 1.8341e-08};
        var x_S311_Na: fmol {init: 2.0379e-09};
        // Constitutive equations
        // C:Cae

        var mu_Cae: J_per_mol;
        var v_Cae: fmol_per_sec;
        v_Cae = f_22;
        ode(x_Cae, t) = v_Cae;
        mu_Cae = R*T*ln(K_Cae*x_Cae);

        // C:Cai

        var mu_Cai: J_per_mol;
        var v_Cai: fmol_per_sec;
        v_Cai = f_21;
        ode(x_Cai, t) = v_Cai;
        mu_Cai = R*T*ln(K_Cai*x_Cai);

        // C:Ke

        var mu_Ke: J_per_mol;
        var v_Ke: fmol_per_sec;
        v_Ke = f_7;
        ode(x_Ke, t) = v_Ke;
        mu_Ke = R*T*ln(K_Ke*x_Ke);

        // C:Ki

        var mu_Ki: J_per_mol;
        var v_Ki: fmol_per_sec;
        v_Ki = f_6;
        ode(x_Ki, t) = v_Ki;
        mu_Ki = R*T*ln(K_Ki*x_Ki);

        // C:Nae

        var mu_Nae: J_per_mol;
        var v_Nae: fmol_per_sec;

        e_1 = mu_Nae;
        v_Nae = f_1;
        ode(x_Nae, t) = v_Nae;
        mu_Nae = R*T*ln(K_Nae*x_Nae);

        // C:Nai

        var mu_Nai: J_per_mol;
        var v_Nai: fmol_per_sec;
        v_Nai = f_15;
        ode(x_Nai, t) = v_Nai;
        mu_Nai = R*T*ln(K_Nai*x_Nai);

        // C:mem

        var V_mem: volt;
        var I_mem: fA;
        I_mem = f_11;
        ode(q_mem, t) = I_mem;
        V_mem = q_mem/C_m;

//            f_11 = f_4+f_10+f_12+f_18+f_23+f_30;
//            f4    v_Na
//            f10   v_K
//            f_12	vLCC
//            f_18	v_NaK
//            f_23	NCX
//            f_30	I_stim

        vLCC = f_60-f_62-f_63;

        f_60 = 
        v_f2_001*(z_rf2-z_ff2) 
        +v_f3_011*(z_rf3-z_ff3)
        +v_f1_101*(z_rf1-z_ff1)
        +v_f2_101*(z_rf2-z_ff2)
        +v_f3_111*(z_rf3-z_ff3)
        +v_f1_001*(z_rf1-z_ff1)
        +v_d021*(z_rd-z_fd)
        +v_d011*(z_rd-z_fd)
        +v_d001*(z_rd-z_fd)
        +v_f3_110*(z_rf3-z_ff3)
        +v_f2_100*(z_rf2-z_ff2)
        +v_f1_100*(z_rf1-z_ff1)
        +v_f3_010*(z_rf3-z_ff3)
        +v_f2_000*(z_rf2-z_ff2)
        +v_f1_000*(z_rf1-z_ff1)
        +v_d020*(z_rd-z_fd)
        +v_d010*(z_rd-z_fd)
        +v_d000*(z_rd-z_fd)

        f_62 = zCa*(v_LCC_Ca1) + zCa*(v_LCC_Ca2);

        f_63 = zK*(v_LCC_K1) +zK*(v_LCC_K2);

        // C:CMDN

        var mu_CMDN: J_per_mol;
        var v_CMDN: fmol_per_sec;

        e_36 = mu_CMDN;
        v_CMDN = -f_36;
        ode(x_CMDN, t) = v_CMDN;
        mu_CMDN = R*T*ln(K_CMDN*x_CMDN);

        // C:CMDNCa

        var mu_CMDNCa: J_per_mol;
        var v_CMDNCa: fmol_per_sec;

        e_42 = mu_CMDNCa;
        v_CMDNCa = f_42;
        ode(x_CMDNCa, t) = v_CMDNCa;
        mu_CMDNCa = R*T*ln(K_CMDNCa*x_CMDNCa);

        // C:TRPN

        var mu_TRPN: J_per_mol;
        var v_TRPN: fmol_per_sec;

        e_35 = mu_TRPN;
        v_TRPN = -f_35;
        ode(x_TRPN, t) = v_TRPN;
        mu_TRPN = R*T*ln(K_TRPN*x_TRPN);

        // C:TRPNCa

        var mu_TRPNCa: J_per_mol;
        var v_TRPNCa: fmol_per_sec;

        e_41 = mu_TRPNCa;
        v_TRPNCa = f_41;
        ode(x_TRPNCa, t) = v_TRPNCa;
        mu_TRPNCa = R*T*ln(K_TRPNCa*x_TRPNCa);

        // Re:b1

        var Af_b1: J_per_mol;
        var Ar_b1: J_per_mol;
        var v_b1: fmol_per_sec;

        Af_b1 = e_37;
        Ar_b1 = e_39;
        f_37 = v_b1;
        f_39 = v_b1;
        v_b1 = kappa_b1*(exp(Af_b1/(R*T))-exp(Ar_b1/(R*T)));

        // Re:b2

        var Af_b2: J_per_mol;
        var Ar_b2: J_per_mol;
        var v_b2: fmol_per_sec;

        Af_b2 = e_38;
        Ar_b2 = e_40;
        f_38 = v_b2;
        f_40 = v_b2;
        v_b2 = kappa_b2*(exp(Af_b2/(R*T))-exp(Ar_b2/(R*T)));

        // C:P1_NCX

        var mu_P1_NCX: J_per_mol;
        var v_P1_NCX: fmol_per_sec;

        e_72 = mu_P1_NCX;
        v_P1_NCX = f_72;
        ode(x_P1_NCX, t) = v_P1_NCX;
        mu_P1_NCX = R*T*ln(K_1_NCX*x_P1_NCX);

        // C:P2_NCX

        var mu_P2_NCX: J_per_mol;
        var v_P2_NCX: fmol_per_sec;

        e_79 = mu_P2_NCX;
        v_P2_NCX = f_79;
        ode(x_P2_NCX, t) = v_P2_NCX;
        mu_P2_NCX = R*T*ln(K_2_NCX*x_P2_NCX);

        // C:P3_NCX

        var mu_P3_NCX: J_per_mol;
        var v_P3_NCX: fmol_per_sec;

        e_95 = mu_P3_NCX;
        v_P3_NCX = f_95;
        ode(x_P3_NCX, t) = v_P3_NCX;
        mu_P3_NCX = R*T*ln(K_3_NCX*x_P3_NCX);

        // C:P4_NCX

        var mu_P4_NCX: J_per_mol;
        var v_P4_NCX: fmol_per_sec;

        e_91 = mu_P4_NCX;
        v_P4_NCX = f_91;
        ode(x_P4_NCX, t) = v_P4_NCX;
        mu_P4_NCX = R*T*ln(K_4_NCX*x_P4_NCX);

        // C:P5_NCX

        var mu_P5_NCX: J_per_mol;
        var v_P5_NCX: fmol_per_sec;

        e_90 = mu_P5_NCX;
        v_P5_NCX = f_90;
        ode(x_P5_NCX, t) = v_P5_NCX;
        mu_P5_NCX = R*T*ln(K_5_NCX*x_P5_NCX);

        // C:P6_NCX

        var mu_P6_NCX: J_per_mol;
        var v_P6_NCX: fmol_per_sec;

        e_67 = mu_P6_NCX;
        v_P6_NCX = f_67;
        ode(x_P6_NCX, t) = v_P6_NCX;
        mu_P6_NCX = R*T*ln(K_6_NCX*x_P6_NCX);

        // Re:r1

        var Af_r1: J_per_mol;
        var Ar_r1: J_per_mol;
        var v_r1: fmol_per_sec;

        Af_r1 = e_75;
        Ar_r1 = e_76;
        f_75 = v_r1;
        f_76 = v_r1;
        v_r1 = kappa_1_NCX*(exp(Af_r1/(R*T))-exp(Ar_r1/(R*T)));

        // Re:r2

        var Af_r2: J_per_mol;
        var Ar_r2: J_per_mol;
        var v_r2: fmol_per_sec;

        Af_r2 = e_81;
        Ar_r2 = e_93;
        f_81 = v_r2;
        f_93 = v_r2;
        v_r2 = kappa_2_NCX*(exp(Af_r2/(R*T))-exp(Ar_r2/(R*T)));

        // Re:r3

        var Af_r3: J_per_mol;
        var Ar_r3: J_per_mol;
        var v_r3: fmol_per_sec;

        Af_r3 = e_94;
        Ar_r3 = e_92;
        f_94 = v_r3;
        f_92 = v_r3;
        v_r3 = kappa_3_NCX*(exp(Af_r3/(R*T))-exp(Ar_r3/(R*T)));

        // Re:r4

        var Af_r4: J_per_mol;
        var Ar_r4: J_per_mol;
        var v_r4: fmol_per_sec;

        Af_r4 = e_88;
        Ar_r4 = e_87;
        f_88 = v_r4;
        f_87 = v_r4;
        v_r4 = kappa_4_NCX*(exp(Af_r4/(R*T))-exp(Ar_r4/(R*T)));

        // Re:r5

        var Af_r5: J_per_mol;
        var Ar_r5: J_per_mol;
        var v_r5: fmol_per_sec;

        Af_r5 = e_84;
        Ar_r5 = e_83;
        f_84 = v_r5;
        f_83 = v_r5;
        v_r5 = kappa_5_NCX*(exp(Af_r5/(R*T))-exp(Ar_r5/(R*T)));

        // Re:r6

        var Af_r6: J_per_mol;
        var Ar_r6: J_per_mol;
        var v_r6: fmol_per_sec;

        Af_r6 = e_69;
        Ar_r6 = e_70;
        f_69 = v_r6;
        f_70 = v_r6;
        v_r6 = kappa_6_NCX*(exp(Af_r6/(R*T))-exp(Ar_r6/(R*T)));

        // C:H

        var mu_H: J_per_mol;
        var v_H: fmol_per_sec;

        v_H = f_109;
        ode(x_H, t) = 0{fmol_per_sec};
        mu_H = R*T*ln(K_H*x_H);

        // C:MgADP

        var mu_MgADP: J_per_mol;
        var v_MgADP: fmol_per_sec;
        v_MgADP = f_105;
        ode(x_MgADP, t) = 0{fmol_per_sec};
        mu_MgADP = R*T*ln(K_MgADP*x_MgADP);

        // C:MgATP

        var mu_MgATP: J_per_mol;
        var v_MgATP: fmol_per_sec;
        v_MgATP = f_111;
        ode(x_MgATP, t) = 0{fmol_per_sec};
        mu_MgATP = R*T*ln(K_MgATP*x_MgATP);

        // C:P7

        var mu_P7: J_per_mol;
        var v_P7: fmol_per_sec;
        v_P7 = f_104;
        ode(x_P7, t) = v_P7;
        mu_P7 = R*T*ln(K_7_NaK*x_P7);

        // C:P_i

        var mu_P_i: J_per_mol;
        var v_P_i: fmol_per_sec;
        v_P_i = f_110;
        ode(x_P_i, t) = 0{fmol_per_sec};
        mu_P_i = R*T*ln(K_P*x_P_i);

        // Re:R13

        var Af_R13: J_per_mol;
        var Ar_R13: J_per_mol;
        var v_R13: fmol_per_sec;

        Af_R13 = e_107;
        Ar_R13 = e_108;
        f_107 = v_R13;
        f_108 = v_R13;
        v_R13 = kappa_13_NaK*(exp(Af_R13/(R*T))-exp(Ar_R13/(R*T)));

        // Re:R15

        var Af_R15: J_per_mol;
        var Ar_R15: J_per_mol;
        var v_R15: fmol_per_sec;

        Af_R15 = e_112;
        Ar_R15 = e_101;
        f_112 = v_R15;
        f_101 = v_R15;
        v_R15 = kappa_15_NaK*(exp(Af_R15/(R*T))-exp(Ar_R15/(R*T)));

        // Re:R6

        var Af_R6: J_per_mol;
        var Ar_R6: J_per_mol;
        var v_R6: fmol_per_sec;

        Af_R6 = e_100;
        Ar_R6 = e_102;
        f_100 = v_R6;
        f_102 = v_R6;
        v_R6 = kappa_6_NaK*(exp(Af_R6/(R*T))-exp(Ar_R6/(R*T)));

        // Re:R7

        var Af_R7: J_per_mol;
        var Ar_R7: J_per_mol;
        var v_R7: fmol_per_sec;

        Af_R7 = e_106;
        Ar_R7 = e_114;
        f_106 = v_R7;
        f_114 = v_R7;
        v_R7 = kappa_7_NaK*(exp(Af_R7/(R*T))-exp(Ar_R7/(R*T)));

        // Sf:I_stim

        f_121 = I_stim;

        // C:P14

        var mu_P14: J_per_mol;
        var v_P14: fmol_per_sec;
        v_P14 = f_172;
        ode(x_P14, t) = v_P14;
        mu_P14 = R*T*ln(K_14_NaK*x_P14);

        // C:P15

        var mu_P15: J_per_mol;
        var v_P15: fmol_per_sec;
        v_P15 = f_176;
        ode(x_P15, t) = v_P15;
        mu_P15 = R*T*ln(K_15_NaK*x_P15);

        // Re:R14

        var Af_R14: J_per_mol;
        var Ar_R14: J_per_mol;
        var v_R14: fmol_per_sec;

        Af_R14 = e_174;
        Ar_R14 = e_175;
        f_174 = v_R14;
        f_175 = v_R14;
        v_R14 = kappa_14_NaK*(exp(Af_R14/(R*T))-exp(Ar_R14/(R*T)));

        // C:P1

        var mu_P1: J_per_mol;
        var v_P1: fmol_per_sec;
        v_P1 = f_178;
        ode(x_P1, t) = v_P1;
        mu_P1 = R*T*ln(K_1_NaK*x_P1);

        // C:P2

        var mu_P2: J_per_mol;
        var v_P2: fmol_per_sec;
        v_P2 = f_187;
        ode(x_P2, t) = v_P2;
        mu_P2 = R*T*ln(K_2_NaK*x_P2);

        // C:P3

        var mu_P3: J_per_mol;
        var v_P3: fmol_per_sec;
        v_P3 = f_182;
        ode(x_P3, t) = v_P3;
        mu_P3 = R*T*ln(K_3_NaK*x_P3);

        // C:P4

        var mu_P4: J_per_mol;
        var v_P4: fmol_per_sec;
        v_P4 = f_192;
        ode(x_P4, t) = v_P4;
        mu_P4 = R*T*ln(K_4_NaK*x_P4);

        // C:P5

        var mu_P5: J_per_mol;
        var v_P5: fmol_per_sec;
        v_P5 = f_199;
        ode(x_P5, t) = v_P5;
        mu_P5 = R*T*ln(K_5_NaK*x_P5);

        // C:P6

        var mu_P6: J_per_mol;
        var v_P6: fmol_per_sec;
        v_P6 = f_205;
        ode(x_P6, t) = v_P6;
        mu_P6 = R*T*ln(K_6_NaK*x_P6);

        // Re:R1

        var Af_R1: J_per_mol;
        var Ar_R1: J_per_mol;
        var v_R1: fmol_per_sec;

        Af_R1 = e_183;
        Ar_R1 = e_184;
        f_183 = v_R1;
        f_184 = v_R1;
        v_R1 = kappa_1_NaK*(exp(Af_R1/(R*T))-exp(Ar_R1/(R*T)));

        // Re:R2

        var Af_R2: J_per_mol;
        var Ar_R2: J_per_mol;
        var v_R2: fmol_per_sec;

        Af_R2 = e_179;
        Ar_R2 = e_180;
        f_179 = v_R2;
        f_180 = v_R2;
        v_R2 = kappa_2_NaK*(exp(Af_R2/(R*T))-exp(Ar_R2/(R*T)));

        // Re:R3

        var Af_R3: J_per_mol;
        var Ar_R3: J_per_mol;
        var v_R3: fmol_per_sec;

        Af_R3 = e_190;
        Ar_R3 = e_191;
        f_190 = v_R3;
        f_191 = v_R3;
        v_R3 = kappa_3_NaK*(exp(Af_R3/(R*T))-exp(Ar_R3/(R*T)));

        // Re:R4

        var Af_R4: J_per_mol;
        var Ar_R4: J_per_mol;
        var v_R4: fmol_per_sec;

        Af_R4 = e_194;
        Ar_R4 = e_195;
        f_194 = v_R4;
        f_195 = v_R4;
        v_R4 = kappa_4_NaK*(exp(Af_R4/(R*T))-exp(Ar_R4/(R*T)));

        // Re:R5

        var Af_R5: J_per_mol;
        var Ar_R5: J_per_mol;
        var v_R5: fmol_per_sec;

        Af_R5 = e_201;
        Ar_R5 = e_203;
        f_201 = v_R5;
        f_203 = v_R5;
        v_R5 = kappa_5_NaK*(exp(Af_R5/(R*T))-exp(Ar_R5/(R*T)));

        // C:P10

        var mu_P10: J_per_mol;
        var v_P10: fmol_per_sec;
        v_P10 = f_218;
        ode(x_P10, t) = v_P10;
        mu_P10 = R*T*ln(K_10_NaK*x_P10);

        // C:P11

        var mu_P11: J_per_mol;
        var v_P11: fmol_per_sec;
        v_P11 = f_223;
        ode(x_P11, t) = v_P11;
        mu_P11 = R*T*ln(K_11_NaK*x_P11);

        // C:P12

        var mu_P12: J_per_mol;
        var v_P12: fmol_per_sec;
        v_P12 = f_231;
        ode(x_P12, t) = v_P12;
        mu_P12 = R*T*ln(K_12_NaK*x_P12);

        // C:P13

        var mu_P13: J_per_mol;
        var v_P13: fmol_per_sec;
        v_P13 = f_227;
        ode(x_P13, t) = v_P13;
        mu_P13 = R*T*ln(K_13_NaK*x_P13);

        // C:P8

        var mu_P8: J_per_mol;
        var v_P8: fmol_per_sec;
        v_P8 = f_207;
        ode(x_P8, t) = v_P8;
        mu_P8 = R*T*ln(K_8_NaK*x_P8);

        // C:P9

        var mu_P9: J_per_mol;
        var v_P9: fmol_per_sec;
        v_P9 = f_213;
        ode(x_P9, t) = v_P9;
        mu_P9 = R*T*ln(K_9_NaK*x_P9);

        // Re:R10

        var Af_R10: J_per_mol;
        var Ar_R10: J_per_mol;
        var v_R10: fmol_per_sec;

        Af_R10 = e_219;
        Ar_R10 = e_220;
        f_219 = v_R10;
        f_220 = v_R10;
        v_R10 = kappa_10_NaK*(exp(Af_R10/(R*T))-exp(Ar_R10/(R*T)));

        // Re:R11

        var Af_R11: J_per_mol;
        var Ar_R11: J_per_mol;
        var v_R11: fmol_per_sec;

        Af_R11 = e_229;
        Ar_R11 = e_230;
        f_229 = v_R11;
        f_230 = v_R11;
        v_R11 = kappa_11_NaK*(exp(Af_R11/(R*T))-exp(Ar_R11/(R*T)));

        // Re:R12

        var Af_R12: J_per_mol;
        var Ar_R12: J_per_mol;
        var v_R12: fmol_per_sec;

        Af_R12 = e_225;
        Ar_R12 = e_226;
        f_225 = v_R12;
        f_226 = v_R12;
        v_R12 = kappa_12_NaK*(exp(Af_R12/(R*T))-exp(Ar_R12/(R*T)));

        // Re:R8

        var Af_R8: J_per_mol;
        var Ar_R8: J_per_mol;
        var v_R8: fmol_per_sec;

        Af_R8 = e_208;
        Ar_R8 = e_209;
        f_208 = v_R8;
        f_209 = v_R8;
        v_R8 = kappa_8_NaK*(exp(Af_R8/(R*T))-exp(Ar_R8/(R*T)));

        // Re:R9

        var Af_R9: J_per_mol;
        var Ar_R9: J_per_mol;
        var v_R9: fmol_per_sec;

        Af_R9 = e_214;
        Ar_R9 = e_234;
        f_214 = v_R9;
        f_234 = v_R9;
        v_R9 = kappa_9_NaK*(exp(Af_R9/(R*T))-exp(Ar_R9/(R*T)));

        // Re_GHK:Na

        var Af_Na: J_per_mol;
        var Ar_Na: J_per_mol;
        var Am_Na: J_per_mol;
        var v_Na: fmol_per_sec;

        Af_Na = e_237;
        Ar_Na = e_240;
        Am_Na = e_239;
        f_237 = v_Na;
        f_240 = v_Na;
        f_239 = 0{fmol_per_sec};

        v_Na = sel
            case Am_Na == 0{J_per_mol}:
                kappa_Na*(exp(Af_Na/(R*T))-exp(Ar_Na/(R*T)));
            otherwise:
                kappa_Na*Am_Na/(R*T)/(exp(Am_Na/(R*T))-1{dimensionless})*(exp(Af_Na/(R*T))-exp(Ar_Na/(R*T)));
        endsel;

        // Re_GHK:K

        var Af_K: J_per_mol;
        var Ar_K: J_per_mol;
        var Am_K: J_per_mol;
        var v_K: fmol_per_sec;

        Af_K = e_276;
        Ar_K = e_279;
        Am_K = e_278;
        f_276 = v_K;
        f_279 = v_K;
        f_278 = 0{fmol_per_sec};

        v_K = sel
            case Am_K == 0{J_per_mol}:
                kappa_K*(exp(Af_K/(R*T))-exp(Ar_K/(R*T)));
            otherwise:
                kappa_K*Am_K/(R*T)/(exp(Am_K/(R*T))-1{dimensionless})*(exp(Af_K/(R*T))-exp(Ar_K/(R*T)));
        endsel;

        // Re_GHK:Kp

        var Af_Kp: J_per_mol;
        var Ar_Kp: J_per_mol;
        var Am_Kp: J_per_mol;
        var v_Kp: fmol_per_sec;

        Af_Kp = e_297;
        Ar_Kp = e_300;
        Am_Kp = e_299;
        f_297 = v_Kp;
        f_300 = v_Kp;
        f_299 = 0{fmol_per_sec};

        v_Kp = sel
            case Am_Kp == 0{J_per_mol}:
                kappa_Kp*(exp(Af_Kp/(R*T))-exp(Ar_Kp/(R*T)));
            otherwise:
                kappa_Kp*Am_Kp/(R*T)/(exp(Am_Kp/(R*T))-1{dimensionless})*(exp(Af_Kp/(R*T))-exp(Ar_Kp/(R*T)));
        endsel;

        // C:C_Kp

        var mu_C_Kp: J_per_mol;
        var v_C_Kp: fmol_per_sec;
        v_C_Kp = f_304;
        ode(x_C_Kp, t) = v_C_Kp;
        mu_C_Kp = R*T*ln(K_C_Kp*x_C_Kp);

        // C:O_Kp

        var mu_O_Kp: J_per_mol;
        var v_O_Kp: fmol_per_sec;
        v_O_Kp = f_305;
        ode(x_O_Kp, t) = v_O_Kp;
        mu_O_Kp = R*T*ln(K_O_Kp*x_O_Kp);

        // Re_GHK:K1

        var Af_K1: J_per_mol;
        var Ar_K1: J_per_mol;
        var Am_K1: J_per_mol;
        var v_K1: fmol_per_sec;

        Af_K1 = e_310;
        Ar_K1 = e_313;
        Am_K1 = e_312;
        f_310 = v_K1;
        f_313 = v_K1;
        f_312 = 0{fmol_per_sec};

        v_K1 = sel
            case Am_K1 == 0{J_per_mol}:
                kappa_K1*(exp(Af_K1/(R*T))-exp(Ar_K1/(R*T)));
            otherwise:
                kappa_K1*Am_K1/(R*T)/(exp(Am_K1/(R*T))-1{dimensionless})*(exp(Af_K1/(R*T))-exp(Ar_K1/(R*T)));
        endsel;

        // C:C_K1

        var mu_C_K1: J_per_mol;
        var v_C_K1: fmol_per_sec;
        v_C_K1 = f_320;
        ode(x_C_K1, t) = v_C_K1;
        mu_C_K1 = R*T*ln(K_C_K1*x_C_K1);

        // C:O_K1

        var mu_O_K1: J_per_mol;
        var v_O_K1: fmol_per_sec;
        v_O_K1 = f_321;
        ode(x_O_K1, t) = v_O_K1;
        mu_O_K1 = R*T*ln(K_O_K1*x_O_K1);

        // Re_GHK:LCC_Ca1

        var Af_LCC_Ca1: J_per_mol;
        var Ar_LCC_Ca1: J_per_mol;
        var Am_LCC_Ca1: J_per_mol;
        var v_LCC_Ca1: fmol_per_sec;

        Af_LCC_Ca1 = e_326;
        Ar_LCC_Ca1 = e_329;
        Am_LCC_Ca1 = e_328;
        f_326 = v_LCC_Ca1;
        f_329 = v_LCC_Ca1;
        f_328 = 0{fmol_per_sec};

        v_LCC_Ca1 = sel
            case Am_LCC_Ca1 == 0{J_per_mol}:
                kappa_LCC_Ca1*(exp(Af_LCC_Ca1/(R*T))-exp(Ar_LCC_Ca1/(R*T)));
            otherwise:
                kappa_LCC_Ca1*Am_LCC_Ca1/(R*T)/(exp(Am_LCC_Ca1/(R*T))-1{dimensionless})*(exp(Af_LCC_Ca1/(R*T))-exp(Ar_LCC_Ca1/(R*T)));
        endsel;

        // Re_GHK:LCC_Ca2

        var Af_LCC_Ca2: J_per_mol;
        var Ar_LCC_Ca2: J_per_mol;
        var Am_LCC_Ca2: J_per_mol;
        var v_LCC_Ca2: fmol_per_sec;

        Af_LCC_Ca2 = e_335;
        Ar_LCC_Ca2 = e_338;
        Am_LCC_Ca2 = e_337;
        f_335 = v_LCC_Ca2;
        f_338 = v_LCC_Ca2;
        f_337 = 0{fmol_per_sec};

        v_LCC_Ca2 = sel
            case Am_LCC_Ca2 == 0{J_per_mol}:
                kappa_LCC_Ca2*(exp(Af_LCC_Ca2/(R*T))-exp(Ar_LCC_Ca2/(R*T)));
            otherwise:
                kappa_LCC_Ca2*Am_LCC_Ca2/(R*T)/(exp(Am_LCC_Ca2/(R*T))-1{dimensionless})*(exp(Af_LCC_Ca2/(R*T))-exp(Ar_LCC_Ca2/(R*T)));
        endsel;

        // Re_GHK:LCC_K1

        var Af_LCC_K1: J_per_mol;
        var Ar_LCC_K1: J_per_mol;
        var Am_LCC_K1: J_per_mol;
        var v_LCC_K1: fmol_per_sec;

        Af_LCC_K1 = e_344;
        Ar_LCC_K1 = e_347;
        Am_LCC_K1 = e_346;
        f_344 = v_LCC_K1;
        f_347 = v_LCC_K1;
        f_346 = 0{fmol_per_sec};

        v_LCC_K1 = sel
            case Am_LCC_K1 == 0{J_per_mol}:
                kappa_LCC_K1*(exp(Af_LCC_K1/(R*T))-exp(Ar_LCC_K1/(R*T)));
            otherwise:
                kappa_LCC_K1*Am_LCC_K1/(R*T)/(exp(Am_LCC_K1/(R*T))-1{dimensionless})*(exp(Af_LCC_K1/(R*T))-exp(Ar_LCC_K1/(R*T)));
        endsel;

        // Re_GHK:LCC_K2

        var Af_LCC_K2: J_per_mol;
        var Ar_LCC_K2: J_per_mol;
        var Am_LCC_K2: J_per_mol;
        var v_LCC_K2: fmol_per_sec;

        Af_LCC_K2 = e_353;
        Ar_LCC_K2 = e_356;
        Am_LCC_K2 = e_355;
        f_353 = v_LCC_K2;
        f_356 = v_LCC_K2;
        f_355 = 0{fmol_per_sec};

        v_LCC_K2 = sel
            case Am_LCC_K2 == 0{J_per_mol}:
                kappa_LCC_K2*(exp(Af_LCC_K2/(R*T))-exp(Ar_LCC_K2/(R*T)));
            otherwise:
                kappa_LCC_K2*Am_LCC_K2/(R*T)/(exp(Am_LCC_K2/(R*T))-1{dimensionless})*(exp(Af_LCC_K2/(R*T))-exp(Ar_LCC_K2/(R*T)));
        endsel;

        // C:S00_K

        var mu_S00_K: J_per_mol;
        var v_S00_K: fmol_per_sec;
        v_S00_K = f_402;
        ode(x_S00_K, t) = v_S00_K;
        mu_S00_K = R*T*ln(K_00_K*x_S00_K);

        // C:S01_K

        var mu_S01_K: J_per_mol;
        var v_S01_K: fmol_per_sec;
        v_S01_K = f_409;
        ode(x_S01_K, t) = v_S01_K;
        mu_S01_K = R*T*ln(K_01_K*x_S01_K);

        // C:S10_K

        var mu_S10_K: J_per_mol;
        var v_S10_K: fmol_per_sec;
        v_S10_K = f_403;
        ode(x_S10_K, t) = v_S10_K;
        mu_S10_K = R*T*ln(K_10_K*x_S10_K);

        // C:S11_K

        var mu_S11_K: J_per_mol;
        var v_S11_K: fmol_per_sec;
        v_S11_K = f_412;
        ode(x_S11_K, t) = v_S11_K;
        mu_S11_K = R*T*ln(K_11_K*x_S11_K);

        // C:S20_K

        var mu_S20_K: J_per_mol;
        var v_S20_K: fmol_per_sec;
        v_S20_K = f_408;
        ode(x_S20_K, t) = v_S20_K;
        mu_S20_K = R*T*ln(K_20_K*x_S20_K);

        // C:S21_K

        var mu_S21_K: J_per_mol;
        var v_S21_K: fmol_per_sec;
        v_S21_K = f_415;
        ode(x_S21_K, t) = v_S21_K;
        mu_S21_K = R*T*ln(K_21_K*x_S21_K);

        // Re:gKp

        var Af_gKp: J_per_mol;
        var Ar_gKp: J_per_mol;
        var v_gKp: fmol_per_sec;

        Af_gKp = e_422;
        Ar_gKp = e_423;
        f_422 = v_gKp;
        f_423 = v_gKp;
        v_gKp = kappa_gKp*(exp(Af_gKp/(R*T))-exp(Ar_gKp/(R*T)));

        // Re:gK1

        var Af_gK1: J_per_mol;
        var Ar_gK1: J_per_mol;
        var v_gK1: fmol_per_sec;

        Af_gK1 = e_428;
        Ar_gK1 = e_429;
        f_428 = v_gK1;
        f_429 = v_gK1;
        v_gK1 = kappa_gK1*(exp(Af_gK1/(R*T))-exp(Ar_gK1/(R*T)));

        // Re:fCa00

        var Af_fCa00: J_per_mol;
        var Ar_fCa00: J_per_mol;
        var v_fCa00: fmol_per_sec;

        Af_fCa00 = e_434;
        Ar_fCa00 = e_435;
        f_434 = v_fCa00;
        f_435 = v_fCa00;
        v_fCa00 = kappa_fCa00*(exp(Af_fCa00/(R*T))-exp(Ar_fCa00/(R*T)));

        // Re:fCa01

        var Af_fCa01: J_per_mol;
        var Ar_fCa01: J_per_mol;
        var v_fCa01: fmol_per_sec;

        Af_fCa01 = e_438;
        Ar_fCa01 = e_439;
        f_438 = v_fCa01;
        f_439 = v_fCa01;
        v_fCa01 = kappa_fCa01*(exp(Af_fCa01/(R*T))-exp(Ar_fCa01/(R*T)));

        // Re:fCa02

        var Af_fCa02: J_per_mol;
        var Ar_fCa02: J_per_mol;
        var v_fCa02: fmol_per_sec;

        Af_fCa02 = e_442;
        Ar_fCa02 = e_443;
        f_442 = v_fCa02;
        f_443 = v_fCa02;
        v_fCa02 = kappa_fCa02*(exp(Af_fCa02/(R*T))-exp(Ar_fCa02/(R*T)));

        // Re:fCa10

        var Af_fCa10: J_per_mol;
        var Ar_fCa10: J_per_mol;
        var v_fCa10: fmol_per_sec;

        Af_fCa10 = e_446;
        Ar_fCa10 = e_447;
        f_446 = v_fCa10;
        f_447 = v_fCa10;
        v_fCa10 = kappa_fCa10*(exp(Af_fCa10/(R*T))-exp(Ar_fCa10/(R*T)));

        // Re:fCa11

        var Af_fCa11: J_per_mol;
        var Ar_fCa11: J_per_mol;
        var v_fCa11: fmol_per_sec;

        Af_fCa11 = e_450;
        Ar_fCa11 = e_451;
        f_450 = v_fCa11;
        f_451 = v_fCa11;
        v_fCa11 = kappa_fCa11*(exp(Af_fCa11/(R*T))-exp(Ar_fCa11/(R*T)));

        // Re:fCa12

        var Af_fCa12: J_per_mol;
        var Ar_fCa12: J_per_mol;
        var v_fCa12: fmol_per_sec;

        Af_fCa12 = e_454;
        Ar_fCa12 = e_455;
        f_454 = v_fCa12;
        f_455 = v_fCa12;
        v_fCa12 = kappa_fCa12*(exp(Af_fCa12/(R*T))-exp(Ar_fCa12/(R*T)));

        // C:S000_LCC

        var mu_S000_LCC: J_per_mol;
        var v_S000_LCC: fmol_per_sec;

0_LCC;
        v_S000_LCC = f_475;
        ode(x_S000_LCC, t) = v_S000_LCC;
        mu_S000_LCC = R*T*ln(K_000_LCC*x_S000_LCC);

        // C:S010_LCC

        var mu_S010_LCC: J_per_mol;
        var v_S010_LCC: fmol_per_sec;

0_LCC;
        v_S010_LCC = f_472;
        ode(x_S010_LCC, t) = v_S010_LCC;
        mu_S010_LCC = R*T*ln(K_010_LCC*x_S010_LCC);

        // C:S020_LCC

        var mu_S020_LCC: J_per_mol;
        var v_S020_LCC: fmol_per_sec;

0_LCC;
        v_S020_LCC = f_458;
        ode(x_S020_LCC, t) = v_S020_LCC;
        mu_S020_LCC = R*T*ln(K_020_LCC*x_S020_LCC);

        // C:S100_LCC

        var mu_S100_LCC: J_per_mol;
        var v_S100_LCC: fmol_per_sec;

0_LCC;
        v_S100_LCC = f_473;
        ode(x_S100_LCC, t) = v_S100_LCC;
        mu_S100_LCC = R*T*ln(K_100_LCC*x_S100_LCC);

        // C:S110_LCC

        var mu_S110_LCC: J_per_mol;
        var v_S110_LCC: fmol_per_sec;

0_LCC;
        v_S110_LCC = f_474;
        ode(x_S110_LCC, t) = v_S110_LCC;
        mu_S110_LCC = R*T*ln(K_110_LCC*x_S110_LCC);

        // C:S120_LCC

        var mu_S120_LCC: J_per_mol;
        var v_S120_LCC: fmol_per_sec;

0_LCC;
        v_S120_LCC = f_460;
        ode(x_S120_LCC, t) = v_S120_LCC;
        mu_S120_LCC = R*T*ln(K_120_LCC*x_S120_LCC);

        // C:S001_LCC

        var mu_S001_LCC: J_per_mol;
        var v_S001_LCC: fmol_per_sec;

1_LCC;
        v_S001_LCC = f_499;
        ode(x_S001_LCC, t) = v_S001_LCC;
        mu_S001_LCC = R*T*ln(K_001_LCC*x_S001_LCC);

        // C:S011_LCC

        var mu_S011_LCC: J_per_mol;
        var v_S011_LCC: fmol_per_sec;

1_LCC;
        v_S011_LCC = f_496;
        ode(x_S011_LCC, t) = v_S011_LCC;
        mu_S011_LCC = R*T*ln(K_011_LCC*x_S011_LCC);

        // C:S021_LCC

        var mu_S021_LCC: J_per_mol;
        var v_S021_LCC: fmol_per_sec;

1_LCC;
        v_S021_LCC = f_482;
        ode(x_S021_LCC, t) = v_S021_LCC;
        mu_S021_LCC = R*T*ln(K_021_LCC*x_S021_LCC);

        // C:S101_LCC

        var mu_S101_LCC: J_per_mol;
        var v_S101_LCC: fmol_per_sec;

1_LCC;
        v_S101_LCC = f_497;
        ode(x_S101_LCC, t) = v_S101_LCC;
        mu_S101_LCC = R*T*ln(K_101_LCC*x_S101_LCC);

        // C:S111_LCC

        var mu_S111_LCC: J_per_mol;
        var v_S111_LCC: fmol_per_sec;

1_LCC;
        v_S111_LCC = f_498;
        ode(x_S111_LCC, t) = v_S111_LCC;
        mu_S111_LCC = R*T*ln(K_111_LCC*x_S111_LCC);

        // C:S121_LCC

        var mu_S121_LCC: J_per_mol;
        var v_S121_LCC: fmol_per_sec;

1_LCC;
        v_S121_LCC = f_484;
        ode(x_S121_LCC, t) = v_S121_LCC;
        mu_S121_LCC = R*T*ln(K_121_LCC*x_S121_LCC);

        // C:S000_Na

        var mu_S000_Na: J_per_mol;
        var v_S000_Na: fmol_per_sec;

0_Na;
        v_S000_Na = f_506;
        ode(x_S000_Na, t) = v_S000_Na;
        mu_S000_Na = R*T*ln(K_000_Na*x_S000_Na);

        // C:S010_Na

        var mu_S010_Na: J_per_mol;
        var v_S010_Na: fmol_per_sec;

0_Na;
        v_S010_Na = f_516;
        ode(x_S010_Na, t) = v_S010_Na;
        mu_S010_Na = R*T*ln(K_010_Na*x_S010_Na);

        // C:S100_Na

        var mu_S100_Na: J_per_mol;
        var v_S100_Na: fmol_per_sec;

0_Na;
        v_S100_Na = f_507;
        ode(x_S100_Na, t) = v_S100_Na;
        mu_S100_Na = R*T*ln(K_100_Na*x_S100_Na);

        // C:S110_Na

        var mu_S110_Na: J_per_mol;
        var v_S110_Na: fmol_per_sec;

0_Na;
        v_S110_Na = f_519;
        ode(x_S110_Na, t) = v_S110_Na;
        mu_S110_Na = R*T*ln(K_110_Na*x_S110_Na);

        // C:S200_Na

        var mu_S200_Na: J_per_mol;
        var v_S200_Na: fmol_per_sec;

0_Na;
        v_S200_Na = f_512;
        ode(x_S200_Na, t) = v_S200_Na;
        mu_S200_Na = R*T*ln(K_200_Na*x_S200_Na);

        // C:S210_Na

        var mu_S210_Na: J_per_mol;
        var v_S210_Na: fmol_per_sec;

0_Na;
        v_S210_Na = f_522;
        ode(x_S210_Na, t) = v_S210_Na;
        mu_S210_Na = R*T*ln(K_210_Na*x_S210_Na);

        // C:S300_Na

        var mu_S300_Na: J_per_mol;
        var v_S300_Na: fmol_per_sec;

0_Na;
        v_S300_Na = f_515;
        ode(x_S300_Na, t) = v_S300_Na;
        mu_S300_Na = R*T*ln(K_300_Na*x_S300_Na);

        // C:S310_Na

        var mu_S310_Na: J_per_mol;
        var v_S310_Na: fmol_per_sec;

0_Na;
        v_S310_Na = f_524;
        ode(x_S310_Na, t) = v_S310_Na;
        mu_S310_Na = R*T*ln(K_310_Na*x_S310_Na);

        // C:S001_Na

        var mu_S001_Na: J_per_mol;
        var v_S001_Na: fmol_per_sec;

1_Na;
        v_S001_Na = f_534;
        ode(x_S001_Na, t) = v_S001_Na;
        mu_S001_Na = R*T*ln(K_001_Na*x_S001_Na);

        // C:S011_Na

        var mu_S011_Na: J_per_mol;
        var v_S011_Na: fmol_per_sec;

1_Na;
        v_S011_Na = f_544;
        ode(x_S011_Na, t) = v_S011_Na;
        mu_S011_Na = R*T*ln(K_011_Na*x_S011_Na);

        // C:S101_Na

        var mu_S101_Na: J_per_mol;
        var v_S101_Na: fmol_per_sec;

1_Na;
        v_S101_Na = f_535;
        ode(x_S101_Na, t) = v_S101_Na;
        mu_S101_Na = R*T*ln(K_101_Na*x_S101_Na);

        // C:S111_Na

        var mu_S111_Na: J_per_mol;
        var v_S111_Na: fmol_per_sec;

1_Na;
        v_S111_Na = f_547;
        ode(x_S111_Na, t) = v_S111_Na;
        mu_S111_Na = R*T*ln(K_111_Na*x_S111_Na);

        // C:S201_Na

        var mu_S201_Na: J_per_mol;
        var v_S201_Na: fmol_per_sec;

1_Na;
        v_S201_Na = f_540;
        ode(x_S201_Na, t) = v_S201_Na;
        mu_S201_Na = R*T*ln(K_201_Na*x_S201_Na);

        // C:S211_Na

        var mu_S211_Na: J_per_mol;
        var v_S211_Na: fmol_per_sec;

1_Na;
        v_S211_Na = f_550;
        ode(x_S211_Na, t) = v_S211_Na;
        mu_S211_Na = R*T*ln(K_211_Na*x_S211_Na);

        // C:S301_Na

        var mu_S301_Na: J_per_mol;
        var v_S301_Na: fmol_per_sec;

1_Na;
        v_S301_Na = f_543;
        ode(x_S301_Na, t) = v_S301_Na;
        mu_S301_Na = R*T*ln(K_301_Na*x_S301_Na);

        // C:S311_Na

        var mu_S311_Na: J_per_mol;
        var v_S311_Na: fmol_per_sec;

1_Na;
        v_S311_Na = f_552;
        ode(x_S311_Na, t) = v_S311_Na;
        mu_S311_Na = R*T*ln(K_311_Na*x_S311_Na);

        // Re:j000

        var Af_j000: J_per_mol;
        var Ar_j000: J_per_mol;
        var v_j000: fmol_per_sec;

        Af_j000 = e_562;
        Ar_j000 = e_563;
        f_562 = v_j000;
        f_563 = v_j000;
        v_j000 = kappa_j000*(exp(Af_j000/(R*T))-exp(Ar_j000/(R*T)));

        // Re:j010

        var Af_j010: J_per_mol;
        var Ar_j010: J_per_mol;
        var v_j010: fmol_per_sec;

        Af_j010 = e_568;
        Ar_j010 = e_569;
        f_568 = v_j010;
        f_569 = v_j010;
        v_j010 = kappa_j010*(exp(Af_j010/(R*T))-exp(Ar_j010/(R*T)));

        // Re:j100

        var Af_j100: J_per_mol;
        var Ar_j100: J_per_mol;
        var v_j100: fmol_per_sec;

        Af_j100 = e_574;
        Ar_j100 = e_575;
        f_574 = v_j100;
        f_575 = v_j100;
        v_j100 = kappa_j100*(exp(Af_j100/(R*T))-exp(Ar_j100/(R*T)));

        // Re:j110

        var Af_j110: J_per_mol;
        var Ar_j110: J_per_mol;
        var v_j110: fmol_per_sec;

        Af_j110 = e_580;
        Ar_j110 = e_581;
        f_580 = v_j110;
        f_581 = v_j110;
        v_j110 = kappa_j110*(exp(Af_j110/(R*T))-exp(Ar_j110/(R*T)));

        // Re:j200

        var Af_j200: J_per_mol;
        var Ar_j200: J_per_mol;
        var v_j200: fmol_per_sec;

        Af_j200 = e_586;
        Ar_j200 = e_587;
        f_586 = v_j200;
        f_587 = v_j200;
        v_j200 = kappa_j200*(exp(Af_j200/(R*T))-exp(Ar_j200/(R*T)));

        // Re:j210

        var Af_j210: J_per_mol;
        var Ar_j210: J_per_mol;
        var v_j210: fmol_per_sec;

        Af_j210 = e_592;
        Ar_j210 = e_593;
        f_592 = v_j210;
        f_593 = v_j210;
        v_j210 = kappa_j210*(exp(Af_j210/(R*T))-exp(Ar_j210/(R*T)));

        // Re:j300

        var Af_j300: J_per_mol;
        var Ar_j300: J_per_mol;
        var v_j300: fmol_per_sec;

        Af_j300 = e_598;
        Ar_j300 = e_599;
        f_598 = v_j300;
        f_599 = v_j300;
        v_j300 = kappa_j300*(exp(Af_j300/(R*T))-exp(Ar_j300/(R*T)));

        // Re:j310

        var Af_j310: J_per_mol;
        var Ar_j310: J_per_mol;
        var v_j310: fmol_per_sec;

        Af_j310 = e_604;
        Ar_j310 = e_605;
        f_604 = v_j310;
        f_605 = v_j310;
        v_j310 = kappa_j310*(exp(Af_j310/(R*T))-exp(Ar_j310/(R*T)));

        // Re:X00

        var Af_X00: J_per_mol;
        var Ar_X00: J_per_mol;
        var v_X00: fmol_per_sec;

        Af_X00 = e_610;
        Ar_X00 = e_611;
        f_610 = v_X00;
        f_611 = v_X00;
        v_X00 = kappa_X00*(exp(Af_X00/(R*T))-exp(Ar_X00/(R*T)));

        // Re:X01

        var Af_X01: J_per_mol;
        var Ar_X01: J_per_mol;
        var v_X01: fmol_per_sec;

        Af_X01 = e_616;
        Ar_X01 = e_617;
        f_616 = v_X01;
        f_617 = v_X01;
        v_X01 = kappa_X01*(exp(Af_X01/(R*T))-exp(Ar_X01/(R*T)));

        // Re:X10

        var Af_X10: J_per_mol;
        var Ar_X10: J_per_mol;
        var v_X10: fmol_per_sec;

        Af_X10 = e_622;
        Ar_X10 = e_623;
        f_622 = v_X10;
        f_623 = v_X10;
        v_X10 = kappa_X10*(exp(Af_X10/(R*T))-exp(Ar_X10/(R*T)));

        // Re:X11

        var Af_X11: J_per_mol;
        var Ar_X11: J_per_mol;
        var v_X11: fmol_per_sec;

        Af_X11 = e_628;
        Ar_X11 = e_629;
        f_628 = v_X11;
        f_629 = v_X11;
        v_X11 = kappa_X11*(exp(Af_X11/(R*T))-exp(Ar_X11/(R*T)));

        // Re:i00

        var Af_i00: J_per_mol;
        var Ar_i00: J_per_mol;
        var v_i00: fmol_per_sec;

        Af_i00 = e_634;
        Ar_i00 = e_635;
        f_634 = v_i00;
        f_635 = v_i00;
        v_i00 = kappa_i00*(exp(Af_i00/(R*T))-exp(Ar_i00/(R*T)));

        // Re:i10

        var Af_i10: J_per_mol;
        var Ar_i10: J_per_mol;
        var v_i10: fmol_per_sec;

        Af_i10 = e_640;
        Ar_i10 = e_641;
        f_640 = v_i10;
        f_641 = v_i10;
        v_i10 = kappa_i10*(exp(Af_i10/(R*T))-exp(Ar_i10/(R*T)));

        // Re:i20

        var Af_i20: J_per_mol;
        var Ar_i20: J_per_mol;
        var v_i20: fmol_per_sec;

        Af_i20 = e_646;
        Ar_i20 = e_647;
        f_646 = v_i20;
        f_647 = v_i20;
        v_i20 = kappa_i20*(exp(Af_i20/(R*T))-exp(Ar_i20/(R*T)));

        // Re:d000

        var Af_d000: J_per_mol;
        var Ar_d000: J_per_mol;
        var v_d000: fmol_per_sec;

        Af_d000 = e_652;
        Ar_d000 = e_653;
        f_652 = v_d000;
        f_653 = v_d000;
        v_d000 = kappa_d000*(exp(Af_d000/(R*T))-exp(Ar_d000/(R*T)));

        // Re:d010

        var Af_d010: J_per_mol;
        var Ar_d010: J_per_mol;
        var v_d010: fmol_per_sec;

        Af_d010 = e_658;
        Ar_d010 = e_659;
        f_658 = v_d010;
        f_659 = v_d010;
        v_d010 = kappa_d010*(exp(Af_d010/(R*T))-exp(Ar_d010/(R*T)));

        // Re:d020

        var Af_d020: J_per_mol;
        var Ar_d020: J_per_mol;
        var v_d020: fmol_per_sec;

        Af_d020 = e_664;
        Ar_d020 = e_665;
        f_664 = v_d020;
        f_665 = v_d020;
        v_d020 = kappa_d020*(exp(Af_d020/(R*T))-exp(Ar_d020/(R*T)));

        // Re:f1_000

        var Af_f1_000: J_per_mol;
        var Ar_f1_000: J_per_mol;
        var v_f1_000: fmol_per_sec;

        Af_f1_000 = e_670;
        Ar_f1_000 = e_671;
        f_670 = v_f1_000;
        f_671 = v_f1_000;
        v_f1_000 = kappa_f1_000*(exp(Af_f1_000/(R*T))-exp(Ar_f1_000/(R*T)));

        // Re:f2_000

        var Af_f2_000: J_per_mol;
        var Ar_f2_000: J_per_mol;
        var v_f2_000: fmol_per_sec;

        Af_f2_000 = e_676;
        Ar_f2_000 = e_677;
        f_676 = v_f2_000;
        f_677 = v_f2_000;
        v_f2_000 = kappa_f2_000*(exp(Af_f2_000/(R*T))-exp(Ar_f2_000/(R*T)));

        // Re:f3_010

        var Af_f3_010: J_per_mol;
        var Ar_f3_010: J_per_mol;
        var v_f3_010: fmol_per_sec;

        Af_f3_010 = e_682;
        Ar_f3_010 = e_683;
        f_682 = v_f3_010;
        f_683 = v_f3_010;
        v_f3_010 = kappa_f3_010*(exp(Af_f3_010/(R*T))-exp(Ar_f3_010/(R*T)));

        // Re:f1_100

        var Af_f1_100: J_per_mol;
        var Ar_f1_100: J_per_mol;
        var v_f1_100: fmol_per_sec;

        Af_f1_100 = e_688;
        Ar_f1_100 = e_689;
        f_688 = v_f1_100;
        f_689 = v_f1_100;
        v_f1_100 = kappa_f1_100*(exp(Af_f1_100/(R*T))-exp(Ar_f1_100/(R*T)));

        // Re:f2_100

        var Af_f2_100: J_per_mol;
        var Ar_f2_100: J_per_mol;
        var v_f2_100: fmol_per_sec;

        Af_f2_100 = e_694;
        Ar_f2_100 = e_695;
        f_694 = v_f2_100;
        f_695 = v_f2_100;
        v_f2_100 = kappa_f2_100*(exp(Af_f2_100/(R*T))-exp(Ar_f2_100/(R*T)));

        // Re:f3_110

        var Af_f3_110: J_per_mol;
        var Ar_f3_110: J_per_mol;
        var v_f3_110: fmol_per_sec;

        Af_f3_110 = e_700;
        Ar_f3_110 = e_701;
        f_700 = v_f3_110;
        f_701 = v_f3_110;
        v_f3_110 = kappa_f3_110*(exp(Af_f3_110/(R*T))-exp(Ar_f3_110/(R*T)));

        // Re:d001

        var Af_d001: J_per_mol;
        var Ar_d001: J_per_mol;
        var v_d001: fmol_per_sec;

        Af_d001 = e_706;
        Ar_d001 = e_707;
        f_706 = v_d001;
        f_707 = v_d001;
        v_d001 = kappa_d001*(exp(Af_d001/(R*T))-exp(Ar_d001/(R*T)));

        // Re:d011

        var Af_d011: J_per_mol;
        var Ar_d011: J_per_mol;
        var v_d011: fmol_per_sec;

        Af_d011 = e_712;
        Ar_d011 = e_713;
        f_712 = v_d011;
        f_713 = v_d011;
        v_d011 = kappa_d011*(exp(Af_d011/(R*T))-exp(Ar_d011/(R*T)));

        // Re:d021

        var Af_d021: J_per_mol;
        var Ar_d021: J_per_mol;
        var v_d021: fmol_per_sec;

        Af_d021 = e_718;
        Ar_d021 = e_719;
        f_718 = v_d021;
        f_719 = v_d021;
        v_d021 = kappa_d021*(exp(Af_d021/(R*T))-exp(Ar_d021/(R*T)));

        // Re:f1_001

        var Af_f1_001: J_per_mol;
        var Ar_f1_001: J_per_mol;
        var v_f1_001: fmol_per_sec;

        Af_f1_001 = e_724;
        Ar_f1_001 = e_725;
        f_724 = v_f1_001;
        f_725 = v_f1_001;
        v_f1_001 = kappa_f1_001*(exp(Af_f1_001/(R*T))-exp(Ar_f1_001/(R*T)));

        // Re:f2_001

        var Af_f2_001: J_per_mol;
        var Ar_f2_001: J_per_mol;
        var v_f2_001: fmol_per_sec;

        Af_f2_001 = e_730;
        Ar_f2_001 = e_731;
        f_730 = v_f2_001;
        f_731 = v_f2_001;
        v_f2_001 = kappa_f2_001*(exp(Af_f2_001/(R*T))-exp(Ar_f2_001/(R*T)));

        // Re:f3_011

        var Af_f3_011: J_per_mol;
        var Ar_f3_011: J_per_mol;
        var v_f3_011: fmol_per_sec;

        Af_f3_011 = e_736;
        Ar_f3_011 = e_737;
        f_736 = v_f3_011;
        f_737 = v_f3_011;
        v_f3_011 = kappa_f3_011*(exp(Af_f3_011/(R*T))-exp(Ar_f3_011/(R*T)));

        // Re:f1_101

        var Af_f1_101: J_per_mol;
        var Ar_f1_101: J_per_mol;
        var v_f1_101: fmol_per_sec;

        Af_f1_101 = e_742;
        Ar_f1_101 = e_743;
        f_742 = v_f1_101;
        f_743 = v_f1_101;
        v_f1_101 = kappa_f1_101*(exp(Af_f1_101/(R*T))-exp(Ar_f1_101/(R*T)));

        // Re:f2_101

        var Af_f2_101: J_per_mol;
        var Ar_f2_101: J_per_mol;
        var v_f2_101: fmol_per_sec;

        Af_f2_101 = e_748;
        Ar_f2_101 = e_749;
        f_748 = v_f2_101;
        f_749 = v_f2_101;
        v_f2_101 = kappa_f2_101*(exp(Af_f2_101/(R*T))-exp(Ar_f2_101/(R*T)));

        // Re:f3_111

        var Af_f3_111: J_per_mol;
        var Ar_f3_111: J_per_mol;
        var v_f3_111: fmol_per_sec;

        Af_f3_111 = e_754;
        Ar_f3_111 = e_755;
        f_754 = v_f3_111;
        f_755 = v_f3_111;
        v_f3_111 = kappa_f3_111*(exp(Af_f3_111/(R*T))-exp(Ar_f3_111/(R*T)));

        // Re:h000

        var Af_h000: J_per_mol;
        var Ar_h000: J_per_mol;
        var v_h000: fmol_per_sec;

        Af_h000 = e_760;
        Ar_h000 = e_761;
        f_760 = v_h000;
        f_761 = v_h000;
        v_h000 = kappa_h000*(exp(Af_h000/(R*T))-exp(Ar_h000/(R*T)));

        // Re:h100

        var Af_h100: J_per_mol;
        var Ar_h100: J_per_mol;
        var v_h100: fmol_per_sec;

        Af_h100 = e_766;
        Ar_h100 = e_767;
        f_766 = v_h100;
        f_767 = v_h100;
        v_h100 = kappa_h100*(exp(Af_h100/(R*T))-exp(Ar_h100/(R*T)));

        // Re:h200

        var Af_h200: J_per_mol;
        var Ar_h200: J_per_mol;
        var v_h200: fmol_per_sec;

        Af_h200 = e_772;
        Ar_h200 = e_773;
        f_772 = v_h200;
        f_773 = v_h200;
        v_h200 = kappa_h200*(exp(Af_h200/(R*T))-exp(Ar_h200/(R*T)));

        // Re:h300

        var Af_h300: J_per_mol;
        var Ar_h300: J_per_mol;
        var v_h300: fmol_per_sec;

        Af_h300 = e_778;
        Ar_h300 = e_779;
        f_778 = v_h300;
        f_779 = v_h300;
        v_h300 = kappa_h300*(exp(Af_h300/(R*T))-exp(Ar_h300/(R*T)));

        // Re:m000

        var Af_m000: J_per_mol;
        var Ar_m000: J_per_mol;
        var v_m000: fmol_per_sec;

        Af_m000 = e_784;
        Ar_m000 = e_785;
        f_784 = v_m000;
        f_785 = v_m000;
        v_m000 = kappa_m000*(exp(Af_m000/(R*T))-exp(Ar_m000/(R*T)));

        // Re:m100

        var Af_m100: J_per_mol;
        var Ar_m100: J_per_mol;
        var v_m100: fmol_per_sec;

        Af_m100 = e_790;
        Ar_m100 = e_791;
        f_790 = v_m100;
        f_791 = v_m100;
        v_m100 = kappa_m100*(exp(Af_m100/(R*T))-exp(Ar_m100/(R*T)));

        // Re:m200

        var Af_m200: J_per_mol;
        var Ar_m200: J_per_mol;
        var v_m200: fmol_per_sec;

        Af_m200 = e_796;
        Ar_m200 = e_797;
        f_796 = v_m200;
        f_797 = v_m200;
        v_m200 = kappa_m200*(exp(Af_m200/(R*T))-exp(Ar_m200/(R*T)));

        // Re:m010

        var Af_m010: J_per_mol;
        var Ar_m010: J_per_mol;
        var v_m010: fmol_per_sec;

        Af_m010 = e_802;
        Ar_m010 = e_803;
        f_802 = v_m010;
        f_803 = v_m010;
        v_m010 = kappa_m010*(exp(Af_m010/(R*T))-exp(Ar_m010/(R*T)));

        // Re:m110

        var Af_m110: J_per_mol;
        var Ar_m110: J_per_mol;
        var v_m110: fmol_per_sec;

        Af_m110 = e_808;
        Ar_m110 = e_809;
        f_808 = v_m110;
        f_809 = v_m110;
        v_m110 = kappa_m110*(exp(Af_m110/(R*T))-exp(Ar_m110/(R*T)));

        // Re:m210

        var Af_m210: J_per_mol;
        var Ar_m210: J_per_mol;
        var v_m210: fmol_per_sec;

        Af_m210 = e_814;
        Ar_m210 = e_815;
        f_814 = v_m210;
        f_815 = v_m210;
        v_m210 = kappa_m210*(exp(Af_m210/(R*T))-exp(Ar_m210/(R*T)));

        // Re:h001

        var Af_h001: J_per_mol;
        var Ar_h001: J_per_mol;
        var v_h001: fmol_per_sec;

        Af_h001 = e_820;
        Ar_h001 = e_821;
        f_820 = v_h001;
        f_821 = v_h001;
        v_h001 = kappa_h001*(exp(Af_h001/(R*T))-exp(Ar_h001/(R*T)));

        // Re:h101

        var Af_h101: J_per_mol;
        var Ar_h101: J_per_mol;
        var v_h101: fmol_per_sec;

        Af_h101 = e_826;
        Ar_h101 = e_827;
        f_826 = v_h101;
        f_827 = v_h101;
        v_h101 = kappa_h101*(exp(Af_h101/(R*T))-exp(Ar_h101/(R*T)));

        // Re:h201

        var Af_h201: J_per_mol;
        var Ar_h201: J_per_mol;
        var v_h201: fmol_per_sec;

        Af_h201 = e_832;
        Ar_h201 = e_833;
        f_832 = v_h201;
        f_833 = v_h201;
        v_h201 = kappa_h201*(exp(Af_h201/(R*T))-exp(Ar_h201/(R*T)));

        // Re:h301

        var Af_h301: J_per_mol;
        var Ar_h301: J_per_mol;
        var v_h301: fmol_per_sec;

        Af_h301 = e_838;
        Ar_h301 = e_839;
        f_838 = v_h301;
        f_839 = v_h301;
        v_h301 = kappa_h301*(exp(Af_h301/(R*T))-exp(Ar_h301/(R*T)));

        // Re:m001

        var Af_m001: J_per_mol;
        var Ar_m001: J_per_mol;
        var v_m001: fmol_per_sec;

        Af_m001 = e_844;
        Ar_m001 = e_845;
        f_844 = v_m001;
        f_845 = v_m001;
        v_m001 = kappa_m001*(exp(Af_m001/(R*T))-exp(Ar_m001/(R*T)));

        // Re:m101

        var Af_m101: J_per_mol;
        var Ar_m101: J_per_mol;
        var v_m101: fmol_per_sec;

        Af_m101 = e_850;
        Ar_m101 = e_851;
        f_850 = v_m101;
        f_851 = v_m101;
        v_m101 = kappa_m101*(exp(Af_m101/(R*T))-exp(Ar_m101/(R*T)));

        // Re:m201

        var Af_m201: J_per_mol;
        var Ar_m201: J_per_mol;
        var v_m201: fmol_per_sec;

        Af_m201 = e_856;
        Ar_m201 = e_857;
        f_856 = v_m201;
        f_857 = v_m201;
        v_m201 = kappa_m201*(exp(Af_m201/(R*T))-exp(Ar_m201/(R*T)));

        // Re:m011

        var Af_m011: J_per_mol;
        var Ar_m011: J_per_mol;
        var v_m011: fmol_per_sec;

        Af_m011 = e_862;
        Ar_m011 = e_863;
        f_862 = v_m011;
        f_863 = v_m011;
        v_m011 = kappa_m011*(exp(Af_m011/(R*T))-exp(Ar_m011/(R*T)));

        // Re:m111

        var Af_m111: J_per_mol;
        var Ar_m111: J_per_mol;
        var v_m111: fmol_per_sec;

        Af_m111 = e_868;
        Ar_m111 = e_869;
        f_868 = v_m111;
        f_869 = v_m111;
        v_m111 = kappa_m111*(exp(Af_m111/(R*T))-exp(Ar_m111/(R*T)));

        // Re:m211

        var Af_m211: J_per_mol;
        var Ar_m211: J_per_mol;
        var v_m211: fmol_per_sec;

        Af_m211 = e_874;
        Ar_m211 = e_875;
        f_874 = v_m211;
        f_875 = v_m211;
        v_m211 = kappa_m211*(exp(Af_m211/(R*T))-exp(Ar_m211/(R*T)));

    enddef;
enddef;