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 120 entries in the algebraic variable array.
C There are a total of 46 entries in each of the rate and state variable arrays.
C There are a total of 99 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(94) is Vmyo in component cell_geometry (uL).
C CONSTS(95) is Vnsr in component cell_geometry (uL).
C CONSTS(90) is Vsr in component cell_geometry (uL).
C CONSTS(96) is Vjsr in component cell_geometry (uL).
C CONSTS(97) is Vss_sr in component cell_geometry (uL).
C CONSTS(98) is Vss_CaL in component cell_geometry (uL).
C CONSTS(92) is Acap in component cell_geometry (cm2).
C CONSTS(93) 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(87) is vcell in component cell_geometry (uL).
C CONSTS(91) 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(74) is caiont in component membrane (uA_per_uF).
C ALGBRC(91) is naiont in component membrane (uA_per_uF).
C ALGBRC(115) is kiont in component membrane (uA_per_uF).
C ALGBRC(119) is clont in component membrane (uA_per_uF).
C ALGBRC(84) is INa in component INa (uA_per_uF).
C ALGBRC(37) is ICaL in component ICaL (uA_per_uF).
C ALGBRC(107) is IK1 in component IK1 (uA_per_uF).
C ALGBRC(97) is IKp in component IKp (uA_per_uF).
C ALGBRC(101) is IKs in component IKs (uA_per_uF).
C ALGBRC(112) is IKr in component IKr (uA_per_uF).
C ALGBRC(72) is IpCa in component IpCa (uA_per_uF).
C ALGBRC(73) is ICab in component ICab (uA_per_uF).
C ALGBRC(70) is INaCa in component INaCa (uA_per_uF).
C ALGBRC(68) is INaCa_ss_sr in component INaCa (uA_per_uF).
C ALGBRC(51) is INaK in component INaK (uA_per_uF).
C ALGBRC(114) is Ito1 in component Ito1 (uA_per_uF).
C ALGBRC(118) is Ito2 in component Ito2 (uA_per_uF).
C ALGBRC(103) is IClb in component IClb (uA_per_uF).
C ALGBRC(88) is INaL in component INaL (uA_per_uF).
C ALGBRC(78) is INab in component INab (uA_per_uF).
C ALGBRC(24) is gKs in component IKs (mS_per_uF).
C CONSTS(14) is gKs_max 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 ALGBRC(28) is O2 in component IKs (dimensionless).
C ALGBRC(31) is OKs in component IKs (dimensionless).
C ALGBRC(32) is alpha in component IKs (per_ms).
C ALGBRC(35) is beta in component IKs (per_ms).
C ALGBRC(40) is delta in component IKs (per_ms).
C ALGBRC(38) is gamma in component IKs (per_ms).
C ALGBRC(43) is eta in component IKs (per_ms).
C CONSTS(86) is theta in component IKs (per_ms).
C ALGBRC(46) is omega in component IKs (per_ms).
C ALGBRC(49) is psi in component IKs (per_ms).
C ALGBRC(99) is EKs in component reversal_potentials (mV).
C STATES(18) is Ca_i in component Ca (mM).
C ALGBRC(34) is ICaL_max in component ICaL (uA_per_uF).
C CONSTS(15) 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 ALGBRC(2) is OI_star in component ICaL (dimensionless).
C ALGBRC(14) is ACT_tau in component ICaL (ms).
C ALGBRC(39) is IV_tau in component ICaL (ms).
C ALGBRC(54) is IV_tau_star in component ICaL (ms).
C ALGBRC(25) is ACT_infinity in component ICaL (dimensionless).
C ALGBRC(36) is IV_infinity in component ICaL (dimensionless).
C ALGBRC(47) is IV_infinity_star in component ICaL (dimensionless).
C ALGBRC(52) is IV_beta_star in component ICaL (dimensionless).
C ALGBRC(50) is IV_beta_infinity_star in component ICaL (dimensionless).
C ALGBRC(29) is alpha in component ICaL (per_ms).
C ALGBRC(33) is beta in component ICaL (per_ms).
C ALGBRC(60) is delta in component ICaL (per_ms).
C ALGBRC(41) is x in component ICaL (per_ms).
C ALGBRC(44) is y in component ICaL (per_ms).
C CONSTS(16) is theta in component ICaL (per_ms).
C ALGBRC(62) is delta_I in component ICaL (per_ms).
C ALGBRC(56) is x_star in component ICaL (per_ms).
C ALGBRC(58) is y_star in component ICaL (per_ms).
C CONSTS(17) is theta_I in component ICaL (per_ms).
C CONSTS(18) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(19) is gamma_Cao in component ICaL (dimensionless).
C STATES(26) is Ca_ss_CaL in component Ca (mM).
C CONSTS(20) is g_Na in component INa (mS_per_uF).
C ALGBRC(81) 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(21) is kmnai in component INaK (mM).
C CONSTS(22) is kmko in component INaK (mM).
C CONSTS(23) is ibarnak in component INaK (uA_per_uF).
C CONSTS(24) is sigma in component INaK (dimensionless).
C CONSTS(25) is V_half in component INaK (mV).
C ALGBRC(42) is phi in component INaK (dimensionless).
C ALGBRC(45) is fv in component INaK (dimensionless).
C CONSTS(99) is PK in component INaK (dimensionless).
C ALGBRC(48) is PNa in component INaK (dimensionless).
C STATES(30) is Na_i in component Na (mM).
C ALGBRC(69) is INaCa_cai in component INaCa (uA_per_uF).
C CONSTS(26) is KmCa_act in component INaCa (mM).
C CONSTS(27) is Vmax in component INaCa (uA_per_uF).
C CONSTS(28) is ksat in component INaCa (dimensionless).
C CONSTS(29) is eta in component INaCa (dimensionless).
C CONSTS(30) is KmNai in component INaCa (mM).
C CONSTS(31) is KmNao in component INaCa (mM).
C CONSTS(32) is KmCai in component INaCa (mM).
C CONSTS(33) is KmCao in component INaCa (mM).
C STATES(31) is Ca_ss_sr in component Ca (mM).
C STATES(32) is Na_ss_sr in component Na (mM).
C ALGBRC(53) is numerator in component INaCa (mM4_uA_per_uF).
C ALGBRC(63) is num_ss in component INaCa (mM4_uA_per_uF).
C ALGBRC(55) is denom_1 in component INaCa (dimensionless).
C ALGBRC(57) is denom_2 in component INaCa (dimensionless).
C ALGBRC(59) is denom_3 in component INaCa (mM4).
C ALGBRC(61) is denom_4 in component INaCa (mM4).
C ALGBRC(64) is denom_ss_1 in component INaCa (dimensionless).
C ALGBRC(65) is denom_ss_2 in component INaCa (dimensionless).
C ALGBRC(66) is denom_ss_3 in component INaCa (mM4).
C ALGBRC(67) is denom_ss_4 in component INaCa (mM4).
C CONSTS(34) is gKp in component IKp (mS_per_uF).
C ALGBRC(71) is Kp in component IKp (dimensionless).
C ALGBRC(94) is EK in component reversal_potentials (mV).
C CONSTS(35) is gpCa in component IpCa (uA_per_uF).
C CONSTS(36) is KmpCa in component IpCa (mM).
C CONSTS(37) is PCab in component ICab (L_per_F_ms).
C CONSTS(38) is gamma_Ca_i in component ICab (dimensionless).
C CONSTS(39) is gamma_Ca_o in component ICab (dimensionless).
C CONSTS(40) is PNab in component INab (L_per_F_ms).
C ALGBRC(76) is phi in component INab (dimensionless).
C CONSTS(41) is gClb in component IClb (mS_per_uF).
C ALGBRC(102) is ECl in component reversal_potentials (mV).
C CONSTS(42) is gNaL in component INaL (mS_per_uF).
C STATES(33) is mL in component INaL_mL_gate (dimensionless).
C STATES(34) 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(43) is tau_hL in component INaL_hL_gate (ms).
C ALGBRC(7) is hL_infinity in component INaL_hL_gate (dimensionless).
C CONSTS(44) is prnak in component reversal_potentials (dimensionless).
C STATES(35) is K_i in component K (mM).
C STATES(36) is Cl_i in component Cl (mM).
C CONSTS(45) is g_K1_max in component IK1 (mS_per_uF).
C CONSTS(88) is gK1 in component IK1 (mS_per_uF).
C ALGBRC(106) is K1 in component IK1_K1_gate (dimensionless).
C ALGBRC(104) is alpha_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(105) is beta_k1 in component IK1_K1_gate (dimensionless).
C ALGBRC(108) is CT_Na_Cl in component CT_Na_Cl (mM_per_ms).
C CONSTS(46) is CT_Na_Cl_max in component CT_Na_Cl (mM_per_ms).
C ALGBRC(109) is CT_K_Cl in component CT_K_Cl (mM_per_ms).
C CONSTS(47) is CT_K_Cl_max in component CT_K_Cl (mM_per_ms).
C CONSTS(89) is gKr in component IKr (mS_per_uF).
C CONSTS(48) is gKr_max in component IKr (mS_per_uF).
C ALGBRC(111) is r in component IKr_r_gate (dimensionless).
C STATES(37) 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(49) is gto1 in component Ito1 (mS_per_uF).
C ALGBRC(113) is rto1 in component Ito1 (dimensionless).
C STATES(38) is a in component Ito1_a_gate (dimensionless).
C STATES(39) is i1f in component Ito1_i1f_gate (dimensionless).
C STATES(40) 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(50) is PCl in component Ito2 (L_per_F_ms).
C CONSTS(51) is zCl in component Ito2 (dimensionless).
C ALGBRC(116) is Ito2_max in component Ito2 (uA_per_uF).
C CONSTS(52) is Ito2_max_scaling_factor in component Ito2 (dimensionless).
C ALGBRC(117) is KCa_ito2 in component Ito2 (dimensionless).
C CONSTS(53) is kCa_ito2 in component Ito2 (mM_per_ms).
C STATES(41) is Irel in component Irel (mM_per_ms).
C STATES(42) 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(54) is tau_i2f in component Ito2_i2f_gate (ms).
C ALGBRC(100) is rel_infinity in component Irel (mM_per_ms).
C ALGBRC(98) is tau_rel in component Irel (ms).
C CONSTS(55) is Krel_infinity in component Irel (mM).
C ALGBRC(95) is alpha_rel in component Irel (mM_per_uA_per_uF_per_ms).
C ALGBRC(92) is beta_tau in component Irel (ms).
C CONSTS(56) is beta_0 in component Irel (ms).
C CONSTS(57) is kappa in component Irel (mM_per_uA_per_uF_per_ms_per_ms).
C CONSTS(58) is h_rel in component Irel (dimensionless).
C CONSTS(59) is h_beta in component Irel (dimensionless).
C ALGBRC(89) is delta_beta_CaMK in component Irel (dimensionless).
C CONSTS(60) is delta_beta_0 in component Irel (dimensionless).
C CONSTS(61) is K_beta in component Irel (mM).
C CONSTS(62) is Krel_tau in component Irel (mM).
C STATES(43) is Ca_JSR in component Ca (mM).
C ALGBRC(85) is CaMK_active in component CaMK_active (dimensionless).
C ALGBRC(96) is Iup in component Iup (mM_per_ms).
C CONSTS(63) is kmup in component Iup (mM).
C ALGBRC(93) is delta_kmPLB in component Iup (mM).
C CONSTS(64) is delta_kmPLB_bar in component Iup (mM).
C CONSTS(65) is iupbar in component Iup (mM_per_ms).
C ALGBRC(90) is delta_iupCaMK in component Iup (dimensionless).
C CONSTS(66) is delta_iupCaMK_bar in component Iup (dimensionless).
C CONSTS(67) is kmCaMK in component Iup (mM).
C ALGBRC(75) is Ileak in component Ileak (mM_per_ms).
C CONSTS(68) is nsrbar in component Ileak (mM).
C STATES(44) is Ca_NSR in component Ca (mM).
C ALGBRC(77) is Itr in component Itr (mM_per_ms).
C CONSTS(69) is tautr in component Itr (ms).
C CONSTS(70) is tau_diff in component Ca (ms).
C CONSTS(71) is tau_diff_ss in component Ca (ms).
C ALGBRC(79) is Idiff in component Ca (mM_per_ms).
C ALGBRC(83) is Idiff_ss in component Ca (mM_per_ms).
C ALGBRC(82) is bmyo in component Ca (dimensionless).
C ALGBRC(86) is bss_sr in component Ca (dimensionless).
C ALGBRC(87) is bss_cal in component Ca (dimensionless).
C ALGBRC(80) is bcsqn in component Ca (dimensionless).
C CONSTS(72) is cmdn_bar in component Ca (mM).
C CONSTS(73) is trpn_bar in component Ca (mM).
C CONSTS(74) is km_cmdn in component Ca (mM).
C CONSTS(75) is km_trpn in component Ca (mM).
C CONSTS(76) is kmcsqn in component Ca (mM).
C CONSTS(77) is csqnbar in component Ca (mM).
C CONSTS(78) is BSRmax in component Ca (mM).
C CONSTS(79) is KmBSR in component Ca (mM).
C CONSTS(80) is BSLmax in component Ca (mM).
C CONSTS(81) is KmBSL in component Ca (mM).
C ALGBRC(110) is Idiff_Na in component Na (mM_per_ms).
C STATES(45) is Cl_ss in component Cl (mM).
C ALGBRC(120) is Idiff_Cl in component Cl (mM_per_ms).
C CONSTS(82) is Km in component CaMK_active (mM).
C STATES(46) is CaMK_trap in component CaMK_active (dimensionless).
C CONSTS(83) is alpha_CaMK in component CaMK_active (per_ms).
C CONSTS(84) is beta_CaMK in component CaMK_active (per_ms).
C CONSTS(85) 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(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(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(33) is d/dt mL in component INaL_mL_gate (dimensionless).
C RATES(34) is d/dt hL in component INaL_hL_gate (dimensionless).
C RATES(37) is d/dt xr in component IKr_xr_gate (dimensionless).
C RATES(38) is d/dt a in component Ito1_a_gate (dimensionless).
C RATES(39) is d/dt i1f in component Ito1_i1f_gate (dimensionless).
C RATES(40) is d/dt i1s in component Ito1_i1s_gate (dimensionless).
C RATES(42) is d/dt i2f in component Ito2_i2f_gate (dimensionless).
C RATES(41) is d/dt Irel in component Irel (mM_per_ms).
C RATES(18) is d/dt Ca_i in component Ca (mM).
C RATES(31) is d/dt Ca_ss_sr in component Ca (mM).
C RATES(26) is d/dt Ca_ss_CaL in component Ca (mM).
C RATES(44) is d/dt Ca_NSR in component Ca (mM).
C RATES(43) is d/dt Ca_JSR in component Ca (mM).
C RATES(30) is d/dt Na_i in component Na (mM).
C RATES(32) is d/dt Na_ss_sr in component Na (mM).
C RATES(36) is d/dt Cl_i in component Cl (mM).
C RATES(45) is d/dt Cl_ss in component Cl (mM).
C RATES(35) is d/dt K_i in component K (mM).
C RATES(46) is d/dt CaMK_trap in component CaMK_active (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 96485
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.4947322902196
CONSTS(10) = 0
CONSTS(11) = 1000
CONSTS(12) = 0.5
CONSTS(13) = -80
CONSTS(14) = 0.0826
STATES(2) = 0.476679784307081
STATES(3) = 0.255181042718489
STATES(4) = 0.0512272866698557
STATES(5) = 0.0045705849736195
STATES(6) = 0.000152923243503507
STATES(7) = 0.132814156809568
STATES(8) = 0.053324572164618
STATES(9) = 0.00713656414955778
STATES(10) = 0.000318368292970148
STATES(11) = 0.0138769647384953
STATES(12) = 0.00371438476499172
STATES(13) = 0.000248553788631871
STATES(14) = 0.000644426171899078
STATES(15) = 8.62495047666861e-5
STATES(16) = 1.12344645485893e-5
STATES(17) = 1.05292905767805e-5
STATES(18) = 8.38749860095482e-5
CONSTS(15) = 0.00015552
STATES(19) = 0.99996982093287
STATES(20) = 2.26882961412491e-5
STATES(21) = 7.49117189549873e-6
STATES(22) = 1.23556404128518e-10
STATES(23) = 5.89853003114132e-11
STATES(24) = 8.01952813127322e-17
STATES(25) = 1.9237001017021e-9
CONSTS(16) = 1
CONSTS(17) = 1e-6
CONSTS(18) = 1
CONSTS(19) = 0.341
STATES(26) = 0.000119239795531548
CONSTS(20) = 9.075
STATES(27) = 0.00101841101248412
STATES(28) = 0.991354143386887
STATES(29) = 0.99709036251243
CONSTS(21) = 2.6
CONSTS(22) = 1.5
CONSTS(23) = 1.4
CONSTS(24) = 1
CONSTS(25) = -92
STATES(30) = 8.01031392415197
CONSTS(26) = 0.000125
CONSTS(27) = 4.5
CONSTS(28) = 0.32
CONSTS(29) = 0.27
CONSTS(30) = 12.3
CONSTS(31) = 87.5
CONSTS(32) = 0.0036
CONSTS(33) = 1.3
STATES(31) = 0.00011296707102434
STATES(32) = 8.01033434859049
CONSTS(34) = 0.00276
CONSTS(35) = 0.0575
CONSTS(36) = 0.0005
CONSTS(37) = 1.99508e-7
CONSTS(38) = 1
CONSTS(39) = 0.341
CONSTS(40) = 3.2e-9
CONSTS(41) = 0.000225
CONSTS(42) = 0.0065
STATES(33) = 0.00101841101248412
STATES(34) = 0.319611809402137
CONSTS(43) = 600
CONSTS(44) = 0.01833
STATES(35) = 144.907134343314
STATES(36) = 20.6683603724432
CONSTS(45) = 0.5
CONSTS(46) = 2.46108e-5
CONSTS(47) = 1.77e-5
CONSTS(48) = 0.0138542
STATES(37) = 4.06959064407486e-6
CONSTS(49) = 0.497458
STATES(38) = 1.76776038804228e-5
STATES(39) = 0.997976503910799
STATES(40) = 0.947581184798815
CONSTS(50) = 9e-7
CONSTS(51) = -1
CONSTS(52) = 1
CONSTS(53) = 0.4
STATES(41) = 0.000714227102027745
STATES(42) = 0.996044894628276
CONSTS(54) = 8
CONSTS(55) = 1
CONSTS(56) = 4.75
CONSTS(57) = 0.1125
CONSTS(58) = 8
CONSTS(59) = 10
CONSTS(60) = 1
CONSTS(61) = 0.28
CONSTS(62) = 0.0123
STATES(43) = 0.993914988616979
CONSTS(63) = 0.00092
CONSTS(64) = 0.00017
CONSTS(65) = 0.004375
CONSTS(66) = 0.75
CONSTS(67) = 0.15
CONSTS(68) = 15
STATES(44) = 1.01475649943057
CONSTS(69) = 25
CONSTS(70) = 0.2
CONSTS(71) = 2
CONSTS(72) = 0.05
CONSTS(73) = 0.07
CONSTS(74) = 0.00238
CONSTS(75) = 0.0005
CONSTS(76) = 0.8
CONSTS(77) = 10
CONSTS(78) = 0.047
CONSTS(79) = 0.00087
CONSTS(80) = 1.124
CONSTS(81) = 0.0087
STATES(45) = 20.6683605553303
CONSTS(82) = 0.0015
STATES(46) = 0.0114911034744777
CONSTS(83) = 0.05
CONSTS(84) = 0.00068
CONSTS(85) = 0.05
CONSTS(86) = 0.0895383
CONSTS(87) = 1000.00* 3.14159265358979*CONSTS(9)*CONSTS(9)*CONSTS(8)
CONSTS(88) = CONSTS(45)* (CONSTS(4)/5.40000) ** (1.0 / 2)
CONSTS(89) = CONSTS(48)* (CONSTS(4)/5.40000) ** (1.0 / 2)
CONSTS(90) = CONSTS(87)*0.0600000
CONSTS(91) = 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(9)+ 2.00000* 3.14159265358979*CONSTS(9)*CONSTS(8)
CONSTS(92) = CONSTS(91)*2.00000
CONSTS(93) = ( 1.00000*1.00000*CONSTS(92))/( CONSTS(1)*1.00000)
CONSTS(94) = CONSTS(87)*0.678000
CONSTS(95) = CONSTS(87)*0.0552000
CONSTS(96) = CONSTS(87)*0.00480000
CONSTS(97) = CONSTS(87)*0.0200000
CONSTS(98) = CONSTS(87)*0.00200000
CONSTS(99) = CONSTS(4)/(CONSTS(4)+CONSTS(22))
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(34) = (ALGBRC(7) - STATES(34))/CONSTS(43)
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(33) = ALGBRC(6)*(1.00000 - STATES(33)) - ALGBRC(18)*STATES(33)
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(37) = (ALGBRC(8) - STATES(37))/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(39) = ALGBRC(21)*(1.00000 - STATES(39)) - ALGBRC(10)*STATES(39)
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(40) = ALGBRC(22)*(1.00000 - STATES(40)) - ALGBRC(11)*STATES(40)
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(42) = (ALGBRC(27) - STATES(42))/CONSTS(54)
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(38) = (ALGBRC(30) - STATES(38))/ALGBRC(26)
ALGBRC(32) = 0.0148646*EXP(( 0.0298773*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(35) = 0.0839863*EXP(( - 0.0554611*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(2) = - 4.00000*ALGBRC(32)*STATES(2)+ ALGBRC(35)*STATES(3)
ALGBRC(40) = 0.00311733*EXP(( - 0.426255*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(38) = 0.0146007*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(3) = - ( 3.00000*ALGBRC(32)+ALGBRC(35)+ALGBRC(38))*STATES(3)+ 4.00000*ALGBRC(32)*STATES(2)+ 2.00000*ALGBRC(35)*STATES(4)+ ALGBRC(40)*STATES(7)
RATES(4) = - ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38))*STATES(4)+ 3.00000*ALGBRC(32)*STATES(3)+ 3.00000*ALGBRC(35)*STATES(5)+ ALGBRC(40)*STATES(8)
RATES(5) = - (ALGBRC(32)+ 3.00000*ALGBRC(35)+ 3.00000*ALGBRC(38))*STATES(5)+ 2.00000*ALGBRC(32)*STATES(4)+ 4.00000*ALGBRC(35)*STATES(6)+ ALGBRC(40)*STATES(9)
RATES(6) = - ( 4.00000*ALGBRC(35)+ 4.00000*ALGBRC(38))*STATES(6)+ ALGBRC(32)*STATES(5)+ ALGBRC(40)*STATES(10)
RATES(7) = - ( 3.00000*ALGBRC(32)+ALGBRC(40))*STATES(7)+ ALGBRC(35)*STATES(8)+ ALGBRC(38)*STATES(3)
RATES(8) = - ( 2.00000*ALGBRC(32)+ALGBRC(35)+ALGBRC(38)+ALGBRC(40))*STATES(8)+ 3.00000*ALGBRC(32)*STATES(7)+ 2.00000*ALGBRC(35)*STATES(9)+ 2.00000*ALGBRC(38)*STATES(4)+ 2.00000*ALGBRC(40)*STATES(11)
RATES(9) = - (ALGBRC(32)+ 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38)+ALGBRC(40))*STATES(9)+ 2.00000*ALGBRC(32)*STATES(8)+ 3.00000*ALGBRC(35)*STATES(10)+ 3.00000*ALGBRC(38)*STATES(5)+ 2.00000*ALGBRC(40)*STATES(12)
RATES(10) = - ( 3.00000*ALGBRC(35)+ 3.00000*ALGBRC(38)+ALGBRC(40))*STATES(10)+ ALGBRC(32)*STATES(9)+ 4.00000*ALGBRC(38)*STATES(6)+ 2.00000*ALGBRC(40)*STATES(13)
RATES(11) = - ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(40))*STATES(11)+ ALGBRC(35)*STATES(12)+ ALGBRC(38)*STATES(8)
RATES(12) = - (ALGBRC(32)+ALGBRC(35)+ALGBRC(38)+ 2.00000*ALGBRC(40))*STATES(12)+ 2.00000*ALGBRC(32)*STATES(11)+ 2.00000*ALGBRC(35)*STATES(13)+ 2.00000*ALGBRC(38)*STATES(9)+ 3.00000*ALGBRC(40)*STATES(14)
RATES(13) = - ( 2.00000*ALGBRC(35)+ 2.00000*ALGBRC(38)+ 2.00000*ALGBRC(40))*STATES(13)+ ALGBRC(32)*STATES(12)+ 3.00000*ALGBRC(38)*STATES(10)+ 3.00000*ALGBRC(40)*STATES(15)
RATES(14) = - (ALGBRC(32)+ 3.00000*ALGBRC(40))*STATES(14)+ ALGBRC(35)*STATES(15)+ ALGBRC(38)*STATES(12)
RATES(15) = - (ALGBRC(35)+ALGBRC(38)+ 3.00000*ALGBRC(40))*STATES(15)+ ALGBRC(32)*STATES(14)+ 2.00000*ALGBRC(38)*STATES(13)+ 4.00000*ALGBRC(40)*STATES(16)
ALGBRC(43) = 0.0773199*EXP(( - 0.0647261*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(16) = - ( 4.00000*ALGBRC(40)+CONSTS(86))*STATES(16)+ ALGBRC(38)*STATES(15)+ ALGBRC(43)*STATES(17)
ALGBRC(28) = 1.00000 - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17))
ALGBRC(46) = 0.794055*EXP(( - 0.0801738*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(49) = 0.586382*EXP(( 0.282056*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
RATES(17) = - (ALGBRC(43)+ALGBRC(49))*STATES(17)+ ALGBRC(46)*ALGBRC(28)+ CONSTS(86)*STATES(16)
ALGBRC(14) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
ALGBRC(25) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
ALGBRC(29) = ALGBRC(25)/ALGBRC(14)
ALGBRC(33) = (1.00000 - ALGBRC(25))/ALGBRC(14)
ALGBRC(60) = 3.00000/(1.00000+0.00300000/STATES(26) ** 4.00000)
ALGBRC(39) = 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(36) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
ALGBRC(41) = ALGBRC(36)/ALGBRC(39)
ALGBRC(44) = (1.00000 - ALGBRC(36))/ALGBRC(39)
RATES(19) = - (ALGBRC(29)+ALGBRC(60)+ALGBRC(44))*STATES(19)+ ALGBRC(33)*STATES(20)+ CONSTS(16)*STATES(21)+ ALGBRC(41)*STATES(23)
RATES(20) = - (ALGBRC(33)+ALGBRC(60)+ALGBRC(44))*STATES(20)+ ALGBRC(29)*STATES(19)+ CONSTS(16)*STATES(22)+ ALGBRC(41)*STATES(24)
ALGBRC(50) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(26) ** 4.00000)
ALGBRC(52) = 1.00000/( ALGBRC(50)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
ALGBRC(54) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(52))
ALGBRC(47) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
ALGBRC(56) = ALGBRC(47)/ALGBRC(54)
ALGBRC(58) = (1.00000 - ALGBRC(47))/ALGBRC(54)
RATES(21) = - (ALGBRC(29)+CONSTS(16)+ALGBRC(58))*STATES(21)+ ALGBRC(60)*STATES(19)+ ALGBRC(33)*STATES(22)+ ALGBRC(56)*STATES(25)
ALGBRC(2) = 1.00000 - (STATES(19)+STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25))
RATES(22) = - (ALGBRC(33)+CONSTS(16)+ALGBRC(58))*STATES(22)+ ALGBRC(60)*STATES(20)+ ALGBRC(29)*STATES(21)+ ALGBRC(56)*ALGBRC(2)
ALGBRC(62) = ( CONSTS(17)*ALGBRC(41)*ALGBRC(58)*ALGBRC(60))/( ALGBRC(44)*ALGBRC(56)*CONSTS(16))
RATES(23) = - (ALGBRC(29)+ALGBRC(62)+ALGBRC(41))*STATES(23)+ ALGBRC(44)*STATES(19)+ CONSTS(17)*STATES(25)+ ALGBRC(33)*STATES(24)
RATES(24) = - (ALGBRC(33)+ALGBRC(62)+ALGBRC(41))*STATES(24)+ ALGBRC(44)*STATES(20)+ CONSTS(17)*ALGBRC(2)+ ALGBRC(29)*STATES(23)
RATES(25) = - (ALGBRC(29)+CONSTS(17)+ALGBRC(56))*STATES(25)+ ALGBRC(62)*STATES(23)+ ALGBRC(58)*STATES(21)+ ALGBRC(33)*ALGBRC(2)
ALGBRC(77) = (STATES(44) - STATES(43))/CONSTS(69)
ALGBRC(80) = 1.00000/(1.00000+( CONSTS(76)*CONSTS(77))/STATES(43)+CONSTS(76) ** 2.00000)
RATES(43) = ALGBRC(80)*(ALGBRC(77) - STATES(41))
ALGBRC(63) = 0.200000*CONSTS(27)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*STATES(31)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(64) = 1.00000+CONSTS(26)/STATES(31) ** 2.00000
ALGBRC(65) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(66) = CONSTS(33)*STATES(32) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(31)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(31)/CONSTS(32))
ALGBRC(67) = CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(30) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(31)
ALGBRC(68) = ALGBRC(63)/( ALGBRC(64)*ALGBRC(65)*(ALGBRC(66)+ALGBRC(67)))
ALGBRC(79) = (STATES(31) - STATES(18))/CONSTS(70)
ALGBRC(83) = (STATES(31) - STATES(26))/CONSTS(71)
ALGBRC(86) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(31) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(31) ** 2.00000)
RATES(31) = - ALGBRC(86)*((ALGBRC(79)+ALGBRC(83)) - (( 2.00000*ALGBRC(68)*CONSTS(93))/( 2.00000*CONSTS(97))+( STATES(41)*CONSTS(96))/CONSTS(97)))
ALGBRC(34) = ( (( CONSTS(15)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(18)*STATES(26)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(19)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(37) = ALGBRC(34)*(STATES(20)+STATES(22))
ALGBRC(87) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(26) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(26) ** 2.00000)
RATES(26) = - ALGBRC(87)*(( ALGBRC(37)*CONSTS(93))/( 2.00000*CONSTS(98)) - ( ALGBRC(83)*CONSTS(97))/CONSTS(98))
ALGBRC(85) = ( CONSTS(85)*(1.00000 - STATES(46)))/(1.00000+CONSTS(82)/STATES(31))+STATES(46)
RATES(46) = CONSTS(83)*ALGBRC(85)*(ALGBRC(85) - STATES(46)) - CONSTS(84)*STATES(46)
ALGBRC(72) = ( CONSTS(35)*STATES(18))/(CONSTS(36)+STATES(18))
ALGBRC(73) = ( (( CONSTS(37)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(38)*STATES(18)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(39)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(53) = 0.800000*CONSTS(27)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*STATES(18)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(55) = 1.00000+CONSTS(26)/STATES(18) ** 2.00000
ALGBRC(57) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(59) = CONSTS(33)*STATES(30) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(18)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(18)/CONSTS(32))
ALGBRC(61) = CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(30) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(18)
ALGBRC(69) = ALGBRC(53)/( ALGBRC(55)*ALGBRC(57)*(ALGBRC(59)+ALGBRC(61)))
ALGBRC(70) = ALGBRC(69)+ALGBRC(68)
ALGBRC(93) = ( CONSTS(64)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
ALGBRC(90) = ( CONSTS(66)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
ALGBRC(96) = ( (ALGBRC(90)+1.00000)*CONSTS(65)*STATES(18))/((STATES(18)+CONSTS(63)) - ALGBRC(93))
ALGBRC(75) = (0.00437500/CONSTS(68))*STATES(44)
ALGBRC(82) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(74))/STATES(18)+CONSTS(74) ** 2.00000+( CONSTS(75)*CONSTS(73))/STATES(18)+CONSTS(75) ** 2.00000)
RATES(18) = ALGBRC(82)*(( - ((ALGBRC(73)+ALGBRC(72)) - 2.00000*ALGBRC(70))*CONSTS(93))/( CONSTS(94)*2.00000)+( (ALGBRC(75) - ALGBRC(96))*CONSTS(95))/CONSTS(94)+( ALGBRC(79)*CONSTS(97))/CONSTS(94))
RATES(44) = ALGBRC(96) - (ALGBRC(75)+( ALGBRC(77)*CONSTS(96))/CONSTS(95))
ALGBRC(89) = CONSTS(60)/(1.00000+CONSTS(61)/( 1.00000*ALGBRC(85)) ** CONSTS(59))
ALGBRC(92) = CONSTS(56)*(1.00000+ALGBRC(89))
ALGBRC(95) = ALGBRC(92)*CONSTS(57)
ALGBRC(100) = ( ALGBRC(37)*ALGBRC(95))/(1.00000+CONSTS(55)/STATES(43) ** CONSTS(58))
ALGBRC(98) = ALGBRC(92)/(1.00000+CONSTS(62)/STATES(43))
RATES(41) = - (ALGBRC(100)+STATES(41))/ALGBRC(98)
ALGBRC(81) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
ALGBRC(84) = CONSTS(20)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(81))
ALGBRC(42) = ( CONSTS(24)*(STATES(1) - CONSTS(25))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
ALGBRC(48) = STATES(30)/(STATES(30)+CONSTS(21)) ** 3.00000
ALGBRC(51) = CONSTS(23)*ALGBRC(45)*CONSTS(99)*ALGBRC(48)
ALGBRC(88) = CONSTS(42)*STATES(33) ** 3.00000*STATES(34)*(STATES(1) - ALGBRC(81))
ALGBRC(76) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(78) = ( CONSTS(1)*CONSTS(40)*ALGBRC(76)*( STATES(30)*EXP(ALGBRC(76)) - CONSTS(6)))/(EXP(ALGBRC(76)) - 1.00000)
ALGBRC(102) = (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(36))
ALGBRC(108) = ( CONSTS(46)*ALGBRC(81) - ALGBRC(102) ** 4.00000)/(ALGBRC(81) - ALGBRC(102) ** 4.00000+87.8251 ** 4.00000)
ALGBRC(110) = (STATES(32) - STATES(30))/CONSTS(70)
RATES(30) = - (( ( 3.00000*ALGBRC(70)+ 3.00000*ALGBRC(51)+ALGBRC(84)+ALGBRC(88)+ALGBRC(78))*CONSTS(93))/CONSTS(94) - (ALGBRC(108)+( ALGBRC(110)*CONSTS(97))/CONSTS(94)))
RATES(32) = - (( 3.00000*ALGBRC(68)*CONSTS(93))/CONSTS(97)+ALGBRC(110))
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(94) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(35))
ALGBRC(104) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(94)) - 59.2150)))
ALGBRC(105) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(94))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(94)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(94))+4.75300)))
ALGBRC(106) = ALGBRC(104)/(ALGBRC(104)+ALGBRC(105))
ALGBRC(107) = CONSTS(88)*ALGBRC(106)*(STATES(1) - ALGBRC(94))
ALGBRC(71) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
ALGBRC(97) = CONSTS(34)*ALGBRC(71)*(STATES(1) - ALGBRC(94))
ALGBRC(24) = CONSTS(14)*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(18) ** 1.40000))
ALGBRC(31) = STATES(17)+ALGBRC(28)
ALGBRC(99) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(44)*CONSTS(6))/(STATES(35)+ CONSTS(44)*STATES(30)))
ALGBRC(101) = ALGBRC(24)*ALGBRC(31)*(STATES(1) - ALGBRC(99))
ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
ALGBRC(112) = CONSTS(89)*STATES(37)*ALGBRC(111)*(STATES(1) - ALGBRC(94))
ALGBRC(113) = EXP(STATES(1)/550.000)
ALGBRC(114) = CONSTS(49)*STATES(38) ** 3.00000*STATES(39)*STATES(40)*ALGBRC(113)*(STATES(1) - ALGBRC(94))
ALGBRC(109) = ( CONSTS(47)*(ALGBRC(94) - ALGBRC(102)))/((ALGBRC(94)+87.8251) - ALGBRC(102))
RATES(35) = - (( ((ALGBRC(101)+ALGBRC(112)+ALGBRC(107)+ALGBRC(114)+ALGBRC(97)+ALGBRC(13)) - 2.00000*ALGBRC(51))*CONSTS(93))/CONSTS(94) - ALGBRC(109))
ALGBRC(74) = (ALGBRC(37)+ALGBRC(73)+ALGBRC(72)) - 2.00000*(ALGBRC(70)+ALGBRC(68))
ALGBRC(91) = ALGBRC(84)+ 3.00000*(ALGBRC(70)+ALGBRC(68))+ 3.00000*ALGBRC(51)+ALGBRC(88)+ALGBRC(78)
ALGBRC(115) = ALGBRC(112)+ALGBRC(101)+ALGBRC(107)+ALGBRC(97)+ - 2.00000*ALGBRC(51)+ALGBRC(114)+ALGBRC(13)
ALGBRC(116) = ( (( CONSTS(50)*CONSTS(51) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(36) - CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(117) = 1.00000 - 1.00000/(1.00000+STATES(41)/CONSTS(53) ** 2.00000)
ALGBRC(118) = ALGBRC(116)*CONSTS(52)*STATES(42)*ALGBRC(117)
ALGBRC(103) = CONSTS(41)*(STATES(1) - ALGBRC(102))
ALGBRC(119) = ALGBRC(103)+ALGBRC(118)
RATES(1) = - (ALGBRC(91)+ALGBRC(115)+ALGBRC(74)+ALGBRC(119))
ALGBRC(120) = (STATES(45) - STATES(36))/CONSTS(70)
RATES(36) = - (( ALGBRC(103)*CONSTS(93))/( - 1.00000*CONSTS(94)) - (ALGBRC(108)+ALGBRC(109)+( ALGBRC(120)*CONSTS(97))/CONSTS(94)))
RATES(45) = - (( ALGBRC(118)*CONSTS(93))/( - 1.00000*CONSTS(97))+ALGBRC(120))
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(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(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.0148646*EXP(( 0.0298773*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(35) = 0.0839863*EXP(( - 0.0554611*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(40) = 0.00311733*EXP(( - 0.426255*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(38) = 0.0146007*EXP(( 0.244650*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(43) = 0.0773199*EXP(( - 0.0647261*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(28) = 1.00000 - (STATES(2)+STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)+STATES(10)+STATES(11)+STATES(12)+STATES(13)+STATES(14)+STATES(15)+STATES(16)+STATES(17))
ALGBRC(46) = 0.794055*EXP(( - 0.0801738*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(49) = 0.586382*EXP(( 0.282056*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(14) = 0.590000+( 0.800000*EXP( 0.0520000*(STATES(1)+13.0000)))/(1.00000+EXP( 0.132000*(STATES(1)+13.0000)))
ALGBRC(25) = 1.00000/(1.00000+EXP(- (STATES(1) - 13.5600)/9.45000))
ALGBRC(29) = ALGBRC(25)/ALGBRC(14)
ALGBRC(33) = (1.00000 - ALGBRC(25))/ALGBRC(14)
ALGBRC(60) = 3.00000/(1.00000+0.00300000/STATES(26) ** 4.00000)
ALGBRC(39) = 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(36) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.250000)/1.25000
ALGBRC(41) = ALGBRC(36)/ALGBRC(39)
ALGBRC(44) = (1.00000 - ALGBRC(36))/ALGBRC(39)
ALGBRC(50) = 25.0000 - 17.5000/(1.00000+0.00300000/STATES(26) ** 4.00000)
ALGBRC(52) = 1.00000/( ALGBRC(50)*(1.00000+EXP(- (STATES(1)+0.213000)/10.8070)))
ALGBRC(54) = 1.00000/(1.00000/( 24.8280*(1.00000+EXP((STATES(1)+49.1000)/10.3490)))+ALGBRC(52))
ALGBRC(47) = (1.00000/(1.00000+EXP((STATES(1)+17.5000)/3.00000))+0.000100000)/1.00010
ALGBRC(56) = ALGBRC(47)/ALGBRC(54)
ALGBRC(58) = (1.00000 - ALGBRC(47))/ALGBRC(54)
ALGBRC(2) = 1.00000 - (STATES(19)+STATES(20)+STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25))
ALGBRC(62) = ( CONSTS(17)*ALGBRC(41)*ALGBRC(58)*ALGBRC(60))/( ALGBRC(44)*ALGBRC(56)*CONSTS(16))
ALGBRC(77) = (STATES(44) - STATES(43))/CONSTS(69)
ALGBRC(80) = 1.00000/(1.00000+( CONSTS(76)*CONSTS(77))/STATES(43)+CONSTS(76) ** 2.00000)
ALGBRC(63) = 0.200000*CONSTS(27)*( STATES(32) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*STATES(31)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(64) = 1.00000+CONSTS(26)/STATES(31) ** 2.00000
ALGBRC(65) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(66) = CONSTS(33)*STATES(32) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(31)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(31)/CONSTS(32))
ALGBRC(67) = CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(32)/CONSTS(30) ** 3.00000)+ STATES(32) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(31)
ALGBRC(68) = ALGBRC(63)/( ALGBRC(64)*ALGBRC(65)*(ALGBRC(66)+ALGBRC(67)))
ALGBRC(79) = (STATES(31) - STATES(18))/CONSTS(70)
ALGBRC(83) = (STATES(31) - STATES(26))/CONSTS(71)
ALGBRC(86) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(31) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(31) ** 2.00000)
ALGBRC(34) = ( (( CONSTS(15)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(18)*STATES(26)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(19)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(37) = ALGBRC(34)*(STATES(20)+STATES(22))
ALGBRC(87) = 1.00000/(1.00000+( CONSTS(78)*CONSTS(79))/CONSTS(79)+STATES(26) ** 2.00000+( CONSTS(80)*CONSTS(81))/CONSTS(81)+STATES(26) ** 2.00000)
ALGBRC(85) = ( CONSTS(85)*(1.00000 - STATES(46)))/(1.00000+CONSTS(82)/STATES(31))+STATES(46)
ALGBRC(72) = ( CONSTS(35)*STATES(18))/(CONSTS(36)+STATES(18))
ALGBRC(73) = ( (( CONSTS(37)*4.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*( CONSTS(38)*STATES(18)*EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(39)*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - 1.00000)
ALGBRC(53) = 0.800000*CONSTS(27)*( STATES(30) ** 3.00000*CONSTS(5)*EXP(( CONSTS(29)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))) - CONSTS(6) ** 3.00000*STATES(18)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(55) = 1.00000+CONSTS(26)/STATES(18) ** 2.00000
ALGBRC(57) = 1.00000+ CONSTS(28)*EXP(( (CONSTS(29) - 1.00000)*STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))
ALGBRC(59) = CONSTS(33)*STATES(30) ** 3.00000+ CONSTS(31) ** 3.00000*STATES(18)+ CONSTS(30) ** 3.00000*CONSTS(5)*(1.00000+STATES(18)/CONSTS(32))
ALGBRC(61) = CONSTS(32)*CONSTS(6) ** 3.00000*(1.00000+STATES(30)/CONSTS(30) ** 3.00000)+ STATES(30) ** 3.00000*CONSTS(5)+ CONSTS(6) ** 3.00000*STATES(18)
ALGBRC(69) = ALGBRC(53)/( ALGBRC(55)*ALGBRC(57)*(ALGBRC(59)+ALGBRC(61)))
ALGBRC(70) = ALGBRC(69)+ALGBRC(68)
ALGBRC(93) = ( CONSTS(64)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
ALGBRC(90) = ( CONSTS(66)*ALGBRC(85)*1.00000)/(CONSTS(67)+ ALGBRC(85)*1.00000)
ALGBRC(96) = ( (ALGBRC(90)+1.00000)*CONSTS(65)*STATES(18))/((STATES(18)+CONSTS(63)) - ALGBRC(93))
ALGBRC(75) = (0.00437500/CONSTS(68))*STATES(44)
ALGBRC(82) = 1.00000/(1.00000+( CONSTS(72)*CONSTS(74))/STATES(18)+CONSTS(74) ** 2.00000+( CONSTS(75)*CONSTS(73))/STATES(18)+CONSTS(75) ** 2.00000)
ALGBRC(89) = CONSTS(60)/(1.00000+CONSTS(61)/( 1.00000*ALGBRC(85)) ** CONSTS(59))
ALGBRC(92) = CONSTS(56)*(1.00000+ALGBRC(89))
ALGBRC(95) = ALGBRC(92)*CONSTS(57)
ALGBRC(100) = ( ALGBRC(37)*ALGBRC(95))/(1.00000+CONSTS(55)/STATES(43) ** CONSTS(58))
ALGBRC(98) = ALGBRC(92)/(1.00000+CONSTS(62)/STATES(43))
ALGBRC(81) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(6)/STATES(30))
ALGBRC(84) = CONSTS(20)*STATES(27) ** 3.00000*STATES(28)*STATES(29)*(STATES(1) - ALGBRC(81))
ALGBRC(42) = ( CONSTS(24)*(STATES(1) - CONSTS(25))*CONSTS(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(45) = 1.00000/(1.00000+EXP(- ALGBRC(42)))
ALGBRC(48) = STATES(30)/(STATES(30)+CONSTS(21)) ** 3.00000
ALGBRC(51) = CONSTS(23)*ALGBRC(45)*CONSTS(99)*ALGBRC(48)
ALGBRC(88) = CONSTS(42)*STATES(33) ** 3.00000*STATES(34)*(STATES(1) - ALGBRC(81))
ALGBRC(76) = ( CONSTS(1)*STATES(1))/( CONSTS(3)*CONSTS(2))
ALGBRC(78) = ( CONSTS(1)*CONSTS(40)*ALGBRC(76)*( STATES(30)*EXP(ALGBRC(76)) - CONSTS(6)))/(EXP(ALGBRC(76)) - 1.00000)
ALGBRC(102) = (( - CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(7)/STATES(36))
ALGBRC(108) = ( CONSTS(46)*ALGBRC(81) - ALGBRC(102) ** 4.00000)/(ALGBRC(81) - ALGBRC(102) ** 4.00000+87.8251 ** 4.00000)
ALGBRC(110) = (STATES(32) - STATES(30))/CONSTS(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(94) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log(CONSTS(4)/STATES(35))
ALGBRC(104) = 1.02000/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(94)) - 59.2150)))
ALGBRC(105) = ( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(94))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(94)) - 594.310)))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(94))+4.75300)))
ALGBRC(106) = ALGBRC(104)/(ALGBRC(104)+ALGBRC(105))
ALGBRC(107) = CONSTS(88)*ALGBRC(106)*(STATES(1) - ALGBRC(94))
ALGBRC(71) = 1.00000/(1.00000+EXP((7.48800 - STATES(1))/5.98000))
ALGBRC(97) = CONSTS(34)*ALGBRC(71)*(STATES(1) - ALGBRC(94))
ALGBRC(24) = CONSTS(14)*(1.00000+0.600000/(1.00000+3.80000e-05/STATES(18) ** 1.40000))
ALGBRC(31) = STATES(17)+ALGBRC(28)
ALGBRC(99) = (( CONSTS(3)*CONSTS(2))/CONSTS(1))*log((CONSTS(4)+ CONSTS(44)*CONSTS(6))/(STATES(35)+ CONSTS(44)*STATES(30)))
ALGBRC(101) = ALGBRC(24)*ALGBRC(31)*(STATES(1) - ALGBRC(99))
ALGBRC(111) = 1.00000/(1.00000+EXP((STATES(1)+10.0000)/15.4000))
ALGBRC(112) = CONSTS(89)*STATES(37)*ALGBRC(111)*(STATES(1) - ALGBRC(94))
ALGBRC(113) = EXP(STATES(1)/550.000)
ALGBRC(114) = CONSTS(49)*STATES(38) ** 3.00000*STATES(39)*STATES(40)*ALGBRC(113)*(STATES(1) - ALGBRC(94))
ALGBRC(109) = ( CONSTS(47)*(ALGBRC(94) - ALGBRC(102)))/((ALGBRC(94)+87.8251) - ALGBRC(102))
ALGBRC(74) = (ALGBRC(37)+ALGBRC(73)+ALGBRC(72)) - 2.00000*(ALGBRC(70)+ALGBRC(68))
ALGBRC(91) = ALGBRC(84)+ 3.00000*(ALGBRC(70)+ALGBRC(68))+ 3.00000*ALGBRC(51)+ALGBRC(88)+ALGBRC(78)
ALGBRC(115) = ALGBRC(112)+ALGBRC(101)+ALGBRC(107)+ALGBRC(97)+ - 2.00000*ALGBRC(51)+ALGBRC(114)+ALGBRC(13)
ALGBRC(116) = ( (( CONSTS(50)*CONSTS(51) ** 2.00000*STATES(1)*CONSTS(1) ** 2.00000)/( CONSTS(3)*CONSTS(2)))*(STATES(36) - CONSTS(7)*EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2)))))/(1.00000 - EXP(( STATES(1)*CONSTS(1))/( CONSTS(3)*CONSTS(2))))
ALGBRC(117) = 1.00000 - 1.00000/(1.00000+STATES(41)/CONSTS(53) ** 2.00000)
ALGBRC(118) = ALGBRC(116)*CONSTS(52)*STATES(42)*ALGBRC(117)
ALGBRC(103) = CONSTS(41)*(STATES(1) - ALGBRC(102))
ALGBRC(119) = ALGBRC(103)+ALGBRC(118)
ALGBRC(120) = (STATES(45) - STATES(36))/CONSTS(70)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END