Generated Code

The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)

The raw code is available.

C
C There are a total of 312 entries in the algebraic variable array.
C There are a total of 25 entries in each of the rate and state variable arrays.
C There are a total of 342 entries in the constant variable array.
C
C
C VOI is time in component environment (minute).
C ALGBRC(1) is pH_cy in component environment (dimensionless).
C CONSTS(35) is addbuffer in component environment (dimensionless).
C CONSTS(36) is pHstat in component environment (dimensionless).
C CONSTS(3) is Par_97 in component parameters (dimensionless).
C CONSTS(4) is Par_98 in component parameters (dimensionless).
C STATES(1) is pH_calc in component differential_equations (dimensionless).
C CONSTS(1) is R in component global_parameters (kilojoule_per_kelvin_per_mole).
C CONSTS(2) is T1 in component global_parameters (kelvin).
C CONSTS(131) is T in component global_parameters (kelvin).
C CONSTS(126) is I in component global_parameters (molar).
C CONSTS(125) is Par_90 in component parameters (molar).
C CONSTS(130) is Par_94 in component parameters (kelvin).
C CONSTS(5) is Par_1 in component parameters (molar_per_minute).
C CONSTS(6) is Par_2 in component parameters (dimensionless).
C CONSTS(7) is Par_3 in component parameters (molar).
C CONSTS(37) is Par_4 in component parameters (molar).
C CONSTS(57) is Par_5 in component parameters (molar).
C CONSTS(8) is Par_6 in component parameters (molar).
C CONSTS(58) is Par_7 in component parameters (molar).
C CONSTS(9) is Par_8 in component parameters (molar).
C CONSTS(59) is Par_9 in component parameters (molar).
C CONSTS(60) is Par_10 in component parameters (molar).
C CONSTS(10) is Par_11 in component parameters (molar).
C CONSTS(61) is Par_12 in component parameters (molar).
C CONSTS(62) is Par_13 in component parameters (molar).
C CONSTS(11) is Par_14 in component parameters (molar).
C CONSTS(63) is Par_15 in component parameters (molar).
C CONSTS(12) is Par_16 in component parameters (molar).
C CONSTS(13) is Par_17 in component parameters (dimensionless).
C CONSTS(14) is Par_18 in component parameters (dimensionless).
C CONSTS(64) is Par_19 in component parameters (molar_per_minute).
C CONSTS(65) is Par_20 in component parameters (molar).
C CONSTS(66) is Par_21 in component parameters (molar).
C CONSTS(67) is Par_22 in component parameters (molar_per_minute).
C CONSTS(68) is Par_23 in component parameters (molar).
C CONSTS(69) is Par_24 in component parameters (molar).
C CONSTS(70) is Par_25 in component parameters (molar_per_minute).
C CONSTS(71) is Par_26 in component parameters (molar).
C CONSTS(72) is Par_27 in component parameters (molar).
C CONSTS(73) is Par_28 in component parameters (molar).
C CONSTS(74) is Par_29 in component parameters (molar).
C CONSTS(75) is Par_30 in component parameters (molar).
C CONSTS(15) is Par_31 in component parameters (molar).
C CONSTS(76) is Par_32 in component parameters (molar).
C CONSTS(16) is Par_33 in component parameters (molar).
C CONSTS(17) is Par_34 in component parameters (molar).
C CONSTS(18) is Par_35 in component parameters (molar).
C CONSTS(19) is Par_36 in component parameters (dimensionless).
C CONSTS(20) is Par_37 in component parameters (dimensionless).
C CONSTS(21) is Par_38 in component parameters (dimensionless).
C CONSTS(77) is Par_39 in component parameters (molar_per_minute).
C CONSTS(78) is Par_40 in component parameters (molar).
C CONSTS(79) is Par_41 in component parameters (molar).
C CONSTS(80) is Par_42 in component parameters (molar).
C CONSTS(81) is Par_43 in component parameters (molar_per_minute).
C CONSTS(82) is Par_44 in component parameters (molar).
C CONSTS(83) is Par_45 in component parameters (molar).
C CONSTS(84) is Par_46 in component parameters (molar_per_minute).
C CONSTS(85) is Par_47 in component parameters (molar).
C CONSTS(86) is Par_48 in component parameters (molar).
C CONSTS(87) is Par_49 in component parameters (molar).
C CONSTS(88) is Par_50 in component parameters (molar).
C CONSTS(89) is Par_51 in component parameters (molar_per_minute).
C CONSTS(90) is Par_52 in component parameters (molar).
C CONSTS(91) is Par_53 in component parameters (molar).
C CONSTS(92) is Par_54 in component parameters (molar).
C CONSTS(93) is Par_55 in component parameters (molar).
C CONSTS(94) is Par_56 in component parameters (molar).
C CONSTS(95) is Par_57 in component parameters (molar_per_minute).
C CONSTS(96) is Par_58 in component parameters (molar).
C CONSTS(97) is Par_59 in component parameters (molar).
C CONSTS(98) is Par_60 in component parameters (molar).
C CONSTS(99) is Par_61 in component parameters (molar).
C CONSTS(100) is Par_62 in component parameters (molar_per_minute).
C CONSTS(101) is Par_63 in component parameters (molar).
C CONSTS(102) is Par_64 in component parameters (molar).
C CONSTS(103) is Par_65 in component parameters (molar_per_minute).
C CONSTS(104) is Par_66 in component parameters (molar).
C CONSTS(105) is Par_67 in component parameters (molar).
C CONSTS(106) is Par_68 in component parameters (molar_per_minute).
C CONSTS(107) is Par_69 in component parameters (molar).
C CONSTS(108) is Par_70 in component parameters (molar).
C CONSTS(109) is Par_71 in component parameters (molar).
C CONSTS(110) is Par_72 in component parameters (molar).
C CONSTS(111) is Par_73 in component parameters (molar_per_minute).
C CONSTS(112) is Par_74 in component parameters (molar).
C CONSTS(113) is Par_75 in component parameters (molar).
C CONSTS(114) is Par_76 in component parameters (molar).
C CONSTS(115) is Par_77 in component parameters (molar).
C CONSTS(116) is Par_78 in component parameters (molar_per_minute).
C CONSTS(117) is Par_79 in component parameters (molar).
C CONSTS(118) is Par_80 in component parameters (molar).
C CONSTS(119) is Par_81 in component parameters (molar).
C CONSTS(22) is Par_82 in component parameters (molar).
C CONSTS(120) is Par_83 in component parameters (molar).
C CONSTS(121) is Par_84 in component parameters (molar_per_minute).
C CONSTS(122) is Par_85 in component parameters (molar).
C CONSTS(123) is Par_86 in component parameters (molar).
C CONSTS(124) is Par_87 in component parameters (molar).
C CONSTS(23) is Par_88 in component parameters (molar_per_minute).
C CONSTS(24) is Par_89 in component parameters (molar).
C CONSTS(127) is Par_91 in component parameters (molar).
C CONSTS(128) is Par_92 in component parameters (molar).
C CONSTS(129) is Par_93 in component parameters (molar).
C CONSTS(25) is Par_95 in component parameters (molar).
C CONSTS(132) is Par_96 in component parameters (molar).
C CONSTS(133) is Par_99 in component parameters (dimensionless).
C CONSTS(26) is Par_100 in component parameters (dimensionless).
C CONSTS(38) is mgT in component equilibrium_constants (molar).
C CONSTS(134) is k in component equilibrium_constants (molar).
C CONSTS(135) is c0 in component equilibrium_constants (molar).
C CONSTS(136) is RT2dadT in component correction_factors (kilojoule_half_liter_per_3_half_mole).
C CONSTS(137) is B in component correction_factors (per_half_molar).
C CONSTS(138) is Icorr in component correction_factors (kilojoule_per_mole).
C CONSTS(139) is I1 in component correction_factors (molar).
C CONSTS(140) is alphadebye in component correction_factors (per_half_molar).
C CONSTS(141) is IcorrpKa in component correction_factors (kilojoule_per_mole).
C CONSTS(142) is TcorrpKa in component correction_factors (mole_per_kilojoule).
C CONSTS(143) is RTalpha in component correction_factors (kilojoule_half_liter_per_3_half_mole).
C CONSTS(144) is IcorrdeltaGpof in component correction_factors (kilojoule_per_mole).
C CONSTS(145) is pKak_Pi in component correction_factors (dimensionless).
C CONSTS(146) is deltaH1o_Pi in component correction_factors (kilojoule_per_mole).
C CONSTS(147) is deltaHmgo_Pi in component correction_factors (kilojoule_per_mole).
C CONSTS(148) is deltaH1_Pi in component correction_factors (kilojoule_per_mole).
C CONSTS(149) is deltaHmg_Pi in component correction_factors (kilojoule_per_mole).
C CONSTS(150) is pKa1_Pi in component correction_factors (dimensionless).
C CONSTS(151) is pKamg_Pi in component correction_factors (dimensionless).
C ALGBRC(2) is P_Pi in component correction_factors (dimensionless).
C ALGBRC(3) is HPi2 in component correction_factors (dimensionless).
C ALGBRC(4) is H2Pi1 in component correction_factors (dimensionless).
C ALGBRC(5) is kPi in component correction_factors (dimensionless).
C ALGBRC(6) is mgPi in component correction_factors (dimensionless).
C ALGBRC(7) is Navg_Pi in component correction_factors (dimensionless).
C ALGBRC(8) is dNavgPidH in component correction_factors (per_molar).
C ALGBRC(9) is dNavgPidmg in component correction_factors (per_molar).
C ALGBRC(10) is dmgPidmg in component correction_factors (per_molar).
C ALGBRC(11) is dmgPidpH in component correction_factors (dimensionless).
C CONSTS(152) is NH_HPi2 in component correction_factors (dimensionless).
C CONSTS(153) is deltaGof_HPi2 in component correction_factors (kilojoule_per_mole).
C ALGBRC(12) is deltaGpof_HPi2 in component correction_factors (kilojoule_per_mole).
C CONSTS(154) is deltaH1o_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(155) is deltaHmgo_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(156) is deltaHko_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(157) is deltaH1_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(158) is deltaHmg_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(159) is deltaHk_ATP in component correction_factors (kilojoule_per_mole).
C CONSTS(160) is pKa1_ATP in component correction_factors (dimensionless).
C CONSTS(161) is pKamg_ATP in component correction_factors (dimensionless).
C CONSTS(162) is pKak_ATP in component correction_factors (dimensionless).
C ALGBRC(13) is P_ATP in component correction_factors (dimensionless).
C ALGBRC(14) is ATP4 in component correction_factors (dimensionless).
C ALGBRC(15) is HATP3 in component correction_factors (dimensionless).
C ALGBRC(16) is mgATP2 in component correction_factors (dimensionless).
C ALGBRC(17) is kATP in component correction_factors (dimensionless).
C ALGBRC(18) is Navg_ATP in component correction_factors (dimensionless).
C ALGBRC(19) is dNavgATPdH in component correction_factors (per_molar).
C ALGBRC(20) is dNavgATPdmg in component correction_factors (per_molar).
C ALGBRC(21) is dmgATP2dmg in component correction_factors (per_molar).
C ALGBRC(22) is dmgATP2dpH in component correction_factors (dimensionless).
C CONSTS(163) is NH_ATP4 in component correction_factors (dimensionless).
C CONSTS(164) is deltaGof_ATP4 in component correction_factors (kilojoule_per_mole).
C ALGBRC(23) is deltaGpof_ATP4 in component correction_factors (kilojoule_per_mole).
C CONSTS(165) is pKak_ADP in component correction_factors (dimensionless).
C CONSTS(166) is deltaH1o_ADP in component correction_factors (kilojoule_per_mole).
C CONSTS(167) is deltaHmgo_ADP in component correction_factors (kilojoule_per_mole).
C CONSTS(168) is deltaH1_ADP in component correction_factors (kilojoule_per_mole).
C CONSTS(169) is deltaHmg_ADP in component correction_factors (kilojoule_per_mole).
C CONSTS(170) is pKa1_ADP in component correction_factors (dimensionless).
C CONSTS(171) is pKamg_ADP in component correction_factors (dimensionless).
C ALGBRC(24) is P_ADP in component correction_factors (dimensionless).
C ALGBRC(25) is ADP3 in component correction_factors (dimensionless).
C ALGBRC(26) is HADP2 in component correction_factors (dimensionless).
C ALGBRC(27) is mgADP in component correction_factors (dimensionless).
C ALGBRC(28) is kADP in component correction_factors (dimensionless).
C ALGBRC(29) is Navg_ADP in component correction_factors (dimensionless).
C ALGBRC(30) is dNavgADPdH in component correction_factors (per_molar).
C ALGBRC(32) is dNavgADPdmg in component correction_factors (per_molar).
C ALGBRC(33) is dmgADPdmg in component correction_factors (per_molar).
C ALGBRC(34) is dmgADPdpH in component correction_factors (dimensionless).
C CONSTS(172) is NH_ADP3 in component correction_factors (dimensionless).
C CONSTS(173) is deltaGof_ADP3 in component correction_factors (kilojoule_per_mole).
C ALGBRC(35) is deltaGpof_ADP3 in component correction_factors (kilojoule_per_mole).
C CONSTS(174) is deltaH1o_AMP in component correction_factors (kilojoule_per_mole).
C CONSTS(175) is deltaHmgo_AMP in component correction_factors (kilojoule_per_mole).
C CONSTS(176) is deltaH1_AMP in component correction_factors (kilojoule_per_mole).
C CONSTS(177) is deltaHmg_AMP in component correction_factors (kilojoule_per_mole).
C CONSTS(178) is pKa1_AMP in component correction_factors (dimensionless).
C CONSTS(179) is pKamg_AMP in component correction_factors (dimensionless).
C ALGBRC(36) is P_AMP in component correction_factors (dimensionless).
C ALGBRC(37) is AMP2 in component correction_factors (dimensionless).
C ALGBRC(38) is HAMP1 in component correction_factors (dimensionless).
C ALGBRC(39) is mgAMP in component correction_factors (dimensionless).
C ALGBRC(40) is Navg_AMP in component correction_factors (dimensionless).
C ALGBRC(41) is dNavgAMPdH in component correction_factors (per_molar).
C ALGBRC(42) is dNavgAMPdmg in component correction_factors (per_molar).
C ALGBRC(43) is dmgAMPdmg in component correction_factors (per_molar).
C ALGBRC(44) is dmgAMPdpH in component correction_factors (dimensionless).
C CONSTS(180) is NH_AMP2 in component correction_factors (dimensionless).
C CONSTS(181) is deltaGof_AMP2 in component correction_factors (kilojoule_per_mole).
C ALGBRC(45) is deltaGpof_AMP2 in component correction_factors (kilojoule_per_mole).
C CONSTS(182) is pKak_PCR in component correction_factors (dimensionless).
C CONSTS(183) is deltaH1o_PCR in component correction_factors (kilojoule_per_mole).
C CONSTS(184) is deltaHmgo_PCR in component correction_factors (kilojoule_per_mole).
C CONSTS(185) is deltaH1_PCR in component correction_factors (kilojoule_per_mole).
C CONSTS(186) is deltaHmg_PCR in component correction_factors (kilojoule_per_mole).
C CONSTS(187) is pKa1_PCR in component correction_factors (dimensionless).
C CONSTS(188) is pKamg_PCR in component correction_factors (dimensionless).
C ALGBRC(46) is P_PCR in component correction_factors (dimensionless).
C ALGBRC(47) is HPCR in component correction_factors (dimensionless).
C ALGBRC(49) is H2PCR in component correction_factors (dimensionless).
C ALGBRC(48) is kPCR in component correction_factors (dimensionless).
C ALGBRC(50) is mgPCR in component correction_factors (dimensionless).
C ALGBRC(51) is Navg_PCR in component correction_factors (dimensionless).
C ALGBRC(52) is dNavgPCRdH in component correction_factors (per_molar).
C ALGBRC(53) is dNavgPCRdmg in component correction_factors (per_molar).
C ALGBRC(54) is dmgPCRdmg in component correction_factors (per_molar).
C ALGBRC(55) is dmgPCRdpH in component correction_factors (dimensionless).
C CONSTS(189) is NH_HPCR in component correction_factors (dimensionless).
C CONSTS(190) is pKa1_CR in component correction_factors (dimensionless).
C ALGBRC(56) is P_CR in component correction_factors (dimensionless).
C ALGBRC(57) is HCR in component correction_factors (dimensionless).
C ALGBRC(58) is H2CR in component correction_factors (dimensionless).
C ALGBRC(59) is Navg_CR in component correction_factors (dimensionless).
C ALGBRC(60) is dNavgCRdH in component correction_factors (per_molar).
C CONSTS(27) is dNavgCRdmg in component correction_factors (per_molar).
C CONSTS(191) is NH_HCR in component correction_factors (dimensionless).
C CONSTS(192) is deltaH1o_G1P in component correction_factors (kilojoule_per_mole).
C CONSTS(193) is deltaHmgo_G1P in component correction_factors (kilojoule_per_mole).
C CONSTS(194) is deltaH1_G1P in component correction_factors (kilojoule_per_mole).
C CONSTS(195) is deltaHmg_G1P in component correction_factors (kilojoule_per_mole).
C CONSTS(196) is pKa1_G1P in component correction_factors (dimensionless).
C CONSTS(197) is pKamg_G1P in component correction_factors (dimensionless).
C ALGBRC(61) is P_G1P in component correction_factors (dimensionless).
C ALGBRC(62) is UG1P in component correction_factors (dimensionless).
C ALGBRC(63) is HG1P in component correction_factors (dimensionless).
C ALGBRC(64) is mgG1P in component correction_factors (dimensionless).
C ALGBRC(65) is Navg_G1P in component correction_factors (dimensionless).
C ALGBRC(66) is dNavgG1PdH in component correction_factors (per_molar).
C ALGBRC(67) is dNavgG1Pdmg in component correction_factors (per_molar).
C ALGBRC(68) is dmgG1Pdmg in component correction_factors (per_molar).
C ALGBRC(69) is dmgG1PdpH in component correction_factors (dimensionless).
C CONSTS(198) is NH_UG1P in component correction_factors (dimensionless).
C CONSTS(199) is deltaGof_UG1P in component correction_factors (kilojoule_per_mole).
C ALGBRC(70) is deltaGpof_UG1P in component correction_factors (kilojoule_per_mole).
C CONSTS(200) is pKa1_G6P in component correction_factors (dimensionless).
C ALGBRC(71) is P_G6P in component correction_factors (dimensionless).
C ALGBRC(72) is UG6P in component correction_factors (dimensionless).
C ALGBRC(73) is HG6P in component correction_factors (dimensionless).
C ALGBRC(74) is Navg_G6P in component correction_factors (dimensionless).
C ALGBRC(75) is dNavgG6PdH in component correction_factors (per_molar).
C CONSTS(28) is dNavgG6Pdmg in component correction_factors (per_molar).
C CONSTS(201) is NH_UG6P in component correction_factors (dimensionless).
C CONSTS(202) is deltaGof_UG6P in component correction_factors (kilojoule_per_mole).
C ALGBRC(76) is deltaGpof_UG6P in component correction_factors (kilojoule_per_mole).
C CONSTS(203) is pKa1_F6P in component correction_factors (dimensionless).
C ALGBRC(77) is P_F6P in component correction_factors (dimensionless).
C ALGBRC(78) is UF6P in component correction_factors (dimensionless).
C ALGBRC(79) is HF6P in component correction_factors (dimensionless).
C ALGBRC(80) is Navg_F6P in component correction_factors (dimensionless).
C ALGBRC(81) is dNavgF6PdH in component correction_factors (per_molar).
C CONSTS(29) is dNavgF6Pdmg in component correction_factors (per_molar).
C CONSTS(204) is NH_UF6P in component correction_factors (dimensionless).
C CONSTS(205) is deltaGof_UF6P in component correction_factors (kilojoule_per_mole).
C ALGBRC(82) is deltaGpof_UF6P in component correction_factors (kilojoule_per_mole).
C CONSTS(206) is pKa1_FDP in component correction_factors (dimensionless).
C CONSTS(207) is pKa2_FDP in component correction_factors (dimensionless).
C CONSTS(208) is pKamg_FDP in component correction_factors (dimensionless).
C ALGBRC(83) is P_FDP in component correction_factors (dimensionless).
C ALGBRC(84) is UFDP in component correction_factors (dimensionless).
C ALGBRC(85) is HFDP in component correction_factors (dimensionless).
C ALGBRC(86) is H2FDP in component correction_factors (dimensionless).
C ALGBRC(87) is mgFDP in component correction_factors (dimensionless).
C ALGBRC(88) is Navg_FDP in component correction_factors (dimensionless).
C ALGBRC(89) is dNavgFDPdH in component correction_factors (per_molar).
C ALGBRC(90) is dNavgFDPdmg in component correction_factors (per_molar).
C ALGBRC(91) is dmgFDPdmg in component correction_factors (per_molar).
C ALGBRC(92) is dmgFDPdpH in component correction_factors (dimensionless).
C CONSTS(209) is NH_UFDP in component correction_factors (dimensionless).
C CONSTS(210) is deltaGof_UFDP in component correction_factors (kilojoule_per_mole).
C ALGBRC(93) is deltaGpof_UFDP in component correction_factors (kilojoule_per_mole).
C CONSTS(211) is pKa1_GAP in component correction_factors (dimensionless).
C ALGBRC(94) is P_GAP in component correction_factors (dimensionless).
C ALGBRC(95) is UGAP in component correction_factors (dimensionless).
C ALGBRC(96) is HGAP in component correction_factors (dimensionless).
C ALGBRC(97) is Navg_GAP in component correction_factors (dimensionless).
C ALGBRC(98) is dNavgGAPdH in component correction_factors (per_molar).
C CONSTS(30) is dNavgGAPdmg in component correction_factors (per_molar).
C CONSTS(212) is NH_UGAP in component correction_factors (dimensionless).
C CONSTS(213) is deltaGof_UGAP in component correction_factors (kilojoule_per_mole).
C ALGBRC(99) is deltaGpof_UGAP in component correction_factors (kilojoule_per_mole).
C CONSTS(214) is pKamg_G3P in component correction_factors (dimensionless).
C CONSTS(215) is deltaH1o_G3P in component correction_factors (kilojoule_per_mole).
C CONSTS(216) is deltaH1_G3P in component correction_factors (kilojoule_per_mole).
C CONSTS(217) is pKa1_G3P in component correction_factors (dimensionless).
C ALGBRC(100) is P_G3P in component correction_factors (dimensionless).
C ALGBRC(101) is UG3P in component correction_factors (dimensionless).
C ALGBRC(102) is HG3P in component correction_factors (dimensionless).
C ALGBRC(103) is mgG3P in component correction_factors (dimensionless).
C ALGBRC(104) is Navg_G3P in component correction_factors (dimensionless).
C ALGBRC(105) is dNavgG3PdH in component correction_factors (per_molar).
C ALGBRC(106) is dNavgG3Pdmg in component correction_factors (per_molar).
C ALGBRC(107) is dmgG3Pdmg in component correction_factors (per_molar).
C ALGBRC(108) is dmgG3PdpH in component correction_factors (dimensionless).
C CONSTS(218) is NH_UG3P in component correction_factors (dimensionless).
C CONSTS(219) is deltaGof_UG3P in component correction_factors (kilojoule_per_mole).
C ALGBRC(109) is deltaGpof_UG3P in component correction_factors (kilojoule_per_mole).
C CONSTS(220) is pKa1_DHAP in component correction_factors (dimensionless).
C CONSTS(221) is pKamg_DHAP in component correction_factors (dimensionless).
C ALGBRC(110) is P_DHAP in component correction_factors (dimensionless).
C ALGBRC(111) is UDHAP in component correction_factors (dimensionless).
C ALGBRC(112) is HDHAP in component correction_factors (dimensionless).
C ALGBRC(113) is mgDHAP in component correction_factors (dimensionless).
C ALGBRC(114) is Navg_DHAP in component correction_factors (dimensionless).
C ALGBRC(115) is dNavgDHAPdH in component correction_factors (per_molar).
C ALGBRC(116) is dNavgDHAPdmg in component correction_factors (per_molar).
C ALGBRC(117) is dmgDHAPdmg in component correction_factors (per_molar).
C ALGBRC(118) is dmgDHAPdpH in component correction_factors (dimensionless).
C CONSTS(222) is NH_UDHAP in component correction_factors (dimensionless).
C CONSTS(223) is deltaGof_UDHAP in component correction_factors (kilojoule_per_mole).
C ALGBRC(119) is deltaGpof_UDHAP in component correction_factors (kilojoule_per_mole).
C CONSTS(224) is pKa1_13DPG in component correction_factors (dimensionless).
C ALGBRC(120) is P_13DPG in component correction_factors (dimensionless).
C ALGBRC(121) is U13DPG in component correction_factors (dimensionless).
C ALGBRC(122) is H13DPG in component correction_factors (dimensionless).
C ALGBRC(123) is Navg_13DPG in component correction_factors (dimensionless).
C ALGBRC(124) is dNavg13DPGdH in component correction_factors (per_molar).
C CONSTS(31) is dNavg13DPGdmg in component correction_factors (per_molar).
C CONSTS(225) is NH_U13DPG in component correction_factors (dimensionless).
C CONSTS(226) is deltaGof_U13DPG in component correction_factors (kilojoule_per_mole).
C ALGBRC(125) is deltaGpof_U13DPG in component correction_factors (kilojoule_per_mole).
C CONSTS(227) is pKa1_3PG in component correction_factors (dimensionless).
C ALGBRC(126) is P_3PG in component correction_factors (dimensionless).
C ALGBRC(127) is U3PG in component correction_factors (dimensionless).
C ALGBRC(128) is H3PG in component correction_factors (dimensionless).
C ALGBRC(129) is Navg_3PG in component correction_factors (dimensionless).
C ALGBRC(130) is dNavg3PGdH in component correction_factors (per_molar).
C CONSTS(32) is dNavg3PGdmg in component correction_factors (per_molar).
C CONSTS(228) is NH_U3PG in component correction_factors (dimensionless).
C CONSTS(229) is deltaGof_U3PG in component correction_factors (kilojoule_per_mole).
C ALGBRC(131) is deltaGpof_U3PG in component correction_factors (kilojoule_per_mole).
C CONSTS(230) is pKa1_2PG in component correction_factors (dimensionless).
C CONSTS(231) is pKamg_2PG in component correction_factors (dimensionless).
C CONSTS(232) is pKak_2PG in component correction_factors (dimensionless).
C ALGBRC(132) is P_2PG in component correction_factors (dimensionless).
C ALGBRC(133) is U2PG in component correction_factors (dimensionless).
C ALGBRC(134) is H2PG in component correction_factors (dimensionless).
C ALGBRC(136) is mg2PG in component correction_factors (dimensionless).
C ALGBRC(135) is k2PG in component correction_factors (dimensionless).
C ALGBRC(137) is Navg_2PG in component correction_factors (dimensionless).
C ALGBRC(138) is dNavg2PGdH in component correction_factors (per_molar).
C ALGBRC(139) is dNavg2PGdmg in component correction_factors (per_molar).
C ALGBRC(140) is dmg2PGdmg in component correction_factors (per_molar).
C ALGBRC(141) is dmg2PGdpH in component correction_factors (dimensionless).
C CONSTS(233) is NH_U2PG in component correction_factors (dimensionless).
C CONSTS(234) is deltaGof_U2PG in component correction_factors (kilojoule_per_mole).
C ALGBRC(142) is deltaGpof_U2PG in component correction_factors (kilojoule_per_mole).
C CONSTS(235) is pKa1_PEP in component correction_factors (dimensionless).
C CONSTS(236) is pKamg_PEP in component correction_factors (dimensionless).
C CONSTS(237) is pKak_PEP in component correction_factors (dimensionless).
C ALGBRC(143) is P_PEP in component correction_factors (dimensionless).
C ALGBRC(144) is UPEP in component correction_factors (dimensionless).
C ALGBRC(145) is HPEP in component correction_factors (dimensionless).
C ALGBRC(146) is kPEP in component correction_factors (dimensionless).
C ALGBRC(147) is mgPEP in component correction_factors (dimensionless).
C ALGBRC(148) is Navg_PEP in component correction_factors (dimensionless).
C ALGBRC(149) is dNavgPEPdH in component correction_factors (per_molar).
C ALGBRC(150) is dNavgPEPdmg in component correction_factors (per_molar).
C ALGBRC(151) is dmgPEPdmg in component correction_factors (per_molar).
C ALGBRC(152) is dmgPEPdpH in component correction_factors (dimensionless).
C CONSTS(238) is NH_UPEP in component correction_factors (dimensionless).
C CONSTS(239) is deltaGof_UPEP in component correction_factors (kilojoule_per_mole).
C ALGBRC(153) is deltaGpof_UPEP in component correction_factors (kilojoule_per_mole).
C CONSTS(240) is pKa1_PYR in component correction_factors (dimensionless).
C ALGBRC(154) is P_PYR in component correction_factors (dimensionless).
C ALGBRC(155) is UPYR in component correction_factors (dimensionless).
C ALGBRC(156) is HPYR in component correction_factors (dimensionless).
C ALGBRC(157) is Navg_PYR in component correction_factors (dimensionless).
C ALGBRC(158) is dNavgPYRdH in component correction_factors (per_molar).
C CONSTS(33) is dNavgPYRdmg in component correction_factors (per_molar).
C CONSTS(241) is NH_UPYR in component correction_factors (dimensionless).
C CONSTS(242) is deltaGof_UPYR in component correction_factors (kilojoule_per_mole).
C ALGBRC(159) is deltaGpof_UPYR in component correction_factors (kilojoule_per_mole).
C CONSTS(243) is pKamg_LAC in component correction_factors (dimensionless).
C CONSTS(244) is deltaH1o_LAC in component correction_factors (kilojoule_per_mole).
C CONSTS(245) is deltaH1_LAC in component correction_factors (kilojoule_per_mole).
C CONSTS(246) is pKa1_LAC in component correction_factors (dimensionless).
C ALGBRC(160) is P_LAC in component correction_factors (dimensionless).
C ALGBRC(161) is ULAC in component correction_factors (dimensionless).
C ALGBRC(162) is HLAC in component correction_factors (dimensionless).
C ALGBRC(163) is mgLAC in component correction_factors (dimensionless).
C ALGBRC(164) is Navg_LAC in component correction_factors (dimensionless).
C ALGBRC(165) is dNavgLACdH in component correction_factors (per_molar).
C ALGBRC(166) is dNavgLACdmg in component correction_factors (per_molar).
C ALGBRC(167) is dmgLACdmg in component correction_factors (per_molar).
C ALGBRC(168) is dmgLACdpH in component correction_factors (dimensionless).
C CONSTS(247) is NH_ULAC in component correction_factors (dimensionless).
C CONSTS(248) is deltaGof_ULAC in component correction_factors (kilojoule_per_mole).
C ALGBRC(169) is deltaGpof_ULAC in component correction_factors (kilojoule_per_mole).
C CONSTS(249) is dNH_GLY in component correction_factors (dimensionless).
C ALGBRC(170) is deltaGpo_GLY in component correction_factors (kilojoule_per_mole).
C CONSTS(250) is NH_NAD in component correction_factors (dimensionless).
C CONSTS(251) is deltaGof_NAD in component correction_factors (kilojoule_per_mole).
C ALGBRC(171) is deltaGpof_NAD in component correction_factors (kilojoule_per_mole).
C CONSTS(252) is NH_NADH in component correction_factors (dimensionless).
C CONSTS(253) is deltaGof_NADH in component correction_factors (kilojoule_per_mole).
C ALGBRC(172) is deltaGpof_NADH in component correction_factors (kilojoule_per_mole).
C CONSTS(254) is NH_H2O in component correction_factors (dimensionless).
C CONSTS(255) is deltaGof_H2O in component correction_factors (kilojoule_per_mole).
C ALGBRC(173) is deltaGpof_H2O in component correction_factors (kilojoule_per_mole).
C CONSTS(256) is NH_H in component correction_factors (dimensionless).
C CONSTS(257) is deltaGof_H in component correction_factors (kilojoule_per_mole).
C ALGBRC(174) is deltaGpof_H in component correction_factors (kilojoule_per_mole).
C ALGBRC(175) is deltaH_CK in component correction_factors (dimensionless).
C CONSTS(258) is Kref_CK in component correction_factors (dimensionless).
C CONSTS(259) is deltaHo_CKo in component correction_factors (kilojoule_per_mole).
C CONSTS(260) is deltaH1_CK in component correction_factors (kilojoule_per_mole).
C CONSTS(261) is Kref_CKI in component correction_factors (dimensionless).
C CONSTS(262) is Kref_CKT in component correction_factors (dimensionless).
C CONSTS(263) is deltaGpo_CK in component correction_factors (kilojoule_per_mole).
C ALGBRC(177) is Kapp_CK in component correction_factors (dimensionless).
C ALGBRC(179) is deltaH_ADK in component correction_factors (dimensionless).
C ALGBRC(180) is deltaGpo_ADK in component correction_factors (kilojoule_per_mole).
C ALGBRC(181) is Kapp_ADK in component correction_factors (dimensionless).
C ALGBRC(182) is deltaH_GP in component correction_factors (dimensionless).
C ALGBRC(183) is deltaGpo_GP in component correction_factors (kilojoule_per_mole).
C ALGBRC(184) is Kapp_GP in component correction_factors (dimensionless).
C ALGBRC(185) is deltaH_PGLM in component correction_factors (dimensionless).
C ALGBRC(186) is deltaGpo_PGLM in component correction_factors (kilojoule_per_mole).
C ALGBRC(187) is Kapp_PGLM in component correction_factors (dimensionless).
C ALGBRC(188) is deltaH_PGI in component correction_factors (dimensionless).
C ALGBRC(189) is deltaGpo_PGI in component correction_factors (kilojoule_per_mole).
C ALGBRC(190) is Kapp_PGI in component correction_factors (dimensionless).
C ALGBRC(191) is deltaH_PFK in component correction_factors (dimensionless).
C ALGBRC(192) is deltaGpo_PFK in component correction_factors (kilojoule_per_mole).
C ALGBRC(193) is Kapp_PFK in component correction_factors (dimensionless).
C ALGBRC(194) is deltaH_ALD in component correction_factors (dimensionless).
C ALGBRC(195) is deltaGpo_ALD in component correction_factors (kilojoule_per_mole).
C ALGBRC(196) is Kapp_ALD in component correction_factors (molar).
C ALGBRC(197) is deltaH_TPI in component correction_factors (dimensionless).
C ALGBRC(198) is deltaGpo_TPI in component correction_factors (kilojoule_per_mole).
C ALGBRC(199) is Kapp_TPI in component correction_factors (dimensionless).
C ALGBRC(200) is deltaH_GAPDH in component correction_factors (dimensionless).
C ALGBRC(201) is deltaGpo_GAPDH in component correction_factors (kilojoule_per_mole).
C ALGBRC(202) is Kapp_GAPDH in component correction_factors (per_molar).
C ALGBRC(203) is deltaH_G3PDH in component correction_factors (dimensionless).
C ALGBRC(204) is deltaGpo_G3PDH in component correction_factors (kilojoule_per_mole).
C ALGBRC(205) is Kapp_G3PDH in component correction_factors (dimensionless).
C ALGBRC(206) is deltaH_PGK in component correction_factors (dimensionless).
C ALGBRC(207) is deltaGpo_PGK in component correction_factors (kilojoule_per_mole).
C ALGBRC(208) is Kapp_PGK in component correction_factors (dimensionless).
C ALGBRC(209) is deltaH_PGM in component correction_factors (dimensionless).
C ALGBRC(210) is deltaGpo_PGM in component correction_factors (kilojoule_per_mole).
C ALGBRC(211) is Kapp_PGM in component correction_factors (dimensionless).
C ALGBRC(212) is deltaH_ENOL in component correction_factors (dimensionless).
C ALGBRC(213) is deltaGpo_ENOL in component correction_factors (kilojoule_per_mole).
C ALGBRC(214) is Kapp_ENOL in component correction_factors (dimensionless).
C ALGBRC(215) is deltaH_PK in component correction_factors (dimensionless).
C ALGBRC(216) is deltaGpo_PK in component correction_factors (kilojoule_per_mole).
C ALGBRC(217) is Kapp_PK in component correction_factors (dimensionless).
C ALGBRC(218) is deltaH_LDH in component correction_factors (dimensionless).
C ALGBRC(219) is deltaGpo_LDH in component correction_factors (kilojoule_per_mole).
C ALGBRC(220) is Kapp_LDH in component correction_factors (dimensionless).
C ALGBRC(31) is deltaH_ATPase in component correction_factors (dimensionless).
C ALGBRC(176) is deltaGpo_ATPase in component correction_factors (kilojoule_per_mole).
C ALGBRC(178) is Kapp_ATPase in component correction_factors (dimensionless).
C STATES(2) is Mg in component differential_equations (molar).
C CONSTS(264) is Vfgly in component glycogen_phosphorylase (molar_per_minute).
C CONSTS(39) is expno in component glycogen_phosphorylase (dimensionless).
C ALGBRC(221) is fracA in component glycogen_phosphorylase (dimensionless).
C CONSTS(40) is KgpA_glyf in component glycogen_phosphorylase (molar).
C CONSTS(265) is KgpA_pi in component glycogen_phosphorylase (molar).
C CONSTS(266) is KgpA_igly in component glycogen_phosphorylase (molar).
C CONSTS(41) is KgpA_ipi in component glycogen_phosphorylase (molar).
C CONSTS(267) is KgpA_glyb in component glycogen_phosphorylase (molar).
C CONSTS(42) is KgpA_g1p in component glycogen_phosphorylase (molar).
C CONSTS(268) is KgpA_ig1p in component glycogen_phosphorylase (molar).
C ALGBRC(222) is Dglya in component glycogen_phosphorylase (dimensionless).
C ALGBRC(223) is pa in component glycogen_phosphorylase (dimensionless).
C ALGBRC(224) is VbglyA in component glycogen_phosphorylase (molar_per_minute).
C ALGBRC(225) is glyAF in component glycogen_phosphorylase (per_minute).
C ALGBRC(226) is glyAR in component glycogen_phosphorylase (per_minute).
C ALGBRC(227) is flux_GPa in component glycogen_phosphorylase (molar_per_minute).
C STATES(3) is G1P in component differential_equations (molar).
C STATES(4) is Pi in component differential_equations (molar).
C STATES(5) is Gly in component differential_equations (molar).
C ALGBRC(228) is fracB in component glycogen_phosphorylase_B (dimensionless).
C CONSTS(269) is KgpB_pi in component glycogen_phosphorylase_B (molar).
C CONSTS(43) is KgpB_ipi in component glycogen_phosphorylase_B (molar).
C CONSTS(270) is KgpB_iglyf in component glycogen_phosphorylase_B (molar).
C CONSTS(271) is KgpB_g1p in component glycogen_phosphorylase_B (molar).
C CONSTS(44) is KgpB_ig1p in component glycogen_phosphorylase_B (molar).
C CONSTS(272) is KgpB_iglyb in component glycogen_phosphorylase_B (molar).
C CONSTS(45) is Kgp_amp in component glycogen_phosphorylase_B (molar).
C CONSTS(46) is interactioncoeff in component glycogen_phosphorylase_B (dimensionless).
C CONSTS(47) is nH in component glycogen_phosphorylase_B (dimensionless).
C ALGBRC(229) is M in component glycogen_phosphorylase_B (dimensionless).
C ALGBRC(230) is Dglyb in component glycogen_phosphorylase_B (dimensionless).
C ALGBRC(231) is pb in component glycogen_phosphorylase_B (dimensionless).
C ALGBRC(232) is VbglyB in component glycogen_phosphorylase_B (molar_per_minute).
C ALGBRC(233) is glyBF in component glycogen_phosphorylase_B (per_minute).
C ALGBRC(234) is glyBR in component glycogen_phosphorylase_B (per_minute).
C ALGBRC(235) is flux_GPb in component glycogen_phosphorylase_B (molar_per_minute).
C STATES(6) is AMP in component differential_equations (molar).
C CONSTS(273) is Vffpglm in component PGLM (molar_per_minute).
C CONSTS(274) is Kpglm_g1p in component PGLM (molar).
C CONSTS(275) is Kpglm_g6p in component PGLM (molar).
C ALGBRC(236) is Vfpglm in component PGLM (molar_per_minute).
C ALGBRC(238) is Vbpglm in component PGLM (molar_per_minute).
C ALGBRC(239) is v_PGLM in component PGLM (molar_per_minute).
C STATES(7) is G6P in component differential_equations (molar).
C CONSTS(276) is Vbbpgi in component PGI (molar_per_minute).
C CONSTS(277) is Kpgi_g6p in component PGI (molar).
C CONSTS(278) is Kpgi_f6p in component PGI (molar).
C ALGBRC(240) is Vbpgi in component PGI (molar_per_minute).
C ALGBRC(241) is Vfpgi in component PGI (molar_per_minute).
C ALGBRC(242) is v_PGI in component PGI (molar_per_minute).
C STATES(8) is F6P in component differential_equations (molar).
C CONSTS(279) is Vffpfk in component PFK (molar_per_minute).
C CONSTS(280) is Kpfk_f6p in component PFK (molar).
C CONSTS(281) is Kpfk_f6pT in component PFK (molar).
C CONSTS(282) is Kpfk_atp in component PFK (molar).
C CONSTS(284) is Kpfk_atpT in component PFK (molar).
C CONSTS(283) is Kpfk_fbp in component PFK (molar).
C CONSTS(48) is Kpfk_fbpT in component PFK (molar).
C CONSTS(285) is Kpfk_adp in component PFK (molar).
C CONSTS(49) is Kpfk_adpT in component PFK (molar).
C CONSTS(50) is Kpfki in component PFK (molar).
C CONSTS(51) is Kmpfk in component PFK (molar).
C CONSTS(52) is d in component PFK (dimensionless).
C CONSTS(53) is e_ in component PFK (dimensionless).
C CONSTS(54) is Lo in component PFK (dimensionless).
C ALGBRC(243) is Vfpfk in component PFK (molar_per_minute).
C ALGBRC(245) is Vbpfk in component PFK (molar_per_minute).
C ALGBRC(246) is L in component PFK (dimensionless).
C CONSTS(286) is alpha in component PFK (dimensionless).
C ALGBRC(247) is Delta in component PFK (dimensionless).
C ALGBRC(248) is Deltap in component PFK (dimensionless).
C ALGBRC(249) is v_PFK in component PFK (molar_per_minute).
C STATES(9) is FBP in component differential_equations (molar).
C STATES(10) is ADP in component differential_equations (molar).
C STATES(11) is ATP in component differential_equations (molar).
C CONSTS(287) is Vffald in component ALD (molar_per_minute).
C CONSTS(288) is Kald_fbp in component ALD (molar).
C CONSTS(289) is Kald_dhap in component ALD (molar).
C CONSTS(290) is Kald_gap in component ALD (molar).
C ALGBRC(250) is Vfald in component ALD (molar_per_minute).
C ALGBRC(252) is Vbald in component ALD (molar_per_minute).
C ALGBRC(253) is v_ALD in component ALD (molar_per_minute).
C STATES(12) is DHAP in component differential_equations (molar).
C STATES(13) is GAP in component differential_equations (molar).
C CONSTS(291) is Vfftpi in component TPI (molar_per_minute).
C CONSTS(292) is Ktpi_gap in component TPI (molar).
C CONSTS(293) is Ktpi_dhap in component TPI (molar).
C CONSTS(294) is Vftpi in component TPI (molar_per_minute).
C ALGBRC(254) is Vbtpi in component TPI (molar_per_minute).
C ALGBRC(255) is v_TPI in component TPI (molar_per_minute).
C CONSTS(295) is Vbbg3pdh in component G3PDH (molar_per_minute).
C CONSTS(296) is Kg3pdh_g3p in component G3PDH (molar).
C CONSTS(297) is Kg3pdh_nad in component G3PDH (molar).
C CONSTS(298) is Kg3pdh_dhap in component G3PDH (molar).
C CONSTS(299) is Kg3pdh_nadh in component G3PDH (molar).
C ALGBRC(256) is Dg3pdh in component G3PDH (dimensionless).
C CONSTS(300) is Vbg3pdh in component G3PDH (molar_per_minute).
C ALGBRC(257) is Vfg3pdh in component G3PDH (molar_per_minute).
C ALGBRC(258) is v_G3PDH in component G3PDH (molar_per_minute).
C STATES(14) is G3P in component differential_equations (molar).
C STATES(15) is NAD in component differential_equations (molar).
C STATES(16) is NADH in component differential_equations (molar).
C CONSTS(301) is Vffgad in component GAPDH (molar_per_minute).
C CONSTS(302) is Kgapdh_gap in component GAPDH (molar).
C CONSTS(303) is Kgapdh_nad in component GAPDH (molar).
C CONSTS(304) is Kgapdh_pi in component GAPDH (molar).
C CONSTS(305) is Kgapdh_bpg in component GAPDH (molar).
C CONSTS(306) is Kgapdh_nadh in component GAPDH (molar).
C ALGBRC(259) is Dgap in component GAPDH (dimensionless).
C ALGBRC(260) is Vfgad in component GAPDH (molar_per_minute).
C ALGBRC(261) is Vbgad in component GAPDH (molar_per_minute).
C ALGBRC(262) is v_GAPDH in component GAPDH (molar_per_minute).
C STATES(17) is BPG in component differential_equations (molar).
C CONSTS(307) is Vbbpgk in component PGK (molar_per_minute).
C CONSTS(308) is Kpgk_bpg in component PGK (molar).
C CONSTS(309) is Kpgk_adp in component PGK (molar).
C CONSTS(310) is Kpgk_3pg in component PGK (molar).
C CONSTS(311) is Kpgk_atp in component PGK (molar).
C CONSTS(312) is Vbpgk in component PGK (molar_per_minute).
C ALGBRC(263) is Vfpgk in component PGK (molar_per_minute).
C ALGBRC(265) is D_PGK in component PGK (dimensionless).
C ALGBRC(266) is v_PGK in component PGK (molar_per_minute).
C STATES(18) is P3G in component differential_equations (molar).
C CONSTS(313) is Vffpgm in component PGM (molar_per_minute).
C CONSTS(314) is Kpgm_3pg in component PGM (molar).
C CONSTS(315) is Kpgm_2pg in component PGM (molar).
C ALGBRC(267) is Vfpgm in component PGM (molar_per_minute).
C ALGBRC(269) is Vbpgm in component PGM (molar_per_minute).
C ALGBRC(270) is v_PGM in component PGM (molar_per_minute).
C STATES(19) is P2G in component differential_equations (molar).
C CONSTS(316) is Vffen in component ENOL (molar_per_minute).
C CONSTS(317) is Ken_2pg in component ENOL (molar).
C CONSTS(318) is Ken_pep in component ENOL (molar).
C CONSTS(319) is Vfen in component ENOL (molar_per_minute).
C ALGBRC(271) is Vben in component ENOL (molar_per_minute).
C ALGBRC(273) is v_ENOL in component ENOL (molar_per_minute).
C STATES(20) is PEP in component differential_equations (molar).
C CONSTS(320) is Vffpk in component PK (molar_per_minute).
C CONSTS(321) is Kpk_pep in component PK (molar).
C CONSTS(322) is Kpk_adp in component PK (molar).
C CONSTS(323) is Kpk_pyr in component PK (molar).
C CONSTS(324) is Kpk_atp in component PK (molar).
C ALGBRC(274) is Vfpk in component PK (molar_per_minute).
C ALGBRC(275) is Vbpk in component PK (molar_per_minute).
C ALGBRC(276) is v_PK in component PK (molar_per_minute).
C STATES(21) is PYR in component differential_equations (molar).
C CONSTS(325) is Vffldh in component LDH (molar_per_minute).
C CONSTS(326) is Kldh_pyr in component LDH (molar).
C CONSTS(327) is Kldh_nadh in component LDH (molar).
C CONSTS(328) is Kldh_lac in component LDH (molar).
C CONSTS(329) is Kldh_nad in component LDH (molar).
C ALGBRC(277) is Vfldh in component LDH (molar_per_minute).
C ALGBRC(278) is Vbldh in component LDH (molar_per_minute).
C ALGBRC(279) is v_LDH in component LDH (molar_per_minute).
C STATES(22) is LAC in component differential_equations (molar).
C ALGBRC(280) is VmaxATPase in component ATPase (molar_per_minute).
C CONSTS(55) is Katp_ATPase in component ATPase (molar).
C ALGBRC(285) is ATPase in component ATPase (molar_per_minute).
C CONSTS(330) is VforCK in component creatine_kinase (molar_per_minute).
C CONSTS(331) is Kck_pcr in component creatine_kinase (molar).
C CONSTS(332) is Kck_iatp in component creatine_kinase (molar).
C CONSTS(333) is Kck_iadp in component creatine_kinase (molar).
C CONSTS(56) is Kck_ipcr in component creatine_kinase (molar).
C CONSTS(334) is Kck_cr in component creatine_kinase (molar).
C ALGBRC(286) is VrevCK in component creatine_kinase (molar_per_minute).
C ALGBRC(287) is CK in component creatine_kinase (molar_per_minute).
C STATES(23) is Cr in component differential_equations (molar).
C STATES(24) is PCr in component differential_equations (molar).
C CONSTS(335) is Vfadk in component adenylate_kinase (molar_per_minute).
C CONSTS(336) is Kadk_amp in component adenylate_kinase (molar).
C CONSTS(337) is Kadk_atp in component adenylate_kinase (molar).
C CONSTS(338) is Kadk_adp in component adenylate_kinase (molar).
C ALGBRC(288) is Vbadk in component adenylate_kinase (molar_per_minute).
C ALGBRC(291) is ADK in component adenylate_kinase (molar_per_minute).
C CONSTS(339) is carnosine in component buffer_capacity (molar).
C CONSTS(340) is tris in component buffer_capacity (molar).
C CONSTS(341) is acetate in component buffer_capacity (molar).
C ALGBRC(292) is bufcapfixed in component buffer_capacity (molar).
C ALGBRC(293) is bufcapmetab in component buffer_capacity (molar).
C ALGBRC(294) is protons_consumed in component buffer_capacity (molar_per_minute).
C ALGBRC(289) is CKprtflux in component buffer_capacity (molar_per_minute).
C ALGBRC(281) is glycprtflux in component buffer_capacity (molar_per_minute).
C ALGBRC(295) is pHODEterm1 in component buffer_capacity (per_minute).
C ALGBRC(296) is pHODEterm2 in component buffer_capacity (per_molar).
C ALGBRC(297) is denom_mgODE in component buffer_capacity (dimensionless).
C ALGBRC(298) is RHSterm1_mgODE in component buffer_capacity (molar).
C ALGBRC(299) is denomMgpHODE in component buffer_capacity (dimensionless).
C ALGBRC(312) is RHSterm2_mgODE in component differential_equations (molar_per_minute).
C CONSTS(34) is fixmg in component differential_equations (dimensionless).
C CONSTS(342) is fixpH in component differential_equations (dimensionless).
C STATES(25) is protonload in component differential_equations (molar).
C ALGBRC(301) is dATPdt in component differential_equations (molar_per_minute).
C ALGBRC(302) is dADPdt in component differential_equations (molar_per_minute).
C ALGBRC(303) is dAMPdt in component differential_equations (molar_per_minute).
C ALGBRC(307) is dDHAPdt in component differential_equations (molar_per_minute).
C ALGBRC(306) is dFBPdt in component differential_equations (molar_per_minute).
C ALGBRC(305) is dG1Pdt in component differential_equations (molar_per_minute).
C ALGBRC(308) is dG3Pdt in component differential_equations (molar_per_minute).
C ALGBRC(311) is dLACdt in component differential_equations (molar_per_minute).
C ALGBRC(309) is dP2Gdt in component differential_equations (molar_per_minute).
C ALGBRC(300) is dPCrdt in component differential_equations (molar_per_minute).
C ALGBRC(310) is dPEPdt in component differential_equations (molar_per_minute).
C ALGBRC(304) is dPidt in component differential_equations (molar_per_minute).
C ALGBRC(290) is dCrdt in component differential_equations (molar_per_minute).
C ALGBRC(282) is dNADdt in component differential_equations (molar_per_minute).
C ALGBRC(283) is dNADHdt in component differential_equations (molar_per_minute).
C ALGBRC(237) is dGlydt in component differential_equations (molar_per_minute).
C ALGBRC(244) is dG6Pdt in component differential_equations (molar_per_minute).
C ALGBRC(251) is dF6Pdt in component differential_equations (molar_per_minute).
C ALGBRC(264) is dGAPdt in component differential_equations (molar_per_minute).
C ALGBRC(268) is dBGPdt in component differential_equations (molar_per_minute).
C ALGBRC(272) is dP3Gdt in component differential_equations (molar_per_minute).
C ALGBRC(284) is dPYRdt in component differential_equations (molar_per_minute).
C RATES(24) is d/dt PCr in component differential_equations (molar).
C RATES(23) is d/dt Cr in component differential_equations (molar).
C RATES(15) is d/dt NAD in component differential_equations (molar).
C RATES(16) is d/dt NADH in component differential_equations (molar).
C RATES(11) is d/dt ATP in component differential_equations (molar).
C RATES(10) is d/dt ADP in component differential_equations (molar).
C RATES(6) is d/dt AMP in component differential_equations (molar).
C RATES(4) is d/dt Pi in component differential_equations (molar).
C RATES(5) is d/dt Gly in component differential_equations (molar).
C RATES(3) is d/dt G1P in component differential_equations (molar).
C RATES(7) is d/dt G6P in component differential_equations (molar).
C RATES(8) is d/dt F6P in component differential_equations (molar).
C RATES(9) is d/dt FBP in component differential_equations (molar).
C RATES(12) is d/dt DHAP in component differential_equations (molar).
C RATES(14) is d/dt G3P in component differential_equations (molar).
C RATES(13) is d/dt GAP in component differential_equations (molar).
C RATES(17) is d/dt BPG in component differential_equations (molar).
C RATES(18) is d/dt P3G in component differential_equations (molar).
C RATES(19) is d/dt P2G in component differential_equations (molar).
C RATES(20) is d/dt PEP in component differential_equations (molar).
C RATES(21) is d/dt PYR in component differential_equations (molar).
C RATES(22) is d/dt LAC in component differential_equations (molar).
C RATES(2) is d/dt Mg in component differential_equations (molar).
C RATES(1) is d/dt pH_calc in component differential_equations (dimensionless).
C RATES(25) is d/dt protonload in component differential_equations (molar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = 7.8
      CONSTS(1) = 8.314e-3
      CONSTS(2) = 298.15
      STATES(2) = 5.132658807e-4
      STATES(3) = 1e-9
      STATES(4) = 0.03
      STATES(5) = 0.04
      STATES(6) = 1e-9
      STATES(7) = 1e-9
      STATES(8) = 1e-9
      STATES(9) = 1e-9
      STATES(10) = 1e-9
      STATES(11) = 0.005
      STATES(12) = 1e-9
      STATES(13) = 1e-9
      STATES(14) = 1e-9
      STATES(15) = 0.0005
      STATES(16) = 1e-9
      STATES(17) = 1e-9
      STATES(18) = 1e-9
      STATES(19) = 1e-9
      STATES(20) = 1e-9
      STATES(21) = 1e-9
      STATES(22) = 1e-9
      STATES(23) = 0.029999999
      STATES(24) = 1e-9
      STATES(25) = 0
      CONSTS(3) = 1.00000
      CONSTS(4) = 7.40000
      CONSTS(5) = 0.0500000
      CONSTS(6) = 0.400000
      CONSTS(7) = 0.00170000
      CONSTS(8) = 0.00470000
      CONSTS(9) = 0.00270000
      CONSTS(10) = 0.00460000
      CONSTS(11) = 0.00740000
      CONSTS(12) = 0.00266055
      CONSTS(13) = 0.0200000
      CONSTS(14) = 1.75052
      CONSTS(15) = 0.00402000
      CONSTS(16) = 0.00270000
      CONSTS(17) = 0.000870000
      CONSTS(18) = 6.00000e-05
      CONSTS(19) = 0.0100000
      CONSTS(20) = 0.0100000
      CONSTS(21) = 13.0000
      CONSTS(22) = 0.00390000
      CONSTS(23) = 0.00000
      CONSTS(24) = 0.000100000
      CONSTS(25) = 0.00500000
      CONSTS(26) = 29.0000
      CONSTS(27) = 0.00000
      CONSTS(28) = 0.00000
      CONSTS(29) = 0.00000
      CONSTS(30) = 0.00000
      CONSTS(31) = 0.00000
      CONSTS(32) = 0.00000
      CONSTS(33) = 0.00000
      CONSTS(34) = 1.00000
      CONSTS(35) = CONSTS(3)
      CONSTS(36) = CONSTS(4)
      CONSTS(37) = 0.00400000
      CONSTS(38) = CONSTS(25)
      CONSTS(39) = CONSTS(26)
      CONSTS(40) = CONSTS(7)
      CONSTS(41) = CONSTS(8)
      CONSTS(42) = CONSTS(9)
      CONSTS(43) = CONSTS(10)
      CONSTS(44) = CONSTS(11)
      CONSTS(45) = CONSTS(12)
      CONSTS(46) = CONSTS(13)
      CONSTS(47) = CONSTS(14)
      CONSTS(48) = CONSTS(15)
      CONSTS(49) = CONSTS(16)
      CONSTS(50) = CONSTS(17)
      CONSTS(51) = CONSTS(18)
      CONSTS(52) = CONSTS(19)
      CONSTS(53) = CONSTS(20)
      CONSTS(54) = CONSTS(21)
      CONSTS(55) = CONSTS(24)
      CONSTS(56) = CONSTS(22)
      CONSTS(57) = 0.00200000/1.50000
      CONSTS(58) = 0.000150000
      CONSTS(59) = 0.0101000
      CONSTS(60) = 0.000200000
      CONSTS(61) = 0.0150000
      CONSTS(62) = 0.00150000
      CONSTS(63) = 0.00440000
      CONSTS(64) = 0.480000
      CONSTS(65) = 6.30000e-05
      CONSTS(66) = 3.00000e-05
      CONSTS(67) = 0.880000
      CONSTS(68) = 0.000480000
      CONSTS(69) = 0.000119000
      CONSTS(70) = 0.0560000
      CONSTS(71) = 0.000180000
      CONSTS(72) = 0.0200000
      CONSTS(73) = 8.00000e-05
      CONSTS(74) = 0.000250000
      CONSTS(75) = 0.00402000
      CONSTS(76) = 0.00270000
      CONSTS(77) = 0.0106591
      CONSTS(78) = 5.00000e-05
      CONSTS(79) = 0.00200000
      CONSTS(80) = 0.00100000
      CONSTS(81) = 12.0000
      CONSTS(82) = 0.000320000
      CONSTS(83) = 0.000610000
      CONSTS(84) = 0.0825000
      CONSTS(85) = 0.000180000
      CONSTS(86) = 1.20000e-05
      CONSTS(87) = 0.000220000
      CONSTS(88) = 8.00000e-06
      CONSTS(89) = 1.26500
      CONSTS(90) = 2.50000e-06
      CONSTS(91) = 9.00000e-05
      CONSTS(92) = 0.000290000
      CONSTS(93) = 8.00000e-07
      CONSTS(94) = 3.30000e-06
      CONSTS(95) = 1.12000
      CONSTS(96) = 0.00200000
      CONSTS(97) = 8.00000e-06
      CONSTS(98) = 0.00120000
      CONSTS(99) = 0.000350000
      CONSTS(100) = 1.12000
      CONSTS(101) = 0.000200000
      CONSTS(102) = 1.40000e-05
      CONSTS(103) = 0.192000
      CONSTS(104) = 0.000100000
      CONSTS(105) = 0.000370000
      CONSTS(106) = 1.44000
      CONSTS(107) = 8.00000e-05
      CONSTS(108) = 0.000300000
      CONSTS(109) = 0.00705000
      CONSTS(110) = 0.00113000
      CONSTS(111) = 1.92000
      CONSTS(112) = 0.000335000
      CONSTS(113) = 2.00000e-06
      CONSTS(114) = 0.0170000
      CONSTS(115) = 0.000849000
      CONSTS(116) = 0.500000
      CONSTS(117) = 0.00111000
      CONSTS(118) = 0.00350000
      CONSTS(119) = 0.000135000
      CONSTS(120) = 0.00380000
      CONSTS(121) = 0.880000
      CONSTS(122) = 0.000320000
      CONSTS(123) = 0.000270000
      CONSTS(124) = 0.000350000
      CONSTS(125) = 0.100000
      CONSTS(126) = CONSTS(125)
      CONSTS(127) = 0.0150000
      CONSTS(128) = 0.0250000
      CONSTS(129) = 0.0100000
      CONSTS(130) = 303.150
      CONSTS(131) = CONSTS(130)
      CONSTS(132) = 0.0800000
      CONSTS(133) = 1.00000
      CONSTS(134) = CONSTS(132)
      CONSTS(135) = 1.00000
      CONSTS(136) = 1.47750
      CONSTS(137) = 1.60000
      CONSTS(138) = ( CONSTS(136)* CONSTS(126) ** (1.0 / 2))/(1.00000+ CONSTS(137)* CONSTS(126) ** (1.0 / 2))
      CONSTS(139) = 0.100000
      CONSTS(140) = 1.17582
      CONSTS(141) = ( 1.00000*CONSTS(140)*( CONSTS(139) ** (1.0 / 2)/(1.00000+ CONSTS(137)* CONSTS(139) ** (1.0 / 2)) -  CONSTS(126) ** (1.0 / 2)/(1.00000+ CONSTS(137)* CONSTS(126) ** (1.0 / 2))))/log(10.0000)
      CONSTS(142) = (1.00000/CONSTS(131) - 1.00000/CONSTS(2))/( log(10.0000)*CONSTS(1))
      CONSTS(143) = 2.91482
      CONSTS(144) = ( CONSTS(143)* CONSTS(126) ** (1.0 / 2))/(1.00000+ CONSTS(137)* CONSTS(126) ** (1.0 / 2))
      CONSTS(145) = 0.500000
      CONSTS(146) = 3.00000
      CONSTS(147) = - 2.90000
      CONSTS(148) = CONSTS(146)+ CONSTS(138)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(149) = CONSTS(147)+ CONSTS(138)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(150) = 6.75000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(148)
      CONSTS(151) = 1.65000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(149)
      CONSTS(152) = 1.00000
      CONSTS(153) = - 1096.10
      CONSTS(154) = - 5.00000
      CONSTS(155) = - 18.0000
      CONSTS(156) = - 1.00000
      CONSTS(157) = CONSTS(154)+ CONSTS(138)*((4.00000 ** 2.00000+1.00000 ** 2.00000) - 3.00000 ** 2.00000)
      CONSTS(158) = CONSTS(155)+ CONSTS(138)*((4.00000 ** 2.00000+2.00000 ** 2.00000) - 2.00000 ** 2.00000)
      CONSTS(159) = CONSTS(156)+ CONSTS(138)*((4.00000 ** 2.00000+1.00000 ** 2.00000) - 3.00000 ** 2.00000)
      CONSTS(160) = 6.48000+ (CONSTS(141)/1.00000)*((4.00000 ** 2.00000+1.00000 ** 2.00000) - 3.00000 ** 2.00000)+ CONSTS(142)*CONSTS(157)
      CONSTS(161) = 4.19000+ (CONSTS(141)/1.00000)*((4.00000 ** 2.00000+2.00000 ** 2.00000) - 2.00000 ** 2.00000)+ CONSTS(142)*CONSTS(158)
      CONSTS(162) = 1.17000+ (CONSTS(141)/1.00000)*((4.00000 ** 2.00000+1.00000 ** 2.00000) - 3.00000 ** 2.00000)+ CONSTS(142)*CONSTS(159)
      CONSTS(163) = 12.0000
      CONSTS(164) = - 2768.10
      CONSTS(165) = 1.00000
      CONSTS(166) = - 3.00000
      CONSTS(167) = - 15.0000
      CONSTS(168) = CONSTS(166)+ CONSTS(138)*((3.00000 ** 2.00000+1.00000 ** 2.00000) - 2.00000 ** 2.00000)
      CONSTS(169) = CONSTS(167)+ CONSTS(138)*((3.00000 ** 2.00000+2.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(170) = 6.38000+ (CONSTS(141)/1.00000)*((3.00000 ** 2.00000+1.00000 ** 2.00000) - 2.00000 ** 2.00000)+ CONSTS(142)*CONSTS(168)
      CONSTS(171) = 3.25000+ (CONSTS(141)/1.00000)*((3.00000 ** 2.00000+2.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(169)
      CONSTS(172) = 12.0000
      CONSTS(173) = - 1906.13
      CONSTS(174) = - 3.00000
      CONSTS(175) = - 7.50000
      CONSTS(176) = CONSTS(174)+ CONSTS(138)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(177) = CONSTS(175)+ CONSTS(138)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(178) = 6.29000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(176)
      CONSTS(179) = 1.92000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)+ CONSTS(142)*CONSTS(177)
      CONSTS(180) = 12.0000
      CONSTS(181) = - 1040.45
      CONSTS(182) = 0.310000
      CONSTS(183) = 2.66000
      CONSTS(184) = 8.19000
      CONSTS(185) = CONSTS(183)+ CONSTS(138)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(186) = CONSTS(184)+ CONSTS(138)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(187) = 4.50000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(185)
      CONSTS(188) = 1.60000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)+ CONSTS(142)*CONSTS(186)
      CONSTS(189) = 8.00000
      CONSTS(190) = 2.30000
      CONSTS(191) = 9.00000
      CONSTS(192) = - 1.70000
      CONSTS(193) = - 12.0000
      CONSTS(194) = CONSTS(192)+ CONSTS(138)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(195) = CONSTS(193)+ CONSTS(138)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(196) = 6.09000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(194)
      CONSTS(197) = 2.48000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+2.00000 ** 2.00000) - 0.00000 ** 2.00000)+ CONSTS(142)*CONSTS(195)
      CONSTS(198) = 11.0000
      CONSTS(199) = - 1756.87
      CONSTS(200) = 6.11000
      CONSTS(201) = 11.0000
      CONSTS(202) = - 1763.94
      CONSTS(203) = 5.89000
      CONSTS(204) = 11.0000
      CONSTS(205) = - 1760.80
      CONSTS(206) = 6.40000
      CONSTS(207) = 5.92000
      CONSTS(208) = 2.70000
      CONSTS(209) = 10.0000
      CONSTS(210) = - 2601.40
      CONSTS(211) = 6.45000
      CONSTS(212) = 5.00000
      CONSTS(213) = - 1288.60
      CONSTS(214) = 1.63000
      CONSTS(215) = - 3.10000
      CONSTS(216) = CONSTS(215)+ CONSTS(138)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)
      CONSTS(217) = 6.22000+ (CONSTS(141)/1.00000)*((2.00000 ** 2.00000+1.00000 ** 2.00000) - 1.00000 ** 2.00000)+ CONSTS(142)*CONSTS(216)
      CONSTS(218) = 7.00000
      CONSTS(219) = - 1339.25
      CONSTS(220) = 5.90000
      CONSTS(221) = 1.57000
      CONSTS(222) = 5.00000
      CONSTS(223) = - 1296.26
      CONSTS(224) = 7.50000
      CONSTS(225) = 4.00000
      CONSTS(226) = - 2356.14
      CONSTS(227) = 6.21000
      CONSTS(228) = 4.00000
      CONSTS(229) = - 1502.54
      CONSTS(230) = 7.00000
      CONSTS(231) = 2.45000
      CONSTS(232) = 1.18000
      CONSTS(233) = 4.00000
      CONSTS(234) = - 1496.38
      CONSTS(235) = 6.35000
      CONSTS(236) = 2.26000
      CONSTS(237) = 1.08000
      CONSTS(238) = 2.00000
      CONSTS(239) = - 1263.65
      CONSTS(240) = 2.49000
      CONSTS(241) = 3.00000
      CONSTS(242) = - 472.270
      CONSTS(243) = 0.980000
      CONSTS(244) = - 0.330000
      CONSTS(245) = CONSTS(244)+ CONSTS(138)*((1.00000 ** 2.00000+1.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(246) = 3.67000+ (CONSTS(141)/1.00000)*((1.00000 ** 2.00000+1.00000 ** 2.00000) - 0.00000 ** 2.00000)+ CONSTS(142)*CONSTS(245)
      CONSTS(247) = 5.00000
      CONSTS(248) = - 516.720
      CONSTS(249) = - 10.0000
      CONSTS(250) = 26.0000
      CONSTS(251) = 0.00000
      CONSTS(252) = 27.0000
      CONSTS(253) = 22.6500
      CONSTS(254) = 2.00000
      CONSTS(255) = - 237.190
      CONSTS(256) = 1.00000
      CONSTS(257) = 0.00000
      CONSTS(258) = 2.58000e+08
      CONSTS(259) = - 17.5500
      CONSTS(260) = CONSTS(259)+ CONSTS(138)*(((2.00000 ** 2.00000+3.00000 ** 2.00000+1.00000 ** 2.00000) - 4.00000 ** 2.00000) - 0.00000 ** 2.00000)
      CONSTS(261) = EXP(log(CONSTS(258))+( CONSTS(140)* CONSTS(126) ** (1.0 / 2)*(((2.00000 ** 2.00000+3.00000 ** 2.00000+1.00000 ** 2.00000) - 4.00000 ** 2.00000) - 0.00000 ** 2.00000))/(1.00000+ CONSTS(137)* CONSTS(126) ** (1.0 / 2)))
      CONSTS(262) = 10.0000 ** arbitrary_log(CONSTS(261), 10) -  CONSTS(142)*CONSTS(260)
      CONSTS(263) =  - CONSTS(1)*CONSTS(131)*log(CONSTS(262))
      CONSTS(264) = CONSTS(5)
      CONSTS(265) = CONSTS(37)
      CONSTS(266) = CONSTS(57)
      CONSTS(267) = CONSTS(58)
      CONSTS(268) = CONSTS(59)
      CONSTS(269) = CONSTS(60)
      CONSTS(270) = CONSTS(61)
      CONSTS(271) = CONSTS(62)
      CONSTS(272) = CONSTS(63)
      CONSTS(273) = CONSTS(64)
      CONSTS(274) = CONSTS(65)
      CONSTS(275) = CONSTS(66)
      CONSTS(276) = CONSTS(67)
      CONSTS(277) = CONSTS(68)
      CONSTS(278) = CONSTS(69)
      CONSTS(279) = CONSTS(70)
      CONSTS(280) = CONSTS(71)
      CONSTS(281) = CONSTS(72)
      CONSTS(282) = CONSTS(73)
      CONSTS(283) = CONSTS(75)
      CONSTS(284) = CONSTS(74)
      CONSTS(285) = CONSTS(76)
      CONSTS(286) = ( CONSTS(280)*CONSTS(282))/( CONSTS(281)*CONSTS(284))
      CONSTS(287) = CONSTS(77)
      CONSTS(288) = CONSTS(78)
      CONSTS(289) = CONSTS(79)
      CONSTS(290) = CONSTS(80)
      CONSTS(291) = CONSTS(81)
      CONSTS(292) = CONSTS(82)
      CONSTS(293) = CONSTS(83)
      CONSTS(294) = CONSTS(291)
      CONSTS(295) = CONSTS(84)
      CONSTS(296) = CONSTS(85)
      CONSTS(297) = CONSTS(86)
      CONSTS(298) = CONSTS(87)
      CONSTS(299) = CONSTS(88)
      CONSTS(300) = CONSTS(295)
      CONSTS(301) = CONSTS(89)
      CONSTS(302) = CONSTS(90)
      CONSTS(303) = CONSTS(91)
      CONSTS(304) = CONSTS(92)
      CONSTS(305) = CONSTS(93)
      CONSTS(306) = CONSTS(94)
      CONSTS(307) = CONSTS(95)
      CONSTS(308) = CONSTS(96)
      CONSTS(309) = CONSTS(97)
      CONSTS(310) = CONSTS(98)
      CONSTS(311) = CONSTS(99)
      CONSTS(312) = CONSTS(307)
      CONSTS(313) = CONSTS(100)
      CONSTS(314) = CONSTS(101)
      CONSTS(315) = CONSTS(102)
      CONSTS(316) = CONSTS(103)
      CONSTS(317) = CONSTS(104)
      CONSTS(318) = CONSTS(105)
      CONSTS(319) = CONSTS(316)
      CONSTS(320) = CONSTS(106)
      CONSTS(321) = CONSTS(107)
      CONSTS(322) = CONSTS(108)
      CONSTS(323) = CONSTS(109)
      CONSTS(324) = CONSTS(110)
      CONSTS(325) = CONSTS(111)
      CONSTS(326) = CONSTS(112)
      CONSTS(327) = CONSTS(113)
      CONSTS(328) = CONSTS(114)
      CONSTS(329) = CONSTS(115)
      CONSTS(330) = CONSTS(116)
      CONSTS(331) = CONSTS(117)
      CONSTS(332) = CONSTS(118)
      CONSTS(333) = CONSTS(119)
      CONSTS(334) = CONSTS(120)
      CONSTS(335) = CONSTS(121)
      CONSTS(336) = CONSTS(122)
      CONSTS(337) = CONSTS(123)
      CONSTS(338) = CONSTS(124)
      CONSTS(339) = CONSTS(128)
      CONSTS(340) = CONSTS(127)
      CONSTS(341) = CONSTS(129)
      CONSTS(342) = CONSTS(133)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(221) = TERNRY(CONSTS(39) .NE. 45.0000, CONSTS(6), TERNRY(VOI.LT.40.0000, 0.00100000, TERNRY(VOI.GE.40.0000.AND.VOI.LT.80.0000, 0.00400000, TERNRY(VOI.GE.80.0000.AND.VOI.LT.100.000, 0.0100000, TERNRY(VOI.GE.100.000, 0.0400000, 0.0/0.0)
      ALGBRC(222) = 1.00000+STATES(5)/CONSTS(40)+STATES(4)/CONSTS(265)+( STATES(5)*STATES(4))/( CONSTS(40)*CONSTS(41))+STATES(5)/CONSTS(267)+STATES(3)/CONSTS(42)+( STATES(5)*STATES(3))/( CONSTS(268)*CONSTS(267))
      ALGBRC(1) = TERNRY(VOI.LE.1.00000.OR.VOI.GT.1.00000.AND.CONSTS(35).EQ.0.00000, STATES(1), CONSTS(36))
      ALGBRC(223) = 1.40400/(1.00000+10.0000 ** 5.94000 - ALGBRC(1)+10.0000 ** ALGBRC(1) - 7.29000)
      ALGBRC(225) = (( ALGBRC(223)*CONSTS(264)*STATES(4))/( CONSTS(266)*CONSTS(265)))/ALGBRC(222)
      ALGBRC(2) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(150)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(151)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(145)
      ALGBRC(61) = 1.00000+ 10.0000 ** - ALGBRC(1)*10.0000 ** CONSTS(196)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(197)
      ALGBRC(12) = (CONSTS(153)+ CONSTS(152)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(152))
      ALGBRC(70) = (CONSTS(199)+ CONSTS(198)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(198))
      ALGBRC(170) = 655.700+ CONSTS(249)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)
      ALGBRC(183) = (ALGBRC(170)+ALGBRC(70)) - ALGBRC(12)
      ALGBRC(184) = ( EXP(- ALGBRC(183)/( CONSTS(1)*CONSTS(131)))*ALGBRC(61))/ALGBRC(2)
      ALGBRC(224) = ( ALGBRC(223)*CONSTS(264)*CONSTS(267)*CONSTS(268))/( CONSTS(266)*CONSTS(265)*ALGBRC(184))
      ALGBRC(226) = (( ALGBRC(224)*STATES(5))/( CONSTS(267)*CONSTS(268)))/ALGBRC(222)
      ALGBRC(227) =  ALGBRC(221)*( STATES(5)*ALGBRC(225) -  STATES(3)*ALGBRC(226))
      ALGBRC(228) = 1.00000 - ALGBRC(221)
      ALGBRC(229) = (STATES(6)/CONSTS(45) ** CONSTS(47)/CONSTS(46))/(1.00000+STATES(6)/CONSTS(45) ** CONSTS(47)/CONSTS(46))
      ALGBRC(230) = 1.00000+STATES(5)/CONSTS(270)+STATES(4)/CONSTS(43)+STATES(5)/CONSTS(272)+STATES(3)/CONSTS(44)+( STATES(5)*STATES(4))/( CONSTS(270)*CONSTS(269))+( STATES(5)*STATES(3))/( CONSTS(271)*CONSTS(272))
      ALGBRC(231) = 1.75000/(1.00000+10.0000 ** 6.12000 - ALGBRC(1)+10.0000 ** ALGBRC(1) - 7.03000)
      ALGBRC(233) = (( ALGBRC(231)*ALGBRC(229)*CONSTS(264)*STATES(4))/( CONSTS(270)*CONSTS(269)))/ALGBRC(230)
      ALGBRC(232) = ( ALGBRC(231)*CONSTS(264)*CONSTS(271)*CONSTS(272))/( CONSTS(270)*CONSTS(269)*ALGBRC(184))
      ALGBRC(234) = (( ALGBRC(229)*ALGBRC(232)*STATES(5))/( CONSTS(271)*CONSTS(272)))/ALGBRC(230)
      ALGBRC(235) =  ALGBRC(228)*( STATES(5)*ALGBRC(233) -  STATES(3)*ALGBRC(234))
      ALGBRC(237) = - (ALGBRC(227)+ALGBRC(235))
      RATES(5) = ALGBRC(237)
      ALGBRC(236) = ( CONSTS(273)*1.32900)/(1.00000+10.0000 ** - ALGBRC(1)+6.64000+10.0000 ** ALGBRC(1) - 8.36000)
      ALGBRC(71) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(200)
      ALGBRC(76) = (CONSTS(202)+ CONSTS(201)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(201))
      ALGBRC(186) = ALGBRC(76) - ALGBRC(70)
      ALGBRC(187) = ( EXP(- ALGBRC(186)/( CONSTS(1)*CONSTS(131)))*ALGBRC(71))/ALGBRC(61)
      ALGBRC(238) = ( ALGBRC(236)*CONSTS(275))/( CONSTS(274)*ALGBRC(187))
      ALGBRC(239) = (( ALGBRC(236)*STATES(3))/CONSTS(274) - ( ALGBRC(238)*STATES(7))/CONSTS(275))/(1.00000+STATES(3)/CONSTS(274)+STATES(7)/CONSTS(275))
      ALGBRC(240) = CONSTS(276)/(1.00000+10.0000 ** - ALGBRC(1)+6.94000+10.0000 ** ALGBRC(1) - 9.35000)
      ALGBRC(77) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(203)
      ALGBRC(82) = (CONSTS(205)+ CONSTS(204)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(204))
      ALGBRC(189) = ALGBRC(82) - ALGBRC(76)
      ALGBRC(190) = ( EXP(- ALGBRC(189)/( CONSTS(1)*CONSTS(131)))*ALGBRC(77))/ALGBRC(71)
      ALGBRC(241) =  (( ALGBRC(240)*CONSTS(277))/CONSTS(278))*ALGBRC(190)
      ALGBRC(242) = (( ALGBRC(241)*STATES(7))/CONSTS(277) - ( ALGBRC(240)*STATES(8))/CONSTS(278))/(1.00000+STATES(8)/CONSTS(278)+STATES(7)/CONSTS(277))
      ALGBRC(244) = ALGBRC(239) - ALGBRC(242)
      RATES(7) = ALGBRC(244)
      ALGBRC(243) = CONSTS(279)/(1.00000+ALGBRC(1)/6.80000 ** - 30.0000)
      ALGBRC(13) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(160)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(161)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(162)
      ALGBRC(24) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(170)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(171)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(165)
      ALGBRC(83) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(206)+10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(208)
      ALGBRC(23) = (CONSTS(164)+ CONSTS(163)*CONSTS(1)*CONSTS(131)*log(10.0000)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(163))
      ALGBRC(35) = (CONSTS(173)+ CONSTS(172)*CONSTS(1)*CONSTS(131)*log(10.0000)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(172))
      ALGBRC(93) = (CONSTS(210)+ CONSTS(209)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(209))
      ALGBRC(174) = CONSTS(257)+ CONSTS(256)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)
      ALGBRC(192) = ((ALGBRC(93)+ALGBRC(35)+ALGBRC(174)) - ALGBRC(82)) - ALGBRC(23)
      ALGBRC(193) = ( EXP(- ALGBRC(192)/( CONSTS(1)*CONSTS(131)))*ALGBRC(83)*ALGBRC(24))/( ALGBRC(77)*ALGBRC(13)*10.0000 ** - ALGBRC(1))
      ALGBRC(245) = ( ALGBRC(243)*CONSTS(283)*CONSTS(285))/( CONSTS(280)*CONSTS(282)*ALGBRC(193))
      ALGBRC(246) =  CONSTS(54)*( ((1.00000+STATES(11)/CONSTS(50))/(1.00000+( CONSTS(52)*STATES(11))/CONSTS(50)))*(1.00000+( CONSTS(53)*STATES(6))/CONSTS(51)))/(1.00000+STATES(6)/CONSTS(51)) ** 4.00000
      ALGBRC(247) =  (1.00000+STATES(8)/CONSTS(280))*(1.00000+STATES(11)/CONSTS(282))+STATES(9)/CONSTS(283)+ (STATES(10)/CONSTS(285))*(1.00000+STATES(9)/CONSTS(283))
      ALGBRC(248) =  (1.00000+STATES(8)/CONSTS(281))*(1.00000+STATES(11)/CONSTS(284))+STATES(9)/CONSTS(48)+ (STATES(10)/CONSTS(49))*(1.00000+STATES(9)/CONSTS(48))
      ALGBRC(249) = ( ((( ALGBRC(243)*STATES(8)*STATES(11))/( CONSTS(280)*CONSTS(282)))/ALGBRC(247) - (( ALGBRC(245)*STATES(10)*STATES(9))/( CONSTS(285)*CONSTS(283)))/ALGBRC(247))*(1.00000+ CONSTS(286)*ALGBRC(246)*ALGBRC(248)/ALGBRC(247) ** 3.00000))/(1.00000+ ALGBRC(246)*ALGBRC(248)/ALGBRC(247) ** 4.00000)
      ALGBRC(251) = ALGBRC(242) - ALGBRC(249)
      RATES(8) = ALGBRC(251)
      ALGBRC(250) = ( CONSTS(287)*1.01300)/(1.00000+10.0000 ** - ALGBRC(1)+5.32000+10.0000 ** ALGBRC(1) - 9.15000)
      ALGBRC(94) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(211)
      ALGBRC(110) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(220)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(221)
      ALGBRC(99) = (CONSTS(213)+ CONSTS(212)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(212))
      ALGBRC(119) = (CONSTS(223)+ CONSTS(222)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(222))
      ALGBRC(195) = (ALGBRC(119)+ALGBRC(99)) - ALGBRC(93)
      ALGBRC(196) = ( 1.00000*EXP(- ALGBRC(195)/( CONSTS(1)*CONSTS(131)))*ALGBRC(94)*ALGBRC(83))/ALGBRC(110)
      ALGBRC(252) = ( ALGBRC(250)*CONSTS(290)*CONSTS(289))/( CONSTS(288)*ALGBRC(196))
      ALGBRC(253) = (( ALGBRC(250)*STATES(9))/CONSTS(288) - ( ALGBRC(252)*STATES(13)*STATES(12))/( CONSTS(290)*CONSTS(289)))/(1.00000+STATES(9)/CONSTS(288)+STATES(13)/CONSTS(290)+STATES(12)/CONSTS(289))
      ALGBRC(198) = ALGBRC(119) - ALGBRC(99)
      ALGBRC(199) = ( EXP(- ALGBRC(198)/( CONSTS(1)*CONSTS(131)))*ALGBRC(110))/ALGBRC(94)
      ALGBRC(254) = ( CONSTS(294)*CONSTS(293))/( CONSTS(292)*ALGBRC(199))
      ALGBRC(255) = (( CONSTS(294)*STATES(13))/CONSTS(292) - ( ALGBRC(254)*STATES(12))/CONSTS(293))/(1.00000+STATES(13)/CONSTS(292)+STATES(12)/CONSTS(293))
      ALGBRC(259) = 1.00000+STATES(4)/CONSTS(304)+STATES(13)/CONSTS(302)+STATES(15)/CONSTS(303)+( STATES(13)*STATES(15))/( CONSTS(302)*CONSTS(303))+( STATES(13)*STATES(15)*STATES(4))/( CONSTS(302)*CONSTS(303)*CONSTS(304))+STATES(17)/CONSTS(305)+STATES(16)/CONSTS(306)+( STATES(17)*STATES(16))/( CONSTS(306)*CONSTS(305))
      ALGBRC(260) =  CONSTS(301)*0.000700000*EXP( ALGBRC(1)*0.897900)
      ALGBRC(120) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(224)
      ALGBRC(125) = (CONSTS(226)+ CONSTS(225)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(225))
      ALGBRC(171) = (CONSTS(251)+ CONSTS(250)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(250))
      ALGBRC(172) = (CONSTS(253)+ CONSTS(252)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(252))
      ALGBRC(201) = (((ALGBRC(125)+ALGBRC(172)+ALGBRC(174)) - ALGBRC(12)) - ALGBRC(99)) - ALGBRC(171)
      ALGBRC(202) = ( EXP(- ALGBRC(201)/( CONSTS(1)*CONSTS(131)))*ALGBRC(120))/( ALGBRC(2)*ALGBRC(94)*10.0000 ** - ALGBRC(1)*1.00000)
      ALGBRC(261) = ( ALGBRC(260)*CONSTS(305)*CONSTS(306))/( CONSTS(302)*CONSTS(304)*CONSTS(303)*ALGBRC(202))
      ALGBRC(262) = (( ALGBRC(260)*STATES(13)*STATES(15)*STATES(4))/( CONSTS(303)*CONSTS(302)*CONSTS(304)) - ( ALGBRC(261)*STATES(17)*STATES(16))/( CONSTS(305)*CONSTS(306)))/ALGBRC(259)
      ALGBRC(264) = (ALGBRC(253) - ALGBRC(255)) - ALGBRC(262)
      RATES(13) = ALGBRC(264)
      ALGBRC(126) = 1.00000+10.0000 ** - ALGBRC(1)+6.21000
      ALGBRC(131) = (CONSTS(229)+ CONSTS(228)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(228))
      ALGBRC(207) = ((ALGBRC(23)+ALGBRC(131)) - ALGBRC(125)) - ALGBRC(35)
      ALGBRC(208) = ( EXP(- ALGBRC(207)/( CONSTS(1)*CONSTS(131)))*ALGBRC(13)*ALGBRC(126))/( ALGBRC(120)*ALGBRC(24))
      ALGBRC(263) =  (( CONSTS(312)*CONSTS(308)*CONSTS(309))/( CONSTS(310)*CONSTS(311)))*ALGBRC(208)
      ALGBRC(265) = 1.00000+STATES(10)/CONSTS(309)+STATES(17)/CONSTS(308)+( STATES(17)*STATES(10))/( CONSTS(308)*CONSTS(309))+STATES(18)/CONSTS(310)+STATES(11)/CONSTS(311)+( STATES(18)*STATES(11))/( CONSTS(310)*CONSTS(311))
      ALGBRC(266) = (( ALGBRC(263)*STATES(17)*STATES(10))/( CONSTS(309)*CONSTS(308)) - ( CONSTS(312)*STATES(11)*STATES(18))/( CONSTS(311)*CONSTS(310)))/ALGBRC(265)
      ALGBRC(268) = ALGBRC(262) - ALGBRC(266)
      RATES(17) = ALGBRC(268)
      ALGBRC(267) = ( CONSTS(313)*0.989000)/(1.00000+10.0000 ** - ALGBRC(1)+5.62000+10.0000 ** ALGBRC(1) - 8.74000)
      ALGBRC(132) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(230)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(231)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(232)
      ALGBRC(142) = (CONSTS(234)+ CONSTS(233)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(233))
      ALGBRC(210) = ALGBRC(142) - ALGBRC(131)
      ALGBRC(211) = ( EXP(- ALGBRC(210)/( CONSTS(1)*CONSTS(131)))*ALGBRC(132))/ALGBRC(126)
      ALGBRC(269) = ( ALGBRC(267)*CONSTS(315))/( CONSTS(314)*ALGBRC(211))
      ALGBRC(270) = (( ALGBRC(267)*STATES(18))/CONSTS(314) - ( ALGBRC(269)*STATES(19))/CONSTS(315))/(1.00000+STATES(18)/CONSTS(314)+STATES(19)/CONSTS(315))
      ALGBRC(272) = ALGBRC(266) - ALGBRC(270)
      RATES(18) = ALGBRC(272)
      ALGBRC(256) =  (1.00000+STATES(14)/CONSTS(296)+STATES(16)/CONSTS(299))*(1.00000+STATES(12)/CONSTS(298)+STATES(15)/CONSTS(299))
      ALGBRC(100) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(217)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(214)
      ALGBRC(109) = (CONSTS(219)+ CONSTS(218)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(218))
      ALGBRC(204) = ((ALGBRC(174)+ALGBRC(172)+ALGBRC(119)) - ALGBRC(171)) - ALGBRC(109)
      ALGBRC(205) = ( EXP(- ALGBRC(204)/( CONSTS(1)*CONSTS(131)))*ALGBRC(110))/( ALGBRC(100)*10.0000 ** - ALGBRC(1))
      ALGBRC(257) = ( CONSTS(300)*CONSTS(296)*CONSTS(297)*ALGBRC(205))/( CONSTS(298)*CONSTS(299))
      ALGBRC(258) = (( ALGBRC(257)*STATES(14)*STATES(15))/( CONSTS(296)*CONSTS(297)) - ( CONSTS(300)*STATES(12)*STATES(16))/( CONSTS(298)*CONSTS(299)))/ALGBRC(256)
      ALGBRC(277) =  CONSTS(325)*( - 0.113400*ALGBRC(1)+1.60690)
      ALGBRC(154) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(240)
      ALGBRC(160) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(246)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(243)
      ALGBRC(159) = (CONSTS(242)+ CONSTS(241)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(241))
      ALGBRC(169) = (CONSTS(248)+ CONSTS(247)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(247))
      ALGBRC(219) = (((ALGBRC(169)+ALGBRC(171)) - ALGBRC(159)) - ALGBRC(172)) - ALGBRC(174)
      ALGBRC(220) = ( EXP(- ALGBRC(219)/( CONSTS(1)*CONSTS(131)))*ALGBRC(160)*10.0000 ** - ALGBRC(1))/ALGBRC(154)
      ALGBRC(278) = ( ALGBRC(277)*CONSTS(328)*CONSTS(329))/( CONSTS(326)*CONSTS(327)*ALGBRC(220))
      ALGBRC(279) = (( ALGBRC(277)*STATES(21)*STATES(16))/( CONSTS(326)*CONSTS(327)) - ( ALGBRC(278)*STATES(22)*STATES(15))/( CONSTS(328)*CONSTS(329)))/(1.00000+STATES(21)/CONSTS(326)+STATES(16)/CONSTS(327)+( STATES(21)*STATES(16))/( CONSTS(326)*CONSTS(327))+STATES(22)/CONSTS(328)+STATES(15)/CONSTS(329)+( STATES(22)*STATES(15))/( CONSTS(328)*CONSTS(329)))
      ALGBRC(282) = (- ALGBRC(262) - ALGBRC(258))+ALGBRC(279)
      RATES(15) = ALGBRC(282)
      ALGBRC(283) = (ALGBRC(262)+ALGBRC(258)) - ALGBRC(279)
      RATES(16) = ALGBRC(283)
      ALGBRC(274) = ( CONSTS(320)*1.05000)/(1.00000+10.0000 ** - ALGBRC(1)+5.58000+10.0000 ** ALGBRC(1) - 8.79000)
      ALGBRC(143) = 1.00000+10.0000 ** CONSTS(235) - ALGBRC(1)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(236)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(237)
      ALGBRC(153) = (CONSTS(239)+ CONSTS(238)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(238))
      ALGBRC(216) = (((ALGBRC(159)+ALGBRC(23)) - ALGBRC(174)) - ALGBRC(153)) - ALGBRC(35)
      ALGBRC(217) = ( EXP(- ALGBRC(216)/( CONSTS(1)*CONSTS(131)))*ALGBRC(154)*ALGBRC(13)*10.0000 ** - ALGBRC(1))/( ALGBRC(143)*ALGBRC(24))
      ALGBRC(275) = ( ALGBRC(274)*CONSTS(323)*CONSTS(324))/( CONSTS(321)*CONSTS(322)*ALGBRC(217))
      ALGBRC(276) = (( ALGBRC(274)*STATES(20)*STATES(10))/( CONSTS(321)*CONSTS(322)) - ( ALGBRC(275)*STATES(21)*STATES(11))/( CONSTS(323)*CONSTS(324)))/(1.00000+STATES(20)/CONSTS(321)+STATES(10)/CONSTS(322)+( STATES(20)*STATES(10))/( CONSTS(321)*CONSTS(322))+STATES(11)/CONSTS(324)+STATES(21)/CONSTS(323)+( STATES(21)*STATES(11))/( CONSTS(323)*CONSTS(324)))
      ALGBRC(284) = ALGBRC(276) - ALGBRC(279)
      RATES(21) = ALGBRC(284)
      ALGBRC(46) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(187)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(188)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(182)
      ALGBRC(56) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(190)
      ALGBRC(177) = ( EXP(- CONSTS(263)/( CONSTS(1)*CONSTS(131)))*10.0000 ** - ALGBRC(1)*ALGBRC(13)*ALGBRC(56))/( ALGBRC(46)*ALGBRC(24))
      ALGBRC(286) = ( (CONSTS(330)/ALGBRC(177))*CONSTS(332)*CONSTS(334))/( CONSTS(333)*CONSTS(331))
      ALGBRC(287) = (( ALGBRC(286)*STATES(11)*STATES(23))/( CONSTS(332)*CONSTS(334)) - ( CONSTS(330)*STATES(10)*STATES(24))/( CONSTS(333)*CONSTS(331)))/(1.00000+STATES(10)/CONSTS(333)+STATES(24)/CONSTS(56)+( STATES(24)*STATES(10))/( CONSTS(333)*CONSTS(331))+STATES(11)/CONSTS(332)+( STATES(23)*STATES(11))/( CONSTS(334)*CONSTS(332)))
      ALGBRC(290) = - ALGBRC(287)
      RATES(23) = ALGBRC(290)
      ALGBRC(14) = 1.00000/ALGBRC(13)
      ALGBRC(15) =  10.0000 ** - ALGBRC(1)+CONSTS(160)*ALGBRC(14)
      ALGBRC(16) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(161)*ALGBRC(14)
      ALGBRC(17) =  (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(162)*ALGBRC(14)
      ALGBRC(18) =  0.00000*ALGBRC(14)+ 1.00000*ALGBRC(15)+ 0.00000*ALGBRC(16)+ 0.00000*ALGBRC(17)
      ALGBRC(25) = 1.00000/ALGBRC(24)
      ALGBRC(26) =  10.0000 ** - ALGBRC(1)+CONSTS(170)*ALGBRC(25)
      ALGBRC(27) =  (( ALGBRC(25)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(171)
      ALGBRC(28) =  (( ALGBRC(25)*CONSTS(134))/CONSTS(135))*10.0000 ** CONSTS(165)
      ALGBRC(29) =  0.00000*ALGBRC(25)+ 1.00000*ALGBRC(26)+ 0.00000*ALGBRC(27)+ 0.00000*ALGBRC(28)
      ALGBRC(47) = 1.00000/ALGBRC(46)
      ALGBRC(49) =  10.0000 ** - ALGBRC(1)+CONSTS(187)*ALGBRC(47)
      ALGBRC(51) = ALGBRC(49)
      ALGBRC(57) = 1.00000/ALGBRC(56)
      ALGBRC(58) =  ALGBRC(57)*10.0000 ** - ALGBRC(1)+CONSTS(190)
      ALGBRC(59) =  0.00000*ALGBRC(57)+ 1.00000*ALGBRC(58)
      ALGBRC(175) = (((ALGBRC(18)+ALGBRC(59)) - ALGBRC(51)) - ALGBRC(29))+(((CONSTS(163)+CONSTS(191)) - CONSTS(189)) - CONSTS(172))
      ALGBRC(36) = 1.00000+10.0000 ** CONSTS(178) - ALGBRC(1)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(179)
      ALGBRC(37) = 1.00000/ALGBRC(36)
      ALGBRC(38) =  ALGBRC(37)*10.0000 ** CONSTS(178) - ALGBRC(1)
      ALGBRC(39) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(179)*ALGBRC(37)
      ALGBRC(40) =  0.00000*ALGBRC(37)+ALGBRC(38)+ 0.00000*ALGBRC(39)
      ALGBRC(179) = (( 2.00000*ALGBRC(29) - ALGBRC(18)) - ALGBRC(40))+(( 2.00000*CONSTS(172) - CONSTS(163)) - CONSTS(180))
      ALGBRC(3) = 1.00000/ALGBRC(2)
      ALGBRC(4) =  10.0000 ** - ALGBRC(1)+CONSTS(150)*ALGBRC(3)
      ALGBRC(7) =  1.00000*ALGBRC(4)
      ALGBRC(62) = 1.00000/ALGBRC(61)
      ALGBRC(63) =  ALGBRC(62)*10.0000 ** - ALGBRC(1)+CONSTS(196)
      ALGBRC(65) = ALGBRC(63)
      ALGBRC(182) = (ALGBRC(65) - ALGBRC(7))+(1.00000 - CONSTS(152))
      ALGBRC(72) = 1.00000/ALGBRC(71)
      ALGBRC(73) =  ALGBRC(72)*10.0000 ** - ALGBRC(1)+CONSTS(200)
      ALGBRC(74) = ALGBRC(73)
      ALGBRC(185) = (ALGBRC(74) - ALGBRC(65))+(CONSTS(201) - CONSTS(198))
      ALGBRC(78) = 1.00000/ALGBRC(77)
      ALGBRC(79) =  ALGBRC(78)*10.0000 ** - ALGBRC(1)+CONSTS(203)
      ALGBRC(80) = ALGBRC(79)
      ALGBRC(188) = (ALGBRC(80) - ALGBRC(74))+(CONSTS(204) - CONSTS(201))
      ALGBRC(84) = 1.00000/ALGBRC(83)
      ALGBRC(85) =  ALGBRC(84)*10.0000 ** - ALGBRC(1)+CONSTS(206)
      ALGBRC(86) =  ALGBRC(84)*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)
      ALGBRC(88) = ALGBRC(85)+ 2.00000*ALGBRC(86)
      ALGBRC(191) = (((((ALGBRC(29)+ALGBRC(88)) - ALGBRC(80)) - ALGBRC(18))+CONSTS(172)+CONSTS(209)) - CONSTS(204)) - CONSTS(163)
      ALGBRC(95) = 1.00000/ALGBRC(94)
      ALGBRC(96) =  ALGBRC(95)*10.0000 ** - ALGBRC(1)+CONSTS(211)
      ALGBRC(97) = ALGBRC(96)
      ALGBRC(111) = 1.00000/ALGBRC(110)
      ALGBRC(112) =  ALGBRC(111)*10.0000 ** - ALGBRC(1)+CONSTS(220)
      ALGBRC(114) = ALGBRC(112)
      ALGBRC(194) = ((ALGBRC(114)+ALGBRC(97)) - ALGBRC(88))+((CONSTS(222)+CONSTS(212)) - CONSTS(209))
      ALGBRC(197) = (ALGBRC(114) - ALGBRC(97))+(CONSTS(222) - CONSTS(212))
      ALGBRC(121) = 1.00000/ALGBRC(120)
      ALGBRC(122) =  ALGBRC(121)*10.0000 ** - ALGBRC(1)+CONSTS(224)
      ALGBRC(123) = ALGBRC(122)
      ALGBRC(200) = ((ALGBRC(123) - ALGBRC(97)) - ALGBRC(7))+((((CONSTS(225)+CONSTS(252)) - CONSTS(212)) - CONSTS(152)) - CONSTS(250))
      ALGBRC(101) = 1.00000/ALGBRC(100)
      ALGBRC(102) =  ALGBRC(101)*10.0000 ** - ALGBRC(1)+CONSTS(217)
      ALGBRC(104) = ALGBRC(102)
      ALGBRC(203) = (ALGBRC(114) - ALGBRC(104))+(((CONSTS(222)+CONSTS(252)) - CONSTS(250)) - CONSTS(218))
      ALGBRC(127) = 1.00000/ALGBRC(126)
      ALGBRC(128) =  ALGBRC(127)*10.0000 ** - ALGBRC(1)+6.21000
      ALGBRC(129) = ALGBRC(128)
      ALGBRC(206) = (((ALGBRC(129)+ALGBRC(18)) - ALGBRC(123)) - ALGBRC(29))+(((CONSTS(228)+CONSTS(163)) - CONSTS(225)) - CONSTS(172))
      ALGBRC(133) = 1.00000/ALGBRC(132)
      ALGBRC(134) =  ALGBRC(133)*10.0000 ** - ALGBRC(1)+CONSTS(230)
      ALGBRC(137) = ALGBRC(134)
      ALGBRC(209) = (ALGBRC(137) - ALGBRC(129))+(CONSTS(233) - CONSTS(228))
      ALGBRC(144) = 1.00000/ALGBRC(143)
      ALGBRC(145) =  ALGBRC(144)*10.0000 ** CONSTS(235) - ALGBRC(1)
      ALGBRC(148) = ALGBRC(145)
      ALGBRC(212) = (ALGBRC(148) - ALGBRC(137))+((CONSTS(254)+CONSTS(238)) - CONSTS(233))
      ALGBRC(155) = 1.00000/ALGBRC(154)
      ALGBRC(156) =  ALGBRC(155)*10.0000 ** - ALGBRC(1)+CONSTS(240)
      ALGBRC(157) = ALGBRC(156)
      ALGBRC(215) = (((ALGBRC(157)+ALGBRC(18)) - ALGBRC(148)) - ALGBRC(29))+(((CONSTS(241)+CONSTS(163)) - CONSTS(238)) - CONSTS(172))
      ALGBRC(161) = 1.00000/ALGBRC(160)
      ALGBRC(162) =  ALGBRC(161)*10.0000 ** - ALGBRC(1)+CONSTS(246)
      ALGBRC(164) = ALGBRC(162)
      ALGBRC(218) = (ALGBRC(164) - ALGBRC(157))+(((CONSTS(247)+CONSTS(250)) - CONSTS(241)) - CONSTS(252))
      ALGBRC(173) = (CONSTS(255)+ CONSTS(254)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(0.00000 - CONSTS(254))
      ALGBRC(213) = (ALGBRC(173)+ALGBRC(153)) - ALGBRC(142)
      ALGBRC(214) = ( EXP(- ALGBRC(213)/( CONSTS(1)*CONSTS(131)))*ALGBRC(143))/ALGBRC(132)
      ALGBRC(271) = ( CONSTS(319)*CONSTS(318))/( CONSTS(317)*ALGBRC(214))
      ALGBRC(273) = (( CONSTS(319)*STATES(19))/CONSTS(317) - ( ALGBRC(271)*STATES(20))/CONSTS(318))/(1.00000+STATES(20)/CONSTS(318)+STATES(19)/CONSTS(317))
      ALGBRC(45) = (CONSTS(181)+ CONSTS(180)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(180))
      ALGBRC(180) = ( 2.00000*ALGBRC(35) - ALGBRC(23)) - ALGBRC(45)
      ALGBRC(181) = ( EXP(- ALGBRC(180)/( CONSTS(1)*CONSTS(131)))*ALGBRC(24) ** 2.00000)/( ALGBRC(13)*ALGBRC(36))
      ALGBRC(288) = ( CONSTS(335)*CONSTS(338) ** 2.00000)/( CONSTS(336)*CONSTS(337)*ALGBRC(181))
      ALGBRC(291) = (( CONSTS(335)*STATES(11)*STATES(6))/( CONSTS(337)*CONSTS(336)) -  ALGBRC(288)*STATES(10)/CONSTS(338) ** 2.00000)/(1.00000+STATES(11)/CONSTS(337)+STATES(6)/CONSTS(336)+( STATES(11)*STATES(6))/( CONSTS(337)*CONSTS(336))+( 2.00000*STATES(10))/CONSTS(338)+STATES(10) ** 2.00000/CONSTS(338) ** 2.00000)
      ALGBRC(294) =  ALGBRC(175)*- ALGBRC(287)+ ALGBRC(179)*ALGBRC(291)+ ALGBRC(185)*ALGBRC(239)+ ALGBRC(182)*(ALGBRC(227)+ALGBRC(235))+ ALGBRC(188)*ALGBRC(242)+ ALGBRC(191)*ALGBRC(249)+ ALGBRC(194)*ALGBRC(253)+ ALGBRC(197)*ALGBRC(255)+ ALGBRC(200)*ALGBRC(262)+ ALGBRC(206)*ALGBRC(266)+ ALGBRC(209)*ALGBRC(270)+ ALGBRC(212)*ALGBRC(273)+ ALGBRC(215)*ALGBRC(276)+ ALGBRC(218)*ALGBRC(279)+ ALGBRC(203)*ALGBRC(258)
      RATES(25) = - ALGBRC(294)
      ALGBRC(300) = ALGBRC(287)
      RATES(24) = ALGBRC(300)
      ALGBRC(280) = TERNRY(CONSTS(39) .NE. 5.00000.AND.CONSTS(39) .NE. 45.0000, CONSTS(23), TERNRY(CONSTS(39) .NE. 5.00000.AND.CONSTS(39).EQ.45.0000.AND.VOI.GE.110.000, CONSTS(23), TERNRY(CONSTS(39).EQ.5.00000.AND.VOI.GT.30.0000, CONSTS(23), 0.00000)
      ALGBRC(285) = ( ALGBRC(280)*STATES(11))/(CONSTS(55)+STATES(11))
      ALGBRC(301) = (((- ALGBRC(287) - ALGBRC(291)) - ALGBRC(249))+ALGBRC(266)+ALGBRC(276)) - ALGBRC(285)
      RATES(11) = ALGBRC(301)
      ALGBRC(302) = (((ALGBRC(287)+ 2.00000*ALGBRC(291)+ALGBRC(249)) - ALGBRC(266)) - ALGBRC(276))+ALGBRC(285)
      RATES(10) = ALGBRC(302)
      ALGBRC(303) = - ALGBRC(291)
      RATES(6) = ALGBRC(303)
      ALGBRC(304) = (- (ALGBRC(227)+ALGBRC(235)) - ALGBRC(262))+ALGBRC(285)
      RATES(4) = ALGBRC(304)
      ALGBRC(305) = (ALGBRC(227)+ALGBRC(235)) - ALGBRC(239)
      RATES(3) = ALGBRC(305)
      ALGBRC(306) = ALGBRC(249) - ALGBRC(253)
      RATES(9) = ALGBRC(306)
      ALGBRC(307) = ALGBRC(253)+ALGBRC(255)+ALGBRC(258)
      RATES(12) = ALGBRC(307)
      ALGBRC(308) = - ALGBRC(258)
      RATES(14) = ALGBRC(308)
      ALGBRC(309) = ALGBRC(270) - ALGBRC(273)
      RATES(19) = ALGBRC(309)
      ALGBRC(310) = ALGBRC(273) - ALGBRC(276)
      RATES(20) = ALGBRC(310)
      ALGBRC(311) = ALGBRC(279)
      RATES(22) = ALGBRC(311)
      ALGBRC(292) = ( log(10.0000)*CONSTS(339)*10.0000 ** - ALGBRC(1) - 6.87000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 6.87000 ** 2.00000+( log(10.0000)*CONSTS(340)*10.0000 ** - ALGBRC(1) - 8.30000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 8.30000 ** 2.00000+( log(10.0000)*CONSTS(341)*10.0000 ** - ALGBRC(1) - 4.80000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 4.80000 ** 2.00000
      ALGBRC(8) = ( 10.0000 ** CONSTS(150)*(ALGBRC(2) - 10.0000 ** - ALGBRC(1)+CONSTS(150)))/( CONSTS(135)*ALGBRC(2) ** 2.00000)
      ALGBRC(19) = ( 10.0000 ** CONSTS(160)*(ALGBRC(13) - 10.0000 ** - ALGBRC(1)+CONSTS(160)))/( CONSTS(135)*ALGBRC(13) ** 2.00000)
      ALGBRC(30) = ( 10.0000 ** CONSTS(170)*(ALGBRC(24) - 10.0000 ** - ALGBRC(1)+CONSTS(170)))/( CONSTS(135)*ALGBRC(24) ** 2.00000)
      ALGBRC(41) = ( 10.0000 ** CONSTS(178)*(ALGBRC(36) - 10.0000 ** - ALGBRC(1)+CONSTS(178)))/( CONSTS(135)*ALGBRC(36) ** 2.00000)
      ALGBRC(52) = ( 10.0000 ** CONSTS(187)*(ALGBRC(46) - 10.0000 ** - ALGBRC(1)+CONSTS(187)))/( CONSTS(135)*ALGBRC(46) ** 2.00000)
      ALGBRC(60) = ( 10.0000 ** CONSTS(190)*(ALGBRC(56) - 10.0000 ** - ALGBRC(1)+CONSTS(190)))/( CONSTS(135)*ALGBRC(56) ** 2.00000)
      ALGBRC(66) = ( 10.0000 ** CONSTS(196)*(ALGBRC(61) - 10.0000 ** - ALGBRC(1)+CONSTS(196)))/( CONSTS(135)*ALGBRC(61) ** 2.00000)
      ALGBRC(75) = ( 10.0000 ** CONSTS(200)*(ALGBRC(71) - 10.0000 ** - ALGBRC(1)+CONSTS(200)))/( CONSTS(135)*ALGBRC(71) ** 2.00000)
      ALGBRC(81) = ( 10.0000 ** CONSTS(203)*(ALGBRC(77) - 10.0000 ** - ALGBRC(1)+CONSTS(203)))/( CONSTS(135)*ALGBRC(77) ** 2.00000)
      ALGBRC(89) = ( ALGBRC(83)*(10.0000 ** CONSTS(206)+ 2.00000*10.0000 ** - ALGBRC(1)+CONSTS(206)+CONSTS(207)) -  (10.0000 ** - ALGBRC(1)+CONSTS(206)+ 2.00000*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207))*(10.0000 ** CONSTS(206)+ 2.00000*10.0000 ** - ALGBRC(1)+CONSTS(206)+CONSTS(207)))/( CONSTS(135)*ALGBRC(83) ** 2.00000)
      ALGBRC(98) = ( 10.0000 ** CONSTS(211)*(ALGBRC(94) - 10.0000 ** - ALGBRC(1)+CONSTS(211)))/( CONSTS(135)*ALGBRC(94) ** 2.00000)
      ALGBRC(105) = ( 10.0000 ** CONSTS(217)*(ALGBRC(100) - 10.0000 ** - ALGBRC(1)+CONSTS(217)))/( CONSTS(135)*ALGBRC(100) ** 2.00000)
      ALGBRC(115) = ( 10.0000 ** CONSTS(220)*(ALGBRC(110) - 10.0000 ** - ALGBRC(1)+CONSTS(220)))/( CONSTS(135)*ALGBRC(110) ** 2.00000)
      ALGBRC(124) = ( 10.0000 ** CONSTS(224)*(ALGBRC(120) - 10.0000 ** - ALGBRC(1)+CONSTS(224)))/( CONSTS(135)*ALGBRC(120) ** 2.00000)
      ALGBRC(130) = ( 10.0000 ** CONSTS(227)*(ALGBRC(126) - 10.0000 ** - ALGBRC(1)+CONSTS(227)))/( CONSTS(135)*ALGBRC(126) ** 2.00000)
      ALGBRC(138) = ( 10.0000 ** CONSTS(230)*(ALGBRC(132) - 10.0000 ** - ALGBRC(1)+CONSTS(230)))/( CONSTS(135)*ALGBRC(132) ** 2.00000)
      ALGBRC(149) = ( 10.0000 ** CONSTS(235)*(ALGBRC(143) - 10.0000 ** - ALGBRC(1)+CONSTS(235)))/( CONSTS(135)*ALGBRC(143) ** 2.00000)
      ALGBRC(158) = ( 10.0000 ** CONSTS(240)*(ALGBRC(154) - 10.0000 ** - ALGBRC(1)+CONSTS(240)))/( CONSTS(135)*ALGBRC(154) ** 2.00000)
      ALGBRC(165) = ( 10.0000 ** CONSTS(246)*(ALGBRC(160) - 10.0000 ** - ALGBRC(1)+CONSTS(246)))/( CONSTS(135)*ALGBRC(160) ** 2.00000)
      ALGBRC(293) =  log(10.0000)*10.0000 ** - ALGBRC(1)*CONSTS(135)*(1.00000+ ALGBRC(8)*STATES(4)+ ALGBRC(19)*STATES(11)+ ALGBRC(30)*STATES(10)+ ALGBRC(41)*STATES(6)+ ALGBRC(52)*STATES(24)+ ALGBRC(60)*STATES(23)+ ALGBRC(66)*STATES(3)+ ALGBRC(75)*STATES(7)+ ALGBRC(81)*STATES(8)+ ALGBRC(89)*STATES(9)+ ALGBRC(98)*STATES(13)+ ALGBRC(115)*STATES(12)+ ALGBRC(105)*STATES(14)+ ALGBRC(124)*STATES(17)+ ALGBRC(130)*STATES(18)+ ALGBRC(138)*STATES(19)+ ALGBRC(149)*STATES(20)+ ALGBRC(158)*STATES(21)+ ALGBRC(165)*STATES(22))
      ALGBRC(295) = ALGBRC(294)/(ALGBRC(292)+ALGBRC(293))
      ALGBRC(9) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(150)+CONSTS(151))/ALGBRC(2) ** 2.00000
      ALGBRC(20) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(160)+CONSTS(161))/ALGBRC(13) ** 2.00000
      ALGBRC(32) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(170)+CONSTS(171))/ALGBRC(24) ** 2.00000
      ALGBRC(42) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(178)+CONSTS(179))/ALGBRC(36) ** 2.00000
      ALGBRC(53) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(187)+CONSTS(188))/ALGBRC(46) ** 2.00000
      ALGBRC(67) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(196)+CONSTS(197))/ALGBRC(61) ** 2.00000
      ALGBRC(90) = ( - (10.0000 ** - ALGBRC(1)+CONSTS(206)+ 2.00000*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207))*10.0000 ** CONSTS(208))/( CONSTS(135)*ALGBRC(83) ** 2.00000)
      ALGBRC(106) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(217)+CONSTS(214))/ALGBRC(100) ** 2.00000
      ALGBRC(116) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(220)+CONSTS(221))/ALGBRC(110) ** 2.00000
      ALGBRC(139) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(230)+CONSTS(231))/ALGBRC(132) ** 2.00000
      ALGBRC(150) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(235)+CONSTS(236))/ALGBRC(143) ** 2.00000
      ALGBRC(166) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(246)+CONSTS(243))/ALGBRC(160) ** 2.00000
      ALGBRC(296) = ( ALGBRC(9)*STATES(4)+ ALGBRC(20)*STATES(11)+ ALGBRC(32)*STATES(10)+ ALGBRC(42)*STATES(6)+ ALGBRC(53)*STATES(24)+ CONSTS(27)*STATES(23)+ ALGBRC(67)*STATES(3)+ CONSTS(28)*STATES(7)+ CONSTS(29)*STATES(8)+ ALGBRC(90)*STATES(9)+ CONSTS(30)*STATES(13)+ ALGBRC(116)*STATES(12)+ ALGBRC(106)*STATES(14)+ CONSTS(31)*STATES(17)+ CONSTS(32)*STATES(18)+ ALGBRC(139)*STATES(19)+ ALGBRC(150)*STATES(20)+ CONSTS(33)*STATES(21)+ ALGBRC(166)*STATES(22))/(ALGBRC(292)+ALGBRC(293))
      ALGBRC(11) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(150)+CONSTS(151)*log(10.0000))/ALGBRC(2) ** 2.00000
      ALGBRC(22) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(160)+CONSTS(161)*log(10.0000))/ALGBRC(13) ** 2.00000
      ALGBRC(34) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(170)+CONSTS(171)*log(10.0000))/ALGBRC(24) ** 2.00000
      ALGBRC(44) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(178)+CONSTS(179)*log(10.0000))/ALGBRC(36) ** 2.00000
      ALGBRC(55) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(187)+CONSTS(188)*log(10.0000))/ALGBRC(46) ** 2.00000
      ALGBRC(69) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(196)+CONSTS(197)*log(10.0000))/ALGBRC(61) ** 2.00000
      ALGBRC(92) = ( (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(208)*( ALGBRC(83)*10.0000 ** - ALGBRC(1)+CONSTS(206)*log(10.0000) -  10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)*2.00000*log(10.0000)))/ALGBRC(83) ** 2.00000
      ALGBRC(108) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(217)+CONSTS(214)*log(10.0000))/ALGBRC(100) ** 2.00000
      ALGBRC(118) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(220)+CONSTS(221)*log(10.0000))/ALGBRC(110) ** 2.00000
      ALGBRC(141) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(230)+CONSTS(231)*log(10.0000))/ALGBRC(132) ** 2.00000
      ALGBRC(152) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(235)+CONSTS(236)*log(10.0000))/ALGBRC(143) ** 2.00000
      ALGBRC(168) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(246)+CONSTS(243)*log(10.0000))/ALGBRC(160) ** 2.00000
      ALGBRC(10) = (( ALGBRC(2)*10.0000 ** CONSTS(151))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(151))/CONSTS(135))/ALGBRC(2) ** 2.00000
      ALGBRC(21) = (( ALGBRC(13)*10.0000 ** CONSTS(161))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(161))/CONSTS(135))/ALGBRC(13) ** 2.00000
      ALGBRC(33) = (( ALGBRC(24)*10.0000 ** CONSTS(171))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(171))/CONSTS(135))/ALGBRC(24) ** 2.00000
      ALGBRC(43) = (( ALGBRC(36)*10.0000 ** CONSTS(179))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(179))/CONSTS(135))/ALGBRC(36) ** 2.00000
      ALGBRC(54) = (( ALGBRC(46)*10.0000 ** CONSTS(188))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(188))/CONSTS(135))/ALGBRC(46) ** 2.00000
      ALGBRC(68) = (( ALGBRC(61)*10.0000 ** CONSTS(197))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(197))/CONSTS(135))/ALGBRC(61) ** 2.00000
      ALGBRC(91) = (( ALGBRC(83)*10.0000 ** CONSTS(208))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(208))/CONSTS(135))/ALGBRC(83) ** 2.00000
      ALGBRC(107) = (( ALGBRC(100)*10.0000 ** CONSTS(214))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(214))/CONSTS(135))/ALGBRC(100) ** 2.00000
      ALGBRC(117) = (( ALGBRC(110)*10.0000 ** CONSTS(221))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(221))/CONSTS(135))/ALGBRC(110) ** 2.00000
      ALGBRC(140) = (( ALGBRC(132)*10.0000 ** CONSTS(231))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(231))/CONSTS(135))/ALGBRC(132) ** 2.00000
      ALGBRC(151) = (( ALGBRC(143)*10.0000 ** CONSTS(236))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(236))/CONSTS(135))/ALGBRC(143) ** 2.00000
      ALGBRC(167) = (( ALGBRC(160)*10.0000 ** CONSTS(243))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(243))/CONSTS(135))/ALGBRC(160) ** 2.00000
      ALGBRC(297) = - 1.00000 - ( ALGBRC(21)*STATES(11)+ ALGBRC(33)*STATES(10)+ ALGBRC(43)*STATES(6)+ ALGBRC(10)*STATES(4)+ ALGBRC(54)*STATES(24)+ ALGBRC(68)*STATES(3)+ ALGBRC(91)*STATES(9)+ ALGBRC(107)*STATES(14)+ ALGBRC(117)*STATES(12)+ ALGBRC(151)*STATES(20)+ ALGBRC(140)*STATES(19)+ ALGBRC(167)*STATES(22))
      ALGBRC(298) = ( ALGBRC(22)*STATES(11)+ ALGBRC(34)*STATES(10)+ ALGBRC(44)*STATES(6)+ ALGBRC(11)*STATES(4)+ ALGBRC(55)*STATES(24)+ ALGBRC(69)*STATES(3)+ ALGBRC(92)*STATES(9)+ ALGBRC(108)*STATES(14)+ ALGBRC(118)*STATES(12)+ ALGBRC(152)*STATES(20)+ ALGBRC(141)*STATES(19)+ ALGBRC(168)*STATES(22))/ALGBRC(297)
      ALGBRC(6) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(151)*ALGBRC(3)
      ALGBRC(50) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(188)*ALGBRC(47)
      ALGBRC(64) =  (( ALGBRC(62)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(197)
      ALGBRC(87) =  (( ALGBRC(84)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(208)
      ALGBRC(103) =  (( ALGBRC(101)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(214)
      ALGBRC(113) =  (( ALGBRC(111)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(221)
      ALGBRC(136) =  (( ALGBRC(133)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(231)
      ALGBRC(147) =  (( ALGBRC(144)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(236)
      ALGBRC(163) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(243)*ALGBRC(161)
      ALGBRC(312) = ( ALGBRC(16)*ALGBRC(301)+ ALGBRC(27)*ALGBRC(302)+ ALGBRC(39)*ALGBRC(303)+ ALGBRC(6)*ALGBRC(304)+ ALGBRC(50)*ALGBRC(300)+ ALGBRC(64)*ALGBRC(305)+ ALGBRC(87)*ALGBRC(306)+ ALGBRC(103)*ALGBRC(308)+ ALGBRC(113)*ALGBRC(307)+ ALGBRC(147)*ALGBRC(310)+ ALGBRC(136)*ALGBRC(309)+ ALGBRC(163)*ALGBRC(311))/ALGBRC(297)
      RATES(2) = TERNRY(VOI.LE.1.00000, ( CONSTS(34)*(ALGBRC(312)+ ALGBRC(298)*ALGBRC(295)))/(1.00000 -  ALGBRC(298)*ALGBRC(296)),  CONSTS(34)*ALGBRC(312))
      RATES(1) = ( CONSTS(342)*(ALGBRC(295)+ ALGBRC(312)*ALGBRC(296)))/(1.00000 -  ALGBRC(298)*ALGBRC(296))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(221) = TERNRY(CONSTS(39) .NE. 45.0000, CONSTS(6), TERNRY(VOI.LT.40.0000, 0.00100000, TERNRY(VOI.GE.40.0000.AND.VOI.LT.80.0000, 0.00400000, TERNRY(VOI.GE.80.0000.AND.VOI.LT.100.000, 0.0100000, TERNRY(VOI.GE.100.000, 0.0400000, 0.0/0.0)
      ALGBRC(222) = 1.00000+STATES(5)/CONSTS(40)+STATES(4)/CONSTS(265)+( STATES(5)*STATES(4))/( CONSTS(40)*CONSTS(41))+STATES(5)/CONSTS(267)+STATES(3)/CONSTS(42)+( STATES(5)*STATES(3))/( CONSTS(268)*CONSTS(267))
      ALGBRC(1) = TERNRY(VOI.LE.1.00000.OR.VOI.GT.1.00000.AND.CONSTS(35).EQ.0.00000, STATES(1), CONSTS(36))
      ALGBRC(223) = 1.40400/(1.00000+10.0000 ** 5.94000 - ALGBRC(1)+10.0000 ** ALGBRC(1) - 7.29000)
      ALGBRC(225) = (( ALGBRC(223)*CONSTS(264)*STATES(4))/( CONSTS(266)*CONSTS(265)))/ALGBRC(222)
      ALGBRC(2) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(150)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(151)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(145)
      ALGBRC(61) = 1.00000+ 10.0000 ** - ALGBRC(1)*10.0000 ** CONSTS(196)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(197)
      ALGBRC(12) = (CONSTS(153)+ CONSTS(152)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(152))
      ALGBRC(70) = (CONSTS(199)+ CONSTS(198)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(198))
      ALGBRC(170) = 655.700+ CONSTS(249)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)
      ALGBRC(183) = (ALGBRC(170)+ALGBRC(70)) - ALGBRC(12)
      ALGBRC(184) = ( EXP(- ALGBRC(183)/( CONSTS(1)*CONSTS(131)))*ALGBRC(61))/ALGBRC(2)
      ALGBRC(224) = ( ALGBRC(223)*CONSTS(264)*CONSTS(267)*CONSTS(268))/( CONSTS(266)*CONSTS(265)*ALGBRC(184))
      ALGBRC(226) = (( ALGBRC(224)*STATES(5))/( CONSTS(267)*CONSTS(268)))/ALGBRC(222)
      ALGBRC(227) =  ALGBRC(221)*( STATES(5)*ALGBRC(225) -  STATES(3)*ALGBRC(226))
      ALGBRC(228) = 1.00000 - ALGBRC(221)
      ALGBRC(229) = (STATES(6)/CONSTS(45) ** CONSTS(47)/CONSTS(46))/(1.00000+STATES(6)/CONSTS(45) ** CONSTS(47)/CONSTS(46))
      ALGBRC(230) = 1.00000+STATES(5)/CONSTS(270)+STATES(4)/CONSTS(43)+STATES(5)/CONSTS(272)+STATES(3)/CONSTS(44)+( STATES(5)*STATES(4))/( CONSTS(270)*CONSTS(269))+( STATES(5)*STATES(3))/( CONSTS(271)*CONSTS(272))
      ALGBRC(231) = 1.75000/(1.00000+10.0000 ** 6.12000 - ALGBRC(1)+10.0000 ** ALGBRC(1) - 7.03000)
      ALGBRC(233) = (( ALGBRC(231)*ALGBRC(229)*CONSTS(264)*STATES(4))/( CONSTS(270)*CONSTS(269)))/ALGBRC(230)
      ALGBRC(232) = ( ALGBRC(231)*CONSTS(264)*CONSTS(271)*CONSTS(272))/( CONSTS(270)*CONSTS(269)*ALGBRC(184))
      ALGBRC(234) = (( ALGBRC(229)*ALGBRC(232)*STATES(5))/( CONSTS(271)*CONSTS(272)))/ALGBRC(230)
      ALGBRC(235) =  ALGBRC(228)*( STATES(5)*ALGBRC(233) -  STATES(3)*ALGBRC(234))
      ALGBRC(237) = - (ALGBRC(227)+ALGBRC(235))
      ALGBRC(236) = ( CONSTS(273)*1.32900)/(1.00000+10.0000 ** - ALGBRC(1)+6.64000+10.0000 ** ALGBRC(1) - 8.36000)
      ALGBRC(71) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(200)
      ALGBRC(76) = (CONSTS(202)+ CONSTS(201)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(201))
      ALGBRC(186) = ALGBRC(76) - ALGBRC(70)
      ALGBRC(187) = ( EXP(- ALGBRC(186)/( CONSTS(1)*CONSTS(131)))*ALGBRC(71))/ALGBRC(61)
      ALGBRC(238) = ( ALGBRC(236)*CONSTS(275))/( CONSTS(274)*ALGBRC(187))
      ALGBRC(239) = (( ALGBRC(236)*STATES(3))/CONSTS(274) - ( ALGBRC(238)*STATES(7))/CONSTS(275))/(1.00000+STATES(3)/CONSTS(274)+STATES(7)/CONSTS(275))
      ALGBRC(240) = CONSTS(276)/(1.00000+10.0000 ** - ALGBRC(1)+6.94000+10.0000 ** ALGBRC(1) - 9.35000)
      ALGBRC(77) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(203)
      ALGBRC(82) = (CONSTS(205)+ CONSTS(204)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(204))
      ALGBRC(189) = ALGBRC(82) - ALGBRC(76)
      ALGBRC(190) = ( EXP(- ALGBRC(189)/( CONSTS(1)*CONSTS(131)))*ALGBRC(77))/ALGBRC(71)
      ALGBRC(241) =  (( ALGBRC(240)*CONSTS(277))/CONSTS(278))*ALGBRC(190)
      ALGBRC(242) = (( ALGBRC(241)*STATES(7))/CONSTS(277) - ( ALGBRC(240)*STATES(8))/CONSTS(278))/(1.00000+STATES(8)/CONSTS(278)+STATES(7)/CONSTS(277))
      ALGBRC(244) = ALGBRC(239) - ALGBRC(242)
      ALGBRC(243) = CONSTS(279)/(1.00000+ALGBRC(1)/6.80000 ** - 30.0000)
      ALGBRC(13) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(160)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(161)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(162)
      ALGBRC(24) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(170)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(171)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(165)
      ALGBRC(83) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(206)+10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(208)
      ALGBRC(23) = (CONSTS(164)+ CONSTS(163)*CONSTS(1)*CONSTS(131)*log(10.0000)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(163))
      ALGBRC(35) = (CONSTS(173)+ CONSTS(172)*CONSTS(1)*CONSTS(131)*log(10.0000)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(172))
      ALGBRC(93) = (CONSTS(210)+ CONSTS(209)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(209))
      ALGBRC(174) = CONSTS(257)+ CONSTS(256)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)
      ALGBRC(192) = ((ALGBRC(93)+ALGBRC(35)+ALGBRC(174)) - ALGBRC(82)) - ALGBRC(23)
      ALGBRC(193) = ( EXP(- ALGBRC(192)/( CONSTS(1)*CONSTS(131)))*ALGBRC(83)*ALGBRC(24))/( ALGBRC(77)*ALGBRC(13)*10.0000 ** - ALGBRC(1))
      ALGBRC(245) = ( ALGBRC(243)*CONSTS(283)*CONSTS(285))/( CONSTS(280)*CONSTS(282)*ALGBRC(193))
      ALGBRC(246) =  CONSTS(54)*( ((1.00000+STATES(11)/CONSTS(50))/(1.00000+( CONSTS(52)*STATES(11))/CONSTS(50)))*(1.00000+( CONSTS(53)*STATES(6))/CONSTS(51)))/(1.00000+STATES(6)/CONSTS(51)) ** 4.00000
      ALGBRC(247) =  (1.00000+STATES(8)/CONSTS(280))*(1.00000+STATES(11)/CONSTS(282))+STATES(9)/CONSTS(283)+ (STATES(10)/CONSTS(285))*(1.00000+STATES(9)/CONSTS(283))
      ALGBRC(248) =  (1.00000+STATES(8)/CONSTS(281))*(1.00000+STATES(11)/CONSTS(284))+STATES(9)/CONSTS(48)+ (STATES(10)/CONSTS(49))*(1.00000+STATES(9)/CONSTS(48))
      ALGBRC(249) = ( ((( ALGBRC(243)*STATES(8)*STATES(11))/( CONSTS(280)*CONSTS(282)))/ALGBRC(247) - (( ALGBRC(245)*STATES(10)*STATES(9))/( CONSTS(285)*CONSTS(283)))/ALGBRC(247))*(1.00000+ CONSTS(286)*ALGBRC(246)*ALGBRC(248)/ALGBRC(247) ** 3.00000))/(1.00000+ ALGBRC(246)*ALGBRC(248)/ALGBRC(247) ** 4.00000)
      ALGBRC(251) = ALGBRC(242) - ALGBRC(249)
      ALGBRC(250) = ( CONSTS(287)*1.01300)/(1.00000+10.0000 ** - ALGBRC(1)+5.32000+10.0000 ** ALGBRC(1) - 9.15000)
      ALGBRC(94) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(211)
      ALGBRC(110) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(220)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(221)
      ALGBRC(99) = (CONSTS(213)+ CONSTS(212)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(212))
      ALGBRC(119) = (CONSTS(223)+ CONSTS(222)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(222))
      ALGBRC(195) = (ALGBRC(119)+ALGBRC(99)) - ALGBRC(93)
      ALGBRC(196) = ( 1.00000*EXP(- ALGBRC(195)/( CONSTS(1)*CONSTS(131)))*ALGBRC(94)*ALGBRC(83))/ALGBRC(110)
      ALGBRC(252) = ( ALGBRC(250)*CONSTS(290)*CONSTS(289))/( CONSTS(288)*ALGBRC(196))
      ALGBRC(253) = (( ALGBRC(250)*STATES(9))/CONSTS(288) - ( ALGBRC(252)*STATES(13)*STATES(12))/( CONSTS(290)*CONSTS(289)))/(1.00000+STATES(9)/CONSTS(288)+STATES(13)/CONSTS(290)+STATES(12)/CONSTS(289))
      ALGBRC(198) = ALGBRC(119) - ALGBRC(99)
      ALGBRC(199) = ( EXP(- ALGBRC(198)/( CONSTS(1)*CONSTS(131)))*ALGBRC(110))/ALGBRC(94)
      ALGBRC(254) = ( CONSTS(294)*CONSTS(293))/( CONSTS(292)*ALGBRC(199))
      ALGBRC(255) = (( CONSTS(294)*STATES(13))/CONSTS(292) - ( ALGBRC(254)*STATES(12))/CONSTS(293))/(1.00000+STATES(13)/CONSTS(292)+STATES(12)/CONSTS(293))
      ALGBRC(259) = 1.00000+STATES(4)/CONSTS(304)+STATES(13)/CONSTS(302)+STATES(15)/CONSTS(303)+( STATES(13)*STATES(15))/( CONSTS(302)*CONSTS(303))+( STATES(13)*STATES(15)*STATES(4))/( CONSTS(302)*CONSTS(303)*CONSTS(304))+STATES(17)/CONSTS(305)+STATES(16)/CONSTS(306)+( STATES(17)*STATES(16))/( CONSTS(306)*CONSTS(305))
      ALGBRC(260) =  CONSTS(301)*0.000700000*EXP( ALGBRC(1)*0.897900)
      ALGBRC(120) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(224)
      ALGBRC(125) = (CONSTS(226)+ CONSTS(225)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(16.0000 - CONSTS(225))
      ALGBRC(171) = (CONSTS(251)+ CONSTS(250)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(250))
      ALGBRC(172) = (CONSTS(253)+ CONSTS(252)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(252))
      ALGBRC(201) = (((ALGBRC(125)+ALGBRC(172)+ALGBRC(174)) - ALGBRC(12)) - ALGBRC(99)) - ALGBRC(171)
      ALGBRC(202) = ( EXP(- ALGBRC(201)/( CONSTS(1)*CONSTS(131)))*ALGBRC(120))/( ALGBRC(2)*ALGBRC(94)*10.0000 ** - ALGBRC(1)*1.00000)
      ALGBRC(261) = ( ALGBRC(260)*CONSTS(305)*CONSTS(306))/( CONSTS(302)*CONSTS(304)*CONSTS(303)*ALGBRC(202))
      ALGBRC(262) = (( ALGBRC(260)*STATES(13)*STATES(15)*STATES(4))/( CONSTS(303)*CONSTS(302)*CONSTS(304)) - ( ALGBRC(261)*STATES(17)*STATES(16))/( CONSTS(305)*CONSTS(306)))/ALGBRC(259)
      ALGBRC(264) = (ALGBRC(253) - ALGBRC(255)) - ALGBRC(262)
      ALGBRC(126) = 1.00000+10.0000 ** - ALGBRC(1)+6.21000
      ALGBRC(131) = (CONSTS(229)+ CONSTS(228)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(228))
      ALGBRC(207) = ((ALGBRC(23)+ALGBRC(131)) - ALGBRC(125)) - ALGBRC(35)
      ALGBRC(208) = ( EXP(- ALGBRC(207)/( CONSTS(1)*CONSTS(131)))*ALGBRC(13)*ALGBRC(126))/( ALGBRC(120)*ALGBRC(24))
      ALGBRC(263) =  (( CONSTS(312)*CONSTS(308)*CONSTS(309))/( CONSTS(310)*CONSTS(311)))*ALGBRC(208)
      ALGBRC(265) = 1.00000+STATES(10)/CONSTS(309)+STATES(17)/CONSTS(308)+( STATES(17)*STATES(10))/( CONSTS(308)*CONSTS(309))+STATES(18)/CONSTS(310)+STATES(11)/CONSTS(311)+( STATES(18)*STATES(11))/( CONSTS(310)*CONSTS(311))
      ALGBRC(266) = (( ALGBRC(263)*STATES(17)*STATES(10))/( CONSTS(309)*CONSTS(308)) - ( CONSTS(312)*STATES(11)*STATES(18))/( CONSTS(311)*CONSTS(310)))/ALGBRC(265)
      ALGBRC(268) = ALGBRC(262) - ALGBRC(266)
      ALGBRC(267) = ( CONSTS(313)*0.989000)/(1.00000+10.0000 ** - ALGBRC(1)+5.62000+10.0000 ** ALGBRC(1) - 8.74000)
      ALGBRC(132) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(230)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(231)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(232)
      ALGBRC(142) = (CONSTS(234)+ CONSTS(233)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(233))
      ALGBRC(210) = ALGBRC(142) - ALGBRC(131)
      ALGBRC(211) = ( EXP(- ALGBRC(210)/( CONSTS(1)*CONSTS(131)))*ALGBRC(132))/ALGBRC(126)
      ALGBRC(269) = ( ALGBRC(267)*CONSTS(315))/( CONSTS(314)*ALGBRC(211))
      ALGBRC(270) = (( ALGBRC(267)*STATES(18))/CONSTS(314) - ( ALGBRC(269)*STATES(19))/CONSTS(315))/(1.00000+STATES(18)/CONSTS(314)+STATES(19)/CONSTS(315))
      ALGBRC(272) = ALGBRC(266) - ALGBRC(270)
      ALGBRC(256) =  (1.00000+STATES(14)/CONSTS(296)+STATES(16)/CONSTS(299))*(1.00000+STATES(12)/CONSTS(298)+STATES(15)/CONSTS(299))
      ALGBRC(100) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(217)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(214)
      ALGBRC(109) = (CONSTS(219)+ CONSTS(218)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(218))
      ALGBRC(204) = ((ALGBRC(174)+ALGBRC(172)+ALGBRC(119)) - ALGBRC(171)) - ALGBRC(109)
      ALGBRC(205) = ( EXP(- ALGBRC(204)/( CONSTS(1)*CONSTS(131)))*ALGBRC(110))/( ALGBRC(100)*10.0000 ** - ALGBRC(1))
      ALGBRC(257) = ( CONSTS(300)*CONSTS(296)*CONSTS(297)*ALGBRC(205))/( CONSTS(298)*CONSTS(299))
      ALGBRC(258) = (( ALGBRC(257)*STATES(14)*STATES(15))/( CONSTS(296)*CONSTS(297)) - ( CONSTS(300)*STATES(12)*STATES(16))/( CONSTS(298)*CONSTS(299)))/ALGBRC(256)
      ALGBRC(277) =  CONSTS(325)*( - 0.113400*ALGBRC(1)+1.60690)
      ALGBRC(154) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(240)
      ALGBRC(160) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(246)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(243)
      ALGBRC(159) = (CONSTS(242)+ CONSTS(241)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(241))
      ALGBRC(169) = (CONSTS(248)+ CONSTS(247)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(1.00000 - CONSTS(247))
      ALGBRC(219) = (((ALGBRC(169)+ALGBRC(171)) - ALGBRC(159)) - ALGBRC(172)) - ALGBRC(174)
      ALGBRC(220) = ( EXP(- ALGBRC(219)/( CONSTS(1)*CONSTS(131)))*ALGBRC(160)*10.0000 ** - ALGBRC(1))/ALGBRC(154)
      ALGBRC(278) = ( ALGBRC(277)*CONSTS(328)*CONSTS(329))/( CONSTS(326)*CONSTS(327)*ALGBRC(220))
      ALGBRC(279) = (( ALGBRC(277)*STATES(21)*STATES(16))/( CONSTS(326)*CONSTS(327)) - ( ALGBRC(278)*STATES(22)*STATES(15))/( CONSTS(328)*CONSTS(329)))/(1.00000+STATES(21)/CONSTS(326)+STATES(16)/CONSTS(327)+( STATES(21)*STATES(16))/( CONSTS(326)*CONSTS(327))+STATES(22)/CONSTS(328)+STATES(15)/CONSTS(329)+( STATES(22)*STATES(15))/( CONSTS(328)*CONSTS(329)))
      ALGBRC(282) = (- ALGBRC(262) - ALGBRC(258))+ALGBRC(279)
      ALGBRC(283) = (ALGBRC(262)+ALGBRC(258)) - ALGBRC(279)
      ALGBRC(274) = ( CONSTS(320)*1.05000)/(1.00000+10.0000 ** - ALGBRC(1)+5.58000+10.0000 ** ALGBRC(1) - 8.79000)
      ALGBRC(143) = 1.00000+10.0000 ** CONSTS(235) - ALGBRC(1)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(236)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(237)
      ALGBRC(153) = (CONSTS(239)+ CONSTS(238)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(9.00000 - CONSTS(238))
      ALGBRC(216) = (((ALGBRC(159)+ALGBRC(23)) - ALGBRC(174)) - ALGBRC(153)) - ALGBRC(35)
      ALGBRC(217) = ( EXP(- ALGBRC(216)/( CONSTS(1)*CONSTS(131)))*ALGBRC(154)*ALGBRC(13)*10.0000 ** - ALGBRC(1))/( ALGBRC(143)*ALGBRC(24))
      ALGBRC(275) = ( ALGBRC(274)*CONSTS(323)*CONSTS(324))/( CONSTS(321)*CONSTS(322)*ALGBRC(217))
      ALGBRC(276) = (( ALGBRC(274)*STATES(20)*STATES(10))/( CONSTS(321)*CONSTS(322)) - ( ALGBRC(275)*STATES(21)*STATES(11))/( CONSTS(323)*CONSTS(324)))/(1.00000+STATES(20)/CONSTS(321)+STATES(10)/CONSTS(322)+( STATES(20)*STATES(10))/( CONSTS(321)*CONSTS(322))+STATES(11)/CONSTS(324)+STATES(21)/CONSTS(323)+( STATES(21)*STATES(11))/( CONSTS(323)*CONSTS(324)))
      ALGBRC(284) = ALGBRC(276) - ALGBRC(279)
      ALGBRC(46) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(187)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(188)+ (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(182)
      ALGBRC(56) = 1.00000+10.0000 ** - ALGBRC(1)+CONSTS(190)
      ALGBRC(177) = ( EXP(- CONSTS(263)/( CONSTS(1)*CONSTS(131)))*10.0000 ** - ALGBRC(1)*ALGBRC(13)*ALGBRC(56))/( ALGBRC(46)*ALGBRC(24))
      ALGBRC(286) = ( (CONSTS(330)/ALGBRC(177))*CONSTS(332)*CONSTS(334))/( CONSTS(333)*CONSTS(331))
      ALGBRC(287) = (( ALGBRC(286)*STATES(11)*STATES(23))/( CONSTS(332)*CONSTS(334)) - ( CONSTS(330)*STATES(10)*STATES(24))/( CONSTS(333)*CONSTS(331)))/(1.00000+STATES(10)/CONSTS(333)+STATES(24)/CONSTS(56)+( STATES(24)*STATES(10))/( CONSTS(333)*CONSTS(331))+STATES(11)/CONSTS(332)+( STATES(23)*STATES(11))/( CONSTS(334)*CONSTS(332)))
      ALGBRC(290) = - ALGBRC(287)
      ALGBRC(14) = 1.00000/ALGBRC(13)
      ALGBRC(15) =  10.0000 ** - ALGBRC(1)+CONSTS(160)*ALGBRC(14)
      ALGBRC(16) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(161)*ALGBRC(14)
      ALGBRC(17) =  (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(162)*ALGBRC(14)
      ALGBRC(18) =  0.00000*ALGBRC(14)+ 1.00000*ALGBRC(15)+ 0.00000*ALGBRC(16)+ 0.00000*ALGBRC(17)
      ALGBRC(25) = 1.00000/ALGBRC(24)
      ALGBRC(26) =  10.0000 ** - ALGBRC(1)+CONSTS(170)*ALGBRC(25)
      ALGBRC(27) =  (( ALGBRC(25)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(171)
      ALGBRC(28) =  (( ALGBRC(25)*CONSTS(134))/CONSTS(135))*10.0000 ** CONSTS(165)
      ALGBRC(29) =  0.00000*ALGBRC(25)+ 1.00000*ALGBRC(26)+ 0.00000*ALGBRC(27)+ 0.00000*ALGBRC(28)
      ALGBRC(47) = 1.00000/ALGBRC(46)
      ALGBRC(49) =  10.0000 ** - ALGBRC(1)+CONSTS(187)*ALGBRC(47)
      ALGBRC(51) = ALGBRC(49)
      ALGBRC(57) = 1.00000/ALGBRC(56)
      ALGBRC(58) =  ALGBRC(57)*10.0000 ** - ALGBRC(1)+CONSTS(190)
      ALGBRC(59) =  0.00000*ALGBRC(57)+ 1.00000*ALGBRC(58)
      ALGBRC(175) = (((ALGBRC(18)+ALGBRC(59)) - ALGBRC(51)) - ALGBRC(29))+(((CONSTS(163)+CONSTS(191)) - CONSTS(189)) - CONSTS(172))
      ALGBRC(36) = 1.00000+10.0000 ** CONSTS(178) - ALGBRC(1)+ (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(179)
      ALGBRC(37) = 1.00000/ALGBRC(36)
      ALGBRC(38) =  ALGBRC(37)*10.0000 ** CONSTS(178) - ALGBRC(1)
      ALGBRC(39) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(179)*ALGBRC(37)
      ALGBRC(40) =  0.00000*ALGBRC(37)+ALGBRC(38)+ 0.00000*ALGBRC(39)
      ALGBRC(179) = (( 2.00000*ALGBRC(29) - ALGBRC(18)) - ALGBRC(40))+(( 2.00000*CONSTS(172) - CONSTS(163)) - CONSTS(180))
      ALGBRC(3) = 1.00000/ALGBRC(2)
      ALGBRC(4) =  10.0000 ** - ALGBRC(1)+CONSTS(150)*ALGBRC(3)
      ALGBRC(7) =  1.00000*ALGBRC(4)
      ALGBRC(62) = 1.00000/ALGBRC(61)
      ALGBRC(63) =  ALGBRC(62)*10.0000 ** - ALGBRC(1)+CONSTS(196)
      ALGBRC(65) = ALGBRC(63)
      ALGBRC(182) = (ALGBRC(65) - ALGBRC(7))+(1.00000 - CONSTS(152))
      ALGBRC(72) = 1.00000/ALGBRC(71)
      ALGBRC(73) =  ALGBRC(72)*10.0000 ** - ALGBRC(1)+CONSTS(200)
      ALGBRC(74) = ALGBRC(73)
      ALGBRC(185) = (ALGBRC(74) - ALGBRC(65))+(CONSTS(201) - CONSTS(198))
      ALGBRC(78) = 1.00000/ALGBRC(77)
      ALGBRC(79) =  ALGBRC(78)*10.0000 ** - ALGBRC(1)+CONSTS(203)
      ALGBRC(80) = ALGBRC(79)
      ALGBRC(188) = (ALGBRC(80) - ALGBRC(74))+(CONSTS(204) - CONSTS(201))
      ALGBRC(84) = 1.00000/ALGBRC(83)
      ALGBRC(85) =  ALGBRC(84)*10.0000 ** - ALGBRC(1)+CONSTS(206)
      ALGBRC(86) =  ALGBRC(84)*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)
      ALGBRC(88) = ALGBRC(85)+ 2.00000*ALGBRC(86)
      ALGBRC(191) = (((((ALGBRC(29)+ALGBRC(88)) - ALGBRC(80)) - ALGBRC(18))+CONSTS(172)+CONSTS(209)) - CONSTS(204)) - CONSTS(163)
      ALGBRC(95) = 1.00000/ALGBRC(94)
      ALGBRC(96) =  ALGBRC(95)*10.0000 ** - ALGBRC(1)+CONSTS(211)
      ALGBRC(97) = ALGBRC(96)
      ALGBRC(111) = 1.00000/ALGBRC(110)
      ALGBRC(112) =  ALGBRC(111)*10.0000 ** - ALGBRC(1)+CONSTS(220)
      ALGBRC(114) = ALGBRC(112)
      ALGBRC(194) = ((ALGBRC(114)+ALGBRC(97)) - ALGBRC(88))+((CONSTS(222)+CONSTS(212)) - CONSTS(209))
      ALGBRC(197) = (ALGBRC(114) - ALGBRC(97))+(CONSTS(222) - CONSTS(212))
      ALGBRC(121) = 1.00000/ALGBRC(120)
      ALGBRC(122) =  ALGBRC(121)*10.0000 ** - ALGBRC(1)+CONSTS(224)
      ALGBRC(123) = ALGBRC(122)
      ALGBRC(200) = ((ALGBRC(123) - ALGBRC(97)) - ALGBRC(7))+((((CONSTS(225)+CONSTS(252)) - CONSTS(212)) - CONSTS(152)) - CONSTS(250))
      ALGBRC(101) = 1.00000/ALGBRC(100)
      ALGBRC(102) =  ALGBRC(101)*10.0000 ** - ALGBRC(1)+CONSTS(217)
      ALGBRC(104) = ALGBRC(102)
      ALGBRC(203) = (ALGBRC(114) - ALGBRC(104))+(((CONSTS(222)+CONSTS(252)) - CONSTS(250)) - CONSTS(218))
      ALGBRC(127) = 1.00000/ALGBRC(126)
      ALGBRC(128) =  ALGBRC(127)*10.0000 ** - ALGBRC(1)+6.21000
      ALGBRC(129) = ALGBRC(128)
      ALGBRC(206) = (((ALGBRC(129)+ALGBRC(18)) - ALGBRC(123)) - ALGBRC(29))+(((CONSTS(228)+CONSTS(163)) - CONSTS(225)) - CONSTS(172))
      ALGBRC(133) = 1.00000/ALGBRC(132)
      ALGBRC(134) =  ALGBRC(133)*10.0000 ** - ALGBRC(1)+CONSTS(230)
      ALGBRC(137) = ALGBRC(134)
      ALGBRC(209) = (ALGBRC(137) - ALGBRC(129))+(CONSTS(233) - CONSTS(228))
      ALGBRC(144) = 1.00000/ALGBRC(143)
      ALGBRC(145) =  ALGBRC(144)*10.0000 ** CONSTS(235) - ALGBRC(1)
      ALGBRC(148) = ALGBRC(145)
      ALGBRC(212) = (ALGBRC(148) - ALGBRC(137))+((CONSTS(254)+CONSTS(238)) - CONSTS(233))
      ALGBRC(155) = 1.00000/ALGBRC(154)
      ALGBRC(156) =  ALGBRC(155)*10.0000 ** - ALGBRC(1)+CONSTS(240)
      ALGBRC(157) = ALGBRC(156)
      ALGBRC(215) = (((ALGBRC(157)+ALGBRC(18)) - ALGBRC(148)) - ALGBRC(29))+(((CONSTS(241)+CONSTS(163)) - CONSTS(238)) - CONSTS(172))
      ALGBRC(161) = 1.00000/ALGBRC(160)
      ALGBRC(162) =  ALGBRC(161)*10.0000 ** - ALGBRC(1)+CONSTS(246)
      ALGBRC(164) = ALGBRC(162)
      ALGBRC(218) = (ALGBRC(164) - ALGBRC(157))+(((CONSTS(247)+CONSTS(250)) - CONSTS(241)) - CONSTS(252))
      ALGBRC(173) = (CONSTS(255)+ CONSTS(254)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(0.00000 - CONSTS(254))
      ALGBRC(213) = (ALGBRC(173)+ALGBRC(153)) - ALGBRC(142)
      ALGBRC(214) = ( EXP(- ALGBRC(213)/( CONSTS(1)*CONSTS(131)))*ALGBRC(143))/ALGBRC(132)
      ALGBRC(271) = ( CONSTS(319)*CONSTS(318))/( CONSTS(317)*ALGBRC(214))
      ALGBRC(273) = (( CONSTS(319)*STATES(19))/CONSTS(317) - ( ALGBRC(271)*STATES(20))/CONSTS(318))/(1.00000+STATES(20)/CONSTS(318)+STATES(19)/CONSTS(317))
      ALGBRC(45) = (CONSTS(181)+ CONSTS(180)*log(10.0000)*CONSTS(1)*CONSTS(131)*ALGBRC(1)) -  CONSTS(144)*(4.00000 - CONSTS(180))
      ALGBRC(180) = ( 2.00000*ALGBRC(35) - ALGBRC(23)) - ALGBRC(45)
      ALGBRC(181) = ( EXP(- ALGBRC(180)/( CONSTS(1)*CONSTS(131)))*ALGBRC(24) ** 2.00000)/( ALGBRC(13)*ALGBRC(36))
      ALGBRC(288) = ( CONSTS(335)*CONSTS(338) ** 2.00000)/( CONSTS(336)*CONSTS(337)*ALGBRC(181))
      ALGBRC(291) = (( CONSTS(335)*STATES(11)*STATES(6))/( CONSTS(337)*CONSTS(336)) -  ALGBRC(288)*STATES(10)/CONSTS(338) ** 2.00000)/(1.00000+STATES(11)/CONSTS(337)+STATES(6)/CONSTS(336)+( STATES(11)*STATES(6))/( CONSTS(337)*CONSTS(336))+( 2.00000*STATES(10))/CONSTS(338)+STATES(10) ** 2.00000/CONSTS(338) ** 2.00000)
      ALGBRC(294) =  ALGBRC(175)*- ALGBRC(287)+ ALGBRC(179)*ALGBRC(291)+ ALGBRC(185)*ALGBRC(239)+ ALGBRC(182)*(ALGBRC(227)+ALGBRC(235))+ ALGBRC(188)*ALGBRC(242)+ ALGBRC(191)*ALGBRC(249)+ ALGBRC(194)*ALGBRC(253)+ ALGBRC(197)*ALGBRC(255)+ ALGBRC(200)*ALGBRC(262)+ ALGBRC(206)*ALGBRC(266)+ ALGBRC(209)*ALGBRC(270)+ ALGBRC(212)*ALGBRC(273)+ ALGBRC(215)*ALGBRC(276)+ ALGBRC(218)*ALGBRC(279)+ ALGBRC(203)*ALGBRC(258)
      ALGBRC(300) = ALGBRC(287)
      ALGBRC(280) = TERNRY(CONSTS(39) .NE. 5.00000.AND.CONSTS(39) .NE. 45.0000, CONSTS(23), TERNRY(CONSTS(39) .NE. 5.00000.AND.CONSTS(39).EQ.45.0000.AND.VOI.GE.110.000, CONSTS(23), TERNRY(CONSTS(39).EQ.5.00000.AND.VOI.GT.30.0000, CONSTS(23), 0.00000)
      ALGBRC(285) = ( ALGBRC(280)*STATES(11))/(CONSTS(55)+STATES(11))
      ALGBRC(301) = (((- ALGBRC(287) - ALGBRC(291)) - ALGBRC(249))+ALGBRC(266)+ALGBRC(276)) - ALGBRC(285)
      ALGBRC(302) = (((ALGBRC(287)+ 2.00000*ALGBRC(291)+ALGBRC(249)) - ALGBRC(266)) - ALGBRC(276))+ALGBRC(285)
      ALGBRC(303) = - ALGBRC(291)
      ALGBRC(304) = (- (ALGBRC(227)+ALGBRC(235)) - ALGBRC(262))+ALGBRC(285)
      ALGBRC(305) = (ALGBRC(227)+ALGBRC(235)) - ALGBRC(239)
      ALGBRC(306) = ALGBRC(249) - ALGBRC(253)
      ALGBRC(307) = ALGBRC(253)+ALGBRC(255)+ALGBRC(258)
      ALGBRC(308) = - ALGBRC(258)
      ALGBRC(309) = ALGBRC(270) - ALGBRC(273)
      ALGBRC(310) = ALGBRC(273) - ALGBRC(276)
      ALGBRC(311) = ALGBRC(279)
      ALGBRC(292) = ( log(10.0000)*CONSTS(339)*10.0000 ** - ALGBRC(1) - 6.87000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 6.87000 ** 2.00000+( log(10.0000)*CONSTS(340)*10.0000 ** - ALGBRC(1) - 8.30000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 8.30000 ** 2.00000+( log(10.0000)*CONSTS(341)*10.0000 ** - ALGBRC(1) - 4.80000)/10.0000 ** - ALGBRC(1)+10.0000 ** - 4.80000 ** 2.00000
      ALGBRC(8) = ( 10.0000 ** CONSTS(150)*(ALGBRC(2) - 10.0000 ** - ALGBRC(1)+CONSTS(150)))/( CONSTS(135)*ALGBRC(2) ** 2.00000)
      ALGBRC(19) = ( 10.0000 ** CONSTS(160)*(ALGBRC(13) - 10.0000 ** - ALGBRC(1)+CONSTS(160)))/( CONSTS(135)*ALGBRC(13) ** 2.00000)
      ALGBRC(30) = ( 10.0000 ** CONSTS(170)*(ALGBRC(24) - 10.0000 ** - ALGBRC(1)+CONSTS(170)))/( CONSTS(135)*ALGBRC(24) ** 2.00000)
      ALGBRC(41) = ( 10.0000 ** CONSTS(178)*(ALGBRC(36) - 10.0000 ** - ALGBRC(1)+CONSTS(178)))/( CONSTS(135)*ALGBRC(36) ** 2.00000)
      ALGBRC(52) = ( 10.0000 ** CONSTS(187)*(ALGBRC(46) - 10.0000 ** - ALGBRC(1)+CONSTS(187)))/( CONSTS(135)*ALGBRC(46) ** 2.00000)
      ALGBRC(60) = ( 10.0000 ** CONSTS(190)*(ALGBRC(56) - 10.0000 ** - ALGBRC(1)+CONSTS(190)))/( CONSTS(135)*ALGBRC(56) ** 2.00000)
      ALGBRC(66) = ( 10.0000 ** CONSTS(196)*(ALGBRC(61) - 10.0000 ** - ALGBRC(1)+CONSTS(196)))/( CONSTS(135)*ALGBRC(61) ** 2.00000)
      ALGBRC(75) = ( 10.0000 ** CONSTS(200)*(ALGBRC(71) - 10.0000 ** - ALGBRC(1)+CONSTS(200)))/( CONSTS(135)*ALGBRC(71) ** 2.00000)
      ALGBRC(81) = ( 10.0000 ** CONSTS(203)*(ALGBRC(77) - 10.0000 ** - ALGBRC(1)+CONSTS(203)))/( CONSTS(135)*ALGBRC(77) ** 2.00000)
      ALGBRC(89) = ( ALGBRC(83)*(10.0000 ** CONSTS(206)+ 2.00000*10.0000 ** - ALGBRC(1)+CONSTS(206)+CONSTS(207)) -  (10.0000 ** - ALGBRC(1)+CONSTS(206)+ 2.00000*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207))*(10.0000 ** CONSTS(206)+ 2.00000*10.0000 ** - ALGBRC(1)+CONSTS(206)+CONSTS(207)))/( CONSTS(135)*ALGBRC(83) ** 2.00000)
      ALGBRC(98) = ( 10.0000 ** CONSTS(211)*(ALGBRC(94) - 10.0000 ** - ALGBRC(1)+CONSTS(211)))/( CONSTS(135)*ALGBRC(94) ** 2.00000)
      ALGBRC(105) = ( 10.0000 ** CONSTS(217)*(ALGBRC(100) - 10.0000 ** - ALGBRC(1)+CONSTS(217)))/( CONSTS(135)*ALGBRC(100) ** 2.00000)
      ALGBRC(115) = ( 10.0000 ** CONSTS(220)*(ALGBRC(110) - 10.0000 ** - ALGBRC(1)+CONSTS(220)))/( CONSTS(135)*ALGBRC(110) ** 2.00000)
      ALGBRC(124) = ( 10.0000 ** CONSTS(224)*(ALGBRC(120) - 10.0000 ** - ALGBRC(1)+CONSTS(224)))/( CONSTS(135)*ALGBRC(120) ** 2.00000)
      ALGBRC(130) = ( 10.0000 ** CONSTS(227)*(ALGBRC(126) - 10.0000 ** - ALGBRC(1)+CONSTS(227)))/( CONSTS(135)*ALGBRC(126) ** 2.00000)
      ALGBRC(138) = ( 10.0000 ** CONSTS(230)*(ALGBRC(132) - 10.0000 ** - ALGBRC(1)+CONSTS(230)))/( CONSTS(135)*ALGBRC(132) ** 2.00000)
      ALGBRC(149) = ( 10.0000 ** CONSTS(235)*(ALGBRC(143) - 10.0000 ** - ALGBRC(1)+CONSTS(235)))/( CONSTS(135)*ALGBRC(143) ** 2.00000)
      ALGBRC(158) = ( 10.0000 ** CONSTS(240)*(ALGBRC(154) - 10.0000 ** - ALGBRC(1)+CONSTS(240)))/( CONSTS(135)*ALGBRC(154) ** 2.00000)
      ALGBRC(165) = ( 10.0000 ** CONSTS(246)*(ALGBRC(160) - 10.0000 ** - ALGBRC(1)+CONSTS(246)))/( CONSTS(135)*ALGBRC(160) ** 2.00000)
      ALGBRC(293) =  log(10.0000)*10.0000 ** - ALGBRC(1)*CONSTS(135)*(1.00000+ ALGBRC(8)*STATES(4)+ ALGBRC(19)*STATES(11)+ ALGBRC(30)*STATES(10)+ ALGBRC(41)*STATES(6)+ ALGBRC(52)*STATES(24)+ ALGBRC(60)*STATES(23)+ ALGBRC(66)*STATES(3)+ ALGBRC(75)*STATES(7)+ ALGBRC(81)*STATES(8)+ ALGBRC(89)*STATES(9)+ ALGBRC(98)*STATES(13)+ ALGBRC(115)*STATES(12)+ ALGBRC(105)*STATES(14)+ ALGBRC(124)*STATES(17)+ ALGBRC(130)*STATES(18)+ ALGBRC(138)*STATES(19)+ ALGBRC(149)*STATES(20)+ ALGBRC(158)*STATES(21)+ ALGBRC(165)*STATES(22))
      ALGBRC(295) = ALGBRC(294)/(ALGBRC(292)+ALGBRC(293))
      ALGBRC(9) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(150)+CONSTS(151))/ALGBRC(2) ** 2.00000
      ALGBRC(20) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(160)+CONSTS(161))/ALGBRC(13) ** 2.00000
      ALGBRC(32) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(170)+CONSTS(171))/ALGBRC(24) ** 2.00000
      ALGBRC(42) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(178)+CONSTS(179))/ALGBRC(36) ** 2.00000
      ALGBRC(53) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(187)+CONSTS(188))/ALGBRC(46) ** 2.00000
      ALGBRC(67) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(196)+CONSTS(197))/ALGBRC(61) ** 2.00000
      ALGBRC(90) = ( - (10.0000 ** - ALGBRC(1)+CONSTS(206)+ 2.00000*10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207))*10.0000 ** CONSTS(208))/( CONSTS(135)*ALGBRC(83) ** 2.00000)
      ALGBRC(106) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(217)+CONSTS(214))/ALGBRC(100) ** 2.00000
      ALGBRC(116) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(220)+CONSTS(221))/ALGBRC(110) ** 2.00000
      ALGBRC(139) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(230)+CONSTS(231))/ALGBRC(132) ** 2.00000
      ALGBRC(150) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(235)+CONSTS(236))/ALGBRC(143) ** 2.00000
      ALGBRC(166) = ( (- 10.0000 ** - ALGBRC(1)/CONSTS(135))*10.0000 ** CONSTS(246)+CONSTS(243))/ALGBRC(160) ** 2.00000
      ALGBRC(296) = ( ALGBRC(9)*STATES(4)+ ALGBRC(20)*STATES(11)+ ALGBRC(32)*STATES(10)+ ALGBRC(42)*STATES(6)+ ALGBRC(53)*STATES(24)+ CONSTS(27)*STATES(23)+ ALGBRC(67)*STATES(3)+ CONSTS(28)*STATES(7)+ CONSTS(29)*STATES(8)+ ALGBRC(90)*STATES(9)+ CONSTS(30)*STATES(13)+ ALGBRC(116)*STATES(12)+ ALGBRC(106)*STATES(14)+ CONSTS(31)*STATES(17)+ CONSTS(32)*STATES(18)+ ALGBRC(139)*STATES(19)+ ALGBRC(150)*STATES(20)+ CONSTS(33)*STATES(21)+ ALGBRC(166)*STATES(22))/(ALGBRC(292)+ALGBRC(293))
      ALGBRC(11) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(150)+CONSTS(151)*log(10.0000))/ALGBRC(2) ** 2.00000
      ALGBRC(22) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(160)+CONSTS(161)*log(10.0000))/ALGBRC(13) ** 2.00000
      ALGBRC(34) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(170)+CONSTS(171)*log(10.0000))/ALGBRC(24) ** 2.00000
      ALGBRC(44) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(178)+CONSTS(179)*log(10.0000))/ALGBRC(36) ** 2.00000
      ALGBRC(55) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(187)+CONSTS(188)*log(10.0000))/ALGBRC(46) ** 2.00000
      ALGBRC(69) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(196)+CONSTS(197)*log(10.0000))/ALGBRC(61) ** 2.00000
      ALGBRC(92) = ( (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(208)*( ALGBRC(83)*10.0000 ** - ALGBRC(1)+CONSTS(206)*log(10.0000) -  10.0000 **  - 2.00000*ALGBRC(1)+CONSTS(206)+CONSTS(207)*2.00000*log(10.0000)))/ALGBRC(83) ** 2.00000
      ALGBRC(108) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(217)+CONSTS(214)*log(10.0000))/ALGBRC(100) ** 2.00000
      ALGBRC(118) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(220)+CONSTS(221)*log(10.0000))/ALGBRC(110) ** 2.00000
      ALGBRC(141) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(230)+CONSTS(231)*log(10.0000))/ALGBRC(132) ** 2.00000
      ALGBRC(152) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(235)+CONSTS(236)*log(10.0000))/ALGBRC(143) ** 2.00000
      ALGBRC(168) = ( (STATES(2)/CONSTS(135))*10.0000 ** - ALGBRC(1)+CONSTS(246)+CONSTS(243)*log(10.0000))/ALGBRC(160) ** 2.00000
      ALGBRC(10) = (( ALGBRC(2)*10.0000 ** CONSTS(151))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(151))/CONSTS(135))/ALGBRC(2) ** 2.00000
      ALGBRC(21) = (( ALGBRC(13)*10.0000 ** CONSTS(161))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(161))/CONSTS(135))/ALGBRC(13) ** 2.00000
      ALGBRC(33) = (( ALGBRC(24)*10.0000 ** CONSTS(171))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(171))/CONSTS(135))/ALGBRC(24) ** 2.00000
      ALGBRC(43) = (( ALGBRC(36)*10.0000 ** CONSTS(179))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(179))/CONSTS(135))/ALGBRC(36) ** 2.00000
      ALGBRC(54) = (( ALGBRC(46)*10.0000 ** CONSTS(188))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(188))/CONSTS(135))/ALGBRC(46) ** 2.00000
      ALGBRC(68) = (( ALGBRC(61)*10.0000 ** CONSTS(197))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(197))/CONSTS(135))/ALGBRC(61) ** 2.00000
      ALGBRC(91) = (( ALGBRC(83)*10.0000 ** CONSTS(208))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(208))/CONSTS(135))/ALGBRC(83) ** 2.00000
      ALGBRC(107) = (( ALGBRC(100)*10.0000 ** CONSTS(214))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(214))/CONSTS(135))/ALGBRC(100) ** 2.00000
      ALGBRC(117) = (( ALGBRC(110)*10.0000 ** CONSTS(221))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(221))/CONSTS(135))/ALGBRC(110) ** 2.00000
      ALGBRC(140) = (( ALGBRC(132)*10.0000 ** CONSTS(231))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(231))/CONSTS(135))/ALGBRC(132) ** 2.00000
      ALGBRC(151) = (( ALGBRC(143)*10.0000 ** CONSTS(236))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(236))/CONSTS(135))/ALGBRC(143) ** 2.00000
      ALGBRC(167) = (( ALGBRC(160)*10.0000 ** CONSTS(243))/CONSTS(135) - ( (STATES(2)/CONSTS(135))*10.0000 **  2.00000*CONSTS(243))/CONSTS(135))/ALGBRC(160) ** 2.00000
      ALGBRC(297) = - 1.00000 - ( ALGBRC(21)*STATES(11)+ ALGBRC(33)*STATES(10)+ ALGBRC(43)*STATES(6)+ ALGBRC(10)*STATES(4)+ ALGBRC(54)*STATES(24)+ ALGBRC(68)*STATES(3)+ ALGBRC(91)*STATES(9)+ ALGBRC(107)*STATES(14)+ ALGBRC(117)*STATES(12)+ ALGBRC(151)*STATES(20)+ ALGBRC(140)*STATES(19)+ ALGBRC(167)*STATES(22))
      ALGBRC(298) = ( ALGBRC(22)*STATES(11)+ ALGBRC(34)*STATES(10)+ ALGBRC(44)*STATES(6)+ ALGBRC(11)*STATES(4)+ ALGBRC(55)*STATES(24)+ ALGBRC(69)*STATES(3)+ ALGBRC(92)*STATES(9)+ ALGBRC(108)*STATES(14)+ ALGBRC(118)*STATES(12)+ ALGBRC(152)*STATES(20)+ ALGBRC(141)*STATES(19)+ ALGBRC(168)*STATES(22))/ALGBRC(297)
      ALGBRC(6) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(151)*ALGBRC(3)
      ALGBRC(50) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(188)*ALGBRC(47)
      ALGBRC(64) =  (( ALGBRC(62)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(197)
      ALGBRC(87) =  (( ALGBRC(84)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(208)
      ALGBRC(103) =  (( ALGBRC(101)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(214)
      ALGBRC(113) =  (( ALGBRC(111)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(221)
      ALGBRC(136) =  (( ALGBRC(133)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(231)
      ALGBRC(147) =  (( ALGBRC(144)*STATES(2))/CONSTS(135))*10.0000 ** CONSTS(236)
      ALGBRC(163) =  (STATES(2)/CONSTS(135))*10.0000 ** CONSTS(243)*ALGBRC(161)
      ALGBRC(312) = ( ALGBRC(16)*ALGBRC(301)+ ALGBRC(27)*ALGBRC(302)+ ALGBRC(39)*ALGBRC(303)+ ALGBRC(6)*ALGBRC(304)+ ALGBRC(50)*ALGBRC(300)+ ALGBRC(64)*ALGBRC(305)+ ALGBRC(87)*ALGBRC(306)+ ALGBRC(103)*ALGBRC(308)+ ALGBRC(113)*ALGBRC(307)+ ALGBRC(147)*ALGBRC(310)+ ALGBRC(136)*ALGBRC(309)+ ALGBRC(163)*ALGBRC(311))/ALGBRC(297)
      ALGBRC(5) =  (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(145)*ALGBRC(3)
      ALGBRC(31) = ((ALGBRC(29)+ALGBRC(7)) - ALGBRC(18))+(((CONSTS(172)+CONSTS(152)) - CONSTS(163)) - CONSTS(254))
      ALGBRC(48) =  (CONSTS(134)/CONSTS(135))*10.0000 ** CONSTS(182)*ALGBRC(47)
      ALGBRC(135) =  (( ALGBRC(133)*CONSTS(134))/CONSTS(135))*10.0000 ** CONSTS(232)
      ALGBRC(146) =  (( ALGBRC(144)*CONSTS(134))/CONSTS(135))*10.0000 ** CONSTS(237)
      ALGBRC(176) = ((ALGBRC(35)+ALGBRC(12)+ALGBRC(174)) - ALGBRC(173)) - ALGBRC(23)
      ALGBRC(178) = ( EXP(- ALGBRC(176)/( CONSTS(1)*CONSTS(131)))*ALGBRC(24)*ALGBRC(2))/( ALGBRC(13)*10.0000 ** - ALGBRC(1))
      ALGBRC(281) =  ALGBRC(185)*ALGBRC(239)+ ALGBRC(182)*(ALGBRC(227)+ALGBRC(235))+ ALGBRC(188)*ALGBRC(242)+ ALGBRC(191)*ALGBRC(249)+ ALGBRC(194)*ALGBRC(253)+ ALGBRC(197)*ALGBRC(255)+ ALGBRC(200)*ALGBRC(262)+ ALGBRC(206)*ALGBRC(266)+ ALGBRC(209)*ALGBRC(270)+ ALGBRC(212)*ALGBRC(273)+ ALGBRC(215)*ALGBRC(276)+ ALGBRC(218)*ALGBRC(279)+ ALGBRC(203)*ALGBRC(258)
      ALGBRC(289) =  ALGBRC(175)*- ALGBRC(287)
      ALGBRC(299) = 1.00000 -  ALGBRC(298)*ALGBRC(296)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END