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