C
C There are a total of 0 entries in the algebraic variable array.
C There are a total of 0 entries in each of the rate and state variable arrays.
C There are a total of 156 entries in the constant variable array.
C
C
C CONSTS(1) is elementary_charge in component codata_2006_electromagnetic (coulomb).
C CONSTS(2) is magnetic_flux_quantum in component codata_2006_electromagnetic (weber).
C CONSTS(3) is conductance_quantum in component codata_2006_electromagnetic (siemens).
C CONSTS(4) is inverse_of_conductance_quantum in component codata_2006_electromagnetic (ohm).
C CONSTS(5) is Josephson_constant in component codata_2006_electromagnetic (hertz_per_volt).
C CONSTS(6) is von_Klitzing_constant in component codata_2006_electromagnetic (ohm).
C CONSTS(7) is Bohr_magneton in component codata_2006_electromagnetic (joule_per_tesla).
C CONSTS(8) is nuclear_magneton in component codata_2006_electromagnetic (joule_per_tesla).
C CONSTS(9) is speed_of_light_in_vacuum in component codata_2006_universal (metre_per_second).
C CONSTS(10) is magnetic_constant in component codata_2006_universal (newton_per_square_ampere).
C CONSTS(11) is electric_constant in component codata_2006_universal (farad_per_metre).
C CONSTS(12) is characteristic_impedance_of_vacuum in component codata_2006_universal (ohm).
C CONSTS(13) is Newtonian_constant_of_gravitation in component codata_2006_universal (cubic_metre_per_kilogram_per_square_second).
C CONSTS(14) is Planck_constant in component codata_2006_universal (joule_second).
C CONSTS(15) is Planck_constant_over_2_pi in component codata_2006_universal (joule_second).
C CONSTS(16) is Planck_mass in component codata_2006_universal (kilogram).
C CONSTS(17) is Planck_temperature in component codata_2006_universal (kelvin).
C CONSTS(18) is Planck_length in component codata_2006_universal (metre).
C CONSTS(19) is Planck_time in component codata_2006_universal (second).
C CONSTS(20) is fine_structure_constant in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(21) is inverse_fine_structure_constant in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(22) is Rydberg_constant in component codata_2006_atomic_and_nuclear (per_metre).
C CONSTS(23) is Bohr_radius in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(24) is Hartree_energy in component codata_2006_atomic_and_nuclear (joule).
C CONSTS(25) is quantum_of_circulation in component codata_2006_atomic_and_nuclear (square_metre_per_second).
C CONSTS(26) is Fermi_coupling_constant in component codata_2006_atomic_and_nuclear (per_square_gigaelectronvolt).
C CONSTS(27) is weak_mixing_angle in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(28) is electron_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(29) is electron_muon_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(30) is electron_tau_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(31) is electron_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(32) is electron_neutron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(33) is electron_deuteron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(34) is electron_alpha_particle_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(35) is electron_charge_to_mass_quotient in component codata_2006_atomic_and_nuclear (coulomb_per_kilogram).
C CONSTS(36) is electron_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(37) is Compton_wavelength in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(38) is classic_electron_radius in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(39) is Thomson_cross_section in component codata_2006_atomic_and_nuclear (square_metre).
C CONSTS(40) is electron_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(41) is electron_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(42) is electron_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(43) is electron_magnetic_moment_anomaly in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(44) is electron_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(45) is electron_muon_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(46) is electron_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(47) is electron_to_shielded_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(48) is electron_neutron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(49) is electron_deuteron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(50) is electron_to_shielded_helion_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(51) is electron_gyromagnetic_ratio in component codata_2006_atomic_and_nuclear (per_second_per_tesla).
C CONSTS(52) is muon_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(53) is muon_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(54) is muon_tau_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(55) is muon_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(56) is muon_neutron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(57) is muon_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(58) is muon_Compton_wavelength in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(59) is muon_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(60) is muon_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(61) is muon_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(62) is muon_magnetic_moment_anomaly in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(63) is muon_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(64) is muon_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(65) is tau_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(66) is tau_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(67) is tau_muon_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(68) is tau_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(69) is tau_neutron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(70) is tau_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(71) is tau_Compton_wavelength in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(72) is proton_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(73) is proton_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(74) is proton_muon_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(75) is proton_tau_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(76) is proton_neutron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(77) is proton_charge_to_mass_quotient in component codata_2006_atomic_and_nuclear (coulomb_per_kilogram).
C CONSTS(78) is proton_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(79) is proton_Compton_wavelength in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(80) is proton_rms_electron_radius in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(81) is proton_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(82) is proton_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(83) is proton_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(84) is proton_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(85) is proton_neutron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(86) is shielded_proton_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(87) is shielded_proton_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(88) is shielded_proton_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(89) is proton_magnetic_shielding_correction in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(90) is proton_gyromagnetic_ratio in component codata_2006_atomic_and_nuclear (per_second_per_tesla).
C CONSTS(91) is shielded_proton_gyromagnetic_ratio in component codata_2006_atomic_and_nuclear (per_second_per_tesla).
C CONSTS(92) is neutron_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(93) is neutron_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(94) is neutron_muon_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(95) is neutron_tau_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(96) is neutron_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(97) is neutron_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(98) is neutron_Compton_wavelength in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(99) is neutron_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(100) is neutron_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(101) is neutron_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(102) is neutron_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(103) is neutron_electron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(104) is neutron_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(105) is neutron_to_shielded_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(106) is neutron_gyromagnetic_ratio in component codata_2006_atomic_and_nuclear (per_second_per_tesla).
C CONSTS(107) is deuteron_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(108) is deuteron_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(109) is deuteron_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(110) is deuteron_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(111) is deuteron_rms_charge_radius in component codata_2006_atomic_and_nuclear (metre).
C CONSTS(112) is deuteron_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(113) is deuteron_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(114) is deuteron_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(115) is deuteron_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(116) is deuteron_electron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(117) is deuteron_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(118) is deuteron_neutron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(119) is triton_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(120) is triton_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(121) is triton_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(122) is triton_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(123) is triton_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(124) is triton_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(125) is triton_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(126) is triton_g_factor in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(127) is triton_electron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(128) is triton_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(129) is triton_neutron_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(130) is helion_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(131) is helion_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(132) is helion_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(133) is helion_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(134) is shielded_helion_magnetic_moment in component codata_2006_atomic_and_nuclear (joule_per_tesla).
C CONSTS(135) is shielded_helion_magnetic_moment_to_Bohr_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(136) is shielded_helion_magnetic_moment_to_nuclear_magneton_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(137) is shielded_helion_to_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(138) is shielded_helion_to_shielded_proton_magnetic_moment_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(139) is shielded_helion_gyromagnetic_ratio in component codata_2006_atomic_and_nuclear (per_second_per_tesla).
C CONSTS(140) is alpha_particle_mass in component codata_2006_atomic_and_nuclear (kilogram).
C CONSTS(141) is alpha_particle_to_electron_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(142) is alpha_particle_to_proton_mass_ratio in component codata_2006_atomic_and_nuclear (dimensionless).
C CONSTS(143) is alpha_particle_molar_mass in component codata_2006_atomic_and_nuclear (kilogram_per_mole).
C CONSTS(144) is Avogadro_constant in component codata_2006_physicochemical (per_mole).
C CONSTS(145) is atomic_mass_constant in component codata_2006_physicochemical (kilogram).
C CONSTS(146) is Faraday_constant in component codata_2006_physicochemical (coulomb_per_mole).
C CONSTS(147) is molar_Planck_constant in component codata_2006_physicochemical (joule_second_per_mole).
C CONSTS(148) is molar_gas_constant in component codata_2006_physicochemical (joule_per_mole_per_kelvin).
C CONSTS(149) is Boltzmann_constant in component codata_2006_physicochemical (joule_per_kelvin).
C CONSTS(150) is molar_volume_of_ideal_gas in component codata_2006_physicochemical (cubic_metre_per_mole).
C CONSTS(151) is Sackur_Tetrode_constant in component codata_2006_physicochemical (dimensionless).
C CONSTS(152) is Stefan_Boltzmann_constant in component codata_2006_physicochemical (watt_per_square_meter_per_quartic_kelvin).
C CONSTS(153) is first_radiation_constant in component codata_2006_physicochemical (watt_per_square_metre).
C CONSTS(154) is first_radiation_constant_for_spectral_radiance in component codata_2006_physicochemical (watt_per_square_metre_per_steradian).
C CONSTS(155) is second_radiation_constant in component codata_2006_physicochemical (metre_kelvin).
C CONSTS(156) is Wien_displacement_law_constant in component codata_2006_physicochemical (metre_kelvin).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 1.602176487e-19
      CONSTS(2) = 2.067833667e-15
      CONSTS(3) = 7.7480917004e-5
      CONSTS(4) = 12906.4037787
      CONSTS(5) = 483597.891e9
      CONSTS(6) = 25812.807557
      CONSTS(7) = 927.400915e-26
      CONSTS(8) = 5.05078324e-27
      CONSTS(9) = 299792458
      CONSTS(10) = 1.2566370614359172953850573533118011536788677597500423283899778369231265625145e-6
      CONSTS(11) = 8.8541878176203898505365630317107502606083701665994498081024171524053950954600e-12
      CONSTS(12) = 376.73031346177065546819840042031930826862350835241865523207463829670726922131
      CONSTS(13) = 6.67428e-11
      CONSTS(14) = 6.62606896e-34
      CONSTS(15) = 1.054571628e-34
      CONSTS(16) = 6.62606896e-34
      CONSTS(17) = 1.416785e32
      CONSTS(18) = 1.616252e-35
      CONSTS(19) = 5.39124e-44
      CONSTS(20) = 7.2973525376e-3
      CONSTS(21) = 137.035999679
      CONSTS(22) = 10973731.568527
      CONSTS(23) = 0.52917720859e-10
      CONSTS(24) = 4.35974394e-18
      CONSTS(25) = 3.6369475199e-4
      CONSTS(26) = 1.6637e-5
      CONSTS(27) = 0.22255
      CONSTS(28) = 9.10938215e-31
      CONSTS(29) = 4.83633171e-3
      CONSTS(30) = 2.87564e-4
      CONSTS(31) = 5.4461702177e-4
      CONSTS(32) = 5.4386734459e-4
      CONSTS(33) = 2.7244371093e-4
      CONSTS(34) = 1.37093355570e-4
      CONSTS(35) = -1.758820150e11
      CONSTS(36) = 5.4857990943e-7
      CONSTS(37) = 2.4263102175e-12
      CONSTS(38) = 2.4263102175e-12
      CONSTS(39) = 0.6652458558e-28
      CONSTS(40) = -928.476377e-26
      CONSTS(41) = -1.00115965218111
      CONSTS(42) = -1838.28197092
      CONSTS(43) = 1.15965218111e-3
      CONSTS(44) = -2.0023193043622
      CONSTS(45) = 206.7669877
      CONSTS(46) = -658.2106848
      CONSTS(47) = -658.2275971
      CONSTS(48) = 960.92050
      CONSTS(49) = -2143.923498
      CONSTS(50) = 864.058257
      CONSTS(51) = 1.760859770e11
      CONSTS(52) = 1.88353130e-28
      CONSTS(53) = 206.7682823
      CONSTS(54) = 5.94592e-2
      CONSTS(55) = 0.1126095261
      CONSTS(56) = 0.1124545167
      CONSTS(57) = 0.1134289256e-3
      CONSTS(58) = 11.73444104e-15
      CONSTS(59) = -4.49044786e-26
      CONSTS(60) = -4.84197049e-3
      CONSTS(61) = -8.89059705
      CONSTS(62) = 1.16592069e-3
      CONSTS(63) = -2.0023318414
      CONSTS(64) = -3.183345137
      CONSTS(65) = 3.16777e-27
      CONSTS(66) = 3477.48
      CONSTS(67) = 16.8183
      CONSTS(68) = 1.89390
      CONSTS(69) = 1.89129
      CONSTS(70) = 1.90768e-3
      CONSTS(71) = 0.69772e-15
      CONSTS(72) = 1.672621637e-27
      CONSTS(73) = 1836.15267247
      CONSTS(74) = 8.88024339
      CONSTS(75) = 0.528012
      CONSTS(76) = 0.99862347824
      CONSTS(77) = 9.57883392e7
      CONSTS(78) = 1.00727646677e-3
      CONSTS(79) = 1.3214098446e-15
      CONSTS(80) = 0.8768e-15
      CONSTS(81) = 1.410606662e-26
      CONSTS(82) = 1.521032209e-3
      CONSTS(83) = 2.792847356
      CONSTS(84) = 5.585694713
      CONSTS(85) = -1.45989806
      CONSTS(86) = 1.410570419e-26
      CONSTS(87) = 1.520993128e-3
      CONSTS(88) = 2.792775598
      CONSTS(89) = 25.694e-6
      CONSTS(90) = 2.675222099e8
      CONSTS(91) = 2.675153362e8
      CONSTS(92) = 1.674927211e-27
      CONSTS(93) = 1838.6836605
      CONSTS(94) = 8.89248409
      CONSTS(95) = 0.528740
      CONSTS(96) = 1.00137841918
      CONSTS(97) = 1.00866491597e-3
      CONSTS(98) = 1.3195908951e-15
      CONSTS(99) = -0.96623641e-26
      CONSTS(100) = -1.04187563e-3
      CONSTS(101) = -1.91304273
      CONSTS(102) = -3.82608545
      CONSTS(103) = 1.04066882e-3
      CONSTS(104) = -0.68497934
      CONSTS(105) = -0.68499694
      CONSTS(106) = 1.83247185e8
      CONSTS(107) = 3.34358320e-27
      CONSTS(108) = 3670.4829654
      CONSTS(109) = 1.99900750108
      CONSTS(110) = 2.013553212724e-3
      CONSTS(111) = 2.1402e-15
      CONSTS(112) = 0.433073465e-26
      CONSTS(113) = -0.4669754556e-3
      CONSTS(114) = 0.8574382308
      CONSTS(115) = 0.8574382308
      CONSTS(116) = -4.664345537e-4
      CONSTS(117) = 0.3070122070
      CONSTS(118) = -0.44820652
      CONSTS(119) = 5.00735588e-27
      CONSTS(120) = 5496.9215269
      CONSTS(121) = 2.9937170309
      CONSTS(122) = 3.0155007134e-3
      CONSTS(123) = 1.504609361e-26
      CONSTS(124) = 1.622393657e-3
      CONSTS(125) = 2.978962448
      CONSTS(126) = 5.957924896
      CONSTS(127) = -1.620514423e-3
      CONSTS(128) = 1.066639908
      CONSTS(129) = -1.55718553
      CONSTS(130) = 5.00641192e-27
      CONSTS(131) = 5495.8852765
      CONSTS(132) = 2.9931526713
      CONSTS(133) = 3.0149322473e-3
      CONSTS(134) = -1.074552982e-26
      CONSTS(135) = -1.158671471e-3
      CONSTS(136) = -2.127497718
      CONSTS(137) = -0.761766558
      CONSTS(138) = -0.7617861313
      CONSTS(139) = 2.037894730e8
      CONSTS(140) = 6.64465620e-27
      CONSTS(141) = 7294.2995365
      CONSTS(142) = 3.97259968951
      CONSTS(143) = 4.001506179127e-3
      CONSTS(144) = 6.02214179e23
      CONSTS(145) = 1.660538782e-27
      CONSTS(146) = 96485.3399
      CONSTS(147) = 3.9903126821e-10
      CONSTS(148) = 8.314472
      CONSTS(149) = 1.3806504e-23
      CONSTS(150) = 22.413996e-3
      CONSTS(151) = -1.1648677
      CONSTS(152) = 5.670400e-8
      CONSTS(153) = 3.74177118e-16
      CONSTS(154) = 1.191042759e-16
      CONSTS(155) = 1.4387752e-2
      CONSTS(156) = 2.8977685e-3
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      RETURN
      END