C
C There are a total of 107 entries in the algebraic variable array.
C There are a total of 43 entries in each of the rate and state variable arrays.
C There are a total of 85 entries in the constant variable array.
C
C
C CONSTS(59) is Aj_nj in component cell (um2).
C CONSTS(1) is BCa in component calcium (mM).
C CONSTS(2) is CSQN in component calcium (mM).
C STATES(1) is CaSR1 in component calcium (mM).
C STATES(2) is CaSR2 in component calcium (mM).
C STATES(3) is CaSR3 in component calcium (mM).
C STATES(4) is CaSR4 in component calcium (mM).
C STATES(5) is Cai1 in component calcium (mM).
C STATES(6) is Cai2 in component calcium (mM).
C STATES(7) is Cai3 in component calcium (mM).
C STATES(8) is Cai4 in component calcium (mM).
C STATES(9) is Cass in component calcium (mM).
C CONSTS(3) is DCa in component calcium (m2_per_s_times_1e_minus_12).
C CONSTS(4) is DCaBm in component calcium (m2_per_s_times_1e_minus_12).
C CONSTS(5) is DCaSR in component calcium (m2_per_s_times_1e_minus_12).
C CONSTS(6) is F in component phys (C_per_mol).
C ALGBRC(59) is ICaL in component ical (pA).
C ALGBRC(64) is ICaP in component icap (pA).
C ALGBRC(82) is ICab in component icab (pA).
C ALGBRC(72) is INaCa in component inaca (pA).
C ALGBRC(76) is JCa1 in component calcium (kat_times_1e_minus_12).
C ALGBRC(77) is JCa2 in component calcium (kat_times_1e_minus_12).
C ALGBRC(78) is JCa3 in component calcium (kat_times_1e_minus_12).
C ALGBRC(49) is JCa4 in component calcium (kat_times_1e_minus_12).
C ALGBRC(94) is JCass in component calcium (kat_times_1e_minus_12).
C ALGBRC(73) is JSRCa1 in component calcium (kat_times_1e_minus_12).
C ALGBRC(74) is JSRCa2 in component calcium (kat_times_1e_minus_12).
C ALGBRC(75) is JSRCa3 in component calcium (kat_times_1e_minus_12).
C ALGBRC(92) is JSRCa4 in component calcium (kat_times_1e_minus_12).
C ALGBRC(1) is JSRCaleak1 in component calcium (kat_times_1e_minus_12).
C ALGBRC(2) is JSRCaleak2 in component calcium (kat_times_1e_minus_12).
C ALGBRC(3) is JSRCaleak3 in component calcium (kat_times_1e_minus_12).
C ALGBRC(4) is JSRCaleakss in component calcium (kat_times_1e_minus_12).
C ALGBRC(65) is J_SERCASR1 in component serca (kat_times_1e_minus_12).
C ALGBRC(66) is J_SERCASR2 in component serca (kat_times_1e_minus_12).
C ALGBRC(67) is J_SERCASR3 in component serca (kat_times_1e_minus_12).
C ALGBRC(89) is J_SERCASRss in component serca (kat_times_1e_minus_12).
C ALGBRC(69) is J_bulkSERCA1 in component serca (kat_times_1e_minus_12).
C ALGBRC(70) is J_bulkSERCA2 in component serca (kat_times_1e_minus_12).
C ALGBRC(71) is J_bulkSERCA3 in component serca (kat_times_1e_minus_12).
C ALGBRC(91) is J_bulkSERCAss in component serca (kat_times_1e_minus_12).
C ALGBRC(32) is Jj_nj in component calcium (kat_times_1e_minus_12).
C ALGBRC(61) is Jrel1 in component ryr (kat_times_1e_minus_12).
C ALGBRC(62) is Jrel2 in component ryr (kat_times_1e_minus_12).
C ALGBRC(63) is Jrel3 in component ryr (kat_times_1e_minus_12).
C ALGBRC(87) is Jrelss in component ryr (kat_times_1e_minus_12).
C CONSTS(7) is KdBCa in component calcium (mM).
C CONSTS(8) is KdCSQN in component calcium (mM).
C CONSTS(9) is KdSLhigh in component calcium (mM).
C CONSTS(10) is KdSLlow in component calcium (mM).
C CONSTS(11) is SLhigh in component calcium (mM).
C CONSTS(12) is SLlow in component calcium (mM).
C CONSTS(73) is VSR1 in component cell (nL).
C CONSTS(78) is VSR2 in component cell (nL).
C CONSTS(81) is VSR3 in component cell (nL).
C CONSTS(84) is VSR4 in component cell (nL).
C CONSTS(60) is Vnonjunct1 in component cell (nL).
C CONSTS(72) is Vnonjunct2 in component cell (nL).
C CONSTS(77) is Vnonjunct3 in component cell (nL).
C CONSTS(80) is Vnonjunct4 in component cell (nL).
C CONSTS(13) is Vss in component cell (nL).
C ALGBRC(29) is calcium_CaSR1_beta in component calcium (dimensionless).
C ALGBRC(30) is calcium_CaSR2_beta in component calcium (dimensionless).
C ALGBRC(31) is calcium_CaSR3_beta in component calcium (dimensionless).
C ALGBRC(50) is calcium_CaSR4_beta in component calcium (dimensionless).
C ALGBRC(51) is calcium_Cai1_beta in component calcium (dimensionless).
C ALGBRC(46) is calcium_Cai1_gamma in component calcium (dimensionless).
C ALGBRC(52) is calcium_Cai2_beta in component calcium (dimensionless).
C ALGBRC(47) is calcium_Cai2_gamma in component calcium (dimensionless).
C ALGBRC(53) is calcium_Cai3_beta in component calcium (dimensionless).
C ALGBRC(48) is calcium_Cai3_gamma in component calcium (dimensionless).
C ALGBRC(60) is calcium_Cai4_beta in component calcium (dimensionless).
C ALGBRC(54) is calcium_Cai4_gamma in component calcium (dimensionless).
C ALGBRC(55) is calcium_Cass_beta in component calcium (dimensionless).
C ALGBRC(83) is calcium_Cass_i_tot in component calcium (pA).
C CONSTS(14) is dx in component cell (um).
C CONSTS(15) is kSRleak in component calcium (hertz).
C VOI is time in component engine (second).
C CONSTS(61) is xj_nj in component cell (um).
C CONSTS(16) is Cm in component cell (nF).
C CONSTS(85) is Vcytosol in component cell (nL).
C CONSTS(83) is Vnonjunct_Nai in component cell (nL).
C CONSTS(17) is lcell in component cell (um).
C CONSTS(18) is pi in component cell (dimensionless).
C CONSTS(19) is rjunct in component cell (um).
C CONSTS(62) is xj_nj_Nai in component cell (um).
C CONSTS(20) is Cao in component extra (mM).
C CONSTS(21) is Ko in component extra (mM).
C CONSTS(22) is Nao in component extra (mM).
C ALGBRC(81) is ECa in component nernst (mV).
C STATES(10) is V in component membrane (mV).
C CONSTS(23) is gCab in component icab (nS).
C CONSTS(24) is ECa_app in component ical (mV).
C STATES(11) is d in component ical (hertz).
C STATES(12) is f1 in component ical (hertz).
C STATES(13) is f2 in component ical (hertz).
C ALGBRC(5) is f_inf in component ical (dimensionless).
C STATES(14) is fca in component ical (hertz).
C CONSTS(25) is gCaL in component ical (nS).
C ALGBRC(6) is ical_d_inf in component ical (dimensionless).
C ALGBRC(33) is ical_d_tau in component ical (second).
C ALGBRC(34) is ical_f1_tau in component ical (second).
C ALGBRC(35) is ical_f2_tau in component ical (second).
C ALGBRC(7) is ical_fca_inf in component ical (dimensionless).
C CONSTS(26) is ical_fca_tau in component ical (second).
C CONSTS(27) is kCa in component ical (mM).
C CONSTS(28) is kCan in component ical (dimensionless).
C CONSTS(29) is ICaPmax in component icap (pA).
C CONSTS(30) is kCaP in component icap (mM).
C ALGBRC(84) is EK in component nernst (mV).
C ALGBRC(97) is ENa in component nernst (mV).
C ALGBRC(99) is If in component if (pA).
C ALGBRC(86) is IfK in component if (pA).
C ALGBRC(98) is IfNa in component if (pA).
C CONSTS(31) is gIf in component if (nS).
C ALGBRC(8) is if_y_inf in component if (dimensionless).
C ALGBRC(36) is if_y_tau in component if (second).
C STATES(15) is y in component if (hertz).
C CONSTS(74) is FRT in component phys (per_mV).
C ALGBRC(88) is IK1 in component ik1 (pA).
C CONSTS(63) is gK1 in component ik1 (nS).
C ALGBRC(90) is IKr in component ikr (pA).
C CONSTS(32) is gKr in component ikr (nS).
C ALGBRC(9) is ikr_pa_inf in component ikr (dimensionless).
C ALGBRC(37) is ikr_pa_tau in component ikr (second).
C STATES(16) is pa in component ikr (hertz).
C ALGBRC(68) is pi in component ikr (dimensionless).
C ALGBRC(93) is IKs in component iks (pA).
C CONSTS(33) is gKs in component iks (nS).
C ALGBRC(10) is iks_n_inf in component iks (dimensionless).
C ALGBRC(38) is iks_n_tau in component iks (second).
C STATES(17) is n in component iks (hertz).
C ALGBRC(95) is IKur in component ikur (pA).
C CONSTS(64) is gKur in component ikur (nS).
C ALGBRC(11) is ikur_r_inf in component ikur (dimensionless).
C ALGBRC(39) is ikur_r_tau in component ikur (second).
C ALGBRC(12) is ikur_s_inf in component ikur (dimensionless).
C ALGBRC(40) is ikur_s_tau in component ikur (second).
C STATES(18) is r in component ikur (hertz).
C STATES(19) is s in component ikur (hertz).
C ALGBRC(100) is INa in component ina (pA).
C CONSTS(34) is PNa in component ina (m3_per_s_times_1e_minus_12).
C STATES(20) is h1 in component ina (hertz).
C STATES(21) is h2 in component ina (hertz).
C ALGBRC(13) is h_inf in component ina (dimensionless).
C ALGBRC(41) is ina_h1_tau in component ina (second).
C ALGBRC(42) is ina_h2_tau in component ina (second).
C ALGBRC(14) is ina_m_inf in component ina (dimensionless).
C ALGBRC(43) is ina_m_tau in component ina (second).
C STATES(22) is m in component ina (hertz).
C ALGBRC(101) is INab in component inab (pA).
C CONSTS(35) is gNab in component inab (nS).
C STATES(23) is Nass in component sodium (mM).
C CONSTS(36) is dNaCa in component inaca (m12_per_mol4).
C CONSTS(37) is fCaNCX in component inaca (dimensionless).
C CONSTS(38) is gam in component inaca (dimensionless).
C CONSTS(39) is kNaCa in component inaca (m12_A_per_mol4_times_1e_minus_12).
C ALGBRC(80) is INaK in component inak (pA).
C CONSTS(40) is INaKmax in component inak (pA).
C ALGBRC(79) is Nass15 in component inak (dimensionless).
C CONSTS(41) is kNaKK in component inak (mM).
C CONSTS(42) is kNaKNa in component inak (dimensionless).
C ALGBRC(96) is It in component it (pA).
C CONSTS(65) is gt in component it (nS).
C ALGBRC(15) is it_r_inf in component it (dimensionless).
C ALGBRC(44) is it_r_tau in component it (second).
C ALGBRC(16) is it_s_inf in component it (dimensionless).
C ALGBRC(45) is it_s_tau in component it (second).
C STATES(24) is r in component it (hertz).
C STATES(25) is s in component it (hertz).
C ALGBRC(102) is i_ion in component membrane (pA).
C ALGBRC(104) is i_stim in component stimulus (pA).
C STATES(26) is Ki in component potassium (mM).
C CONSTS(66) is RTF in component phys (mV).
C CONSTS(43) is R in component phys (mJ_per_mol_per_K).
C CONSTS(44) is T in component phys (kelvin).
C ALGBRC(106) is i_tot in component potassium (pA).
C ALGBRC(56) is SRCa1 in component ryr (dimensionless).
C ALGBRC(57) is SRCa2 in component ryr (dimensionless).
C ALGBRC(58) is SRCa3 in component ryr (dimensionless).
C ALGBRC(85) is SRCass in component ryr (dimensionless).
C STATES(27) is a1 in component ryr (hertz).
C STATES(28) is a2 in component ryr (hertz).
C STATES(29) is a3 in component ryr (hertz).
C ALGBRC(17) is ainf1 in component ryr (dimensionless).
C ALGBRC(18) is ainf2 in component ryr (dimensionless).
C ALGBRC(19) is ainf3 in component ryr (dimensionless).
C ALGBRC(20) is ainfss in component ryr (dimensionless).
C STATES(30) is ass in component ryr (hertz).
C STATES(31) is c1 in component ryr (hertz).
C STATES(32) is c2 in component ryr (hertz).
C STATES(33) is c3 in component ryr (hertz).
C ALGBRC(21) is cinf1 in component ryr (dimensionless).
C ALGBRC(22) is cinf2 in component ryr (dimensionless).
C ALGBRC(23) is cinf3 in component ryr (dimensionless).
C ALGBRC(24) is cinfss in component ryr (dimensionless).
C STATES(34) is css in component ryr (hertz).
C CONSTS(75) is nu1 in component ryr (m3_per_s_times_1e_minus_12).
C CONSTS(79) is nu2 in component ryr (m3_per_s_times_1e_minus_12).
C CONSTS(82) is nu3 in component ryr (m3_per_s_times_1e_minus_12).
C CONSTS(67) is nuss in component ryr (m3_per_s_times_1e_minus_12).
C STATES(35) is o1 in component ryr (hertz).
C STATES(36) is o2 in component ryr (hertz).
C STATES(37) is o3 in component ryr (hertz).
C ALGBRC(25) is oinf1 in component ryr (dimensionless).
C ALGBRC(26) is oinf2 in component ryr (dimensionless).
C ALGBRC(27) is oinf3 in component ryr (dimensionless).
C ALGBRC(28) is oinfss in component ryr (dimensionless).
C STATES(38) is oss in component ryr (hertz).
C CONSTS(45) is tau_act in component ryr (second).
C CONSTS(46) is tau_actss in component ryr (second).
C CONSTS(47) is tau_adapt in component ryr (second).
C CONSTS(48) is tau_inact in component ryr (second).
C CONSTS(49) is tau_inactss in component ryr (second).
C CONSTS(50) is SERCAKmf in component serca (mM).
C CONSTS(51) is SERCAKmr in component serca (mM).
C STATES(39) is a1 in component serca (mM).
C STATES(40) is a2 in component serca (mM).
C STATES(41) is a3 in component serca (mM).
C STATES(42) is ass in component serca (mM).
C CONSTS(52) is cpumps in component serca (mM).
C CONSTS(68) is k1 in component serca (m6_per_s_per_mol2).
C CONSTS(76) is k2 in component serca (hertz).
C CONSTS(69) is k3 in component serca (m6_per_s_per_mol2).
C CONSTS(53) is k4 in component serca (hertz).
C CONSTS(70) is BNa in component sodium (mM).
C CONSTS(54) is DNa in component sodium (m2_per_s_times_1e_minus_12).
C ALGBRC(103) is JNa in component sodium (kat_times_1e_minus_12).
C CONSTS(55) is KdBNa in component sodium (mM).
C STATES(43) is Nai in component sodium (mM).
C ALGBRC(105) is betaNass in component sodium (dimensionless).
C ALGBRC(107) is i_ss in component sodium (pA).
C CONSTS(71) is amplitude in component stimulus (pA).
C CONSTS(56) is duration in component stimulus (second).
C CONSTS(57) is offset in component stimulus (second).
C CONSTS(58) is period in component stimulus (second).
C RATES(1) is d/dt CaSR1 in component calcium (mM).
C RATES(2) is d/dt CaSR2 in component calcium (mM).
C RATES(3) is d/dt CaSR3 in component calcium (mM).
C RATES(4) is d/dt CaSR4 in component calcium (mM).
C RATES(5) is d/dt Cai1 in component calcium (mM).
C RATES(6) is d/dt Cai2 in component calcium (mM).
C RATES(7) is d/dt Cai3 in component calcium (mM).
C RATES(8) is d/dt Cai4 in component calcium (mM).
C RATES(9) is d/dt Cass in component calcium (mM).
C RATES(11) is d/dt d in component ical (hertz).
C RATES(12) is d/dt f1 in component ical (hertz).
C RATES(13) is d/dt f2 in component ical (hertz).
C RATES(14) is d/dt fca in component ical (hertz).
C RATES(15) is d/dt y in component if (hertz).
C RATES(16) is d/dt pa in component ikr (hertz).
C RATES(17) is d/dt n in component iks (hertz).
C RATES(18) is d/dt r in component ikur (hertz).
C RATES(19) is d/dt s in component ikur (hertz).
C RATES(20) is d/dt h1 in component ina (hertz).
C RATES(21) is d/dt h2 in component ina (hertz).
C RATES(22) is d/dt m in component ina (hertz).
C RATES(24) is d/dt r in component it (hertz).
C RATES(25) is d/dt s in component it (hertz).
C RATES(10) is d/dt V in component membrane (mV).
C RATES(26) is d/dt Ki in component potassium (mM).
C RATES(27) is d/dt a1 in component ryr (hertz).
C RATES(28) is d/dt a2 in component ryr (hertz).
C RATES(29) is d/dt a3 in component ryr (hertz).
C RATES(30) is d/dt ass in component ryr (hertz).
C RATES(31) is d/dt c1 in component ryr (hertz).
C RATES(32) is d/dt c2 in component ryr (hertz).
C RATES(33) is d/dt c3 in component ryr (hertz).
C RATES(34) is d/dt css in component ryr (hertz).
C RATES(35) is d/dt o1 in component ryr (hertz).
C RATES(36) is d/dt o2 in component ryr (hertz).
C RATES(37) is d/dt o3 in component ryr (hertz).
C RATES(38) is d/dt oss in component ryr (hertz).
C RATES(39) is d/dt a1 in component serca (mM).
C RATES(40) is d/dt a2 in component serca (mM).
C RATES(41) is d/dt a3 in component serca (mM).
C RATES(42) is d/dt ass in component serca (mM).
C RATES(43) is d/dt Nai in component sodium (mM).
C RATES(23) is d/dt Nass in component sodium (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.024
      CONSTS(2) = 6.7
      STATES(1) = 0.6189225
      STATES(2) = 0.6076289
      STATES(3) = 0.5905266
      STATES(4) = 0.5738108
      STATES(5) = 1.35496500000000013e-04
      STATES(6) = 1.38142100000000014e-04
      STATES(7) = 1.44208699999999994e-04
      STATES(8) = 1.56184399999999995e-04
      STATES(9) = 1.61937700000000013e-04
      CONSTS(3) = 780.0
      CONSTS(4) = 25.0
      CONSTS(5) = 44.0
      CONSTS(6) = 96487.0
      CONSTS(7) = 0.00238
      CONSTS(8) = 0.8
      CONSTS(9) = 0.013
      CONSTS(10) = 1.1
      CONSTS(11) = 13.0
      CONSTS(12) = 165.0
      CONSTS(13) = 4.99231999999999966e-05
      CONSTS(14) = 1.625
      CONSTS(15) = 0.006
      CONSTS(16) = 0.05
      CONSTS(17) = 122.051
      CONSTS(18) = 3.14159265358979312e+00
      CONSTS(19) = 6.5
      CONSTS(20) = 1.8
      CONSTS(21) = 5.4
      CONSTS(22) = 130.0
      STATES(10) = -75.42786
      CONSTS(23) = 0.0952
      CONSTS(24) = 60.0
      STATES(11) = 1.06091699999999996e-05
      STATES(12) = 0.9988566
      STATES(13) = 0.9988624
      STATES(14) = 0.9744374
      CONSTS(25) = 25.3125
      CONSTS(26) = 0.002
      CONSTS(27) = 0.001
      CONSTS(28) = 2.0
      CONSTS(29) = 2.0
      CONSTS(30) = 0.0005
      CONSTS(31) = 1.0
      STATES(15) = 5.62066499999999969e-02
      CONSTS(32) = 0.5
      STATES(16) = 4.18941700000000008e-05
      CONSTS(33) = 1.0
      STATES(17) = 4.10975100000000003e-03
      STATES(18) = 3.11170299999999984e-04
      STATES(19) = 0.9751094
      CONSTS(34) = 0.0018
      STATES(20) = 0.90391
      STATES(21) = 0.9039673
      STATES(22) = 2.77581199999999990e-03
      CONSTS(35) = 0.060599
      STATES(23) = 8.691504
      CONSTS(36) = 0.0003
      CONSTS(37) = 1.0
      CONSTS(38) = 0.45
      CONSTS(39) = 0.0084
      CONSTS(40) = 70.8253
      CONSTS(41) = 1.0
      CONSTS(42) = 11.0
      STATES(24) = 9.59425800000000026e-04
      STATES(25) = 0.954338
      STATES(26) = 134.6313
      CONSTS(43) = 8314.0
      CONSTS(44) = 306.15
      STATES(27) = 0.1925362
      STATES(28) = 0.2010345
      STATES(29) = 0.2163122
      STATES(30) = 0.2455297
      STATES(31) = 0.9993722
      STATES(32) = 0.9995086
      STATES(33) = 0.9995604
      STATES(34) = 0.9999717
      STATES(35) = 9.47851400000000044e-05
      STATES(36) = 7.76550300000000031e-05
      STATES(37) = 5.67494700000000006e-05
      STATES(38) = 3.97509699999999973e-05
      CONSTS(45) = 0.01875
      CONSTS(46) = 0.005
      CONSTS(47) = 1.0
      CONSTS(48) = 0.0875
      CONSTS(49) = 0.015
      CONSTS(50) = 0.00025
      CONSTS(51) = 1.8
      STATES(39) = 4.63856499999999988e-03
      STATES(40) = 4.51207800000000010e-03
      STATES(41) = 4.32640899999999981e-03
      STATES(42) = 4.25044500000000026e-03
      CONSTS(52) = 0.04
      CONSTS(53) = 7.5
      CONSTS(54) = 0.12
      CONSTS(55) = 10.0
      STATES(43) = 9.28686
      CONSTS(56) = 0.001
      CONSTS(57) = 0.01
      CONSTS(58) = 1.0
      CONSTS(59) =  ( ( ( CONSTS(18)*CONSTS(19))*2.00000)*CONSTS(17))*0.500000
      CONSTS(60) =  ( ( ( ( 1.00000*CONSTS(14) ** 2.00000 -  0.00000*CONSTS(14) ** 2.00000)*CONSTS(18))*CONSTS(17))*0.500000)*1.00000e-06
      CONSTS(61) = 0.0200000/2.00000+CONSTS(14)/2.00000
      CONSTS(62) = 0.0200000/2.00000+ 2.00000*CONSTS(14)
      CONSTS(63) =  3.82500*0.900000
      CONSTS(64) =  0.890000*2.75000
      CONSTS(65) =  1.09000*7.50000
      CONSTS(66) = ( CONSTS(43)*CONSTS(44))/CONSTS(6)
      CONSTS(67) =  625.000*CONSTS(13)
      CONSTS(68) =  1000.00 ** 2.00000*CONSTS(53)
      CONSTS(69) = CONSTS(53)/CONSTS(51) ** 2.00000
      CONSTS(70) =  0.490000*2.31000
      CONSTS(71) = - 2500.00
      CONSTS(72) =  ( ( ( ( 2.00000*CONSTS(14) ** 2.00000 -  1.00000*CONSTS(14) ** 2.00000)*CONSTS(18))*CONSTS(17))*0.500000)*1.00000e-06
      CONSTS(73) =  (( 0.0500000*CONSTS(60))/2.00000)*0.900000
      CONSTS(74) = 1.00000/CONSTS(66)
      CONSTS(75) =  1.00000*CONSTS(60)
      CONSTS(76) =  CONSTS(68)*CONSTS(50) ** 2.00000
      CONSTS(77) =  ( ( ( ( 3.00000*CONSTS(14) ** 2.00000 -  2.00000*CONSTS(14) ** 2.00000)*CONSTS(18))*CONSTS(17))*0.500000)*1.00000e-06
      CONSTS(78) =  (( 0.0500000*CONSTS(72))/2.00000)*0.900000
      CONSTS(79) =  1.00000*CONSTS(72)
      CONSTS(80) =  ( ( ( ( 4.00000*CONSTS(14) ** 2.00000 -  3.00000*CONSTS(14) ** 2.00000)*CONSTS(18))*CONSTS(17))*0.500000)*1.00000e-06
      CONSTS(81) =  (( 0.0500000*CONSTS(77))/2.00000)*0.900000
      CONSTS(82) =  1.00000*CONSTS(77)
      CONSTS(83) = ((CONSTS(60)+CONSTS(72))+CONSTS(77))+CONSTS(80)
      CONSTS(84) =  (( 0.0500000*CONSTS(80))/2.00000)*0.900000
      CONSTS(85) = CONSTS(83)+CONSTS(13)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = 1.00000 - 1.00000/(1.00000+CONSTS(27)/STATES(9) ** CONSTS(28))
      RATES(14) = (ALGBRC(7) - STATES(14))/CONSTS(26)
      ALGBRC(17) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(5)*1000.00 - 0.290000)/0.0820000))
      RATES(27) = (ALGBRC(17) - STATES(27))/CONSTS(47)
      ALGBRC(18) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(6)*1000.00 - 0.290000)/0.0820000))
      RATES(28) = (ALGBRC(18) - STATES(28))/CONSTS(47)
      ALGBRC(19) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(7)*1000.00 - 0.290000)/0.0820000))
      RATES(29) = (ALGBRC(19) - STATES(29))/CONSTS(47)
      ALGBRC(20) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(9)*1000.00 - 0.290000)/0.0820000))
      RATES(30) = (ALGBRC(20) - STATES(30))/CONSTS(47)
      ALGBRC(21) = 1.00000/(1.00000+EXP(( STATES(5)*1000.00 - (STATES(27)+0.0200000))/0.0100000))
      RATES(31) = (ALGBRC(21) - STATES(31))/CONSTS(48)
      ALGBRC(22) = 1.00000/(1.00000+EXP(( STATES(6)*1000.00 - (STATES(28)+0.0200000))/0.0100000))
      RATES(32) = (ALGBRC(22) - STATES(32))/CONSTS(48)
      ALGBRC(23) = 1.00000/(1.00000+EXP(( STATES(7)*1000.00 - (STATES(29)+0.0200000))/0.0100000))
      RATES(33) = (ALGBRC(23) - STATES(33))/CONSTS(48)
      ALGBRC(24) = 1.00000/(1.00000+EXP(( STATES(9)*1000.00 - (STATES(30)+0.0200000))/0.0100000))
      RATES(34) = (ALGBRC(24) - STATES(34))/CONSTS(49)
      ALGBRC(25) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(5)*1000.00 - (STATES(27)+0.220000))/0.0300000))
      RATES(35) = (ALGBRC(25) - STATES(35))/CONSTS(45)
      ALGBRC(26) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(6)*1000.00 - (STATES(28)+0.220000))/0.0300000))
      RATES(36) = (ALGBRC(26) - STATES(36))/CONSTS(45)
      ALGBRC(27) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(7)*1000.00 - (STATES(29)+0.220000))/0.0300000))
      RATES(37) = (ALGBRC(27) - STATES(37))/CONSTS(45)
      ALGBRC(28) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(9)*1000.00 - (STATES(30)+0.220000))/0.0300000))
      RATES(38) = (ALGBRC(28) - STATES(38))/CONSTS(46)
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(10)+9.00000)/- 5.80000))
      ALGBRC(33) =  0.00270000*EXP(- (STATES(10)+35.0000)/30.0000 ** 2.00000)+0.00200000
      RATES(11) = (ALGBRC(6) - STATES(11))/ALGBRC(33)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(10)+27.4000)/7.10000))
      ALGBRC(34) = (( 0.986980*EXP(- (STATES(10)+30.1605)/7.09396 ** 2.00000)+0.0427500/(1.00000+EXP((STATES(10) - 51.6155)/- 80.6133)))+0.0357600/(1.00000+EXP((STATES(10)+29.5727)/13.2176))) - 0.00821000
      RATES(12) = (ALGBRC(5) - STATES(12))/ALGBRC(34)
      ALGBRC(35) =  1.33230*EXP(- (STATES(10)+40.0000)/14.2000 ** 2.00000)+0.0626000
      RATES(13) = (ALGBRC(5) - STATES(13))/ALGBRC(35)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(10)+97.8287)/12.4802))
      ALGBRC(36) = 1.00000/( 0.00332000*EXP(- STATES(10)/16.5410)+ 23.7184*EXP(STATES(10)/16.5410))
      RATES(15) = (ALGBRC(8) - STATES(15))/ALGBRC(36)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(10)+15.0000)/- 6.00000))
      ALGBRC(37) =  0.217180*EXP(- (STATES(10)+20.1376)/22.1996 ** 2.00000)+0.0311800
      RATES(16) = (ALGBRC(9) - STATES(16))/ALGBRC(37)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(10) - 19.9000)/- 12.7000))
      ALGBRC(38) =  0.400000*EXP(- (STATES(10) - 20.0000)/20.0000 ** 2.00000)+0.700000
      RATES(17) = (ALGBRC(10) - STATES(17))/ALGBRC(38)
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(10)+6.00000)/- 8.60000))
      ALGBRC(39) = 0.00900000/(1.00000+EXP((STATES(10)+5.00000)/12.0000))+0.000500000
      RATES(18) = (ALGBRC(11) - STATES(18))/ALGBRC(39)
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(10)+7.50000)/10.0000))
      ALGBRC(40) = 0.590000/(1.00000+EXP((STATES(10)+60.0000)/10.0000))+3.05000
      RATES(19) = (ALGBRC(12) - STATES(19))/ALGBRC(40)
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(10)+63.6000)/5.30000))
      ALGBRC(41) = 0.0300000/(1.00000+EXP((STATES(10)+35.1000)/3.20000))+0.000300000
      RATES(20) = (ALGBRC(13) - STATES(20))/ALGBRC(41)
      ALGBRC(42) = 0.120000/(1.00000+EXP((STATES(10)+35.1000)/3.20000))+0.00300000
      RATES(21) = (ALGBRC(13) - STATES(21))/ALGBRC(42)
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(10)+27.1200)/- 8.21000))
      ALGBRC(43) =  4.20000e-05*EXP(- (STATES(10)+25.5700)/28.8000 ** 2.00000)+2.40000e-05
      RATES(22) = (ALGBRC(14) - STATES(22))/ALGBRC(43)
      ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(10) - 1.00000)/- 11.0000))
      ALGBRC(44) =  0.00350000*EXP(- (STATES(10)+0.00000)/30.0000 ** 2.00000)+0.00150000
      RATES(24) = (ALGBRC(15) - STATES(24))/ALGBRC(44)
      ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(10)+40.5000)/11.5000))
      ALGBRC(45) =  0.0256350*EXP(- (STATES(10)+52.4500)/15.8827 ** 2.00000)+0.0141400
      RATES(25) = (ALGBRC(16) - STATES(25))/ALGBRC(45)
      ALGBRC(32) =  ( (( CONSTS(3)*CONSTS(59))/CONSTS(61))*(STATES(9) - STATES(8)))*1.00000e-06
      ALGBRC(49) = ALGBRC(32)
      ALGBRC(54) = ( CONSTS(1)*CONSTS(7))/STATES(8)+CONSTS(7) ** 2.00000
      ALGBRC(60) = 1.00000/(1.00000+ALGBRC(54))
      RATES(8) = ( ( ALGBRC(60)*(CONSTS(3)+ ALGBRC(54)*CONSTS(4)))*(((STATES(8) -  2.00000*STATES(8))+STATES(7))/CONSTS(14) ** 2.00000+(STATES(8) - STATES(7))/( ( 2.00000*4.00000)*CONSTS(14) ** 2.00000)) -  (( ( ( 2.00000*ALGBRC(60))*ALGBRC(54))*CONSTS(4))/(CONSTS(7)+STATES(8)))*(STATES(8) - STATES(7))/( 2.00000*CONSTS(14)) ** 2.00000)+ (ALGBRC(49)/CONSTS(80))*ALGBRC(60)
      ALGBRC(1) =  ( CONSTS(15)*(STATES(1) - STATES(5)))*CONSTS(60)
      ALGBRC(65) =  ( ( ( - CONSTS(69)*STATES(1) ** 2.00000)*(CONSTS(52) - STATES(39))+ CONSTS(53)*STATES(39))*CONSTS(60))*2.00000
      ALGBRC(56) = 1.00000 - 1.00000/(1.00000+EXP((STATES(1) - 0.300000)/0.100000))
      ALGBRC(61) =  ( ( ( CONSTS(75)*STATES(35))*STATES(31))*ALGBRC(56))*(STATES(1) - STATES(5))
      ALGBRC(73) = (ALGBRC(65) - ALGBRC(1)) - ALGBRC(61)
      ALGBRC(29) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(1)+CONSTS(8) ** 2.00000)
      RATES(1) =  ( ALGBRC(29)*CONSTS(5))*(((STATES(2) -  2.00000*STATES(1))+STATES(1))/CONSTS(14) ** 2.00000+(STATES(2) - STATES(1))/( ( 2.00000*1.00000)*CONSTS(14) ** 2.00000))+ (ALGBRC(73)/CONSTS(73))*ALGBRC(29)
      ALGBRC(2) =  ( CONSTS(15)*(STATES(2) - STATES(6)))*CONSTS(72)
      ALGBRC(66) =  ( ( ( - CONSTS(69)*STATES(2) ** 2.00000)*(CONSTS(52) - STATES(40))+ CONSTS(53)*STATES(40))*CONSTS(72))*2.00000
      ALGBRC(57) = 1.00000 - 1.00000/(1.00000+EXP((STATES(2) - 0.300000)/0.100000))
      ALGBRC(62) =  ( ( ( CONSTS(79)*STATES(36))*STATES(32))*ALGBRC(57))*(STATES(2) - STATES(6))
      ALGBRC(74) = (ALGBRC(66) - ALGBRC(2)) - ALGBRC(62)
      ALGBRC(30) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(2)+CONSTS(8) ** 2.00000)
      RATES(2) =  ( ALGBRC(30)*CONSTS(5))*(((STATES(3) -  2.00000*STATES(2))+STATES(1))/CONSTS(14) ** 2.00000+(STATES(3) - STATES(1))/( ( 2.00000*2.00000)*CONSTS(14) ** 2.00000))+ (ALGBRC(74)/CONSTS(78))*ALGBRC(30)
      ALGBRC(3) =  ( CONSTS(15)*(STATES(3) - STATES(7)))*CONSTS(77)
      ALGBRC(67) =  ( ( ( - CONSTS(69)*STATES(3) ** 2.00000)*(CONSTS(52) - STATES(41))+ CONSTS(53)*STATES(41))*CONSTS(77))*2.00000
      ALGBRC(58) = 1.00000 - 1.00000/(1.00000+EXP((STATES(3) - 0.300000)/0.100000))
      ALGBRC(63) =  ( ( ( CONSTS(82)*STATES(37))*STATES(33))*ALGBRC(58))*(STATES(3) - STATES(7))
      ALGBRC(75) = (ALGBRC(67) - ALGBRC(3)) - ALGBRC(63)
      ALGBRC(31) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(3)+CONSTS(8) ** 2.00000)
      RATES(3) =  ( ALGBRC(31)*CONSTS(5))*(((STATES(4) -  2.00000*STATES(3))+STATES(2))/CONSTS(14) ** 2.00000+(STATES(4) - STATES(2))/( ( 2.00000*3.00000)*CONSTS(14) ** 2.00000))+ (ALGBRC(75)/CONSTS(81))*ALGBRC(31)
      ALGBRC(69) =  ( ( ( CONSTS(68)*STATES(5) ** 2.00000)*(CONSTS(52) - STATES(39)) -  CONSTS(76)*STATES(39))*CONSTS(60))*2.00000
      RATES(39) = ( 0.500000*(- ALGBRC(65)+ALGBRC(69)))/CONSTS(60)
      ALGBRC(70) =  ( ( ( CONSTS(68)*STATES(6) ** 2.00000)*(CONSTS(52) - STATES(40)) -  CONSTS(76)*STATES(40))*CONSTS(72))*2.00000
      RATES(40) = ( 0.500000*(- ALGBRC(66)+ALGBRC(70)))/CONSTS(72)
      ALGBRC(71) =  ( ( ( CONSTS(68)*STATES(7) ** 2.00000)*(CONSTS(52) - STATES(41)) -  CONSTS(76)*STATES(41))*CONSTS(77))*2.00000
      RATES(41) = ( 0.500000*(- ALGBRC(67)+ALGBRC(71)))/CONSTS(77)
      ALGBRC(76) = (- ALGBRC(69)+ALGBRC(1))+ALGBRC(61)
      ALGBRC(46) = ( CONSTS(1)*CONSTS(7))/STATES(5)+CONSTS(7) ** 2.00000
      ALGBRC(51) = 1.00000/(1.00000+ALGBRC(46))
      RATES(5) = ( ( ALGBRC(51)*(CONSTS(3)+ ALGBRC(46)*CONSTS(4)))*(((STATES(6) -  2.00000*STATES(5))+STATES(5))/CONSTS(14) ** 2.00000+(STATES(6) - STATES(5))/( ( 2.00000*1.00000)*CONSTS(14) ** 2.00000)) -  (( ( ( 2.00000*ALGBRC(51))*ALGBRC(46))*CONSTS(4))/(CONSTS(7)+STATES(5)))*(STATES(6) - STATES(5))/( 2.00000*CONSTS(14)) ** 2.00000)+ (ALGBRC(76)/CONSTS(60))*ALGBRC(51)
      ALGBRC(77) = (- ALGBRC(70)+ALGBRC(2))+ALGBRC(62)
      ALGBRC(47) = ( CONSTS(1)*CONSTS(7))/STATES(6)+CONSTS(7) ** 2.00000
      ALGBRC(52) = 1.00000/(1.00000+ALGBRC(47))
      RATES(6) = ( ( ALGBRC(52)*(CONSTS(3)+ ALGBRC(47)*CONSTS(4)))*(((STATES(7) -  2.00000*STATES(6))+STATES(5))/CONSTS(14) ** 2.00000+(STATES(7) - STATES(5))/( ( 2.00000*2.00000)*CONSTS(14) ** 2.00000)) -  (( ( ( 2.00000*ALGBRC(52))*ALGBRC(47))*CONSTS(4))/(CONSTS(7)+STATES(6)))*(STATES(7) - STATES(5))/( 2.00000*CONSTS(14)) ** 2.00000)+ (ALGBRC(77)/CONSTS(72))*ALGBRC(52)
      ALGBRC(78) = (- ALGBRC(71)+ALGBRC(3))+ALGBRC(63)
      ALGBRC(48) = ( CONSTS(1)*CONSTS(7))/STATES(7)+CONSTS(7) ** 2.00000
      ALGBRC(53) = 1.00000/(1.00000+ALGBRC(48))
      RATES(7) = ( ( ALGBRC(53)*(CONSTS(3)+ ALGBRC(48)*CONSTS(4)))*(((STATES(8) -  2.00000*STATES(7))+STATES(6))/CONSTS(14) ** 2.00000+(STATES(8) - STATES(6))/( ( 2.00000*3.00000)*CONSTS(14) ** 2.00000)) -  (( ( ( 2.00000*ALGBRC(53))*ALGBRC(48))*CONSTS(4))/(CONSTS(7)+STATES(7)))*(STATES(8) - STATES(6))/( 2.00000*CONSTS(14)) ** 2.00000)+ (ALGBRC(78)/CONSTS(77))*ALGBRC(53)
      ALGBRC(4) =  ( CONSTS(15)*(STATES(4) - STATES(9)))*CONSTS(13)
      ALGBRC(89) =  ( ( ( - CONSTS(69)*STATES(4) ** 2.00000)*(CONSTS(52) - STATES(42))+ CONSTS(53)*STATES(42))*CONSTS(13))*2.00000
      ALGBRC(85) = 1.00000 - 1.00000/(1.00000+EXP((STATES(4) - 0.300000)/0.100000))
      ALGBRC(87) =  ( ( ( CONSTS(67)*STATES(38))*STATES(34))*ALGBRC(85))*(STATES(4) - STATES(9))
      ALGBRC(92) = (ALGBRC(89) - ALGBRC(4)) - ALGBRC(87)
      ALGBRC(50) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(4)+CONSTS(8) ** 2.00000)
      RATES(4) =  ( ALGBRC(50)*CONSTS(5))*(((STATES(4) -  2.00000*STATES(4))+STATES(3))/CONSTS(14) ** 2.00000+(STATES(4) - STATES(3))/( ( 2.00000*4.00000)*CONSTS(14) ** 2.00000))+ (ALGBRC(92)/CONSTS(84))*ALGBRC(50)
      ALGBRC(91) =  ( ( ( CONSTS(68)*STATES(9) ** 2.00000)*(CONSTS(52) - STATES(42)) -  CONSTS(76)*STATES(42))*CONSTS(13))*2.00000
      RATES(42) = ( 0.500000*(- ALGBRC(89)+ALGBRC(91)))/CONSTS(13)
      ALGBRC(94) = ((- ALGBRC(32)+ALGBRC(4)) - ALGBRC(91))+ALGBRC(87)
      ALGBRC(55) = 1.00000/(((1.00000+( CONSTS(12)*CONSTS(10))/STATES(9)+CONSTS(10) ** 2.00000)+( CONSTS(11)*CONSTS(9))/STATES(9)+CONSTS(9) ** 2.00000)+( CONSTS(1)*CONSTS(7))/STATES(9)+CONSTS(7) ** 2.00000)
      ALGBRC(59) =  ( ( ( ( CONSTS(25)*STATES(11))*STATES(14))*STATES(12))*STATES(13))*(STATES(10) - CONSTS(24))
      ALGBRC(64) = ( CONSTS(29)*STATES(9))/(CONSTS(30)+STATES(9))
      ALGBRC(81) = ( CONSTS(66)*log(CONSTS(20)/STATES(9)))/2.00000
      ALGBRC(82) =  CONSTS(23)*(STATES(10) - ALGBRC(81))
      ALGBRC(72) = ( CONSTS(39)*( ( EXP( ( CONSTS(38)*STATES(10))*CONSTS(74))*STATES(23) ** 3.00000)*CONSTS(20) -  ( ( EXP( ( (CONSTS(38) - 1.00000)*STATES(10))*CONSTS(74))*CONSTS(22) ** 3.00000)*STATES(9))*CONSTS(37)))/(1.00000+ CONSTS(36)*( ( CONSTS(22) ** 3.00000*STATES(9))*CONSTS(37)+ STATES(23) ** 3.00000*CONSTS(20)))
      ALGBRC(83) = ((- ALGBRC(59) - ALGBRC(82)) - ALGBRC(64))+ 2.00000*ALGBRC(72)
      RATES(9) =  ALGBRC(55)*(ALGBRC(94)/CONSTS(13)+ALGBRC(83)/( ( 2.00000*CONSTS(13))*CONSTS(6)))
      ALGBRC(103) =  ( (( CONSTS(54)*CONSTS(59))/CONSTS(62))*(STATES(23) - STATES(43)))*1.00000e-06
      RATES(43) = ALGBRC(103)/CONSTS(83)
      ALGBRC(84) =  CONSTS(66)*log(CONSTS(21)/STATES(26))
      ALGBRC(86) =  ( CONSTS(31)*STATES(15))*( (1.00000 - 0.267700)*(STATES(10) - ALGBRC(84)))
      ALGBRC(97) =  CONSTS(66)*log(CONSTS(22)/STATES(23))
      ALGBRC(98) =  ( CONSTS(31)*STATES(15))*( 0.267700*(STATES(10) - ALGBRC(97)))
      ALGBRC(99) = ALGBRC(86)+ALGBRC(98)
      ALGBRC(88) = ( ( CONSTS(63)* CONSTS(21)*1.00000 ** 0.445700)*(STATES(10) - ALGBRC(84)))/(1.00000+EXP( ( 1.50000*((STATES(10) - ALGBRC(84))+3.60000))*CONSTS(74)))
      ALGBRC(68) = 1.00000/(1.00000+EXP((STATES(10)+55.0000)/24.0000))
      ALGBRC(90) =  ( ( CONSTS(32)*STATES(16))*ALGBRC(68))*(STATES(10) - ALGBRC(84))
      ALGBRC(93) =  ( CONSTS(33)*STATES(17))*(STATES(10) - ALGBRC(84))
      ALGBRC(95) =  ( ( CONSTS(64)*STATES(18))*STATES(19))*(STATES(10) - ALGBRC(84))
      ALGBRC(100) = ( ( ( ( ( ( ( CONSTS(34)*STATES(22) ** 3.00000)*( 0.900000*STATES(20)+ 0.100000*STATES(21)))*CONSTS(22))*STATES(10))*CONSTS(6))*CONSTS(74))*(EXP( (STATES(10) - ALGBRC(97))*CONSTS(74)) - 1.00000))/(EXP( STATES(10)*CONSTS(74)) - 1.00000)
      ALGBRC(101) =  CONSTS(35)*(STATES(10) - ALGBRC(97))
      ALGBRC(79) =  STATES(23)*1.00000 ** 1.50000
      ALGBRC(80) = ( (( (( CONSTS(40)*CONSTS(21))/(CONSTS(21)+CONSTS(41)))*ALGBRC(79))/(ALGBRC(79)+CONSTS(42) ** 1.50000))*(STATES(10)+150.000))/(STATES(10)+200.000)
      ALGBRC(96) =  ( ( CONSTS(65)*STATES(24))*STATES(25))*(STATES(10) - ALGBRC(84))
      ALGBRC(102) = ((((((((((((ALGBRC(100))+ALGBRC(59))+ALGBRC(96))+ALGBRC(95))+ALGBRC(88))+ALGBRC(90))+ALGBRC(93))+ALGBRC(101))+ALGBRC(82))+ALGBRC(80))+ALGBRC(64))+ALGBRC(72))+ALGBRC(99)
      ALGBRC(104) =  (TERNRY((VOI - CONSTS(57)) -  CONSTS(58)*INT((VOI - CONSTS(57))/CONSTS(58)).LT.CONSTS(56), 1.00000, 0.00000))*CONSTS(71)
      RATES(10) = - (ALGBRC(102)+ALGBRC(104))/CONSTS(16)
      ALGBRC(106) = (((((((ALGBRC(96))+ALGBRC(95))+ALGBRC(88))+ALGBRC(90))+ALGBRC(93)) -  2.00000*ALGBRC(80))+ALGBRC(86))+ALGBRC(104)
      RATES(26) = - ALGBRC(106)/( CONSTS(85)*CONSTS(6))
      ALGBRC(105) = 1.00000/(1.00000+( CONSTS(70)*CONSTS(55))/STATES(23)+CONSTS(55) ** 2.00000)
      ALGBRC(107) = ((((ALGBRC(100))+ALGBRC(101))+ 3.00000*ALGBRC(80))+ 3.00000*ALGBRC(72))+ALGBRC(98)
      RATES(23) =  ALGBRC(105)*(- ALGBRC(103)/CONSTS(13) - ALGBRC(107)/( CONSTS(13)*CONSTS(6)))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(7) = 1.00000 - 1.00000/(1.00000+CONSTS(27)/STATES(9) ** CONSTS(28))
      ALGBRC(17) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(5)*1000.00 - 0.290000)/0.0820000))
      ALGBRC(18) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(6)*1000.00 - 0.290000)/0.0820000))
      ALGBRC(19) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(7)*1000.00 - 0.290000)/0.0820000))
      ALGBRC(20) = 0.505000 - 0.427000/(1.00000+EXP(( STATES(9)*1000.00 - 0.290000)/0.0820000))
      ALGBRC(21) = 1.00000/(1.00000+EXP(( STATES(5)*1000.00 - (STATES(27)+0.0200000))/0.0100000))
      ALGBRC(22) = 1.00000/(1.00000+EXP(( STATES(6)*1000.00 - (STATES(28)+0.0200000))/0.0100000))
      ALGBRC(23) = 1.00000/(1.00000+EXP(( STATES(7)*1000.00 - (STATES(29)+0.0200000))/0.0100000))
      ALGBRC(24) = 1.00000/(1.00000+EXP(( STATES(9)*1000.00 - (STATES(30)+0.0200000))/0.0100000))
      ALGBRC(25) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(5)*1000.00 - (STATES(27)+0.220000))/0.0300000))
      ALGBRC(26) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(6)*1000.00 - (STATES(28)+0.220000))/0.0300000))
      ALGBRC(27) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(7)*1000.00 - (STATES(29)+0.220000))/0.0300000))
      ALGBRC(28) = 1.00000 - 1.00000/(1.00000+EXP(( STATES(9)*1000.00 - (STATES(30)+0.220000))/0.0300000))
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(10)+9.00000)/- 5.80000))
      ALGBRC(33) =  0.00270000*EXP(- (STATES(10)+35.0000)/30.0000 ** 2.00000)+0.00200000
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(10)+27.4000)/7.10000))
      ALGBRC(34) = (( 0.986980*EXP(- (STATES(10)+30.1605)/7.09396 ** 2.00000)+0.0427500/(1.00000+EXP((STATES(10) - 51.6155)/- 80.6133)))+0.0357600/(1.00000+EXP((STATES(10)+29.5727)/13.2176))) - 0.00821000
      ALGBRC(35) =  1.33230*EXP(- (STATES(10)+40.0000)/14.2000 ** 2.00000)+0.0626000
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(10)+97.8287)/12.4802))
      ALGBRC(36) = 1.00000/( 0.00332000*EXP(- STATES(10)/16.5410)+ 23.7184*EXP(STATES(10)/16.5410))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(10)+15.0000)/- 6.00000))
      ALGBRC(37) =  0.217180*EXP(- (STATES(10)+20.1376)/22.1996 ** 2.00000)+0.0311800
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(10) - 19.9000)/- 12.7000))
      ALGBRC(38) =  0.400000*EXP(- (STATES(10) - 20.0000)/20.0000 ** 2.00000)+0.700000
      ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(10)+6.00000)/- 8.60000))
      ALGBRC(39) = 0.00900000/(1.00000+EXP((STATES(10)+5.00000)/12.0000))+0.000500000
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(10)+7.50000)/10.0000))
      ALGBRC(40) = 0.590000/(1.00000+EXP((STATES(10)+60.0000)/10.0000))+3.05000
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(10)+63.6000)/5.30000))
      ALGBRC(41) = 0.0300000/(1.00000+EXP((STATES(10)+35.1000)/3.20000))+0.000300000
      ALGBRC(42) = 0.120000/(1.00000+EXP((STATES(10)+35.1000)/3.20000))+0.00300000
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(10)+27.1200)/- 8.21000))
      ALGBRC(43) =  4.20000e-05*EXP(- (STATES(10)+25.5700)/28.8000 ** 2.00000)+2.40000e-05
      ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(10) - 1.00000)/- 11.0000))
      ALGBRC(44) =  0.00350000*EXP(- (STATES(10)+0.00000)/30.0000 ** 2.00000)+0.00150000
      ALGBRC(16) = 1.00000/(1.00000+EXP((STATES(10)+40.5000)/11.5000))
      ALGBRC(45) =  0.0256350*EXP(- (STATES(10)+52.4500)/15.8827 ** 2.00000)+0.0141400
      ALGBRC(32) =  ( (( CONSTS(3)*CONSTS(59))/CONSTS(61))*(STATES(9) - STATES(8)))*1.00000e-06
      ALGBRC(49) = ALGBRC(32)
      ALGBRC(54) = ( CONSTS(1)*CONSTS(7))/STATES(8)+CONSTS(7) ** 2.00000
      ALGBRC(60) = 1.00000/(1.00000+ALGBRC(54))
      ALGBRC(1) =  ( CONSTS(15)*(STATES(1) - STATES(5)))*CONSTS(60)
      ALGBRC(65) =  ( ( ( - CONSTS(69)*STATES(1) ** 2.00000)*(CONSTS(52) - STATES(39))+ CONSTS(53)*STATES(39))*CONSTS(60))*2.00000
      ALGBRC(56) = 1.00000 - 1.00000/(1.00000+EXP((STATES(1) - 0.300000)/0.100000))
      ALGBRC(61) =  ( ( ( CONSTS(75)*STATES(35))*STATES(31))*ALGBRC(56))*(STATES(1) - STATES(5))
      ALGBRC(73) = (ALGBRC(65) - ALGBRC(1)) - ALGBRC(61)
      ALGBRC(29) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(1)+CONSTS(8) ** 2.00000)
      ALGBRC(2) =  ( CONSTS(15)*(STATES(2) - STATES(6)))*CONSTS(72)
      ALGBRC(66) =  ( ( ( - CONSTS(69)*STATES(2) ** 2.00000)*(CONSTS(52) - STATES(40))+ CONSTS(53)*STATES(40))*CONSTS(72))*2.00000
      ALGBRC(57) = 1.00000 - 1.00000/(1.00000+EXP((STATES(2) - 0.300000)/0.100000))
      ALGBRC(62) =  ( ( ( CONSTS(79)*STATES(36))*STATES(32))*ALGBRC(57))*(STATES(2) - STATES(6))
      ALGBRC(74) = (ALGBRC(66) - ALGBRC(2)) - ALGBRC(62)
      ALGBRC(30) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(2)+CONSTS(8) ** 2.00000)
      ALGBRC(3) =  ( CONSTS(15)*(STATES(3) - STATES(7)))*CONSTS(77)
      ALGBRC(67) =  ( ( ( - CONSTS(69)*STATES(3) ** 2.00000)*(CONSTS(52) - STATES(41))+ CONSTS(53)*STATES(41))*CONSTS(77))*2.00000
      ALGBRC(58) = 1.00000 - 1.00000/(1.00000+EXP((STATES(3) - 0.300000)/0.100000))
      ALGBRC(63) =  ( ( ( CONSTS(82)*STATES(37))*STATES(33))*ALGBRC(58))*(STATES(3) - STATES(7))
      ALGBRC(75) = (ALGBRC(67) - ALGBRC(3)) - ALGBRC(63)
      ALGBRC(31) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(3)+CONSTS(8) ** 2.00000)
      ALGBRC(69) =  ( ( ( CONSTS(68)*STATES(5) ** 2.00000)*(CONSTS(52) - STATES(39)) -  CONSTS(76)*STATES(39))*CONSTS(60))*2.00000
      ALGBRC(70) =  ( ( ( CONSTS(68)*STATES(6) ** 2.00000)*(CONSTS(52) - STATES(40)) -  CONSTS(76)*STATES(40))*CONSTS(72))*2.00000
      ALGBRC(71) =  ( ( ( CONSTS(68)*STATES(7) ** 2.00000)*(CONSTS(52) - STATES(41)) -  CONSTS(76)*STATES(41))*CONSTS(77))*2.00000
      ALGBRC(76) = (- ALGBRC(69)+ALGBRC(1))+ALGBRC(61)
      ALGBRC(46) = ( CONSTS(1)*CONSTS(7))/STATES(5)+CONSTS(7) ** 2.00000
      ALGBRC(51) = 1.00000/(1.00000+ALGBRC(46))
      ALGBRC(77) = (- ALGBRC(70)+ALGBRC(2))+ALGBRC(62)
      ALGBRC(47) = ( CONSTS(1)*CONSTS(7))/STATES(6)+CONSTS(7) ** 2.00000
      ALGBRC(52) = 1.00000/(1.00000+ALGBRC(47))
      ALGBRC(78) = (- ALGBRC(71)+ALGBRC(3))+ALGBRC(63)
      ALGBRC(48) = ( CONSTS(1)*CONSTS(7))/STATES(7)+CONSTS(7) ** 2.00000
      ALGBRC(53) = 1.00000/(1.00000+ALGBRC(48))
      ALGBRC(4) =  ( CONSTS(15)*(STATES(4) - STATES(9)))*CONSTS(13)
      ALGBRC(89) =  ( ( ( - CONSTS(69)*STATES(4) ** 2.00000)*(CONSTS(52) - STATES(42))+ CONSTS(53)*STATES(42))*CONSTS(13))*2.00000
      ALGBRC(85) = 1.00000 - 1.00000/(1.00000+EXP((STATES(4) - 0.300000)/0.100000))
      ALGBRC(87) =  ( ( ( CONSTS(67)*STATES(38))*STATES(34))*ALGBRC(85))*(STATES(4) - STATES(9))
      ALGBRC(92) = (ALGBRC(89) - ALGBRC(4)) - ALGBRC(87)
      ALGBRC(50) = 1.00000/(1.00000+( CONSTS(2)*CONSTS(8))/STATES(4)+CONSTS(8) ** 2.00000)
      ALGBRC(91) =  ( ( ( CONSTS(68)*STATES(9) ** 2.00000)*(CONSTS(52) - STATES(42)) -  CONSTS(76)*STATES(42))*CONSTS(13))*2.00000
      ALGBRC(94) = ((- ALGBRC(32)+ALGBRC(4)) - ALGBRC(91))+ALGBRC(87)
      ALGBRC(55) = 1.00000/(((1.00000+( CONSTS(12)*CONSTS(10))/STATES(9)+CONSTS(10) ** 2.00000)+( CONSTS(11)*CONSTS(9))/STATES(9)+CONSTS(9) ** 2.00000)+( CONSTS(1)*CONSTS(7))/STATES(9)+CONSTS(7) ** 2.00000)
      ALGBRC(59) =  ( ( ( ( CONSTS(25)*STATES(11))*STATES(14))*STATES(12))*STATES(13))*(STATES(10) - CONSTS(24))
      ALGBRC(64) = ( CONSTS(29)*STATES(9))/(CONSTS(30)+STATES(9))
      ALGBRC(81) = ( CONSTS(66)*log(CONSTS(20)/STATES(9)))/2.00000
      ALGBRC(82) =  CONSTS(23)*(STATES(10) - ALGBRC(81))
      ALGBRC(72) = ( CONSTS(39)*( ( EXP( ( CONSTS(38)*STATES(10))*CONSTS(74))*STATES(23) ** 3.00000)*CONSTS(20) -  ( ( EXP( ( (CONSTS(38) - 1.00000)*STATES(10))*CONSTS(74))*CONSTS(22) ** 3.00000)*STATES(9))*CONSTS(37)))/(1.00000+ CONSTS(36)*( ( CONSTS(22) ** 3.00000*STATES(9))*CONSTS(37)+ STATES(23) ** 3.00000*CONSTS(20)))
      ALGBRC(83) = ((- ALGBRC(59) - ALGBRC(82)) - ALGBRC(64))+ 2.00000*ALGBRC(72)
      ALGBRC(103) =  ( (( CONSTS(54)*CONSTS(59))/CONSTS(62))*(STATES(23) - STATES(43)))*1.00000e-06
      ALGBRC(84) =  CONSTS(66)*log(CONSTS(21)/STATES(26))
      ALGBRC(86) =  ( CONSTS(31)*STATES(15))*( (1.00000 - 0.267700)*(STATES(10) - ALGBRC(84)))
      ALGBRC(97) =  CONSTS(66)*log(CONSTS(22)/STATES(23))
      ALGBRC(98) =  ( CONSTS(31)*STATES(15))*( 0.267700*(STATES(10) - ALGBRC(97)))
      ALGBRC(99) = ALGBRC(86)+ALGBRC(98)
      ALGBRC(88) = ( ( CONSTS(63)* CONSTS(21)*1.00000 ** 0.445700)*(STATES(10) - ALGBRC(84)))/(1.00000+EXP( ( 1.50000*((STATES(10) - ALGBRC(84))+3.60000))*CONSTS(74)))
      ALGBRC(68) = 1.00000/(1.00000+EXP((STATES(10)+55.0000)/24.0000))
      ALGBRC(90) =  ( ( CONSTS(32)*STATES(16))*ALGBRC(68))*(STATES(10) - ALGBRC(84))
      ALGBRC(93) =  ( CONSTS(33)*STATES(17))*(STATES(10) - ALGBRC(84))
      ALGBRC(95) =  ( ( CONSTS(64)*STATES(18))*STATES(19))*(STATES(10) - ALGBRC(84))
      ALGBRC(100) = ( ( ( ( ( ( ( CONSTS(34)*STATES(22) ** 3.00000)*( 0.900000*STATES(20)+ 0.100000*STATES(21)))*CONSTS(22))*STATES(10))*CONSTS(6))*CONSTS(74))*(EXP( (STATES(10) - ALGBRC(97))*CONSTS(74)) - 1.00000))/(EXP( STATES(10)*CONSTS(74)) - 1.00000)
      ALGBRC(101) =  CONSTS(35)*(STATES(10) - ALGBRC(97))
      ALGBRC(79) =  STATES(23)*1.00000 ** 1.50000
      ALGBRC(80) = ( (( (( CONSTS(40)*CONSTS(21))/(CONSTS(21)+CONSTS(41)))*ALGBRC(79))/(ALGBRC(79)+CONSTS(42) ** 1.50000))*(STATES(10)+150.000))/(STATES(10)+200.000)
      ALGBRC(96) =  ( ( CONSTS(65)*STATES(24))*STATES(25))*(STATES(10) - ALGBRC(84))
      ALGBRC(102) = ((((((((((((ALGBRC(100))+ALGBRC(59))+ALGBRC(96))+ALGBRC(95))+ALGBRC(88))+ALGBRC(90))+ALGBRC(93))+ALGBRC(101))+ALGBRC(82))+ALGBRC(80))+ALGBRC(64))+ALGBRC(72))+ALGBRC(99)
      ALGBRC(104) =  (TERNRY((VOI - CONSTS(57)) -  CONSTS(58)*INT((VOI - CONSTS(57))/CONSTS(58)).LT.CONSTS(56), 1.00000, 0.00000))*CONSTS(71)
      ALGBRC(106) = (((((((ALGBRC(96))+ALGBRC(95))+ALGBRC(88))+ALGBRC(90))+ALGBRC(93)) -  2.00000*ALGBRC(80))+ALGBRC(86))+ALGBRC(104)
      ALGBRC(105) = 1.00000/(1.00000+( CONSTS(70)*CONSTS(55))/STATES(23)+CONSTS(55) ** 2.00000)
      ALGBRC(107) = ((((ALGBRC(100))+ALGBRC(101))+ 3.00000*ALGBRC(80))+ 3.00000*ALGBRC(72))+ALGBRC(98)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END