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