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 128 entries in the algebraic variable array.
C There are a total of 48 entries in each of the rate and state variable arrays.
C There are a total of 94 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C CONSTS(1) is F in component model_parameters (C_per_mole).
C CONSTS(2) is T in component model_parameters (kelvin).
C CONSTS(3) is R in component model_parameters (J_per_kmole_K).
C CONSTS(4) is K_o in component model_parameters (mM).
C CONSTS(5) is Ca_o in component model_parameters (mM).
C CONSTS(6) is Na_o in component model_parameters (mM).
C CONSTS(7) is Cl_o in component model_parameters (mM).
C CONSTS(89) is Vmyo in component cell_geometry (uL).
C CONSTS(90) is Vnsr in component cell_geometry (uL).
C CONSTS(85) is Vsr in component cell_geometry (uL).
C CONSTS(91) is Vjsr in component cell_geometry (uL).
C CONSTS(92) is Vss_sr in component cell_geometry (uL).
C CONSTS(93) is Vss_CaL in component cell_geometry (uL).
C CONSTS(87) is Acap in component cell_geometry (cm2).
C CONSTS(88) is AF in component cell_geometry (uF_mole_per_C).
C CONSTS(8) is l in component cell_geometry (cm).
C CONSTS(9) is a in component cell_geometry (cm).
C CONSTS(82) is vcell in component cell_geometry (uL).
C CONSTS(86) is ageo in component cell_geometry (cm2).
C STATES(1) is Vm in component membrane (mV).
C ALGBRC(13) is i_Stim in component membrane (uA_per_uF).
C CONSTS(10) is stim_offset in component membrane (ms).
C CONSTS(11) is stim_period in component membrane (ms).
C CONSTS(12) is stim_duration in component membrane (ms).
C CONSTS(13) is stim_amplitude in component membrane (uA_per_uF).
C ALGBRC(1) is past in component membrane (ms).
C ALGBRC(111) is caiont in component membrane (uA_per_uF).
C ALGBRC(114) is naiont in component membrane (uA_per_uF).
C ALGBRC(81) is kiont in component membrane (uA_per_uF).
C ALGBRC(72) is clont in component membrane (uA_per_uF).
C ALGBRC(52) is INa in component INa (uA_per_uF).
C ALGBRC(99) is ICaL in component ICaL (uA_per_uF).
C ALGBRC(62) is IK1 in component IK1 (uA_per_uF).
C ALGBRC(55) is IKp in component IKp (uA_per_uF).
C ALGBRC(80) is IKs in component IKs (uA_per_uF).
C ALGBRC(66) is IKr in component IKr (uA_per_uF).
C ALGBRC(87) is IpCa in component IpCa (uA_per_uF).
C ALGBRC(88) is ICab in component ICab (uA_per_uF).
C ALGBRC(86) is INaCa in component INaCa (uA_per_uF).
C ALGBRC(110) is INaCa_ss_sr in component INaCa (uA_per_uF).
C ALGBRC(45) is INaK in component INaK (uA_per_uF).
C ALGBRC(68) is Ito1 in component Ito1 (uA_per_uF).
C ALGBRC(71) is Ito2 in component Ito2 (uA_per_uF).
C ALGBRC(58) is IClb in component IClb (uA_per_uF).
C ALGBRC(53) is INaL in component INaL (uA_per_uF).
C ALGBRC(50) is INab in component INab (uA_per_uF).
C ALGBRC(79) is gKs in component IKs (mS_per_uF).
C STATES(2) is C1 in component IKs (dimensionless).
C STATES(3) is C2 in component IKs (dimensionless).
C STATES(4) is C3 in component IKs (dimensionless).
C STATES(5) is C4 in component IKs (dimensionless).
C STATES(6) is C5 in component IKs (dimensionless).
C STATES(7) is C6 in component IKs (dimensionless).
C STATES(8) is C7 in component IKs (dimensionless).
C STATES(9) is C8 in component IKs (dimensionless).
C STATES(10) is C9 in component IKs (dimensionless).
C STATES(11) is C10 in component IKs (dimensionless).
C STATES(12) is C11 in component IKs (dimensionless).
C STATES(13) is C12 in component IKs (dimensionless).
C STATES(14) is C13 in component IKs (dimensionless).
C STATES(15) is C14 in component IKs (dimensionless).
C STATES(16) is C15 in component IKs (dimensionless).
C STATES(17) is O1 in component IKs (dimensionless).
C STATES(18) is O2 in component IKs (dimensionless).
C ALGBRC(24) is OKs in component IKs (dimensionless).
C ALGBRC(2) is alpha in component IKs (per_ms).
C ALGBRC(14) is beta in component IKs (per_ms).
C ALGBRC(29) is delta in component IKs (per_ms).
C ALGBRC(25) is gamma in component IKs (per_ms).
C ALGBRC(32) is eta in component IKs (per_ms).
C CONSTS(81) is theta in component IKs (per_ms).
C ALGBRC(34) is omega in component IKs (per_ms).
C ALGBRC(36) is psi in component IKs (per_ms).
C ALGBRC(56) is EKs in component reversal_potentials (mV).
C ALGBRC(78) is Ca_i in component Ca_buffering (mM).
C ALGBRC(97) is ICaL_max in component ICaL (uA_per_uF).
C CONSTS(14) is PCa in component ICaL (L_per_F_ms).
C STATES(19) is C in component ICaL (dimensionless).
C STATES(20) is O in component ICaL (dimensionless).
C STATES(21) is C_star in component ICaL (dimensionless).
C STATES(22) is O_star in component ICaL (dimensionless).
C STATES(23) is CI in component ICaL (dimensionless).
C STATES(24) is OI in component ICaL (dimensionless).
C STATES(25) is CI_star in component ICaL (dimensionless).
C STATES(26) is OI_star in component ICaL (dimensionless).
C ALGBRC(28) is ACT_tau in component ICaL (ms).
C ALGBRC(38) is IV_tau in component ICaL (ms).
C ALGBRC(100) is IV_tau_star in component ICaL (ms).
C ALGBRC(31) is ACT_infinity in component ICaL (dimensionless).
C ALGBRC(37) is IV_infinity in component ICaL (dimensionless).
C ALGBRC(41) is IV_infinity_star in component ICaL (dimensionless).
C ALGBRC(98) is IV_beta_star in component ICaL (dimensionless).
C ALGBRC(96) is IV_beta_infinity_star in component ICaL (dimensionless).
C ALGBRC(33) is alpha in component ICaL (per_ms).
C ALGBRC(35) is beta in component ICaL (per_ms).
C ALGBRC(106) is delta in component ICaL (per_ms).
C ALGBRC(39) is x in component ICaL (per_ms).
C ALGBRC(40) is y in component ICaL (per_ms).
C CONSTS(15) is theta in component ICaL (per_ms).
C ALGBRC(108) is delta_I in component ICaL (per_ms).
C ALGBRC(102) is x_star in component ICaL (per_ms).
C ALGBRC(104) is y_star in component ICaL (per_ms).
C CONSTS(16) is theta_I in component ICaL (per_ms).
C CONSTS(17) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(18) is gamma_Cao in component ICaL (dimensionless).
C ALGBRC(95) is Ca_ss_CaL in component Ca_ss_buffering (mM).
C CONSTS(19) is g_Na in component INa (mS_per_uF).
C ALGBRC(51) is ENa in component reversal_potentials (mV).
C STATES(27) is m in component INa_m_gate (dimensionless).
C STATES(28) is h in component INa_h_gate (dimensionless).
C STATES(29) is j in component INa_j_gate (dimensionless).
C ALGBRC(3) is am in component INa_m_gate (per_ms).
C ALGBRC(15) is bm in component INa_m_gate (per_ms).
C ALGBRC(4) is ah in component INa_h_gate (per_ms).
C ALGBRC(16) is bh in component INa_h_gate (per_ms).
C ALGBRC(5) is aj in component INa_j_gate (per_ms).
C ALGBRC(17) is bj in component INa_j_gate (per_ms).
C CONSTS(20) is kmnai in component INaK (mM).
C CONSTS(21) is kmko in component INaK (mM).
C CONSTS(22) is ibarnak in component INaK (uA_per_uF).
C CONSTS(23) is sigma in component INaK (dimensionless).
C CONSTS(24) is V_half in component INaK (mV).
C ALGBRC(42) is phi in component INaK (dimensionless).
C ALGBRC(43) is fv in component INaK (dimensionless).
C CONSTS(94) is PK in component INaK (dimensionless).
C ALGBRC(44) is PNa in component INaK (dimensionless).
C STATES(30) is Na_i in component Na (mM).
C CONSTS(25) is KmCa_act in component INaCa (mM).
C CONSTS(26) is Vmax in component INaCa (uA_per_uF).
C CONSTS(27) is ksat in component INaCa (dimensionless).
C CONSTS(28) is eta in component INaCa (dimensionless).
C CONSTS(29) is KmNai in component INaCa (mM).
C CONSTS(30) is KmNao in component INaCa (mM).
C CONSTS(31) is KmCai in component INaCa (mM).
C CONSTS(32) is KmCao in component INaCa (mM).
C ALGBRC(101) is Ca_ss_sr in component Ca_ss_buffering (mM).
C STATES(31) is Na_ss_sr in component Na (mM).
C ALGBRC(82) is num in component INaCa (mM4_uA_per_uF).
C ALGBRC(103) is num_ss in component INaCa (mM4_uA_per_uF).
C ALGBRC(83) is denom_1 in component INaCa (dimensionless).
C ALGBRC(46) is denom_2 in component INaCa (dimensionless).
C ALGBRC(84) is denom_3 in component INaCa (mM4).
C ALGBRC(85) is denom_4 in component INaCa (mM4).
C ALGBRC(105) is denom_ss_1 in component INaCa (dimensionless).
C ALGBRC(47) is denom_ss_2 in component INaCa (dimensionless).
C ALGBRC(107) is denom_ss_3 in component INaCa (mM4).
C ALGBRC(109) is denom_ss_4 in component INaCa (mM4).
C CONSTS(33) is gKp in component IKp (mS_per_uF).
C ALGBRC(48) is Kp in component IKp (dimensionless).
C ALGBRC(54) is EK in component reversal_potentials (mV).
C CONSTS(34) is gpCa in component IpCa (uA_per_uF).
C CONSTS(35) is KmpCa in component IpCa (mM).
C CONSTS(36) is PCab in component ICab (L_per_F_ms).
C CONSTS(37) is gamma_Ca_i in component ICab (dimensionless).
C CONSTS(38) is gamma_Ca_o in component ICab (dimensionless).
C CONSTS(39) is PNab in component INab (L_per_F_ms).
C ALGBRC(49) is phi in component INab (dimensionless).
C CONSTS(40) is gClb in component IClb (mS_per_uF).
C ALGBRC(57) is ECl in component reversal_potentials (mV).
C CONSTS(41) is gNaL in component INaL (mS_per_uF).
C STATES(32) is mL in component INaL_mL_gate (dimensionless).
C STATES(33) is hL in component INaL_hL_gate (dimensionless).
C ALGBRC(6) is amL in component INaL_mL_gate (per_ms).
C ALGBRC(18) is bmL in component INaL_mL_gate (per_ms).
C CONSTS(42) is tau_hL in component INaL_hL_gate (ms).
C ALGBRC(7) is hL_infinity in component INaL_hL_gate (dimensionless).
C CONSTS(43) is prnak in component reversal_potentials (dimensionless).
C STATES(34) is K_i in component K (mM).
C STATES(35) is Cl_i in component Cl (mM).
C CONSTS(83) is gK1 in component IK1 (mS_per_uF).
C ALGBRC(61) is K1 in component IK1_K1_gate (dimensionless).
C ALGBRC(59) is alpha_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(60) is beta_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(63) is CT_Na_Cl in component CT_Na_Cl (mM_per_ms).
C CONSTS(44) is CT_Na_Cl_max in component CT_Na_Cl (mM_per_ms).
C ALGBRC(64) is CT_K_Cl in component CT_K_Cl (mM_per_ms).
C CONSTS(45) is CT_K_Cl_max in component CT_K_Cl (mM_per_ms).
C CONSTS(84) is gKr in component IKr (mS_per_uF).
C ALGBRC(65) is r in component IKr_r_gate (dimensionless).
C STATES(36) is xr in component IKr_xr_gate (dimensionless).
C ALGBRC(8) is xr_infinity in component IKr_xr_gate (dimensionless).
C ALGBRC(19) is tau_xr in component IKr_xr_gate (ms).
C CONSTS(46) is gto1 in component Ito1 (mS_per_uF).
C ALGBRC(67) is rto1 in component Ito1 (dimensionless).
C STATES(37) is a in component Ito1_a_gate (dimensionless).
C STATES(38) is i1f in component Ito1_i1f_gate (dimensionless).
C STATES(39) is i1s in component Ito1_i1s_gate (dimensionless).
C ALGBRC(9) is alpha_a in component Ito1_a_gate (per_ms).
C ALGBRC(20) is beta_a in component Ito1_a_gate (per_ms).
C ALGBRC(26) is tau_a in component Ito1_a_gate (ms).
C ALGBRC(30) is a_infinity in component Ito1_a_gate (dimensionless).
C ALGBRC(21) is alpha_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(10) is beta_i1f in component Ito1_i1f_gate (per_ms).
C ALGBRC(22) is alpha_i1s in component Ito1_i1s_gate (per_ms).
C ALGBRC(11) is beta_i1s in component Ito1_i1s_gate (per_ms).
C CONSTS(47) is PCl in component Ito2 (L_per_F_ms).
C CONSTS(48) is zCl in component Ito2 (dimensionless).
C ALGBRC(69) is Ito2_max in component Ito2 (uA_per_uF).
C ALGBRC(70) is KCa_ito2 in component Ito2 (dimensionless).
C CONSTS(49) is kCa_ito2 in component Ito2 (mM_per_ms).
C STATES(40) is Irel in component Irel (mM_per_ms).
C STATES(41) is i2f in component Ito2_i2f_gate (dimensionless).
C ALGBRC(23) is alpha_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(12) is beta_i2f in component Ito2_i2f_gate (per_ms).
C ALGBRC(27) is i2f_infinity in component Ito2_i2f_gate (dimensionless).
C CONSTS(50) is tau_i2f in component Ito2_i2f_gate (ms).
C ALGBRC(128) is rel_infinity in component Irel (mM_per_ms).
C ALGBRC(127) is tau_rel in component Irel (ms).
C CONSTS(51) is Krel_infinity in component Irel (mM).
C ALGBRC(125) is alpha_rel in component Irel (mM_per_uA_per_uF_per_ms).
C ALGBRC(123) is beta_tau in component Irel (ms).
C CONSTS(52) is beta_0 in component Irel (ms).
C CONSTS(53) is kappa in component Irel (mM_per_uA_per_uF_per_ms_per_ms).
C CONSTS(54) is h_rel in component Irel (dimensionless).
C CONSTS(55) is h_beta in component Irel (dimensionless).
C ALGBRC(121) is delta_beta_CaMK in component Irel (dimensionless).
C CONSTS(56) is delta_beta_0 in component Irel (dimensionless).
C CONSTS(57) is K_beta in component Irel (mM).
C CONSTS(58) is Krel_tau in component Irel (mM).
C ALGBRC(118) is Ca_JSR in component Ca_JSR_buffering (mM).
C ALGBRC(120) is CaMK_active in component CaMK_active (dimensionless).
C ALGBRC(126) is Iup in component Iup (mM_per_ms).
C CONSTS(59) is kmup in component Iup (mM).
C ALGBRC(124) is delta_kmPLB in component Iup (mM).
C CONSTS(60) is delta_kmPLB_bar in component Iup (mM).
C CONSTS(61) is iupbar in component Iup (mM_per_ms).
C ALGBRC(122) is delta_iupCaMK in component Iup (dimensionless).
C CONSTS(62) is delta_iupCaMK_bar in component Iup (dimensionless).
C CONSTS(63) is kmCaMK in component CaMK_active (mM).
C ALGBRC(74) is Ileak in component Ileak (mM_per_ms).
C CONSTS(64) is nsrbar in component Ileak (mM).
C STATES(42) is Ca_NSR in component Ca (mM).
C ALGBRC(119) is Itr in component Itr (mM_per_ms).
C CONSTS(65) is tautr in component Itr (ms).
C STATES(43) is Ca_JSR_0 in component Ca (mM).
C STATES(44) is Ca_i_0 in component Ca (mM).
C STATES(45) is Ca_ss_sr_0 in component Ca (mM).
C STATES(46) is Ca_ss_CaL_0 in component Ca (mM).
C CONSTS(66) is tau_diff_ss in component Ca (ms).
C CONSTS(67) is tau_diff in component Ca (ms).
C ALGBRC(112) is Idiff in component Ca (mM_per_ms).
C ALGBRC(116) is Idiff_ss in component Ca (mM_per_ms).
C ALGBRC(75) is b_myo in component Ca_buffering (mM).
C ALGBRC(76) is c_myo in component Ca_buffering (mM2).
C ALGBRC(77) is d_myo in component Ca_buffering (mM3).
C CONSTS(68) is cmdn_bar in component Ca_buffering (mM).
C CONSTS(69) is trpn_bar in component Ca_buffering (mM).
C CONSTS(70) is km_cmdn in component Ca_buffering (mM).
C CONSTS(71) is km_trpn in component Ca_buffering (mM).
C ALGBRC(89) is b_myo1 in component Ca_ss_buffering (mM).
C ALGBRC(90) is c_myo1 in component Ca_ss_buffering (mM2).
C ALGBRC(91) is d_myo1 in component Ca_ss_buffering (mM3).
C ALGBRC(92) is b_myo2 in component Ca_ss_buffering (mM).
C ALGBRC(93) is c_myo2 in component Ca_ss_buffering (mM2).
C ALGBRC(94) is d_myo2 in component Ca_ss_buffering (mM3).
C CONSTS(72) is BSL_max in component Ca_ss_buffering (mM).
C CONSTS(73) is BSR_max in component Ca_ss_buffering (mM).
C CONSTS(74) is Km_BSR in component Ca_ss_buffering (mM).
C CONSTS(75) is Km_BSL in component Ca_ss_buffering (mM).
C ALGBRC(115) is b in component Ca_JSR_buffering (mM).
C ALGBRC(117) is c in component Ca_JSR_buffering (mM2).
C CONSTS(76) is csqn_bar in component Ca_JSR_buffering (mM).
C CONSTS(77) is km_csqn in component Ca_JSR_buffering (mM).
C ALGBRC(113) is Idiff_Na in component Na (mM_per_ms).
C STATES(47) is Cl_ss in component Cl (mM).
C ALGBRC(73) is Idiff_Cl in component Cl (mM_per_ms).
C STATES(48) is CaMK_trap in component CaMK_active (dimensionless).
C CONSTS(78) is alpha_CaMK in component CaMK_active (per_ms).
C CONSTS(79) is beta_CaMK in component CaMK_active (per_ms).
C CONSTS(80) is CaMK_0 in component CaMK_active (dimensionless).
C RATES(1) is d/dt Vm in component membrane (mV).
C RATES(2) is d/dt C1 in component IKs (dimensionless).
C RATES(3) is d/dt C2 in component IKs (dimensionless).
C RATES(4) is d/dt C3 in component IKs (dimensionless).
C RATES(5) is d/dt C4 in component IKs (dimensionless).
C RATES(6) is d/dt C5 in component IKs (dimensionless).
C RATES(7) is d/dt C6 in component IKs (dimensionless).
C RATES(8) is d/dt C7 in component IKs (dimensionless).
C RATES(9) is d/dt C8 in component IKs (dimensionless).
C RATES(10) is d/dt C9 in component IKs (dimensionless).
C RATES(11) is d/dt C10 in component IKs (dimensionless).
C RATES(12) is d/dt C11 in component IKs (dimensionless).
C RATES(13) is d/dt C12 in component IKs (dimensionless).
C RATES(14) is d/dt C13 in component IKs (dimensionless).
C RATES(15) is d/dt C14 in component IKs (dimensionless).
C RATES(16) is d/dt C15 in component IKs (dimensionless).
C RATES(17) is d/dt O1 in component IKs (dimensionless).
C RATES(18) is d/dt O2 in component IKs (dimensionless).
C RATES(19) is d/dt C in component ICaL (dimensionless).
C RATES(20) is d/dt O in component ICaL (dimensionless).
C RATES(21) is d/dt C_star in component ICaL (dimensionless).
C RATES(22) is d/dt O_star in component ICaL (dimensionless).
C RATES(23) is d/dt CI in component ICaL (dimensionless).
C RATES(24) is d/dt OI in component ICaL (dimensionless).
C RATES(25) is d/dt CI_star in component ICaL (dimensionless).
C RATES(26) is d/dt OI_star in component ICaL (dimensionless).
C RATES(27) is d/dt m in component INa_m_gate (dimensionless).
C RATES(28) is d/dt h in component INa_h_gate (dimensionless).
C RATES(29) is d/dt j in component INa_j_gate (dimensionless).
C RATES(32) is d/dt mL in component INaL_mL_gate (dimensionless).
C RATES(33) is d/dt hL in component INaL_hL_gate (dimensionless).
C RATES(36) is d/dt xr in component IKr_xr_gate (dimensionless).
C RATES(37) is d/dt a in component Ito1_a_gate (dimensionless).
C RATES(38) is d/dt i1f in component Ito1_i1f_gate (dimensionless).
C RATES(39) is d/dt i1s in component Ito1_i1s_gate (dimensionless).
C RATES(41) is d/dt i2f in component Ito2_i2f_gate (dimensionless).
C RATES(40) is d/dt Irel in component Irel (mM_per_ms).
C RATES(44) is d/dt Ca_i_0 in component Ca (mM).
C RATES(45) is d/dt Ca_ss_sr_0 in component Ca (mM).
C RATES(46) is d/dt Ca_ss_CaL_0 in component Ca (mM).
C RATES(42) is d/dt Ca_NSR in component Ca (mM).
C RATES(43) is d/dt Ca_JSR_0 in component Ca (mM).
C RATES(30) is d/dt Na_i in component Na (mM).
C RATES(31) is d/dt Na_ss_sr in component Na (mM).
C RATES(35) is d/dt Cl_i in component Cl (mM).
C RATES(47) is d/dt Cl_ss in component Cl (mM).
C RATES(34) is d/dt K_i in component K (mM).
C RATES(48) is d/dt CaMK_trap in component CaMK_active (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 96487
CONSTS(2) = 310
CONSTS(3) = 8314
CONSTS(4) = 5.4
CONSTS(5) = 1.8
CONSTS(6) = 140
CONSTS(7) = 100
CONSTS(8) = 0.01
CONSTS(9) = 0.0011
STATES(1) = -87.472528
CONSTS(10) = 0
CONSTS(11) = 1e3
CONSTS(12) = 0.5
CONSTS(13) = -80
STATES(2) = 0.476653873938053
STATES(3) = 0.2551771850026437
STATES(4) = 0.05122852186842027
STATES(5) = 0.004570874421666423
STATES(6) = 0.0001529389211563687
STATES(7) = 0.1328281770413747
STATES(8) = 0.05333228495514583
STATES(9) = 0.007137874613383642
STATES(10) = 0.0003184390970637536
STATES(11) = 0.01388061042363991
STATES(12) = 0.003715500164470033
STATES(13) = 0.0002486371465457659
STATES(14) = 0.0006446803014431032
STATES(15) = 8.628258026770862e-5
STATES(16) = 1.122823524761588e-5
STATES(17) = 1.051906959485266e-5
STATES(18) = 2.372219883411999e-6
CONSTS(14) = 1.5552e-4
STATES(19) = 0.9999613782488198
STATES(20) = 2.27419494597377e-5
STATES(21) = 5.936477396599124e-11
STATES(22) = 1.350122566369778e-15
STATES(23) = 1.587938121237188e-5
STATES(24) = 3.611420243201921e-10
STATES(25) = 0
STATES(26) = 0
CONSTS(15) = 1
CONSTS(16) = 0.000001
CONSTS(17) = 1
CONSTS(18) = 0.341
CONSTS(19) = 9.075
STATES(27) = 0.001022252813581797
STATES(28) = 0.9913109066867
STATES(29) = 0.9970771239754529
CONSTS(20) = 2.6
CONSTS(21) = 1.5
CONSTS(22) = 1.4
CONSTS(23) = 1
CONSTS(24) = -92
STATES(30) = 6.935575
CONSTS(25) = 1.25e-4
CONSTS(26) = 4.5
CONSTS(27) = 0.32
CONSTS(28) = 0.27
CONSTS(29) = 12.3
CONSTS(30) = 87.5
CONSTS(31) = 0.0036
CONSTS(32) = 1.3
STATES(31) = 6.935607
CONSTS(33) = 0.00276
CONSTS(34) = 0.0575
CONSTS(35) = 0.5e-3
CONSTS(36) = 1.995084e-7
CONSTS(37) = 1
CONSTS(38) = 0.341
CONSTS(39) = 0.32e-8
CONSTS(40) = 2.25e-4
CONSTS(41) = 0.0065
STATES(32) = 0.001022252813581797
STATES(33) = 0.3593255076983574
CONSTS(42) = 600
CONSTS(43) = 0.01833
STATES(34) = 145.5689
STATES(35) = 20.26842
CONSTS(44) = 2.46108e-5
CONSTS(45) = 1.77e-5
STATES(36) = 1.235996442041028e-8
CONSTS(46) = 0.4975
STATES(37) = 1.773318268499115e-5
STATES(38) = 0.9979722
STATES(39) = 0.9874435
CONSTS(47) = 9e-7
CONSTS(48) = -1
CONSTS(49) = 0.4
STATES(40) = 0.00107332
STATES(41) = 0.9960352
CONSTS(50) = 8
CONSTS(51) = 1
CONSTS(52) = 4.75
CONSTS(53) = 0.1125
CONSTS(54) = 8
CONSTS(55) = 10
CONSTS(56) = 1
CONSTS(57) = 0.28
CONSTS(58) = 0.0123
CONSTS(59) = 0.00092
CONSTS(60) = 1.7e-4
CONSTS(61) = 0.004375
CONSTS(62) = 0.75
CONSTS(63) = 0.15
CONSTS(64) = 15
STATES(42) = 1.159939
CONSTS(65) = 25
STATES(43) = 6.99469048
STATES(44) = 0.013408
STATES(45) = 0.02414815
STATES(46) = 0.02510283
CONSTS(66) = 2
CONSTS(67) = 0.2
CONSTS(68) = 0.05
CONSTS(69) = 0.07
CONSTS(70) = 0.00238
CONSTS(71) = 0.0005
CONSTS(72) = 1.124
CONSTS(73) = 0.047
CONSTS(74) = 0.00087
CONSTS(75) = 0.0087
CONSTS(76) = 10
CONSTS(77) = 0.8
STATES(47) = 20.26842
STATES(48) = 0.0018861341
CONSTS(78) = 0.05
CONSTS(79) = 0.00068
CONSTS(80) = 0.05
CONSTS(81) = 0.0895380
CONSTS(82) = 1000.00* 3.14159265358979*CONSTS(9)*CONSTS(9)*CONSTS(8)
CONSTS(83) = 0.500000* (CONSTS(4)/5.40000) ** (1.0 / 2)
CONSTS(84) = 0.0138542* (CONSTS(4)/5.40000) ** (1.0 / 2)
CONSTS(85) = CONSTS(82)*0.0600000
CONSTS(86) = 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(9)+ 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(8)
CONSTS(87) = CONSTS(86)*2.00000
CONSTS(88) = ( 1.00000*1.00000*CONSTS(87))/( CONSTS(1)*1.00000)
CONSTS(89) = CONSTS(82)*0.678000
CONSTS(90) = CONSTS(82)*0.0552000
CONSTS(91) = CONSTS(82)*0.00480000
CONSTS(92) = CONSTS(82)*0.0200000
CONSTS(93) = CONSTS(82)*0.00200000
CONSTS(94) = CONSTS(4)/(CONSTS(4)+CONSTS(21))
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
RATES(33) = (ALGBRC(7) - STATES(33))/CONSTS(42)
ALGBRC(2) = 0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(14) = 0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(2) = - 4.00000*ALGBRC(2)*STATES(2)+ ALGBRC(14)*STATES(3)
ALGBRC(3) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(15) = 0.0800000*EXP(- STATES(1)/11.0000)
RATES(27) = ALGBRC(3)*(1.00000 - STATES(27)) - ALGBRC(15)*STATES(27)
ALGBRC(4) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000))
ALGBRC(16) = TERNRY(STATES(1).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))), 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))
RATES(28) = ALGBRC(4)*(1.00000 - STATES(28)) - ALGBRC(16)*STATES(28)
ALGBRC(5) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(1)) - 6.94800e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
ALGBRC(17) = TERNRY(STATES(1).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
RATES(29) = ALGBRC(5)*(1.00000 - STATES(29)) - ALGBRC(17)*STATES(29)
ALGBRC(6) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(18) = 0.0800000*EXP(- STATES(1)/11.0000)
RATES(32) = ALGBRC(6)*(1.00000 - STATES(32)) - ALGBRC(18)*STATES(32)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
ALGBRC(19) = 1.00000/(( 0.000600000*(STATES(1) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(1) - 1.73840)))+( 0.000300000*(STATES(1)+38.3608))/(EXP( 0.152200*(STATES(1)+38.3608)) - 1.00000))
RATES(36) = (ALGBRC(8) - STATES(36))/ALGBRC(19)
ALGBRC(21) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
ALGBRC(10) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
RATES(38) = ALGBRC(21)*(1.00000 - STATES(38)) - ALGBRC(10)*STATES(38)
ALGBRC(22) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
RATES(39) = ALGBRC(22)*(1.00000 - STATES(39)) - ALGBRC(11)*STATES(39)
ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
ALGBRC(12) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
ALGBRC(27) = ALGBRC(23)/(ALGBRC(23)+ALGBRC(12))
RATES(41) = (ALGBRC(27) - STATES(41))/CONSTS(50)
ALGBRC(29) = 0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(25) = 0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(3) = - ( 3.00000*ALGBRC(2)+ALGBRC(14)+ALGBRC(25))*STATES(3)+ 4.00000*ALGBRC(2)*STATES(2)+ 2.00000*ALGBRC(14)*STATES(4)+ ALGBRC(29)*STATES(7)
RATES(4) = - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25))*STATES(4)+ 3.00000*ALGBRC(2)*STATES(3)+ 3.00000*ALGBRC(14)*STATES(5)+ ALGBRC(29)*STATES(8)
RATES(5) = - (ALGBRC(2)+ 3.00000*ALGBRC(14)+ 3.00000*ALGBRC(25))*STATES(5)+ 2.00000*ALGBRC(2)*STATES(4)+ 4.00000*ALGBRC(14)*STATES(6)+ ALGBRC(29)*STATES(9)
RATES(6) = - ( 4.00000*ALGBRC(14)+ 4.00000*ALGBRC(25))*STATES(6)+ ALGBRC(2)*STATES(5)+ ALGBRC(29)*STATES(10)
RATES(7) = - ( 3.00000*ALGBRC(2)+ALGBRC(29))*STATES(7)+ ALGBRC(14)*STATES(8)+ ALGBRC(25)*STATES(3)
RATES(8) = - ( 2.00000*ALGBRC(2)+ALGBRC(14)+ALGBRC(25)+ALGBRC(29))*STATES(8)+ 3.00000*ALGBRC(2)*STATES(7)+ 2.00000*ALGBRC(14)*STATES(9)+ 2.00000*ALGBRC(25)*STATES(4)+ 2.00000*ALGBRC(29)*STATES(11)
RATES(9) = - (ALGBRC(2)+ 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25)+ALGBRC(29))*STATES(9)+ 2.00000*ALGBRC(2)*STATES(8)+ 3.00000*ALGBRC(14)*STATES(10)+ 3.00000*ALGBRC(25)*STATES(5)+ 2.00000*ALGBRC(29)*STATES(12)
RATES(10) = - ( 3.00000*ALGBRC(14)+ 3.00000*ALGBRC(25)+ALGBRC(29))*STATES(10)+ ALGBRC(2)*STATES(9)+ 4.00000*ALGBRC(25)*STATES(6)+ 2.00000*ALGBRC(29)*STATES(13)
RATES(11) = - ( 2.00000*ALGBRC(2)+ 2.00000*ALGBRC(29))*STATES(11)+ ALGBRC(14)*STATES(12)+ ALGBRC(25)*STATES(8)
RATES(12) = - (ALGBRC(2)+ALGBRC(14)+ALGBRC(25)+ 2.00000*ALGBRC(29))*STATES(12)+ 2.00000*ALGBRC(2)*STATES(11)+ 2.00000*ALGBRC(14)*STATES(13)+ 2.00000*ALGBRC(25)*STATES(9)+ 3.00000*ALGBRC(29)*STATES(14)
RATES(13) = - ( 2.00000*ALGBRC(14)+ 2.00000*ALGBRC(25)+ 2.00000*ALGBRC(29))*STATES(13)+ ALGBRC(2)*STATES(12)+ 3.00000*ALGBRC(25)*STATES(10)+ 3.00000*ALGBRC(29)*STATES(15)
RATES(14) = - (ALGBRC(2)+ 3.00000*ALGBRC(29))*STATES(14)+ ALGBRC(14)*STATES(15)+ ALGBRC(25)*STATES(12)
RATES(15) = - (ALGBRC(14)+ALGBRC(25)+ 3.00000*ALGBRC(29))*STATES(15)+ ALGBRC(2)*STATES(14)+ 2.00000*ALGBRC(25)*STATES(13)+ 4.00000*ALGBRC(29)*STATES(16)
ALGBRC(9) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
ALGBRC(20) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
ALGBRC(26) = 1.00000/(ALGBRC(9)+ALGBRC(20))
ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
RATES(37) = (ALGBRC(30) - STATES(37))/ALGBRC(26)
ALGBRC(32) = 0.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(16) = - ( 4.00000*ALGBRC(29)+CONSTS(81))*STATES(16)+ ALGBRC(25)*STATES(15)+ ALGBRC(32)*STATES(17)
ALGBRC(34) = 0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(36) = 0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(17) = - (ALGBRC(32)+ALGBRC(36))*STATES(17)+ ALGBRC(34)*STATES(18)+ CONSTS(81)*STATES(16)
RATES(18) = - ALGBRC(34)*STATES(18)+ ALGBRC(36)*STATES(17)
ALGBRC(57) = (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(35))
ALGBRC(58) = CONSTS(40)*(STATES(1) - ALGBRC(57))
ALGBRC(51) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
ALGBRC(63) = ( CONSTS(44)*ALGBRC(51) - ALGBRC(57) ** 4.00000)/(ALGBRC(51) - ALGBRC(57) ** 4.00000+87.8251 ** 4.00000)
ALGBRC(54) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(34))
ALGBRC(64) = ( CONSTS(45)*(ALGBRC(54) - ALGBRC(57)))/((ALGBRC(54)+87.8251) - ALGBRC(57))
ALGBRC(73) = (STATES(47) - STATES(35))/CONSTS(67)
RATES(35) = - (( ALGBRC(58)*CONSTS(88))/( - 1.00000*CONSTS(89)) - (ALGBRC(63)+ALGBRC(64)+( ALGBRC(73)*CONSTS(92))/CONSTS(89)))
ALGBRC(69) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(35) - CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(70) = 1.00000 - 1.00000/(1.00000+STATES(40)/CONSTS(49) ** 2.00000)
ALGBRC(71) = ALGBRC(69)*STATES(41)*ALGBRC(70)
RATES(47) = - (( ALGBRC(71)*CONSTS(88))/( - 1.00000*CONSTS(92))+ALGBRC(73))
ALGBRC(1) = INT(VOI/CONSTS(11))*CONSTS(11)
ALGBRC(13) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
ALGBRC(59) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(54)) - 59.2150)))
ALGBRC(60) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(54))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(54)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(54))+4.75300)))
ALGBRC(61) = ALGBRC(59)/(ALGBRC(59)+ALGBRC(60))
ALGBRC(62) = CONSTS(83)*ALGBRC(61)*(STATES(1) - ALGBRC(54))
ALGBRC(48) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
ALGBRC(55) = CONSTS(33)*ALGBRC(48)*(STATES(1) - ALGBRC(54))
ALGBRC(75) = ((CONSTS(68)+CONSTS(69)) - STATES(44))+CONSTS(71)+CONSTS(70)
ALGBRC(76) = ( CONSTS(71)*CONSTS(70) - STATES(44)*(CONSTS(71)+CONSTS(70)))+ CONSTS(69)*CONSTS(70)+ CONSTS(68)*CONSTS(71)
ALGBRC(77) = - CONSTS(71)*CONSTS(70)*STATES(44)
ALGBRC(78) = (( 2.00000* (ALGBRC(75) ** 2.00000 - 3.00000*ALGBRC(76)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(75)*ALGBRC(76) - 2.00000*ALGBRC(75) ** 3.00000) - 27.0000*ALGBRC(77))/( 2.00000*ALGBRC(75) ** 2.00000 - 3.00000*ALGBRC(76) ** 1.50000))/3.00000) - ALGBRC(75)/3.00000
ALGBRC(79) = 0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(78) ** 1.40000))
ALGBRC(24) = STATES(17)+STATES(18)
ALGBRC(56) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(34)+ CONSTS(43)*STATES(30)))
ALGBRC(80) = ALGBRC(79)*ALGBRC(24)*(STATES(1) - ALGBRC(56))
ALGBRC(65) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
ALGBRC(66) = CONSTS(84)*STATES(36)*ALGBRC(65)*(STATES(1) - ALGBRC(54))
ALGBRC(42) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(43) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
ALGBRC(44) = STATES(30)/(STATES(30)+CONSTS(20)) ** 3.00000
ALGBRC(45) = CONSTS(22)*ALGBRC(43)*CONSTS(94)*ALGBRC(44)
ALGBRC(67) = EXP(STATES(1)/550.000)
ALGBRC(68) = CONSTS(46)*STATES(37) ** 3.00000*STATES(38)*STATES(39)*ALGBRC(67)*(STATES(1) - ALGBRC(54))
RATES(34) = - (( ((ALGBRC(80)+ALGBRC(66)+ALGBRC(62)+ALGBRC(68)+ALGBRC(55)+ALGBRC(13)) - 2.00000*ALGBRC(45))*CONSTS(88))/CONSTS(89) - ALGBRC(64))
ALGBRC(28) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
ALGBRC(33) = ALGBRC(31)/ALGBRC(28)
ALGBRC(35) = (1.00000 - ALGBRC(31))/ALGBRC(28)
ALGBRC(89) = ((CONSTS(72)+CONSTS(73)) - STATES(46))+CONSTS(74)+CONSTS(75)
ALGBRC(90) = ( CONSTS(75)*CONSTS(74) - STATES(46)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
ALGBRC(91) = - CONSTS(74)*CONSTS(75)*STATES(46)
ALGBRC(95) = (( 2.00000* (ALGBRC(89) ** 2.00000 - 3.00000*ALGBRC(90)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(89)*ALGBRC(90) - 2.00000*ALGBRC(89) ** 3.00000) - 27.0000*ALGBRC(91))/( 2.00000*ALGBRC(89) ** 2.00000 - 3.00000*ALGBRC(90) ** 1.50000))/3.00000) - ALGBRC(89)/3.00000
ALGBRC(106) = 3.00000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
ALGBRC(38) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+1.00000/( 30.5530*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070))))
ALGBRC(37) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
ALGBRC(39) = ALGBRC(37)/ALGBRC(38)
ALGBRC(40) = (1.00000 - ALGBRC(37))/ALGBRC(38)
RATES(19) = - (ALGBRC(33)+ALGBRC(106)+ALGBRC(40))*STATES(19)+ ALGBRC(35)*STATES(20)+ CONSTS(15)*STATES(21)+ ALGBRC(39)*STATES(23)
RATES(20) = - (ALGBRC(35)+ALGBRC(106)+ALGBRC(40))*STATES(20)+ ALGBRC(33)*STATES(19)+ CONSTS(15)*STATES(22)+ ALGBRC(39)*STATES(24)
ALGBRC(96) = 25.0000 - 17.5000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
ALGBRC(98) = 1.00000/( ALGBRC(96)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
ALGBRC(100) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(98))
ALGBRC(41) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
ALGBRC(102) = ALGBRC(41)/ALGBRC(100)
ALGBRC(104) = (1.00000 - ALGBRC(41))/ALGBRC(100)
RATES(21) = - (ALGBRC(33)+CONSTS(15)+ALGBRC(104))*STATES(21)+ ALGBRC(106)*STATES(19)+ ALGBRC(35)*STATES(22)+ ALGBRC(102)*STATES(25)
RATES(22) = - (ALGBRC(35)+CONSTS(15)+ALGBRC(104))*STATES(22)+ ALGBRC(106)*STATES(20)+ ALGBRC(33)*STATES(21)+ ALGBRC(102)*STATES(26)
ALGBRC(108) = ( CONSTS(16)*ALGBRC(39)*ALGBRC(104)*ALGBRC(106))/( ALGBRC(40)*ALGBRC(102)*CONSTS(15))
RATES(23) = - (ALGBRC(33)+ALGBRC(108)+ALGBRC(39))*STATES(23)+ ALGBRC(40)*STATES(19)+ CONSTS(16)*STATES(25)+ ALGBRC(35)*STATES(24)
RATES(24) = - (ALGBRC(35)+ALGBRC(108)+ALGBRC(39))*STATES(24)+ ALGBRC(40)*STATES(20)+ CONSTS(16)*STATES(26)+ ALGBRC(33)*STATES(23)
RATES(25) = - (ALGBRC(33)+CONSTS(16)+ALGBRC(102))*STATES(25)+ ALGBRC(108)*STATES(23)+ ALGBRC(104)*STATES(21)+ ALGBRC(35)*STATES(26)
RATES(26) = - (ALGBRC(35)+CONSTS(16)+ALGBRC(102))*STATES(26)+ ALGBRC(108)*STATES(24)+ ALGBRC(104)*STATES(22)+ ALGBRC(33)*STATES(25)
ALGBRC(52) = CONSTS(19)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(51))
ALGBRC(82) = 0.800000*CONSTS(26)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*ALGBRC(78)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(83) = 1.00000+CONSTS(25)/ALGBRC(78) ** 2.00000
ALGBRC(46) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(84) = CONSTS(32)*STATES(30) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(78)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(78)/CONSTS(31))
ALGBRC(85) = CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(29) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(78)
ALGBRC(86) = ALGBRC(82)/( ALGBRC(83)*ALGBRC(46)*(ALGBRC(84)+ALGBRC(85)))
ALGBRC(53) = CONSTS(41)*STATES(32) ** 3.00000*STATES(33)*(STATES(1) - ALGBRC(51))
ALGBRC(49) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(50) = ( CONSTS(1)*CONSTS(39)*ALGBRC(49)*( STATES(30)*EXP(ALGBRC(49)) - CONSTS(6)))/(EXP(ALGBRC(49)) - 1.00000)
ALGBRC(113) = (STATES(31) - STATES(30))/CONSTS(67)
RATES(30) = - (( ( 3.00000*ALGBRC(86)+ 3.00000*ALGBRC(45)+ALGBRC(52)+ALGBRC(53)+ALGBRC(50))*CONSTS(88))/CONSTS(89) - (ALGBRC(63)+( ALGBRC(113)*CONSTS(92))/CONSTS(89)))
ALGBRC(92) = ((CONSTS(72)+CONSTS(73)) - STATES(45))+CONSTS(74)+CONSTS(75)
ALGBRC(93) = ( CONSTS(75)*CONSTS(74) - STATES(45)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
ALGBRC(94) = - CONSTS(74)*CONSTS(75)*STATES(45)
ALGBRC(101) = (( 2.00000* (ALGBRC(92) ** 2.00000 - 3.00000*ALGBRC(93)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(92)*ALGBRC(93) - 2.00000*ALGBRC(92) ** 3.00000) - 27.0000*ALGBRC(94))/( 2.00000*ALGBRC(92) ** 2.00000 - 3.00000*ALGBRC(93) ** 1.50000))/3.00000) - ALGBRC(92)/3.00000
ALGBRC(103) = 0.200000*CONSTS(26)*( STATES(31) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*ALGBRC(101)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(105) = 1.00000+CONSTS(25)/ALGBRC(101) ** 2.00000
ALGBRC(47) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(107) = CONSTS(32)*STATES(31) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(101)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(101)/CONSTS(31))
ALGBRC(109) = CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(31)/CONSTS(29) ** 3.00000)+ STATES(31) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(101)
ALGBRC(110) = ALGBRC(103)/( ALGBRC(105)*ALGBRC(47)*(ALGBRC(107)+ALGBRC(109)))
RATES(31) = - (( 3.00000*ALGBRC(110)*CONSTS(88))/CONSTS(92)+ALGBRC(113))
ALGBRC(97) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*ALGBRC(95)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(18)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(99) = ALGBRC(97)*(STATES(20)+STATES(22))
ALGBRC(87) = ( CONSTS(34)*ALGBRC(78))/(CONSTS(35)+ALGBRC(78))
ALGBRC(88) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*ALGBRC(78)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(38)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(111) = (ALGBRC(99)+ALGBRC(88)+ALGBRC(87)) - 2.00000*(ALGBRC(86)+ALGBRC(110))
ALGBRC(114) = ALGBRC(52)+ 3.00000*(ALGBRC(86)+ALGBRC(110))+ 3.00000*ALGBRC(45)+ALGBRC(53)+ALGBRC(50)
ALGBRC(81) = ALGBRC(66)+ALGBRC(80)+ALGBRC(62)+ALGBRC(55)+ - 2.00000*ALGBRC(45)+ALGBRC(68)+ALGBRC(13)
ALGBRC(72) = ALGBRC(58)+ALGBRC(71)
RATES(1) = - (ALGBRC(114)+ALGBRC(81)+ALGBRC(111)+ALGBRC(72))
ALGBRC(112) = (ALGBRC(101) - ALGBRC(78))/CONSTS(67)
ALGBRC(116) = (ALGBRC(101) - ALGBRC(95))/CONSTS(66)
RATES(45) = - ((ALGBRC(112)+ALGBRC(116)) - (( 2.00000*ALGBRC(110)*CONSTS(88))/( 2.00000*CONSTS(92))+( STATES(40)*CONSTS(91))/CONSTS(92)))
RATES(46) = - (( ALGBRC(99)*CONSTS(88))/( 2.00000*CONSTS(93)) - ( ALGBRC(116)*CONSTS(92))/CONSTS(93))
ALGBRC(115) = (CONSTS(76)+CONSTS(77)) - STATES(43)
ALGBRC(117) = STATES(43)*CONSTS(77)
ALGBRC(118) = - ALGBRC(115)/2.00000+ (ALGBRC(115) ** 2.00000+ 4.00000*ALGBRC(117)) ** (1.0 / 2)/2.00000
ALGBRC(119) = (STATES(42) - ALGBRC(118))/CONSTS(65)
RATES(43) = ALGBRC(119) - STATES(40)
ALGBRC(120) = ( CONSTS(80)*(1.00000 - STATES(48))*1.00000)/(1.00000+CONSTS(63)/ALGBRC(101))
RATES(48) = CONSTS(78)*ALGBRC(120)*(ALGBRC(120) - STATES(48)) - CONSTS(79)*STATES(48)
ALGBRC(124) = ( CONSTS(60)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
ALGBRC(122) = ( CONSTS(62)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
ALGBRC(126) = ( (ALGBRC(122)+1.00000)*CONSTS(61)*ALGBRC(78))/((ALGBRC(78)+CONSTS(59)) - ALGBRC(124))
ALGBRC(74) = (0.00437500/CONSTS(64))*STATES(42)
RATES(44) = (( - ((ALGBRC(88)+ALGBRC(87)) - 2.00000*ALGBRC(86))*CONSTS(88))/( 2.00000*CONSTS(89))+( (ALGBRC(126) - ALGBRC(74))*CONSTS(90))/CONSTS(89)) - ( ALGBRC(112)*CONSTS(92))/CONSTS(89)
RATES(42) = ALGBRC(126) - (ALGBRC(74)+( ALGBRC(119)*CONSTS(91))/CONSTS(90))
ALGBRC(121) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(120)) ** CONSTS(55))
ALGBRC(123) = CONSTS(52)*(1.00000+ALGBRC(121))
ALGBRC(125) = ALGBRC(123)*CONSTS(53)
ALGBRC(128) = ( ALGBRC(99)*ALGBRC(125))/(1.00000+CONSTS(51)/ALGBRC(118) ** CONSTS(54))
ALGBRC(127) = ALGBRC(123)/(1.00000+CONSTS(58)/ALGBRC(118))
RATES(40) = - (ALGBRC(128)+STATES(40))/ALGBRC(127)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+91.0000)/6.10000))
ALGBRC(2) = 0.0148640*EXP(( 0.0298770*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(14) = 0.0839860*EXP(( - 0.0554610*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(3) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(15) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(4) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000))
ALGBRC(16) = TERNRY(STATES(1).GE.- 40.0000, 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))), 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)))
ALGBRC(5) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 127140.*EXP( 0.244400*STATES(1)) - 6.94800e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
ALGBRC(17) = TERNRY(STATES(1).GE.- 40.0000, ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
ALGBRC(6) = ( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(18) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0850)/4.25000))
ALGBRC(19) = 1.00000/(( 0.000600000*(STATES(1) - 1.73840))/(1.00000 - EXP( - 0.136000*(STATES(1) - 1.73840)))+( 0.000300000*(STATES(1)+38.3608))/(EXP( 0.152200*(STATES(1)+38.3608)) - 1.00000))
ALGBRC(21) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
ALGBRC(10) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
ALGBRC(22) = 1.00000/( 250.000*(1.00000+EXP((STATES(1)+60.0000)/5.00000)))
ALGBRC(11) = 1.00000/( 9.79530*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
ALGBRC(23) = 0.0250000/(1.00000+EXP((STATES(1)+58.0000)/5.00000))
ALGBRC(12) = 1.00000/( 5.00000*(1.00000+EXP((STATES(1)+19.0000)/- 9.00000)))
ALGBRC(27) = ALGBRC(23)/(ALGBRC(23)+ALGBRC(12))
ALGBRC(29) = 0.00311730*EXP(( - 0.426250*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(25) = 0.0146010*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(9) = 1.00000/( 1.20890*(1.00000+EXP((STATES(1) - 18.4099)/- 29.3814)))
ALGBRC(20) = 3.50000/(1.00000+EXP((STATES(1)+100.000)/29.3814))
ALGBRC(26) = 1.00000/(ALGBRC(9)+ALGBRC(20))
ALGBRC(30) = 1.00000/(1.00000+EXP((STATES(1)+9.43700)/- 7.13300))
ALGBRC(32) = 0.0773200*EXP(( - 0.0647260*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(34) = 0.794050*EXP(( - 0.0801740*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(36) = 0.586380*EXP(( 0.282060*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(57) = (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(35))
ALGBRC(58) = CONSTS(40)*(STATES(1) - ALGBRC(57))
ALGBRC(51) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
ALGBRC(63) = ( CONSTS(44)*ALGBRC(51) - ALGBRC(57) ** 4.00000)/(ALGBRC(51) - ALGBRC(57) ** 4.00000+87.8251 ** 4.00000)
ALGBRC(54) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(34))
ALGBRC(64) = ( CONSTS(45)*(ALGBRC(54) - ALGBRC(57)))/((ALGBRC(54)+87.8251) - ALGBRC(57))
ALGBRC(73) = (STATES(47) - STATES(35))/CONSTS(67)
ALGBRC(69) = ( (( CONSTS(47)*CONSTS(48) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(35) - CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(70) = 1.00000 - 1.00000/(1.00000+STATES(40)/CONSTS(49) ** 2.00000)
ALGBRC(71) = ALGBRC(69)*STATES(41)*ALGBRC(70)
ALGBRC(1) = INT(VOI/CONSTS(11))*CONSTS(11)
ALGBRC(13) = TERNRY(VOI - ALGBRC(1).GE.CONSTS(10).AND.VOI - ALGBRC(1).LE.CONSTS(10)+CONSTS(12), CONSTS(13), 0.00000)
ALGBRC(59) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(54)) - 59.2150)))
ALGBRC(60) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(54))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(54)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(54))+4.75300)))
ALGBRC(61) = ALGBRC(59)/(ALGBRC(59)+ALGBRC(60))
ALGBRC(62) = CONSTS(83)*ALGBRC(61)*(STATES(1) - ALGBRC(54))
ALGBRC(48) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
ALGBRC(55) = CONSTS(33)*ALGBRC(48)*(STATES(1) - ALGBRC(54))
ALGBRC(75) = ((CONSTS(68)+CONSTS(69)) - STATES(44))+CONSTS(71)+CONSTS(70)
ALGBRC(76) = ( CONSTS(71)*CONSTS(70) - STATES(44)*(CONSTS(71)+CONSTS(70)))+ CONSTS(69)*CONSTS(70)+ CONSTS(68)*CONSTS(71)
ALGBRC(77) = - CONSTS(71)*CONSTS(70)*STATES(44)
ALGBRC(78) = (( 2.00000* (ALGBRC(75) ** 2.00000 - 3.00000*ALGBRC(76)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(75)*ALGBRC(76) - 2.00000*ALGBRC(75) ** 3.00000) - 27.0000*ALGBRC(77))/( 2.00000*ALGBRC(75) ** 2.00000 - 3.00000*ALGBRC(76) ** 1.50000))/3.00000) - ALGBRC(75)/3.00000
ALGBRC(79) = 0.0826000*(1.00000+0.600000/(1.00000+3.80000e-05/ALGBRC(78) ** 1.40000))
ALGBRC(24) = STATES(17)+STATES(18)
ALGBRC(56) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(43)*CONSTS(6))/(STATES(34)+ CONSTS(43)*STATES(30)))
ALGBRC(80) = ALGBRC(79)*ALGBRC(24)*(STATES(1) - ALGBRC(56))
ALGBRC(65) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
ALGBRC(66) = CONSTS(84)*STATES(36)*ALGBRC(65)*(STATES(1) - ALGBRC(54))
ALGBRC(42) = ( CONSTS(23)*(STATES(1) - CONSTS(24))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(43) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
ALGBRC(44) = STATES(30)/(STATES(30)+CONSTS(20)) ** 3.00000
ALGBRC(45) = CONSTS(22)*ALGBRC(43)*CONSTS(94)*ALGBRC(44)
ALGBRC(67) = EXP(STATES(1)/550.000)
ALGBRC(68) = CONSTS(46)*STATES(37) ** 3.00000*STATES(38)*STATES(39)*ALGBRC(67)*(STATES(1) - ALGBRC(54))
ALGBRC(28) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
ALGBRC(31) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
ALGBRC(33) = ALGBRC(31)/ALGBRC(28)
ALGBRC(35) = (1.00000 - ALGBRC(31))/ALGBRC(28)
ALGBRC(89) = ((CONSTS(72)+CONSTS(73)) - STATES(46))+CONSTS(74)+CONSTS(75)
ALGBRC(90) = ( CONSTS(75)*CONSTS(74) - STATES(46)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
ALGBRC(91) = - CONSTS(74)*CONSTS(75)*STATES(46)
ALGBRC(95) = (( 2.00000* (ALGBRC(89) ** 2.00000 - 3.00000*ALGBRC(90)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(89)*ALGBRC(90) - 2.00000*ALGBRC(89) ** 3.00000) - 27.0000*ALGBRC(91))/( 2.00000*ALGBRC(89) ** 2.00000 - 3.00000*ALGBRC(90) ** 1.50000))/3.00000) - ALGBRC(89)/3.00000
ALGBRC(106) = 3.00000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
ALGBRC(38) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+1.00000/( 30.5530*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070))))
ALGBRC(37) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
ALGBRC(39) = ALGBRC(37)/ALGBRC(38)
ALGBRC(40) = (1.00000 - ALGBRC(37))/ALGBRC(38)
ALGBRC(96) = 25.0000 - 17.5000/(1.00000+0.00300000/ALGBRC(95) ** 4.00000)
ALGBRC(98) = 1.00000/( ALGBRC(96)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
ALGBRC(100) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(98))
ALGBRC(41) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
ALGBRC(102) = ALGBRC(41)/ALGBRC(100)
ALGBRC(104) = (1.00000 - ALGBRC(41))/ALGBRC(100)
ALGBRC(108) = ( CONSTS(16)*ALGBRC(39)*ALGBRC(104)*ALGBRC(106))/( ALGBRC(40)*ALGBRC(102)*CONSTS(15))
ALGBRC(52) = CONSTS(19)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(51))
ALGBRC(82) = 0.800000*CONSTS(26)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*ALGBRC(78)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(83) = 1.00000+CONSTS(25)/ALGBRC(78) ** 2.00000
ALGBRC(46) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(84) = CONSTS(32)*STATES(30) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(78)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(78)/CONSTS(31))
ALGBRC(85) = CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(29) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(78)
ALGBRC(86) = ALGBRC(82)/( ALGBRC(83)*ALGBRC(46)*(ALGBRC(84)+ALGBRC(85)))
ALGBRC(53) = CONSTS(41)*STATES(32) ** 3.00000*STATES(33)*(STATES(1) - ALGBRC(51))
ALGBRC(49) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(50) = ( CONSTS(1)*CONSTS(39)*ALGBRC(49)*( STATES(30)*EXP(ALGBRC(49)) - CONSTS(6)))/(EXP(ALGBRC(49)) - 1.00000)
ALGBRC(113) = (STATES(31) - STATES(30))/CONSTS(67)
ALGBRC(92) = ((CONSTS(72)+CONSTS(73)) - STATES(45))+CONSTS(74)+CONSTS(75)
ALGBRC(93) = ( CONSTS(75)*CONSTS(74) - STATES(45)*(CONSTS(74)+CONSTS(75)))+ CONSTS(73)*CONSTS(75)+ CONSTS(72)*CONSTS(74)
ALGBRC(94) = - CONSTS(74)*CONSTS(75)*STATES(45)
ALGBRC(101) = (( 2.00000* (ALGBRC(92) ** 2.00000 - 3.00000*ALGBRC(93)) ** (1.0 / 2))/3.00000)*cos(acos((( 9.00000*ALGBRC(92)*ALGBRC(93) - 2.00000*ALGBRC(92) ** 3.00000) - 27.0000*ALGBRC(94))/( 2.00000*ALGBRC(92) ** 2.00000 - 3.00000*ALGBRC(93) ** 1.50000))/3.00000) - ALGBRC(92)/3.00000
ALGBRC(103) = 0.200000*CONSTS(26)*( STATES(31) ** 3.00000*CONSTS(5)*EXP(( CONSTS(28)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*ALGBRC(101)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(105) = 1.00000+CONSTS(25)/ALGBRC(101) ** 2.00000
ALGBRC(47) = 1.00000+ CONSTS(27)*EXP(( (CONSTS(28) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(107) = CONSTS(32)*STATES(31) ** 3.00000+ CONSTS(30) ** 3.00000*ALGBRC(101)+ CONSTS(29) ** 3.00000*CONSTS(5)*(1.00000+ALGBRC(101)/CONSTS(31))
ALGBRC(109) = CONSTS(31)*CONSTS(6) ** 3.00000*(1.00000+STATES(31)/CONSTS(29) ** 3.00000)+ STATES(31) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*ALGBRC(101)
ALGBRC(110) = ALGBRC(103)/( ALGBRC(105)*ALGBRC(47)*(ALGBRC(107)+ALGBRC(109)))
ALGBRC(97) = ( (( CONSTS(14)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(17)*ALGBRC(95)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(18)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(99) = ALGBRC(97)*(STATES(20)+STATES(22))
ALGBRC(87) = ( CONSTS(34)*ALGBRC(78))/(CONSTS(35)+ALGBRC(78))
ALGBRC(88) = ( (( CONSTS(36)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(37)*ALGBRC(78)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(38)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(111) = (ALGBRC(99)+ALGBRC(88)+ALGBRC(87)) - 2.00000*(ALGBRC(86)+ALGBRC(110))
ALGBRC(114) = ALGBRC(52)+ 3.00000*(ALGBRC(86)+ALGBRC(110))+ 3.00000*ALGBRC(45)+ALGBRC(53)+ALGBRC(50)
ALGBRC(81) = ALGBRC(66)+ALGBRC(80)+ALGBRC(62)+ALGBRC(55)+ - 2.00000*ALGBRC(45)+ALGBRC(68)+ALGBRC(13)
ALGBRC(72) = ALGBRC(58)+ALGBRC(71)
ALGBRC(112) = (ALGBRC(101) - ALGBRC(78))/CONSTS(67)
ALGBRC(116) = (ALGBRC(101) - ALGBRC(95))/CONSTS(66)
ALGBRC(115) = (CONSTS(76)+CONSTS(77)) - STATES(43)
ALGBRC(117) = STATES(43)*CONSTS(77)
ALGBRC(118) = - ALGBRC(115)/2.00000+ (ALGBRC(115) ** 2.00000+ 4.00000*ALGBRC(117)) ** (1.0 / 2)/2.00000
ALGBRC(119) = (STATES(42) - ALGBRC(118))/CONSTS(65)
ALGBRC(120) = ( CONSTS(80)*(1.00000 - STATES(48))*1.00000)/(1.00000+CONSTS(63)/ALGBRC(101))
ALGBRC(124) = ( CONSTS(60)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
ALGBRC(122) = ( CONSTS(62)*ALGBRC(120)*1.00000)/(CONSTS(63)+ ALGBRC(120)*1.00000)
ALGBRC(126) = ( (ALGBRC(122)+1.00000)*CONSTS(61)*ALGBRC(78))/((ALGBRC(78)+CONSTS(59)) - ALGBRC(124))
ALGBRC(74) = (0.00437500/CONSTS(64))*STATES(42)
ALGBRC(121) = CONSTS(56)/(1.00000+CONSTS(57)/( 1.00000*ALGBRC(120)) ** CONSTS(55))
ALGBRC(123) = CONSTS(52)*(1.00000+ALGBRC(121))
ALGBRC(125) = ALGBRC(123)*CONSTS(53)
ALGBRC(128) = ( ALGBRC(99)*ALGBRC(125))/(1.00000+CONSTS(51)/ALGBRC(118) ** CONSTS(54))
ALGBRC(127) = ALGBRC(123)/(1.00000+CONSTS(58)/ALGBRC(118))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END