C
C There are a total of 80 entries in the algebraic variable array.
C There are a total of 25 entries in each of the rate and state variable arrays.
C There are a total of 160 entries in the constant variable array.
C
C
C ALGBRC(1) is sum_g_syn_e_s_1 in component synaptic_coupling (nanoS).
C ALGBRC(6) is sum_g_syn_e_s_2 in component synaptic_coupling (nanoS).
C ALGBRC(11) is sum_g_syn_e_s_3 in component synaptic_coupling (nanoS).
C ALGBRC(16) is sum_g_syn_e_s_4 in component synaptic_coupling (nanoS).
C ALGBRC(21) is sum_g_syn_e_s_5 in component synaptic_coupling (nanoS).
C CONSTS(1) is g_syn_e_1_2 in component synaptic_coupling (nanoS).
C CONSTS(2) is g_syn_e_1_3 in component synaptic_coupling (nanoS).
C CONSTS(3) is g_syn_e_1_4 in component synaptic_coupling (nanoS).
C CONSTS(4) is g_syn_e_1_5 in component synaptic_coupling (nanoS).
C CONSTS(5) is g_syn_e_2_1 in component synaptic_coupling (nanoS).
C CONSTS(6) is g_syn_e_2_3 in component synaptic_coupling (nanoS).
C CONSTS(7) is g_syn_e_2_4 in component synaptic_coupling (nanoS).
C CONSTS(8) is g_syn_e_2_5 in component synaptic_coupling (nanoS).
C CONSTS(9) is g_syn_e_3_1 in component synaptic_coupling (nanoS).
C CONSTS(10) is g_syn_e_3_2 in component synaptic_coupling (nanoS).
C CONSTS(11) is g_syn_e_3_4 in component synaptic_coupling (nanoS).
C CONSTS(12) is g_syn_e_3_5 in component synaptic_coupling (nanoS).
C CONSTS(13) is g_syn_e_4_1 in component synaptic_coupling (nanoS).
C CONSTS(14) is g_syn_e_4_2 in component synaptic_coupling (nanoS).
C CONSTS(15) is g_syn_e_4_3 in component synaptic_coupling (nanoS).
C CONSTS(16) is g_syn_e_4_5 in component synaptic_coupling (nanoS).
C CONSTS(17) is g_syn_e_5_1 in component synaptic_coupling (nanoS).
C CONSTS(18) is g_syn_e_5_2 in component synaptic_coupling (nanoS).
C CONSTS(19) is g_syn_e_5_3 in component synaptic_coupling (nanoS).
C CONSTS(20) is g_syn_e_5_4 in component synaptic_coupling (nanoS).
C STATES(1) is s in component synaptic_input (dimensionless).
C STATES(2) is s in component synaptic_input (dimensionless).
C STATES(3) is s in component synaptic_input (dimensionless).
C STATES(4) is s in component synaptic_input (dimensionless).
C STATES(5) is s in component synaptic_input (dimensionless).
C VOI is time in component environment (millisecond).
C STATES(6) is V in component membrane (millivolt).
C CONSTS(21) is C in component membrane (picoF).
C CONSTS(22) is i_app in component membrane (picoA).
C ALGBRC(61) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(46) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(51) is i_K in component potassium_current (picoA).
C ALGBRC(66) is i_L in component leakage_current (picoA).
C ALGBRC(71) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(76) is i_syn_e in component synaptic_input (picoA).
C CONSTS(23) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(24) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(26) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(7) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(25) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(26) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(2) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(27) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(27) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(28) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(29) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(30) is g_K in component potassium_current (nanoS).
C CONSTS(31) is E_K in component potassium_current (millivolt).
C STATES(8) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(3) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(28) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(32) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(33) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(34) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(35) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(56) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(9) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(36) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(37) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(4) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(29) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(38) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(39) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(40) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(41) is g_L in component leakage_current (nanoS).
C CONSTS(42) is E_L in component leakage_current (millivolt).
C CONSTS(43) is E_syn_e in component tonic_current (millivolt).
C CONSTS(44) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(5) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(45) is kr in component synaptic_input (dimensionless).
C CONSTS(46) is tau_s in component synaptic_input (millisecond).
C CONSTS(47) is sigma_s in component synaptic_input (millivolt).
C CONSTS(48) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(10) is V in component membrane (millivolt).
C CONSTS(49) is C in component membrane (picoF).
C CONSTS(50) is i_app in component membrane (picoA).
C ALGBRC(62) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(47) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(52) is i_K in component potassium_current (picoA).
C ALGBRC(67) is i_L in component leakage_current (picoA).
C ALGBRC(72) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(77) is i_syn_e in component synaptic_input (picoA).
C CONSTS(51) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(52) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(30) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(11) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(53) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(54) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(7) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(31) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(55) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(56) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(57) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(58) is g_K in component potassium_current (nanoS).
C CONSTS(59) is E_K in component potassium_current (millivolt).
C STATES(12) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(8) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(32) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(60) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(61) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(62) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(63) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(57) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(13) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(64) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(65) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(9) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(33) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(66) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(67) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(68) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(69) is g_L in component leakage_current (nanoS).
C CONSTS(70) is E_L in component leakage_current (millivolt).
C CONSTS(71) is E_syn_e in component tonic_current (millivolt).
C CONSTS(72) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(10) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(73) is kr in component synaptic_input (dimensionless).
C CONSTS(74) is tau_s in component synaptic_input (millisecond).
C CONSTS(75) is sigma_s in component synaptic_input (millivolt).
C CONSTS(76) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(14) is V in component membrane (millivolt).
C CONSTS(77) is C in component membrane (picoF).
C CONSTS(78) is i_app in component membrane (picoA).
C ALGBRC(63) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(48) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(53) is i_K in component potassium_current (picoA).
C ALGBRC(68) is i_L in component leakage_current (picoA).
C ALGBRC(73) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(78) is i_syn_e in component synaptic_input (picoA).
C CONSTS(79) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(80) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(34) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(15) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(81) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(82) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(12) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(35) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(83) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(84) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(85) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(86) is g_K in component potassium_current (nanoS).
C CONSTS(87) is E_K in component potassium_current (millivolt).
C STATES(16) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(13) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(36) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(88) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(89) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(90) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(91) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(58) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(17) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(92) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(93) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(14) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(37) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(94) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(95) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(96) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(97) is g_L in component leakage_current (nanoS).
C CONSTS(98) is E_L in component leakage_current (millivolt).
C CONSTS(99) is E_syn_e in component tonic_current (millivolt).
C CONSTS(100) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(15) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(101) is kr in component synaptic_input (dimensionless).
C CONSTS(102) is tau_s in component synaptic_input (millisecond).
C CONSTS(103) is sigma_s in component synaptic_input (millivolt).
C CONSTS(104) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(18) is V in component membrane (millivolt).
C CONSTS(105) is C in component membrane (picoF).
C CONSTS(106) is i_app in component membrane (picoA).
C ALGBRC(64) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(49) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(54) is i_K in component potassium_current (picoA).
C ALGBRC(69) is i_L in component leakage_current (picoA).
C ALGBRC(74) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(79) is i_syn_e in component synaptic_input (picoA).
C CONSTS(107) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(108) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(38) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(19) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(109) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(110) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(17) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(39) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(111) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(112) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(113) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(114) is g_K in component potassium_current (nanoS).
C CONSTS(115) is E_K in component potassium_current (millivolt).
C STATES(20) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(18) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(40) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(116) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(117) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(118) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(119) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(59) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(21) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(120) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(121) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(19) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(41) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(122) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(123) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(124) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(125) is g_L in component leakage_current (nanoS).
C CONSTS(126) is E_L in component leakage_current (millivolt).
C CONSTS(127) is E_syn_e in component tonic_current (millivolt).
C CONSTS(128) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(20) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(129) is kr in component synaptic_input (dimensionless).
C CONSTS(130) is tau_s in component synaptic_input (millisecond).
C CONSTS(131) is sigma_s in component synaptic_input (millivolt).
C CONSTS(132) is theta_s in component synaptic_input (millivolt).
C VOI is time in component environment (millisecond).
C STATES(22) is V in component membrane (millivolt).
C CONSTS(133) is C in component membrane (picoF).
C CONSTS(134) is i_app in component membrane (picoA).
C ALGBRC(65) is i_NaP in component persistent_sodium_current (picoA).
C ALGBRC(50) is i_Na in component fast_sodium_current (picoA).
C ALGBRC(55) is i_K in component potassium_current (picoA).
C ALGBRC(70) is i_L in component leakage_current (picoA).
C ALGBRC(75) is i_tonic_e in component tonic_current (picoA).
C ALGBRC(80) is i_syn_e in component synaptic_input (picoA).
C CONSTS(135) is E_Na in component fast_sodium_current (millivolt).
C CONSTS(136) is g_Na in component fast_sodium_current (nanoS).
C ALGBRC(42) is m_infinity in component fast_sodium_current_m_gate (dimensionless).
C STATES(23) is n in component fast_sodium_current_n_gate (dimensionless).
C CONSTS(137) is theta_m in component fast_sodium_current_m_gate (millivolt).
C CONSTS(138) is sigma_m in component fast_sodium_current_m_gate (millivolt).
C ALGBRC(22) is n_infinity in component fast_sodium_current_n_gate (dimensionless).
C ALGBRC(43) is tau_n in component fast_sodium_current_n_gate (millisecond).
C CONSTS(139) is tau_n_max in component fast_sodium_current_n_gate (millisecond).
C CONSTS(140) is theta_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(141) is sigma_n in component fast_sodium_current_n_gate (millivolt).
C CONSTS(142) is g_K in component potassium_current (nanoS).
C CONSTS(143) is E_K in component potassium_current (millivolt).
C STATES(24) is n in component potassium_current_n_gate (dimensionless).
C ALGBRC(23) is n_infinity in component potassium_current_n_gate (dimensionless).
C ALGBRC(44) is tau_n in component potassium_current_n_gate (millisecond).
C CONSTS(144) is tau_n_max in component potassium_current_n_gate (millisecond).
C CONSTS(145) is theta_n in component potassium_current_n_gate (millivolt).
C CONSTS(146) is sigma_n in component potassium_current_n_gate (millivolt).
C CONSTS(147) is g_NaP in component persistent_sodium_current (nanoS).
C ALGBRC(60) is m_infinity in component persistent_sodium_current_m_gate (dimensionless).
C STATES(25) is h in component persistent_sodium_current_h_gate (dimensionless).
C CONSTS(148) is theta_m in component persistent_sodium_current_m_gate (millivolt).
C CONSTS(149) is sigma_m in component persistent_sodium_current_m_gate (millivolt).
C ALGBRC(24) is h_infinity in component persistent_sodium_current_h_gate (dimensionless).
C ALGBRC(45) is tau_h in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(150) is tau_h_max in component persistent_sodium_current_h_gate (millisecond).
C CONSTS(151) is theta_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(152) is sigma_h in component persistent_sodium_current_h_gate (millivolt).
C CONSTS(153) is g_L in component leakage_current (nanoS).
C CONSTS(154) is E_L in component leakage_current (millivolt).
C CONSTS(155) is E_syn_e in component tonic_current (millivolt).
C CONSTS(156) is g_tonic_e in component tonic_current (nanoS).
C ALGBRC(25) is s_infinity in component synaptic_input (dimensionless).
C CONSTS(157) is kr in component synaptic_input (dimensionless).
C CONSTS(158) is tau_s in component synaptic_input (millisecond).
C CONSTS(159) is sigma_s in component synaptic_input (millivolt).
C CONSTS(160) is theta_s in component synaptic_input (millivolt).
C RATES(6) is d/dt V in component membrane (millivolt).
C RATES(7) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(8) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(9) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(1) is d/dt s in component synaptic_input (dimensionless).
C RATES(10) is d/dt V in component membrane (millivolt).
C RATES(11) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(12) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(13) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(2) is d/dt s in component synaptic_input (dimensionless).
C RATES(14) is d/dt V in component membrane (millivolt).
C RATES(15) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(16) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(17) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(3) is d/dt s in component synaptic_input (dimensionless).
C RATES(18) is d/dt V in component membrane (millivolt).
C RATES(19) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(20) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(21) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(4) is d/dt s in component synaptic_input (dimensionless).
C RATES(22) is d/dt V in component membrane (millivolt).
C RATES(23) is d/dt n in component fast_sodium_current_n_gate (dimensionless).
C RATES(24) is d/dt n in component potassium_current_n_gate (dimensionless).
C RATES(25) is d/dt h in component persistent_sodium_current_h_gate (dimensionless).
C RATES(5) is d/dt s in component synaptic_input (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      CONSTS(1) = 0.10
      CONSTS(2) = 0.10
      CONSTS(3) = 0.10
      CONSTS(4) = 0.10
      CONSTS(5) = 0.10
      CONSTS(6) = 0.10
      CONSTS(7) = 0.10
      CONSTS(8) = 0.10
      CONSTS(9) = 0.10
      CONSTS(10) = 0.10
      CONSTS(11) = 0.10
      CONSTS(12) = 0.10
      CONSTS(13) = 0.10
      CONSTS(14) = 0.10
      CONSTS(15) = 0.10
      CONSTS(16) = 0.10
      CONSTS(17) = 0.10
      CONSTS(18) = 0.10
      CONSTS(19) = 0.10
      CONSTS(20) = 0.10
      STATES(1) = 1
      STATES(2) = 1
      STATES(3) = 1
      STATES(4) = 1
      STATES(5) = 1
      STATES(6) = -50.0
      CONSTS(21) = 21.0
      CONSTS(22) = 0.0
      CONSTS(23) = 50.0
      CONSTS(24) = 28.0
      STATES(7) = 0.01
      CONSTS(25) = -34.0
      CONSTS(26) = -5.0
      CONSTS(27) = 10.0
      CONSTS(28) = -29.0
      CONSTS(29) = -4.0
      CONSTS(30) = 11.2
      CONSTS(31) = -85.0
      STATES(8) = 0.01
      CONSTS(32) = 10.0
      CONSTS(33) = -29.0
      CONSTS(34) = -4.0
      CONSTS(35) = 2.8
      STATES(9) = 0.46
      CONSTS(36) = -40.0
      CONSTS(37) = -6.0
      CONSTS(38) = 10000.0
      CONSTS(39) = -48.0
      CONSTS(40) = 6.0
      CONSTS(41) = 2.8
      CONSTS(42) = -57.5
      CONSTS(43) = 0.0
      CONSTS(44) = 0.0
      CONSTS(45) = 1.0
      CONSTS(46) = 5.0
      CONSTS(47) = -5.0
      CONSTS(48) = -10.0
      STATES(10) = -50.0
      CONSTS(49) = 21.0
      CONSTS(50) = 0.0
      CONSTS(51) = 50.0
      CONSTS(52) = 28.0
      STATES(11) = 0.01
      CONSTS(53) = -34.0
      CONSTS(54) = -5.0
      CONSTS(55) = 10.0
      CONSTS(56) = -29.0
      CONSTS(57) = -4.0
      CONSTS(58) = 11.2
      CONSTS(59) = -85.0
      STATES(12) = 0.01
      CONSTS(60) = 10.0
      CONSTS(61) = -29.0
      CONSTS(62) = -4.0
      CONSTS(63) = 2.8
      STATES(13) = 0.46
      CONSTS(64) = -40.0
      CONSTS(65) = -6.0
      CONSTS(66) = 10000.0
      CONSTS(67) = -48.0
      CONSTS(68) = 6.0
      CONSTS(69) = 2.8
      CONSTS(70) = -57.5
      CONSTS(71) = 0.0
      CONSTS(72) = 0.0
      CONSTS(73) = 1.0
      CONSTS(74) = 5.0
      CONSTS(75) = -5.0
      CONSTS(76) = -10.0
      STATES(14) = -50.0
      CONSTS(77) = 21.0
      CONSTS(78) = 0.0
      CONSTS(79) = 50.0
      CONSTS(80) = 28.0
      STATES(15) = 0.01
      CONSTS(81) = -34.0
      CONSTS(82) = -5.0
      CONSTS(83) = 10.0
      CONSTS(84) = -29.0
      CONSTS(85) = -4.0
      CONSTS(86) = 11.2
      CONSTS(87) = -85.0
      STATES(16) = 0.01
      CONSTS(88) = 10.0
      CONSTS(89) = -29.0
      CONSTS(90) = -4.0
      CONSTS(91) = 2.8
      STATES(17) = 0.46
      CONSTS(92) = -40.0
      CONSTS(93) = -6.0
      CONSTS(94) = 10000.0
      CONSTS(95) = -48.0
      CONSTS(96) = 6.0
      CONSTS(97) = 2.8
      CONSTS(98) = -57.5
      CONSTS(99) = 0.0
      CONSTS(100) = 0.0
      CONSTS(101) = 1.0
      CONSTS(102) = 5.0
      CONSTS(103) = -5.0
      CONSTS(104) = -10.0
      STATES(18) = -50.0
      CONSTS(105) = 21.0
      CONSTS(106) = 0.0
      CONSTS(107) = 50.0
      CONSTS(108) = 28.0
      STATES(19) = 0.01
      CONSTS(109) = -34.0
      CONSTS(110) = -5.0
      CONSTS(111) = 10.0
      CONSTS(112) = -29.0
      CONSTS(113) = -4.0
      CONSTS(114) = 11.2
      CONSTS(115) = -85.0
      STATES(20) = 0.01
      CONSTS(116) = 10.0
      CONSTS(117) = -29.0
      CONSTS(118) = -4.0
      CONSTS(119) = 2.8
      STATES(21) = 0.46
      CONSTS(120) = -40.0
      CONSTS(121) = -6.0
      CONSTS(122) = 10000.0
      CONSTS(123) = -48.0
      CONSTS(124) = 6.0
      CONSTS(125) = 2.8
      CONSTS(126) = -57.5
      CONSTS(127) = 0.0
      CONSTS(128) = 0.0
      CONSTS(129) = 1.0
      CONSTS(130) = 5.0
      CONSTS(131) = -5.0
      CONSTS(132) = -10.0
      STATES(22) = -50.0
      CONSTS(133) = 21.0
      CONSTS(134) = 0.0
      CONSTS(135) = 50.0
      CONSTS(136) = 28.0
      STATES(23) = 0.01
      CONSTS(137) = -34.0
      CONSTS(138) = -5.0
      CONSTS(139) = 10.0
      CONSTS(140) = -29.0
      CONSTS(141) = -4.0
      CONSTS(142) = 11.2
      CONSTS(143) = -85.0
      STATES(24) = 0.01
      CONSTS(144) = 10.0
      CONSTS(145) = -29.0
      CONSTS(146) = -4.0
      CONSTS(147) = 2.8
      STATES(25) = 0.46
      CONSTS(148) = -40.0
      CONSTS(149) = -6.0
      CONSTS(150) = 10000.0
      CONSTS(151) = -48.0
      CONSTS(152) = 6.0
      CONSTS(153) = 2.8
      CONSTS(154) = -57.5
      CONSTS(155) = 0.0
      CONSTS(156) = 0.0
      CONSTS(157) = 1.0
      CONSTS(158) = 5.0
      CONSTS(159) = -5.0
      CONSTS(160) = -10.0
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(48))/CONSTS(47)))
      RATES(1) = ( (1.00000 - STATES(1))*ALGBRC(5) - - ( CONSTS(45)*STATES(1)))/CONSTS(46)
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(76))/CONSTS(75)))
      RATES(2) = ( (1.00000 - STATES(2))*ALGBRC(10) - - ( CONSTS(73)*STATES(2)))/CONSTS(74)
      ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(104))/CONSTS(103)))
      RATES(3) = ( (1.00000 - STATES(3))*ALGBRC(15) - - ( CONSTS(101)*STATES(3)))/CONSTS(102)
      ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(132))/CONSTS(131)))
      RATES(4) = ( (1.00000 - STATES(4))*ALGBRC(20) - - ( CONSTS(129)*STATES(4)))/CONSTS(130)
      ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(160))/CONSTS(159)))
      RATES(5) = ( (1.00000 - STATES(5))*ALGBRC(25) - - ( CONSTS(157)*STATES(5)))/CONSTS(158)
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(28))/CONSTS(29)))
      ALGBRC(27) = CONSTS(27)/cosh((STATES(6) - CONSTS(28))/( 2.00000*CONSTS(29)))
      RATES(7) = (ALGBRC(2) - STATES(7))/ALGBRC(27)
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(33))/CONSTS(34)))
      ALGBRC(28) = CONSTS(32)/cosh((STATES(6) - CONSTS(33))/( 2.00000*CONSTS(34)))
      RATES(8) = (ALGBRC(3) - STATES(8))/ALGBRC(28)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(39))/CONSTS(40)))
      ALGBRC(29) = CONSTS(38)/cosh((STATES(6) - CONSTS(39))/( 2.00000*CONSTS(40)))
      RATES(9) = (ALGBRC(4) - STATES(9))/ALGBRC(29)
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(56))/CONSTS(57)))
      ALGBRC(31) = CONSTS(55)/cosh((STATES(10) - CONSTS(56))/( 2.00000*CONSTS(57)))
      RATES(11) = (ALGBRC(7) - STATES(11))/ALGBRC(31)
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(61))/CONSTS(62)))
      ALGBRC(32) = CONSTS(60)/cosh((STATES(10) - CONSTS(61))/( 2.00000*CONSTS(62)))
      RATES(12) = (ALGBRC(8) - STATES(12))/ALGBRC(32)
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(67))/CONSTS(68)))
      ALGBRC(33) = CONSTS(66)/cosh((STATES(10) - CONSTS(67))/( 2.00000*CONSTS(68)))
      RATES(13) = (ALGBRC(9) - STATES(13))/ALGBRC(33)
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(84))/CONSTS(85)))
      ALGBRC(35) = CONSTS(83)/cosh((STATES(14) - CONSTS(84))/( 2.00000*CONSTS(85)))
      RATES(15) = (ALGBRC(12) - STATES(15))/ALGBRC(35)
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(89))/CONSTS(90)))
      ALGBRC(36) = CONSTS(88)/cosh((STATES(14) - CONSTS(89))/( 2.00000*CONSTS(90)))
      RATES(16) = (ALGBRC(13) - STATES(16))/ALGBRC(36)
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(95))/CONSTS(96)))
      ALGBRC(37) = CONSTS(94)/cosh((STATES(14) - CONSTS(95))/( 2.00000*CONSTS(96)))
      RATES(17) = (ALGBRC(14) - STATES(17))/ALGBRC(37)
      ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(112))/CONSTS(113)))
      ALGBRC(39) = CONSTS(111)/cosh((STATES(18) - CONSTS(112))/( 2.00000*CONSTS(113)))
      RATES(19) = (ALGBRC(17) - STATES(19))/ALGBRC(39)
      ALGBRC(18) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(117))/CONSTS(118)))
      ALGBRC(40) = CONSTS(116)/cosh((STATES(18) - CONSTS(117))/( 2.00000*CONSTS(118)))
      RATES(20) = (ALGBRC(18) - STATES(20))/ALGBRC(40)
      ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(123))/CONSTS(124)))
      ALGBRC(41) = CONSTS(122)/cosh((STATES(18) - CONSTS(123))/( 2.00000*CONSTS(124)))
      RATES(21) = (ALGBRC(19) - STATES(21))/ALGBRC(41)
      ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(140))/CONSTS(141)))
      ALGBRC(43) = CONSTS(139)/cosh((STATES(22) - CONSTS(140))/( 2.00000*CONSTS(141)))
      RATES(23) = (ALGBRC(22) - STATES(23))/ALGBRC(43)
      ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(145))/CONSTS(146)))
      ALGBRC(44) = CONSTS(144)/cosh((STATES(22) - CONSTS(145))/( 2.00000*CONSTS(146)))
      RATES(24) = (ALGBRC(23) - STATES(24))/ALGBRC(44)
      ALGBRC(24) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(151))/CONSTS(152)))
      ALGBRC(45) = CONSTS(150)/cosh((STATES(22) - CONSTS(151))/( 2.00000*CONSTS(152)))
      RATES(25) = (ALGBRC(24) - STATES(25))/ALGBRC(45)
      ALGBRC(56) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(36))/CONSTS(37)))
      ALGBRC(61) =  CONSTS(35)*ALGBRC(56)*STATES(9)*(STATES(6) - CONSTS(23))
      ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(25))/CONSTS(26)))
      ALGBRC(46) =  CONSTS(24)*ALGBRC(26) ** 3.00000*(1.00000 - STATES(7))*(STATES(6) - CONSTS(23))
      ALGBRC(51) =  CONSTS(30)*STATES(8) ** 4.00000*(STATES(6) - CONSTS(31))
      ALGBRC(66) =  CONSTS(41)*(STATES(6) - CONSTS(42))
      ALGBRC(71) =  CONSTS(44)*(STATES(6) - CONSTS(43))
      ALGBRC(1) =  CONSTS(1)*STATES(1)+ CONSTS(2)*STATES(1)+ CONSTS(3)*STATES(1)+ CONSTS(4)*STATES(1)
      ALGBRC(76) =  ALGBRC(1)*(STATES(6) - CONSTS(43))
      RATES(6) = (- (ALGBRC(61)+ALGBRC(46)+ALGBRC(51)+ALGBRC(66)+ALGBRC(71)+ALGBRC(76))+CONSTS(22))/CONSTS(21)
      ALGBRC(57) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(64))/CONSTS(65)))
      ALGBRC(62) =  CONSTS(63)*ALGBRC(57)*STATES(13)*(STATES(10) - CONSTS(51))
      ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(53))/CONSTS(54)))
      ALGBRC(47) =  CONSTS(52)*ALGBRC(30) ** 3.00000*(1.00000 - STATES(11))*(STATES(10) - CONSTS(51))
      ALGBRC(52) =  CONSTS(58)*STATES(12) ** 4.00000*(STATES(10) - CONSTS(59))
      ALGBRC(67) =  CONSTS(69)*(STATES(10) - CONSTS(70))
      ALGBRC(72) =  CONSTS(72)*(STATES(10) - CONSTS(71))
      ALGBRC(6) =  CONSTS(5)*STATES(2)+ CONSTS(6)*STATES(2)+ CONSTS(7)*STATES(2)+ CONSTS(8)*STATES(2)
      ALGBRC(77) =  ALGBRC(6)*(STATES(10) - CONSTS(71))
      RATES(10) = (- (ALGBRC(62)+ALGBRC(47)+ALGBRC(52)+ALGBRC(67)+ALGBRC(72)+ALGBRC(77))+CONSTS(50))/CONSTS(49)
      ALGBRC(58) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(92))/CONSTS(93)))
      ALGBRC(63) =  CONSTS(91)*ALGBRC(58)*STATES(17)*(STATES(14) - CONSTS(79))
      ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(81))/CONSTS(82)))
      ALGBRC(48) =  CONSTS(80)*ALGBRC(34) ** 3.00000*(1.00000 - STATES(15))*(STATES(14) - CONSTS(79))
      ALGBRC(53) =  CONSTS(86)*STATES(16) ** 4.00000*(STATES(14) - CONSTS(87))
      ALGBRC(68) =  CONSTS(97)*(STATES(14) - CONSTS(98))
      ALGBRC(73) =  CONSTS(100)*(STATES(14) - CONSTS(99))
      ALGBRC(11) =  CONSTS(9)*STATES(3)+ CONSTS(10)*STATES(3)+ CONSTS(11)*STATES(3)+ CONSTS(12)*STATES(3)
      ALGBRC(78) =  ALGBRC(11)*(STATES(14) - CONSTS(99))
      RATES(14) = (- (ALGBRC(63)+ALGBRC(48)+ALGBRC(53)+ALGBRC(68)+ALGBRC(73)+ALGBRC(78))+CONSTS(78))/CONSTS(77)
      ALGBRC(59) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(120))/CONSTS(121)))
      ALGBRC(64) =  CONSTS(119)*ALGBRC(59)*STATES(21)*(STATES(18) - CONSTS(107))
      ALGBRC(38) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(109))/CONSTS(110)))
      ALGBRC(49) =  CONSTS(108)*ALGBRC(38) ** 3.00000*(1.00000 - STATES(19))*(STATES(18) - CONSTS(107))
      ALGBRC(54) =  CONSTS(114)*STATES(20) ** 4.00000*(STATES(18) - CONSTS(115))
      ALGBRC(69) =  CONSTS(125)*(STATES(18) - CONSTS(126))
      ALGBRC(74) =  CONSTS(128)*(STATES(18) - CONSTS(127))
      ALGBRC(16) =  CONSTS(13)*STATES(4)+ CONSTS(14)*STATES(4)+ CONSTS(15)*STATES(4)+ CONSTS(16)*STATES(4)
      ALGBRC(79) =  ALGBRC(16)*(STATES(18) - CONSTS(127))
      RATES(18) = (- (ALGBRC(64)+ALGBRC(49)+ALGBRC(54)+ALGBRC(69)+ALGBRC(74)+ALGBRC(79))+CONSTS(106))/CONSTS(105)
      ALGBRC(60) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(148))/CONSTS(149)))
      ALGBRC(65) =  CONSTS(147)*ALGBRC(60)*STATES(25)*(STATES(22) - CONSTS(135))
      ALGBRC(42) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(137))/CONSTS(138)))
      ALGBRC(50) =  CONSTS(136)*ALGBRC(42) ** 3.00000*(1.00000 - STATES(23))*(STATES(22) - CONSTS(135))
      ALGBRC(55) =  CONSTS(142)*STATES(24) ** 4.00000*(STATES(22) - CONSTS(143))
      ALGBRC(70) =  CONSTS(153)*(STATES(22) - CONSTS(154))
      ALGBRC(75) =  CONSTS(156)*(STATES(22) - CONSTS(155))
      ALGBRC(21) =  CONSTS(17)*STATES(5)+ CONSTS(18)*STATES(5)+ CONSTS(19)*STATES(5)+ CONSTS(20)*STATES(5)
      ALGBRC(80) =  ALGBRC(21)*(STATES(22) - CONSTS(155))
      RATES(22) = (- (ALGBRC(65)+ALGBRC(50)+ALGBRC(55)+ALGBRC(70)+ALGBRC(75)+ALGBRC(80))+CONSTS(134))/CONSTS(133)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(48))/CONSTS(47)))
      ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(76))/CONSTS(75)))
      ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(104))/CONSTS(103)))
      ALGBRC(20) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(132))/CONSTS(131)))
      ALGBRC(25) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(160))/CONSTS(159)))
      ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(28))/CONSTS(29)))
      ALGBRC(27) = CONSTS(27)/cosh((STATES(6) - CONSTS(28))/( 2.00000*CONSTS(29)))
      ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(33))/CONSTS(34)))
      ALGBRC(28) = CONSTS(32)/cosh((STATES(6) - CONSTS(33))/( 2.00000*CONSTS(34)))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(39))/CONSTS(40)))
      ALGBRC(29) = CONSTS(38)/cosh((STATES(6) - CONSTS(39))/( 2.00000*CONSTS(40)))
      ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(56))/CONSTS(57)))
      ALGBRC(31) = CONSTS(55)/cosh((STATES(10) - CONSTS(56))/( 2.00000*CONSTS(57)))
      ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(61))/CONSTS(62)))
      ALGBRC(32) = CONSTS(60)/cosh((STATES(10) - CONSTS(61))/( 2.00000*CONSTS(62)))
      ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(67))/CONSTS(68)))
      ALGBRC(33) = CONSTS(66)/cosh((STATES(10) - CONSTS(67))/( 2.00000*CONSTS(68)))
      ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(84))/CONSTS(85)))
      ALGBRC(35) = CONSTS(83)/cosh((STATES(14) - CONSTS(84))/( 2.00000*CONSTS(85)))
      ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(89))/CONSTS(90)))
      ALGBRC(36) = CONSTS(88)/cosh((STATES(14) - CONSTS(89))/( 2.00000*CONSTS(90)))
      ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(95))/CONSTS(96)))
      ALGBRC(37) = CONSTS(94)/cosh((STATES(14) - CONSTS(95))/( 2.00000*CONSTS(96)))
      ALGBRC(17) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(112))/CONSTS(113)))
      ALGBRC(39) = CONSTS(111)/cosh((STATES(18) - CONSTS(112))/( 2.00000*CONSTS(113)))
      ALGBRC(18) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(117))/CONSTS(118)))
      ALGBRC(40) = CONSTS(116)/cosh((STATES(18) - CONSTS(117))/( 2.00000*CONSTS(118)))
      ALGBRC(19) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(123))/CONSTS(124)))
      ALGBRC(41) = CONSTS(122)/cosh((STATES(18) - CONSTS(123))/( 2.00000*CONSTS(124)))
      ALGBRC(22) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(140))/CONSTS(141)))
      ALGBRC(43) = CONSTS(139)/cosh((STATES(22) - CONSTS(140))/( 2.00000*CONSTS(141)))
      ALGBRC(23) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(145))/CONSTS(146)))
      ALGBRC(44) = CONSTS(144)/cosh((STATES(22) - CONSTS(145))/( 2.00000*CONSTS(146)))
      ALGBRC(24) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(151))/CONSTS(152)))
      ALGBRC(45) = CONSTS(150)/cosh((STATES(22) - CONSTS(151))/( 2.00000*CONSTS(152)))
      ALGBRC(56) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(36))/CONSTS(37)))
      ALGBRC(61) =  CONSTS(35)*ALGBRC(56)*STATES(9)*(STATES(6) - CONSTS(23))
      ALGBRC(26) = 1.00000/(1.00000+EXP((STATES(6) - CONSTS(25))/CONSTS(26)))
      ALGBRC(46) =  CONSTS(24)*ALGBRC(26) ** 3.00000*(1.00000 - STATES(7))*(STATES(6) - CONSTS(23))
      ALGBRC(51) =  CONSTS(30)*STATES(8) ** 4.00000*(STATES(6) - CONSTS(31))
      ALGBRC(66) =  CONSTS(41)*(STATES(6) - CONSTS(42))
      ALGBRC(71) =  CONSTS(44)*(STATES(6) - CONSTS(43))
      ALGBRC(1) =  CONSTS(1)*STATES(1)+ CONSTS(2)*STATES(1)+ CONSTS(3)*STATES(1)+ CONSTS(4)*STATES(1)
      ALGBRC(76) =  ALGBRC(1)*(STATES(6) - CONSTS(43))
      ALGBRC(57) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(64))/CONSTS(65)))
      ALGBRC(62) =  CONSTS(63)*ALGBRC(57)*STATES(13)*(STATES(10) - CONSTS(51))
      ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(10) - CONSTS(53))/CONSTS(54)))
      ALGBRC(47) =  CONSTS(52)*ALGBRC(30) ** 3.00000*(1.00000 - STATES(11))*(STATES(10) - CONSTS(51))
      ALGBRC(52) =  CONSTS(58)*STATES(12) ** 4.00000*(STATES(10) - CONSTS(59))
      ALGBRC(67) =  CONSTS(69)*(STATES(10) - CONSTS(70))
      ALGBRC(72) =  CONSTS(72)*(STATES(10) - CONSTS(71))
      ALGBRC(6) =  CONSTS(5)*STATES(2)+ CONSTS(6)*STATES(2)+ CONSTS(7)*STATES(2)+ CONSTS(8)*STATES(2)
      ALGBRC(77) =  ALGBRC(6)*(STATES(10) - CONSTS(71))
      ALGBRC(58) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(92))/CONSTS(93)))
      ALGBRC(63) =  CONSTS(91)*ALGBRC(58)*STATES(17)*(STATES(14) - CONSTS(79))
      ALGBRC(34) = 1.00000/(1.00000+EXP((STATES(14) - CONSTS(81))/CONSTS(82)))
      ALGBRC(48) =  CONSTS(80)*ALGBRC(34) ** 3.00000*(1.00000 - STATES(15))*(STATES(14) - CONSTS(79))
      ALGBRC(53) =  CONSTS(86)*STATES(16) ** 4.00000*(STATES(14) - CONSTS(87))
      ALGBRC(68) =  CONSTS(97)*(STATES(14) - CONSTS(98))
      ALGBRC(73) =  CONSTS(100)*(STATES(14) - CONSTS(99))
      ALGBRC(11) =  CONSTS(9)*STATES(3)+ CONSTS(10)*STATES(3)+ CONSTS(11)*STATES(3)+ CONSTS(12)*STATES(3)
      ALGBRC(78) =  ALGBRC(11)*(STATES(14) - CONSTS(99))
      ALGBRC(59) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(120))/CONSTS(121)))
      ALGBRC(64) =  CONSTS(119)*ALGBRC(59)*STATES(21)*(STATES(18) - CONSTS(107))
      ALGBRC(38) = 1.00000/(1.00000+EXP((STATES(18) - CONSTS(109))/CONSTS(110)))
      ALGBRC(49) =  CONSTS(108)*ALGBRC(38) ** 3.00000*(1.00000 - STATES(19))*(STATES(18) - CONSTS(107))
      ALGBRC(54) =  CONSTS(114)*STATES(20) ** 4.00000*(STATES(18) - CONSTS(115))
      ALGBRC(69) =  CONSTS(125)*(STATES(18) - CONSTS(126))
      ALGBRC(74) =  CONSTS(128)*(STATES(18) - CONSTS(127))
      ALGBRC(16) =  CONSTS(13)*STATES(4)+ CONSTS(14)*STATES(4)+ CONSTS(15)*STATES(4)+ CONSTS(16)*STATES(4)
      ALGBRC(79) =  ALGBRC(16)*(STATES(18) - CONSTS(127))
      ALGBRC(60) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(148))/CONSTS(149)))
      ALGBRC(65) =  CONSTS(147)*ALGBRC(60)*STATES(25)*(STATES(22) - CONSTS(135))
      ALGBRC(42) = 1.00000/(1.00000+EXP((STATES(22) - CONSTS(137))/CONSTS(138)))
      ALGBRC(50) =  CONSTS(136)*ALGBRC(42) ** 3.00000*(1.00000 - STATES(23))*(STATES(22) - CONSTS(135))
      ALGBRC(55) =  CONSTS(142)*STATES(24) ** 4.00000*(STATES(22) - CONSTS(143))
      ALGBRC(70) =  CONSTS(153)*(STATES(22) - CONSTS(154))
      ALGBRC(75) =  CONSTS(156)*(STATES(22) - CONSTS(155))
      ALGBRC(21) =  CONSTS(17)*STATES(5)+ CONSTS(18)*STATES(5)+ CONSTS(19)*STATES(5)+ CONSTS(20)*STATES(5)
      ALGBRC(80) =  ALGBRC(21)*(STATES(22) - CONSTS(155))
      RETURN
      END