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.00200000
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
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