Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 84 entries in the algebraic variable array.
C There are a total of 39 entries in each of the rate and state variable arrays.
C There are a total of 148 entries in the constant variable array.
C
C
C CONSTS(1) is t_period in component Istim_parameters (ms).
C CONSTS(2) is V_ICCrest in component Istim_parameters (mV).
C CONSTS(3) is V_ICCamp in component Istim_parameters (mV).
C CONSTS(4) is t_ICCpeak in component Istim_parameters (ms).
C CONSTS(5) is t_slope in component Istim_parameters (ms).
C CONSTS(6) is f_1 in component Istim_parameters (ms).
C CONSTS(7) is f_2 in component Istim_parameters (ms).
C CONSTS(8) is Gcouple in component Istim_parameters (nS).
C VOI is time in component time (ms).
C ALGBRC(14) is I_stim in component Periodic_IStim_protocol (pA).
C ALGBRC(9) is V_ICC in component Periodic_IStim_protocol (mV).
C STATES(1) is V in component Membrane_potential (mV).
C ALGBRC(83) is I_ion in component Membrane_potential (pA).
C ALGBRC(38) is I in component Ionic_currents (pA).
C ALGBRC(46) is I in component Ionic_currents (pA).
C ALGBRC(55) is I in component Ionic_currents (pA).
C ALGBRC(62) is I in component Ionic_currents (pA).
C ALGBRC(68) is I in component Ionic_currents (pA).
C ALGBRC(26) is I_NCX in component INCX (pA).
C ALGBRC(30) is I_NaK in component INaK (pA).
C ALGBRC(74) is I in component Ionic_currents (pA).
C ALGBRC(80) is I in component Ionic_currents (pA).
C STATES(2) is Cai in component Cai (mM).
C STATES(3) is Nai in component Nai (mM).
C STATES(4) is Ki in component Ki (mM).
C STATES(5) is O_CaL in component ICaL_channel_states (dimensionless).
C ALGBRC(18) is P_CaT in component ICaT (dimensionless).
C STATES(6) is O4 in component IBK_channel_states (dimensionless).
C ALGBRC(22) is P_Kv in component IKv (dimensionless).
C STATES(7) is O_Na in component INa_channel_states (dimensionless).
C ALGBRC(84) is norm in component ICaL_channel_states (dimensionless).
C ALGBRC(82) is norm in component IBK_channel_states (dimensionless).
C ALGBRC(57) is norm in component INa_channel_states (dimensionless).
C STATES(8) is X in component gating_kinetics (dimensionless).
C STATES(9) is X in component gating_kinetics (dimensionless).
C STATES(10) is X in component gating_kinetics (dimensionless).
C STATES(11) is X in component gating_kinetics (dimensionless).
C CONSTS(9) is g_CaL in component g_parameters (nS).
C CONSTS(10) is g_Kv in component g_parameters (nS).
C CONSTS(11) is g_BK in component g_parameters (nS).
C CONSTS(12) is g_Na in component g_parameters (nS).
C CONSTS(13) is V_init in component initial_conditions (mV).
C CONSTS(14) is Cm in component model_parameters (pF).
C CONSTS(15) is R in component model_parameters (J_per_K_mol).
C CONSTS(16) is T in component model_parameters (kelvin).
C CONSTS(17) is F in component model_parameters (C_per_mmol).
C CONSTS(18) is T0CaL in component model_parameters (kelvin).
C CONSTS(19) is T0CaT in component model_parameters (kelvin).
C CONSTS(20) is Q10Ca in component model_parameters (dimensionless).
C CONSTS(21) is Cao in component model_parameters (mM).
C CONSTS(22) is C0_init in component initial_conditions (dimensionless).
C CONSTS(23) is C1_init in component initial_conditions (dimensionless).
C CONSTS(24) is C2_init in component initial_conditions (dimensionless).
C CONSTS(25) is C3_init in component initial_conditions (dimensionless).
C CONSTS(26) is C0Ca_init in component initial_conditions (dimensionless).
C CONSTS(27) is C1Ca_init in component initial_conditions (dimensionless).
C CONSTS(28) is C2Ca_init in component initial_conditions (dimensionless).
C CONSTS(29) is C3Ca_init in component initial_conditions (dimensionless).
C CONSTS(30) is ICa_init in component initial_conditions (dimensionless).
C CONSTS(31) is IVs_init in component initial_conditions (dimensionless).
C CONSTS(32) is IVf_init in component initial_conditions (dimensionless).
C CONSTS(33) is IVfCa_init in component initial_conditions (dimensionless).
C CONSTS(34) is IVsCa_init in component initial_conditions (dimensionless).
C CONSTS(35) is O_CaL_init in component initial_conditions (dimensionless).
C CONSTS(36) is g_CaT in component model_parameters (nS).
C CONSTS(37) is d_CaT_init in component initial_conditions (dimensionless).
C CONSTS(38) is f_CaT_init in component initial_conditions (dimensionless).
C CONSTS(39) is tau_dCaT in component model_parameters (ms).
C CONSTS(40) is T0K in component model_parameters (kelvin).
C CONSTS(41) is Q10K in component model_parameters (dimensionless).
C CONSTS(42) is Ko in component model_parameters (mM).
C CONSTS(43) is x_Kv_init in component initial_conditions (dimensionless).
C CONSTS(44) is y_Kv_init in component initial_conditions (dimensionless).
C CONSTS(45) is tau_xKv in component model_parameters (ms).
C CONSTS(46) is tau_yKv in component model_parameters (ms).
C CONSTS(47) is BK_C0_init in component initial_conditions (dimensionless).
C CONSTS(48) is BK_C1_init in component initial_conditions (dimensionless).
C CONSTS(49) is BK_C2_init in component initial_conditions (dimensionless).
C CONSTS(50) is BK_C3_init in component initial_conditions (dimensionless).
C CONSTS(51) is BK_C4_init in component initial_conditions (dimensionless).
C CONSTS(52) is O0_init in component initial_conditions (dimensionless).
C CONSTS(53) is O1_init in component initial_conditions (dimensionless).
C CONSTS(54) is O2_init in component initial_conditions (dimensionless).
C CONSTS(55) is O3_init in component initial_conditions (dimensionless).
C CONSTS(56) is O4_init in component initial_conditions (dimensionless).
C CONSTS(57) is k_on in component model_parameters (per_ms).
C CONSTS(58) is k_off_C in component model_parameters (per_ms).
C CONSTS(59) is k_off_O in component model_parameters (per_ms).
C CONSTS(60) is T0Na in component model_parameters (kelvin).
C CONSTS(61) is Q10Na in component model_parameters (dimensionless).
C CONSTS(62) is Nao in component model_parameters (mM).
C CONSTS(63) is Na_C1_init in component initial_conditions (dimensionless).
C CONSTS(64) is Na_C2_init in component initial_conditions (dimensionless).
C CONSTS(65) is Na_C3_init in component initial_conditions (dimensionless).
C CONSTS(66) is I1_init in component initial_conditions (dimensionless).
C CONSTS(67) is I2_init in component initial_conditions (dimensionless).
C CONSTS(68) is O_Na_init in component initial_conditions (dimensionless).
C CONSTS(69) is K_mCa in component model_parameters (mM).
C CONSTS(70) is K_mNai in component model_parameters (mM).
C CONSTS(71) is k_sat in component model_parameters (dimensionless).
C CONSTS(72) is gamma in component model_parameters (dimensionless).
C CONSTS(73) is P_NCX in component model_parameters (pA).
C CONSTS(74) is K_mK in component model_parameters (mM).
C CONSTS(75) is K_mNa in component model_parameters (mM).
C CONSTS(76) is P_NaK in component model_parameters (pA).
C CONSTS(77) is g_NsNa in component model_parameters (nS).
C CONSTS(78) is g_NsK in component model_parameters (nS).
C CONSTS(79) is n_CRT in component model_parameters (dimensionless).
C CONSTS(80) is n_CaM in component model_parameters (dimensionless).
C CONSTS(81) is CRT_total in component model_parameters (mM).
C CONSTS(82) is CaM_total in component model_parameters (mM).
C CONSTS(83) is K_D_CRT in component model_parameters (mM).
C CONSTS(84) is K_D_CaM in component model_parameters (mM4).
C CONSTS(85) is V_cell in component model_parameters (litre).
C CONSTS(86) is Cai_total_init in component initial_conditions (mM).
C CONSTS(87) is Cai_init in component initial_conditions (mM).
C CONSTS(88) is Nai_init in component initial_conditions (mM).
C CONSTS(89) is Ki_init in component initial_conditions (mM).
C ALGBRC(1) is t in component Periodic_IStim_protocol (ms).
C CONSTS(90) is g_CaL in component model_parameters (nS).
C CONSTS(91) is g_Kv in component model_parameters (nS).
C CONSTS(92) is g_BK in component model_parameters (nS).
C CONSTS(93) is g_Na in component model_parameters (nS).
C CONSTS(94) is PO in component INS_Na (dimensionless).
C CONSTS(95) is PO in component INS_K (dimensionless).
C STATES(12) is Cai_total in component Cai (mM).
C ALGBRC(53) is Cai_total_rate in component Cai (mM_per_ms).
C CONSTS(96) is Cvt in component Cai (fM_per_s2mM_per_ms).
C CONSTS(97) is Cvt in component Nai (fM_per_s2mM_per_ms).
C CONSTS(98) is Cvt in component Ki (fM_per_s2mM_per_ms).
C ALGBRC(34) is E in component Nernst_potential (mV).
C STATES(13) is C0 in component ICaL_channel_states (dimensionless).
C STATES(14) is C1 in component ICaL_channel_states (dimensionless).
C STATES(15) is C2 in component ICaL_channel_states (dimensionless).
C STATES(16) is C3 in component ICaL_channel_states (dimensionless).
C STATES(17) is C0Ca in component ICaL_channel_states (dimensionless).
C STATES(18) is C1Ca in component ICaL_channel_states (dimensionless).
C STATES(19) is C2Ca in component ICaL_channel_states (dimensionless).
C STATES(20) is C3Ca in component ICaL_channel_states (dimensionless).
C STATES(21) is ICa in component ICaL_channel_states (dimensionless).
C STATES(22) is IVs in component ICaL_channel_states (dimensionless).
C STATES(23) is IVf in component ICaL_channel_states (dimensionless).
C STATES(24) is IVfCa in component ICaL_channel_states (dimensionless).
C STATES(25) is IVsCa in component ICaL_channel_states (dimensionless).
C ALGBRC(2) is alpha in component ICaL_channel_states (per_ms).
C ALGBRC(10) is beta in component ICaL_channel_states (per_ms).
C ALGBRC(15) is alpha_0 in component ICaL_channel_states (per_ms).
C ALGBRC(19) is alpha_1 in component ICaL_channel_states (per_ms).
C ALGBRC(23) is alpha_2 in component ICaL_channel_states (per_ms).
C ALGBRC(27) is alpha_3 in component ICaL_channel_states (per_ms).
C ALGBRC(32) is beta_0 in component ICaL_channel_states (per_ms).
C ALGBRC(37) is beta_1 in component ICaL_channel_states (per_ms).
C ALGBRC(41) is beta_2 in component ICaL_channel_states (per_ms).
C ALGBRC(45) is beta_3 in component ICaL_channel_states (per_ms).
C CONSTS(148) is delta in component ICaL_channel_states (per_ms).
C ALGBRC(81) is theta in component ICaL_channel_states (per_ms).
C ALGBRC(54) is phi_s in component ICaL_channel_states (per_ms).
C ALGBRC(49) is phi_f in component ICaL_channel_states (per_ms).
C ALGBRC(72) is omega_s in component ICaL_channel_states (per_ms).
C ALGBRC(69) is omega_f in component ICaL_channel_states (per_ms).
C ALGBRC(78) is omega_fs in component ICaL_channel_states (per_ms).
C ALGBRC(75) is omega_sf in component ICaL_channel_states (per_ms).
C ALGBRC(67) is psi_s in component ICaL_channel_states (per_ms).
C ALGBRC(64) is psi_f in component ICaL_channel_states (per_ms).
C ALGBRC(61) is xi_s in component ICaL_channel_states (per_ms).
C ALGBRC(58) is xi_f in component ICaL_channel_states (per_ms).
C CONSTS(134) is phi in component temperature_factor (dimensionless).
C ALGBRC(42) is E in component Nernst_potential (mV).
C ALGBRC(3) is X_inf in component d_inf (dimensionless).
C CONSTS(135) is phi in component temperature_factor (dimensionless).
C ALGBRC(4) is X_inf in component f_inf (dimensionless).
C ALGBRC(11) is tau in component tau_f (ms).
C CONSTS(136) is phi in component temperature_factor (dimensionless).
C ALGBRC(50) is E in component Nernst_potential (mV).
C ALGBRC(5) is X_inf in component x_inf (dimensionless).
C CONSTS(137) is phi in component temperature_factor (dimensionless).
C ALGBRC(6) is X_inf in component y_inf (dimensionless).
C CONSTS(138) is phi in component temperature_factor (dimensionless).
C ALGBRC(59) is E in component Nernst_potential (mV).
C STATES(26) is C0 in component IBK_channel_states (dimensionless).
C STATES(27) is C1 in component IBK_channel_states (dimensionless).
C STATES(28) is C2 in component IBK_channel_states (dimensionless).
C STATES(29) is C3 in component IBK_channel_states (dimensionless).
C STATES(30) is C4 in component IBK_channel_states (dimensionless).
C STATES(31) is O0 in component IBK_channel_states (dimensionless).
C STATES(32) is O1 in component IBK_channel_states (dimensionless).
C STATES(33) is O2 in component IBK_channel_states (dimensionless).
C STATES(34) is O3 in component IBK_channel_states (dimensionless).
C ALGBRC(7) is alpha in component IBK_channel_states (per_ms).
C ALGBRC(12) is beta in component IBK_channel_states (per_ms).
C ALGBRC(16) is k_C0O0 in component IBK_channel_states (per_ms).
C ALGBRC(20) is k_C1O1 in component IBK_channel_states (per_ms).
C ALGBRC(24) is k_C2O2 in component IBK_channel_states (per_ms).
C ALGBRC(28) is k_C3O3 in component IBK_channel_states (per_ms).
C ALGBRC(33) is k_C4O4 in component IBK_channel_states (per_ms).
C ALGBRC(35) is k_O0C0 in component IBK_channel_states (per_ms).
C ALGBRC(39) is k_O1C1 in component IBK_channel_states (per_ms).
C ALGBRC(43) is k_O2C2 in component IBK_channel_states (per_ms).
C ALGBRC(47) is k_O3C3 in component IBK_channel_states (per_ms).
C ALGBRC(51) is k_O4C4 in component IBK_channel_states (per_ms).
C ALGBRC(56) is k_C0C1 in component IBK_channel_states (per_ms).
C ALGBRC(60) is k_C1C2 in component IBK_channel_states (per_ms).
C ALGBRC(63) is k_C2C3 in component IBK_channel_states (per_ms).
C ALGBRC(66) is k_C3C4 in component IBK_channel_states (per_ms).
C CONSTS(139) is k_C4C3 in component IBK_channel_states (per_ms).
C CONSTS(140) is k_C3C2 in component IBK_channel_states (per_ms).
C CONSTS(141) is k_C2C1 in component IBK_channel_states (per_ms).
C CONSTS(142) is k_C1C0 in component IBK_channel_states (per_ms).
C ALGBRC(70) is k_O0O1 in component IBK_channel_states (per_ms).
C ALGBRC(73) is k_O1O2 in component IBK_channel_states (per_ms).
C ALGBRC(76) is k_O2O3 in component IBK_channel_states (per_ms).
C ALGBRC(79) is k_O3O4 in component IBK_channel_states (per_ms).
C CONSTS(143) is k_O4O3 in component IBK_channel_states (per_ms).
C CONSTS(144) is k_O3O2 in component IBK_channel_states (per_ms).
C CONSTS(145) is k_O2O1 in component IBK_channel_states (per_ms).
C CONSTS(146) is k_O1O0 in component IBK_channel_states (per_ms).
C ALGBRC(65) is E in component Nernst_potential (mV).
C STATES(35) is C1 in component INa_channel_states (dimensionless).
C STATES(36) is C2 in component INa_channel_states (dimensionless).
C STATES(37) is C3 in component INa_channel_states (dimensionless).
C STATES(38) is I1 in component INa_channel_states (dimensionless).
C STATES(39) is I2 in component INa_channel_states (dimensionless).
C ALGBRC(8) is k_I2I1 in component INa_channel_states (per_ms).
C ALGBRC(13) is k_I1O in component INa_channel_states (per_ms).
C ALGBRC(17) is k_OC1 in component INa_channel_states (per_ms).
C ALGBRC(21) is k_C1C2 in component INa_channel_states (per_ms).
C ALGBRC(25) is k_C2C3 in component INa_channel_states (per_ms).
C ALGBRC(29) is k_C3C2 in component INa_channel_states (per_ms).
C ALGBRC(31) is k_C2C1 in component INa_channel_states (per_ms).
C ALGBRC(36) is k_C1O in component INa_channel_states (per_ms).
C ALGBRC(40) is k_OI1 in component INa_channel_states (per_ms).
C ALGBRC(44) is k_I1I2 in component INa_channel_states (per_ms).
C ALGBRC(48) is k_I1C1 in component INa_channel_states (per_ms).
C ALGBRC(52) is k_C1I1 in component INa_channel_states (per_ms).
C CONSTS(147) is phi in component temperature_factor (dimensionless).
C ALGBRC(71) is E in component Nernst_potential (mV).
C ALGBRC(77) is E in component Nernst_potential (mV).
C CONSTS(99) is R in component constants (J_per_K_mol).
C CONSTS(100) is F in component constants (C_per_mmol).
C CONSTS(101) is z_Ca in component constants (dimensionless).
C CONSTS(102) is z_Na in component constants (dimensionless).
C CONSTS(103) is z_K in component constants (dimensionless).
C CONSTS(104) is R in component constants (J_per_K_mol).
C CONSTS(105) is F in component constants (C_per_mmol).
C CONSTS(106) is z_Ca in component constants (dimensionless).
C CONSTS(107) is z_Na in component constants (dimensionless).
C CONSTS(108) is z_K in component constants (dimensionless).
C CONSTS(109) is R in component constants (J_per_K_mol).
C CONSTS(110) is F in component constants (C_per_mmol).
C CONSTS(111) is z_K in component constants (dimensionless).
C CONSTS(112) is z_Na in component constants (dimensionless).
C CONSTS(113) is z_Ca in component constants (dimensionless).
C CONSTS(114) is R in component constants (J_per_K_mol).
C CONSTS(115) is F in component constants (C_per_mmol).
C CONSTS(116) is z_K in component constants (dimensionless).
C CONSTS(117) is z_Na in component constants (dimensionless).
C CONSTS(118) is z_Ca in component constants (dimensionless).
C CONSTS(119) is R in component constants (J_per_K_mol).
C CONSTS(120) is F in component constants (C_per_mmol).
C CONSTS(121) is z_Na in component constants (dimensionless).
C CONSTS(122) is z_Ca in component constants (dimensionless).
C CONSTS(123) is z_K in component constants (dimensionless).
C CONSTS(124) is R in component constants (J_per_K_mol).
C CONSTS(125) is F in component constants (C_per_mmol).
C CONSTS(126) is z_Na in component constants (dimensionless).
C CONSTS(127) is z_Ca in component constants (dimensionless).
C CONSTS(128) is z_K in component constants (dimensionless).
C CONSTS(129) is R in component constants (J_per_K_mol).
C CONSTS(130) is F in component constants (C_per_mmol).
C CONSTS(131) is z_K in component constants (dimensionless).
C CONSTS(132) is z_Na in component constants (dimensionless).
C CONSTS(133) is z_Ca in component constants (dimensionless).
C RATES(1) is d/dt V in component Membrane_potential (mV).
C RATES(12) is d/dt Cai_total in component Cai (mM).
C RATES(2) is d/dt Cai in component Cai (mM).
C RATES(3) is d/dt Nai in component Nai (mM).
C RATES(4) is d/dt Ki in component Ki (mM).
C RATES(16) is d/dt C3 in component ICaL_channel_states (dimensionless).
C RATES(15) is d/dt C2 in component ICaL_channel_states (dimensionless).
C RATES(14) is d/dt C1 in component ICaL_channel_states (dimensionless).
C RATES(13) is d/dt C0 in component ICaL_channel_states (dimensionless).
C RATES(20) is d/dt C3Ca in component ICaL_channel_states (dimensionless).
C RATES(19) is d/dt C2Ca in component ICaL_channel_states (dimensionless).
C RATES(18) is d/dt C1Ca in component ICaL_channel_states (dimensionless).
C RATES(17) is d/dt C0Ca in component ICaL_channel_states (dimensionless).
C RATES(5) is d/dt O_CaL in component ICaL_channel_states (dimensionless).
C RATES(21) is d/dt ICa in component ICaL_channel_states (dimensionless).
C RATES(22) is d/dt IVs in component ICaL_channel_states (dimensionless).
C RATES(23) is d/dt IVf in component ICaL_channel_states (dimensionless).
C RATES(25) is d/dt IVsCa in component ICaL_channel_states (dimensionless).
C RATES(24) is d/dt IVfCa in component ICaL_channel_states (dimensionless).
C RATES(8) is d/dt X in component gating_kinetics (dimensionless).
C RATES(9) is d/dt X in component gating_kinetics (dimensionless).
C RATES(10) is d/dt X in component gating_kinetics (dimensionless).
C RATES(11) is d/dt X in component gating_kinetics (dimensionless).
C RATES(30) is d/dt C4 in component IBK_channel_states (dimensionless).
C RATES(29) is d/dt C3 in component IBK_channel_states (dimensionless).
C RATES(28) is d/dt C2 in component IBK_channel_states (dimensionless).
C RATES(27) is d/dt C1 in component IBK_channel_states (dimensionless).
C RATES(26) is d/dt C0 in component IBK_channel_states (dimensionless).
C RATES(6) is d/dt O4 in component IBK_channel_states (dimensionless).
C RATES(34) is d/dt O3 in component IBK_channel_states (dimensionless).
C RATES(33) is d/dt O2 in component IBK_channel_states (dimensionless).
C RATES(32) is d/dt O1 in component IBK_channel_states (dimensionless).
C RATES(31) is d/dt O0 in component IBK_channel_states (dimensionless).
C RATES(37) is d/dt C3 in component INa_channel_states (dimensionless).
C RATES(36) is d/dt C2 in component INa_channel_states (dimensionless).
C RATES(35) is d/dt C1 in component INa_channel_states (dimensionless).
C RATES(7) is d/dt O_Na in component INa_channel_states (dimensionless).
C RATES(39) is d/dt I2 in component INa_channel_states (dimensionless).
C RATES(38) is d/dt I1 in component INa_channel_states (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 10000
CONSTS(2) = -57
CONSTS(3) = 33.5
CONSTS(4) = 300
CONSTS(5) = 600
CONSTS(6) = 12000
CONSTS(7) = 300
CONSTS(8) = 2.6
CONSTS(9) = 1.44
CONSTS(10) = 1.0217
CONSTS(11) = 80
CONSTS(12) = 25.1
CONSTS(13) = -73.5049651455872
CONSTS(14) = 50
CONSTS(15) = 8.314
CONSTS(16) = 310
CONSTS(17) = 96.48534
CONSTS(18) = 310
CONSTS(19) = 297
CONSTS(20) = 2.1
CONSTS(21) = 2
CONSTS(22) = 0.815464741971086
CONSTS(23) = 0.152399266235657
CONSTS(24) = 0.0106805060777161
CONSTS(25) = 0.000332673548872087
CONSTS(26) = 0.0175888495282545
CONSTS(27) = 0.00328711668724504
CONSTS(28) = 0.000230369020877669
CONSTS(29) = 7.1754726923539e-6
CONSTS(30) = 8.38123983500905e-8
CONSTS(31) = 1.1193313274705e-6
CONSTS(32) = 4.0998751301597e-6
CONSTS(33) = 8.84306615061238e-8
CONSTS(34) = 2.41429816075123e-8
CONSTS(35) = 3.88576045134351e-6
CONSTS(36) = 0.0425
CONSTS(37) = 0.0791635737410974
CONSTS(38) = 0.377831534375835
CONSTS(39) = 1.9508
CONSTS(40) = 297
CONSTS(41) = 3.1
CONSTS(42) = 5.4
CONSTS(43) = 0.14714161078933
CONSTS(44) = 0.99994773314105
CONSTS(45) = 4.7803
CONSTS(46) = 763.7564
CONSTS(47) = 0.48379087935899
CONSTS(48) = 0.385183559520031
CONSTS(49) = 0.115002824567753
CONSTS(50) = 0.0152602714149774
CONSTS(51) = 0.000759264410974374
CONSTS(52) = 6.94960798375172e-7
CONSTS(53) = 5.55636826398253e-8
CONSTS(54) = 2.85143702125325e-8
CONSTS(55) = 1.59832806123435e-6
CONSTS(56) = 1.82113764497095e-6
CONSTS(57) = 40633
CONSTS(58) = 11
CONSTS(59) = 1.1
CONSTS(60) = 297
CONSTS(61) = 2.45
CONSTS(62) = 140
CONSTS(63) = 0.0119443135223679
CONSTS(64) = 0.0109545368437155
CONSTS(65) = 0.973782548650071
CONSTS(66) = 0.000126882921013389
CONSTS(67) = 0.00318975045717667
CONSTS(68) = 1.96760342050475e-6
CONSTS(69) = 1.38
CONSTS(70) = 87.5
CONSTS(71) = 0.1
CONSTS(72) = 0.35
CONSTS(73) = 1992.1865
CONSTS(74) = 1
CONSTS(75) = 40
CONSTS(76) = 16.197
CONSTS(77) = 0.022488
CONSTS(78) = 0.017512
CONSTS(79) = 1
CONSTS(80) = 4
CONSTS(81) = 0.034
CONSTS(82) = 0.012
CONSTS(83) = 0.0009
CONSTS(84) = 0.0001
CONSTS(85) = 3.5e-12
CONSTS(86) = 0.004914
CONSTS(87) = 5.38843941249284e-5
CONSTS(88) = 10.5731241425458
CONSTS(89) = 153.604280337996
CONSTS(90) = 1.44
CONSTS(91) = 1.0217
CONSTS(92) = 80
CONSTS(93) = 25.1
CONSTS(94) = 1
CONSTS(95) = 1
CONSTS(96) = 1e-15
CONSTS(97) = 1e-15
CONSTS(98) = 1e-15
CONSTS(99) = 8.314
CONSTS(100) = 96.48534
CONSTS(101) = 2
CONSTS(102) = 1
CONSTS(103) = 1
CONSTS(104) = 8.314
CONSTS(105) = 96.48534
CONSTS(106) = 2
CONSTS(107) = 1
CONSTS(108) = 1
CONSTS(109) = 8.314
CONSTS(110) = 96.48534
CONSTS(111) = 1
CONSTS(112) = 1
CONSTS(113) = 2
CONSTS(114) = 8.314
CONSTS(115) = 96.48534
CONSTS(116) = 1
CONSTS(117) = 1
CONSTS(118) = 2
CONSTS(119) = 8.314
CONSTS(120) = 96.48534
CONSTS(121) = 1
CONSTS(122) = 2
CONSTS(123) = 1
CONSTS(124) = 8.314
CONSTS(125) = 96.48534
CONSTS(126) = 1
CONSTS(127) = 2
CONSTS(128) = 1
CONSTS(129) = 8.314
CONSTS(130) = 96.48534
CONSTS(131) = 1
CONSTS(132) = 1
CONSTS(133) = 2
CONSTS(134) = CONSTS(20) ** (CONSTS(16) - CONSTS(18))/10.0000
CONSTS(135) = CONSTS(20) ** (CONSTS(16) - CONSTS(19))/10.0000
CONSTS(136) = CONSTS(20) ** (CONSTS(16) - CONSTS(19))/10.0000
CONSTS(137) = CONSTS(41) ** (CONSTS(16) - CONSTS(40))/10.0000
CONSTS(138) = CONSTS(41) ** (CONSTS(16) - CONSTS(40))/10.0000
CONSTS(139) = 4.00000*CONSTS(58)
CONSTS(140) = 3.00000*CONSTS(58)
CONSTS(141) = 2.00000*CONSTS(58)
CONSTS(142) = 1.00000*CONSTS(58)
CONSTS(143) = 4.00000*CONSTS(59)
CONSTS(144) = 3.00000*CONSTS(59)
CONSTS(145) = 2.00000*CONSTS(59)
CONSTS(146) = 1.00000*CONSTS(59)
CONSTS(147) = CONSTS(61) ** (CONSTS(16) - CONSTS(60))/10.0000
CONSTS(148) = CONSTS(134)*0.0100000
STATES(1) = CONSTS(13)
STATES(2) = CONSTS(87)
STATES(3) = CONSTS(88)
STATES(4) = CONSTS(89)
STATES(5) = CONSTS(35)
STATES(6) = CONSTS(56)
STATES(7) = CONSTS(68)
STATES(8) = CONSTS(37)
STATES(9) = CONSTS(38)
STATES(10) = CONSTS(43)
STATES(11) = CONSTS(44)
STATES(12) = CONSTS(86)
STATES(13) = CONSTS(22)
STATES(14) = CONSTS(23)
STATES(15) = CONSTS(24)
STATES(16) = CONSTS(25)
STATES(17) = CONSTS(26)
STATES(18) = CONSTS(27)
STATES(19) = CONSTS(28)
STATES(20) = CONSTS(29)
STATES(21) = CONSTS(30)
STATES(22) = CONSTS(31)
STATES(23) = CONSTS(32)
STATES(24) = CONSTS(33)
STATES(25) = CONSTS(33)
STATES(26) = CONSTS(47)
STATES(27) = CONSTS(48)
STATES(28) = CONSTS(49)
STATES(29) = CONSTS(50)
STATES(30) = CONSTS(51)
STATES(31) = CONSTS(52)
STATES(32) = CONSTS(53)
STATES(33) = CONSTS(54)
STATES(34) = CONSTS(55)
STATES(35) = CONSTS(63)
STATES(36) = CONSTS(64)
STATES(37) = CONSTS(65)
STATES(38) = CONSTS(66)
STATES(39) = CONSTS(67)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+60.5000)/5.30000))
RATES(8) = ((ALGBRC(3) - STATES(8))/CONSTS(39))*CONSTS(135)
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+43.0000)/17.3600))
RATES(10) = ((ALGBRC(5) - STATES(10))/CONSTS(45))*CONSTS(137)
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1) - 44.9000)/12.0096))
RATES(11) = ((ALGBRC(6) - STATES(11))/CONSTS(46))*CONSTS(138)
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.5000)/4.00000))
ALGBRC(11) = 0.381170*(8.60000+ 14.7000*EXP(- STATES(1)+50.0000 ** 2.00000/( 900.000*1.00000)))
RATES(9) = ((ALGBRC(4) - STATES(9))/ALGBRC(11))*CONSTS(136)
ALGBRC(34) = (( CONSTS(99)*CONSTS(16))/( CONSTS(101)*CONSTS(100)))*log(CONSTS(21)/STATES(2))
ALGBRC(38) = CONSTS(9)*STATES(5)*(STATES(1) - ALGBRC(34))
ALGBRC(18) = STATES(8)*STATES(9)
ALGBRC(42) = (( CONSTS(104)*CONSTS(16))/( CONSTS(106)*CONSTS(105)))*log(CONSTS(21)/STATES(2))
ALGBRC(46) = CONSTS(36)*ALGBRC(18)*(STATES(1) - ALGBRC(42))
ALGBRC(26) = ( CONSTS(73)*( EXP(( CONSTS(72)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))*STATES(3) ** 3.00000*CONSTS(21) - 2.50000*EXP(( (CONSTS(72) - 1.00000)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))*CONSTS(62) ** 3.00000*STATES(2)))/( (CONSTS(70) ** 3.00000+CONSTS(62) ** 3.00000)*(CONSTS(69)+CONSTS(21))*(1.00000+ CONSTS(71)*EXP(( (CONSTS(72) - 1.00000)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))))
RATES(12) = ( - ((ALGBRC(38)+ALGBRC(46)) - 2.00000*ALGBRC(26))*CONSTS(96))/( 2.00000*CONSTS(17)*CONSTS(85))
ALGBRC(53) = ( - ((ALGBRC(38)+ALGBRC(46)) - 2.00000*ALGBRC(26))*CONSTS(96))/( 2.00000*CONSTS(17)*CONSTS(85))
RATES(2) = ALGBRC(53)/(1.00000+( CONSTS(79)*CONSTS(81)*CONSTS(83)*STATES(2) ** CONSTS(79) - 1.00000)/STATES(2) ** CONSTS(79)+CONSTS(83) ** 2.00000+( CONSTS(80)*CONSTS(82)*CONSTS(84)*STATES(2) ** CONSTS(80) - 1.00000)/STATES(2) ** CONSTS(80)+CONSTS(84) ** 2.00000)
ALGBRC(57) = STATES(7)+STATES(35)+STATES(36)+STATES(37)+STATES(38)+STATES(39)
ALGBRC(25) = CONSTS(147)*0.554320*EXP(- 0.0990740+ 0.0364410*STATES(1))
ALGBRC(29) = CONSTS(147)*0.000525480*EXP(- 0.0691020+ 0.00319450*STATES(1))
RATES(37) = ( - ALGBRC(29)*STATES(37))/ALGBRC(57)+( ALGBRC(25)*STATES(36))/ALGBRC(57)
ALGBRC(21) = CONSTS(147)*3.15660*EXP(0.363520+ 0.0771930*STATES(1))
ALGBRC(31) = CONSTS(147)*1.44960*EXP(- 0.156600+ 0.0583530*STATES(1))
RATES(36) = ( - (ALGBRC(31)+ALGBRC(25))*STATES(36))/ALGBRC(57)+( ALGBRC(21)*STATES(35))/ALGBRC(57)+( ALGBRC(29)*STATES(37))/ALGBRC(57)
ALGBRC(17) = CONSTS(147)*2.39100*EXP(- 13.3350 - 0.252890*STATES(1))
ALGBRC(36) = CONSTS(147)*1.53290*EXP(0.00931930+ 0.0410750*STATES(1))
ALGBRC(48) = CONSTS(147)*1.90460*EXP(- 2.48400+ 0.0204060*STATES(1))
ALGBRC(52) = CONSTS(147)*0.000216880*EXP(- 0.0634380+ 0.00466830*STATES(1))
RATES(35) = ( - (ALGBRC(21)+ALGBRC(36)+ALGBRC(52))*STATES(35))/ALGBRC(57)+( ALGBRC(17)*STATES(7))/ALGBRC(57)+( ALGBRC(31)*STATES(36))/ALGBRC(57)+( ALGBRC(48)*STATES(38))/ALGBRC(57)
ALGBRC(13) = CONSTS(147)*0.120520*EXP(- 9.60280+ 0.0830250*STATES(1))
ALGBRC(40) = CONSTS(147)*1.61640*EXP(0.307630+ 0.00605350*STATES(1))
RATES(7) = ( - (ALGBRC(17)+ALGBRC(40))*STATES(7))/ALGBRC(57)+( ALGBRC(36)*STATES(35))/ALGBRC(57)+( ALGBRC(13)*STATES(38))/ALGBRC(57)
ALGBRC(8) = CONSTS(147)*0.00392390*EXP(2.67930+ 0.00614680*STATES(1))
ALGBRC(44) = CONSTS(147)*0.0277350*EXP(0.0514900 - 0.0468650*STATES(1))
RATES(39) = ( - ALGBRC(8)*STATES(39))/ALGBRC(57)+( ALGBRC(44)*STATES(38))/ALGBRC(57)
RATES(38) = ( - (ALGBRC(13)+ALGBRC(44)+ALGBRC(48))*STATES(38))/ALGBRC(57)+( ALGBRC(8)*STATES(39))/ALGBRC(57)+( ALGBRC(52)*STATES(35))/ALGBRC(57)+( ALGBRC(40)*STATES(7))/ALGBRC(57)
ALGBRC(65) = (( CONSTS(119)*CONSTS(16))/( CONSTS(121)*CONSTS(120)))*log(CONSTS(62)/STATES(3))
ALGBRC(68) = CONSTS(12)*STATES(7)*(STATES(1) - ALGBRC(65))
ALGBRC(30) = ( CONSTS(76)*CONSTS(42)*STATES(3))/( (CONSTS(42)+CONSTS(74))*(CONSTS(75)+STATES(3))*(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))+ 0.0353000*EXP(( - STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))))
ALGBRC(71) = (( CONSTS(124)*CONSTS(16))/( CONSTS(126)*CONSTS(125)))*log(CONSTS(62)/STATES(3))
ALGBRC(74) = CONSTS(77)*CONSTS(94)*(STATES(1) - ALGBRC(71))
RATES(3) = ( - (ALGBRC(68)+ 3.00000*ALGBRC(30)+ 3.00000*ALGBRC(26)+ALGBRC(74))*CONSTS(97))/( CONSTS(17)*CONSTS(85))
ALGBRC(1) = INT(MOD(VOI, CONSTS(1)))
ALGBRC(9) = TERNRY(ALGBRC(1).LT.CONSTS(4), CONSTS(2)+( CONSTS(3)*ALGBRC(1))/CONSTS(7), CONSTS(2)+( CONSTS(3)*(1.00000+EXP(- CONSTS(6)/( 2.00000*CONSTS(5)))))/(1.00000+EXP(((ALGBRC(1) - CONSTS(7)) - 0.500000*CONSTS(6))/CONSTS(5))))
ALGBRC(14) = CONSTS(8)*(STATES(1) - ALGBRC(9))
ALGBRC(22) = STATES(10)*STATES(11)
ALGBRC(50) = (( CONSTS(109)*CONSTS(16))/( CONSTS(111)*CONSTS(110)))*log(CONSTS(42)/STATES(4))
ALGBRC(55) = CONSTS(10)*ALGBRC(22)*(STATES(1) - ALGBRC(50))
ALGBRC(59) = (( CONSTS(114)*CONSTS(16))/( CONSTS(116)*CONSTS(115)))*log(CONSTS(42)/STATES(4))
ALGBRC(62) = CONSTS(11)*STATES(6)*(STATES(1) - ALGBRC(59))
ALGBRC(77) = (( CONSTS(129)*CONSTS(16))/( CONSTS(131)*CONSTS(130)))*log(CONSTS(42)/STATES(4))
ALGBRC(80) = CONSTS(78)*CONSTS(95)*(STATES(1) - ALGBRC(77))
RATES(4) = ( - (((ALGBRC(55)+ALGBRC(62)+ALGBRC(14)) - 2.00000*ALGBRC(30))+ALGBRC(80))*CONSTS(98))/( CONSTS(17)*CONSTS(85))
ALGBRC(82) = STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)+STATES(31)+STATES(32)+STATES(33)+STATES(34)+STATES(6)
ALGBRC(7) = 1.00000*EXP(( 8.47188*STATES(1))/( 1.00000*CONSTS(16)))
ALGBRC(33) = 0.0443240*ALGBRC(7)
ALGBRC(12) = 1.00000*EXP(( - 7.77556*STATES(1))/( 1.00000*CONSTS(16)))
ALGBRC(51) = 0.000106000*ALGBRC(12)
ALGBRC(66) = 1.00000*CONSTS(57)*STATES(2)
RATES(30) = ( - (CONSTS(139)+ALGBRC(33))*STATES(30))/ALGBRC(82)+( ALGBRC(66)*STATES(29))/ALGBRC(82)+( ALGBRC(51)*STATES(6))/ALGBRC(82)
ALGBRC(28) = 0.000781000*ALGBRC(7)
ALGBRC(47) = 0.0953120*ALGBRC(12)
ALGBRC(63) = 2.00000*CONSTS(57)*STATES(2)
RATES(29) = ( - (CONSTS(140)+ALGBRC(28)+ALGBRC(66))*STATES(29))/ALGBRC(82)+( ALGBRC(63)*STATES(28))/ALGBRC(82)+( ALGBRC(47)*STATES(34))/ALGBRC(82)+( CONSTS(139)*STATES(30))/ALGBRC(82)
ALGBRC(24) = 2.81000e-05*ALGBRC(7)
ALGBRC(43) = 32.6594*ALGBRC(12)
ALGBRC(60) = 3.00000*CONSTS(57)*STATES(2)
RATES(28) = ( - (CONSTS(141)+ALGBRC(24)+ALGBRC(63))*STATES(28))/ALGBRC(82)+( ALGBRC(60)*STATES(27))/ALGBRC(82)+( ALGBRC(43)*STATES(33))/ALGBRC(82)+( CONSTS(140)*STATES(29))/ALGBRC(82)
ALGBRC(20) = 0.000869000*ALGBRC(7)
ALGBRC(39) = 144.174*ALGBRC(12)
ALGBRC(56) = 4.00000*CONSTS(57)*STATES(2)
RATES(27) = ( - (CONSTS(142)+ALGBRC(20)+ALGBRC(60))*STATES(27))/ALGBRC(82)+( ALGBRC(56)*STATES(26))/ALGBRC(82)+( ALGBRC(39)*STATES(32))/ALGBRC(82)+( CONSTS(141)*STATES(28))/ALGBRC(82)
ALGBRC(16) = 0.0216200*ALGBRC(7)
ALGBRC(35) = 318.108*ALGBRC(12)
RATES(26) = ( - (ALGBRC(56)+ALGBRC(16))*STATES(26))/ALGBRC(82)+( CONSTS(142)*STATES(27))/ALGBRC(82)+( ALGBRC(35)*STATES(31))/ALGBRC(82)
ALGBRC(79) = 1.00000*CONSTS(57)*STATES(2)
RATES(6) = ( - (CONSTS(143)+ALGBRC(51))*STATES(6))/ALGBRC(82)+( ALGBRC(79)*STATES(34))/ALGBRC(82)+( ALGBRC(33)*STATES(30))/ALGBRC(82)
ALGBRC(76) = 2.00000*CONSTS(57)*STATES(2)
RATES(34) = ( - (CONSTS(144)+ALGBRC(47)+ALGBRC(79))*STATES(34))/ALGBRC(82)+( ALGBRC(76)*STATES(33))/ALGBRC(82)+( ALGBRC(28)*STATES(29))/ALGBRC(82)+( CONSTS(143)*STATES(6))/ALGBRC(82)
ALGBRC(73) = 3.00000*CONSTS(57)*STATES(2)
RATES(33) = ( - (CONSTS(145)+ALGBRC(43)+ALGBRC(76))*STATES(33))/ALGBRC(82)+( ALGBRC(73)*STATES(32))/ALGBRC(82)+( ALGBRC(24)*STATES(28))/ALGBRC(82)+( CONSTS(144)*STATES(34))/ALGBRC(82)
ALGBRC(70) = 4.00000*CONSTS(57)*STATES(2)
RATES(32) = ( - (CONSTS(146)+ALGBRC(39)+ALGBRC(73))*STATES(32))/ALGBRC(82)+( ALGBRC(70)*STATES(31))/ALGBRC(82)+( ALGBRC(20)*STATES(27))/ALGBRC(82)+( CONSTS(145)*STATES(33))/ALGBRC(82)
RATES(31) = ( - (ALGBRC(70)+ALGBRC(35))*STATES(31))/ALGBRC(82)+( CONSTS(146)*STATES(32))/ALGBRC(82)+( ALGBRC(16)*STATES(26))/ALGBRC(82)
ALGBRC(83) = ALGBRC(38)+ALGBRC(46)+ALGBRC(55)+ALGBRC(62)+ALGBRC(68)+ALGBRC(26)+ALGBRC(30)+ALGBRC(74)+ALGBRC(80)
RATES(1) = - (ALGBRC(14)+ALGBRC(83))/CONSTS(14)
ALGBRC(84) = STATES(16)+STATES(15)+STATES(14)+STATES(13)+STATES(20)+STATES(19)+STATES(18)+STATES(17)+STATES(5)+STATES(21)+STATES(22)+STATES(23)+STATES(25)+STATES(24)
ALGBRC(2) = CONSTS(134)*0.731000*EXP(STATES(1)/30.0000)
ALGBRC(23) = 2.00000*ALGBRC(2)
ALGBRC(27) = 1.00000*ALGBRC(2)
ALGBRC(10) = CONSTS(134)*0.214900*EXP(- STATES(1)/40.0000)
ALGBRC(41) = 3.00000*ALGBRC(10)
ALGBRC(45) = 4.00000*ALGBRC(10)
ALGBRC(81) = ( CONSTS(134)*4.00000)/(1.00000+1.00000/STATES(2))
ALGBRC(54) = CONSTS(134)*0.00595600*EXP(- STATES(1)/40.0000)
ALGBRC(49) = CONSTS(134)*0.474200*EXP(STATES(1)/10.0000)
ALGBRC(67) = CONSTS(134)*0.00232000*EXP(- STATES(1)/280.000)
ALGBRC(61) = CONSTS(134)*0.0121300*EXP(STATES(1)/500.000)
ALGBRC(72) = ( ALGBRC(45)*ALGBRC(61)*ALGBRC(54))/( ALGBRC(27)*ALGBRC(67))
ALGBRC(64) = CONSTS(134)*0.0219700*EXP(STATES(1)/500.000)
ALGBRC(58) = CONSTS(134)*0.0140700*EXP(- STATES(1)/300.000)
ALGBRC(69) = ( ALGBRC(45)*ALGBRC(58)*ALGBRC(49))/( ALGBRC(27)*ALGBRC(64))
RATES(16) = (( ALGBRC(23)*STATES(15))/ALGBRC(84) - ( (ALGBRC(27)+ALGBRC(41)+ALGBRC(54)+ALGBRC(49)+ALGBRC(81))*STATES(16))/ALGBRC(84))+( ALGBRC(45)*STATES(5))/ALGBRC(84)+( ALGBRC(69)*STATES(23))/ALGBRC(84)+( ALGBRC(72)*STATES(22))/ALGBRC(84)+( CONSTS(148)*STATES(20))/ALGBRC(84)
ALGBRC(19) = 3.00000*ALGBRC(2)
ALGBRC(37) = 2.00000*ALGBRC(10)
RATES(15) = (( ALGBRC(19)*STATES(14))/ALGBRC(84) - ( (ALGBRC(23)+ALGBRC(37)+ALGBRC(81))*STATES(15))/ALGBRC(84))+( ALGBRC(41)*STATES(16))/ALGBRC(84)+( CONSTS(148)*STATES(19))/ALGBRC(84)
ALGBRC(15) = 4.00000*ALGBRC(2)
ALGBRC(32) = 1.00000*ALGBRC(10)
RATES(14) = (( ALGBRC(15)*STATES(13))/ALGBRC(84) - ( (ALGBRC(19)+ALGBRC(32)+ALGBRC(81))*STATES(14))/ALGBRC(84))+( ALGBRC(37)*STATES(15))/ALGBRC(84)+( CONSTS(148)*STATES(18))/ALGBRC(84)
RATES(13) = ( - (ALGBRC(15)+ALGBRC(81))*STATES(13))/ALGBRC(84)+( ALGBRC(32)*STATES(14))/ALGBRC(84)+( CONSTS(148)*STATES(17))/ALGBRC(84)
RATES(20) = ((( ALGBRC(81)*STATES(16))/ALGBRC(84)+( ALGBRC(23)*STATES(19))/ALGBRC(84)) - ( (ALGBRC(27)+ALGBRC(41)+ALGBRC(49)+ALGBRC(54)+CONSTS(148))*STATES(20))/ALGBRC(84))+( ALGBRC(45)*STATES(21))/ALGBRC(84)+( ALGBRC(69)*STATES(24))/ALGBRC(84)+( ALGBRC(72)*STATES(25))/ALGBRC(84)
RATES(19) = ((( ALGBRC(81)*STATES(15))/ALGBRC(84)+( ALGBRC(19)*STATES(18))/ALGBRC(84)) - ( (ALGBRC(23)+ALGBRC(37)+CONSTS(148))*STATES(19))/ALGBRC(84))+( ALGBRC(41)*STATES(20))/ALGBRC(84)
RATES(18) = ((( ALGBRC(81)*STATES(14))/ALGBRC(84)+( ALGBRC(15)*STATES(17))/ALGBRC(84)) - ( (ALGBRC(19)+ALGBRC(32)+CONSTS(148))*STATES(18))/ALGBRC(84))+( ALGBRC(37)*STATES(19))/ALGBRC(84)
RATES(17) = (( ALGBRC(81)*STATES(13))/ALGBRC(84) - ( (ALGBRC(15)+CONSTS(148))*STATES(17))/ALGBRC(84))+( ALGBRC(32)*STATES(18))/ALGBRC(84)
RATES(5) = (( ALGBRC(27)*STATES(16))/ALGBRC(84) - ( (ALGBRC(45)+ALGBRC(64)+ALGBRC(67)+ALGBRC(81))*STATES(5))/ALGBRC(84))+( ALGBRC(58)*STATES(23))/ALGBRC(84)+( ALGBRC(61)*STATES(22))/ALGBRC(84)+( CONSTS(148)*STATES(21))/ALGBRC(84)
RATES(21) = ((( ALGBRC(81)*STATES(5))/ALGBRC(84)+( ALGBRC(27)*STATES(20))/ALGBRC(84)) - ( (ALGBRC(45)+ALGBRC(64)+ALGBRC(67)+CONSTS(148))*STATES(21))/ALGBRC(84))+( ALGBRC(58)*STATES(24))/ALGBRC(84)+( ALGBRC(61)*STATES(25))/ALGBRC(84)
ALGBRC(78) = ALGBRC(67)
ALGBRC(75) = ( ALGBRC(61)*ALGBRC(64))/ALGBRC(58)
RATES(22) = ((( ALGBRC(54)*STATES(16))/ALGBRC(84)+( ALGBRC(67)*STATES(5))/ALGBRC(84)+( ALGBRC(78)*STATES(23))/ALGBRC(84)) - ( (ALGBRC(72)+ALGBRC(61)+ALGBRC(75)+ALGBRC(81))*STATES(22))/ALGBRC(84))+( CONSTS(148)*STATES(25))/ALGBRC(84)
RATES(23) = ((( ALGBRC(49)*STATES(16))/ALGBRC(84)+( ALGBRC(64)*STATES(5))/ALGBRC(84)) - ( (ALGBRC(69)+ALGBRC(58)+ALGBRC(78)+ALGBRC(81))*STATES(23))/ALGBRC(84))+( ALGBRC(75)*STATES(22))/ALGBRC(84)+( CONSTS(148)*STATES(24))/ALGBRC(84)
RATES(25) = (( ALGBRC(81)*STATES(22))/ALGBRC(84)+( ALGBRC(54)*STATES(20))/ALGBRC(84)+( ALGBRC(67)*STATES(21))/ALGBRC(84)+( ALGBRC(78)*STATES(24))/ALGBRC(84)) - ( (ALGBRC(72)+ALGBRC(61)+ALGBRC(75)+CONSTS(148))*STATES(25))/ALGBRC(84)
RATES(24) = ((( ALGBRC(81)*STATES(23))/ALGBRC(84)+( ALGBRC(49)*STATES(20))/ALGBRC(84)+( ALGBRC(64)*STATES(21))/ALGBRC(84)) - ( (ALGBRC(69)+ALGBRC(58)+ALGBRC(78)+CONSTS(148))*STATES(24))/ALGBRC(84))+( ALGBRC(75)*STATES(25))/ALGBRC(84)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+60.5000)/5.30000))
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+43.0000)/17.3600))
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1) - 44.9000)/12.0096))
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.5000)/4.00000))
ALGBRC(11) = 0.381170*(8.60000+ 14.7000*EXP(- STATES(1)+50.0000 ** 2.00000/( 900.000*1.00000)))
ALGBRC(34) = (( CONSTS(99)*CONSTS(16))/( CONSTS(101)*CONSTS(100)))*log(CONSTS(21)/STATES(2))
ALGBRC(38) = CONSTS(9)*STATES(5)*(STATES(1) - ALGBRC(34))
ALGBRC(18) = STATES(8)*STATES(9)
ALGBRC(42) = (( CONSTS(104)*CONSTS(16))/( CONSTS(106)*CONSTS(105)))*log(CONSTS(21)/STATES(2))
ALGBRC(46) = CONSTS(36)*ALGBRC(18)*(STATES(1) - ALGBRC(42))
ALGBRC(26) = ( CONSTS(73)*( EXP(( CONSTS(72)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))*STATES(3) ** 3.00000*CONSTS(21) - 2.50000*EXP(( (CONSTS(72) - 1.00000)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))*CONSTS(62) ** 3.00000*STATES(2)))/( (CONSTS(70) ** 3.00000+CONSTS(62) ** 3.00000)*(CONSTS(69)+CONSTS(21))*(1.00000+ CONSTS(71)*EXP(( (CONSTS(72) - 1.00000)*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))))
ALGBRC(53) = ( - ((ALGBRC(38)+ALGBRC(46)) - 2.00000*ALGBRC(26))*CONSTS(96))/( 2.00000*CONSTS(17)*CONSTS(85))
ALGBRC(57) = STATES(7)+STATES(35)+STATES(36)+STATES(37)+STATES(38)+STATES(39)
ALGBRC(25) = CONSTS(147)*0.554320*EXP(- 0.0990740+ 0.0364410*STATES(1))
ALGBRC(29) = CONSTS(147)*0.000525480*EXP(- 0.0691020+ 0.00319450*STATES(1))
ALGBRC(21) = CONSTS(147)*3.15660*EXP(0.363520+ 0.0771930*STATES(1))
ALGBRC(31) = CONSTS(147)*1.44960*EXP(- 0.156600+ 0.0583530*STATES(1))
ALGBRC(17) = CONSTS(147)*2.39100*EXP(- 13.3350 - 0.252890*STATES(1))
ALGBRC(36) = CONSTS(147)*1.53290*EXP(0.00931930+ 0.0410750*STATES(1))
ALGBRC(48) = CONSTS(147)*1.90460*EXP(- 2.48400+ 0.0204060*STATES(1))
ALGBRC(52) = CONSTS(147)*0.000216880*EXP(- 0.0634380+ 0.00466830*STATES(1))
ALGBRC(13) = CONSTS(147)*0.120520*EXP(- 9.60280+ 0.0830250*STATES(1))
ALGBRC(40) = CONSTS(147)*1.61640*EXP(0.307630+ 0.00605350*STATES(1))
ALGBRC(8) = CONSTS(147)*0.00392390*EXP(2.67930+ 0.00614680*STATES(1))
ALGBRC(44) = CONSTS(147)*0.0277350*EXP(0.0514900 - 0.0468650*STATES(1))
ALGBRC(65) = (( CONSTS(119)*CONSTS(16))/( CONSTS(121)*CONSTS(120)))*log(CONSTS(62)/STATES(3))
ALGBRC(68) = CONSTS(12)*STATES(7)*(STATES(1) - ALGBRC(65))
ALGBRC(30) = ( CONSTS(76)*CONSTS(42)*STATES(3))/( (CONSTS(42)+CONSTS(74))*(CONSTS(75)+STATES(3))*(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))+ 0.0353000*EXP(( - STATES(1)*CONSTS(17))/( CONSTS(15)*CONSTS(16)))))
ALGBRC(71) = (( CONSTS(124)*CONSTS(16))/( CONSTS(126)*CONSTS(125)))*log(CONSTS(62)/STATES(3))
ALGBRC(74) = CONSTS(77)*CONSTS(94)*(STATES(1) - ALGBRC(71))
ALGBRC(1) = INT(MOD(VOI, CONSTS(1)))
ALGBRC(9) = TERNRY(ALGBRC(1).LT.CONSTS(4), CONSTS(2)+( CONSTS(3)*ALGBRC(1))/CONSTS(7), CONSTS(2)+( CONSTS(3)*(1.00000+EXP(- CONSTS(6)/( 2.00000*CONSTS(5)))))/(1.00000+EXP(((ALGBRC(1) - CONSTS(7)) - 0.500000*CONSTS(6))/CONSTS(5))))
ALGBRC(14) = CONSTS(8)*(STATES(1) - ALGBRC(9))
ALGBRC(22) = STATES(10)*STATES(11)
ALGBRC(50) = (( CONSTS(109)*CONSTS(16))/( CONSTS(111)*CONSTS(110)))*log(CONSTS(42)/STATES(4))
ALGBRC(55) = CONSTS(10)*ALGBRC(22)*(STATES(1) - ALGBRC(50))
ALGBRC(59) = (( CONSTS(114)*CONSTS(16))/( CONSTS(116)*CONSTS(115)))*log(CONSTS(42)/STATES(4))
ALGBRC(62) = CONSTS(11)*STATES(6)*(STATES(1) - ALGBRC(59))
ALGBRC(77) = (( CONSTS(129)*CONSTS(16))/( CONSTS(131)*CONSTS(130)))*log(CONSTS(42)/STATES(4))
ALGBRC(80) = CONSTS(78)*CONSTS(95)*(STATES(1) - ALGBRC(77))
ALGBRC(82) = STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)+STATES(31)+STATES(32)+STATES(33)+STATES(34)+STATES(6)
ALGBRC(7) = 1.00000*EXP(( 8.47188*STATES(1))/( 1.00000*CONSTS(16)))
ALGBRC(33) = 0.0443240*ALGBRC(7)
ALGBRC(12) = 1.00000*EXP(( - 7.77556*STATES(1))/( 1.00000*CONSTS(16)))
ALGBRC(51) = 0.000106000*ALGBRC(12)
ALGBRC(66) = 1.00000*CONSTS(57)*STATES(2)
ALGBRC(28) = 0.000781000*ALGBRC(7)
ALGBRC(47) = 0.0953120*ALGBRC(12)
ALGBRC(63) = 2.00000*CONSTS(57)*STATES(2)
ALGBRC(24) = 2.81000e-05*ALGBRC(7)
ALGBRC(43) = 32.6594*ALGBRC(12)
ALGBRC(60) = 3.00000*CONSTS(57)*STATES(2)
ALGBRC(20) = 0.000869000*ALGBRC(7)
ALGBRC(39) = 144.174*ALGBRC(12)
ALGBRC(56) = 4.00000*CONSTS(57)*STATES(2)
ALGBRC(16) = 0.0216200*ALGBRC(7)
ALGBRC(35) = 318.108*ALGBRC(12)
ALGBRC(79) = 1.00000*CONSTS(57)*STATES(2)
ALGBRC(76) = 2.00000*CONSTS(57)*STATES(2)
ALGBRC(73) = 3.00000*CONSTS(57)*STATES(2)
ALGBRC(70) = 4.00000*CONSTS(57)*STATES(2)
ALGBRC(83) = ALGBRC(38)+ALGBRC(46)+ALGBRC(55)+ALGBRC(62)+ALGBRC(68)+ALGBRC(26)+ALGBRC(30)+ALGBRC(74)+ALGBRC(80)
ALGBRC(84) = STATES(16)+STATES(15)+STATES(14)+STATES(13)+STATES(20)+STATES(19)+STATES(18)+STATES(17)+STATES(5)+STATES(21)+STATES(22)+STATES(23)+STATES(25)+STATES(24)
ALGBRC(2) = CONSTS(134)*0.731000*EXP(STATES(1)/30.0000)
ALGBRC(23) = 2.00000*ALGBRC(2)
ALGBRC(27) = 1.00000*ALGBRC(2)
ALGBRC(10) = CONSTS(134)*0.214900*EXP(- STATES(1)/40.0000)
ALGBRC(41) = 3.00000*ALGBRC(10)
ALGBRC(45) = 4.00000*ALGBRC(10)
ALGBRC(81) = ( CONSTS(134)*4.00000)/(1.00000+1.00000/STATES(2))
ALGBRC(54) = CONSTS(134)*0.00595600*EXP(- STATES(1)/40.0000)
ALGBRC(49) = CONSTS(134)*0.474200*EXP(STATES(1)/10.0000)
ALGBRC(67) = CONSTS(134)*0.00232000*EXP(- STATES(1)/280.000)
ALGBRC(61) = CONSTS(134)*0.0121300*EXP(STATES(1)/500.000)
ALGBRC(72) = ( ALGBRC(45)*ALGBRC(61)*ALGBRC(54))/( ALGBRC(27)*ALGBRC(67))
ALGBRC(64) = CONSTS(134)*0.0219700*EXP(STATES(1)/500.000)
ALGBRC(58) = CONSTS(134)*0.0140700*EXP(- STATES(1)/300.000)
ALGBRC(69) = ( ALGBRC(45)*ALGBRC(58)*ALGBRC(49))/( ALGBRC(27)*ALGBRC(64))
ALGBRC(19) = 3.00000*ALGBRC(2)
ALGBRC(37) = 2.00000*ALGBRC(10)
ALGBRC(15) = 4.00000*ALGBRC(2)
ALGBRC(32) = 1.00000*ALGBRC(10)
ALGBRC(78) = ALGBRC(67)
ALGBRC(75) = ( ALGBRC(61)*ALGBRC(64))/ALGBRC(58)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END