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 107 entries in the algebraic variable array.
C There are a total of 35 entries in each of the rate and state variable arrays.
C There are a total of 89 entries in the constant variable array.
C
C
C CONSTS(1) is I_hold in component interface (pA_per_pF).
C CONSTS(2) is I_test in component interface (pA_per_pF).
C CONSTS(3) is test_start in component interface (msec).
C CONSTS(4) is test_end in component interface (msec).
C ALGBRC(1) is Ist in component interface (pA_per_pF).
C VOI is time in component environment (msec).
C ALGBRC(105) is I_tot in component membrane_potential (pA_per_pF).
C ALGBRC(87) is I_Ca_tot in component membrane_potential (pA_per_pF).
C STATES(1) is v in component membrane_potential (mV).
C CONSTS(74) is Cm in component parameters (uF_per_cm2).
C ALGBRC(64) is ina in component I_Na (pA_per_pF).
C ALGBRC(71) is ical in component I_CaL (pA_per_pF).
C ALGBRC(73) is icat in component I_CaT (pA_per_pF).
C ALGBRC(74) is ib in component I_b (pA_per_pF).
C ALGBRC(75) is iherg in component I_HERG (pA_per_pF).
C ALGBRC(76) is ikcnq1 in component I_KCNQ1 (pA_per_pF).
C ALGBRC(77) is ikcnq4 in component I_KCNQ4 (pA_per_pF).
C ALGBRC(78) is ikcnq5 in component I_KCNQ5 (pA_per_pF).
C ALGBRC(79) is ik1 in component I_K1 (pA_per_pF).
C ALGBRC(80) is ik2 in component I_K2 (pA_per_pF).
C ALGBRC(81) is ika in component I_Ka (pA_per_pF).
C ALGBRC(82) is iBKa in component I_BKa (pA_per_pF).
C ALGBRC(83) is iBKab in component I_BKab (pA_per_pF).
C ALGBRC(84) is ih in component I_h (pA_per_pF).
C ALGBRC(85) is icl in component I_Cl (pA_per_pF).
C ALGBRC(89) is insna in component I_ns (pA_per_pF).
C ALGBRC(86) is insca in component I_ns (pA_per_pF).
C ALGBRC(91) is insk in component I_ns (pA_per_pF).
C ALGBRC(93) is inak in component I_NaK (pA_per_pF).
C ALGBRC(102) is inaca in component I_NaCa (pA_per_pF).
C ALGBRC(106) is J_tot in component Ca_Concentrations (mM_per_msec).
C ALGBRC(88) is J_Ca_mem in component Ca_Concentrations (mM_per_msec).
C STATES(2) is cai in component Ca_Concentrations (mM).
C ALGBRC(101) is jnaca in component J_NaCa (mM_per_msec).
C ALGBRC(103) is jpmca in component J_PMCA (mM_per_msec).
C CONSTS(69) is buff in component parameters (dimensionless).
C CONSTS(70) is AV in component parameters (cm2_per_uL).
C CONSTS(71) is zca in component parameters (dimensionless).
C CONSTS(73) is frdy in component parameters (coulomb_per_mole).
C ALGBRC(90) is jcamem_plot in component Ca_Concentrations (M_per_msec).
C ALGBRC(107) is jpmca_plot in component Ca_Concentrations (M_per_msec).
C ALGBRC(104) is jnaca_plot in component Ca_Concentrations (M_per_msec).
C CONSTS(5) is conversion in component Ca_Concentrations (mM_to_M).
C CONSTS(12) is ki in component parameters (mM).
C CONSTS(60) is nai in component parameters (mM).
C CONSTS(63) is cli in component parameters (mM).
C CONSTS(64) is ko in component parameters (mM).
C CONSTS(13) is cao in component parameters (mM).
C CONSTS(65) is nao in component parameters (mM).
C CONSTS(67) is clo in component parameters (mM).
C CONSTS(14) is mgo in component parameters (mM).
C CONSTS(15) is zna in component parameters (dimensionless).
C CONSTS(16) is zk in component parameters (dimensionless).
C CONSTS(72) is R in component parameters (joule_per_kelvin_per_kilomole).
C CONSTS(75) is temp in component parameters (kelvin).
C CONSTS(17) is gna in component parameters (nS_per_pF).
C CONSTS(18) is gcal in component parameters (nS_per_pF).
C CONSTS(19) is ecal in component parameters (mV).
C CONSTS(20) is kmca in component parameters (mM).
C CONSTS(21) is gcat in component parameters (nS_per_pF).
C CONSTS(22) is ecat in component parameters (mV).
C CONSTS(23) is gkca in component parameters (nS_per_pF).
C CONSTS(24) is gb in component parameters (nS_per_pF).
C CONSTS(25) is gk1 in component parameters (nS_per_pF).
C CONSTS(26) is gk2 in component parameters (nS_per_pF).
C CONSTS(27) is gbka in component parameters (dimensionless).
C CONSTS(28) is gbkab in component parameters (dimensionless).
C CONSTS(29) is gka in component parameters (nS_per_pF).
C CONSTS(30) is gkq1 in component parameters (nS_per_pF).
C CONSTS(31) is gkq4 in component parameters (nS_per_pF).
C CONSTS(32) is gkq5 in component parameters (nS_per_pF).
C CONSTS(33) is gherg in component parameters (nS_per_pF).
C CONSTS(34) is gcl in component parameters (nS_per_pF).
C CONSTS(35) is gh in component parameters (nS_per_pF).
C CONSTS(36) is gns in component parameters (nS_per_pF).
C CONSTS(37) is PnsK in component parameters (dimensionless).
C CONSTS(38) is PnsNa in component parameters (dimensionless).
C CONSTS(39) is PnsCa in component parameters (dimensionless).
C CONSTS(40) is PnsCs in component parameters (dimensionless).
C CONSTS(41) is gnsCa in component parameters (dimensionless).
C CONSTS(42) is gnsNa in component parameters (dimensionless).
C CONSTS(43) is gnsK in component parameters (dimensionless).
C CONSTS(44) is gnsCs in component parameters (dimensionless).
C CONSTS(76) is ginak in component parameters (pA_per_pF).
C CONSTS(81) is nakKmko in component parameters (mM).
C CONSTS(83) is nakKmnai in component parameters (mM).
C CONSTS(77) is PK in component parameters (dimensionless).
C CONSTS(82) is PNa in component parameters (dimensionless).
C CONSTS(45) is Jpmca in component parameters (mM_per_msec).
C CONSTS(46) is Kmpmca in component parameters (mM).
C CONSTS(47) is npmca in component parameters (dimensionless).
C CONSTS(85) is Jnaca in component parameters (mM_per_msec).
C CONSTS(48) is Kmallo in component parameters (mM).
C CONSTS(49) is nallo in component parameters (dimensionless).
C CONSTS(86) is ksat in component parameters (dimensionless).
C CONSTS(87) is xgamma in component parameters (dimensionless).
C CONSTS(50) is Kmnai in component parameters (mM).
C CONSTS(51) is Kmcai in component parameters (mM).
C CONSTS(52) is Kmnao in component parameters (mM).
C CONSTS(53) is Kmcao in component parameters (mM).
C CONSTS(54) is Fmax in component parameters (uN).
C CONSTS(55) is FKm in component parameters (nM).
C CONSTS(56) is Fn in component parameters (dimensionless).
C ALGBRC(29) is vFRT in component parameters (dimensionless).
C CONSTS(78) is ena in component parameters (mV).
C CONSTS(79) is ek in component parameters (mV).
C CONSTS(84) is eh in component parameters (mV).
C CONSTS(80) is ecl in component parameters (mV).
C ALGBRC(59) is enscc in component parameters (mV).
C ALGBRC(2) is wss in component Ca_dependent_Force (dimensionless).
C ALGBRC(30) is wtc in component Ca_dependent_Force (msec).
C CONSTS(6) is conversion in component Ca_dependent_Force (nM_to_mM).
C ALGBRC(3) is Force in component Ca_dependent_Force (uN).
C STATES(3) is w in component Ca_dependent_Force (dimensionless).
C ALGBRC(4) is mss in component I_Na (dimensionless).
C ALGBRC(5) is hss in component I_Na (dimensionless).
C ALGBRC(31) is mtc in component I_Na (msec).
C ALGBRC(32) is htc in component I_Na (msec).
C STATES(4) is m in component I_Na (dimensionless).
C STATES(5) is h in component I_Na (dimensionless).
C ALGBRC(6) is dss in component I_CaL (dimensionless).
C ALGBRC(7) is fss in component I_CaL (dimensionless).
C ALGBRC(69) is fca in component I_CaL (dimensionless).
C ALGBRC(33) is dtc in component I_CaL (msec).
C CONSTS(57) is f1tc in component I_CaL (msec).
C ALGBRC(34) is f2tc in component I_CaL (msec).
C STATES(6) is d in component I_CaL (dimensionless).
C STATES(7) is f1 in component I_CaL (dimensionless).
C STATES(8) is f2 in component I_CaL (dimensionless).
C ALGBRC(8) is bss in component I_CaT (dimensionless).
C ALGBRC(9) is gss in component I_CaT (dimensionless).
C ALGBRC(35) is btc in component I_CaT (msec).
C ALGBRC(36) is gtc in component I_CaT (msec).
C STATES(9) is b in component I_CaT (dimensionless).
C STATES(10) is g in component I_CaT (dimensionless).
C ALGBRC(10) is hnss in component I_HERG (dimensionless).
C ALGBRC(11) is hsss in component I_HERG (dimensionless).
C ALGBRC(37) is hn1tc in component I_HERG (msec).
C ALGBRC(38) is hn2tc in component I_HERG (msec).
C ALGBRC(39) is hstc in component I_HERG (msec).
C STATES(11) is hn1 in component I_HERG (dimensionless).
C STATES(12) is hn2 in component I_HERG (dimensionless).
C STATES(13) is hs in component I_HERG (dimensionless).
C ALGBRC(12) is nq1ss in component I_KCNQ1 (dimensionless).
C ALGBRC(13) is wq1ss in component I_KCNQ1 (dimensionless).
C ALGBRC(14) is sq1ss in component I_KCNQ1 (dimensionless).
C ALGBRC(40) is nq1ftc in component I_KCNQ1 (msec).
C ALGBRC(41) is nq1stc in component I_KCNQ1 (msec).
C ALGBRC(42) is wq1tc in component I_KCNQ1 (msec).
C CONSTS(58) is sq1tc in component I_KCNQ1 (msec).
C STATES(14) is nq1f in component I_KCNQ1 (dimensionless).
C STATES(15) is nq1s in component I_KCNQ1 (dimensionless).
C STATES(16) is wq1 in component I_KCNQ1 (dimensionless).
C STATES(17) is sq1 in component I_KCNQ1 (dimensionless).
C ALGBRC(15) is nq4ss in component I_KCNQ4 (dimensionless).
C ALGBRC(16) is sq4ss in component I_KCNQ4 (dimensionless).
C ALGBRC(43) is nq4tc in component I_KCNQ4 (msec).
C ALGBRC(44) is sq4tc in component I_KCNQ4 (msec).
C STATES(18) is nq4 in component I_KCNQ4 (dimensionless).
C STATES(19) is sq4 in component I_KCNQ4 (dimensionless).
C ALGBRC(17) is nq5ss in component I_KCNQ5 (dimensionless).
C ALGBRC(18) is wq5ss in component I_KCNQ5 (dimensionless).
C ALGBRC(19) is sq5ss in component I_KCNQ5 (dimensionless).
C ALGBRC(45) is nq5ftc in component I_KCNQ5 (msec).
C CONSTS(59) is nq5stc in component I_KCNQ5 (msec).
C ALGBRC(46) is wq5tc in component I_KCNQ5 (msec).
C ALGBRC(47) is sq5tc in component I_KCNQ5 (msec).
C STATES(20) is nq5f in component I_KCNQ5 (dimensionless).
C STATES(21) is nq5s in component I_KCNQ5 (dimensionless).
C STATES(22) is wq5 in component I_KCNQ5 (dimensionless).
C STATES(23) is sq5 in component I_KCNQ5 (dimensionless).
C ALGBRC(20) is qss in component I_K1 (dimensionless).
C ALGBRC(21) is rss in component I_K1 (dimensionless).
C ALGBRC(48) is qtc in component I_K1 (msec).
C ALGBRC(49) is r1tc in component I_K1 (msec).
C ALGBRC(50) is r2tc in component I_K1 (msec).
C STATES(24) is q in component I_K1 (dimensionless).
C STATES(25) is r1 in component I_K1 (dimensionless).
C STATES(26) is r2 in component I_K1 (dimensionless).
C ALGBRC(22) is pss in component I_K2 (dimensionless).
C ALGBRC(23) is kss in component I_K2 (dimensionless).
C ALGBRC(51) is ptc in component I_K2 (msec).
C ALGBRC(52) is k1tc in component I_K2 (msec).
C ALGBRC(53) is k2tc in component I_K2 (msec).
C STATES(27) is p in component I_K2 (dimensionless).
C STATES(28) is k1 in component I_K2 (dimensionless).
C STATES(29) is k2 in component I_K2 (dimensionless).
C ALGBRC(24) is sss in component I_Ka (dimensionless).
C ALGBRC(25) is xss in component I_Ka (dimensionless).
C ALGBRC(54) is stc in component I_Ka (msec).
C ALGBRC(55) is xtc in component I_Ka (msec).
C STATES(30) is s in component I_Ka (dimensionless).
C STATES(31) is x in component I_Ka (dimensionless).
C ALGBRC(26) is xass_z in component I_BKa (dimensionless).
C ALGBRC(56) is xass_vh in component I_BKa (mV).
C CONSTS(7) is conversion in component I_BKa (mM_to_M).
C ALGBRC(60) is xass in component I_BKa (dimensionless).
C ALGBRC(65) is xatc in component I_BKa (msec).
C STATES(32) is xa in component I_BKa (dimensionless).
C ALGBRC(27) is xabss_z in component I_BKab (dimensionless).
C ALGBRC(57) is xabss_vh in component I_BKab (mV).
C CONSTS(8) is conversion in component I_BKab (mM_to_M).
C ALGBRC(61) is xabss in component I_BKab (dimensionless).
C ALGBRC(66) is xabtc in component I_BKab (msec).
C STATES(33) is xab in component I_BKab (dimensionless).
C ALGBRC(28) is yss in component I_h (dimensionless).
C ALGBRC(67) is ytc in component I_h (msec).
C ALGBRC(58) is ya in component I_h (per_msec).
C ALGBRC(62) is yb in component I_h (per_msec).
C STATES(34) is y in component I_h (dimensionless).
C ALGBRC(70) is css in component I_Cl (dimensionless).
C ALGBRC(72) is ctc in component I_Cl (msec).
C ALGBRC(63) is K1cl in component I_Cl (mM).
C ALGBRC(68) is K2cl in component I_Cl (dimensionless).
C STATES(35) is c in component I_Cl (dimensionless).
C CONSTS(61) is fmg in component I_ns (dimensionless).
C CONSTS(68) is gs_nao in component I_ns (dimensionless).
C CONSTS(62) is gs_cao in component I_ns (dimensionless).
C CONSTS(66) is gs_ko in component I_ns (dimensionless).
C CONSTS(9) is tinyamount in component I_ns (mM).
C ALGBRC(92) is fnak in component I_NaK (dimensionless).
C CONSTS(88) is knak in component I_NaK (dimensionless).
C CONSTS(89) is nnak in component I_NaK (dimensionless).
C CONSTS(10) is inaca_sign in component I_NaCa (dimensionless).
C ALGBRC(94) is f1naca in component J_NaCa (dimensionless).
C ALGBRC(95) is f2naca in component J_NaCa (dimensionless).
C ALGBRC(96) is fallo in component J_NaCa (dimensionless).
C ALGBRC(97) is naca_Eup in component J_NaCa (mM4).
C ALGBRC(98) is naca_Ed1 in component J_NaCa (dimensionless).
C ALGBRC(99) is naca_Ed2 in component J_NaCa (mM4).
C ALGBRC(100) is naca_Ed3 in component J_NaCa (mM4).
C CONSTS(11) is jnaca_sign in component J_NaCa (dimensionless).
C RATES(1) is d/dt v in component membrane_potential (mV).
C RATES(2) is d/dt cai in component Ca_Concentrations (mM).
C RATES(3) is d/dt w in component Ca_dependent_Force (dimensionless).
C RATES(4) is d/dt m in component I_Na (dimensionless).
C RATES(5) is d/dt h in component I_Na (dimensionless).
C RATES(6) is d/dt d in component I_CaL (dimensionless).
C RATES(7) is d/dt f1 in component I_CaL (dimensionless).
C RATES(8) is d/dt f2 in component I_CaL (dimensionless).
C RATES(9) is d/dt b in component I_CaT (dimensionless).
C RATES(10) is d/dt g in component I_CaT (dimensionless).
C RATES(11) is d/dt hn1 in component I_HERG (dimensionless).
C RATES(12) is d/dt hn2 in component I_HERG (dimensionless).
C RATES(13) is d/dt hs in component I_HERG (dimensionless).
C RATES(14) is d/dt nq1f in component I_KCNQ1 (dimensionless).
C RATES(15) is d/dt nq1s in component I_KCNQ1 (dimensionless).
C RATES(16) is d/dt wq1 in component I_KCNQ1 (dimensionless).
C RATES(17) is d/dt sq1 in component I_KCNQ1 (dimensionless).
C RATES(18) is d/dt nq4 in component I_KCNQ4 (dimensionless).
C RATES(19) is d/dt sq4 in component I_KCNQ4 (dimensionless).
C RATES(20) is d/dt nq5f in component I_KCNQ5 (dimensionless).
C RATES(21) is d/dt nq5s in component I_KCNQ5 (dimensionless).
C RATES(22) is d/dt wq5 in component I_KCNQ5 (dimensionless).
C RATES(23) is d/dt sq5 in component I_KCNQ5 (dimensionless).
C RATES(24) is d/dt q in component I_K1 (dimensionless).
C RATES(25) is d/dt r1 in component I_K1 (dimensionless).
C RATES(26) is d/dt r2 in component I_K1 (dimensionless).
C RATES(27) is d/dt p in component I_K2 (dimensionless).
C RATES(28) is d/dt k1 in component I_K2 (dimensionless).
C RATES(29) is d/dt k2 in component I_K2 (dimensionless).
C RATES(30) is d/dt s in component I_Ka (dimensionless).
C RATES(31) is d/dt x in component I_Ka (dimensionless).
C RATES(32) is d/dt xa in component I_BKa (dimensionless).
C RATES(33) is d/dt xab in component I_BKab (dimensionless).
C RATES(34) is d/dt y in component I_h (dimensionless).
C RATES(35) is d/dt c in component I_Cl (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0
CONSTS(2) = -0.25
CONSTS(3) = 1000
CONSTS(4) = 50000
STATES(1) = -50.80774403486136
STATES(2) = 0.0001235555354004079
CONSTS(5) = 1000
CONSTS(6) = 1e-6
STATES(3) = 0.2768302028689621
STATES(4) = 0.166998688814229
STATES(5) = 0.3156075507278521
STATES(6) = 0.01605749091924106
STATES(7) = 0.861723329545759
STATES(8) = 0.8617233295451014
STATES(9) = 0.5857399935992883
STATES(10) = 0.02817518341064734
STATES(11) = 0.02498997383730429
STATES(12) = 0.02498997383730429
STATES(13) = 0.5292140214748325
STATES(14) = 0.09043683263784785
STATES(15) = 0.09043683263784785
STATES(16) = 0.9230513956601629
STATES(17) = 0.7426740827665872
STATES(18) = 0.1081102112425531
STATES(19) = 0.6280622663818354
STATES(20) = 0.2618890409031491
STATES(21) = 0.2618890409031491
STATES(22) = 0.9230513956601629
STATES(23) = 0.6280622663818354
STATES(24) = 0.22560249352574
STATES(25) = 0.1261674553968813
STATES(26) = 0.1261674420211442
STATES(27) = 0.1358747732875166
STATES(28) = 0.9944827384537837
STATES(29) = 0.9944827384537837
STATES(30) = 0.04562272513582834
STATES(31) = 0.06162789823439722
CONSTS(7) = 1000.0
STATES(32) = 0.0003575122973592095
CONSTS(8) = 1000.0
STATES(33) = 0.002673927875795617
STATES(34) = 0.001821587846781853
STATES(35) = 0.0006695090454068198
CONSTS(9) = 1e-8
CONSTS(10) = -1
CONSTS(11) = -1
CONSTS(12) = 140.000
CONSTS(13) = 2.50000
CONSTS(14) = 0.500000
CONSTS(15) = 1.00000
CONSTS(16) = 1.00000
CONSTS(17) = 0.00000
CONSTS(18) = 0.600000
CONSTS(19) = 45.0000
CONSTS(20) = 0.00100000
CONSTS(21) = 0.0580000
CONSTS(22) = 42.0000
CONSTS(23) = 0.800000
CONSTS(24) = 0.00000
CONSTS(25) = 0.240000
CONSTS(26) = 0.0320000
CONSTS(27) = 0.200000
CONSTS(28) = 0.100000
CONSTS(29) = 0.160000
CONSTS(30) = 0.00320000
CONSTS(31) = 0.0240000
CONSTS(32) = 0.0160000
CONSTS(33) = 0.0800000
CONSTS(34) = 0.187500
CONSTS(35) = 0.0542000
CONSTS(36) = 0.0123000
CONSTS(37) = 1.30000
CONSTS(38) = 0.900000
CONSTS(39) = 0.890000
CONSTS(40) = 1.00000
CONSTS(41) = 0.500000
CONSTS(42) = 1.00000
CONSTS(43) = 1.19000
CONSTS(44) = 1.60000
CONSTS(45) = 3.50000e-07
CONSTS(46) = 0.000500000
CONSTS(47) = 2.00000
CONSTS(48) = 0.000300000
CONSTS(49) = 4.00000
CONSTS(50) = 30.0000
CONSTS(51) = 0.00700000
CONSTS(52) = 87.5000
CONSTS(53) = 1.30000
CONSTS(54) = 3.00000
CONSTS(55) = 161.301
CONSTS(56) = 3.60205
CONSTS(57) = 12.0000
CONSTS(58) = 50000.0
CONSTS(59) = 1000.00
CONSTS(60) = 4.00000
CONSTS(61) = 0.108043+0.903902/(1.00000+CONSTS(14)/0.281007 ** 1.29834)
CONSTS(62) = ( (1.00000/0.000525000)*0.0300000)/(1.00000+150.000/(CONSTS(13)+CONSTS(9)) ** 2.00000)
CONSTS(63) = 46.0000
CONSTS(64) = 6.00000
CONSTS(65) = 130.000
CONSTS(66) = ( (1.00000/0.0123000)*0.0300000)/(1.00000+150.000/(CONSTS(64)+CONSTS(9)) ** 2.00000)
CONSTS(67) = 130.000
CONSTS(68) = ( (1.00000/0.0123000)*0.0300000)/(1.00000+150.000/(CONSTS(65)+CONSTS(9)) ** 2.00000)
CONSTS(69) = 0.0150000
CONSTS(70) = 4.00000
CONSTS(71) = 2.00000
CONSTS(72) = 8314.00
CONSTS(73) = 96485.0
CONSTS(74) = 1.00000
CONSTS(75) = 308.000
CONSTS(76) = 1.70000
CONSTS(77) = 1.00000
CONSTS(78) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log(CONSTS(65)/CONSTS(60))
CONSTS(79) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log(CONSTS(64)/CONSTS(12))
CONSTS(80) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log(CONSTS(63)/CONSTS(67))
CONSTS(81) = 2.00000
CONSTS(82) = 0.350000
CONSTS(83) = 22.0000
CONSTS(84) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log((CONSTS(64)+ (CONSTS(82)/CONSTS(77))*CONSTS(65))/(CONSTS(12)+ (CONSTS(82)/CONSTS(77))*CONSTS(60)))
CONSTS(85) = 3.50000e-06
CONSTS(86) = 0.270000
CONSTS(87) = 0.350000
CONSTS(88) = 1.00000/(1.00000+CONSTS(81)/CONSTS(64) ** 1.50000)
CONSTS(89) = 1.00000/(1.00000+CONSTS(83)/CONSTS(60) ** 2.00000)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+38.0000)/7.00000))
RATES(7) = (ALGBRC(7) - STATES(7))/CONSTS(57)
ALGBRC(14) = 0.340000+0.660000/(1.00000+EXP((STATES(1)+45.3000)/12.3000))
RATES(17) = (ALGBRC(14) - STATES(17))/CONSTS(58)
ALGBRC(17) = 1.00000/(1.00000+EXP(- (STATES(1)+36.5500)/13.7600))
RATES(21) = (ALGBRC(17) - STATES(21))/CONSTS(59)
ALGBRC(2) = 1.00000/(1.00000+( CONSTS(55)*CONSTS(6))/STATES(2) ** CONSTS(56))
ALGBRC(30) = 4000.00*(0.234845+(1.00000 - 0.234845)/(1.00000+STATES(2)/( CONSTS(55)*CONSTS(6)) ** CONSTS(56)))
RATES(3) = (ALGBRC(2) - STATES(3))/ALGBRC(30)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+35.9584)/9.24013))
ALGBRC(31) = 0.250000+7.00000/(1.00000+EXP((STATES(1)+38.0000)/10.0000))
RATES(4) = (ALGBRC(4) - STATES(4))/ALGBRC(31)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+57.0000)/8.00000))
ALGBRC(32) = 0.900000+1002.85/(1.00000+(STATES(1)+47.5000)/1.50000 ** 2.00000)
RATES(5) = (ALGBRC(5) - STATES(5))/ALGBRC(32)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+22.0000)/7.00000))
ALGBRC(33) = 2.29000+5.70000/(1.00000+(STATES(1)+29.9700)/9.00000 ** 2.00000)
RATES(6) = (ALGBRC(6) - STATES(6))/ALGBRC(33)
ALGBRC(34) = 90.9699*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1)+13.9629)/45.3782))*(1.00000+EXP(- (STATES(1)+9.49866)/3.39450))))
RATES(8) = (ALGBRC(7) - STATES(8))/ALGBRC(34)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+54.2300)/9.88000))
ALGBRC(35) = 0.450000+3.90000/(1.00000+(STATES(1)+66.0000)/26.0000 ** 2.00000)
RATES(9) = (ALGBRC(8) - STATES(9))/ALGBRC(35)
ALGBRC(9) = 0.0200000+0.980000/(1.00000+EXP((STATES(1)+72.9780)/4.64000))
ALGBRC(36) = 150.000*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 417.430)/203.180))*(1.00000+EXP(- (STATES(1)+61.1100)/8.07000))))
RATES(10) = (ALGBRC(9) - STATES(10))/ALGBRC(36)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+16.0000)/9.50000))
ALGBRC(37) = 46.0999+1685.76/( (1.00000+EXP(- (STATES(1)+40.8489)/13.7802))*(1.00000+EXP((STATES(1)+20.6372)/15.1113)))
RATES(11) = (ALGBRC(10) - STATES(11))/ALGBRC(37)
ALGBRC(38) = 475.667+16321.6/( (1.00000+EXP(- (STATES(1)+41.8328)/6.96673))*(1.00000+EXP((STATES(1)+23.2432)/21.2949)))
RATES(12) = (ALGBRC(10) - STATES(12))/ALGBRC(38)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+48.0000)/24.0000))
ALGBRC(39) = 19.7864/(1.00000+(STATES(1)+20.7136)/44.2868 ** 2.00000) - 0.378843
RATES(13) = (ALGBRC(11) - STATES(13))/ALGBRC(39)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+22.0000)/12.4800))
ALGBRC(40) = 395.300/(1.00000+(STATES(1)+38.1000)/33.5900 ** 2.00000)
RATES(14) = (ALGBRC(12) - STATES(14))/ALGBRC(40)
ALGBRC(41) = (5503.00+5345.40/(1.00000+10.0000 ** (- 23.9000 - STATES(1))*- 0.0282700)) - 4590.60/(1.00000+10.0000 ** (STATES(1)+14.1500)*- 0.0357000)
RATES(15) = (ALGBRC(12) - STATES(15))/ALGBRC(41)
ALGBRC(13) = 0.490000+0.510000/(1.00000+EXP((STATES(1)+1.08400)/28.7800))
ALGBRC(42) = 5.44000+29.2000/(1.00000+(STATES(1)+48.0900)/48.8300 ** 2.00000)
RATES(16) = (ALGBRC(13) - STATES(16))/ALGBRC(42)
ALGBRC(15) = 1.00000/(1.00000+EXP(- (STATES(1)+15.0400)/16.9500))
ALGBRC(43) = 10.0000+895.900/(1.00000+EXP((- 18.0100 - STATES(1))/31.0400))
RATES(18) = (ALGBRC(15) - STATES(18))/ALGBRC(43)
ALGBRC(16) = 0.405800/(1.00000+EXP((STATES(1)+86.8400)/15.0500))+0.594200/(1.00000+EXP((STATES(1) - 70.1300)/13.3700))
ALGBRC(44) = 1077.00+185845./(1.00000+(STATES(1) - 39.4400)/7.34400 ** 2.00000)
RATES(19) = (ALGBRC(16) - STATES(19))/ALGBRC(44)
ALGBRC(45) = 37.5100+539.000/(1.00000+(STATES(1)+40.2400)/17.7200 ** 2.00000)
RATES(20) = (ALGBRC(17) - STATES(20))/ALGBRC(45)
ALGBRC(18) = 0.490000+0.510000/(1.00000+EXP((STATES(1)+1.08400)/28.7800))
ALGBRC(46) = 5.44000+29.2000/(1.00000+(STATES(1)+48.0900)/48.8300 ** 2.00000)
RATES(22) = (ALGBRC(18) - STATES(22))/ALGBRC(46)
ALGBRC(19) = 0.405800/(1.00000+EXP((STATES(1)+86.8400)/15.0500))+0.594200/(1.00000+EXP((STATES(1) - 70.1300)/13.3700))
ALGBRC(47) = 1077.00+185845./(1.00000+(STATES(1) - 39.4400)/7.34400 ** 2.00000)
RATES(23) = (ALGBRC(19) - STATES(23))/ALGBRC(47)
ALGBRC(20) = 0.978613/(1.00000+EXP(- (STATES(1)+18.6736)/26.6606))
ALGBRC(48) = 500.000/(1.00000+(STATES(1)+60.7100)/15.7900 ** 2.00000)
RATES(24) = (ALGBRC(20) - STATES(24))/ALGBRC(48)
ALGBRC(21) = 1.00000/(1.00000+EXP((STATES(1)+63.0000)/6.30000))
ALGBRC(49) = 5000.00/(1.00000+(STATES(1)+62.7133)/35.8611 ** 2.00000)
RATES(25) = (ALGBRC(21) - STATES(25))/ALGBRC(49)
ALGBRC(50) = 30000.0+220000./(1.00000+EXP((STATES(1)+22.0000)/4.00000))
RATES(26) = (ALGBRC(21) - STATES(26))/ALGBRC(50)
ALGBRC(22) = 0.948000/(1.00000+EXP(- (STATES(1)+17.9100)/18.4000))
ALGBRC(51) = 100.000/(1.00000+(STATES(1)+64.1000)/28.6700 ** 2.00000)
RATES(27) = (ALGBRC(22) - STATES(27))/ALGBRC(51)
ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(1)+21.2000)/5.70000))
ALGBRC(52) = 1.00000e+06*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 315.000)/50.0000))*(1.00000+EXP(- (STATES(1)+74.9000)/8.00000))))
RATES(28) = (ALGBRC(23) - STATES(28))/ALGBRC(52)
ALGBRC(53) = 2.50000e+06*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 132.868)/25.3992))*(1.00000+EXP(- (STATES(1)+24.9203)/2.67915))))
RATES(29) = (ALGBRC(23) - STATES(29))/ALGBRC(53)
ALGBRC(24) = 1.00000/(1.00000+EXP(- (STATES(1)+27.7900)/7.57000))
ALGBRC(54) = 17.0000/(1.00000+(STATES(1)+20.5232)/35.0000 ** 2.00000)
RATES(30) = (ALGBRC(24) - STATES(30))/ALGBRC(54)
ALGBRC(25) = 0.0200000+0.980000/(1.00000+EXP((STATES(1)+69.5000)/6.00000))
ALGBRC(55) = 7.50000+10.0000/(1.00000+(STATES(1)+34.1765)/120.000 ** 2.00000)
RATES(31) = (ALGBRC(25) - STATES(31))/ALGBRC(55)
ALGBRC(26) = - 0.749234/(1.00000+( STATES(2)*CONSTS(7) - 0.0630535)/0.161942 ** 2.00000)+8.38384/(1.00000+( STATES(2)*CONSTS(7)+1538.29)/739.057 ** 2.00000)
ALGBRC(56) = 5011.47/(1.00000+( STATES(2)*CONSTS(7)+0.237503)/0.000239278 ** 0.422910) - 37.5137
ALGBRC(60) = 1.00000/(1.00000+EXP(( - ALGBRC(26)*CONSTS(73)*(STATES(1) - ALGBRC(56)))/( CONSTS(72)*CONSTS(75))))
ALGBRC(65) = 2.40914/(1.00000+(STATES(1) - 158.779)/- 52.1497 ** 2.00000)
RATES(32) = (ALGBRC(60) - STATES(32))/ALGBRC(65)
ALGBRC(27) = - 0.681249/(1.00000+( STATES(2)*CONSTS(8) - 0.218988)/0.428335 ** 2.00000)+1.40001/(1.00000+( STATES(2)*CONSTS(8)+228.710)/684.946 ** 2.00000)
ALGBRC(57) = 8540.23/(1.00000+( STATES(2)*CONSTS(8)+0.401189)/0.00399115 ** 0.668054) - 109.275
ALGBRC(61) = 1.00000/(1.00000+EXP(( - ALGBRC(27)*CONSTS(73)*(STATES(1) - ALGBRC(57)))/( CONSTS(72)*CONSTS(75))))
ALGBRC(66) = 13.8049/(1.00000+(STATES(1) - 153.019)/66.4952 ** 2.00000)
RATES(33) = (ALGBRC(61) - STATES(33))/ALGBRC(66)
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(1)+105.390)/8.65530))
ALGBRC(58) = 3.50000e-06*EXP( - 0.0497000*STATES(1))
ALGBRC(62) = 0.0400300*EXP( 0.0521100*STATES(1))
ALGBRC(67) = 1.00000/(ALGBRC(58)+ALGBRC(62))
RATES(34) = (ALGBRC(28) - STATES(34))/ALGBRC(67)
ALGBRC(29) = ( STATES(1)*CONSTS(73))/( CONSTS(72)*CONSTS(75))
ALGBRC(63) = 0.000600000*EXP( 2.53000*ALGBRC(29))
ALGBRC(68) = 0.100000*EXP( - 5.00000*ALGBRC(29))
ALGBRC(70) = 1.00000/(1.00000+ ALGBRC(68)*(ALGBRC(63)/STATES(2) ** 2.00000+ALGBRC(63)/STATES(2)+1.00000))
ALGBRC(72) = - 160.000+210.000/(1.00000+EXP((STATES(1)+4.56000)/11.6200))+170.000/(1.00000+EXP(- (STATES(1)+25.5000)/11.6200))
RATES(35) = (ALGBRC(70) - STATES(35))/ALGBRC(72)
ALGBRC(1) = TERNRY(VOI.GT.CONSTS(3).AND.VOI.LT.CONSTS(4), CONSTS(2), CONSTS(1))
ALGBRC(64) = CONSTS(17)*STATES(4)*STATES(4)*STATES(4)*STATES(5)*(STATES(1) - CONSTS(78))
ALGBRC(69) = 1.00000/(1.00000+STATES(2)/CONSTS(20) ** 4.00000)
ALGBRC(71) = CONSTS(18)*ALGBRC(69)*STATES(6)*STATES(6)*( 0.800000*STATES(7)+ 0.200000*STATES(8))*(STATES(1) - CONSTS(19))
ALGBRC(73) = CONSTS(21)*STATES(9)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(22))
ALGBRC(74) = CONSTS(24)*(STATES(1) - CONSTS(79))
ALGBRC(75) = CONSTS(33)*( 0.800000*STATES(11)+ 0.200000*STATES(12))*STATES(13)*(STATES(1) - CONSTS(79))
ALGBRC(76) = CONSTS(30)*( 0.300000*STATES(14)+ 0.700000*STATES(15))*STATES(16)*STATES(17)*(STATES(1) - CONSTS(79))
ALGBRC(77) = CONSTS(31)*STATES(18)*STATES(19)*(STATES(1) - CONSTS(79))
ALGBRC(78) = CONSTS(32)*( 0.200000*STATES(20)+ 0.800000*STATES(21))*STATES(22)*STATES(23)*(STATES(1) - CONSTS(79))
ALGBRC(79) = CONSTS(25)*STATES(24)*STATES(24)*( 0.380000*STATES(25)+ 0.630000*STATES(26))*(STATES(1) - CONSTS(79))
ALGBRC(80) = CONSTS(26)*STATES(27)*STATES(27)*( 0.750000*STATES(28)+ 0.250000*STATES(29))*(STATES(1) - CONSTS(79))
ALGBRC(81) = CONSTS(29)*STATES(30)*STATES(31)*(STATES(1) - CONSTS(79))
ALGBRC(82) = CONSTS(23)*CONSTS(27)*STATES(32)*(STATES(1) - CONSTS(79))
ALGBRC(83) = CONSTS(23)*CONSTS(28)*STATES(33)*(STATES(1) - CONSTS(79))
ALGBRC(84) = CONSTS(35)*STATES(34)*(STATES(1) - CONSTS(84))
ALGBRC(85) = CONSTS(34)*STATES(35)*(STATES(1) - CONSTS(80))
ALGBRC(59) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log(( CONSTS(37)*CONSTS(64)+ CONSTS(38)*CONSTS(65)+( 4.00000*CONSTS(39)*CONSTS(13))/(1.00000+EXP(ALGBRC(29))))/( CONSTS(37)*CONSTS(12)+ CONSTS(38)*CONSTS(60)+( 4.00000*CONSTS(39)*STATES(2))/(1.00000+EXP(ALGBRC(29)))))
ALGBRC(89) = CONSTS(61)*CONSTS(68)*CONSTS(42)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(86) = CONSTS(61)*CONSTS(62)*CONSTS(41)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(91) = CONSTS(61)*CONSTS(66)*CONSTS(43)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(92) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(29))+ 0.00219000*EXP(CONSTS(65)/49.7100)*EXP( - 1.90000*ALGBRC(29)))
ALGBRC(93) = CONSTS(76)*CONSTS(88)*CONSTS(89)*ALGBRC(92)
ALGBRC(96) = 1.00000/(1.00000+CONSTS(48)/STATES(2) ** CONSTS(49))
ALGBRC(94) = EXP( (CONSTS(87) - 1.00000)*ALGBRC(29))
ALGBRC(95) = EXP( CONSTS(87)*ALGBRC(29))
ALGBRC(97) = CONSTS(60) ** 3.00000*CONSTS(13)*ALGBRC(95) - CONSTS(65) ** 3.00000*STATES(2)*ALGBRC(94)
ALGBRC(98) = 1.00000+ CONSTS(86)*ALGBRC(94)
ALGBRC(99) = CONSTS(53)*CONSTS(60) ** 3.00000+ CONSTS(52) ** 3.00000*STATES(2)+ CONSTS(50) ** 3.00000*CONSTS(13)*(1.00000+STATES(2)/CONSTS(51))
ALGBRC(100) = CONSTS(13)*CONSTS(60) ** 3.00000+ CONSTS(65) ** 3.00000*STATES(2)+ CONSTS(65) ** 3.00000*CONSTS(51)*(1.00000+CONSTS(60)/CONSTS(50) ** 3.00000)
ALGBRC(101) = ( CONSTS(11)*CONSTS(85)*ALGBRC(96)*ALGBRC(97))/( ALGBRC(98)*(ALGBRC(99)+ALGBRC(100)))
ALGBRC(102) = (( 0.500000*CONSTS(71)*CONSTS(73))/( CONSTS(70)*CONSTS(74)*CONSTS(69)))*CONSTS(10)*ALGBRC(101)
ALGBRC(105) = ALGBRC(64)+ALGBRC(84)+ALGBRC(102)+ALGBRC(93)+ALGBRC(71)+ALGBRC(73)+ALGBRC(85)+ALGBRC(75)+ALGBRC(76)+ALGBRC(77)+ALGBRC(78)+ALGBRC(79)+ALGBRC(80)+ALGBRC(81)+ALGBRC(82)+ALGBRC(83)+ALGBRC(89)+ALGBRC(91)+ALGBRC(86)+ALGBRC(74)
RATES(1) = - (ALGBRC(105)+ALGBRC(1))
ALGBRC(87) = ALGBRC(71)+ALGBRC(73)+ALGBRC(86)
ALGBRC(88) = (( CONSTS(70)*CONSTS(74)*CONSTS(69))/( CONSTS(71)*CONSTS(73)))*ALGBRC(87)
ALGBRC(103) = CONSTS(45)/(1.00000+CONSTS(46)/STATES(2) ** CONSTS(47))
ALGBRC(106) = ALGBRC(88)+ALGBRC(101)+ALGBRC(103)
RATES(2) = - ALGBRC(106)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+38.0000)/7.00000))
ALGBRC(14) = 0.340000+0.660000/(1.00000+EXP((STATES(1)+45.3000)/12.3000))
ALGBRC(17) = 1.00000/(1.00000+EXP(- (STATES(1)+36.5500)/13.7600))
ALGBRC(2) = 1.00000/(1.00000+( CONSTS(55)*CONSTS(6))/STATES(2) ** CONSTS(56))
ALGBRC(30) = 4000.00*(0.234845+(1.00000 - 0.234845)/(1.00000+STATES(2)/( CONSTS(55)*CONSTS(6)) ** CONSTS(56)))
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+35.9584)/9.24013))
ALGBRC(31) = 0.250000+7.00000/(1.00000+EXP((STATES(1)+38.0000)/10.0000))
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+57.0000)/8.00000))
ALGBRC(32) = 0.900000+1002.85/(1.00000+(STATES(1)+47.5000)/1.50000 ** 2.00000)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1)+22.0000)/7.00000))
ALGBRC(33) = 2.29000+5.70000/(1.00000+(STATES(1)+29.9700)/9.00000 ** 2.00000)
ALGBRC(34) = 90.9699*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1)+13.9629)/45.3782))*(1.00000+EXP(- (STATES(1)+9.49866)/3.39450))))
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+54.2300)/9.88000))
ALGBRC(35) = 0.450000+3.90000/(1.00000+(STATES(1)+66.0000)/26.0000 ** 2.00000)
ALGBRC(9) = 0.0200000+0.980000/(1.00000+EXP((STATES(1)+72.9780)/4.64000))
ALGBRC(36) = 150.000*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 417.430)/203.180))*(1.00000+EXP(- (STATES(1)+61.1100)/8.07000))))
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+16.0000)/9.50000))
ALGBRC(37) = 46.0999+1685.76/( (1.00000+EXP(- (STATES(1)+40.8489)/13.7802))*(1.00000+EXP((STATES(1)+20.6372)/15.1113)))
ALGBRC(38) = 475.667+16321.6/( (1.00000+EXP(- (STATES(1)+41.8328)/6.96673))*(1.00000+EXP((STATES(1)+23.2432)/21.2949)))
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+48.0000)/24.0000))
ALGBRC(39) = 19.7864/(1.00000+(STATES(1)+20.7136)/44.2868 ** 2.00000) - 0.378843
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+22.0000)/12.4800))
ALGBRC(40) = 395.300/(1.00000+(STATES(1)+38.1000)/33.5900 ** 2.00000)
ALGBRC(41) = (5503.00+5345.40/(1.00000+10.0000 ** (- 23.9000 - STATES(1))*- 0.0282700)) - 4590.60/(1.00000+10.0000 ** (STATES(1)+14.1500)*- 0.0357000)
ALGBRC(13) = 0.490000+0.510000/(1.00000+EXP((STATES(1)+1.08400)/28.7800))
ALGBRC(42) = 5.44000+29.2000/(1.00000+(STATES(1)+48.0900)/48.8300 ** 2.00000)
ALGBRC(15) = 1.00000/(1.00000+EXP(- (STATES(1)+15.0400)/16.9500))
ALGBRC(43) = 10.0000+895.900/(1.00000+EXP((- 18.0100 - STATES(1))/31.0400))
ALGBRC(16) = 0.405800/(1.00000+EXP((STATES(1)+86.8400)/15.0500))+0.594200/(1.00000+EXP((STATES(1) - 70.1300)/13.3700))
ALGBRC(44) = 1077.00+185845./(1.00000+(STATES(1) - 39.4400)/7.34400 ** 2.00000)
ALGBRC(45) = 37.5100+539.000/(1.00000+(STATES(1)+40.2400)/17.7200 ** 2.00000)
ALGBRC(18) = 0.490000+0.510000/(1.00000+EXP((STATES(1)+1.08400)/28.7800))
ALGBRC(46) = 5.44000+29.2000/(1.00000+(STATES(1)+48.0900)/48.8300 ** 2.00000)
ALGBRC(19) = 0.405800/(1.00000+EXP((STATES(1)+86.8400)/15.0500))+0.594200/(1.00000+EXP((STATES(1) - 70.1300)/13.3700))
ALGBRC(47) = 1077.00+185845./(1.00000+(STATES(1) - 39.4400)/7.34400 ** 2.00000)
ALGBRC(20) = 0.978613/(1.00000+EXP(- (STATES(1)+18.6736)/26.6606))
ALGBRC(48) = 500.000/(1.00000+(STATES(1)+60.7100)/15.7900 ** 2.00000)
ALGBRC(21) = 1.00000/(1.00000+EXP((STATES(1)+63.0000)/6.30000))
ALGBRC(49) = 5000.00/(1.00000+(STATES(1)+62.7133)/35.8611 ** 2.00000)
ALGBRC(50) = 30000.0+220000./(1.00000+EXP((STATES(1)+22.0000)/4.00000))
ALGBRC(22) = 0.948000/(1.00000+EXP(- (STATES(1)+17.9100)/18.4000))
ALGBRC(51) = 100.000/(1.00000+(STATES(1)+64.1000)/28.6700 ** 2.00000)
ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(1)+21.2000)/5.70000))
ALGBRC(52) = 1.00000e+06*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 315.000)/50.0000))*(1.00000+EXP(- (STATES(1)+74.9000)/8.00000))))
ALGBRC(53) = 2.50000e+06*(1.00000 - 1.00000/( (1.00000+EXP((STATES(1) - 132.868)/25.3992))*(1.00000+EXP(- (STATES(1)+24.9203)/2.67915))))
ALGBRC(24) = 1.00000/(1.00000+EXP(- (STATES(1)+27.7900)/7.57000))
ALGBRC(54) = 17.0000/(1.00000+(STATES(1)+20.5232)/35.0000 ** 2.00000)
ALGBRC(25) = 0.0200000+0.980000/(1.00000+EXP((STATES(1)+69.5000)/6.00000))
ALGBRC(55) = 7.50000+10.0000/(1.00000+(STATES(1)+34.1765)/120.000 ** 2.00000)
ALGBRC(26) = - 0.749234/(1.00000+( STATES(2)*CONSTS(7) - 0.0630535)/0.161942 ** 2.00000)+8.38384/(1.00000+( STATES(2)*CONSTS(7)+1538.29)/739.057 ** 2.00000)
ALGBRC(56) = 5011.47/(1.00000+( STATES(2)*CONSTS(7)+0.237503)/0.000239278 ** 0.422910) - 37.5137
ALGBRC(60) = 1.00000/(1.00000+EXP(( - ALGBRC(26)*CONSTS(73)*(STATES(1) - ALGBRC(56)))/( CONSTS(72)*CONSTS(75))))
ALGBRC(65) = 2.40914/(1.00000+(STATES(1) - 158.779)/- 52.1497 ** 2.00000)
ALGBRC(27) = - 0.681249/(1.00000+( STATES(2)*CONSTS(8) - 0.218988)/0.428335 ** 2.00000)+1.40001/(1.00000+( STATES(2)*CONSTS(8)+228.710)/684.946 ** 2.00000)
ALGBRC(57) = 8540.23/(1.00000+( STATES(2)*CONSTS(8)+0.401189)/0.00399115 ** 0.668054) - 109.275
ALGBRC(61) = 1.00000/(1.00000+EXP(( - ALGBRC(27)*CONSTS(73)*(STATES(1) - ALGBRC(57)))/( CONSTS(72)*CONSTS(75))))
ALGBRC(66) = 13.8049/(1.00000+(STATES(1) - 153.019)/66.4952 ** 2.00000)
ALGBRC(28) = 1.00000/(1.00000+EXP((STATES(1)+105.390)/8.65530))
ALGBRC(58) = 3.50000e-06*EXP( - 0.0497000*STATES(1))
ALGBRC(62) = 0.0400300*EXP( 0.0521100*STATES(1))
ALGBRC(67) = 1.00000/(ALGBRC(58)+ALGBRC(62))
ALGBRC(29) = ( STATES(1)*CONSTS(73))/( CONSTS(72)*CONSTS(75))
ALGBRC(63) = 0.000600000*EXP( 2.53000*ALGBRC(29))
ALGBRC(68) = 0.100000*EXP( - 5.00000*ALGBRC(29))
ALGBRC(70) = 1.00000/(1.00000+ ALGBRC(68)*(ALGBRC(63)/STATES(2) ** 2.00000+ALGBRC(63)/STATES(2)+1.00000))
ALGBRC(72) = - 160.000+210.000/(1.00000+EXP((STATES(1)+4.56000)/11.6200))+170.000/(1.00000+EXP(- (STATES(1)+25.5000)/11.6200))
ALGBRC(1) = TERNRY(VOI.GT.CONSTS(3).AND.VOI.LT.CONSTS(4), CONSTS(2), CONSTS(1))
ALGBRC(64) = CONSTS(17)*STATES(4)*STATES(4)*STATES(4)*STATES(5)*(STATES(1) - CONSTS(78))
ALGBRC(69) = 1.00000/(1.00000+STATES(2)/CONSTS(20) ** 4.00000)
ALGBRC(71) = CONSTS(18)*ALGBRC(69)*STATES(6)*STATES(6)*( 0.800000*STATES(7)+ 0.200000*STATES(8))*(STATES(1) - CONSTS(19))
ALGBRC(73) = CONSTS(21)*STATES(9)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(22))
ALGBRC(74) = CONSTS(24)*(STATES(1) - CONSTS(79))
ALGBRC(75) = CONSTS(33)*( 0.800000*STATES(11)+ 0.200000*STATES(12))*STATES(13)*(STATES(1) - CONSTS(79))
ALGBRC(76) = CONSTS(30)*( 0.300000*STATES(14)+ 0.700000*STATES(15))*STATES(16)*STATES(17)*(STATES(1) - CONSTS(79))
ALGBRC(77) = CONSTS(31)*STATES(18)*STATES(19)*(STATES(1) - CONSTS(79))
ALGBRC(78) = CONSTS(32)*( 0.200000*STATES(20)+ 0.800000*STATES(21))*STATES(22)*STATES(23)*(STATES(1) - CONSTS(79))
ALGBRC(79) = CONSTS(25)*STATES(24)*STATES(24)*( 0.380000*STATES(25)+ 0.630000*STATES(26))*(STATES(1) - CONSTS(79))
ALGBRC(80) = CONSTS(26)*STATES(27)*STATES(27)*( 0.750000*STATES(28)+ 0.250000*STATES(29))*(STATES(1) - CONSTS(79))
ALGBRC(81) = CONSTS(29)*STATES(30)*STATES(31)*(STATES(1) - CONSTS(79))
ALGBRC(82) = CONSTS(23)*CONSTS(27)*STATES(32)*(STATES(1) - CONSTS(79))
ALGBRC(83) = CONSTS(23)*CONSTS(28)*STATES(33)*(STATES(1) - CONSTS(79))
ALGBRC(84) = CONSTS(35)*STATES(34)*(STATES(1) - CONSTS(84))
ALGBRC(85) = CONSTS(34)*STATES(35)*(STATES(1) - CONSTS(80))
ALGBRC(59) = (( CONSTS(72)*CONSTS(75))/CONSTS(73))*log(( CONSTS(37)*CONSTS(64)+ CONSTS(38)*CONSTS(65)+( 4.00000*CONSTS(39)*CONSTS(13))/(1.00000+EXP(ALGBRC(29))))/( CONSTS(37)*CONSTS(12)+ CONSTS(38)*CONSTS(60)+( 4.00000*CONSTS(39)*STATES(2))/(1.00000+EXP(ALGBRC(29)))))
ALGBRC(89) = CONSTS(61)*CONSTS(68)*CONSTS(42)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(86) = CONSTS(61)*CONSTS(62)*CONSTS(41)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(91) = CONSTS(61)*CONSTS(66)*CONSTS(43)*CONSTS(36)*(STATES(1) - ALGBRC(59))
ALGBRC(92) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*ALGBRC(29))+ 0.00219000*EXP(CONSTS(65)/49.7100)*EXP( - 1.90000*ALGBRC(29)))
ALGBRC(93) = CONSTS(76)*CONSTS(88)*CONSTS(89)*ALGBRC(92)
ALGBRC(96) = 1.00000/(1.00000+CONSTS(48)/STATES(2) ** CONSTS(49))
ALGBRC(94) = EXP( (CONSTS(87) - 1.00000)*ALGBRC(29))
ALGBRC(95) = EXP( CONSTS(87)*ALGBRC(29))
ALGBRC(97) = CONSTS(60) ** 3.00000*CONSTS(13)*ALGBRC(95) - CONSTS(65) ** 3.00000*STATES(2)*ALGBRC(94)
ALGBRC(98) = 1.00000+ CONSTS(86)*ALGBRC(94)
ALGBRC(99) = CONSTS(53)*CONSTS(60) ** 3.00000+ CONSTS(52) ** 3.00000*STATES(2)+ CONSTS(50) ** 3.00000*CONSTS(13)*(1.00000+STATES(2)/CONSTS(51))
ALGBRC(100) = CONSTS(13)*CONSTS(60) ** 3.00000+ CONSTS(65) ** 3.00000*STATES(2)+ CONSTS(65) ** 3.00000*CONSTS(51)*(1.00000+CONSTS(60)/CONSTS(50) ** 3.00000)
ALGBRC(101) = ( CONSTS(11)*CONSTS(85)*ALGBRC(96)*ALGBRC(97))/( ALGBRC(98)*(ALGBRC(99)+ALGBRC(100)))
ALGBRC(102) = (( 0.500000*CONSTS(71)*CONSTS(73))/( CONSTS(70)*CONSTS(74)*CONSTS(69)))*CONSTS(10)*ALGBRC(101)
ALGBRC(105) = ALGBRC(64)+ALGBRC(84)+ALGBRC(102)+ALGBRC(93)+ALGBRC(71)+ALGBRC(73)+ALGBRC(85)+ALGBRC(75)+ALGBRC(76)+ALGBRC(77)+ALGBRC(78)+ALGBRC(79)+ALGBRC(80)+ALGBRC(81)+ALGBRC(82)+ALGBRC(83)+ALGBRC(89)+ALGBRC(91)+ALGBRC(86)+ALGBRC(74)
ALGBRC(87) = ALGBRC(71)+ALGBRC(73)+ALGBRC(86)
ALGBRC(88) = (( CONSTS(70)*CONSTS(74)*CONSTS(69))/( CONSTS(71)*CONSTS(73)))*ALGBRC(87)
ALGBRC(103) = CONSTS(45)/(1.00000+CONSTS(46)/STATES(2) ** CONSTS(47))
ALGBRC(106) = ALGBRC(88)+ALGBRC(101)+ALGBRC(103)
ALGBRC(3) = CONSTS(54)*(STATES(3) - 0.234500)
ALGBRC(90) = ALGBRC(88)*CONSTS(5)
ALGBRC(104) = ALGBRC(101)*CONSTS(5)
ALGBRC(107) = ALGBRC(103)*CONSTS(5)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END