C
C There are a total of 90 entries in the algebraic variable array.
C There are a total of 31 entries in each of the rate and state variable arrays.
C There are a total of 133 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is isotonic in component parameters (dimensionless).
C CONSTS(2) is alpha_1 in component parameters (per_micrometre).
C CONSTS(3) is beta_1 in component parameters (millinewton).
C CONSTS(4) is alpha_2 in component parameters (per_micrometre).
C CONSTS(5) is beta_2 in component parameters (millinewton).
C CONSTS(6) is alpha_3 in component parameters (per_micrometre).
C CONSTS(7) is beta_3 in component parameters (millinewton).
C CONSTS(8) is lambda in component parameters (millinewton).
C CONSTS(9) is mu in component parameters (dimensionless).
C CONSTS(10) is k_mu in component parameters (dimensionless).
C CONSTS(11) is kappa in component parameters (dimensionless).
C CONSTS(12) is kappa_0 in component parameters (dimensionless).
C CONSTS(13) is m_0 in component parameters (dimensionless).
C CONSTS(14) is v_max in component parameters (micrometre_per_second).
C CONSTS(15) is a in component parameters (dimensionless).
C CONSTS(16) is d_h in component parameters (dimensionless).
C CONSTS(17) is alpha_P in component parameters (dimensionless).
C ALGBRC(38) is l in component length (micrometre).
C ALGBRC(27) is F_muscle in component force (millinewton).
C CONSTS(18) is F_afterload in component isotonic (millinewton).
C ALGBRC(41) is isotonic_mode in component isotonic (dimensionless).
C CONSTS(19) is l_0 in component length (micrometre).
C CONSTS(20) is stim_period in component membrane (second).
C CONSTS(21) is S_0 in component parameters_izakov_et_al_1991 (micrometre).
C ALGBRC(3) is q_v in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(22) is q_1 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(23) is q_2 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(24) is q_3 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(25) is q_4 in component parameters_izakov_et_al_1991 (per_second).
C CONSTS(26) is x_st in component parameters_izakov_et_al_1991 (dimensionless).
C CONSTS(124) is v_st in component parameters_izakov_et_al_1991 (micrometre_per_second).
C CONSTS(125) is v_1 in component parameters_izakov_et_al_1991 (micrometre_per_second).
C CONSTS(27) is alpha_G in component parameters_izakov_et_al_1991 (dimensionless).
C CONSTS(28) is k_A in component parameters_izakov_et_al_1991 (per_millimolar).
C STATES(1) is v in component CE_velocity (micrometre_per_second).
C CONSTS(29) is alpha_Q in component parameters_izakov_et_al_1991 (dimensionless).
C CONSTS(30) is beta_Q in component parameters_izakov_et_al_1991 (dimensionless).
C ALGBRC(72) is F_CE in component force (millinewton).
C ALGBRC(16) is F_XSE in component force (millinewton).
C ALGBRC(1) is F_SE in component force (millinewton).
C ALGBRC(2) is F_PE in component force (millinewton).
C STATES(2) is N in component crossbridge_kinetics (dimensionless).
C ALGBRC(45) is k_P_vis in component CE_velocity (millinewton_second_per_micrometre).
C ALGBRC(52) is k_S_vis in component PE_velocity (millinewton_second_per_micrometre).
C STATES(3) is w in component PE_velocity (micrometre_per_second).
C STATES(4) is l_1 in component length (micrometre).
C STATES(5) is l_2 in component length (micrometre).
C STATES(6) is l_3 in component length (micrometre).
C ALGBRC(70) is p_v in component average_crossbridge_force (dimensionless).
C ALGBRC(64) is K_kappa in component crossbridge_kinetics (per_second).
C ALGBRC(32) is M_A in component crossbridge_kinetics (dimensionless).
C ALGBRC(34) is n_1 in component crossbridge_kinetics (dimensionless).
C ALGBRC(36) is L_oz in component crossbridge_kinetics (dimensionless).
C ALGBRC(60) is k_p_v in component crossbridge_kinetics (per_second).
C ALGBRC(62) is k_m_v in component crossbridge_kinetics (per_second).
C STATES(7) is A in component intracellular_calcium_concentration (millimolar).
C CONSTS(31) is A_tot in component intracellular_calcium_concentration (millimolar).
C ALGBRC(58) is G_star in component average_crossbridge_force (dimensionless).
C ALGBRC(54) is P_star in component average_crossbridge_force (dimensionless).
C CONSTS(32) is g_1 in component crossbridge_kinetics (per_micrometre).
C CONSTS(33) is g_2 in component crossbridge_kinetics (dimensionless).
C ALGBRC(4) is dl_1_dt in component length (micrometre_per_second).
C ALGBRC(55) is dl_2_dt in component length (micrometre_per_second).
C ALGBRC(56) is dl_3_dt in component length (micrometre_per_second).
C ALGBRC(48) is phi_chi_2 in component CE_velocity (micrometre_per_second).
C CONSTS(34) is stim_start in component membrane (second).
C CONSTS(35) is stim_end in component membrane (second).
C CONSTS(36) is stim_duration in component membrane (second).
C CONSTS(37) is stim_amplitude in component membrane (nanoA).
C ALGBRC(75) is phi_chi in component CE_velocity (micrometre_per_second2).
C ALGBRC(73) is p_prime_v in component average_crossbridge_force (second_per_micrometre).
C CONSTS(38) is alpha_P_lengthening in component CE_velocity (per_micrometre).
C CONSTS(39) is beta_P_lengthening in component CE_velocity (millinewton_second_per_micrometre).
C CONSTS(40) is alpha_P_shortening in component CE_velocity (per_micrometre).
C CONSTS(41) is beta_P_shortening in component CE_velocity (millinewton_second_per_micrometre).
C ALGBRC(43) is alp_p in component CE_velocity (per_micrometre).
C CONSTS(42) is alpha_S_lengthening in component PE_velocity (per_micrometre).
C CONSTS(43) is beta_S_lengthening in component PE_velocity (millinewton_second_per_micrometre).
C CONSTS(44) is alpha_S_shortening in component PE_velocity (per_micrometre).
C CONSTS(45) is beta_S_shortening in component PE_velocity (millinewton_second_per_micrometre).
C ALGBRC(50) is alp_s in component PE_velocity (per_micrometre).
C CONSTS(131) is gamma in component average_crossbridge_force (dimensionless).
C CONSTS(126) is case_1 in component average_crossbridge_force (second_per_micrometre).
C ALGBRC(66) is case_2 in component average_crossbridge_force (second_per_micrometre).
C CONSTS(127) is case_3 in component average_crossbridge_force (second_per_micrometre).
C ALGBRC(68) is case_4 in component average_crossbridge_force (second_per_micrometre).
C STATES(8) is V in component membrane (millivolt).
C CONSTS(46) is R in component membrane (joule_per_kilomole_kelvin).
C CONSTS(47) is T in component membrane (kelvin).
C CONSTS(48) is F in component membrane (coulomb_per_mole).
C CONSTS(49) is Cm in component membrane (microF).
C ALGBRC(53) is i_K1 in component time_independent_potassium_current (nanoA).
C ALGBRC(80) is i_to in component transient_outward_current (nanoA).
C ALGBRC(57) is i_Kr in component rapid_delayed_rectifier_potassium_current (nanoA).
C ALGBRC(59) is i_Ks in component slow_delayed_rectifier_potassium_current (nanoA).
C ALGBRC(69) is i_Ca_L_K_cyt in component L_type_Ca_channel (nanoA).
C ALGBRC(76) is i_Ca_L_K_ds in component L_type_Ca_channel (nanoA).
C ALGBRC(81) is i_NaK in component sodium_potassium_pump (nanoA).
C ALGBRC(61) is i_Na in component fast_sodium_current (nanoA).
C ALGBRC(65) is i_b_Na in component sodium_background_current (nanoA).
C ALGBRC(63) is i_p_Na in component persistent_sodium_current (nanoA).
C ALGBRC(71) is i_Ca_L_Na_cyt in component L_type_Ca_channel (nanoA).
C ALGBRC(77) is i_Ca_L_Na_ds in component L_type_Ca_channel (nanoA).
C ALGBRC(82) is i_NaCa_cyt in component sodium_calcium_exchanger (nanoA).
C ALGBRC(83) is i_NaCa_ds in component sodium_calcium_exchanger (nanoA).
C ALGBRC(67) is i_Ca_L_Ca_cyt in component L_type_Ca_channel (nanoA).
C ALGBRC(74) is i_Ca_L_Ca_ds in component L_type_Ca_channel (nanoA).
C ALGBRC(79) is i_b_Ca in component calcium_background_current (nanoA).
C ALGBRC(40) is i_Stim in component membrane (nanoA).
C ALGBRC(42) is E_Na in component reversal_potentials (millivolt).
C ALGBRC(44) is E_K in component reversal_potentials (millivolt).
C ALGBRC(47) is E_Ks in component reversal_potentials (millivolt).
C ALGBRC(49) is E_Ca in component reversal_potentials (millivolt).
C ALGBRC(51) is E_mh in component reversal_potentials (millivolt).
C CONSTS(50) is P_kna in component reversal_potentials (dimensionless).
C STATES(9) is K_o in component extracellular_potassium_concentration (millimolar).
C CONSTS(51) is Na_o in component extracellular_sodium_concentration (millimolar).
C STATES(10) is K_i in component intracellular_potassium_concentration (millimolar).
C STATES(11) is Na_i in component intracellular_sodium_concentration (millimolar).
C CONSTS(52) is Ca_o in component extracellular_calcium_concentration (millimolar).
C STATES(12) is Ca_i in component intracellular_calcium_concentration (millimolar).
C CONSTS(53) is K_mk1 in component time_independent_potassium_current (millimolar).
C CONSTS(54) is g_K1 in component time_independent_potassium_current (microS).
C CONSTS(55) is g_Kr1 in component rapid_delayed_rectifier_potassium_current (microS).
C CONSTS(56) is g_Kr2 in component rapid_delayed_rectifier_potassium_current (microS).
C STATES(13) is xr1 in component rapid_delayed_rectifier_potassium_current_xr1_gate (dimensionless).
C STATES(14) is xr2 in component rapid_delayed_rectifier_potassium_current_xr2_gate (dimensionless).
C ALGBRC(5) is alpha_xr1 in component rapid_delayed_rectifier_potassium_current_xr1_gate (per_second).
C ALGBRC(17) is beta_xr1 in component rapid_delayed_rectifier_potassium_current_xr1_gate (per_second).
C ALGBRC(6) is alpha_xr2 in component rapid_delayed_rectifier_potassium_current_xr2_gate (per_second).
C ALGBRC(18) is beta_xr2 in component rapid_delayed_rectifier_potassium_current_xr2_gate (per_second).
C CONSTS(57) is g_Ks in component slow_delayed_rectifier_potassium_current (microS).
C STATES(15) is xs in component slow_delayed_rectifier_potassium_current_xs_gate (dimensionless).
C ALGBRC(7) is alpha_xs in component slow_delayed_rectifier_potassium_current_xs_gate (per_second).
C ALGBRC(19) is beta_xs in component slow_delayed_rectifier_potassium_current_xs_gate (per_second).
C ALGBRC(46) is i_KNa in component sodium_activated_potassium_current (nanoA).
C CONSTS(58) is g_K_Na in component sodium_activated_potassium_current (microS).
C CONSTS(59) is K_kna in component sodium_activated_potassium_current (millimolar).
C CONSTS(60) is g_Na in component fast_sodium_current (microS).
C STATES(16) is m in component fast_sodium_current_m_gate (dimensionless).
C STATES(17) is h in component fast_sodium_current_h_gate (dimensionless).
C ALGBRC(20) is alpha_m in component fast_sodium_current_m_gate (per_second).
C ALGBRC(28) is beta_m in component fast_sodium_current_m_gate (per_second).
C CONSTS(61) is delta_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(8) is E0_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(9) is alpha_h in component fast_sodium_current_h_gate (per_second).
C ALGBRC(21) is beta_h in component fast_sodium_current_h_gate (per_second).
C CONSTS(62) is g_pna in component persistent_sodium_current (microS).
C CONSTS(63) is g_bna in component sodium_background_current (microS).
C ALGBRC(78) is i_Ca_L in component L_type_Ca_channel (nanoA).
C CONSTS(64) is P_Ca_L in component L_type_Ca_channel (nanoA_per_millimolar).
C CONSTS(65) is P_CaK in component L_type_Ca_channel (dimensionless).
C CONSTS(66) is P_CaNa in component L_type_Ca_channel (dimensionless).
C STATES(18) is Ca_ds in component intracellular_calcium_concentration (millimolar).
C STATES(19) is d in component L_type_Ca_channel_d_gate (dimensionless).
C STATES(20) is f in component L_type_Ca_channel_f_gate (dimensionless).
C STATES(21) is f2 in component L_type_Ca_channel_f2_gate (dimensionless).
C STATES(22) is f2ds in component L_type_Ca_channel_f2ds_gate (dimensionless).
C CONSTS(67) is Km_f2 in component L_type_Ca_channel (millimolar).
C CONSTS(68) is Km_f2ds in component L_type_Ca_channel (millimolar).
C CONSTS(69) is R_decay in component L_type_Ca_channel (per_second).
C CONSTS(70) is FrICa in component L_type_Ca_channel (dimensionless).
C ALGBRC(22) is alpha_d in component L_type_Ca_channel_d_gate (per_second).
C ALGBRC(29) is beta_d in component L_type_Ca_channel_d_gate (per_second).
C ALGBRC(10) is E0_d in component L_type_Ca_channel_d_gate (millivolt).
C CONSTS(71) is speed_d in component L_type_Ca_channel_d_gate (dimensionless).
C ALGBRC(23) is alpha_f in component L_type_Ca_channel_f_gate (per_second).
C ALGBRC(30) is beta_f in component L_type_Ca_channel_f_gate (per_second).
C CONSTS(72) is speed_f in component L_type_Ca_channel_f_gate (dimensionless).
C CONSTS(73) is delta_f in component L_type_Ca_channel_f_gate (millivolt).
C ALGBRC(11) is E0_f in component L_type_Ca_channel_f_gate (millivolt).
C CONSTS(74) is g_bca in component calcium_background_current (microS).
C CONSTS(75) is g_to in component transient_outward_current (microS).
C CONSTS(76) is g_tos in component transient_outward_current (dimensionless).
C STATES(23) is s in component transient_outward_current_s_gate (dimensionless).
C STATES(24) is r in component transient_outward_current_r_gate (dimensionless).
C ALGBRC(12) is alpha_s in component transient_outward_current_s_gate (per_second).
C ALGBRC(24) is beta_s in component transient_outward_current_s_gate (per_second).
C CONSTS(77) is i_NaK_max in component sodium_potassium_pump (nanoA).
C CONSTS(78) is K_mK in component sodium_potassium_pump (millimolar).
C CONSTS(79) is K_mNa in component sodium_potassium_pump (millimolar).
C ALGBRC(85) is i_NaCa in component sodium_calcium_exchanger (nanoA).
C CONSTS(80) is k_NaCa in component sodium_calcium_exchanger (nanoA).
C CONSTS(81) is n_NaCa in component sodium_calcium_exchanger (dimensionless).
C CONSTS(82) is d_NaCa in component sodium_calcium_exchanger (dimensionless).
C CONSTS(83) is gamma in component sodium_calcium_exchanger (dimensionless).
C CONSTS(84) is FRiNaCa in component sodium_calcium_exchanger (dimensionless).
C ALGBRC(86) is i_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C CONSTS(129) is K_1 in component sarcoplasmic_reticulum_calcium_pump (dimensionless).
C ALGBRC(84) is K_2 in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(85) is K_cyca in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(86) is K_xcs in component sarcoplasmic_reticulum_calcium_pump (dimensionless).
C CONSTS(87) is K_srca in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C CONSTS(88) is alpha_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C CONSTS(89) is beta_up in component sarcoplasmic_reticulum_calcium_pump (millimolar_per_second).
C STATES(25) is Ca_up in component intracellular_calcium_concentration (millimolar).
C CONSTS(90) is flag_ingib in component sarcoplasmic_reticulum_calcium_pump (dimensionless).
C CONSTS(91) is K_inh in component sarcoplasmic_reticulum_calcium_pump (millimolar).
C ALGBRC(87) is i_trans in component calcium_translocation (millimolar_per_second).
C STATES(26) is Ca_rel in component intracellular_calcium_concentration (millimolar).
C CONSTS(92) is a_tr in component calcium_translocation (per_second).
C CONSTS(93) is alpha_CaS in component calcium_translocation (per_millimolar_second).
C CONSTS(94) is beta_CaS in component calcium_translocation (per_second).
C CONSTS(95) is CaS_tot in component calcium_translocation (millimolar).
C CONSTS(130) is beta in component calcium_translocation (millimolar).
C ALGBRC(88) is i_rel in component calcium_release (millimolar_per_second).
C ALGBRC(13) is VoltDep in component calcium_release (dimensionless).
C ALGBRC(31) is RegBindSite in component calcium_release (dimensionless).
C ALGBRC(14) is CaiReg in component calcium_release (dimensionless).
C ALGBRC(25) is CadsReg in component calcium_release (dimensionless).
C ALGBRC(33) is ActRate in component calcium_release (per_second).
C ALGBRC(35) is InactRate in component calcium_release (per_second).
C CONSTS(96) is SRLeak in component calcium_release (per_second).
C CONSTS(97) is K_m_rel in component calcium_release (per_second).
C CONSTS(98) is K_m_Ca_cyt in component calcium_release (millimolar).
C CONSTS(99) is K_m_Ca_ds in component calcium_release (millimolar).
C ALGBRC(39) is PrecFrac in component calcium_release (dimensionless).
C STATES(27) is ActFrac in component calcium_release (dimensionless).
C STATES(28) is ProdFrac in component calcium_release (dimensionless).
C ALGBRC(37) is SpeedRel in component calcium_release (dimensionless).
C CONSTS(132) is V_i in component intracellular_calcium_concentration (micrometre3).
C CONSTS(100) is n_NaK in component intracellular_sodium_concentration (dimensionless).
C CONSTS(101) is K_b in component extracellular_potassium_concentration (millimolar).
C CONSTS(133) is V_e in component intracellular_calcium_concentration (micrometre3).
C CONSTS(102) is radius in component intracellular_calcium_concentration (micrometre).
C CONSTS(103) is length in component intracellular_calcium_concentration (micrometre).
C CONSTS(128) is V_Cell in component intracellular_calcium_concentration (micrometre3).
C CONSTS(104) is V_i_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(105) is V_ds_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(106) is V_rel_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(107) is V_e_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(108) is V_up_ratio in component intracellular_calcium_concentration (dimensionless).
C CONSTS(109) is Kdecay in component intracellular_calcium_concentration (per_second).
C ALGBRC(89) is N_A in component intracellular_calcium_concentration (dimensionless).
C ALGBRC(90) is pi_N_A in component intracellular_calcium_concentration (dimensionless).
C CONSTS(110) is pi_min in component intracellular_calcium_concentration (dimensionless).
C CONSTS(111) is a_on in component intracellular_calcium_concentration (per_millimolar_second).
C CONSTS(112) is a_off in component intracellular_calcium_concentration (per_second).
C STATES(29) is B_1 in component intracellular_calcium_concentration (millimolar).
C CONSTS(113) is B_1_tot in component intracellular_calcium_concentration (millimolar).
C CONSTS(114) is b_1_on in component intracellular_calcium_concentration (per_millimolar_second).
C CONSTS(115) is b_1_off in component intracellular_calcium_concentration (per_second).
C STATES(30) is B_2 in component intracellular_calcium_concentration (millimolar).
C CONSTS(116) is B_2_tot in component intracellular_calcium_concentration (millimolar).
C CONSTS(117) is b_2_on in component intracellular_calcium_concentration (per_millimolar_second).
C CONSTS(118) is b_2_off in component intracellular_calcium_concentration (per_second).
C CONSTS(119) is g_fibro_junct in component fibroblast (microS).
C CONSTS(120) is g_fibro in component fibroblast (microS).
C CONSTS(121) is c_fibro in component fibroblast (microF).
C CONSTS(122) is g_fibro_stretch in component fibroblast (microS).
C CONSTS(123) is E_fibro_stretch in component fibroblast (millivolt).
C STATES(31) is V_fibro in component fibroblast (millivolt).
C ALGBRC(15) is i_fibro in component fibroblast (nanoA).
C ALGBRC(26) is i_fibro_junct in component fibroblast (nanoA).
C RATES(2) is d/dt N in component crossbridge_kinetics (dimensionless).
C RATES(4) is d/dt l_1 in component length (micrometre).
C RATES(5) is d/dt l_2 in component length (micrometre).
C RATES(6) is d/dt l_3 in component length (micrometre).
C RATES(1) is d/dt v in component CE_velocity (micrometre_per_second).
C RATES(3) is d/dt w in component PE_velocity (micrometre_per_second).
C RATES(8) is d/dt V in component membrane (millivolt).
C RATES(13) is d/dt xr1 in component rapid_delayed_rectifier_potassium_current_xr1_gate (dimensionless).
C RATES(14) is d/dt xr2 in component rapid_delayed_rectifier_potassium_current_xr2_gate (dimensionless).
C RATES(15) is d/dt xs in component slow_delayed_rectifier_potassium_current_xs_gate (dimensionless).
C RATES(16) is d/dt m in component fast_sodium_current_m_gate (dimensionless).
C RATES(17) is d/dt h in component fast_sodium_current_h_gate (dimensionless).
C RATES(19) is d/dt d in component L_type_Ca_channel_d_gate (dimensionless).
C RATES(20) is d/dt f in component L_type_Ca_channel_f_gate (dimensionless).
C RATES(21) is d/dt f2 in component L_type_Ca_channel_f2_gate (dimensionless).
C RATES(22) is d/dt f2ds in component L_type_Ca_channel_f2ds_gate (dimensionless).
C RATES(23) is d/dt s in component transient_outward_current_s_gate (dimensionless).
C RATES(24) is d/dt r in component transient_outward_current_r_gate (dimensionless).
C RATES(27) is d/dt ActFrac in component calcium_release (dimensionless).
C RATES(28) is d/dt ProdFrac in component calcium_release (dimensionless).
C RATES(11) is d/dt Na_i in component intracellular_sodium_concentration (millimolar).
C RATES(9) is d/dt K_o in component extracellular_potassium_concentration (millimolar).
C RATES(10) is d/dt K_i in component intracellular_potassium_concentration (millimolar).
C RATES(7) is d/dt A in component intracellular_calcium_concentration (millimolar).
C RATES(29) is d/dt B_1 in component intracellular_calcium_concentration (millimolar).
C RATES(30) is d/dt B_2 in component intracellular_calcium_concentration (millimolar).
C RATES(12) is d/dt Ca_i in component intracellular_calcium_concentration (millimolar).
C RATES(18) is d/dt Ca_ds in component intracellular_calcium_concentration (millimolar).
C RATES(25) is d/dt Ca_up in component intracellular_calcium_concentration (millimolar).
C RATES(26) is d/dt Ca_rel in component intracellular_calcium_concentration (millimolar).
C RATES(31) is d/dt V_fibro in component fibroblast (millivolt).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0
      CONSTS(2) = 14.6
      CONSTS(3) = 0.84
      CONSTS(4) = 14.6
      CONSTS(5) = 0.0018
      CONSTS(6) = 48
      CONSTS(7) = 0.015
      CONSTS(8) = 30
      CONSTS(9) = 3
      CONSTS(10) = 0.6
      CONSTS(11) = 0.705
      CONSTS(12) = 3
      CONSTS(13) = 0.9
      CONSTS(14) = 5.5
      CONSTS(15) = 0.25
      CONSTS(16) = 0.5
      CONSTS(17) = 4
      CONSTS(18) = 2
      CONSTS(19) = 0.525139356105856
      CONSTS(20) = 1
      CONSTS(21) = 1.14
      CONSTS(22) = 17.3
      CONSTS(23) = 259
      CONSTS(24) = 17.3
      CONSTS(25) = 15
      CONSTS(26) = 0.964285
      CONSTS(27) = 1
      CONSTS(28) = 40
      STATES(1) = 0
      CONSTS(29) = 10
      CONSTS(30) = 5
      STATES(2) = 7.455e-8
      STATES(3) = 0
      STATES(4) = 0.436333342969918
      STATES(5) = 0.436333525334166
      STATES(6) = 0.088805830771694
      STATES(7) = 0.00015
      CONSTS(31) = 0.07
      CONSTS(32) = 0.6
      CONSTS(33) = 0.52
      CONSTS(34) = 0.06
      CONSTS(35) = 10000
      CONSTS(36) = 0.0025
      CONSTS(37) = -3
      CONSTS(38) = 16
      CONSTS(39) = 0.0015
      CONSTS(40) = 16
      CONSTS(41) = 0.0015
      CONSTS(42) = 46
      CONSTS(43) = 0
      CONSTS(44) = 39
      CONSTS(45) = 0
      STATES(8) = -93.658148
      CONSTS(46) = 8314.472
      CONSTS(47) = 310
      CONSTS(48) = 96485.3415
      CONSTS(49) = 9.5e-5
      CONSTS(50) = 0.03
      STATES(9) = 3.988
      CONSTS(51) = 140
      STATES(10) = 139.054
      STATES(11) = 5.18787513289509
      CONSTS(52) = 2
      STATES(12) = 6.15e-6
      CONSTS(53) = 10
      CONSTS(54) = 0.5
      CONSTS(55) = 0.0021
      CONSTS(56) = 0.0013
      STATES(13) = 8.88859784542779e-6
      STATES(14) = 1.53745791069154e-7
      CONSTS(57) = 0.0026
      STATES(15) = 0.001
      CONSTS(58) = 0
      CONSTS(59) = 20
      CONSTS(60) = 2.5
      STATES(16) = 0.0015
      STATES(17) = 0.995
      CONSTS(61) = 1e-5
      CONSTS(62) = 0.004
      CONSTS(63) = 0.0006
      CONSTS(64) = 0.1
      CONSTS(65) = 0.002
      CONSTS(66) = 0.01
      STATES(18) = 2.55e-6
      STATES(19) = 0
      STATES(20) = 1
      STATES(21) = 1
      STATES(22) = 0.997
      CONSTS(67) = 100000
      CONSTS(68) = 0.001
      CONSTS(69) = 20
      CONSTS(70) = 1
      CONSTS(71) = 3
      CONSTS(72) = 0.3
      CONSTS(73) = 0.0001
      CONSTS(74) = 0.00025
      CONSTS(75) = 0.006
      CONSTS(76) = 0
      STATES(23) = 0.997044616031121
      STATES(24) = 1.63117173173398e-8
      CONSTS(77) = 0.7
      CONSTS(78) = 1
      CONSTS(79) = 24.2
      CONSTS(80) = 0.0005
      CONSTS(81) = 3
      CONSTS(82) = 0
      CONSTS(83) = 0.5
      CONSTS(84) = 0.001
      CONSTS(85) = 0.00015
      CONSTS(86) = 0.4
      CONSTS(87) = 0.5
      CONSTS(88) = 1
      CONSTS(89) = 0.03
      STATES(25) = 0.994579
      CONSTS(90) = 0
      CONSTS(91) = 4
      STATES(26) = 0.989665
      CONSTS(92) = 15
      CONSTS(93) = 50000
      CONSTS(94) = 32500
      CONSTS(95) = 40
      CONSTS(96) = 0.05
      CONSTS(97) = 10000
      CONSTS(98) = 0.0005
      CONSTS(99) = 0.01
      STATES(27) = 0.001914
      STATES(28) = 0.2854569
      CONSTS(100) = 1.5
      CONSTS(101) = 4
      CONSTS(102) = 12
      CONSTS(103) = 74
      CONSTS(104) = 0.49
      CONSTS(105) = 0.1
      CONSTS(106) = 0.003
      CONSTS(107) = 0.4
      CONSTS(108) = 0.03
      CONSTS(109) = 10
      CONSTS(110) = 0.03
      CONSTS(111) = 70000
      CONSTS(112) = 200
      STATES(29) = 0
      CONSTS(113) = 0.08
      CONSTS(114) = 100000
      CONSTS(115) = 182
      STATES(30) = 0
      CONSTS(116) = 0.1
      CONSTS(117) = 1000
      CONSTS(118) = 3
      CONSTS(119) = 2.9e-4
      CONSTS(120) = 2e-4
      CONSTS(121) = 1e-5
      CONSTS(122) = 0
      CONSTS(123) = 0
      STATES(31) = -20
      CONSTS(124) =  CONSTS(26)*CONSTS(14)
      CONSTS(125) = CONSTS(14)/10.0000
      CONSTS(126) = ( CONSTS(15)*(0.400000+ 0.400000*CONSTS(15)))/( CONSTS(14)* (CONSTS(15)+1.00000)*0.400000 ** 2.00000)
      CONSTS(127) = ( 0.400000*CONSTS(15)+1.00000)/( CONSTS(15)*CONSTS(14))
      CONSTS(128) = ( 3.14159*CONSTS(102)/1000.00 ** 2.00000*CONSTS(103))/1000.00
      CONSTS(129) = ( CONSTS(85)*CONSTS(86))/CONSTS(87)
      CONSTS(130) = CONSTS(94)/CONSTS(93)
      CONSTS(131) = ( CONSTS(15)*CONSTS(16)*CONSTS(125)/CONSTS(14) ** 2.00000)/( 3.00000*CONSTS(15)*CONSTS(16) - ( (CONSTS(15)+1.00000)*CONSTS(125))/CONSTS(14))
      CONSTS(132) =  CONSTS(128)*CONSTS(104)
      CONSTS(133) =  CONSTS(128)*CONSTS(107)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RATES(21) = 1.00000 -  1.00000*(STATES(12)/(CONSTS(67)+STATES(12))+STATES(21))
      RATES(22) =  CONSTS(69)*(1.00000 - (STATES(18)/(CONSTS(68)+STATES(18))+STATES(22)))
      RATES(24) =  333.000*(1.00000/(1.00000+EXP(- (STATES(8)+4.00000)/5.00000)) - STATES(24))
      ALGBRC(4) = STATES(1)
      RATES(4) = ALGBRC(4)
      ALGBRC(5) = 50.0000/(1.00000+EXP(- (STATES(8) - 5.00000)/9.00000))
      ALGBRC(17) =  0.0500000*EXP(- (STATES(8) - 20.0000)/15.0000)
      RATES(13) =  ALGBRC(5)*(1.00000 - STATES(13)) -  ALGBRC(17)*STATES(13)
      ALGBRC(6) = 50.0000/(1.00000+EXP(- (STATES(8) - 5.00000)/9.00000))
      ALGBRC(18) =  0.400000*EXP(- (STATES(8)+30.0000)/30.0000 ** 3.00000)
      RATES(14) =  ALGBRC(6)*(1.00000 - STATES(14)) -  ALGBRC(18)*STATES(14)
      ALGBRC(7) = 14.0000/(1.00000+EXP(- (STATES(8) - 40.0000)/9.00000))
      ALGBRC(19) =  1.00000*EXP(- STATES(8)/45.0000)
      RATES(15) =  ALGBRC(7)*(1.00000 - STATES(15)) -  ALGBRC(19)*STATES(15)
      ALGBRC(9) =  20.0000*EXP( - 0.125000*(STATES(8)+75.0000))
      ALGBRC(21) = 2000.00/(1.00000+ 320.000*EXP( - 0.100000*(STATES(8)+75.0000)))
      RATES(17) =  ALGBRC(9)*(1.00000 - STATES(17)) -  ALGBRC(21)*STATES(17)
      ALGBRC(12) =  0.0330000*EXP(- STATES(8)/17.0000)
      ALGBRC(24) = 33.0000/(1.00000+EXP( - 0.125000*(STATES(8)+10.0000)))
      RATES(23) =  ALGBRC(12)*(1.00000 - STATES(23)) -  ALGBRC(24)*STATES(23)
      ALGBRC(15) =  CONSTS(120)*(STATES(31)+20.0000)+ CONSTS(122)*(STATES(31) - CONSTS(123))
      ALGBRC(26) =  - CONSTS(119)*(STATES(8) - STATES(31))
      RATES(31) = - (ALGBRC(15)+ALGBRC(26))/CONSTS(121)
      ALGBRC(8) = STATES(8)+41.0000
      ALGBRC(20) = TERNRY(ABS(ALGBRC(8)).LT.CONSTS(61), 2000.00, ( 200.000*ALGBRC(8))/(1.00000 - EXP( - 0.100000*ALGBRC(8))))
      ALGBRC(28) =  8000.00*EXP( - 0.0560000*(STATES(8)+66.0000))
      RATES(16) =  ALGBRC(20)*(1.00000 - STATES(16)) -  ALGBRC(28)*STATES(16)
      ALGBRC(10) = (STATES(8)+24.0000) - 5.00000
      ALGBRC(22) = TERNRY(ABS(ALGBRC(10)).LT.0.000100000, 120.000, ( 30.0000*ALGBRC(10))/(1.00000 - EXP(- ALGBRC(10)/4.00000)))
      ALGBRC(29) = TERNRY(ABS(ALGBRC(10)).LT.0.000100000, 120.000, ( 12.0000*ALGBRC(10))/(EXP(ALGBRC(10)/10.0000) - 1.00000))
      RATES(19) =  CONSTS(71)*( ALGBRC(22)*(1.00000 - STATES(19)) -  ALGBRC(29)*STATES(19))
      ALGBRC(11) = STATES(8)+34.0000
      ALGBRC(23) = TERNRY(ABS(ALGBRC(11)).LT.CONSTS(73), 25.0000, ( 6.25000*ALGBRC(11))/(EXP(ALGBRC(11)/4.00000) - 1.00000))
      ALGBRC(30) = 12.0000/(1.00000+EXP(( - 1.00000*(STATES(8)+34.0000))/4.00000))
      RATES(20) =  CONSTS(72)*( ALGBRC(23)*(1.00000 - STATES(20)) -  ALGBRC(30)*STATES(20))
      ALGBRC(14) = STATES(12)/(STATES(12)+CONSTS(98))
      ALGBRC(25) = STATES(18)/(STATES(18)+CONSTS(99))
      ALGBRC(31) = ALGBRC(14)+ (1.00000 - ALGBRC(14))*ALGBRC(25)
      ALGBRC(35) = 60.0000+ 500.000*ALGBRC(31) ** 2.00000
      ALGBRC(37) = TERNRY(STATES(8).LT.- 50.0000, 5.00000, 1.00000)
      RATES(28) =  STATES(27)*ALGBRC(37)*ALGBRC(35) -  ALGBRC(37)*0.600000*STATES(28)
      ALGBRC(33) =  500.000*ALGBRC(31) ** 2.00000
      ALGBRC(39) = (1.00000 - STATES(27)) - STATES(28)
      RATES(27) =  ALGBRC(39)*ALGBRC(37)*ALGBRC(33) -  STATES(27)*ALGBRC(37)*ALGBRC(35)
      ALGBRC(52) = TERNRY(STATES(3).LE.STATES(1),  CONSTS(43)*EXP( CONSTS(42)*(STATES(5) - STATES(4))),  CONSTS(45)*EXP( CONSTS(44)*(STATES(5) - STATES(4))))
      ALGBRC(38) = STATES(5)+STATES(6)
      ALGBRC(16) =  CONSTS(7)*(EXP( CONSTS(6)*STATES(6)) - 1.00000)
      ALGBRC(27) = ALGBRC(16)
      ALGBRC(41) = TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(27).GT.CONSTS(18).AND.ALGBRC(38).LE. CONSTS(19)*(1.00000+0.000100000), 1.00000, 0.00000)
      ALGBRC(48) = TERNRY(ALGBRC(41).EQ.1.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))), ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6))))
      ALGBRC(55) = TERNRY(ALGBRC(52).EQ.0.00000, ALGBRC(48), STATES(3))
      RATES(5) = ALGBRC(55)
      ALGBRC(56) = TERNRY(ALGBRC(41).EQ.1.00000, 0.00000, TERNRY(ALGBRC(41).EQ.0.00000.AND.ALGBRC(52).EQ.0.00000, - ALGBRC(48), TERNRY(ALGBRC(41).EQ.0.00000.AND.ALGBRC(52) .NE. 0.00000, - STATES(3), 0.0/0.0)
      RATES(6) = TERNRY(VOI.GE.CONSTS(34).AND.VOI.LE.CONSTS(35).AND.(VOI - CONSTS(34)) -  INT((VOI - CONSTS(34))/CONSTS(20))*CONSTS(20).LE.CONSTS(36), - ((STATES(5)+STATES(6)) - CONSTS(19))/CONSTS(36), ALGBRC(56))
      ALGBRC(32) = ( STATES(7)/CONSTS(31) ** CONSTS(9)*(1.00000+CONSTS(10) ** CONSTS(9)))/(STATES(7)/CONSTS(31) ** CONSTS(9)+CONSTS(10) ** CONSTS(9))
      ALGBRC(34) = TERNRY( CONSTS(32)*STATES(4)+CONSTS(33).LT.0.00000, 0.00000, TERNRY( CONSTS(32)*STATES(4)+CONSTS(33).LT.1.00000,  CONSTS(32)*STATES(4)+CONSTS(33), 1.00000)
      ALGBRC(36) = TERNRY(STATES(4).GT.0.550000, (STATES(4)+CONSTS(21))/(0.460000+CONSTS(21)),  (CONSTS(21)+0.550000)*1.00000)
      ALGBRC(3) = TERNRY(STATES(1).LE.0.00000, CONSTS(22) - ( CONSTS(23)*STATES(1))/CONSTS(14), TERNRY(STATES(1).LE.CONSTS(124).AND.0.00000.LT.STATES(1), ( (CONSTS(25) - CONSTS(24))*STATES(1))/CONSTS(124)+CONSTS(24), CONSTS(25)/1.00000+( CONSTS(30)*(STATES(1) - CONSTS(124)))/CONSTS(14) ** CONSTS(29))
      ALGBRC(54) = TERNRY(STATES(1).LE.0.00000, ( CONSTS(15)*(1.00000+STATES(1)/CONSTS(14)))/(CONSTS(15) - STATES(1)/CONSTS(14)), (1.00000+CONSTS(16)) - ( CONSTS(16) ** 2.00000*CONSTS(15))/( (( CONSTS(15)*CONSTS(16))/CONSTS(131))*STATES(1)/CONSTS(14) ** 2.00000+( (CONSTS(15)+1.00000)*STATES(1))/CONSTS(14)+ CONSTS(15)*CONSTS(16)))
      ALGBRC(58) = TERNRY(- CONSTS(14).LE.STATES(1).AND.STATES(1).LE.0.00000, 1.00000+( 0.600000*STATES(1))/CONSTS(14), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(125), ALGBRC(54)/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000), ( ALGBRC(54)*EXP( - CONSTS(27)*(STATES(1) - CONSTS(125))/CONSTS(14) ** CONSTS(17)))/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000))
      ALGBRC(60) =  CONSTS(11)*CONSTS(12)*ALGBRC(3)*CONSTS(13)*ALGBRC(58)
      ALGBRC(62) =  CONSTS(12)*ALGBRC(3)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(58))
      ALGBRC(64) =  ALGBRC(60)*ALGBRC(32)*ALGBRC(34)*ALGBRC(36)*(1.00000 - STATES(2)) -  ALGBRC(62)*STATES(2)
      RATES(2) = ALGBRC(64)
      ALGBRC(45) = TERNRY(STATES(1).LE.0.00000,  CONSTS(39)*EXP( CONSTS(38)*STATES(4)),  CONSTS(41)*EXP( CONSTS(40)*STATES(4)))
      ALGBRC(70) = ALGBRC(54)/ALGBRC(58)
      ALGBRC(66) = ( CONSTS(15)*1.00000*(1.00000+ 0.400000*CONSTS(15)+( 1.20000*STATES(1))/CONSTS(14)+ 0.600000*STATES(1)/CONSTS(14) ** 2.00000))/( CONSTS(14)* (CONSTS(15) - STATES(1)/CONSTS(14))*(1.00000+( 0.600000*STATES(1))/CONSTS(14)) ** 2.00000)
      ALGBRC(68) =  (1.00000/CONSTS(14))*EXP( - CONSTS(27)*STATES(1)/CONSTS(14) - CONSTS(125)/CONSTS(14) ** CONSTS(17))*(( 0.400000*CONSTS(15)+1.00000)/CONSTS(15)+ CONSTS(27)*CONSTS(17)*(1.00000+( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14)))*STATES(1)/CONSTS(14) - CONSTS(125)/CONSTS(14) ** CONSTS(17) - 1.00000)
      ALGBRC(73) = TERNRY(STATES(1).LE.- CONSTS(14), CONSTS(126), TERNRY(- CONSTS(14).LT.STATES(1).AND.STATES(1).LE.0.00000, ALGBRC(66), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(125), CONSTS(127), ALGBRC(68))
      ALGBRC(43) = TERNRY(STATES(1).LE.0.00000, CONSTS(38), CONSTS(40))
      ALGBRC(75) = TERNRY(ALGBRC(41).EQ.1.00000, - ( CONSTS(8)*ALGBRC(64)*ALGBRC(70)+ ALGBRC(43)*ALGBRC(45)*STATES(1) ** 2.00000+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(73)+ALGBRC(45)), - ( CONSTS(8)*ALGBRC(64)*ALGBRC(70)+ ALGBRC(43)*ALGBRC(45)*STATES(1) ** 2.00000+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(73)+ALGBRC(45)))
      RATES(1) = TERNRY(ALGBRC(52).EQ.0.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(ALGBRC(48) - STATES(1)) - ( CONSTS(8)*ALGBRC(64)*ALGBRC(70)+ ALGBRC(43)*ALGBRC(45)*STATES(1) ** 2.00000))/( CONSTS(8)*STATES(2)*ALGBRC(73)+ALGBRC(45)), ALGBRC(75))
      ALGBRC(50) = TERNRY(STATES(3).LE.STATES(1), CONSTS(42), CONSTS(44))
      RATES(3) = TERNRY(ALGBRC(41).EQ.1.00000.AND.ALGBRC(52) .NE. 0.00000, (( ALGBRC(52)*(ALGBRC(75) -  ALGBRC(50)*STATES(3) - STATES(1) ** 2.00000) -  CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(STATES(3) - STATES(1))) -  CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/ALGBRC(52), TERNRY(ALGBRC(41).EQ.0.00000.AND.ALGBRC(52) .NE. 0.00000, (ALGBRC(75) -  ALGBRC(50)*STATES(3) - STATES(1) ** 2.00000) - ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*(STATES(3) - STATES(1))+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/ALGBRC(52), TERNRY(ALGBRC(52).EQ.0.00000, 0.00000, 0.0/0.0)
      ALGBRC(44) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log(STATES(9)/STATES(10))
      ALGBRC(53) = ( (( CONSTS(54)*STATES(9))/(STATES(9)+CONSTS(53)))*(STATES(8) - ALGBRC(44)))/(1.00000+EXP(( ((STATES(8) - ALGBRC(44)) - 10.0000)*CONSTS(48)*1.25000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(80) =  CONSTS(75)*(CONSTS(76)+ STATES(23)*(1.00000 - CONSTS(76)))*STATES(24)*(STATES(8) - ALGBRC(44))
      ALGBRC(57) =  (( ( CONSTS(55)*STATES(13)+ CONSTS(56)*STATES(14))*1.00000)/(1.00000+EXP((STATES(8)+9.00000)/22.4000)))*(STATES(8) - ALGBRC(44))
      ALGBRC(47) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log((STATES(9)+ CONSTS(50)*CONSTS(51))/(STATES(10)+ CONSTS(50)*STATES(11)))
      ALGBRC(59) =  CONSTS(57)*STATES(15) ** 2.00000*(STATES(8) - ALGBRC(47))
      ALGBRC(69) =  ((( (1.00000 - CONSTS(70))*CONSTS(65)*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(10)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  STATES(9)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(76) =  ((( CONSTS(70)*CONSTS(65)*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(10)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  STATES(9)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(81) = ( (( CONSTS(77)*STATES(9))/(CONSTS(78)+STATES(9)))*STATES(11))/(CONSTS(79)+STATES(11))
      RATES(9) = ( 1.00000*(((ALGBRC(57)+ALGBRC(59)+ALGBRC(53)+ALGBRC(80)) -  (1.00000/(CONSTS(100) - 1.00000))*ALGBRC(81))+ALGBRC(69)+ALGBRC(76)))/( 1.00000*CONSTS(133)*CONSTS(48)) -  0.700000*(STATES(9) - CONSTS(101))
      RATES(10) =  (- 1.00000/( 1.00000*CONSTS(132)*CONSTS(48)))*((ALGBRC(53)+ALGBRC(57)+ALGBRC(59)+ALGBRC(69)+ALGBRC(76)+ALGBRC(80)) -  (1.00000/(CONSTS(100) - 1.00000))*ALGBRC(81))
      ALGBRC(51) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log((CONSTS(51)+ 0.120000*STATES(9))/(STATES(11)+ 0.120000*STATES(10)))
      ALGBRC(61) =  CONSTS(60)*STATES(16) ** 3.00000*STATES(17)*(STATES(8) - ALGBRC(51))
      ALGBRC(42) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log(CONSTS(51)/STATES(11))
      ALGBRC(65) =  CONSTS(63)*(STATES(8) - ALGBRC(42))
      ALGBRC(63) =  (( CONSTS(62)*1.00000)/(1.00000+EXP(- (STATES(8)+52.0000)/8.00000)))*(STATES(8) - ALGBRC(42))
      ALGBRC(71) =  ((( (1.00000 - CONSTS(70))*CONSTS(66)*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(11)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(51)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(77) =  ((( CONSTS(70)*CONSTS(66)*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(11)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(51)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(82) = ( (1.00000 - CONSTS(84))*CONSTS(80)*( EXP(( CONSTS(83)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*STATES(11) ** CONSTS(81)*CONSTS(52) -  EXP(( (CONSTS(83) - 1.00000)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*CONSTS(51) ** CONSTS(81)*STATES(12)))/( (1.00000+ CONSTS(82)*( STATES(12)*CONSTS(51) ** CONSTS(81)+ CONSTS(52)*STATES(11) ** CONSTS(81)))*(1.00000+STATES(12)/0.00690000))
      ALGBRC(83) = ( CONSTS(84)*CONSTS(80)*( EXP(( CONSTS(83)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*STATES(11) ** CONSTS(81)*CONSTS(52) -  EXP(( (CONSTS(83) - 1.00000)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*CONSTS(51) ** CONSTS(81)*STATES(18)))/( (1.00000+ CONSTS(82)*( STATES(18)*CONSTS(51) ** CONSTS(81)+ CONSTS(52)*STATES(11) ** CONSTS(81)))*(1.00000+STATES(18)/0.00690000))
      ALGBRC(67) =  ((( (1.00000 - CONSTS(70))*4.00000*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47)))))*( STATES(12)*EXP(( 100.000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(52)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(74) =  ((( CONSTS(70)*4.00000*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47)))))*( STATES(12)*EXP(( 100.000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(52)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(49) =  (( 0.500000*CONSTS(46)*CONSTS(47))/CONSTS(48))*log(CONSTS(52)/STATES(12))
      ALGBRC(79) =  CONSTS(74)*(STATES(8) - ALGBRC(49))
      ALGBRC(40) = TERNRY(VOI.GE.CONSTS(34).AND.VOI.LE.CONSTS(35).AND.(VOI - CONSTS(34)) -  INT((VOI - CONSTS(34))/CONSTS(20))*CONSTS(20).LE.CONSTS(36), CONSTS(37), 0.00000)
      RATES(8) =  (- 1.00000/CONSTS(49))*(ALGBRC(40)+ALGBRC(53)+ALGBRC(80)+ALGBRC(57)+ALGBRC(59)+ALGBRC(81)+ALGBRC(61)+ALGBRC(65)+ALGBRC(63)+ALGBRC(71)+ALGBRC(77)+ALGBRC(82)+ALGBRC(83)+ALGBRC(67)+ALGBRC(74)+ALGBRC(69)+ALGBRC(76)+ALGBRC(79))
      RATES(18) = ( - 1.00000*ALGBRC(74)+( 2.00000*ALGBRC(83))/(CONSTS(81) - 2.00000))/( 2.00000*1.00000*CONSTS(105)*CONSTS(132)*CONSTS(48)) -  STATES(18)*CONSTS(109)
      ALGBRC(85) = ALGBRC(82)+ALGBRC(83)
      RATES(11) =  (- 1.00000/( 1.00000*CONSTS(132)*CONSTS(48)))*(ALGBRC(61)+ALGBRC(63)+ALGBRC(65)+ALGBRC(71)+ALGBRC(77)+ (CONSTS(100)/(CONSTS(100) - 1.00000))*ALGBRC(81)+ (CONSTS(81)/(CONSTS(81) - 2.00000))*ALGBRC(85))
      ALGBRC(84) = STATES(12)+ STATES(25)*CONSTS(129)+ CONSTS(85)*CONSTS(86)+CONSTS(85)
      ALGBRC(86) = TERNRY(CONSTS(90).EQ.0.00000,  (STATES(12)/ALGBRC(84))*CONSTS(88) -  (( STATES(25)*CONSTS(129))/ALGBRC(84))*CONSTS(89), ( (STATES(12)/ALGBRC(84))*CONSTS(88))/(1.00000+STATES(25)/CONSTS(91)) -  (( STATES(25)*CONSTS(129))/ALGBRC(84))*CONSTS(89))
      ALGBRC(87) =  CONSTS(92)*(STATES(25) - STATES(26))
      RATES(25) =  (CONSTS(104)/CONSTS(108))*ALGBRC(86) - ALGBRC(87)
      ALGBRC(88) =  ( STATES(27)/(STATES(27)+0.250000) ** 2.00000*CONSTS(97)+CONSTS(96))*STATES(26)
      RATES(26) = ( (CONSTS(108)/CONSTS(106))*ALGBRC(87) - ALGBRC(88))/(1.00000+( CONSTS(130)*CONSTS(95))/STATES(26)+CONSTS(130) ** 2.00000)
      ALGBRC(89) = ( CONSTS(31)*STATES(2))/( ALGBRC(36)*STATES(7))
      ALGBRC(90) = TERNRY(ALGBRC(89).LE.0.00000, CONSTS(110), TERNRY(ALGBRC(89).LE.1.00000, CONSTS(110) ** ALGBRC(89), 1.00000)
      RATES(7) =  CONSTS(111)*(CONSTS(31) - STATES(7))*STATES(12) -  CONSTS(112)*EXP( - CONSTS(28)*STATES(7))*ALGBRC(90)*STATES(7)
      RATES(29) =  CONSTS(114)*(CONSTS(113) - STATES(29))*STATES(12) -  CONSTS(115)*STATES(29)
      RATES(30) =  CONSTS(117)*(CONSTS(116) - STATES(30))*STATES(12) -  CONSTS(118)*STATES(30)
      RATES(12) = (((( (- 1.00000/( 2.00000*1.00000*CONSTS(132)*CONSTS(48)))*((ALGBRC(67)+ALGBRC(79)) -  (2.00000/(CONSTS(81) - 2.00000))*ALGBRC(82))+ STATES(18)*CONSTS(105)*CONSTS(109)+( ALGBRC(88)*CONSTS(106))/CONSTS(104)) - RATES(7)) - RATES(29)) - RATES(30)) - ALGBRC(86)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(4) = STATES(1)
      ALGBRC(5) = 50.0000/(1.00000+EXP(- (STATES(8) - 5.00000)/9.00000))
      ALGBRC(17) =  0.0500000*EXP(- (STATES(8) - 20.0000)/15.0000)
      ALGBRC(6) = 50.0000/(1.00000+EXP(- (STATES(8) - 5.00000)/9.00000))
      ALGBRC(18) =  0.400000*EXP(- (STATES(8)+30.0000)/30.0000 ** 3.00000)
      ALGBRC(7) = 14.0000/(1.00000+EXP(- (STATES(8) - 40.0000)/9.00000))
      ALGBRC(19) =  1.00000*EXP(- STATES(8)/45.0000)
      ALGBRC(9) =  20.0000*EXP( - 0.125000*(STATES(8)+75.0000))
      ALGBRC(21) = 2000.00/(1.00000+ 320.000*EXP( - 0.100000*(STATES(8)+75.0000)))
      ALGBRC(12) =  0.0330000*EXP(- STATES(8)/17.0000)
      ALGBRC(24) = 33.0000/(1.00000+EXP( - 0.125000*(STATES(8)+10.0000)))
      ALGBRC(15) =  CONSTS(120)*(STATES(31)+20.0000)+ CONSTS(122)*(STATES(31) - CONSTS(123))
      ALGBRC(26) =  - CONSTS(119)*(STATES(8) - STATES(31))
      ALGBRC(8) = STATES(8)+41.0000
      ALGBRC(20) = TERNRY(ABS(ALGBRC(8)).LT.CONSTS(61), 2000.00, ( 200.000*ALGBRC(8))/(1.00000 - EXP( - 0.100000*ALGBRC(8))))
      ALGBRC(28) =  8000.00*EXP( - 0.0560000*(STATES(8)+66.0000))
      ALGBRC(10) = (STATES(8)+24.0000) - 5.00000
      ALGBRC(22) = TERNRY(ABS(ALGBRC(10)).LT.0.000100000, 120.000, ( 30.0000*ALGBRC(10))/(1.00000 - EXP(- ALGBRC(10)/4.00000)))
      ALGBRC(29) = TERNRY(ABS(ALGBRC(10)).LT.0.000100000, 120.000, ( 12.0000*ALGBRC(10))/(EXP(ALGBRC(10)/10.0000) - 1.00000))
      ALGBRC(11) = STATES(8)+34.0000
      ALGBRC(23) = TERNRY(ABS(ALGBRC(11)).LT.CONSTS(73), 25.0000, ( 6.25000*ALGBRC(11))/(EXP(ALGBRC(11)/4.00000) - 1.00000))
      ALGBRC(30) = 12.0000/(1.00000+EXP(( - 1.00000*(STATES(8)+34.0000))/4.00000))
      ALGBRC(14) = STATES(12)/(STATES(12)+CONSTS(98))
      ALGBRC(25) = STATES(18)/(STATES(18)+CONSTS(99))
      ALGBRC(31) = ALGBRC(14)+ (1.00000 - ALGBRC(14))*ALGBRC(25)
      ALGBRC(35) = 60.0000+ 500.000*ALGBRC(31) ** 2.00000
      ALGBRC(37) = TERNRY(STATES(8).LT.- 50.0000, 5.00000, 1.00000)
      ALGBRC(33) =  500.000*ALGBRC(31) ** 2.00000
      ALGBRC(39) = (1.00000 - STATES(27)) - STATES(28)
      ALGBRC(52) = TERNRY(STATES(3).LE.STATES(1),  CONSTS(43)*EXP( CONSTS(42)*(STATES(5) - STATES(4))),  CONSTS(45)*EXP( CONSTS(44)*(STATES(5) - STATES(4))))
      ALGBRC(38) = STATES(5)+STATES(6)
      ALGBRC(16) =  CONSTS(7)*(EXP( CONSTS(6)*STATES(6)) - 1.00000)
      ALGBRC(27) = ALGBRC(16)
      ALGBRC(41) = TERNRY(CONSTS(1).EQ.1.00000.AND.ALGBRC(27).GT.CONSTS(18).AND.ALGBRC(38).LE. CONSTS(19)*(1.00000+0.000100000), 1.00000, 0.00000)
      ALGBRC(48) = TERNRY(ALGBRC(41).EQ.1.00000, ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))), ( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))*STATES(1))/( CONSTS(2)*CONSTS(3)*EXP( CONSTS(2)*(STATES(5) - STATES(4)))+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6))))
      ALGBRC(55) = TERNRY(ALGBRC(52).EQ.0.00000, ALGBRC(48), STATES(3))
      ALGBRC(56) = TERNRY(ALGBRC(41).EQ.1.00000, 0.00000, TERNRY(ALGBRC(41).EQ.0.00000.AND.ALGBRC(52).EQ.0.00000, - ALGBRC(48), TERNRY(ALGBRC(41).EQ.0.00000.AND.ALGBRC(52) .NE. 0.00000, - STATES(3), 0.0/0.0)
      ALGBRC(32) = ( STATES(7)/CONSTS(31) ** CONSTS(9)*(1.00000+CONSTS(10) ** CONSTS(9)))/(STATES(7)/CONSTS(31) ** CONSTS(9)+CONSTS(10) ** CONSTS(9))
      ALGBRC(34) = TERNRY( CONSTS(32)*STATES(4)+CONSTS(33).LT.0.00000, 0.00000, TERNRY( CONSTS(32)*STATES(4)+CONSTS(33).LT.1.00000,  CONSTS(32)*STATES(4)+CONSTS(33), 1.00000)
      ALGBRC(36) = TERNRY(STATES(4).GT.0.550000, (STATES(4)+CONSTS(21))/(0.460000+CONSTS(21)),  (CONSTS(21)+0.550000)*1.00000)
      ALGBRC(3) = TERNRY(STATES(1).LE.0.00000, CONSTS(22) - ( CONSTS(23)*STATES(1))/CONSTS(14), TERNRY(STATES(1).LE.CONSTS(124).AND.0.00000.LT.STATES(1), ( (CONSTS(25) - CONSTS(24))*STATES(1))/CONSTS(124)+CONSTS(24), CONSTS(25)/1.00000+( CONSTS(30)*(STATES(1) - CONSTS(124)))/CONSTS(14) ** CONSTS(29))
      ALGBRC(54) = TERNRY(STATES(1).LE.0.00000, ( CONSTS(15)*(1.00000+STATES(1)/CONSTS(14)))/(CONSTS(15) - STATES(1)/CONSTS(14)), (1.00000+CONSTS(16)) - ( CONSTS(16) ** 2.00000*CONSTS(15))/( (( CONSTS(15)*CONSTS(16))/CONSTS(131))*STATES(1)/CONSTS(14) ** 2.00000+( (CONSTS(15)+1.00000)*STATES(1))/CONSTS(14)+ CONSTS(15)*CONSTS(16)))
      ALGBRC(58) = TERNRY(- CONSTS(14).LE.STATES(1).AND.STATES(1).LE.0.00000, 1.00000+( 0.600000*STATES(1))/CONSTS(14), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(125), ALGBRC(54)/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000), ( ALGBRC(54)*EXP( - CONSTS(27)*(STATES(1) - CONSTS(125))/CONSTS(14) ** CONSTS(17)))/(( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14))+1.00000))
      ALGBRC(60) =  CONSTS(11)*CONSTS(12)*ALGBRC(3)*CONSTS(13)*ALGBRC(58)
      ALGBRC(62) =  CONSTS(12)*ALGBRC(3)*(1.00000 -  CONSTS(11)*CONSTS(13)*ALGBRC(58))
      ALGBRC(64) =  ALGBRC(60)*ALGBRC(32)*ALGBRC(34)*ALGBRC(36)*(1.00000 - STATES(2)) -  ALGBRC(62)*STATES(2)
      ALGBRC(45) = TERNRY(STATES(1).LE.0.00000,  CONSTS(39)*EXP( CONSTS(38)*STATES(4)),  CONSTS(41)*EXP( CONSTS(40)*STATES(4)))
      ALGBRC(70) = ALGBRC(54)/ALGBRC(58)
      ALGBRC(66) = ( CONSTS(15)*1.00000*(1.00000+ 0.400000*CONSTS(15)+( 1.20000*STATES(1))/CONSTS(14)+ 0.600000*STATES(1)/CONSTS(14) ** 2.00000))/( CONSTS(14)* (CONSTS(15) - STATES(1)/CONSTS(14))*(1.00000+( 0.600000*STATES(1))/CONSTS(14)) ** 2.00000)
      ALGBRC(68) =  (1.00000/CONSTS(14))*EXP( - CONSTS(27)*STATES(1)/CONSTS(14) - CONSTS(125)/CONSTS(14) ** CONSTS(17))*(( 0.400000*CONSTS(15)+1.00000)/CONSTS(15)+ CONSTS(27)*CONSTS(17)*(1.00000+( ( 0.400000*CONSTS(15)+1.00000)*STATES(1))/( CONSTS(15)*CONSTS(14)))*STATES(1)/CONSTS(14) - CONSTS(125)/CONSTS(14) ** CONSTS(17) - 1.00000)
      ALGBRC(73) = TERNRY(STATES(1).LE.- CONSTS(14), CONSTS(126), TERNRY(- CONSTS(14).LT.STATES(1).AND.STATES(1).LE.0.00000, ALGBRC(66), TERNRY(0.00000.LT.STATES(1).AND.STATES(1).LE.CONSTS(125), CONSTS(127), ALGBRC(68))
      ALGBRC(43) = TERNRY(STATES(1).LE.0.00000, CONSTS(38), CONSTS(40))
      ALGBRC(75) = TERNRY(ALGBRC(41).EQ.1.00000, - ( CONSTS(8)*ALGBRC(64)*ALGBRC(70)+ ALGBRC(43)*ALGBRC(45)*STATES(1) ** 2.00000+ CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(73)+ALGBRC(45)), - ( CONSTS(8)*ALGBRC(64)*ALGBRC(70)+ ALGBRC(43)*ALGBRC(45)*STATES(1) ** 2.00000+ ( CONSTS(4)*CONSTS(5)*EXP( CONSTS(4)*STATES(5))+ CONSTS(6)*CONSTS(7)*EXP( CONSTS(6)*STATES(6)))*STATES(3))/( CONSTS(8)*STATES(2)*ALGBRC(73)+ALGBRC(45)))
      ALGBRC(50) = TERNRY(STATES(3).LE.STATES(1), CONSTS(42), CONSTS(44))
      ALGBRC(44) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log(STATES(9)/STATES(10))
      ALGBRC(53) = ( (( CONSTS(54)*STATES(9))/(STATES(9)+CONSTS(53)))*(STATES(8) - ALGBRC(44)))/(1.00000+EXP(( ((STATES(8) - ALGBRC(44)) - 10.0000)*CONSTS(48)*1.25000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(80) =  CONSTS(75)*(CONSTS(76)+ STATES(23)*(1.00000 - CONSTS(76)))*STATES(24)*(STATES(8) - ALGBRC(44))
      ALGBRC(57) =  (( ( CONSTS(55)*STATES(13)+ CONSTS(56)*STATES(14))*1.00000)/(1.00000+EXP((STATES(8)+9.00000)/22.4000)))*(STATES(8) - ALGBRC(44))
      ALGBRC(47) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log((STATES(9)+ CONSTS(50)*CONSTS(51))/(STATES(10)+ CONSTS(50)*STATES(11)))
      ALGBRC(59) =  CONSTS(57)*STATES(15) ** 2.00000*(STATES(8) - ALGBRC(47))
      ALGBRC(69) =  ((( (1.00000 - CONSTS(70))*CONSTS(65)*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(10)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  STATES(9)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(76) =  ((( CONSTS(70)*CONSTS(65)*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(10)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  STATES(9)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(81) = ( (( CONSTS(77)*STATES(9))/(CONSTS(78)+STATES(9)))*STATES(11))/(CONSTS(79)+STATES(11))
      ALGBRC(51) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log((CONSTS(51)+ 0.120000*STATES(9))/(STATES(11)+ 0.120000*STATES(10)))
      ALGBRC(61) =  CONSTS(60)*STATES(16) ** 3.00000*STATES(17)*(STATES(8) - ALGBRC(51))
      ALGBRC(42) =  (( CONSTS(46)*CONSTS(47))/CONSTS(48))*log(CONSTS(51)/STATES(11))
      ALGBRC(65) =  CONSTS(63)*(STATES(8) - ALGBRC(42))
      ALGBRC(63) =  (( CONSTS(62)*1.00000)/(1.00000+EXP(- (STATES(8)+52.0000)/8.00000)))*(STATES(8) - ALGBRC(42))
      ALGBRC(71) =  ((( (1.00000 - CONSTS(70))*CONSTS(66)*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(11)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(51)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(77) =  ((( CONSTS(70)*CONSTS(66)*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))))*( STATES(11)*EXP(( 50.0000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(51)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47))))
      ALGBRC(82) = ( (1.00000 - CONSTS(84))*CONSTS(80)*( EXP(( CONSTS(83)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*STATES(11) ** CONSTS(81)*CONSTS(52) -  EXP(( (CONSTS(83) - 1.00000)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*CONSTS(51) ** CONSTS(81)*STATES(12)))/( (1.00000+ CONSTS(82)*( STATES(12)*CONSTS(51) ** CONSTS(81)+ CONSTS(52)*STATES(11) ** CONSTS(81)))*(1.00000+STATES(12)/0.00690000))
      ALGBRC(83) = ( CONSTS(84)*CONSTS(80)*( EXP(( CONSTS(83)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*STATES(11) ** CONSTS(81)*CONSTS(52) -  EXP(( (CONSTS(83) - 1.00000)*(CONSTS(81) - 2.00000)*STATES(8)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))*CONSTS(51) ** CONSTS(81)*STATES(18)))/( (1.00000+ CONSTS(82)*( STATES(18)*CONSTS(51) ** CONSTS(81)+ CONSTS(52)*STATES(11) ** CONSTS(81)))*(1.00000+STATES(18)/0.00690000))
      ALGBRC(67) =  ((( (1.00000 - CONSTS(70))*4.00000*CONSTS(64)*STATES(19)*STATES(20)*STATES(21)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47)))))*( STATES(12)*EXP(( 100.000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(52)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(74) =  ((( CONSTS(70)*4.00000*CONSTS(64)*STATES(19)*STATES(20)*STATES(22)*(STATES(8) - 50.0000)*CONSTS(48))/( CONSTS(46)*CONSTS(47)))/(1.00000 - EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47)))))*( STATES(12)*EXP(( 100.000*CONSTS(48))/( CONSTS(46)*CONSTS(47))) -  CONSTS(52)*EXP(( - (STATES(8) - 50.0000)*CONSTS(48)*2.00000)/( CONSTS(46)*CONSTS(47))))
      ALGBRC(49) =  (( 0.500000*CONSTS(46)*CONSTS(47))/CONSTS(48))*log(CONSTS(52)/STATES(12))
      ALGBRC(79) =  CONSTS(74)*(STATES(8) - ALGBRC(49))
      ALGBRC(40) = TERNRY(VOI.GE.CONSTS(34).AND.VOI.LE.CONSTS(35).AND.(VOI - CONSTS(34)) -  INT((VOI - CONSTS(34))/CONSTS(20))*CONSTS(20).LE.CONSTS(36), CONSTS(37), 0.00000)
      ALGBRC(85) = ALGBRC(82)+ALGBRC(83)
      ALGBRC(84) = STATES(12)+ STATES(25)*CONSTS(129)+ CONSTS(85)*CONSTS(86)+CONSTS(85)
      ALGBRC(86) = TERNRY(CONSTS(90).EQ.0.00000,  (STATES(12)/ALGBRC(84))*CONSTS(88) -  (( STATES(25)*CONSTS(129))/ALGBRC(84))*CONSTS(89), ( (STATES(12)/ALGBRC(84))*CONSTS(88))/(1.00000+STATES(25)/CONSTS(91)) -  (( STATES(25)*CONSTS(129))/ALGBRC(84))*CONSTS(89))
      ALGBRC(87) =  CONSTS(92)*(STATES(25) - STATES(26))
      ALGBRC(88) =  ( STATES(27)/(STATES(27)+0.250000) ** 2.00000*CONSTS(97)+CONSTS(96))*STATES(26)
      ALGBRC(89) = ( CONSTS(31)*STATES(2))/( ALGBRC(36)*STATES(7))
      ALGBRC(90) = TERNRY(ALGBRC(89).LE.0.00000, CONSTS(110), TERNRY(ALGBRC(89).LE.1.00000, CONSTS(110) ** ALGBRC(89), 1.00000)
      ALGBRC(1) =  CONSTS(3)*(EXP( CONSTS(2)*(STATES(5) - STATES(4))) - 1.00000)
      ALGBRC(2) =  CONSTS(5)*(EXP( CONSTS(4)*STATES(5)) - 1.00000)
      ALGBRC(13) = EXP( 0.0800000*(STATES(8) - 40.0000))
      ALGBRC(46) =  (( CONSTS(58)*STATES(11))/(STATES(11)+CONSTS(59)))*(STATES(8) - ALGBRC(44))
      ALGBRC(72) =  CONSTS(8)*ALGBRC(70)*STATES(2)
      ALGBRC(78) = ALGBRC(67)+ALGBRC(69)+ALGBRC(71)+ALGBRC(74)+ALGBRC(76)+ALGBRC(77)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END