C
C There are a total of 89 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 151 entries in the constant variable array.
C
C
C CONSTS(1) is Buf_C in component cai (mM).
C STATES(1) is Cai in component cai (mM).
C ALGBRC(1) is Cai_bufc in component cai (dimensionless).
C CONSTS(2) is Cm in component geom (pF).
C CONSTS(3) is F in component phys (C_per_mmol).
C CONSTS(4) is Kbuf_C in component cai (mM).
C CONSTS(145) is Vc in component geom (um3).
C ALGBRC(66) is i_CaL_Ca in component ical (A_per_F).
C ALGBRC(72) is i_CaT in component icat (A_per_F).
C ALGBRC(84) is i_NaCa in component inaca (A_per_F).
C ALGBRC(86) is i_PCa in component ipca (A_per_F).
C ALGBRC(63) is i_b_Ca in component ibca (A_per_F).
C ALGBRC(82) is i_leak in component ileak (mM_per_ms).
C ALGBRC(87) is i_rel in component irel (mM_per_ms).
C ALGBRC(89) is i_up in component iup (mM_per_ms).
C VOI is time in component engine (ms).
C CONSTS(87) is Buf_SR in component casr (mM).
C STATES(2) is Ca_SR in component casr (mM).
C ALGBRC(17) is Ca_SR_bufSR in component casr (dimensionless).
C CONSTS(5) is Kbuf_SR in component casr (mM).
C CONSTS(127) is V_SR in component geom (um3).
C CONSTS(6) is pace in component engine (dimensionless).
C CONSTS(7) is Cao in component extra (mM).
C ALGBRC(33) is E_Ca in component erev (mV).
C ALGBRC(47) is E_K in component erev (mV).
C ALGBRC(61) is E_Na in component erev (mV).
C STATES(3) is Ki in component ki (mM).
C CONSTS(8) is Ko in component extra (mM).
C STATES(4) is Nai in component nai (mM).
C CONSTS(9) is Nao in component extra (mM).
C CONSTS(88) is RTF in component phys (mV).
C CONSTS(10) is VSR_tenT in component geom (um3).
C CONSTS(11) is V_tot in component geom (um3).
C CONSTS(86) is V_tot_tenT in component geom (um3).
C CONSTS(12) is Vc_tenT in component geom (um3).
C STATES(5) is V in component membrane (mV).
C CONSTS(89) is g_b_Ca in component ibca (mS_per_uF).
C CONSTS(90) is g_b_Na in component ibna (mS_per_uF).
C ALGBRC(64) is i_b_Na in component ibna (A_per_F).
C CONSTS(128) is FFRT in component phys (s4_A2_per_g_per_m2_per_mol_times_1e3).
C CONSTS(91) is FRT in component phys (per_mV).
C ALGBRC(4) is alpha_fCa in component ical (dimensionless).
C ALGBRC(20) is beta_fCa in component ical (dimensionless).
C STATES(6) is d in component ical (dimensionless).
C CONSTS(13) is d1 in component ical (mS_per_uF).
C CONSTS(14) is d2 in component ical (mV).
C CONSTS(92) is d3 in component ical (mS_per_uF).
C CONSTS(129) is d4 in component ical (mV).
C CONSTS(15) is d5 in component ical (dimensionless).
C CONSTS(93) is d6 in component ical (mV).
C STATES(7) is f in component ical (dimensionless).
C CONSTS(16) is f1 in component ical (mS_per_uF).
C CONSTS(94) is f2 in component ical (mV).
C CONSTS(95) is f3 in component ical (mS_per_uF).
C CONSTS(130) is f4 in component ical (mV).
C CONSTS(17) is f5 in component ical (dimensionless).
C CONSTS(18) is f6 in component ical (mV).
C STATES(8) is fCa in component ical (dimensionless).
C ALGBRC(50) is fCa_inf in component ical (dimensionless).
C ALGBRC(36) is gamma_fCa in component ical (dimensionless).
C ALGBRC(71) is i_CaL in component ical (A_per_F).
C ALGBRC(68) is i_CaL_K in component ical (A_per_F).
C ALGBRC(70) is i_CaL_Na in component ical (A_per_F).
C ALGBRC(65) is ibarca in component ical (A_per_F).
C ALGBRC(67) is ibark in component ical (A_per_F).
C ALGBRC(69) is ibarna in component ical (A_per_F).
C ALGBRC(2) is ical_d_a in component ical (mS_per_uF).
C ALGBRC(18) is ical_d_b in component ical (mS_per_uF).
C ALGBRC(34) is ical_d_inf in component ical (dimensionless).
C ALGBRC(48) is ical_d_tau in component ical (ms).
C ALGBRC(3) is ical_f_a in component ical (mS_per_uF).
C ALGBRC(19) is ical_f_b in component ical (mS_per_uF).
C ALGBRC(35) is ical_f_inf in component ical (dimensionless).
C ALGBRC(49) is ical_f_tau in component ical (ms).
C ALGBRC(62) is k_fca in component ical (dimensionless).
C CONSTS(19) is p_CaL in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(146) is p_CaL_Ca in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(147) is p_CaL_K in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(148) is p_CaL_Na in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(20) is p_CaL_shannonCa in component ical (dimensionless).
C CONSTS(131) is p_CaL_shannonCap in component ical (dimensionless).
C CONSTS(21) is p_CaL_shannonK in component ical (dimensionless).
C CONSTS(132) is p_CaL_shannonKp in component ical (dimensionless).
C CONSTS(22) is p_CaL_shannonNa in component ical (dimensionless).
C CONSTS(133) is p_CaL_shannonNap in component ical (dimensionless).
C CONSTS(96) is p_CaL_shannonTot in component ical (dimensionless).
C CONSTS(23) is scale in component ical (dimensionless).
C CONSTS(24) is tau_fCa in component ical (ms).
C CONSTS(25) is taud_const in component ical (ms).
C CONSTS(26) is tauf_const in component ical (ms).
C STATES(9) is d in component icat (dimensionless).
C STATES(10) is f in component icat (dimensionless).
C CONSTS(27) is g_CaT in component icat (mS_per_uF).
C ALGBRC(5) is icat_d_inf in component icat (dimensionless).
C ALGBRC(21) is icat_d_tau in component icat (ms).
C ALGBRC(6) is icat_f_inf in component icat (dimensionless).
C ALGBRC(22) is icat_f_tau in component icat (ms).
C CONSTS(97) is Na_frac in component ifunny (dimensionless).
C CONSTS(28) is NatoK_ratio in component ifunny (dimensionless).
C STATES(11) is Xf in component ifunny (dimensionless).
C CONSTS(29) is g_f in component ifunny (mS_per_uF).
C ALGBRC(75) is i_f in component ifunny (A_per_F).
C ALGBRC(73) is i_fK in component ifunny (A_per_F).
C ALGBRC(74) is i_fNa in component ifunny (A_per_F).
C ALGBRC(7) is ifunny_Xf_a in component ifunny (mS_per_uF).
C ALGBRC(23) is ifunny_Xf_b in component ifunny (mS_per_uF).
C ALGBRC(37) is ifunny_Xf_inf in component ifunny (dimensionless).
C ALGBRC(51) is ifunny_Xf_tau in component ifunny (ms).
C CONSTS(30) is xF1 in component ifunny (mS_per_uF).
C CONSTS(98) is xF2 in component ifunny (mV).
C CONSTS(99) is xF3 in component ifunny (mS_per_uF).
C CONSTS(134) is xF4 in component ifunny (mV).
C CONSTS(31) is xF5 in component ifunny (dimensionless).
C CONSTS(32) is xF6 in component ifunny (mV).
C CONSTS(33) is xF_const in component ifunny (ms).
C CONSTS(34) is g_K1 in component ik1 (mS_per_uF).
C ALGBRC(79) is i_K1 in component ik1 (A_per_F).
C ALGBRC(76) is ik1_inf_a in component ik1 (mS_per_uF).
C ALGBRC(77) is ik1_inf_b in component ik1 (mS_per_uF).
C ALGBRC(78) is inf in component ik1 (dimensionless).
C CONSTS(35) is xK11 in component ik1 (mS_per_uF).
C CONSTS(36) is xK12 in component ik1 (mV).
C CONSTS(37) is xK13 in component ik1 (mV).
C CONSTS(38) is xK14 in component ik1 (mV).
C CONSTS(39) is xK15 in component ik1 (mV).
C STATES(12) is Xr1 in component ikr (dimensionless).
C CONSTS(40) is Xr1_1 in component ikr (mS_per_uF).
C CONSTS(41) is Xr1_2 in component ikr (mV).
C CONSTS(100) is Xr1_3 in component ikr (mS_per_uF).
C CONSTS(135) is Xr1_4 in component ikr (mV).
C CONSTS(42) is Xr1_5 in component ikr (dimensionless).
C CONSTS(101) is Xr1_6 in component ikr (mV).
C STATES(13) is Xr2 in component ikr (dimensionless).
C CONSTS(43) is Xr2_1 in component ikr (mS_per_uF).
C CONSTS(102) is Xr2_2 in component ikr (mV).
C CONSTS(103) is Xr2_3 in component ikr (mS_per_uF).
C CONSTS(136) is Xr2_4 in component ikr (mV).
C CONSTS(44) is Xr2_5 in component ikr (dimensionless).
C CONSTS(45) is Xr2_6 in component ikr (mV).
C CONSTS(46) is g_Kr in component ikr (mS_per_uF).
C ALGBRC(80) is i_Kr in component ikr (A_per_F).
C ALGBRC(8) is ikr_Xr1_a in component ikr (mS_per_uF).
C ALGBRC(24) is ikr_Xr1_b in component ikr (mS_per_uF).
C ALGBRC(38) is ikr_Xr1_inf in component ikr (dimensionless).
C ALGBRC(52) is ikr_Xr1_tau in component ikr (ms).
C ALGBRC(9) is ikr_Xr2_a in component ikr (mS_per_uF).
C ALGBRC(25) is ikr_Xr2_b in component ikr (mS_per_uF).
C ALGBRC(39) is ikr_Xr2_inf in component ikr (dimensionless).
C ALGBRC(53) is ikr_Xr2_tau in component ikr (ms).
C CONSTS(47) is tau_1_offset in component ikr (ms).
C CONSTS(48) is tau_2_offset in component ikr (ms).
C STATES(14) is Xs in component iks (dimensionless).
C CONSTS(49) is g_Ks in component iks (mS_per_uF).
C ALGBRC(81) is i_Ks in component iks (A_per_F).
C ALGBRC(10) is iks_Xs_a in component iks (mS_per_uF).
C ALGBRC(26) is iks_Xs_b in component iks (mS_per_uF).
C ALGBRC(40) is iks_Xs_inf in component iks (dimensionless).
C ALGBRC(54) is iks_Xs_tau in component iks (ms).
C CONSTS(50) is ks1 in component iks (mS_per_uF).
C CONSTS(51) is ks2 in component iks (mV).
C CONSTS(104) is ks3 in component iks (mS_per_uF).
C CONSTS(137) is ks4 in component iks (mV).
C CONSTS(52) is ks5 in component iks (dimensionless).
C CONSTS(105) is ks6 in component iks (mV).
C CONSTS(53) is tauks_const in component iks (ms).
C CONSTS(106) is V_leak in component ileak (mS_per_uF).
C CONSTS(54) is g_Na in component ina (mS_per_uF).
C STATES(15) is h in component ina (dimensionless).
C CONSTS(55) is h1 in component ina (mS_per_uF).
C CONSTS(107) is h2 in component ina (mV).
C CONSTS(108) is h3 in component ina (mS_per_uF).
C CONSTS(138) is h4 in component ina (mV).
C CONSTS(56) is h5 in component ina (dimensionless).
C CONSTS(57) is h6 in component ina (mV).
C ALGBRC(83) is i_Na in component ina (A_per_F).
C ALGBRC(11) is ina_h_a in component ina (mS_per_uF).
C ALGBRC(27) is ina_h_b in component ina (mS_per_uF).
C ALGBRC(41) is ina_h_inf in component ina (dimensionless).
C ALGBRC(55) is ina_h_tau in component ina (ms).
C ALGBRC(12) is ina_j_a in component ina (mS_per_uF).
C ALGBRC(28) is ina_j_b in component ina (mS_per_uF).
C ALGBRC(42) is ina_j_inf in component ina (dimensionless).
C ALGBRC(56) is ina_j_tau in component ina (ms).
C ALGBRC(13) is ina_m_a in component ina (mS_per_uF).
C ALGBRC(29) is ina_m_b in component ina (mS_per_uF).
C ALGBRC(43) is ina_m_inf in component ina (dimensionless).
C ALGBRC(57) is ina_m_tau in component ina (ms).
C STATES(16) is j in component ina (dimensionless).
C CONSTS(58) is j1 in component ina (mS_per_uF).
C CONSTS(109) is j2 in component ina (mV).
C CONSTS(139) is j3 in component ina (mS_per_uF).
C CONSTS(149) is j4 in component ina (mV).
C CONSTS(110) is j5 in component ina (dimensionless).
C CONSTS(140) is j6 in component ina (mV).
C STATES(17) is m in component ina (dimensionless).
C CONSTS(59) is m1 in component ina (mS_per_uF).
C CONSTS(60) is m2 in component ina (mV).
C CONSTS(111) is m3 in component ina (mS_per_uF).
C CONSTS(141) is m4 in component ina (mV).
C CONSTS(61) is m5 in component ina (dimensionless).
C CONSTS(112) is m6 in component ina (mV).
C CONSTS(62) is tau_h_const in component ina (ms).
C CONSTS(63) is tau_j_const in component ina (ms).
C CONSTS(64) is tau_m_const in component ina (ms).
C CONSTS(65) is KmCa in component inaca (mM).
C CONSTS(66) is KmNai in component inaca (mM).
C CONSTS(67) is Ksat in component inaca (dimensionless).
C CONSTS(113) is alpha in component inaca (dimensionless).
C CONSTS(114) is gamma in component inaca (dimensionless).
C CONSTS(115) is kNaCa in component inaca (A_per_F).
C CONSTS(68) is Km_K in component inak (mM).
C CONSTS(69) is Km_Na in component inak (mM).
C CONSTS(116) is PNaK in component inak (A_per_F).
C ALGBRC(85) is i_NaK in component inak (A_per_F).
C CONSTS(70) is KPCa in component ipca (mM).
C CONSTS(117) is g_PCa in component ipca (A_per_F).
C STATES(18) is I in component irel (dimensionless).
C CONSTS(71) is MaxSR in component irel (dimensionless).
C CONSTS(72) is MinSR in component irel (dimensionless).
C STATES(19) is O in component irel (dimensionless).
C STATES(20) is R in component irel (dimensionless).
C ALGBRC(14) is RI in component irel (dimensionless).
C CONSTS(73) is ec50SR in component irel (mM).
C ALGBRC(30) is kCaSR in component irel (dimensionless).
C CONSTS(118) is kiCa in component irel (per_mM_per_ms).
C ALGBRC(44) is kiSRCa in component irel (per_mM_per_ms).
C CONSTS(119) is kim in component irel (mS_per_uF).
C CONSTS(142) is koCa in component irel (per_mM2_per_ms).
C ALGBRC(58) is koSRCa in component irel (per_mM2_per_ms).
C CONSTS(120) is kom in component irel (mS_per_uF).
C CONSTS(74) is ks in component irel (mS_per_uF).
C CONSTS(75) is g_to in component ito (mS_per_uF).
C ALGBRC(88) is i_to in component ito (A_per_F).
C ALGBRC(15) is ito_r_a in component ito (mS_per_uF).
C ALGBRC(31) is ito_r_b in component ito (mS_per_uF).
C ALGBRC(45) is ito_r_inf in component ito (dimensionless).
C ALGBRC(59) is ito_r_tau in component ito (ms).
C ALGBRC(16) is ito_s_a in component ito (mS_per_uF).
C ALGBRC(32) is ito_s_b in component ito (mS_per_uF).
C ALGBRC(46) is ito_s_inf in component ito (dimensionless).
C ALGBRC(60) is ito_s_tau in component ito (ms).
C STATES(21) is r in component ito (dimensionless).
C CONSTS(76) is r1 in component ito (mS_per_uF).
C CONSTS(77) is r2 in component ito (mV).
C CONSTS(121) is r3 in component ito (mS_per_uF).
C CONSTS(143) is r4 in component ito (mV).
C CONSTS(78) is r5 in component ito (dimensionless).
C CONSTS(122) is r6 in component ito (mV).
C STATES(22) is s in component ito (dimensionless).
C CONSTS(79) is s1 in component ito (mS_per_uF).
C CONSTS(123) is s2 in component ito (mV).
C CONSTS(124) is s3 in component ito (mS_per_uF).
C CONSTS(144) is s4 in component ito (mV).
C CONSTS(80) is s5 in component ito (dimensionless).
C CONSTS(81) is s6 in component ito (mV).
C CONSTS(82) is tau_r_const in component ito (ms).
C CONSTS(83) is tau_s_const in component ito (ms).
C CONSTS(125) is Kup in component iup (mM).
C CONSTS(126) is VmaxUp in component iup (mM_per_ms).
C CONSTS(151) is i_stim in component stimulus (A_per_F).
C CONSTS(84) is R in component phys (J_per_mol_per_K).
C CONSTS(85) is T in component phys (kelvin).
C CONSTS(150) is amplitude in component stimulus (A_per_F).
C RATES(1) is d/dt Cai in component cai (mM).
C RATES(2) is d/dt Ca_SR in component casr (mM).
C RATES(6) is d/dt d in component ical (dimensionless).
C RATES(7) is d/dt f in component ical (dimensionless).
C RATES(8) is d/dt fCa in component ical (dimensionless).
C RATES(9) is d/dt d in component icat (dimensionless).
C RATES(10) is d/dt f in component icat (dimensionless).
C RATES(11) is d/dt Xf in component ifunny (dimensionless).
C RATES(12) is d/dt Xr1 in component ikr (dimensionless).
C RATES(13) is d/dt Xr2 in component ikr (dimensionless).
C RATES(14) is d/dt Xs in component iks (dimensionless).
C RATES(15) is d/dt h in component ina (dimensionless).
C RATES(16) is d/dt j in component ina (dimensionless).
C RATES(17) is d/dt m in component ina (dimensionless).
C RATES(18) is d/dt I in component irel (dimensionless).
C RATES(19) is d/dt O in component irel (dimensionless).
C RATES(20) is d/dt R in component irel (dimensionless).
C RATES(21) is d/dt r in component ito (dimensionless).
C RATES(22) is d/dt s in component ito (dimensionless).
C RATES(3) is d/dt Ki in component ki (mM).
C RATES(5) is d/dt V in component membrane (mV).
C RATES(4) is d/dt Nai in component nai (mM).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.06
      STATES(1) = 2.19191642424964390e-04
      CONSTS(2) = 60.0
      CONSTS(3) = 9.64853415000000041e+01
      CONSTS(4) = 0.0006
      STATES(2) = 3.35086796732326109e-01
      CONSTS(5) = 0.3
      CONSTS(6) = 0.0
      CONSTS(7) = 1.8
      STATES(3) = 1.04748824394112106e+02
      CONSTS(8) = 5.4
      STATES(4) = 7.16928091250999167
      CONSTS(9) = 140.0
      CONSTS(10) = 1094.0
      CONSTS(11) = 3960.0
      CONSTS(12) = 16404.0
      STATES(5) = -7.55966016388546791e+01
      STATES(6) = 3.94925342652924281e-04
      CONSTS(13) = 1.29662941897219994e+01
      CONSTS(14) = 7.07914596471100044
      CONSTS(15) = 4.49094155069999987e-02
      STATES(7) = 1.70990105585540286e-01
      CONSTS(16) = 5.12589825999999987e-04
      CONSTS(17) = 1.93121122351431995e+03
      CONSTS(18) = 5.73002749969900016
      STATES(8) = 8.77798946134088598e-01
      CONSTS(19) = 3.08027691378999990e-01
      CONSTS(20) = 0.00054
      CONSTS(21) = 2.7e-07
      CONSTS(22) = 1.5e-08
      CONSTS(23) = 1.2
      CONSTS(24) = 2.0
      CONSTS(25) = 1.65824694683000007
      CONSTS(26) = 1.00462559171102995e+02
      STATES(9) = 2.70195573471577175e-04
      STATES(10) = 7.56032904368393432e-01
      CONSTS(27) = 0.185
      CONSTS(28) = 0.491
      STATES(11) = 6.40338504912615469e-03
      CONSTS(29) = 0.0435
      CONSTS(30) = 5.78970000000000002e-07
      CONSTS(31) = 2.00866502378844016e+04
      CONSTS(32) = 1.02023528452800001e+01
      CONSTS(33) = 2.39452913465299986e+01
      CONSTS(34) = 1.33785777797606004e-01
      CONSTS(35) = 4.77994972217041014e-01
      CONSTS(36) = 2.72427558793486995e+01
      CONSTS(37) = 4.92502331781412028
      CONSTS(38) = 8.72223760006881932
      CONSTS(39) = 5.66361974998243980e+01
      STATES(12) = 3.09767485715433222e-01
      CONSTS(40) = 5.74885237435000026e-03
      CONSTS(41) = 1.36234926362576001e+01
      CONSTS(42) = 4.76305711818360011e-02
      STATES(13) = 4.50577185148518577e-01
      CONSTS(43) = 1.24566405268270002e-02
      CONSTS(44) = 3.73426331501040991e+01
      CONSTS(45) = 2.20919642353902006e+01
      CONSTS(46) = 0.218025
      CONSTS(47) = 50.0
      CONSTS(48) = 0.0
      STATES(14) = 1.53788281650948710e-01
      CONSTS(49) = 0.0077
      CONSTS(50) = 1.16558447999999992e-03
      CONSTS(51) = 6.67268386758935958e+04
      CONSTS(52) = 2.80458908250000027e-01
      CONSTS(53) = 4.74115000000000034e-06
      CONSTS(54) = 9.72061340924100037
      STATES(15) = 7.39543607812429227e-01
      CONSTS(55) = 3.62659886399999999e-03
      CONSTS(56) = 9.66329497711473959e+03
      CONSTS(57) = 7.39550356461299963
      STATES(16) = 1.24515982574504899e-01
      CONSTS(58) = 5.12257182000000044e-04
      STATES(17) = 2.97549962926413614e-02
      CONSTS(59) = 1.08045846384818006e+02
      CONSTS(60) = 1.31070157339409992e+01
      CONSTS(61) = 2.32691436700000007e-03
      CONSTS(62) = 1.67331502516000014e-01
      CONSTS(63) = 9.51088724962000032e-01
      CONSTS(64) = 3.19775803839999970e-02
      CONSTS(65) = 1.38
      CONSTS(66) = 87.5
      CONSTS(67) = 0.1
      CONSTS(68) = 1.0
      CONSTS(69) = 40.0
      CONSTS(70) = 0.0005
      STATES(18) = 1.42153622323011597e-02
      CONSTS(71) = 15.0
      CONSTS(72) = 1.0
      STATES(19) = 1.65045105312396393e-04
      STATES(20) = 1.13120363433751106e-02
      CONSTS(73) = 0.45
      CONSTS(74) = 12.5
      CONSTS(75) = 1.17833333333299997e-01
      STATES(21) = 2.67597833344160611e-04
      CONSTS(76) = 5.53614181712999975e-02
      CONSTS(77) = 1.16842023429669002e+01
      CONSTS(78) = 3.98918108037750008
      STATES(22) = 7.46802810614006107e-01
      CONSTS(79) = 3.44230944300000013e-04
      CONSTS(80) = 1.86760536909694991e+02
      CONSTS(81) = 8.18093387332270083
      CONSTS(82) = 6.96758421171499998e-01
      CONSTS(83) = 1.12244577239468999e+01
      CONSTS(84) = 8.314472
      CONSTS(85) = 310.0
      CONSTS(86) = CONSTS(12)+CONSTS(10)
      CONSTS(87) =  10.0000*1.20000
      CONSTS(88) = ( CONSTS(84)*CONSTS(85))/CONSTS(3)
      CONSTS(89) =  0.000592000*0.620000
      CONSTS(90) =  0.000290000*1.50000
      CONSTS(91) = CONSTS(3)/( CONSTS(84)*CONSTS(85))
      CONSTS(92) =  CONSTS(15)*CONSTS(13)
      CONSTS(93) = - 6.90988
      CONSTS(94) = - 49.5057
      CONSTS(95) =  CONSTS(17)*CONSTS(16)
      CONSTS(96) = (CONSTS(20)+CONSTS(22))+CONSTS(21)
      CONSTS(97) = CONSTS(28)/(CONSTS(28)+1.00000)
      CONSTS(98) = - 14.5897
      CONSTS(99) =  CONSTS(31)*CONSTS(30)
      CONSTS(100) =  CONSTS(42)*CONSTS(40)
      CONSTS(101) = - 7.06809
      CONSTS(102) = - 25.9945
      CONSTS(103) =  CONSTS(44)*CONSTS(43)
      CONSTS(104) =  CONSTS(52)*CONSTS(50)
      CONSTS(105) = - 18.8670
      CONSTS(106) =  8.00000e-05*0.0200000
      CONSTS(107) = - 19.8394
      CONSTS(108) =  CONSTS(56)*CONSTS(55)
      CONSTS(109) = - 66.5838
      CONSTS(110) = CONSTS(56)
      CONSTS(111) =  CONSTS(61)*CONSTS(59)
      CONSTS(112) = - 7.91773
      CONSTS(113) =  2.50000*1.10000
      CONSTS(114) =  0.350000*2.00000
      CONSTS(115) =  1000.00*1.10000
      CONSTS(116) =  1.36200*1.81800
      CONSTS(117) =  0.0250000*10.5000
      CONSTS(118) =  54.0000*0.342500
      CONSTS(119) =  0.00100000*0.557100
      CONSTS(120) =  1.50000*0.142900
      CONSTS(121) =  CONSTS(78)*CONSTS(76)
      CONSTS(122) = - 11.0471
      CONSTS(123) = - 17.6345
      CONSTS(124) =  CONSTS(80)*CONSTS(79)
      CONSTS(125) =  0.000250000*0.702000
      CONSTS(126) =  0.000425000*0.260000
      CONSTS(127) =  CONSTS(11)*(CONSTS(10)/CONSTS(86))
      CONSTS(128) =  CONSTS(3)*CONSTS(91)
      CONSTS(129) = 1.00000/(1.00000/CONSTS(14)+1.00000/CONSTS(93))
      CONSTS(130) = 1.00000/(1.00000/CONSTS(94)+1.00000/CONSTS(18))
      CONSTS(131) = CONSTS(20)/CONSTS(96)
      CONSTS(132) = CONSTS(21)/CONSTS(96)
      CONSTS(133) = CONSTS(22)/CONSTS(96)
      CONSTS(134) = 1.00000/(1.00000/CONSTS(98)+1.00000/CONSTS(32))
      CONSTS(135) = 1.00000/(1.00000/CONSTS(41)+1.00000/CONSTS(101))
      CONSTS(136) = 1.00000/(1.00000/CONSTS(102)+1.00000/CONSTS(45))
      CONSTS(137) = 1.00000/(1.00000/CONSTS(51)+1.00000/CONSTS(105))
      CONSTS(138) = 1.00000/(1.00000/CONSTS(107)+1.00000/CONSTS(57))
      CONSTS(139) =  CONSTS(110)*CONSTS(58)
      CONSTS(140) = CONSTS(57)
      CONSTS(141) = 1.00000/(1.00000/CONSTS(60)+1.00000/CONSTS(112))
      CONSTS(142) =  56320.0*11.4302
      CONSTS(143) = 1.00000/(1.00000/CONSTS(77)+1.00000/CONSTS(122))
      CONSTS(144) = 1.00000/(1.00000/CONSTS(123)+1.00000/CONSTS(81))
      CONSTS(145) =  CONSTS(11)*(CONSTS(12)/CONSTS(86))
      CONSTS(146) =  CONSTS(131)*CONSTS(19)
      CONSTS(147) =  CONSTS(132)*CONSTS(19)
      CONSTS(148) =  CONSTS(133)*CONSTS(19)
      CONSTS(149) = 1.00000/(1.00000/CONSTS(109)+1.00000/CONSTS(140))
      CONSTS(150) = - 3.00000
      CONSTS(151) =  CONSTS(6)*CONSTS(150)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(5)+26.3000)/- 6.00000))
      ALGBRC(21) = 1.00000/( 1.06800*EXP((STATES(5)+26.3000)/30.0000)+ 1.06800*EXP((STATES(5)+26.3000)/- 30.0000))
      RATES(9) = (ALGBRC(5) - STATES(9))/ALGBRC(21)
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(5)+61.7000)/5.60000))
      ALGBRC(22) = 1.00000/( 0.0153000*EXP(- (STATES(5)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(5)+61.7000)/15.3800))
      RATES(10) = (ALGBRC(6) - STATES(10))/ALGBRC(22)
      ALGBRC(2) =  CONSTS(13)*EXP(STATES(5)/CONSTS(14))
      ALGBRC(18) =  CONSTS(92)*EXP(STATES(5)/CONSTS(129))
      ALGBRC(34) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(18))
      ALGBRC(48) = 1.00000/(ALGBRC(2)+ALGBRC(18))+CONSTS(25)
      RATES(6) = (ALGBRC(34) - STATES(6))/ALGBRC(48)
      ALGBRC(3) =  CONSTS(16)*EXP(STATES(5)/CONSTS(94))
      ALGBRC(19) =  CONSTS(95)*EXP(STATES(5)/CONSTS(130))
      ALGBRC(35) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(19))
      ALGBRC(49) = 1.00000/(ALGBRC(3)+ALGBRC(19))+CONSTS(26)
      RATES(7) = (ALGBRC(35) - STATES(7))/ALGBRC(49)
      ALGBRC(7) =  CONSTS(30)*EXP(STATES(5)/CONSTS(98))
      ALGBRC(23) =  CONSTS(99)*EXP(STATES(5)/CONSTS(134))
      ALGBRC(37) = ALGBRC(7)/(ALGBRC(7)+ALGBRC(23))
      ALGBRC(51) = 1.00000/(ALGBRC(7)+ALGBRC(23))+CONSTS(33)
      RATES(11) = (ALGBRC(37) - STATES(11))/ALGBRC(51)
      ALGBRC(8) =  CONSTS(40)*EXP(STATES(5)/CONSTS(41))
      ALGBRC(24) =  CONSTS(100)*EXP(STATES(5)/CONSTS(135))
      ALGBRC(38) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(24))
      ALGBRC(52) = 1.00000/(ALGBRC(8)+ALGBRC(24))+CONSTS(47)
      RATES(12) = (ALGBRC(38) - STATES(12))/ALGBRC(52)
      ALGBRC(9) =  CONSTS(43)*EXP(STATES(5)/CONSTS(102))
      ALGBRC(25) =  CONSTS(103)*EXP(STATES(5)/CONSTS(136))
      ALGBRC(39) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(25))
      ALGBRC(53) = 1.00000/(ALGBRC(9)+ALGBRC(25))+CONSTS(48)
      RATES(13) = (ALGBRC(39) - STATES(13))/ALGBRC(53)
      ALGBRC(10) =  CONSTS(50)*EXP(STATES(5)/CONSTS(51))
      ALGBRC(26) =  CONSTS(104)*EXP(STATES(5)/CONSTS(137))
      ALGBRC(40) = ALGBRC(10)/(ALGBRC(10)+ALGBRC(26))
      ALGBRC(54) = 1.00000/(ALGBRC(10)+ALGBRC(26))+CONSTS(53)
      RATES(14) = (ALGBRC(40) - STATES(14))/ALGBRC(54)
      ALGBRC(11) =  CONSTS(55)*EXP(STATES(5)/CONSTS(107))
      ALGBRC(27) =  CONSTS(108)*EXP(STATES(5)/CONSTS(138))
      ALGBRC(41) = ALGBRC(11)/(ALGBRC(11)+ALGBRC(27))
      ALGBRC(55) = 1.00000/(ALGBRC(11)+ALGBRC(27))+CONSTS(62)
      RATES(15) = (ALGBRC(41) - STATES(15))/ALGBRC(55)
      ALGBRC(12) =  CONSTS(58)*EXP(STATES(5)/CONSTS(109))
      ALGBRC(28) =  CONSTS(139)*EXP(STATES(5)/CONSTS(149))
      ALGBRC(42) = ALGBRC(12)/(ALGBRC(12)+ALGBRC(28))
      ALGBRC(56) = 1.00000/(ALGBRC(12)+ALGBRC(28))+CONSTS(63)
      RATES(16) = (ALGBRC(42) - STATES(16))/ALGBRC(56)
      ALGBRC(13) =  CONSTS(59)*EXP(STATES(5)/CONSTS(60))
      ALGBRC(29) =  CONSTS(111)*EXP(STATES(5)/CONSTS(141))
      ALGBRC(43) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(29))
      ALGBRC(57) = 1.00000/(ALGBRC(13)+ALGBRC(29))+CONSTS(64)
      RATES(17) = (ALGBRC(43) - STATES(17))/ALGBRC(57)
      ALGBRC(14) = ((1.00000 - STATES(20)) - STATES(19)) - STATES(18)
      ALGBRC(30) = CONSTS(71) - (CONSTS(71) - CONSTS(72))/(1.00000+CONSTS(73)/STATES(2) ** 2.50000)
      ALGBRC(44) =  CONSTS(118)*ALGBRC(30)
      ALGBRC(58) = CONSTS(142)/ALGBRC(30)
      RATES(18) = (( ( ALGBRC(44)*STATES(1))*STATES(19) -  CONSTS(119)*STATES(18)) -  CONSTS(120)*STATES(18))+ ( ALGBRC(58)*STATES(1) ** 2.00000)*ALGBRC(14)
      RATES(19) = (( ( ALGBRC(58)*STATES(1) ** 2.00000)*STATES(20) -  CONSTS(120)*STATES(19)) -  ( ALGBRC(44)*STATES(1))*STATES(19))+ CONSTS(119)*STATES(18)
      RATES(20) = (( CONSTS(119)*ALGBRC(14) -  ( ALGBRC(44)*STATES(1))*STATES(20)) -  ( ALGBRC(58)*STATES(1) ** 2.00000)*STATES(20))+ CONSTS(120)*STATES(19)
      ALGBRC(15) =  CONSTS(76)*EXP(STATES(5)/CONSTS(77))
      ALGBRC(31) =  CONSTS(121)*EXP(STATES(5)/CONSTS(143))
      ALGBRC(45) = ALGBRC(15)/(ALGBRC(15)+ALGBRC(31))
      ALGBRC(59) = 1.00000/(ALGBRC(15)+ALGBRC(31))+CONSTS(82)
      RATES(21) = (ALGBRC(45) - STATES(21))/ALGBRC(59)
      ALGBRC(16) =  CONSTS(79)*EXP(STATES(5)/CONSTS(123))
      ALGBRC(32) =  CONSTS(124)*EXP(STATES(5)/CONSTS(144))
      ALGBRC(46) = ALGBRC(16)/(ALGBRC(16)+ALGBRC(32))
      ALGBRC(60) = 1.00000/(ALGBRC(16)+ALGBRC(32))+CONSTS(83)
      RATES(22) = (ALGBRC(46) - STATES(22))/ALGBRC(60)
      ALGBRC(4) = 1.00000/(1.00000+( CONSTS(23)*STATES(1))/0.000325000 ** 8.00000)
      ALGBRC(20) = 0.100000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000500000)/0.000100000))
      ALGBRC(36) = 0.200000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000750000)/0.000800000))
      ALGBRC(50) = (((ALGBRC(4)+ALGBRC(20))+ALGBRC(36))+0.230000)/1.46000
      ALGBRC(62) = TERNRY(ALGBRC(50).GT.STATES(8).AND.STATES(5).GT.- 60.0000, 0.00000, 1.00000)
      RATES(8) = ( ALGBRC(62)*(ALGBRC(50) - STATES(8)))/CONSTS(24)
      ALGBRC(84) = ( CONSTS(115)*( ( EXP( ( CONSTS(114)*STATES(5))*CONSTS(91))*STATES(4) ** 3.00000)*CONSTS(7) -  ( ( EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))*CONSTS(9) ** 3.00000)*STATES(1))*CONSTS(113)))/( ( (CONSTS(66) ** 3.00000+CONSTS(9) ** 3.00000)*(CONSTS(65)+CONSTS(7)))*(1.00000+ CONSTS(67)*EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))))
      ALGBRC(61) =  CONSTS(88)*log(CONSTS(9)/STATES(4))
      ALGBRC(64) =  CONSTS(90)*(STATES(5) - ALGBRC(61))
      ALGBRC(69) = ( ( ( CONSTS(148)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(4))*EXP( STATES(5)*CONSTS(91)) -  0.750000*CONSTS(9)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
      ALGBRC(70) =  ( ( ALGBRC(69)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(74) =  ( ( CONSTS(97)*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(61))
      ALGBRC(83) =  ( ( ( CONSTS(54)*STATES(17) ** 3.00000)*STATES(15))*STATES(16))*(STATES(5) - ALGBRC(61))
      ALGBRC(85) = ( ( CONSTS(116)*CONSTS(8))*STATES(4))/( ( (CONSTS(8)+CONSTS(68))*(STATES(4)+CONSTS(69)))*((1.00000+ 0.124500*EXP( ( - 0.100000*STATES(5))*CONSTS(91)))+ 0.0353000*EXP( - STATES(5)*CONSTS(91))))
      RATES(4) =  (- CONSTS(2)/( CONSTS(3)*CONSTS(145)))*((((((ALGBRC(83))+ALGBRC(64))+ALGBRC(74))+ 3.00000*ALGBRC(85))+ 3.00000*ALGBRC(84))+ALGBRC(70))
      ALGBRC(67) = ( ( ( CONSTS(147)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(3))*EXP( STATES(5)*CONSTS(91)) -  0.750000*CONSTS(8)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
      ALGBRC(68) =  ( ( ALGBRC(67)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(47) =  CONSTS(88)*log(CONSTS(8)/STATES(3))
      ALGBRC(73) =  ( ( (1.00000 - CONSTS(97))*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(47))
      ALGBRC(76) =  CONSTS(35)*EXP((STATES(5)+CONSTS(37))/CONSTS(36))
      ALGBRC(77) =  1.00000*EXP((STATES(5)+CONSTS(39))/CONSTS(38))
      ALGBRC(78) = ALGBRC(76)/(ALGBRC(76)+ALGBRC(77))
      ALGBRC(79) =  ( ( CONSTS(34)* (CONSTS(8)/5.40000) ** (1.0 / 2))*ALGBRC(78))*(STATES(5) - ALGBRC(47))
      ALGBRC(80) =  ( ( ( CONSTS(46)* (CONSTS(8)/5.40000) ** (1.0 / 2))*STATES(12))*STATES(13))*(STATES(5) - ALGBRC(47))
      ALGBRC(81) =  ( CONSTS(49)*STATES(14) ** 2.00000)*(STATES(5) - ALGBRC(47))
      ALGBRC(88) =  ( ( CONSTS(75)*STATES(21))*STATES(22))*(STATES(5) - ALGBRC(47))
      RATES(3) =  (- CONSTS(2)/( CONSTS(3)*CONSTS(145)))*(((((((ALGBRC(79))+ALGBRC(88))+ALGBRC(80))+ALGBRC(81))+ALGBRC(73)) -  2.00000*ALGBRC(85))+ALGBRC(68))
      ALGBRC(33) =  ( 0.500000*CONSTS(88))*log(CONSTS(7)/STATES(1))
      ALGBRC(72) =  ( ( CONSTS(27)*STATES(9))*STATES(10))*(STATES(5) - ALGBRC(33))
      ALGBRC(86) = ( CONSTS(117)*STATES(1))/(STATES(1)+CONSTS(70))
      ALGBRC(63) =  CONSTS(89)*(STATES(5) - ALGBRC(33))
      ALGBRC(65) = ( ( ( ( CONSTS(146)*4.00000)*STATES(5))*CONSTS(128))*( ( 0.341000*STATES(1))*EXP( ( 2.00000*STATES(5))*CONSTS(91)) -  0.341000*CONSTS(7)))/(EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 1.00000)
      ALGBRC(66) =  ( ( ALGBRC(65)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(71) = (ALGBRC(66)+ALGBRC(70))+ALGBRC(68)
      ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
      RATES(5) = - ((((((((((((((ALGBRC(79))+ALGBRC(88))+ALGBRC(80))+ALGBRC(81))+ALGBRC(71))+ALGBRC(72))+ALGBRC(85))+ALGBRC(83))+ALGBRC(84))+ALGBRC(86))+ALGBRC(75))+ALGBRC(64))+ALGBRC(63))+CONSTS(151))
      ALGBRC(1) = 1.00000/(1.00000+( CONSTS(1)*CONSTS(4))/STATES(1)+CONSTS(4) ** 2.00000)
      ALGBRC(82) =  (STATES(2) - STATES(1))*CONSTS(106)
      ALGBRC(87) =  ( ( CONSTS(74)*STATES(19))*(STATES(2) - STATES(1)))*(CONSTS(127)/CONSTS(145))
      ALGBRC(89) = CONSTS(126)/(1.00000+CONSTS(125) ** 2.00000/STATES(1) ** 2.00000)
      RATES(1) =  ALGBRC(1)*(((- ALGBRC(89)+ALGBRC(82))+ALGBRC(87)) -  (CONSTS(2)/( ( 2.00000*CONSTS(145))*CONSTS(3)))*(((((ALGBRC(66))+ALGBRC(72))+ALGBRC(63))+ALGBRC(86)) -  2.00000*ALGBRC(84)))
      ALGBRC(17) = 1.00000/(1.00000+( CONSTS(87)*CONSTS(5))/STATES(2)+CONSTS(5) ** 2.00000)
      RATES(2) =  (( ALGBRC(17)*CONSTS(145))/CONSTS(127))*((ALGBRC(89) - ALGBRC(87)) - ALGBRC(82))
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(5)+26.3000)/- 6.00000))
      ALGBRC(21) = 1.00000/( 1.06800*EXP((STATES(5)+26.3000)/30.0000)+ 1.06800*EXP((STATES(5)+26.3000)/- 30.0000))
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(5)+61.7000)/5.60000))
      ALGBRC(22) = 1.00000/( 0.0153000*EXP(- (STATES(5)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(5)+61.7000)/15.3800))
      ALGBRC(2) =  CONSTS(13)*EXP(STATES(5)/CONSTS(14))
      ALGBRC(18) =  CONSTS(92)*EXP(STATES(5)/CONSTS(129))
      ALGBRC(34) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(18))
      ALGBRC(48) = 1.00000/(ALGBRC(2)+ALGBRC(18))+CONSTS(25)
      ALGBRC(3) =  CONSTS(16)*EXP(STATES(5)/CONSTS(94))
      ALGBRC(19) =  CONSTS(95)*EXP(STATES(5)/CONSTS(130))
      ALGBRC(35) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(19))
      ALGBRC(49) = 1.00000/(ALGBRC(3)+ALGBRC(19))+CONSTS(26)
      ALGBRC(7) =  CONSTS(30)*EXP(STATES(5)/CONSTS(98))
      ALGBRC(23) =  CONSTS(99)*EXP(STATES(5)/CONSTS(134))
      ALGBRC(37) = ALGBRC(7)/(ALGBRC(7)+ALGBRC(23))
      ALGBRC(51) = 1.00000/(ALGBRC(7)+ALGBRC(23))+CONSTS(33)
      ALGBRC(8) =  CONSTS(40)*EXP(STATES(5)/CONSTS(41))
      ALGBRC(24) =  CONSTS(100)*EXP(STATES(5)/CONSTS(135))
      ALGBRC(38) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(24))
      ALGBRC(52) = 1.00000/(ALGBRC(8)+ALGBRC(24))+CONSTS(47)
      ALGBRC(9) =  CONSTS(43)*EXP(STATES(5)/CONSTS(102))
      ALGBRC(25) =  CONSTS(103)*EXP(STATES(5)/CONSTS(136))
      ALGBRC(39) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(25))
      ALGBRC(53) = 1.00000/(ALGBRC(9)+ALGBRC(25))+CONSTS(48)
      ALGBRC(10) =  CONSTS(50)*EXP(STATES(5)/CONSTS(51))
      ALGBRC(26) =  CONSTS(104)*EXP(STATES(5)/CONSTS(137))
      ALGBRC(40) = ALGBRC(10)/(ALGBRC(10)+ALGBRC(26))
      ALGBRC(54) = 1.00000/(ALGBRC(10)+ALGBRC(26))+CONSTS(53)
      ALGBRC(11) =  CONSTS(55)*EXP(STATES(5)/CONSTS(107))
      ALGBRC(27) =  CONSTS(108)*EXP(STATES(5)/CONSTS(138))
      ALGBRC(41) = ALGBRC(11)/(ALGBRC(11)+ALGBRC(27))
      ALGBRC(55) = 1.00000/(ALGBRC(11)+ALGBRC(27))+CONSTS(62)
      ALGBRC(12) =  CONSTS(58)*EXP(STATES(5)/CONSTS(109))
      ALGBRC(28) =  CONSTS(139)*EXP(STATES(5)/CONSTS(149))
      ALGBRC(42) = ALGBRC(12)/(ALGBRC(12)+ALGBRC(28))
      ALGBRC(56) = 1.00000/(ALGBRC(12)+ALGBRC(28))+CONSTS(63)
      ALGBRC(13) =  CONSTS(59)*EXP(STATES(5)/CONSTS(60))
      ALGBRC(29) =  CONSTS(111)*EXP(STATES(5)/CONSTS(141))
      ALGBRC(43) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(29))
      ALGBRC(57) = 1.00000/(ALGBRC(13)+ALGBRC(29))+CONSTS(64)
      ALGBRC(14) = ((1.00000 - STATES(20)) - STATES(19)) - STATES(18)
      ALGBRC(30) = CONSTS(71) - (CONSTS(71) - CONSTS(72))/(1.00000+CONSTS(73)/STATES(2) ** 2.50000)
      ALGBRC(44) =  CONSTS(118)*ALGBRC(30)
      ALGBRC(58) = CONSTS(142)/ALGBRC(30)
      ALGBRC(15) =  CONSTS(76)*EXP(STATES(5)/CONSTS(77))
      ALGBRC(31) =  CONSTS(121)*EXP(STATES(5)/CONSTS(143))
      ALGBRC(45) = ALGBRC(15)/(ALGBRC(15)+ALGBRC(31))
      ALGBRC(59) = 1.00000/(ALGBRC(15)+ALGBRC(31))+CONSTS(82)
      ALGBRC(16) =  CONSTS(79)*EXP(STATES(5)/CONSTS(123))
      ALGBRC(32) =  CONSTS(124)*EXP(STATES(5)/CONSTS(144))
      ALGBRC(46) = ALGBRC(16)/(ALGBRC(16)+ALGBRC(32))
      ALGBRC(60) = 1.00000/(ALGBRC(16)+ALGBRC(32))+CONSTS(83)
      ALGBRC(4) = 1.00000/(1.00000+( CONSTS(23)*STATES(1))/0.000325000 ** 8.00000)
      ALGBRC(20) = 0.100000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000500000)/0.000100000))
      ALGBRC(36) = 0.200000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000750000)/0.000800000))
      ALGBRC(50) = (((ALGBRC(4)+ALGBRC(20))+ALGBRC(36))+0.230000)/1.46000
      ALGBRC(62) = TERNRY(ALGBRC(50).GT.STATES(8).AND.STATES(5).GT.- 60.0000, 0.00000, 1.00000)
      ALGBRC(84) = ( CONSTS(115)*( ( EXP( ( CONSTS(114)*STATES(5))*CONSTS(91))*STATES(4) ** 3.00000)*CONSTS(7) -  ( ( EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))*CONSTS(9) ** 3.00000)*STATES(1))*CONSTS(113)))/( ( (CONSTS(66) ** 3.00000+CONSTS(9) ** 3.00000)*(CONSTS(65)+CONSTS(7)))*(1.00000+ CONSTS(67)*EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))))
      ALGBRC(61) =  CONSTS(88)*log(CONSTS(9)/STATES(4))
      ALGBRC(64) =  CONSTS(90)*(STATES(5) - ALGBRC(61))
      ALGBRC(69) = ( ( ( CONSTS(148)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(4))*EXP( STATES(5)*CONSTS(91)) -  0.750000*CONSTS(9)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
      ALGBRC(70) =  ( ( ALGBRC(69)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(74) =  ( ( CONSTS(97)*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(61))
      ALGBRC(83) =  ( ( ( CONSTS(54)*STATES(17) ** 3.00000)*STATES(15))*STATES(16))*(STATES(5) - ALGBRC(61))
      ALGBRC(85) = ( ( CONSTS(116)*CONSTS(8))*STATES(4))/( ( (CONSTS(8)+CONSTS(68))*(STATES(4)+CONSTS(69)))*((1.00000+ 0.124500*EXP( ( - 0.100000*STATES(5))*CONSTS(91)))+ 0.0353000*EXP( - STATES(5)*CONSTS(91))))
      ALGBRC(67) = ( ( ( CONSTS(147)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(3))*EXP( STATES(5)*CONSTS(91)) -  0.750000*CONSTS(8)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
      ALGBRC(68) =  ( ( ALGBRC(67)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(47) =  CONSTS(88)*log(CONSTS(8)/STATES(3))
      ALGBRC(73) =  ( ( (1.00000 - CONSTS(97))*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(47))
      ALGBRC(76) =  CONSTS(35)*EXP((STATES(5)+CONSTS(37))/CONSTS(36))
      ALGBRC(77) =  1.00000*EXP((STATES(5)+CONSTS(39))/CONSTS(38))
      ALGBRC(78) = ALGBRC(76)/(ALGBRC(76)+ALGBRC(77))
      ALGBRC(79) =  ( ( CONSTS(34)* (CONSTS(8)/5.40000) ** (1.0 / 2))*ALGBRC(78))*(STATES(5) - ALGBRC(47))
      ALGBRC(80) =  ( ( ( CONSTS(46)* (CONSTS(8)/5.40000) ** (1.0 / 2))*STATES(12))*STATES(13))*(STATES(5) - ALGBRC(47))
      ALGBRC(81) =  ( CONSTS(49)*STATES(14) ** 2.00000)*(STATES(5) - ALGBRC(47))
      ALGBRC(88) =  ( ( CONSTS(75)*STATES(21))*STATES(22))*(STATES(5) - ALGBRC(47))
      ALGBRC(33) =  ( 0.500000*CONSTS(88))*log(CONSTS(7)/STATES(1))
      ALGBRC(72) =  ( ( CONSTS(27)*STATES(9))*STATES(10))*(STATES(5) - ALGBRC(33))
      ALGBRC(86) = ( CONSTS(117)*STATES(1))/(STATES(1)+CONSTS(70))
      ALGBRC(63) =  CONSTS(89)*(STATES(5) - ALGBRC(33))
      ALGBRC(65) = ( ( ( ( CONSTS(146)*4.00000)*STATES(5))*CONSTS(128))*( ( 0.341000*STATES(1))*EXP( ( 2.00000*STATES(5))*CONSTS(91)) -  0.341000*CONSTS(7)))/(EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 1.00000)
      ALGBRC(66) =  ( ( ALGBRC(65)*STATES(6))*STATES(7))*STATES(8)
      ALGBRC(71) = (ALGBRC(66)+ALGBRC(70))+ALGBRC(68)
      ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
      ALGBRC(1) = 1.00000/(1.00000+( CONSTS(1)*CONSTS(4))/STATES(1)+CONSTS(4) ** 2.00000)
      ALGBRC(82) =  (STATES(2) - STATES(1))*CONSTS(106)
      ALGBRC(87) =  ( ( CONSTS(74)*STATES(19))*(STATES(2) - STATES(1)))*(CONSTS(127)/CONSTS(145))
      ALGBRC(89) = CONSTS(126)/(1.00000+CONSTS(125) ** 2.00000/STATES(1) ** 2.00000)
      ALGBRC(17) = 1.00000/(1.00000+( CONSTS(87)*CONSTS(5))/STATES(2)+CONSTS(5) ** 2.00000)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END