Generated Code
The following is f77 code generated by the CellML API from this CellML file. (Back to language selection)
The raw code is available.
C
C There are a total of 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