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 251 entries in the algebraic variable array.
C There are a total of 57 entries in each of the rate and state variable arrays.
C There are a total of 232 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C CONSTS(1) is celltype in component environment (dimensionless).
C CONSTS(2) is nao in component extracellular (millimolar).
C CONSTS(3) is cao in component extracellular (millimolar).
C CONSTS(4) is ko in component extracellular (millimolar).
C CONSTS(5) is R in component physical_constants (joule_per_kilomole_kelvin).
C CONSTS(6) is T in component physical_constants (kelvin).
C CONSTS(7) is F in component physical_constants (coulomb_per_mole).
C CONSTS(8) is zna in component physical_constants (dimensionless).
C CONSTS(9) is zca in component physical_constants (dimensionless).
C CONSTS(10) is zk in component physical_constants (dimensionless).
C CONSTS(11) is L in component cell_geometry (centimeter).
C CONSTS(12) is rad in component cell_geometry (centimeter).
C CONSTS(185) is vcell in component cell_geometry (microliter).
C CONSTS(201) is Ageo in component cell_geometry (centimeter_squared).
C CONSTS(203) is Acap in component cell_geometry (centimeter_squared).
C CONSTS(204) is vmyo in component cell_geometry (microliter).
C CONSTS(205) is vnsr in component cell_geometry (microliter).
C CONSTS(206) is vjsr in component cell_geometry (microliter).
C CONSTS(207) is vss in component cell_geometry (microliter).
C CONSTS(208) is vsr in component cell_geometry (microliter).
C STATES(1) is v in component membrane (millivolt).
C ALGBRC(23) is vffrt in component membrane (coulomb_per_mole).
C ALGBRC(32) is vfrt in component membrane (dimensionless).
C ALGBRC(55) is INa in component INa (microA_per_microF).
C ALGBRC(57) is INaL in component INaL (microA_per_microF).
C ALGBRC(63) is Ito in component Ito (microA_per_microF).
C ALGBRC(144) is ICaL in component ICaL (microA_per_microF).
C ALGBRC(146) is ICaNa in component ICaL (microA_per_microF).
C ALGBRC(147) is ICaK in component ICaL (microA_per_microF).
C ALGBRC(148) is IKr in component IKr (microA_per_microF).
C ALGBRC(150) is IKs in component IKs (microA_per_microF).
C ALGBRC(152) is IK1 in component IK1 (microA_per_microF).
C ALGBRC(184) is INaCa_i in component INaCa_i (microA_per_microF).
C ALGBRC(214) is INaCa_ss in component INaCa_i (microA_per_microF).
C ALGBRC(233) is INaK in component INaK (microA_per_microF).
C ALGBRC(236) is INab in component INab (microA_per_microF).
C ALGBRC(235) is IKb in component IKb (microA_per_microF).
C ALGBRC(240) is IpCa in component IpCa (microA_per_microF).
C ALGBRC(238) is ICab in component ICab (microA_per_microF).
C CONSTS(13) is pstim in component membrane (dimensionless).
C ALGBRC(9) is Istim in component membrane (microA_per_microF).
C CONSTS(14) is i_Stim_Start in component membrane (millisecond).
C CONSTS(15) is i_Stim_End in component membrane (millisecond).
C CONSTS(16) is i_Stim_Amplitude in component membrane (microA_per_microF).
C CONSTS(17) is i_Stim_Period in component membrane (millisecond).
C CONSTS(18) is i_Stim_PulseDuration in component membrane (millisecond).
C CONSTS(19) is KmCaMK in component CaMK (millimolar).
C CONSTS(20) is aCaMK in component CaMK (per_millimolar_per_millisecond).
C CONSTS(21) is bCaMK in component CaMK (per_millisecond).
C CONSTS(22) is CaMKo in component CaMK (dimensionless).
C CONSTS(23) is KmCaM in component CaMK (millimolar).
C ALGBRC(38) is CaMKb in component CaMK (millimolar).
C ALGBRC(40) is CaMKa in component CaMK (millimolar).
C STATES(2) is CaMKt in component CaMK (millimolar).
C STATES(3) is cass in component intracellular_ions (millimolar).
C CONSTS(24) is cmdnmax_b in component intracellular_ions (millimolar).
C CONSTS(174) is cmdnmax in component intracellular_ions (millimolar).
C CONSTS(25) is kmcmdn in component intracellular_ions (millimolar).
C CONSTS(26) is trpnmax in component intracellular_ions (millimolar).
C CONSTS(27) is kmtrpn in component intracellular_ions (millimolar).
C CONSTS(28) is BSRmax in component intracellular_ions (millimolar).
C CONSTS(29) is KmBSR in component intracellular_ions (millimolar).
C CONSTS(30) is BSLmax in component intracellular_ions (millimolar).
C CONSTS(31) is KmBSL in component intracellular_ions (millimolar).
C CONSTS(32) is csqnmax in component intracellular_ions (millimolar).
C CONSTS(33) is kmcsqn in component intracellular_ions (millimolar).
C STATES(4) is nai in component intracellular_ions (millimolar).
C STATES(5) is nass in component intracellular_ions (millimolar).
C STATES(6) is ki in component intracellular_ions (millimolar).
C STATES(7) is kss in component intracellular_ions (millimolar).
C CONSTS(34) is cansr in component intracellular_ions (millimolar).
C CONSTS(35) is cajsr in component intracellular_ions (millimolar).
C STATES(8) is casr in component intracellular_ions (millimolar).
C STATES(9) is cai in component intracellular_ions (millimolar).
C ALGBRC(239) is JdiffNa in component diff (millimolar_per_millisecond).
C ALGBRC(241) is Jdiff in component diff (millimolar_per_millisecond).
C ALGBRC(251) is Jup in component SERCA (millimolar_per_millisecond).
C ALGBRC(237) is JdiffK in component diff (millimolar_per_millisecond).
C ALGBRC(246) is Jrel in component ryr (millimolar_per_millisecond).
C CONSTS(175) is Jtr in component trans_flux (millimolar_per_millisecond).
C ALGBRC(250) is Jleak in component SERCA (millimolar_per_millisecond).
C ALGBRC(42) is Bcai in component intracellular_ions (dimensionless).
C ALGBRC(46) is Bcasr in component intracellular_ions (dimensionless).
C ALGBRC(44) is Bcass in component intracellular_ions (dimensionless).
C CONSTS(36) is cm in component intracellular_ions (microF_per_centimeter_squared).
C CONSTS(37) is PKNa in component reversal_potentials (dimensionless).
C ALGBRC(49) is ENa in component reversal_potentials (millivolt).
C ALGBRC(50) is EK in component reversal_potentials (millivolt).
C ALGBRC(51) is EKs in component reversal_potentials (millivolt).
C CONSTS(38) is EKshift in component reversal_potentials (millivolt).
C CONSTS(39) is btj in component INa (dimensionless).
C CONSTS(40) is bGNa in component INa (dimensionless).
C ALGBRC(1) is mss in component INa (dimensionless).
C ALGBRC(16) is tm in component INa (millisecond).
C CONSTS(41) is mssV1 in component INa (millivolt).
C CONSTS(42) is mssV2 in component INa (millivolt).
C CONSTS(43) is mtV1 in component INa (millivolt).
C CONSTS(44) is mtV2 in component INa (millivolt).
C CONSTS(45) is mtD1 in component INa (dimensionless).
C CONSTS(46) is mtD2 in component INa (dimensionless).
C CONSTS(47) is mtV3 in component INa (millivolt).
C CONSTS(48) is mtV4 in component INa (millivolt).
C STATES(10) is m in component INa (dimensionless).
C ALGBRC(2) is hss in component INa (dimensionless).
C ALGBRC(17) is thf in component INa (millisecond).
C ALGBRC(18) is ths in component INa (millisecond).
C CONSTS(49) is hssV1 in component INa (millivolt).
C CONSTS(50) is hssV2 in component INa (millivolt).
C CONSTS(176) is Ahs in component INa (dimensionless).
C CONSTS(51) is Ahf in component INa (dimensionless).
C STATES(11) is hf in component INa (dimensionless).
C STATES(12) is hs in component INa (dimensionless).
C ALGBRC(52) is h in component INa (dimensionless).
C CONSTS(52) is GNa in component INa (milliS_per_microF).
C ALGBRC(19) is jss in component INa (dimensionless).
C ALGBRC(27) is tj in component INa (millisecond).
C STATES(13) is j in component INa (dimensionless).
C ALGBRC(28) is hssp in component INa (dimensionless).
C ALGBRC(34) is thsp in component INa (millisecond).
C STATES(14) is hsp in component INa (dimensionless).
C ALGBRC(53) is hp in component INa (dimensionless).
C ALGBRC(35) is tjp in component INa (millisecond).
C STATES(15) is jp in component INa (dimensionless).
C ALGBRC(54) is fINap in component INa (dimensionless).
C CONSTS(53) is bGnal in component INaL (dimensionless).
C CONSTS(54) is bthL in component INaL (dimensionless).
C ALGBRC(29) is mLss in component INaL (dimensionless).
C ALGBRC(36) is tmL in component INaL (millisecond).
C STATES(16) is mL in component INaL (dimensionless).
C CONSTS(177) is thL in component INaL (millisecond).
C ALGBRC(3) is hLss in component INaL (dimensionless).
C STATES(17) is hL in component INaL (dimensionless).
C ALGBRC(4) is hLssp in component INaL (dimensionless).
C CONSTS(199) is thLp in component INaL (millisecond).
C STATES(18) is hLp in component INaL (dimensionless).
C CONSTS(55) is GNaL_b in component INaL (milliS_per_microF).
C CONSTS(178) is GNaL in component INaL (milliS_per_microF).
C ALGBRC(56) is fINaLp in component INaL (dimensionless).
C CONSTS(56) is bGto in component Ito (dimensionless).
C CONSTS(57) is Gto_b in component Ito (milliS_per_microF).
C ALGBRC(5) is ass in component Ito (dimensionless).
C ALGBRC(20) is ta in component Ito (millisecond).
C STATES(19) is a in component Ito (dimensionless).
C ALGBRC(6) is iss in component Ito (dimensionless).
C ALGBRC(21) is delta_epi in component Ito (dimensionless).
C ALGBRC(30) is tiF_b in component Ito (millisecond).
C ALGBRC(37) is tiS_b in component Ito (millisecond).
C ALGBRC(39) is tiF in component Ito (millisecond).
C ALGBRC(41) is tiS in component Ito (millisecond).
C ALGBRC(58) is AiF in component Ito (dimensionless).
C ALGBRC(59) is AiS in component Ito (dimensionless).
C STATES(20) is iF in component Ito (dimensionless).
C STATES(21) is iS in component Ito (dimensionless).
C ALGBRC(60) is i in component Ito (dimensionless).
C ALGBRC(31) is assp in component Ito (dimensionless).
C STATES(22) is ap in component Ito (dimensionless).
C ALGBRC(43) is dti_develop in component Ito (dimensionless).
C ALGBRC(45) is dti_recover in component Ito (dimensionless).
C ALGBRC(47) is tiFp in component Ito (millisecond).
C ALGBRC(48) is tiSp in component Ito (millisecond).
C STATES(23) is iFp in component Ito (dimensionless).
C STATES(24) is iSp in component Ito (dimensionless).
C ALGBRC(61) is ip in component Ito (dimensionless).
C CONSTS(179) is Gto in component Ito (milliS_per_microF).
C ALGBRC(62) is fItop in component Ito (dimensionless).
C CONSTS(58) is EKshift in component Ito (millivolt).
C CONSTS(180) is r_down in component ICaL (dimensionless).
C ALGBRC(65) is r_up in component ICaL (dimensionless).
C CONSTS(59) is undo_CDI in component ICaL (dimensionless).
C STATES(25) is nca in component ICaL (dimensionless).
C ALGBRC(7) is jncass in component ICaL (dimensionless).
C CONSTS(60) is tjnca in component ICaL (millisecond).
C STATES(26) is jnca in component ICaL (dimensionless).
C ALGBRC(8) is km2n in component ICaL (per_millisecond).
C CONSTS(61) is Kmn in component ICaL (millimolar).
C ALGBRC(22) is anca in component ICaL (dimensionless).
C CONSTS(62) is k2n in component ICaL (per_millisecond).
C CONSTS(63) is kmn in component ICaL (millimolar).
C ALGBRC(66) is dss in component ICaL (dimensionless).
C ALGBRC(67) is td in component ICaL (millisecond).
C ALGBRC(68) is alpha in component ICaL (dimensionless).
C ALGBRC(69) is beta in component ICaL (dimensionless).
C ALGBRC(70) is jcass_new in component ICaL (dimensionless).
C ALGBRC(71) is jcass_VD in component ICaL (dimensionless).
C ALGBRC(72) is jcass_CD in component ICaL (dimensionless).
C ALGBRC(73) is jcass_VDp in component ICaL (dimensionless).
C ALGBRC(74) is jcass_CDp in component ICaL (dimensionless).
C ALGBRC(75) is tjca_new in component ICaL (dimensionless).
C ALGBRC(76) is tjca_VD in component ICaL (dimensionless).
C ALGBRC(77) is tjca_CD in component ICaL (dimensionless).
C ALGBRC(78) is tjca_VDp in component ICaL (dimensionless).
C ALGBRC(79) is tjca_CDp in component ICaL (dimensionless).
C ALGBRC(80) is psi_VD in component ICaL (dimensionless).
C ALGBRC(81) is psi_VDp in component ICaL (dimensionless).
C ALGBRC(82) is psi_CD in component ICaL (dimensionless).
C ALGBRC(83) is psi_CDp in component ICaL (dimensionless).
C ALGBRC(84) is omega_VD in component ICaL (dimensionless).
C ALGBRC(85) is omega_VDp in component ICaL (dimensionless).
C ALGBRC(86) is omega_CD in component ICaL (dimensionless).
C ALGBRC(87) is omega_CDp in component ICaL (dimensionless).
C ALGBRC(88) is f1ss_0 in component ICaL (dimensionless).
C ALGBRC(89) is tf1_0 in component ICaL (dimensionless).
C CONSTS(181) is ktaup in component ICaL (dimensionless).
C ALGBRC(90) is gamma_VD in component ICaL (dimensionless).
C ALGBRC(91) is delta_VD in component ICaL (dimensionless).
C ALGBRC(94) is gamma_VDp in component ICaL (dimensionless).
C ALGBRC(95) is delta_VDp in component ICaL (dimensionless).
C CONSTS(64) is kCDI in component ICaL (dimensionless).
C ALGBRC(96) is gamma_CD in component ICaL (dimensionless).
C ALGBRC(97) is delta_CD in component ICaL (dimensionless).
C ALGBRC(100) is gamma_CDp in component ICaL (dimensionless).
C ALGBRC(101) is delta_CDp in component ICaL (dimensionless).
C ALGBRC(92) is tf1_VD in component ICaL (dimensionless).
C ALGBRC(98) is tf1_CD in component ICaL (dimensionless).
C ALGBRC(93) is f1ss_VD in component ICaL (dimensionless).
C ALGBRC(99) is f1ss_CD in component ICaL (dimensionless).
C ALGBRC(102) is tf2_new in component ICaL (dimensionless).
C ALGBRC(104) is tf2_VD in component ICaL (dimensionless).
C ALGBRC(107) is tf2_CD in component ICaL (dimensionless).
C ALGBRC(103) is tf2_VDp in component ICaL (dimensionless).
C ALGBRC(110) is tf2_CDp in component ICaL (dimensionless).
C ALGBRC(105) is theta_VD in component ICaL (dimensionless).
C ALGBRC(113) is theta_CD in component ICaL (dimensionless).
C ALGBRC(106) is theta_VDp in component ICaL (dimensionless).
C ALGBRC(114) is theta_CDp in component ICaL (dimensionless).
C ALGBRC(108) is eta_VD in component ICaL (dimensionless).
C ALGBRC(109) is eta_VDp in component ICaL (dimensionless).
C ALGBRC(115) is eta_CD in component ICaL (dimensionless).
C ALGBRC(116) is eta_CDp in component ICaL (dimensionless).
C ALGBRC(111) is tf2post_VD in component ICaL (dimensionless).
C ALGBRC(117) is tf2post_CD in component ICaL (dimensionless).
C ALGBRC(112) is f2ss_VD in component ICaL (dimensionless).
C ALGBRC(118) is f2ss_CD in component ICaL (dimensionless).
C ALGBRC(119) is PhiCaL in component ICaL (coulomb_per_metre_3).
C ALGBRC(120) is PhiCaNa in component ICaL (coulomb_per_metre_3).
C ALGBRC(121) is PhiCaK in component ICaL (coulomb_per_metre_3).
C CONSTS(209) is PCa in component ICaL (dimensionless).
C CONSTS(210) is PCap in component ICaL (dimensionless).
C CONSTS(211) is PCaNa in component ICaL (dimensionless).
C CONSTS(212) is PCaK in component ICaL (dimensionless).
C CONSTS(213) is PCaNap in component ICaL (dimensionless).
C CONSTS(214) is PCaKp in component ICaL (dimensionless).
C CONSTS(65) is PCa_b in component ICaL (dimensionless).
C STATES(27) is I1k in component ICaL (dimensionless).
C STATES(28) is I2k in component ICaL (dimensionless).
C STATES(29) is Ck in component ICaL (dimensionless).
C STATES(30) is I1kp in component ICaL (dimensionless).
C STATES(31) is I2kp in component ICaL (dimensionless).
C STATES(32) is Ckp in component ICaL (dimensionless).
C STATES(33) is I1Cak in component ICaL (dimensionless).
C STATES(34) is I2Cak in component ICaL (dimensionless).
C STATES(35) is CCak in component ICaL (dimensionless).
C STATES(36) is I1Cakp in component ICaL (dimensionless).
C STATES(37) is I2Cakp in component ICaL (dimensionless).
C STATES(38) is CCakp in component ICaL (dimensionless).
C STATES(39) is Ok in component ICaL (dimensionless).
C STATES(40) is Okp in component ICaL (dimensionless).
C ALGBRC(122) is OCak in component ICaL (dimensionless).
C ALGBRC(123) is OCakp in component ICaL (dimensionless).
C ALGBRC(124) is ICaL_VD in component ICaL (dimensionless).
C ALGBRC(125) is ICaL_VDp in component ICaL (dimensionless).
C ALGBRC(127) is ICaL_CD in component ICaL (dimensionless).
C ALGBRC(128) is ICaL_CDp in component ICaL (dimensionless).
C ALGBRC(130) is ICaNa_VD in component ICaL (dimensionless).
C ALGBRC(131) is ICaNa_VDp in component ICaL (dimensionless).
C ALGBRC(132) is ICaNa_CD in component ICaL (dimensionless).
C ALGBRC(133) is ICaNa_CDp in component ICaL (dimensionless).
C ALGBRC(134) is ICaK_VD in component ICaL (dimensionless).
C ALGBRC(135) is ICaK_VDp in component ICaL (dimensionless).
C ALGBRC(136) is ICaK_CD in component ICaL (dimensionless).
C ALGBRC(137) is ICaK_CDp in component ICaL (dimensionless).
C ALGBRC(138) is ICaLnp in component ICaL (dimensionless).
C ALGBRC(139) is ICaLp in component ICaL (dimensionless).
C ALGBRC(126) is ICaLVD in component ICaL (dimensionless).
C ALGBRC(129) is ICaLCD in component ICaL (dimensionless).
C ALGBRC(140) is ICaNanp in component ICaL (dimensionless).
C ALGBRC(141) is ICaNap in component ICaL (dimensionless).
C ALGBRC(142) is ICaKnp in component ICaL (dimensionless).
C ALGBRC(143) is ICaKp in component ICaL (dimensionless).
C ALGBRC(64) is fICaLp in component ICaL (dimensionless).
C CONSTS(66) is bGCaL in component ICaL (dimensionless).
C ALGBRC(145) is gICaL in component ICaL (microF).
C CONSTS(67) is GKr_b in component IKr (milliS_per_microF).
C STATES(41) is IC1 in component IKr (dimensionless).
C STATES(42) is IC2 in component IKr (dimensionless).
C STATES(43) is C1 in component IKr (dimensionless).
C STATES(44) is C2 in component IKr (dimensionless).
C STATES(45) is O in component IKr (dimensionless).
C STATES(46) is IO in component IKr (dimensionless).
C STATES(47) is IObound in component IKr (dimensionless).
C STATES(48) is Obound in component IKr (dimensionless).
C STATES(49) is Cbound in component IKr (dimensionless).
C STATES(50) is D in component IKr (dimensionless).
C CONSTS(182) is GKr in component IKr (milliS_per_microF).
C CONSTS(200) is GKr_total in component IKr (milliS_per_microF).
C CONSTS(68) is A1 in component IKr (per_millisecond).
C CONSTS(69) is B1 in component IKr (per_millivolt).
C CONSTS(70) is q1 in component IKr (dimensionless).
C CONSTS(71) is A2 in component IKr (per_millisecond).
C CONSTS(72) is B2 in component IKr (per_millivolt).
C CONSTS(73) is q2 in component IKr (dimensionless).
C CONSTS(74) is A3 in component IKr (per_millisecond).
C CONSTS(75) is B3 in component IKr (per_millivolt).
C CONSTS(76) is q3 in component IKr (dimensionless).
C CONSTS(77) is A4 in component IKr (per_millisecond).
C CONSTS(78) is B4 in component IKr (per_millivolt).
C CONSTS(79) is q4 in component IKr (dimensionless).
C CONSTS(80) is A11 in component IKr (per_millisecond).
C CONSTS(81) is B11 in component IKr (per_millivolt).
C CONSTS(82) is q11 in component IKr (dimensionless).
C CONSTS(83) is A21 in component IKr (per_millisecond).
C CONSTS(84) is B21 in component IKr (per_millivolt).
C CONSTS(85) is q21 in component IKr (dimensionless).
C CONSTS(86) is A31 in component IKr (per_millisecond).
C CONSTS(87) is B31 in component IKr (per_millivolt).
C CONSTS(88) is q31 in component IKr (dimensionless).
C CONSTS(89) is A41 in component IKr (per_millisecond).
C CONSTS(90) is B41 in component IKr (per_millivolt).
C CONSTS(91) is q41 in component IKr (dimensionless).
C CONSTS(92) is A51 in component IKr (per_millisecond).
C CONSTS(93) is B51 in component IKr (per_millivolt).
C CONSTS(94) is q51 in component IKr (dimensionless).
C CONSTS(95) is A52 in component IKr (per_millisecond).
C CONSTS(96) is B52 in component IKr (per_millivolt).
C CONSTS(97) is q52 in component IKr (dimensionless).
C CONSTS(98) is A53 in component IKr (per_millisecond).
C CONSTS(99) is B53 in component IKr (per_millivolt).
C CONSTS(100) is q53 in component IKr (dimensionless).
C CONSTS(101) is A61 in component IKr (per_millisecond).
C CONSTS(102) is B61 in component IKr (per_millivolt).
C CONSTS(103) is q61 in component IKr (dimensionless).
C CONSTS(104) is A62 in component IKr (per_millisecond).
C CONSTS(105) is B62 in component IKr (per_millivolt).
C CONSTS(106) is q62 in component IKr (dimensionless).
C CONSTS(107) is A63 in component IKr (per_millisecond).
C CONSTS(108) is B63 in component IKr (per_millivolt).
C CONSTS(109) is q63 in component IKr (dimensionless).
C CONSTS(110) is Kmax in component IKr (dimensionless).
C CONSTS(111) is Ku in component IKr (per_millisecond).
C CONSTS(112) is n in component IKr (dimensionless).
C CONSTS(113) is halfmax in component IKr (dimensionless).
C CONSTS(114) is Kt in component IKr (per_millisecond).
C CONSTS(115) is Vhalf in component IKr (millivolt).
C CONSTS(116) is Temp in component IKr (kelvin).
C CONSTS(117) is bGKr in component IKr (dimensionless).
C CONSTS(118) is GKs_b in component IKs (milliS_per_microF).
C CONSTS(119) is bGKs in component IKs (milliS_per_microF).
C CONSTS(120) is EKshift in component IKs (millivolt).
C CONSTS(183) is GKs in component IKs (milliS_per_microF).
C ALGBRC(10) is xs1ss in component IKs (dimensionless).
C ALGBRC(24) is xs2ss in component IKs (dimensionless).
C ALGBRC(25) is txs1 in component IKs (millisecond).
C CONSTS(121) is txs1_max in component IKs (millisecond).
C STATES(51) is xs1 in component IKs (dimensionless).
C STATES(52) is xs2 in component IKs (dimensionless).
C ALGBRC(149) is KsCa in component IKs (dimensionless).
C ALGBRC(33) is txs2 in component IKs (millisecond).
C CONSTS(184) is GK1 in component IK1 (milliS_per_microF).
C CONSTS(122) is GK1_b in component IK1 (milliS_per_microF).
C CONSTS(123) is bGK1 in component IK1 (dimensionless).
C CONSTS(124) is EKshift in component IK1 (millivolt).
C ALGBRC(11) is xk1ss in component IK1 (dimensionless).
C ALGBRC(26) is txk1 in component IK1 (millisecond).
C STATES(53) is xk1 in component IK1 (dimensionless).
C ALGBRC(151) is rk1 in component IK1 (millisecond).
C CONSTS(125) is kslope_rk1 in component IK1 (dimensionless).
C CONSTS(126) is kna1 in component INaCa_i (per_millisecond).
C CONSTS(127) is kna2 in component INaCa_i (per_millisecond).
C CONSTS(128) is kna3 in component INaCa_i (per_millisecond).
C CONSTS(129) is kasymm in component INaCa_i (dimensionless).
C CONSTS(130) is wna in component INaCa_i (dimensionless).
C CONSTS(131) is wca in component INaCa_i (dimensionless).
C CONSTS(132) is wnaca in component INaCa_i (dimensionless).
C CONSTS(133) is kcaon in component INaCa_i (per_millisecond).
C CONSTS(134) is kcaoff in component INaCa_i (per_millisecond).
C CONSTS(135) is qna in component INaCa_i (dimensionless).
C CONSTS(136) is qca in component INaCa_i (dimensionless).
C ALGBRC(154) is hna in component INaCa_i (dimensionless).
C ALGBRC(153) is hca in component INaCa_i (dimensionless).
C CONSTS(137) is KmCaAct in component INaCa_i (millimolar).
C CONSTS(138) is Gncx_b in component INaCa_i (milliS_per_microF).
C CONSTS(139) is bGncx in component INaCa_i (dimensionless).
C CONSTS(221) is Gncx in component INaCa_i (milliS_per_microF).
C ALGBRC(155) is h1_i in component INaCa_i (dimensionless).
C ALGBRC(156) is h2_i in component INaCa_i (dimensionless).
C ALGBRC(157) is h3_i in component INaCa_i (dimensionless).
C ALGBRC(158) is h4_i in component INaCa_i (dimensionless).
C ALGBRC(159) is h5_i in component INaCa_i (dimensionless).
C ALGBRC(160) is h6_i in component INaCa_i (dimensionless).
C ALGBRC(161) is h7_i in component INaCa_i (dimensionless).
C ALGBRC(162) is h8_i in component INaCa_i (dimensionless).
C ALGBRC(163) is h9_i in component INaCa_i (dimensionless).
C CONSTS(215) is h10_i in component INaCa_i (dimensionless).
C CONSTS(216) is h11_i in component INaCa_i (dimensionless).
C CONSTS(217) is h12_i in component INaCa_i (dimensionless).
C CONSTS(218) is k1_i in component INaCa_i (dimensionless).
C CONSTS(219) is k2_i in component INaCa_i (dimensionless).
C ALGBRC(164) is k3p_i in component INaCa_i (dimensionless).
C ALGBRC(165) is k3pp_i in component INaCa_i (dimensionless).
C ALGBRC(166) is k3_i in component INaCa_i (dimensionless).
C ALGBRC(169) is k4_i in component INaCa_i (dimensionless).
C ALGBRC(167) is k4p_i in component INaCa_i (dimensionless).
C ALGBRC(168) is k4pp_i in component INaCa_i (dimensionless).
C CONSTS(220) is k5_i in component INaCa_i (dimensionless).
C ALGBRC(170) is k6_i in component INaCa_i (dimensionless).
C ALGBRC(171) is k7_i in component INaCa_i (dimensionless).
C ALGBRC(172) is k8_i in component INaCa_i (dimensionless).
C ALGBRC(173) is x1_i in component INaCa_i (dimensionless).
C ALGBRC(174) is x2_i in component INaCa_i (dimensionless).
C ALGBRC(175) is x3_i in component INaCa_i (dimensionless).
C ALGBRC(176) is x4_i in component INaCa_i (dimensionless).
C ALGBRC(177) is E1_i in component INaCa_i (dimensionless).
C ALGBRC(178) is E2_i in component INaCa_i (dimensionless).
C ALGBRC(179) is E3_i in component INaCa_i (dimensionless).
C ALGBRC(180) is E4_i in component INaCa_i (dimensionless).
C ALGBRC(181) is allo_i in component INaCa_i (dimensionless).
C ALGBRC(182) is JncxNa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(183) is JncxCa_i in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(185) is h1_ss in component INaCa_i (dimensionless).
C ALGBRC(186) is h2_ss in component INaCa_i (dimensionless).
C ALGBRC(187) is h3_ss in component INaCa_i (dimensionless).
C ALGBRC(188) is h4_ss in component INaCa_i (dimensionless).
C ALGBRC(189) is h5_ss in component INaCa_i (dimensionless).
C ALGBRC(190) is h6_ss in component INaCa_i (dimensionless).
C ALGBRC(191) is h7_ss in component INaCa_i (dimensionless).
C ALGBRC(192) is h8_ss in component INaCa_i (dimensionless).
C ALGBRC(193) is h9_ss in component INaCa_i (dimensionless).
C CONSTS(222) is h10_ss in component INaCa_i (dimensionless).
C CONSTS(223) is h11_ss in component INaCa_i (dimensionless).
C CONSTS(224) is h12_ss in component INaCa_i (dimensionless).
C CONSTS(225) is k1_ss in component INaCa_i (dimensionless).
C CONSTS(226) is k2_ss in component INaCa_i (dimensionless).
C ALGBRC(194) is k3p_ss in component INaCa_i (dimensionless).
C ALGBRC(195) is k3pp_ss in component INaCa_i (dimensionless).
C ALGBRC(196) is k3_ss in component INaCa_i (dimensionless).
C ALGBRC(199) is k4_ss in component INaCa_i (dimensionless).
C ALGBRC(197) is k4p_ss in component INaCa_i (dimensionless).
C ALGBRC(198) is k4pp_ss in component INaCa_i (dimensionless).
C CONSTS(227) is k5_ss in component INaCa_i (dimensionless).
C ALGBRC(200) is k6_ss in component INaCa_i (dimensionless).
C ALGBRC(201) is k7_ss in component INaCa_i (dimensionless).
C ALGBRC(202) is k8_ss in component INaCa_i (dimensionless).
C ALGBRC(203) is x1_ss in component INaCa_i (dimensionless).
C ALGBRC(204) is x2_ss in component INaCa_i (dimensionless).
C ALGBRC(205) is x3_ss in component INaCa_i (dimensionless).
C ALGBRC(206) is x4_ss in component INaCa_i (dimensionless).
C ALGBRC(207) is E1_ss in component INaCa_i (dimensionless).
C ALGBRC(208) is E2_ss in component INaCa_i (dimensionless).
C ALGBRC(209) is E3_ss in component INaCa_i (dimensionless).
C ALGBRC(210) is E4_ss in component INaCa_i (dimensionless).
C ALGBRC(211) is allo_ss in component INaCa_i (dimensionless).
C ALGBRC(212) is JncxNa_ss in component INaCa_i (millimolar_per_millisecond).
C ALGBRC(213) is JncxCa_ss in component INaCa_i (millimolar_per_millisecond).
C CONSTS(140) is k1p in component INaK (per_millisecond).
C CONSTS(141) is k1m in component INaK (per_millisecond).
C CONSTS(142) is k2p in component INaK (per_millisecond).
C CONSTS(143) is k2m in component INaK (per_millisecond).
C CONSTS(144) is k3p in component INaK (per_millisecond).
C CONSTS(145) is k3m in component INaK (per_millisecond).
C CONSTS(146) is k4p in component INaK (per_millisecond).
C CONSTS(147) is k4m in component INaK (per_millisecond).
C CONSTS(148) is Knai0 in component INaK (millimolar).
C CONSTS(149) is Knao0 in component INaK (millimolar).
C CONSTS(150) is delta in component INaK (millivolt).
C CONSTS(151) is Kki in component INaK (per_millisecond).
C CONSTS(152) is Kko in component INaK (per_millisecond).
C CONSTS(153) is MgADP in component INaK (millimolar).
C CONSTS(154) is MgATP in component INaK (millimolar).
C CONSTS(155) is Kmgatp in component INaK (millimolar).
C CONSTS(156) is H in component INaK (millimolar).
C CONSTS(157) is eP in component INaK (dimensionless).
C CONSTS(158) is Khp in component INaK (millimolar).
C CONSTS(159) is Knap in component INaK (millimolar).
C CONSTS(160) is Kxkur in component INaK (millimolar).
C CONSTS(161) is Pnak_b in component INaK (milliS_per_microF).
C CONSTS(231) is Pnak in component INaK (milliS_per_microF).
C CONSTS(162) is bGnak in component INaK (dimensionless).
C ALGBRC(215) is Knai in component INaK (millimolar).
C ALGBRC(216) is Knao in component INaK (millimolar).
C ALGBRC(217) is P in component INaK (dimensionless).
C ALGBRC(218) is a1 in component INaK (dimensionless).
C CONSTS(228) is b1 in component INaK (dimensionless).
C CONSTS(229) is a2 in component INaK (dimensionless).
C ALGBRC(219) is b2 in component INaK (dimensionless).
C ALGBRC(220) is a3 in component INaK (dimensionless).
C ALGBRC(221) is b3 in component INaK (dimensionless).
C CONSTS(230) is a4 in component INaK (dimensionless).
C ALGBRC(222) is b4 in component INaK (dimensionless).
C ALGBRC(223) is x1 in component INaK (dimensionless).
C ALGBRC(224) is x2 in component INaK (dimensionless).
C ALGBRC(225) is x3 in component INaK (dimensionless).
C ALGBRC(226) is x4 in component INaK (dimensionless).
C ALGBRC(227) is E1 in component INaK (dimensionless).
C ALGBRC(228) is E2 in component INaK (dimensionless).
C ALGBRC(229) is E3 in component INaK (dimensionless).
C ALGBRC(230) is E4 in component INaK (dimensionless).
C ALGBRC(231) is JnakNa in component INaK (millimolar_per_millisecond).
C ALGBRC(232) is JnakK in component INaK (millimolar_per_millisecond).
C ALGBRC(234) is xkb in component IKb (dimensionless).
C CONSTS(163) is GKb_b in component IKb (milliS_per_microF).
C CONSTS(186) is GKb in component IKb (milliS_per_microF).
C CONSTS(164) is PNab in component INab (milliS_per_microF).
C CONSTS(165) is PCab in component ICab (milliS_per_microF).
C CONSTS(166) is undo_ICab in component ICab (dimensionless).
C CONSTS(167) is GpCa in component IpCa (milliS_per_microF).
C CONSTS(168) is KmCap in component IpCa (millimolar).
C CONSTS(169) is RyRa1 in component ryr (micromolar).
C CONSTS(170) is RyRa2 in component ryr (micromolar).
C CONSTS(187) is RyRohalf in component ryr (micromolar).
C CONSTS(188) is RyRchalf in component ryr (micromolar).
C ALGBRC(242) is RyRSRCass in component ryr (dimensionless).
C ALGBRC(12) is RyRainfss in component ryr (micromolar).
C CONSTS(189) is RyRtauadapt in component ryr (millisecond).
C STATES(54) is RyRa in component ryr (micromolar).
C ALGBRC(13) is RyRoinfss in component ryr (dimensionless).
C CONSTS(190) is RyRtauact in component ryr (millisecond).
C ALGBRC(14) is RyRcinfss in component ryr (dimensionless).
C STATES(55) is RyRo in component ryr (dimensionless).
C CONSTS(191) is RyRtauinact in component ryr (millisecond).
C STATES(56) is RyRc in component ryr (dimensionless).
C CONSTS(202) is RyRtauinactp in component ryr (millisecond).
C STATES(57) is RyRcp in component ryr (dimensionless).
C ALGBRC(243) is fJrelp in component ryr (dimensionless).
C ALGBRC(244) is Jrelnp in component ryr (millimolar_per_millisecond).
C ALGBRC(245) is Jrelp in component ryr (millimolar_per_millisecond).
C CONSTS(171) is g_irel_max in component ryr (millimolar_per_second).
C CONSTS(192) is g_irel_max_M in component ryr (per_second).
C CONSTS(193) is g_irel_max_p in component ryr (per_second).
C CONSTS(194) is upScale in component SERCA (dimensionless).
C ALGBRC(247) is Jupnp in component SERCA (millimolar_per_millisecond).
C ALGBRC(248) is Jupp in component SERCA (millimolar_per_millisecond).
C ALGBRC(249) is fJupp in component SERCA (dimensionless).
C CONSTS(172) is Jup_b in component SERCA (dimensionless).
C CONSTS(173) is cJup in component SERCA (dimensionless).
C CONSTS(195) is Vmax_SRCaP in component SERCA (millimolar_per_millisecond).
C CONSTS(196) is Kmf in component SERCA (millimolar).
C CONSTS(197) is Kmr in component SERCA (millimolar).
C CONSTS(198) is hillSRCaP in component SERCA (millimolar).
C ALGBRC(15) is Jup2 in component SERCA (millimolar_per_millisecond).
C RATES(1) is d/dt v in component membrane (millivolt).
C RATES(2) is d/dt CaMKt in component CaMK (millimolar).
C RATES(4) is d/dt nai in component intracellular_ions (millimolar).
C RATES(5) is d/dt nass in component intracellular_ions (millimolar).
C RATES(6) is d/dt ki in component intracellular_ions (millimolar).
C RATES(7) is d/dt kss in component intracellular_ions (millimolar).
C RATES(9) is d/dt cai in component intracellular_ions (millimolar).
C RATES(3) is d/dt cass in component intracellular_ions (millimolar).
C RATES(8) is d/dt casr in component intracellular_ions (millimolar).
C RATES(10) is d/dt m in component INa (dimensionless).
C RATES(11) is d/dt hf in component INa (dimensionless).
C RATES(12) is d/dt hs in component INa (dimensionless).
C RATES(13) is d/dt j in component INa (dimensionless).
C RATES(14) is d/dt hsp in component INa (dimensionless).
C RATES(15) is d/dt jp in component INa (dimensionless).
C RATES(16) is d/dt mL in component INaL (dimensionless).
C RATES(17) is d/dt hL in component INaL (dimensionless).
C RATES(18) is d/dt hLp in component INaL (dimensionless).
C RATES(19) is d/dt a in component Ito (dimensionless).
C RATES(20) is d/dt iF in component Ito (dimensionless).
C RATES(21) is d/dt iS in component Ito (dimensionless).
C RATES(22) is d/dt ap in component Ito (dimensionless).
C RATES(23) is d/dt iFp in component Ito (dimensionless).
C RATES(24) is d/dt iSp in component Ito (dimensionless).
C RATES(26) is d/dt jnca in component ICaL (dimensionless).
C RATES(25) is d/dt nca in component ICaL (dimensionless).
C RATES(39) is d/dt Ok in component ICaL (dimensionless).
C RATES(28) is d/dt I2k in component ICaL (dimensionless).
C RATES(27) is d/dt I1k in component ICaL (dimensionless).
C RATES(29) is d/dt Ck in component ICaL (dimensionless).
C RATES(40) is d/dt Okp in component ICaL (dimensionless).
C RATES(31) is d/dt I2kp in component ICaL (dimensionless).
C RATES(30) is d/dt I1kp in component ICaL (dimensionless).
C RATES(32) is d/dt Ckp in component ICaL (dimensionless).
C RATES(34) is d/dt I2Cak in component ICaL (dimensionless).
C RATES(33) is d/dt I1Cak in component ICaL (dimensionless).
C RATES(35) is d/dt CCak in component ICaL (dimensionless).
C RATES(37) is d/dt I2Cakp in component ICaL (dimensionless).
C RATES(36) is d/dt I1Cakp in component ICaL (dimensionless).
C RATES(38) is d/dt CCakp in component ICaL (dimensionless).
C RATES(41) is d/dt IC1 in component IKr (dimensionless).
C RATES(42) is d/dt IC2 in component IKr (dimensionless).
C RATES(43) is d/dt C1 in component IKr (dimensionless).
C RATES(44) is d/dt C2 in component IKr (dimensionless).
C RATES(45) is d/dt O in component IKr (dimensionless).
C RATES(46) is d/dt IO in component IKr (dimensionless).
C RATES(47) is d/dt IObound in component IKr (dimensionless).
C RATES(48) is d/dt Obound in component IKr (dimensionless).
C RATES(49) is d/dt Cbound in component IKr (dimensionless).
C RATES(50) is d/dt D in component IKr (dimensionless).
C RATES(51) is d/dt xs1 in component IKs (dimensionless).
C RATES(52) is d/dt xs2 in component IKs (dimensionless).
C RATES(53) is d/dt xk1 in component IK1 (dimensionless).
C RATES(54) is d/dt RyRa in component ryr (micromolar).
C RATES(55) is d/dt RyRo in component ryr (dimensionless).
C RATES(56) is d/dt RyRc in component ryr (dimensionless).
C RATES(57) is d/dt RyRcp in component ryr (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0
CONSTS(2) = 144
CONSTS(3) = 2.7
CONSTS(4) = 5.4
CONSTS(5) = 8314.0
CONSTS(6) = 310.0
CONSTS(7) = 96485.0
CONSTS(8) = 1.0
CONSTS(9) = 2.0
CONSTS(10) = 1.0
CONSTS(11) = 0.01
CONSTS(12) = 0.0011
STATES(1) = -87
CONSTS(13) = 3
CONSTS(14) = 0
CONSTS(15) = 999000
CONSTS(16) = -53
CONSTS(17) = 1000
CONSTS(18) = 1
CONSTS(19) = 0.15
CONSTS(20) = 0.05
CONSTS(21) = 0.00068
CONSTS(22) = 0.05
CONSTS(23) = 0.0015
STATES(2) = 0
STATES(3) = 1.0e-04
CONSTS(24) = 0.05
CONSTS(25) = 0.00238
CONSTS(26) = 0.07
CONSTS(27) = 0.0005
CONSTS(28) = 0.047
CONSTS(29) = 0.00087
CONSTS(30) = 1.124
CONSTS(31) = 0.0087
CONSTS(32) = 1.0
CONSTS(33) = 0.8
STATES(4) = 7
STATES(5) = 7
STATES(6) = 145
STATES(7) = 145
CONSTS(34) = 1.50490908825974
CONSTS(35) = 0.821364973867864
STATES(8) = 1.2
STATES(9) = 1.0e-04
CONSTS(36) = 1
CONSTS(37) = 0.01833
CONSTS(38) = 8
CONSTS(39) = 1
CONSTS(40) = 0.27
CONSTS(41) = 39.57
CONSTS(42) = 9.871
CONSTS(43) = 11.64
CONSTS(44) = 34.77
CONSTS(45) = 6.765
CONSTS(46) = 8.552
CONSTS(47) = 77.42
CONSTS(48) = 5.955
STATES(10) = 0
CONSTS(49) = 78.5
CONSTS(50) = 6.22
CONSTS(51) = 0.99
STATES(11) = 1
STATES(12) = 1
CONSTS(52) = 75
STATES(13) = 1
STATES(14) = 1
STATES(15) = 1
CONSTS(53) = 2.8
CONSTS(54) = 1
STATES(16) = 0
STATES(17) = 1
STATES(18) = 1
CONSTS(55) = 0.0075
CONSTS(56) = 1
CONSTS(57) = 0.02
STATES(19) = 0
STATES(20) = 1
STATES(21) = 1
STATES(22) = 0
STATES(23) = 1
STATES(24) = 1
CONSTS(58) = 8
CONSTS(59) = 0
STATES(25) = 0
CONSTS(60) = 1
STATES(26) = 1
CONSTS(61) = 0.05
CONSTS(62) = 1000
CONSTS(63) = 0.05
CONSTS(64) = 9
CONSTS(65) = 0.0001
STATES(27) = 0
STATES(28) = 0
STATES(29) = 1
STATES(30) = 0
STATES(31) = 0
STATES(32) = 1
STATES(33) = 0
STATES(34) = 0
STATES(35) = 0
STATES(36) = 0
STATES(37) = 0
STATES(38) = 0
STATES(39) = 0
STATES(40) = 0
CONSTS(66) = 1
CONSTS(67) = 0.046
STATES(41) = 0.999637
STATES(42) = 6.83208e-05
STATES(43) = 1.80145e-08
STATES(44) = 8.26619e-05
STATES(45) = 0.00015551
STATES(46) = 5.67623e-05
STATES(47) = 0
STATES(48) = 0
STATES(49) = 0
STATES(50) = 0
CONSTS(68) = 0.0264
CONSTS(69) = 4.631e-05
CONSTS(70) = 4.843
CONSTS(71) = 4.986e-06
CONSTS(72) = -0.004226
CONSTS(73) = 4.23
CONSTS(74) = 0.001214
CONSTS(75) = 0.008516
CONSTS(76) = 4.962
CONSTS(77) = 1.854e-05
CONSTS(78) = -0.04641
CONSTS(79) = 3.769
CONSTS(80) = 0.0007868
CONSTS(81) = 1.535e-08
CONSTS(82) = 4.942
CONSTS(83) = 5.455e-06
CONSTS(84) = -0.1688
CONSTS(85) = 4.156
CONSTS(86) = 0.005509
CONSTS(87) = 7.771e-09
CONSTS(88) = 4.22
CONSTS(89) = 0.001416
CONSTS(90) = -0.02877
CONSTS(91) = 1.459
CONSTS(92) = 0.4492
CONSTS(93) = 0.008595
CONSTS(94) = 5.0
CONSTS(95) = 0.3181
CONSTS(96) = 3.613e-08
CONSTS(97) = 4.663
CONSTS(98) = 0.149
CONSTS(99) = 0.004668
CONSTS(100) = 2.412
CONSTS(101) = 0.01241
CONSTS(102) = 0.1725
CONSTS(103) = 5.568
CONSTS(104) = 0.3226
CONSTS(105) = -0.0006575
CONSTS(106) = 5.0
CONSTS(107) = 0.008978
CONSTS(108) = -0.02215
CONSTS(109) = 5.682
CONSTS(110) = 0
CONSTS(111) = 0
CONSTS(112) = 1.0
CONSTS(113) = 1
CONSTS(114) = 3.5e-5
CONSTS(115) = 1.0
CONSTS(116) = 310.0
CONSTS(117) = 1.2
CONSTS(118) = 0.0034
CONSTS(119) = 2
CONSTS(120) = 8
CONSTS(121) = 817.3
STATES(51) = 0
STATES(52) = 0
CONSTS(122) = 0.1908
CONSTS(123) = 0.71
CONSTS(124) = 8
STATES(53) = 1
CONSTS(125) = 1.09
CONSTS(126) = 15.0
CONSTS(127) = 5.0
CONSTS(128) = 88.12
CONSTS(129) = 12.5
CONSTS(130) = 60000
CONSTS(131) = 60000
CONSTS(132) = 5000
CONSTS(133) = 1500000
CONSTS(134) = 5000
CONSTS(135) = 0.5224
CONSTS(136) = 0.1670
CONSTS(137) = 150e-06
CONSTS(138) = 0.0008
CONSTS(139) = 2.4
CONSTS(140) = 949.5
CONSTS(141) = 182.4
CONSTS(142) = 687.2
CONSTS(143) = 39.4
CONSTS(144) = 1899.0
CONSTS(145) = 79300.0
CONSTS(146) = 639.0
CONSTS(147) = 40.0
CONSTS(148) = 9.073
CONSTS(149) = 27.78
CONSTS(150) = -0.1550
CONSTS(151) = 0.5
CONSTS(152) = 0.3582
CONSTS(153) = 0.05
CONSTS(154) = 9.8
CONSTS(155) = 1.698e-7
CONSTS(156) = 1e-7
CONSTS(157) = 4.2
CONSTS(158) = 1.698e-7
CONSTS(159) = 224.0
CONSTS(160) = 292.0
CONSTS(161) = 30
CONSTS(162) = 2
CONSTS(163) = 0.003
CONSTS(164) = 3.75e-10
CONSTS(165) = 2.5e-8
CONSTS(166) = 0
CONSTS(167) = 0.0005
CONSTS(168) = 0.0005
CONSTS(169) = 0.05
CONSTS(170) = 0.03
STATES(54) = 0.03
STATES(55) = 0
STATES(56) = 1
STATES(57) = 1
CONSTS(171) = 0.02
CONSTS(172) = 1
CONSTS(173) = 3.13
CONSTS(174) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(24)*1.20000, CONSTS(24))
CONSTS(175) = 0.00000
CONSTS(176) = 1.00000 - CONSTS(51)
CONSTS(177) = 200.000*CONSTS(54)
CONSTS(178) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(55)*CONSTS(53)*0.700000, CONSTS(55)*CONSTS(53))
CONSTS(179) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(56)*CONSTS(57)*4.00000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(56)*CONSTS(57)*4.00000, CONSTS(57)*CONSTS(56))
CONSTS(180) = 0.100000*(1.00000 - CONSTS(59))
CONSTS(181) = 2.50000
CONSTS(182) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(67)*CONSTS(117)*1.10000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(67)*CONSTS(117)*0.800000, CONSTS(67)*CONSTS(117))
CONSTS(183) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(119)*CONSTS(118)*1.40000, CONSTS(119)*CONSTS(118))
CONSTS(184) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(123)*CONSTS(122)*1.20000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(123)*CONSTS(122)*1.30000, CONSTS(123)*CONSTS(122))
CONSTS(185) = 1000.00*3.14160*CONSTS(12)*CONSTS(12)*CONSTS(11)
CONSTS(186) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(163)*0.600000, CONSTS(163))
CONSTS(187) = 0.120000 - (CONSTS(169) - CONSTS(170)/2.00000)
CONSTS(188) = 0.100000 - (CONSTS(169) - CONSTS(170)/2.00000)
CONSTS(189) = 1000.00
CONSTS(190) = 1.87500/1.87500
CONSTS(191) = ( 2.00000*87.5000)/10.0000
CONSTS(192) = CONSTS(171)*1.70000
CONSTS(193) = CONSTS(171)*1.25000
CONSTS(194) = TERNRY(CONSTS(1).EQ.1.00000, 1.30000, 1.00000)
CONSTS(195) = 1.00000*0.00531140
CONSTS(196) = 0.000246000
CONSTS(197) = 1.70000
CONSTS(198) = 1.78700
CONSTS(231) = 0.00000
CONSTS(199) = 3.00000*CONSTS(177)
CONSTS(200) = CONSTS(182)* (CONSTS(4)/5.40000) ** (1.0 / 2)
CONSTS(201) = 2.00000*3.14160*CONSTS(12)*CONSTS(12)+ 2.00000*3.14160*CONSTS(12)*CONSTS(11)
CONSTS(202) = CONSTS(191)*1.25000
CONSTS(203) = 2.00000*CONSTS(201)
CONSTS(204) = 0.680000*CONSTS(185)
CONSTS(205) = 0.0552000*CONSTS(185)
CONSTS(206) = 0.00480000*CONSTS(185)
CONSTS(207) = 0.0200000*CONSTS(185)
CONSTS(208) = 0.950000*(CONSTS(205)+CONSTS(206))
CONSTS(209) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(65)*0.900000*1.40000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(65)*0.900000*2.00000, CONSTS(65)*0.900000)
CONSTS(210) = 1.10000*CONSTS(209)
CONSTS(211) = 0.00125000*CONSTS(209)
CONSTS(212) = 0.000357400*CONSTS(209)
CONSTS(213) = 0.00125000*CONSTS(210)
CONSTS(214) = 0.000357400*CONSTS(210)
CONSTS(215) = CONSTS(129)+1.00000+ (CONSTS(2)/CONSTS(126))*(1.00000+CONSTS(2)/CONSTS(127))
CONSTS(216) = ( CONSTS(2)*CONSTS(2))/( CONSTS(215)*CONSTS(126)*CONSTS(127))
CONSTS(217) = 1.00000/CONSTS(215)
CONSTS(218) = CONSTS(217)*CONSTS(3)*CONSTS(133)
CONSTS(219) = CONSTS(134)
CONSTS(220) = CONSTS(134)
CONSTS(221) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(139)*CONSTS(138)*1.20000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(139)*CONSTS(138)*1.40000, CONSTS(139)*CONSTS(138))
CONSTS(222) = CONSTS(129)+1.00000+ (CONSTS(2)/CONSTS(126))*(1.00000+CONSTS(2)/CONSTS(127))
CONSTS(223) = ( CONSTS(2)*CONSTS(2))/( CONSTS(222)*CONSTS(126)*CONSTS(127))
CONSTS(224) = 1.00000/CONSTS(222)
CONSTS(225) = CONSTS(224)*CONSTS(3)*CONSTS(133)
CONSTS(226) = CONSTS(134)
CONSTS(227) = CONSTS(134)
CONSTS(228) = CONSTS(141)*CONSTS(153)
CONSTS(229) = CONSTS(142)
CONSTS(230) = (( CONSTS(146)*CONSTS(154))/CONSTS(155))/(1.00000+CONSTS(154)/CONSTS(155))
CONSTS(231) = TERNRY(CONSTS(1).EQ.1.00000, CONSTS(162)*CONSTS(161)*0.900000, TERNRY(CONSTS(1).EQ.2.00000, CONSTS(162)*CONSTS(161)*0.700000, CONSTS(162)*CONSTS(161))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(50) = CONSTS(231)
RATES(41) = (- ( CONSTS(80)*EXP( CONSTS(81)*STATES(1))*STATES(41)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(82)))/10.0000) - CONSTS(83)*EXP( CONSTS(84)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(85)))/10.0000))+ CONSTS(92)*EXP( CONSTS(93)*STATES(1))*STATES(43)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(94)))/10.0000)) - CONSTS(101)*EXP( CONSTS(102)*STATES(1))*STATES(41)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(103)))/10.0000)
RATES(42) = ((( CONSTS(80)*EXP( CONSTS(81)*STATES(1))*STATES(41)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(82)))/10.0000) - CONSTS(83)*EXP( CONSTS(84)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(85)))/10.0000)) - ( CONSTS(74)*EXP( CONSTS(75)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(76)))/10.0000) - CONSTS(77)*EXP( CONSTS(78)*STATES(1))*STATES(46)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(79)))/10.0000)))+ CONSTS(95)*EXP( CONSTS(96)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(97)))/10.0000)) - CONSTS(104)*EXP( CONSTS(105)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(106)))/10.0000)
RATES(43) = - ( CONSTS(68)*EXP( CONSTS(69)*STATES(1))*STATES(43)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(70)))/10.0000) - CONSTS(71)*EXP( CONSTS(72)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(73)))/10.0000)) - ( CONSTS(92)*EXP( CONSTS(93)*STATES(1))*STATES(43)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(94)))/10.0000) - CONSTS(101)*EXP( CONSTS(102)*STATES(1))*STATES(41)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(103)))/10.0000))
RATES(44) = (( CONSTS(68)*EXP( CONSTS(69)*STATES(1))*STATES(43)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(70)))/10.0000) - CONSTS(71)*EXP( CONSTS(72)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(73)))/10.0000)) - ( CONSTS(86)*EXP( CONSTS(87)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(88)))/10.0000) - CONSTS(89)*EXP( CONSTS(90)*STATES(1))*STATES(45)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(91)))/10.0000))) - ( CONSTS(95)*EXP( CONSTS(96)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(97)))/10.0000) - CONSTS(104)*EXP( CONSTS(105)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(106)))/10.0000))
RATES(45) = (( CONSTS(86)*EXP( CONSTS(87)*STATES(1))*STATES(44)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(88)))/10.0000) - CONSTS(89)*EXP( CONSTS(90)*STATES(1))*STATES(45)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(91)))/10.0000)) - ( CONSTS(98)*EXP( CONSTS(99)*STATES(1))*STATES(45)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(100)))/10.0000) - CONSTS(107)*EXP( CONSTS(108)*STATES(1))*STATES(46)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(109)))/10.0000))) - ( (( CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*log(STATES(50))))/(EXP( CONSTS(112)*log(STATES(50)))+CONSTS(113)))*STATES(45) - CONSTS(111)*STATES(48))
RATES(46) = ((( CONSTS(74)*EXP( CONSTS(75)*STATES(1))*STATES(42)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(76)))/10.0000) - CONSTS(77)*EXP( CONSTS(78)*STATES(1))*STATES(46)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(79)))/10.0000))+ CONSTS(98)*EXP( CONSTS(99)*STATES(1))*STATES(45)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(100)))/10.0000)) - CONSTS(107)*EXP( CONSTS(108)*STATES(1))*STATES(46)*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(109)))/10.0000)) - ( (( CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*log(STATES(50))))/(EXP( CONSTS(112)*log(STATES(50)))+CONSTS(113)))*STATES(46) - (( CONSTS(111)*CONSTS(98)*EXP( CONSTS(99)*STATES(1))*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(100)))/10.0000))/( CONSTS(107)*EXP( CONSTS(108)*STATES(1))*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(109)))/10.0000)))*STATES(47))
RATES(47) = (( (( CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*log(STATES(50))))/(EXP( CONSTS(112)*log(STATES(50)))+CONSTS(113)))*STATES(46) - (( CONSTS(111)*CONSTS(98)*EXP( CONSTS(99)*STATES(1))*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(100)))/10.0000))/( CONSTS(107)*EXP( CONSTS(108)*STATES(1))*EXP(( (CONSTS(116) - 293.000)*log(CONSTS(109)))/10.0000)))*STATES(47))+ (CONSTS(114)/(1.00000+EXP(- (STATES(1) - CONSTS(115))/6.78900)))*STATES(49)) - CONSTS(114)*STATES(47)
RATES(48) = (( (( CONSTS(110)*CONSTS(111)*EXP( CONSTS(112)*log(STATES(50))))/(EXP( CONSTS(112)*log(STATES(50)))+CONSTS(113)))*STATES(45) - CONSTS(111)*STATES(48))+ (CONSTS(114)/(1.00000+EXP(- (STATES(1) - CONSTS(115))/6.78900)))*STATES(49)) - CONSTS(114)*STATES(48)
RATES(49) = - ( (CONSTS(114)/(1.00000+EXP(- (STATES(1) - CONSTS(115))/6.78900)))*STATES(49) - CONSTS(114)*STATES(48)) - ( (CONSTS(114)/(1.00000+EXP(- (STATES(1) - CONSTS(115))/6.78900)))*STATES(49) - CONSTS(114)*STATES(47))
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
RATES(17) = (ALGBRC(3) - STATES(17))/CONSTS(177)
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
RATES(18) = (ALGBRC(4) - STATES(18))/CONSTS(199)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+19.5800+25.0000)/3.69600))
RATES(26) = (ALGBRC(7) - STATES(26))/CONSTS(60)
ALGBRC(12) = CONSTS(169) - CONSTS(170)/(1.00000+EXP(( 1000.00*STATES(3) - 0.0430000)/0.00820000))
RATES(54) = (ALGBRC(12) - STATES(54))/CONSTS(189)
ALGBRC(13) = 1.00000 - 1.00000/(1.00000+EXP(( 1000.00*STATES(3) - (STATES(54)+CONSTS(187)))/0.00300000))
RATES(55) = (ALGBRC(13) - STATES(55))/CONSTS(190)
ALGBRC(14) = 1.00000/(1.00000+EXP(( 1000.00*STATES(3) - (STATES(54)+CONSTS(188)))/0.00100000))
RATES(56) = (ALGBRC(14) - STATES(56))/CONSTS(191)
RATES(57) = (ALGBRC(14) - STATES(57))/CONSTS(202)
ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
ALGBRC(16) = 1.00000/( CONSTS(45)*EXP((STATES(1)+CONSTS(43))/CONSTS(44))+ CONSTS(46)*EXP(- (STATES(1)+CONSTS(47))/CONSTS(48)))
RATES(10) = (ALGBRC(1) - STATES(10))/ALGBRC(16)
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(49))/CONSTS(50)))
ALGBRC(17) = 1.00000/( 3.68600e-06*EXP(- (STATES(1)+3.88750)/7.85790)+ 16.0000*EXP((STATES(1) - 0.496300)/9.18430))+0.0750000
RATES(11) = (ALGBRC(2) - STATES(11))/ALGBRC(17)
ALGBRC(18) = 1.00000/( 0.00979400*EXP(- (STATES(1)+17.9500)/28.0500)+ 0.334300*EXP((STATES(1)+5.73000)/56.6600))
RATES(12) = (ALGBRC(2) - STATES(12))/ALGBRC(18)
ALGBRC(5) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 14.3400)/14.8200))
ALGBRC(20) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+CONSTS(58)+100.000)/29.3814)))
RATES(19) = (ALGBRC(5) - STATES(19))/ALGBRC(20)
ALGBRC(8) = STATES(26)*150.000
ALGBRC(22) = (1.00000 - STATES(25))/1.00000+CONSTS(61)/STATES(3) ** 4.00000
RATES(25) = ALGBRC(22)*CONSTS(62) - STATES(25)*ALGBRC(8)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000+CONSTS(120))/8.93200))
ALGBRC(25) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800+CONSTS(120))/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000+CONSTS(120))/230.000))
RATES(51) = (ALGBRC(10) - STATES(51))/ALGBRC(25)
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(4)+144.590+CONSTS(124))/( 1.56920*CONSTS(4)+3.81150)))
ALGBRC(26) = 122.200/(EXP(- (STATES(1)+CONSTS(124)+127.200)/20.3600)+EXP((STATES(1)+CONSTS(124)+236.800)/69.3300))
RATES(53) = (ALGBRC(11) - STATES(53))/ALGBRC(26)
ALGBRC(19) = ALGBRC(2)
ALGBRC(27) = (4.85900+1.00000/( 0.862800*EXP(- (STATES(1)+116.726)/7.60050)+ 1.10960*EXP((STATES(1)+6.27190)/9.03580)))*CONSTS(39)
RATES(13) = (ALGBRC(19) - STATES(13))/ALGBRC(27)
ALGBRC(31) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 24.3400)/14.8200))
RATES(22) = (ALGBRC(31) - STATES(22))/ALGBRC(20)
ALGBRC(24) = ALGBRC(10)
ALGBRC(33) = 1.00000/( 0.0100000*EXP(((STATES(1) - 50.0000)+CONSTS(120))/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400+CONSTS(120))/31.0000))
RATES(52) = (ALGBRC(24) - STATES(52))/ALGBRC(33)
ALGBRC(38) = ( CONSTS(22)*(1.00000 - STATES(2)))/(1.00000+CONSTS(23)/STATES(3))
RATES(2) = CONSTS(20)*ALGBRC(38)*(ALGBRC(38)+STATES(2)) - CONSTS(21)*STATES(2)
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(1)+84.7000)/6.22000))
ALGBRC(34) = 3.00000*ALGBRC(18)
RATES(14) = (ALGBRC(28) - STATES(14))/ALGBRC(34)
ALGBRC(35) = 1.46000*ALGBRC(27)
RATES(15) = (ALGBRC(19) - STATES(15))/ALGBRC(35)
ALGBRC(29) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
ALGBRC(36) = ALGBRC(16)
RATES(16) = (ALGBRC(29) - STATES(16))/ALGBRC(36)
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(58)+43.9400)/5.71100))
ALGBRC(21) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+CONSTS(58)+70.0000)/5.00000)), 1.00000)
ALGBRC(30) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+CONSTS(58)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+CONSTS(58)+50.0000)/16.5900))
ALGBRC(39) = ALGBRC(30)*ALGBRC(21)
RATES(20) = (ALGBRC(6) - STATES(20))/ALGBRC(39)
ALGBRC(37) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+CONSTS(58)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+CONSTS(58)+114.100)/8.07900))
ALGBRC(41) = ALGBRC(37)*ALGBRC(21)
RATES(21) = (ALGBRC(6) - STATES(21))/ALGBRC(41)
ALGBRC(43) = 1.35400+0.000100000/(EXP(((STATES(1)+CONSTS(58)) - 167.400)/15.8900)+EXP(- ((STATES(1)+CONSTS(58)) - 12.2300)/0.215400))
ALGBRC(45) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+CONSTS(58)+70.0000)/20.0000))
ALGBRC(47) = ALGBRC(43)*ALGBRC(45)*ALGBRC(39)
RATES(23) = (ALGBRC(6) - STATES(23))/ALGBRC(47)
ALGBRC(48) = ALGBRC(43)*ALGBRC(45)*ALGBRC(41)
RATES(24) = (ALGBRC(6) - STATES(24))/ALGBRC(48)
ALGBRC(65) = (( CONSTS(180)*STATES(25))/(1.00000 - STATES(25)))*(1.00000 - CONSTS(59))
ALGBRC(66) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000)/4.23000))
ALGBRC(67) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
ALGBRC(68) = ALGBRC(66)/ALGBRC(67)
ALGBRC(69) = (1.00000 - ALGBRC(66))/ALGBRC(67)
ALGBRC(70) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
ALGBRC(71) = ALGBRC(70)
ALGBRC(75) = 35.0000+ 350.000*EXP(- STATES(1)+20.0000 ** 2.00000/( 2.00000*100.000))
ALGBRC(76) = ALGBRC(75)
ALGBRC(80) = ALGBRC(71)/ALGBRC(76)
ALGBRC(84) = (1.00000 - ALGBRC(71))/ALGBRC(76)
RATES(29) = ((( ALGBRC(69)*STATES(39)+ ALGBRC(80)*STATES(28)) - (ALGBRC(84)+ALGBRC(68))*STATES(29)) - ALGBRC(65)*STATES(29))+ CONSTS(180)*STATES(35)
ALGBRC(73) = ALGBRC(70)
ALGBRC(78) = ALGBRC(75)
ALGBRC(81) = ALGBRC(73)/ALGBRC(78)
ALGBRC(85) = (1.00000 - ALGBRC(73))/ALGBRC(78)
RATES(32) = ((( ALGBRC(69)*STATES(40)+ ALGBRC(81)*STATES(31)) - (ALGBRC(85)+ALGBRC(68))*STATES(32)) - ALGBRC(65)*STATES(32))+ CONSTS(180)*STATES(38)
ALGBRC(88) = 0.800000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))+0.200000
ALGBRC(89) = 1.00000*(70.0000+1.20000/( 0.00450000*EXP((STATES(1)+20.0000)/- 50.0000)+ 0.00450000*EXP((STATES(1)+30.0000)/10.0000)))
ALGBRC(90) = (1.00000 - ALGBRC(88))/ALGBRC(89)
ALGBRC(91) = ALGBRC(88)/ALGBRC(89)
ALGBRC(102) = 1.00000*(100.000+0.00000/( 0.00350000*EXP((STATES(1)+5.00000)/- 84.0000)+ 0.00350000*EXP((STATES(1)+5.00000)/4.00000)))
ALGBRC(104) = ALGBRC(102)
ALGBRC(105) = (( ALGBRC(68)*ALGBRC(90)*ALGBRC(80))/ALGBRC(104))/( ALGBRC(68)*ALGBRC(90)*ALGBRC(80)+ ALGBRC(69)*ALGBRC(91)*ALGBRC(84))
ALGBRC(108) = 1.00000/ALGBRC(104) - ALGBRC(105)
RATES(28) = ((( ALGBRC(108)*STATES(27)+ ALGBRC(84)*STATES(29)) - (ALGBRC(105)+ALGBRC(80))*STATES(28)) - ALGBRC(65)*STATES(28))+ CONSTS(180)*STATES(34)
RATES(27) = ((( ALGBRC(105)*STATES(28)+ ALGBRC(90)*STATES(39)) - (ALGBRC(108)+ALGBRC(91))*STATES(27)) - ALGBRC(65)*STATES(27))+ CONSTS(180)*STATES(33)
ALGBRC(94) = ALGBRC(90)/CONSTS(181)
ALGBRC(95) = ALGBRC(91)/CONSTS(181)
ALGBRC(103) = ALGBRC(102)*CONSTS(181)
ALGBRC(106) = (( ALGBRC(68)*ALGBRC(94)*ALGBRC(81))/ALGBRC(103))/( ALGBRC(68)*ALGBRC(94)*ALGBRC(81)+ ALGBRC(69)*ALGBRC(95)*ALGBRC(85))
ALGBRC(109) = 1.00000/ALGBRC(103) - ALGBRC(106)
RATES(31) = ((( ALGBRC(109)*STATES(30)+ ALGBRC(85)*STATES(32)) - (ALGBRC(106)+ALGBRC(81))*STATES(31)) - ALGBRC(65)*STATES(31))+ CONSTS(180)*STATES(37)
RATES(30) = ((( ALGBRC(106)*STATES(31)+ ALGBRC(94)*STATES(40)) - (ALGBRC(109)+ALGBRC(95))*STATES(30)) - ALGBRC(65)*STATES(30))+ CONSTS(180)*STATES(36)
ALGBRC(72) = ALGBRC(70)
ALGBRC(77) = ALGBRC(75)
ALGBRC(82) = ALGBRC(72)/ALGBRC(77)
ALGBRC(86) = (1.00000 - ALGBRC(72))/ALGBRC(77)
ALGBRC(96) = ALGBRC(90)*CONSTS(64)
ALGBRC(97) = ALGBRC(91)*CONSTS(64)
ALGBRC(107) = ALGBRC(104)/CONSTS(64)
ALGBRC(113) = (( ALGBRC(68)*ALGBRC(96)*ALGBRC(82))/ALGBRC(107))/( ALGBRC(68)*ALGBRC(96)*ALGBRC(82)+ ALGBRC(69)*ALGBRC(97)*ALGBRC(86))
ALGBRC(115) = 1.00000/ALGBRC(107) - ALGBRC(113)
RATES(34) = ((( ALGBRC(115)*STATES(33)+ ALGBRC(86)*STATES(35)) - (ALGBRC(113)+ALGBRC(82))*STATES(34))+ ALGBRC(65)*STATES(28)) - CONSTS(180)*STATES(34)
ALGBRC(74) = ALGBRC(70)
ALGBRC(79) = ALGBRC(75)
ALGBRC(83) = ALGBRC(74)/ALGBRC(79)
ALGBRC(87) = (1.00000 - ALGBRC(74))/ALGBRC(79)
ALGBRC(100) = ALGBRC(94)*CONSTS(64)
ALGBRC(101) = ALGBRC(95)*CONSTS(64)
ALGBRC(110) = (ALGBRC(104)/CONSTS(64))*CONSTS(181)
ALGBRC(114) = (( ALGBRC(68)*ALGBRC(100)*ALGBRC(83))/ALGBRC(110))/( ALGBRC(68)*ALGBRC(100)*ALGBRC(83)+ ALGBRC(69)*ALGBRC(101)*ALGBRC(87))
ALGBRC(116) = 1.00000/ALGBRC(110) - ALGBRC(114)
RATES(37) = ((( ALGBRC(116)*STATES(36)+ ALGBRC(87)*STATES(38)) - (ALGBRC(114)+ALGBRC(83))*STATES(37))+ ALGBRC(65)*STATES(31)) - CONSTS(180)*STATES(37)
ALGBRC(122) = ((((((1.00000 - STATES(35)) - STATES(33)) - STATES(34)) - STATES(29)) - STATES(27)) - STATES(28)) - STATES(39)
RATES(39) = ((( ALGBRC(68)*STATES(29)+ ALGBRC(91)*STATES(27)) - (ALGBRC(69)+ALGBRC(90))*STATES(39)) - ALGBRC(65)*STATES(39))+ CONSTS(180)*ALGBRC(122)
RATES(33) = ((( ALGBRC(113)*STATES(34)+ ALGBRC(96)*ALGBRC(122)) - (ALGBRC(115)+ALGBRC(97))*STATES(33))+ ALGBRC(65)*STATES(27)) - CONSTS(180)*STATES(33)
RATES(35) = ((( ALGBRC(69)*ALGBRC(122)+ ALGBRC(82)*STATES(34)) - (ALGBRC(86)+ALGBRC(68))*STATES(35))+ ALGBRC(65)*STATES(29)) - CONSTS(180)*STATES(35)
ALGBRC(123) = ((((((1.00000 - STATES(38)) - STATES(36)) - STATES(37)) - STATES(32)) - STATES(30)) - STATES(31)) - STATES(40)
RATES(40) = ((( ALGBRC(68)*STATES(32)+ ALGBRC(95)*STATES(30)) - (ALGBRC(69)+ALGBRC(94))*STATES(40)) - ALGBRC(65)*STATES(40))+ CONSTS(180)*ALGBRC(123)
RATES(36) = ((( ALGBRC(114)*STATES(37)+ ALGBRC(100)*ALGBRC(123)) - (ALGBRC(116)+ALGBRC(101))*STATES(36))+ ALGBRC(65)*STATES(30)) - CONSTS(180)*STATES(36)
RATES(38) = ((( ALGBRC(69)*ALGBRC(123)+ ALGBRC(83)*STATES(37)) - (ALGBRC(87)+ALGBRC(68))*STATES(38))+ ALGBRC(65)*STATES(32)) - CONSTS(180)*STATES(38)
ALGBRC(50) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(4)/STATES(6))
ALGBRC(58) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(58)) - 213.600)/151.200))
ALGBRC(59) = 1.00000 - ALGBRC(58)
ALGBRC(60) = ALGBRC(58)*STATES(20)+ ALGBRC(59)*STATES(21)
ALGBRC(61) = ALGBRC(58)*STATES(23)+ ALGBRC(59)*STATES(24)
ALGBRC(40) = ALGBRC(38)+STATES(2)
ALGBRC(62) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(63) = CONSTS(179)*(STATES(1) - ALGBRC(50))*( (1.00000 - ALGBRC(62))*STATES(19)*ALGBRC(60)+ ALGBRC(62)*STATES(22)*ALGBRC(61))
ALGBRC(148) = CONSTS(200)*STATES(45)*(STATES(1) - ALGBRC(50))
ALGBRC(51) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log((CONSTS(4)+ CONSTS(37)*CONSTS(2))/(STATES(6)+ CONSTS(37)*STATES(4)))
ALGBRC(149) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(9) ** 1.40000)
ALGBRC(150) = CONSTS(183)*ALGBRC(149)*STATES(51)*STATES(52)*(STATES(1) - ALGBRC(51))
ALGBRC(151) = 1.00000/(1.00000+EXP((((STATES(1)+105.800) - 2.60000*CONSTS(4))+CONSTS(124))/( CONSTS(125)*9.49300)))
ALGBRC(152) = CONSTS(184)* CONSTS(4) ** (1.0 / 2)*ALGBRC(151)*STATES(53)*(STATES(1) - ALGBRC(50))
ALGBRC(216) = CONSTS(149)*EXP(( (1.00000 - CONSTS(150))*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
ALGBRC(220) = ( CONSTS(144)*CONSTS(4)/CONSTS(152) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(216) ** 3.00000+1.00000+CONSTS(4)/CONSTS(152) ** 2.00000) - 1.00000)
ALGBRC(217) = CONSTS(157)/(1.00000+CONSTS(156)/CONSTS(158)+STATES(4)/CONSTS(159)+STATES(6)/CONSTS(160))
ALGBRC(221) = ( CONSTS(145)*ALGBRC(217)*CONSTS(156))/(1.00000+CONSTS(154)/CONSTS(155))
ALGBRC(215) = CONSTS(148)*EXP(( CONSTS(150)*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
ALGBRC(218) = ( CONSTS(140)*STATES(4)/ALGBRC(215) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(215) ** 3.00000+1.00000+STATES(6)/CONSTS(151) ** 2.00000) - 1.00000)
ALGBRC(219) = ( CONSTS(143)*CONSTS(2)/ALGBRC(216) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(216) ** 3.00000+1.00000+CONSTS(4)/CONSTS(152) ** 2.00000) - 1.00000)
ALGBRC(222) = ( CONSTS(147)*STATES(6)/CONSTS(151) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(215) ** 3.00000+1.00000+STATES(6)/CONSTS(151) ** 2.00000) - 1.00000)
ALGBRC(223) = CONSTS(230)*ALGBRC(218)*CONSTS(229)+ ALGBRC(219)*ALGBRC(222)*ALGBRC(221)+ CONSTS(229)*ALGBRC(222)*ALGBRC(221)+ ALGBRC(221)*ALGBRC(218)*CONSTS(229)
ALGBRC(224) = ALGBRC(219)*CONSTS(228)*ALGBRC(222)+ ALGBRC(218)*CONSTS(229)*ALGBRC(220)+ ALGBRC(220)*CONSTS(228)*ALGBRC(222)+ CONSTS(229)*ALGBRC(220)*ALGBRC(222)
ALGBRC(225) = CONSTS(229)*ALGBRC(220)*CONSTS(230)+ ALGBRC(221)*ALGBRC(219)*CONSTS(228)+ ALGBRC(219)*CONSTS(228)*CONSTS(230)+ ALGBRC(220)*CONSTS(230)*CONSTS(228)
ALGBRC(226) = ALGBRC(222)*ALGBRC(221)*ALGBRC(219)+ ALGBRC(220)*CONSTS(230)*ALGBRC(218)+ ALGBRC(219)*CONSTS(230)*ALGBRC(218)+ ALGBRC(221)*ALGBRC(219)*ALGBRC(218)
ALGBRC(227) = ALGBRC(223)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(228) = ALGBRC(224)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(231) = 3.00000*( ALGBRC(227)*ALGBRC(220) - ALGBRC(228)*ALGBRC(221))
ALGBRC(229) = ALGBRC(225)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(230) = ALGBRC(226)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(232) = 2.00000*( ALGBRC(230)*CONSTS(228) - ALGBRC(229)*ALGBRC(218))
ALGBRC(233) = CONSTS(231)*( CONSTS(8)*ALGBRC(231)+ CONSTS(10)*ALGBRC(232))
ALGBRC(234) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
ALGBRC(235) = CONSTS(186)*ALGBRC(234)*(STATES(1) - ALGBRC(50))
ALGBRC(9) = TERNRY(VOI.GE.CONSTS(14).AND.VOI.LE.CONSTS(15).AND.(VOI - CONSTS(14)) - INT((VOI - CONSTS(14))/CONSTS(17))*CONSTS(17).LE.CONSTS(18), CONSTS(16), 0.00000)
ALGBRC(237) = (STATES(7) - STATES(6))/2.00000
RATES(6) = ( - ((ALGBRC(63)+ALGBRC(148)+ALGBRC(150)+ALGBRC(152)+ALGBRC(235)+ALGBRC(9)) - 2.00000*ALGBRC(233))*CONSTS(36)*CONSTS(203))/( CONSTS(7)*CONSTS(204))+( ALGBRC(237)*CONSTS(207))/CONSTS(204)
ALGBRC(23) = ( STATES(1)*CONSTS(7)*CONSTS(7))/( CONSTS(5)*CONSTS(6))
ALGBRC(32) = ( STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6))
ALGBRC(121) = ( 1.00000*ALGBRC(23)*( 0.750000*STATES(7)*EXP( 1.00000*ALGBRC(32)) - 0.750000*CONSTS(4)))/(EXP( 1.00000*ALGBRC(32)) - 1.00000)
ALGBRC(134) = CONSTS(212)*ALGBRC(121)*STATES(39)
ALGBRC(136) = CONSTS(212)*ALGBRC(121)*ALGBRC(122)
ALGBRC(142) = ALGBRC(134)+ALGBRC(136)
ALGBRC(135) = CONSTS(214)*ALGBRC(121)*STATES(40)
ALGBRC(137) = CONSTS(214)*ALGBRC(121)*ALGBRC(123)
ALGBRC(143) = ALGBRC(135)+ALGBRC(137)
ALGBRC(64) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(147) = ( ALGBRC(143)*ALGBRC(64)+ ALGBRC(142)*(1.00000 - ALGBRC(64)))*CONSTS(66)
RATES(7) = ( - ALGBRC(147)*CONSTS(36)*CONSTS(203))/( CONSTS(7)*CONSTS(207)) - ALGBRC(237)
ALGBRC(49) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(2)/STATES(4))
ALGBRC(52) = CONSTS(51)*STATES(11)+ CONSTS(176)*STATES(12)
ALGBRC(53) = CONSTS(51)*STATES(11)+ CONSTS(176)*STATES(14)
ALGBRC(54) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(55) = CONSTS(52)*CONSTS(40)*(STATES(1) - ALGBRC(49))*STATES(10) ** 3.00000*( (1.00000 - ALGBRC(54))*ALGBRC(52)*STATES(13)+ ALGBRC(54)*ALGBRC(53)*STATES(15))
ALGBRC(56) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(57) = CONSTS(178)*(STATES(1) - ALGBRC(49))*STATES(16)*( (1.00000 - ALGBRC(56))*STATES(17)+ ALGBRC(56)*STATES(18))
ALGBRC(181) = 1.00000/(1.00000+CONSTS(137)/STATES(9) ** 2.00000)
ALGBRC(154) = EXP(( CONSTS(135)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
ALGBRC(161) = 1.00000+ (CONSTS(2)/CONSTS(128))*(1.00000+1.00000/ALGBRC(154))
ALGBRC(162) = CONSTS(2)/( CONSTS(128)*ALGBRC(154)*ALGBRC(161))
ALGBRC(165) = ALGBRC(162)*CONSTS(132)
ALGBRC(155) = 1.00000+ (STATES(4)/CONSTS(128))*(1.00000+ALGBRC(154))
ALGBRC(156) = ( STATES(4)*ALGBRC(154))/( CONSTS(128)*ALGBRC(155))
ALGBRC(168) = ALGBRC(156)*CONSTS(132)
ALGBRC(158) = 1.00000+ (STATES(4)/CONSTS(126))*(1.00000+STATES(4)/CONSTS(127))
ALGBRC(159) = ( STATES(4)*STATES(4))/( ALGBRC(158)*CONSTS(126)*CONSTS(127))
ALGBRC(171) = ALGBRC(159)*ALGBRC(156)*CONSTS(130)
ALGBRC(172) = ALGBRC(162)*CONSTS(216)*CONSTS(130)
ALGBRC(163) = 1.00000/ALGBRC(161)
ALGBRC(164) = ALGBRC(163)*CONSTS(131)
ALGBRC(166) = ALGBRC(164)+ALGBRC(165)
ALGBRC(153) = EXP(( CONSTS(136)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
ALGBRC(157) = 1.00000/ALGBRC(155)
ALGBRC(167) = ( ALGBRC(157)*CONSTS(131))/ALGBRC(153)
ALGBRC(169) = ALGBRC(167)+ALGBRC(168)
ALGBRC(160) = 1.00000/ALGBRC(158)
ALGBRC(170) = ALGBRC(160)*STATES(9)*CONSTS(133)
ALGBRC(173) = CONSTS(219)*ALGBRC(169)*(ALGBRC(171)+ALGBRC(170))+ CONSTS(220)*ALGBRC(171)*(CONSTS(219)+ALGBRC(166))
ALGBRC(174) = CONSTS(218)*ALGBRC(171)*(ALGBRC(169)+CONSTS(220))+ ALGBRC(169)*ALGBRC(170)*(CONSTS(218)+ALGBRC(172))
ALGBRC(175) = CONSTS(218)*ALGBRC(166)*(ALGBRC(171)+ALGBRC(170))+ ALGBRC(172)*ALGBRC(170)*(CONSTS(219)+ALGBRC(166))
ALGBRC(176) = CONSTS(219)*ALGBRC(172)*(ALGBRC(169)+CONSTS(220))+ ALGBRC(166)*CONSTS(220)*(CONSTS(218)+ALGBRC(172))
ALGBRC(177) = ALGBRC(173)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(178) = ALGBRC(174)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(179) = ALGBRC(175)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(180) = ALGBRC(176)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(182) = ( 3.00000*( ALGBRC(180)*ALGBRC(171) - ALGBRC(177)*ALGBRC(172))+ ALGBRC(179)*ALGBRC(168)) - ALGBRC(178)*ALGBRC(165)
ALGBRC(183) = ALGBRC(178)*CONSTS(219) - ALGBRC(177)*CONSTS(218)
ALGBRC(184) = 0.800000*CONSTS(221)*ALGBRC(181)*( CONSTS(8)*ALGBRC(182)+ CONSTS(9)*ALGBRC(183))
ALGBRC(236) = ( CONSTS(164)*ALGBRC(23)*( STATES(4)*EXP(ALGBRC(32)) - CONSTS(2)))/(EXP(ALGBRC(32)) - 1.00000)
ALGBRC(239) = (STATES(5) - STATES(4))/2.00000
RATES(4) = ( - (ALGBRC(55)+ALGBRC(57)+ 3.00000*ALGBRC(184)+ 3.00000*ALGBRC(233)+ALGBRC(236))*CONSTS(203)*CONSTS(36))/( CONSTS(7)*CONSTS(204))+( ALGBRC(239)*CONSTS(207))/CONSTS(204)
ALGBRC(120) = ( 1.00000*ALGBRC(23)*( 0.750000*STATES(5)*EXP( 1.00000*ALGBRC(32)) - 0.750000*CONSTS(2)))/(EXP( 1.00000*ALGBRC(32)) - 1.00000)
ALGBRC(130) = CONSTS(211)*ALGBRC(120)*STATES(39)
ALGBRC(132) = CONSTS(211)*ALGBRC(120)*ALGBRC(122)
ALGBRC(140) = ALGBRC(130)+ALGBRC(132)
ALGBRC(131) = CONSTS(213)*ALGBRC(120)*STATES(40)
ALGBRC(133) = CONSTS(213)*ALGBRC(120)*ALGBRC(123)
ALGBRC(141) = ALGBRC(131)+ALGBRC(133)
ALGBRC(146) = ( ALGBRC(141)*ALGBRC(64)+ ALGBRC(140)*(1.00000 - ALGBRC(64)))*CONSTS(66)
ALGBRC(211) = 1.00000/(1.00000+CONSTS(137)/STATES(3) ** 2.00000)
ALGBRC(191) = 1.00000+ (CONSTS(2)/CONSTS(128))*(1.00000+1.00000/ALGBRC(154))
ALGBRC(192) = CONSTS(2)/( CONSTS(128)*ALGBRC(154)*ALGBRC(191))
ALGBRC(195) = ALGBRC(192)*CONSTS(132)
ALGBRC(185) = 1.00000+ (STATES(5)/CONSTS(128))*(1.00000+ALGBRC(154))
ALGBRC(186) = ( STATES(5)*ALGBRC(154))/( CONSTS(128)*ALGBRC(185))
ALGBRC(198) = ALGBRC(186)*CONSTS(132)
ALGBRC(188) = 1.00000+ (STATES(5)/CONSTS(126))*(1.00000+STATES(5)/CONSTS(127))
ALGBRC(189) = ( STATES(5)*STATES(5))/( ALGBRC(188)*CONSTS(126)*CONSTS(127))
ALGBRC(201) = ALGBRC(189)*ALGBRC(186)*CONSTS(130)
ALGBRC(202) = ALGBRC(192)*CONSTS(223)*CONSTS(130)
ALGBRC(193) = 1.00000/ALGBRC(191)
ALGBRC(194) = ALGBRC(193)*CONSTS(131)
ALGBRC(196) = ALGBRC(194)+ALGBRC(195)
ALGBRC(187) = 1.00000/ALGBRC(185)
ALGBRC(197) = ( ALGBRC(187)*CONSTS(131))/ALGBRC(153)
ALGBRC(199) = ALGBRC(197)+ALGBRC(198)
ALGBRC(190) = 1.00000/ALGBRC(188)
ALGBRC(200) = ALGBRC(190)*STATES(3)*CONSTS(133)
ALGBRC(203) = CONSTS(226)*ALGBRC(199)*(ALGBRC(201)+ALGBRC(200))+ CONSTS(227)*ALGBRC(201)*(CONSTS(226)+ALGBRC(196))
ALGBRC(204) = CONSTS(225)*ALGBRC(201)*(ALGBRC(199)+CONSTS(227))+ ALGBRC(199)*ALGBRC(200)*(CONSTS(225)+ALGBRC(202))
ALGBRC(205) = CONSTS(225)*ALGBRC(196)*(ALGBRC(201)+ALGBRC(200))+ ALGBRC(202)*ALGBRC(200)*(CONSTS(226)+ALGBRC(196))
ALGBRC(206) = CONSTS(226)*ALGBRC(202)*(ALGBRC(199)+CONSTS(227))+ ALGBRC(196)*CONSTS(227)*(CONSTS(225)+ALGBRC(202))
ALGBRC(207) = ALGBRC(203)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(208) = ALGBRC(204)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(209) = ALGBRC(205)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(210) = ALGBRC(206)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(212) = ( 3.00000*( ALGBRC(210)*ALGBRC(201) - ALGBRC(207)*ALGBRC(202))+ ALGBRC(209)*ALGBRC(198)) - ALGBRC(208)*ALGBRC(195)
ALGBRC(213) = ALGBRC(208)*CONSTS(226) - ALGBRC(207)*CONSTS(225)
ALGBRC(214) = 0.200000*CONSTS(221)*ALGBRC(211)*( CONSTS(8)*ALGBRC(212)+ CONSTS(9)*ALGBRC(213))
RATES(5) = ( - (ALGBRC(146)+ 3.00000*ALGBRC(214))*CONSTS(36)*CONSTS(203))/( CONSTS(7)*CONSTS(207)) - ALGBRC(239)
ALGBRC(119) = ( 4.00000*ALGBRC(23)*( 1.20000*STATES(3)*EXP( 2.00000*ALGBRC(32)) - 0.341000*CONSTS(3)))/(EXP( 2.00000*ALGBRC(32)) - 1.00000)
ALGBRC(124) = CONSTS(209)*ALGBRC(119)*STATES(39)
ALGBRC(127) = CONSTS(209)*ALGBRC(119)*ALGBRC(122)
ALGBRC(138) = ALGBRC(124)+ALGBRC(127)
ALGBRC(125) = CONSTS(210)*ALGBRC(119)*STATES(40)
ALGBRC(128) = CONSTS(210)*ALGBRC(119)*ALGBRC(123)
ALGBRC(139) = ALGBRC(125)+ALGBRC(128)
ALGBRC(144) = ( ALGBRC(139)*ALGBRC(64)+ ALGBRC(138)*(1.00000 - ALGBRC(64)))*CONSTS(66)
ALGBRC(240) = ( CONSTS(167)*STATES(9))/(CONSTS(168)+STATES(9))
ALGBRC(238) = ( (1.00000 - CONSTS(166))*CONSTS(165)*16.0000*ALGBRC(23)*( 1.20000*STATES(9)*EXP( 2.00000*ALGBRC(32)) - 0.341000*CONSTS(3)))/(EXP( 2.00000*ALGBRC(32)) - 1.00000)
RATES(1) = - (ALGBRC(55)+ALGBRC(57)+ALGBRC(63)+ALGBRC(144)+ALGBRC(146)+ALGBRC(147)+ALGBRC(148)+ALGBRC(150)+ALGBRC(152)+ALGBRC(184)+ALGBRC(214)+ALGBRC(233)+ALGBRC(236)+ALGBRC(235)+ALGBRC(240)+ALGBRC(238)+ALGBRC(9))
ALGBRC(241) = ( (STATES(3) - STATES(9))*1.70000)/0.200000
ALGBRC(243) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(242) = 1.00000 - 1.00000/(1.00000+EXP((STATES(8) - 0.300000)/0.100000))
ALGBRC(244) = TERNRY(CONSTS(1).EQ.2.00000, CONSTS(192)*ALGBRC(242)*STATES(55)*STATES(56)*(STATES(8) - STATES(3)), CONSTS(171)*ALGBRC(242)*STATES(55)*STATES(56)*(STATES(8) - STATES(3)))
ALGBRC(245) = TERNRY(CONSTS(1).EQ.2.00000, CONSTS(193)*1.70000*ALGBRC(242)*STATES(55)*STATES(57)*(STATES(8) - STATES(3)), CONSTS(193)*ALGBRC(242)*STATES(55)*STATES(57)*(STATES(8) - STATES(3)))
ALGBRC(246) = (1.00000 - ALGBRC(243))*ALGBRC(244)+ ALGBRC(243)*ALGBRC(245)
ALGBRC(44) = 1.00000/(1.00000+( CONSTS(28)*CONSTS(29))/CONSTS(29)+STATES(3) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(3) ** 2.00000)
RATES(3) = ALGBRC(44)*((( - (ALGBRC(144) - 2.00000*ALGBRC(214))*CONSTS(203))/( 2.00000*CONSTS(7)*CONSTS(207))+( ALGBRC(246)*CONSTS(208))/CONSTS(207)) - ALGBRC(241))
ALGBRC(247) = ( CONSTS(194)*0.00437500*STATES(9))/(STATES(9)+0.000920000)
ALGBRC(248) = ( CONSTS(194)*2.75000*0.00437500*STATES(9))/((STATES(9)+0.000920000) - 0.000170000)
ALGBRC(249) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(251) = CONSTS(173)*( (1.00000 - ALGBRC(249))*ALGBRC(247)+ ALGBRC(249)*ALGBRC(248))
ALGBRC(250) = ( 0.0123000*STATES(8))/15.0000
ALGBRC(42) = 1.00000/(1.00000+( CONSTS(174)*CONSTS(25))/CONSTS(25)+STATES(9) ** 2.00000+( CONSTS(26)*CONSTS(27))/CONSTS(27)+STATES(9) ** 2.00000)
RATES(9) = ALGBRC(42)*((( - ((ALGBRC(240)+ALGBRC(238)) - 2.00000*ALGBRC(184))*CONSTS(203))/( 2.00000*CONSTS(7)*CONSTS(204)) - ( ALGBRC(251)*CONSTS(208))/CONSTS(204))+( ALGBRC(250)*CONSTS(208))/CONSTS(204)+( ALGBRC(241)*CONSTS(207))/CONSTS(204))
ALGBRC(46) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(33))/CONSTS(33)+STATES(8) ** 2.00000)
RATES(8) = ALGBRC(46)*((ALGBRC(251) - ALGBRC(250)) - ALGBRC(246))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+87.6100)/7.48800))
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+93.8100)/7.48800))
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+19.5800+25.0000)/3.69600))
ALGBRC(12) = CONSTS(169) - CONSTS(170)/(1.00000+EXP(( 1000.00*STATES(3) - 0.0430000)/0.00820000))
ALGBRC(13) = 1.00000 - 1.00000/(1.00000+EXP(( 1000.00*STATES(3) - (STATES(54)+CONSTS(187)))/0.00300000))
ALGBRC(14) = 1.00000/(1.00000+EXP(( 1000.00*STATES(3) - (STATES(54)+CONSTS(188)))/0.00100000))
ALGBRC(1) = 1.00000/(1.00000+EXP(- (STATES(1)+CONSTS(41))/CONSTS(42)))
ALGBRC(16) = 1.00000/( CONSTS(45)*EXP((STATES(1)+CONSTS(43))/CONSTS(44))+ CONSTS(46)*EXP(- (STATES(1)+CONSTS(47))/CONSTS(48)))
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(49))/CONSTS(50)))
ALGBRC(17) = 1.00000/( 3.68600e-06*EXP(- (STATES(1)+3.88750)/7.85790)+ 16.0000*EXP((STATES(1) - 0.496300)/9.18430))+0.0750000
ALGBRC(18) = 1.00000/( 0.00979400*EXP(- (STATES(1)+17.9500)/28.0500)+ 0.334300*EXP((STATES(1)+5.73000)/56.6600))
ALGBRC(5) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 14.3400)/14.8200))
ALGBRC(20) = 1.05150/(1.00000/( 1.20890*(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 18.4099)/29.3814)))+3.50000/(1.00000+EXP((STATES(1)+CONSTS(58)+100.000)/29.3814)))
ALGBRC(8) = STATES(26)*150.000
ALGBRC(22) = (1.00000 - STATES(25))/1.00000+CONSTS(61)/STATES(3) ** 4.00000
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+11.6000+CONSTS(120))/8.93200))
ALGBRC(25) = 817.300+1.00000/( 0.000232600*EXP((STATES(1)+48.2800+CONSTS(120))/17.8000)+ 0.00129200*EXP(- (STATES(1)+210.000+CONSTS(120))/230.000))
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+ 2.55380*CONSTS(4)+144.590+CONSTS(124))/( 1.56920*CONSTS(4)+3.81150)))
ALGBRC(26) = 122.200/(EXP(- (STATES(1)+CONSTS(124)+127.200)/20.3600)+EXP((STATES(1)+CONSTS(124)+236.800)/69.3300))
ALGBRC(19) = ALGBRC(2)
ALGBRC(27) = (4.85900+1.00000/( 0.862800*EXP(- (STATES(1)+116.726)/7.60050)+ 1.10960*EXP((STATES(1)+6.27190)/9.03580)))*CONSTS(39)
ALGBRC(31) = 1.00000/(1.00000+EXP(- ((STATES(1)+CONSTS(58)) - 24.3400)/14.8200))
ALGBRC(24) = ALGBRC(10)
ALGBRC(33) = 1.00000/( 0.0100000*EXP(((STATES(1) - 50.0000)+CONSTS(120))/20.0000)+ 0.0193000*EXP(- (STATES(1)+66.5400+CONSTS(120))/31.0000))
ALGBRC(38) = ( CONSTS(22)*(1.00000 - STATES(2)))/(1.00000+CONSTS(23)/STATES(3))
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(1)+84.7000)/6.22000))
ALGBRC(34) = 3.00000*ALGBRC(18)
ALGBRC(35) = 1.46000*ALGBRC(27)
ALGBRC(29) = 1.00000/(1.00000+EXP(- (STATES(1)+42.8500)/5.26400))
ALGBRC(36) = ALGBRC(16)
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+CONSTS(58)+43.9400)/5.71100))
ALGBRC(21) = TERNRY(CONSTS(1).EQ.1.00000, 1.00000 - 0.950000/(1.00000+EXP((STATES(1)+CONSTS(58)+70.0000)/5.00000)), 1.00000)
ALGBRC(30) = 4.56200+1.00000/( 0.393300*EXP(- (STATES(1)+CONSTS(58)+100.000)/100.000)+ 0.0800400*EXP((STATES(1)+CONSTS(58)+50.0000)/16.5900))
ALGBRC(39) = ALGBRC(30)*ALGBRC(21)
ALGBRC(37) = 23.6200+1.00000/( 0.00141600*EXP(- (STATES(1)+CONSTS(58)+96.5200)/59.0500)+ 1.78000e-08*EXP((STATES(1)+CONSTS(58)+114.100)/8.07900))
ALGBRC(41) = ALGBRC(37)*ALGBRC(21)
ALGBRC(43) = 1.35400+0.000100000/(EXP(((STATES(1)+CONSTS(58)) - 167.400)/15.8900)+EXP(- ((STATES(1)+CONSTS(58)) - 12.2300)/0.215400))
ALGBRC(45) = 1.00000 - 0.500000/(1.00000+EXP((STATES(1)+CONSTS(58)+70.0000)/20.0000))
ALGBRC(47) = ALGBRC(43)*ALGBRC(45)*ALGBRC(39)
ALGBRC(48) = ALGBRC(43)*ALGBRC(45)*ALGBRC(41)
ALGBRC(65) = (( CONSTS(180)*STATES(25))/(1.00000 - STATES(25)))*(1.00000 - CONSTS(59))
ALGBRC(66) = 1.00000/(1.00000+EXP(- (STATES(1)+3.94000)/4.23000))
ALGBRC(67) = 0.600000+1.00000/(EXP( - 0.0500000*(STATES(1)+6.00000))+EXP( 0.0900000*(STATES(1)+14.0000)))
ALGBRC(68) = ALGBRC(66)/ALGBRC(67)
ALGBRC(69) = (1.00000 - ALGBRC(66))/ALGBRC(67)
ALGBRC(70) = 1.00000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))
ALGBRC(71) = ALGBRC(70)
ALGBRC(75) = 35.0000+ 350.000*EXP(- STATES(1)+20.0000 ** 2.00000/( 2.00000*100.000))
ALGBRC(76) = ALGBRC(75)
ALGBRC(80) = ALGBRC(71)/ALGBRC(76)
ALGBRC(84) = (1.00000 - ALGBRC(71))/ALGBRC(76)
ALGBRC(73) = ALGBRC(70)
ALGBRC(78) = ALGBRC(75)
ALGBRC(81) = ALGBRC(73)/ALGBRC(78)
ALGBRC(85) = (1.00000 - ALGBRC(73))/ALGBRC(78)
ALGBRC(88) = 0.800000/(1.00000+EXP((STATES(1)+19.5800)/3.69600))+0.200000
ALGBRC(89) = 1.00000*(70.0000+1.20000/( 0.00450000*EXP((STATES(1)+20.0000)/- 50.0000)+ 0.00450000*EXP((STATES(1)+30.0000)/10.0000)))
ALGBRC(90) = (1.00000 - ALGBRC(88))/ALGBRC(89)
ALGBRC(91) = ALGBRC(88)/ALGBRC(89)
ALGBRC(102) = 1.00000*(100.000+0.00000/( 0.00350000*EXP((STATES(1)+5.00000)/- 84.0000)+ 0.00350000*EXP((STATES(1)+5.00000)/4.00000)))
ALGBRC(104) = ALGBRC(102)
ALGBRC(105) = (( ALGBRC(68)*ALGBRC(90)*ALGBRC(80))/ALGBRC(104))/( ALGBRC(68)*ALGBRC(90)*ALGBRC(80)+ ALGBRC(69)*ALGBRC(91)*ALGBRC(84))
ALGBRC(108) = 1.00000/ALGBRC(104) - ALGBRC(105)
ALGBRC(94) = ALGBRC(90)/CONSTS(181)
ALGBRC(95) = ALGBRC(91)/CONSTS(181)
ALGBRC(103) = ALGBRC(102)*CONSTS(181)
ALGBRC(106) = (( ALGBRC(68)*ALGBRC(94)*ALGBRC(81))/ALGBRC(103))/( ALGBRC(68)*ALGBRC(94)*ALGBRC(81)+ ALGBRC(69)*ALGBRC(95)*ALGBRC(85))
ALGBRC(109) = 1.00000/ALGBRC(103) - ALGBRC(106)
ALGBRC(72) = ALGBRC(70)
ALGBRC(77) = ALGBRC(75)
ALGBRC(82) = ALGBRC(72)/ALGBRC(77)
ALGBRC(86) = (1.00000 - ALGBRC(72))/ALGBRC(77)
ALGBRC(96) = ALGBRC(90)*CONSTS(64)
ALGBRC(97) = ALGBRC(91)*CONSTS(64)
ALGBRC(107) = ALGBRC(104)/CONSTS(64)
ALGBRC(113) = (( ALGBRC(68)*ALGBRC(96)*ALGBRC(82))/ALGBRC(107))/( ALGBRC(68)*ALGBRC(96)*ALGBRC(82)+ ALGBRC(69)*ALGBRC(97)*ALGBRC(86))
ALGBRC(115) = 1.00000/ALGBRC(107) - ALGBRC(113)
ALGBRC(74) = ALGBRC(70)
ALGBRC(79) = ALGBRC(75)
ALGBRC(83) = ALGBRC(74)/ALGBRC(79)
ALGBRC(87) = (1.00000 - ALGBRC(74))/ALGBRC(79)
ALGBRC(100) = ALGBRC(94)*CONSTS(64)
ALGBRC(101) = ALGBRC(95)*CONSTS(64)
ALGBRC(110) = (ALGBRC(104)/CONSTS(64))*CONSTS(181)
ALGBRC(114) = (( ALGBRC(68)*ALGBRC(100)*ALGBRC(83))/ALGBRC(110))/( ALGBRC(68)*ALGBRC(100)*ALGBRC(83)+ ALGBRC(69)*ALGBRC(101)*ALGBRC(87))
ALGBRC(116) = 1.00000/ALGBRC(110) - ALGBRC(114)
ALGBRC(122) = ((((((1.00000 - STATES(35)) - STATES(33)) - STATES(34)) - STATES(29)) - STATES(27)) - STATES(28)) - STATES(39)
ALGBRC(123) = ((((((1.00000 - STATES(38)) - STATES(36)) - STATES(37)) - STATES(32)) - STATES(30)) - STATES(31)) - STATES(40)
ALGBRC(50) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(4)/STATES(6))
ALGBRC(58) = 1.00000/(1.00000+EXP(((STATES(1)+CONSTS(58)) - 213.600)/151.200))
ALGBRC(59) = 1.00000 - ALGBRC(58)
ALGBRC(60) = ALGBRC(58)*STATES(20)+ ALGBRC(59)*STATES(21)
ALGBRC(61) = ALGBRC(58)*STATES(23)+ ALGBRC(59)*STATES(24)
ALGBRC(40) = ALGBRC(38)+STATES(2)
ALGBRC(62) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(63) = CONSTS(179)*(STATES(1) - ALGBRC(50))*( (1.00000 - ALGBRC(62))*STATES(19)*ALGBRC(60)+ ALGBRC(62)*STATES(22)*ALGBRC(61))
ALGBRC(148) = CONSTS(200)*STATES(45)*(STATES(1) - ALGBRC(50))
ALGBRC(51) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log((CONSTS(4)+ CONSTS(37)*CONSTS(2))/(STATES(6)+ CONSTS(37)*STATES(4)))
ALGBRC(149) = 1.00000+0.600000/(1.00000+3.80000e-05/STATES(9) ** 1.40000)
ALGBRC(150) = CONSTS(183)*ALGBRC(149)*STATES(51)*STATES(52)*(STATES(1) - ALGBRC(51))
ALGBRC(151) = 1.00000/(1.00000+EXP((((STATES(1)+105.800) - 2.60000*CONSTS(4))+CONSTS(124))/( CONSTS(125)*9.49300)))
ALGBRC(152) = CONSTS(184)* CONSTS(4) ** (1.0 / 2)*ALGBRC(151)*STATES(53)*(STATES(1) - ALGBRC(50))
ALGBRC(216) = CONSTS(149)*EXP(( (1.00000 - CONSTS(150))*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
ALGBRC(220) = ( CONSTS(144)*CONSTS(4)/CONSTS(152) ** 2.00000)/((1.00000+CONSTS(2)/ALGBRC(216) ** 3.00000+1.00000+CONSTS(4)/CONSTS(152) ** 2.00000) - 1.00000)
ALGBRC(217) = CONSTS(157)/(1.00000+CONSTS(156)/CONSTS(158)+STATES(4)/CONSTS(159)+STATES(6)/CONSTS(160))
ALGBRC(221) = ( CONSTS(145)*ALGBRC(217)*CONSTS(156))/(1.00000+CONSTS(154)/CONSTS(155))
ALGBRC(215) = CONSTS(148)*EXP(( CONSTS(150)*STATES(1)*CONSTS(7))/( 3.00000*CONSTS(5)*CONSTS(6)))
ALGBRC(218) = ( CONSTS(140)*STATES(4)/ALGBRC(215) ** 3.00000)/((1.00000+STATES(4)/ALGBRC(215) ** 3.00000+1.00000+STATES(6)/CONSTS(151) ** 2.00000) - 1.00000)
ALGBRC(219) = ( CONSTS(143)*CONSTS(2)/ALGBRC(216) ** 3.00000)/((1.00000+CONSTS(2)/ALGBRC(216) ** 3.00000+1.00000+CONSTS(4)/CONSTS(152) ** 2.00000) - 1.00000)
ALGBRC(222) = ( CONSTS(147)*STATES(6)/CONSTS(151) ** 2.00000)/((1.00000+STATES(4)/ALGBRC(215) ** 3.00000+1.00000+STATES(6)/CONSTS(151) ** 2.00000) - 1.00000)
ALGBRC(223) = CONSTS(230)*ALGBRC(218)*CONSTS(229)+ ALGBRC(219)*ALGBRC(222)*ALGBRC(221)+ CONSTS(229)*ALGBRC(222)*ALGBRC(221)+ ALGBRC(221)*ALGBRC(218)*CONSTS(229)
ALGBRC(224) = ALGBRC(219)*CONSTS(228)*ALGBRC(222)+ ALGBRC(218)*CONSTS(229)*ALGBRC(220)+ ALGBRC(220)*CONSTS(228)*ALGBRC(222)+ CONSTS(229)*ALGBRC(220)*ALGBRC(222)
ALGBRC(225) = CONSTS(229)*ALGBRC(220)*CONSTS(230)+ ALGBRC(221)*ALGBRC(219)*CONSTS(228)+ ALGBRC(219)*CONSTS(228)*CONSTS(230)+ ALGBRC(220)*CONSTS(230)*CONSTS(228)
ALGBRC(226) = ALGBRC(222)*ALGBRC(221)*ALGBRC(219)+ ALGBRC(220)*CONSTS(230)*ALGBRC(218)+ ALGBRC(219)*CONSTS(230)*ALGBRC(218)+ ALGBRC(221)*ALGBRC(219)*ALGBRC(218)
ALGBRC(227) = ALGBRC(223)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(228) = ALGBRC(224)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(231) = 3.00000*( ALGBRC(227)*ALGBRC(220) - ALGBRC(228)*ALGBRC(221))
ALGBRC(229) = ALGBRC(225)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(230) = ALGBRC(226)/(ALGBRC(223)+ALGBRC(224)+ALGBRC(225)+ALGBRC(226))
ALGBRC(232) = 2.00000*( ALGBRC(230)*CONSTS(228) - ALGBRC(229)*ALGBRC(218))
ALGBRC(233) = CONSTS(231)*( CONSTS(8)*ALGBRC(231)+ CONSTS(10)*ALGBRC(232))
ALGBRC(234) = 1.00000/(1.00000+EXP(- (STATES(1) - 14.4800)/18.3400))
ALGBRC(235) = CONSTS(186)*ALGBRC(234)*(STATES(1) - ALGBRC(50))
ALGBRC(9) = TERNRY(VOI.GE.CONSTS(14).AND.VOI.LE.CONSTS(15).AND.(VOI - CONSTS(14)) - INT((VOI - CONSTS(14))/CONSTS(17))*CONSTS(17).LE.CONSTS(18), CONSTS(16), 0.00000)
ALGBRC(237) = (STATES(7) - STATES(6))/2.00000
ALGBRC(23) = ( STATES(1)*CONSTS(7)*CONSTS(7))/( CONSTS(5)*CONSTS(6))
ALGBRC(32) = ( STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6))
ALGBRC(121) = ( 1.00000*ALGBRC(23)*( 0.750000*STATES(7)*EXP( 1.00000*ALGBRC(32)) - 0.750000*CONSTS(4)))/(EXP( 1.00000*ALGBRC(32)) - 1.00000)
ALGBRC(134) = CONSTS(212)*ALGBRC(121)*STATES(39)
ALGBRC(136) = CONSTS(212)*ALGBRC(121)*ALGBRC(122)
ALGBRC(142) = ALGBRC(134)+ALGBRC(136)
ALGBRC(135) = CONSTS(214)*ALGBRC(121)*STATES(40)
ALGBRC(137) = CONSTS(214)*ALGBRC(121)*ALGBRC(123)
ALGBRC(143) = ALGBRC(135)+ALGBRC(137)
ALGBRC(64) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(147) = ( ALGBRC(143)*ALGBRC(64)+ ALGBRC(142)*(1.00000 - ALGBRC(64)))*CONSTS(66)
ALGBRC(49) = (( CONSTS(5)*CONSTS(6))/CONSTS(7))*log(CONSTS(2)/STATES(4))
ALGBRC(52) = CONSTS(51)*STATES(11)+ CONSTS(176)*STATES(12)
ALGBRC(53) = CONSTS(51)*STATES(11)+ CONSTS(176)*STATES(14)
ALGBRC(54) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(55) = CONSTS(52)*CONSTS(40)*(STATES(1) - ALGBRC(49))*STATES(10) ** 3.00000*( (1.00000 - ALGBRC(54))*ALGBRC(52)*STATES(13)+ ALGBRC(54)*ALGBRC(53)*STATES(15))
ALGBRC(56) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(57) = CONSTS(178)*(STATES(1) - ALGBRC(49))*STATES(16)*( (1.00000 - ALGBRC(56))*STATES(17)+ ALGBRC(56)*STATES(18))
ALGBRC(181) = 1.00000/(1.00000+CONSTS(137)/STATES(9) ** 2.00000)
ALGBRC(154) = EXP(( CONSTS(135)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
ALGBRC(161) = 1.00000+ (CONSTS(2)/CONSTS(128))*(1.00000+1.00000/ALGBRC(154))
ALGBRC(162) = CONSTS(2)/( CONSTS(128)*ALGBRC(154)*ALGBRC(161))
ALGBRC(165) = ALGBRC(162)*CONSTS(132)
ALGBRC(155) = 1.00000+ (STATES(4)/CONSTS(128))*(1.00000+ALGBRC(154))
ALGBRC(156) = ( STATES(4)*ALGBRC(154))/( CONSTS(128)*ALGBRC(155))
ALGBRC(168) = ALGBRC(156)*CONSTS(132)
ALGBRC(158) = 1.00000+ (STATES(4)/CONSTS(126))*(1.00000+STATES(4)/CONSTS(127))
ALGBRC(159) = ( STATES(4)*STATES(4))/( ALGBRC(158)*CONSTS(126)*CONSTS(127))
ALGBRC(171) = ALGBRC(159)*ALGBRC(156)*CONSTS(130)
ALGBRC(172) = ALGBRC(162)*CONSTS(216)*CONSTS(130)
ALGBRC(163) = 1.00000/ALGBRC(161)
ALGBRC(164) = ALGBRC(163)*CONSTS(131)
ALGBRC(166) = ALGBRC(164)+ALGBRC(165)
ALGBRC(153) = EXP(( CONSTS(136)*STATES(1)*CONSTS(7))/( CONSTS(5)*CONSTS(6)))
ALGBRC(157) = 1.00000/ALGBRC(155)
ALGBRC(167) = ( ALGBRC(157)*CONSTS(131))/ALGBRC(153)
ALGBRC(169) = ALGBRC(167)+ALGBRC(168)
ALGBRC(160) = 1.00000/ALGBRC(158)
ALGBRC(170) = ALGBRC(160)*STATES(9)*CONSTS(133)
ALGBRC(173) = CONSTS(219)*ALGBRC(169)*(ALGBRC(171)+ALGBRC(170))+ CONSTS(220)*ALGBRC(171)*(CONSTS(219)+ALGBRC(166))
ALGBRC(174) = CONSTS(218)*ALGBRC(171)*(ALGBRC(169)+CONSTS(220))+ ALGBRC(169)*ALGBRC(170)*(CONSTS(218)+ALGBRC(172))
ALGBRC(175) = CONSTS(218)*ALGBRC(166)*(ALGBRC(171)+ALGBRC(170))+ ALGBRC(172)*ALGBRC(170)*(CONSTS(219)+ALGBRC(166))
ALGBRC(176) = CONSTS(219)*ALGBRC(172)*(ALGBRC(169)+CONSTS(220))+ ALGBRC(166)*CONSTS(220)*(CONSTS(218)+ALGBRC(172))
ALGBRC(177) = ALGBRC(173)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(178) = ALGBRC(174)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(179) = ALGBRC(175)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(180) = ALGBRC(176)/(ALGBRC(173)+ALGBRC(174)+ALGBRC(175)+ALGBRC(176))
ALGBRC(182) = ( 3.00000*( ALGBRC(180)*ALGBRC(171) - ALGBRC(177)*ALGBRC(172))+ ALGBRC(179)*ALGBRC(168)) - ALGBRC(178)*ALGBRC(165)
ALGBRC(183) = ALGBRC(178)*CONSTS(219) - ALGBRC(177)*CONSTS(218)
ALGBRC(184) = 0.800000*CONSTS(221)*ALGBRC(181)*( CONSTS(8)*ALGBRC(182)+ CONSTS(9)*ALGBRC(183))
ALGBRC(236) = ( CONSTS(164)*ALGBRC(23)*( STATES(4)*EXP(ALGBRC(32)) - CONSTS(2)))/(EXP(ALGBRC(32)) - 1.00000)
ALGBRC(239) = (STATES(5) - STATES(4))/2.00000
ALGBRC(120) = ( 1.00000*ALGBRC(23)*( 0.750000*STATES(5)*EXP( 1.00000*ALGBRC(32)) - 0.750000*CONSTS(2)))/(EXP( 1.00000*ALGBRC(32)) - 1.00000)
ALGBRC(130) = CONSTS(211)*ALGBRC(120)*STATES(39)
ALGBRC(132) = CONSTS(211)*ALGBRC(120)*ALGBRC(122)
ALGBRC(140) = ALGBRC(130)+ALGBRC(132)
ALGBRC(131) = CONSTS(213)*ALGBRC(120)*STATES(40)
ALGBRC(133) = CONSTS(213)*ALGBRC(120)*ALGBRC(123)
ALGBRC(141) = ALGBRC(131)+ALGBRC(133)
ALGBRC(146) = ( ALGBRC(141)*ALGBRC(64)+ ALGBRC(140)*(1.00000 - ALGBRC(64)))*CONSTS(66)
ALGBRC(211) = 1.00000/(1.00000+CONSTS(137)/STATES(3) ** 2.00000)
ALGBRC(191) = 1.00000+ (CONSTS(2)/CONSTS(128))*(1.00000+1.00000/ALGBRC(154))
ALGBRC(192) = CONSTS(2)/( CONSTS(128)*ALGBRC(154)*ALGBRC(191))
ALGBRC(195) = ALGBRC(192)*CONSTS(132)
ALGBRC(185) = 1.00000+ (STATES(5)/CONSTS(128))*(1.00000+ALGBRC(154))
ALGBRC(186) = ( STATES(5)*ALGBRC(154))/( CONSTS(128)*ALGBRC(185))
ALGBRC(198) = ALGBRC(186)*CONSTS(132)
ALGBRC(188) = 1.00000+ (STATES(5)/CONSTS(126))*(1.00000+STATES(5)/CONSTS(127))
ALGBRC(189) = ( STATES(5)*STATES(5))/( ALGBRC(188)*CONSTS(126)*CONSTS(127))
ALGBRC(201) = ALGBRC(189)*ALGBRC(186)*CONSTS(130)
ALGBRC(202) = ALGBRC(192)*CONSTS(223)*CONSTS(130)
ALGBRC(193) = 1.00000/ALGBRC(191)
ALGBRC(194) = ALGBRC(193)*CONSTS(131)
ALGBRC(196) = ALGBRC(194)+ALGBRC(195)
ALGBRC(187) = 1.00000/ALGBRC(185)
ALGBRC(197) = ( ALGBRC(187)*CONSTS(131))/ALGBRC(153)
ALGBRC(199) = ALGBRC(197)+ALGBRC(198)
ALGBRC(190) = 1.00000/ALGBRC(188)
ALGBRC(200) = ALGBRC(190)*STATES(3)*CONSTS(133)
ALGBRC(203) = CONSTS(226)*ALGBRC(199)*(ALGBRC(201)+ALGBRC(200))+ CONSTS(227)*ALGBRC(201)*(CONSTS(226)+ALGBRC(196))
ALGBRC(204) = CONSTS(225)*ALGBRC(201)*(ALGBRC(199)+CONSTS(227))+ ALGBRC(199)*ALGBRC(200)*(CONSTS(225)+ALGBRC(202))
ALGBRC(205) = CONSTS(225)*ALGBRC(196)*(ALGBRC(201)+ALGBRC(200))+ ALGBRC(202)*ALGBRC(200)*(CONSTS(226)+ALGBRC(196))
ALGBRC(206) = CONSTS(226)*ALGBRC(202)*(ALGBRC(199)+CONSTS(227))+ ALGBRC(196)*CONSTS(227)*(CONSTS(225)+ALGBRC(202))
ALGBRC(207) = ALGBRC(203)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(208) = ALGBRC(204)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(209) = ALGBRC(205)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(210) = ALGBRC(206)/(ALGBRC(203)+ALGBRC(204)+ALGBRC(205)+ALGBRC(206))
ALGBRC(212) = ( 3.00000*( ALGBRC(210)*ALGBRC(201) - ALGBRC(207)*ALGBRC(202))+ ALGBRC(209)*ALGBRC(198)) - ALGBRC(208)*ALGBRC(195)
ALGBRC(213) = ALGBRC(208)*CONSTS(226) - ALGBRC(207)*CONSTS(225)
ALGBRC(214) = 0.200000*CONSTS(221)*ALGBRC(211)*( CONSTS(8)*ALGBRC(212)+ CONSTS(9)*ALGBRC(213))
ALGBRC(119) = ( 4.00000*ALGBRC(23)*( 1.20000*STATES(3)*EXP( 2.00000*ALGBRC(32)) - 0.341000*CONSTS(3)))/(EXP( 2.00000*ALGBRC(32)) - 1.00000)
ALGBRC(124) = CONSTS(209)*ALGBRC(119)*STATES(39)
ALGBRC(127) = CONSTS(209)*ALGBRC(119)*ALGBRC(122)
ALGBRC(138) = ALGBRC(124)+ALGBRC(127)
ALGBRC(125) = CONSTS(210)*ALGBRC(119)*STATES(40)
ALGBRC(128) = CONSTS(210)*ALGBRC(119)*ALGBRC(123)
ALGBRC(139) = ALGBRC(125)+ALGBRC(128)
ALGBRC(144) = ( ALGBRC(139)*ALGBRC(64)+ ALGBRC(138)*(1.00000 - ALGBRC(64)))*CONSTS(66)
ALGBRC(240) = ( CONSTS(167)*STATES(9))/(CONSTS(168)+STATES(9))
ALGBRC(238) = ( (1.00000 - CONSTS(166))*CONSTS(165)*16.0000*ALGBRC(23)*( 1.20000*STATES(9)*EXP( 2.00000*ALGBRC(32)) - 0.341000*CONSTS(3)))/(EXP( 2.00000*ALGBRC(32)) - 1.00000)
ALGBRC(241) = ( (STATES(3) - STATES(9))*1.70000)/0.200000
ALGBRC(243) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(242) = 1.00000 - 1.00000/(1.00000+EXP((STATES(8) - 0.300000)/0.100000))
ALGBRC(244) = TERNRY(CONSTS(1).EQ.2.00000, CONSTS(192)*ALGBRC(242)*STATES(55)*STATES(56)*(STATES(8) - STATES(3)), CONSTS(171)*ALGBRC(242)*STATES(55)*STATES(56)*(STATES(8) - STATES(3)))
ALGBRC(245) = TERNRY(CONSTS(1).EQ.2.00000, CONSTS(193)*1.70000*ALGBRC(242)*STATES(55)*STATES(57)*(STATES(8) - STATES(3)), CONSTS(193)*ALGBRC(242)*STATES(55)*STATES(57)*(STATES(8) - STATES(3)))
ALGBRC(246) = (1.00000 - ALGBRC(243))*ALGBRC(244)+ ALGBRC(243)*ALGBRC(245)
ALGBRC(44) = 1.00000/(1.00000+( CONSTS(28)*CONSTS(29))/CONSTS(29)+STATES(3) ** 2.00000+( CONSTS(30)*CONSTS(31))/CONSTS(31)+STATES(3) ** 2.00000)
ALGBRC(247) = ( CONSTS(194)*0.00437500*STATES(9))/(STATES(9)+0.000920000)
ALGBRC(248) = ( CONSTS(194)*2.75000*0.00437500*STATES(9))/((STATES(9)+0.000920000) - 0.000170000)
ALGBRC(249) = 1.00000/(1.00000+CONSTS(19)/ALGBRC(40))
ALGBRC(251) = CONSTS(173)*( (1.00000 - ALGBRC(249))*ALGBRC(247)+ ALGBRC(249)*ALGBRC(248))
ALGBRC(250) = ( 0.0123000*STATES(8))/15.0000
ALGBRC(42) = 1.00000/(1.00000+( CONSTS(174)*CONSTS(25))/CONSTS(25)+STATES(9) ** 2.00000+( CONSTS(26)*CONSTS(27))/CONSTS(27)+STATES(9) ** 2.00000)
ALGBRC(46) = 1.00000/(1.00000+( CONSTS(32)*CONSTS(33))/CONSTS(33)+STATES(8) ** 2.00000)
ALGBRC(15) = ( CONSTS(195)*(STATES(9)/CONSTS(196) ** CONSTS(198) - STATES(8)/CONSTS(197) ** CONSTS(198)))/(1.00000+STATES(9)/CONSTS(196) ** CONSTS(198)+STATES(8)/CONSTS(197) ** CONSTS(198))
ALGBRC(92) = 1.00000/(ALGBRC(90)+ALGBRC(91))
ALGBRC(93) = ALGBRC(90)/(ALGBRC(90)+ALGBRC(91))
ALGBRC(98) = 1.00000/(ALGBRC(96)+ALGBRC(97))
ALGBRC(99) = ALGBRC(96)/(ALGBRC(96)+ALGBRC(97))
ALGBRC(111) = 1.00000/(ALGBRC(108)+ALGBRC(105))
ALGBRC(112) = ALGBRC(108)/(ALGBRC(108)+ALGBRC(105))
ALGBRC(117) = 1.00000/(ALGBRC(115)+ALGBRC(113))
ALGBRC(118) = ALGBRC(115)/(ALGBRC(115)+ALGBRC(113))
ALGBRC(126) = ALGBRC(124)*(1.00000 - ALGBRC(64))+ ALGBRC(125)*ALGBRC(64)
ALGBRC(129) = ALGBRC(127)*(1.00000 - ALGBRC(64))+ ALGBRC(128)*ALGBRC(64)
ALGBRC(145) = ALGBRC(144)/ALGBRC(119)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END