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 65 entries in the algebraic variable array.
C There are a total of 34 entries in each of the rate and state variable arrays.
C There are a total of 155 entries in the constant variable array.
C
C
C CONSTS(1) is t_ss in component Clamp_parameters (ms).
C CONSTS(2) is V_actHolding in component Clamp_parameters (mV).
C CONSTS(3) is V_actTest in component Clamp_parameters (mV).
C CONSTS(4) is Cai in component Clamp_parameters (mM).
C CONSTS(5) is Ki in component Clamp_parameters (mM).
C CONSTS(6) is Nai in component Clamp_parameters (mM).
C VOI is time in component time (ms).
C ALGBRC(1) is V in component Patch_clamp_protocol (mV).
C ALGBRC(51) is I_ion in component Clamped_current_Xi (pA).
C ALGBRC(23) is I in component Ionic_currents (pA).
C ALGBRC(27) is I in component Ionic_currents (pA).
C ALGBRC(31) is I in component Ionic_currents (pA).
C ALGBRC(35) is I in component Ionic_currents (pA).
C ALGBRC(39) is I in component Ionic_currents (pA).
C ALGBRC(15) is I_NCX in component INCX (pA).
C ALGBRC(19) is I_NaK in component INaK (pA).
C ALGBRC(43) is I in component Ionic_currents (pA).
C ALGBRC(47) is I in component Ionic_currents (pA).
C STATES(1) is O_CaL in component ICaL_channel_states (dimensionless).
C ALGBRC(2) is P_CaT in component ICaT (dimensionless).
C STATES(2) is O4 in component IBK_channel_states (dimensionless).
C ALGBRC(10) is P_Kv in component IKv (dimensionless).
C STATES(3) is O_Na in component INa_channel_states (dimensionless).
C STATES(4) is X in component gating_kinetics (dimensionless).
C STATES(5) is X in component gating_kinetics (dimensionless).
C STATES(6) is X in component gating_kinetics (dimensionless).
C STATES(7) is X in component gating_kinetics (dimensionless).
C CONSTS(7) is g_CaL in component model_parameters (nS).
C CONSTS(8) is R in component model_parameters (J_per_K_mol).
C CONSTS(9) is T in component model_parameters (kelvin).
C CONSTS(10) is F in component model_parameters (C_per_mmol).
C CONSTS(11) is T0CaL in component model_parameters (kelvin).
C CONSTS(12) is T0CaT in component model_parameters (kelvin).
C CONSTS(13) is Q10Ca in component model_parameters (dimensionless).
C CONSTS(14) is Cao in component model_parameters (mM).
C CONSTS(15) is C0_init in component initial_conditions (dimensionless).
C CONSTS(16) is C1_init in component initial_conditions (dimensionless).
C CONSTS(17) is C2_init in component initial_conditions (dimensionless).
C CONSTS(18) is C3_init in component initial_conditions (dimensionless).
C CONSTS(19) is C0Ca_init in component initial_conditions (dimensionless).
C CONSTS(20) is C1Ca_init in component initial_conditions (dimensionless).
C CONSTS(21) is C2Ca_init in component initial_conditions (dimensionless).
C CONSTS(22) is C3Ca_init in component initial_conditions (dimensionless).
C CONSTS(23) is ICa_init in component initial_conditions (dimensionless).
C CONSTS(24) is IVs_init in component initial_conditions (dimensionless).
C CONSTS(25) is IVf_init in component initial_conditions (dimensionless).
C CONSTS(26) is IVfCa_init in component initial_conditions (dimensionless).
C CONSTS(27) is IVsCa_init in component initial_conditions (dimensionless).
C CONSTS(28) is O_CaL_init in component initial_conditions (dimensionless).
C CONSTS(29) is g_CaT in component model_parameters (nS).
C CONSTS(30) is d_CaT_init in component initial_conditions (dimensionless).
C CONSTS(31) is f_CaT_init in component initial_conditions (dimensionless).
C CONSTS(32) is tau_dCaT in component model_parameters (ms).
C CONSTS(33) is g_Kv in component model_parameters (nS).
C CONSTS(34) is T0K in component model_parameters (kelvin).
C CONSTS(35) is Q10K in component model_parameters (dimensionless).
C CONSTS(36) is Ko in component model_parameters (mM).
C CONSTS(37) is x_Kv_init in component initial_conditions (dimensionless).
C CONSTS(38) is y_Kv_init in component initial_conditions (dimensionless).
C CONSTS(39) is tau_xKv in component model_parameters (ms).
C CONSTS(40) is tau_yKv in component model_parameters (ms).
C CONSTS(41) is g_BK in component model_parameters (nS).
C CONSTS(42) is BK_C0_init in component initial_conditions (dimensionless).
C CONSTS(43) is BK_C1_init in component initial_conditions (dimensionless).
C CONSTS(44) is BK_C2_init in component initial_conditions (dimensionless).
C CONSTS(45) is BK_C3_init in component initial_conditions (dimensionless).
C CONSTS(46) is BK_C4_init in component initial_conditions (dimensionless).
C CONSTS(47) is O0_init in component initial_conditions (dimensionless).
C CONSTS(48) is O1_init in component initial_conditions (dimensionless).
C CONSTS(49) is O2_init in component initial_conditions (dimensionless).
C CONSTS(50) is O3_init in component initial_conditions (dimensionless).
C CONSTS(51) is O4_init in component initial_conditions (dimensionless).
C CONSTS(52) is k_on in component model_parameters (per_ms).
C CONSTS(53) is k_off_C in component model_parameters (per_ms).
C CONSTS(54) is k_off_O in component model_parameters (per_ms).
C CONSTS(55) is g_Na in component model_parameters (nS).
C CONSTS(56) is T0Na in component model_parameters (kelvin).
C CONSTS(57) is Q10Na in component model_parameters (dimensionless).
C CONSTS(58) is Nao in component model_parameters (mM).
C CONSTS(59) is Na_C1_init in component initial_conditions (dimensionless).
C CONSTS(60) is Na_C2_init in component initial_conditions (dimensionless).
C CONSTS(61) is Na_C3_init in component initial_conditions (dimensionless).
C CONSTS(62) is I1_init in component initial_conditions (dimensionless).
C CONSTS(63) is I2_init in component initial_conditions (dimensionless).
C CONSTS(64) is O_Na_init in component initial_conditions (dimensionless).
C CONSTS(65) is K_mCa in component model_parameters (mM).
C CONSTS(66) is K_mNai in component model_parameters (mM).
C CONSTS(67) is k_sat in component model_parameters (dimensionless).
C CONSTS(68) is gamma in component model_parameters (dimensionless).
C CONSTS(69) is P_NCX in component model_parameters (pA).
C CONSTS(70) is K_mK in component model_parameters (mM).
C CONSTS(71) is K_mNa in component model_parameters (mM).
C CONSTS(72) is P_NaK in component model_parameters (pA).
C CONSTS(73) is g_NsNa in component model_parameters (nS).
C CONSTS(74) is g_NsK in component model_parameters (nS).
C CONSTS(75) is Cm in component model_parameters (pF).
C CONSTS(76) is V_cell in component model_parameters (litre).
C CONSTS(77) is n_CRT in component model_parameters (dimensionless).
C CONSTS(78) is n_CaM in component model_parameters (dimensionless).
C CONSTS(79) is CRT_total in component model_parameters (mM).
C CONSTS(80) is CaM_total in component model_parameters (mM).
C CONSTS(81) is K_D_CRT in component model_parameters (mM).
C CONSTS(82) is K_D_CaM in component model_parameters (mM4).
C CONSTS(83) is V_init in component initial_conditions (mV).
C CONSTS(84) is Cai_total_init in component initial_conditions (mM).
C CONSTS(85) is Cai_init in component initial_conditions (mM).
C CONSTS(86) is Ki_init in component initial_conditions (mM).
C CONSTS(87) is Nai_init in component initial_conditions (mM).
C ALGBRC(55) is norm in component IBK_channel_states (dimensionless).
C ALGBRC(56) is norm in component INa_channel_states (dimensionless).
C CONSTS(88) is PO in component INS_Na (dimensionless).
C CONSTS(89) is PO in component INS_K (dimensionless).
C CONSTS(126) is E in component Nernst_potential (mV).
C STATES(8) is C0 in component ICaL_channel_states (dimensionless).
C STATES(9) is C1 in component ICaL_channel_states (dimensionless).
C STATES(10) is C2 in component ICaL_channel_states (dimensionless).
C STATES(11) is C3 in component ICaL_channel_states (dimensionless).
C STATES(12) is C0Ca in component ICaL_channel_states (dimensionless).
C STATES(13) is C1Ca in component ICaL_channel_states (dimensionless).
C STATES(14) is C2Ca in component ICaL_channel_states (dimensionless).
C STATES(15) is C3Ca in component ICaL_channel_states (dimensionless).
C STATES(16) is ICa in component ICaL_channel_states (dimensionless).
C STATES(17) is IVs in component ICaL_channel_states (dimensionless).
C STATES(18) is IVf in component ICaL_channel_states (dimensionless).
C STATES(19) is IVfCa in component ICaL_channel_states (dimensionless).
C STATES(20) is IVsCa in component ICaL_channel_states (dimensionless).
C ALGBRC(3) is alpha in component ICaL_channel_states (per_ms).
C ALGBRC(11) is beta in component ICaL_channel_states (per_ms).
C ALGBRC(16) is alpha_0 in component ICaL_channel_states (per_ms).
C ALGBRC(20) is alpha_1 in component ICaL_channel_states (per_ms).
C ALGBRC(24) is alpha_2 in component ICaL_channel_states (per_ms).
C ALGBRC(28) is alpha_3 in component ICaL_channel_states (per_ms).
C ALGBRC(33) is beta_0 in component ICaL_channel_states (per_ms).
C ALGBRC(38) is beta_1 in component ICaL_channel_states (per_ms).
C ALGBRC(42) is beta_2 in component ICaL_channel_states (per_ms).
C ALGBRC(46) is beta_3 in component ICaL_channel_states (per_ms).
C CONSTS(127) is delta in component ICaL_channel_states (per_ms).
C CONSTS(128) is theta in component ICaL_channel_states (per_ms).
C ALGBRC(54) is phi_s in component ICaL_channel_states (per_ms).
C ALGBRC(50) is phi_f in component ICaL_channel_states (per_ms).
C ALGBRC(62) is omega_s in component ICaL_channel_states (per_ms).
C ALGBRC(61) is omega_f in component ICaL_channel_states (per_ms).
C ALGBRC(64) is omega_fs in component ICaL_channel_states (per_ms).
C ALGBRC(63) is omega_sf in component ICaL_channel_states (per_ms).
C ALGBRC(60) is psi_s in component ICaL_channel_states (per_ms).
C ALGBRC(59) is psi_f in component ICaL_channel_states (per_ms).
C ALGBRC(58) is xi_s in component ICaL_channel_states (per_ms).
C ALGBRC(57) is xi_f in component ICaL_channel_states (per_ms).
C ALGBRC(65) is norm in component ICaL_channel_states (dimensionless).
C CONSTS(125) is phi in component temperature_factor (dimensionless).
C CONSTS(151) is E in component Nernst_potential (mV).
C ALGBRC(4) is X_inf in component d_inf (dimensionless).
C CONSTS(129) is phi in component temperature_factor (dimensionless).
C ALGBRC(5) is X_inf in component f_inf (dimensionless).
C ALGBRC(12) is tau in component tau_f (ms).
C CONSTS(130) is phi in component temperature_factor (dimensionless).
C CONSTS(152) is E in component Nernst_potential (mV).
C ALGBRC(6) is X_inf in component x_inf (dimensionless).
C CONSTS(131) is phi in component temperature_factor (dimensionless).
C ALGBRC(7) is X_inf in component y_inf (dimensionless).
C CONSTS(132) is phi in component temperature_factor (dimensionless).
C CONSTS(133) is E in component Nernst_potential (mV).
C STATES(21) is C0 in component IBK_channel_states (dimensionless).
C STATES(22) is C1 in component IBK_channel_states (dimensionless).
C STATES(23) is C2 in component IBK_channel_states (dimensionless).
C STATES(24) is C3 in component IBK_channel_states (dimensionless).
C STATES(25) is C4 in component IBK_channel_states (dimensionless).
C STATES(26) is O0 in component IBK_channel_states (dimensionless).
C STATES(27) is O1 in component IBK_channel_states (dimensionless).
C STATES(28) is O2 in component IBK_channel_states (dimensionless).
C STATES(29) is O3 in component IBK_channel_states (dimensionless).
C ALGBRC(8) is alpha in component IBK_channel_states (per_ms).
C ALGBRC(13) is beta in component IBK_channel_states (per_ms).
C ALGBRC(17) is k_C0O0 in component IBK_channel_states (per_ms).
C ALGBRC(21) is k_C1O1 in component IBK_channel_states (per_ms).
C ALGBRC(25) is k_C2O2 in component IBK_channel_states (per_ms).
C ALGBRC(29) is k_C3O3 in component IBK_channel_states (per_ms).
C ALGBRC(34) is k_C4O4 in component IBK_channel_states (per_ms).
C ALGBRC(36) is k_O0C0 in component IBK_channel_states (per_ms).
C ALGBRC(40) is k_O1C1 in component IBK_channel_states (per_ms).
C ALGBRC(44) is k_O2C2 in component IBK_channel_states (per_ms).
C ALGBRC(48) is k_O3C3 in component IBK_channel_states (per_ms).
C ALGBRC(52) is k_O4C4 in component IBK_channel_states (per_ms).
C CONSTS(134) is k_C0C1 in component IBK_channel_states (per_ms).
C CONSTS(135) is k_C1C2 in component IBK_channel_states (per_ms).
C CONSTS(136) is k_C2C3 in component IBK_channel_states (per_ms).
C CONSTS(137) is k_C3C4 in component IBK_channel_states (per_ms).
C CONSTS(138) is k_C4C3 in component IBK_channel_states (per_ms).
C CONSTS(139) is k_C3C2 in component IBK_channel_states (per_ms).
C CONSTS(140) is k_C2C1 in component IBK_channel_states (per_ms).
C CONSTS(141) is k_C1C0 in component IBK_channel_states (per_ms).
C CONSTS(142) is k_O0O1 in component IBK_channel_states (per_ms).
C CONSTS(143) is k_O1O2 in component IBK_channel_states (per_ms).
C CONSTS(144) is k_O2O3 in component IBK_channel_states (per_ms).
C CONSTS(145) is k_O3O4 in component IBK_channel_states (per_ms).
C CONSTS(146) is k_O4O3 in component IBK_channel_states (per_ms).
C CONSTS(147) is k_O3O2 in component IBK_channel_states (per_ms).
C CONSTS(148) is k_O2O1 in component IBK_channel_states (per_ms).
C CONSTS(149) is k_O1O0 in component IBK_channel_states (per_ms).
C CONSTS(150) is E in component Nernst_potential (mV).
C STATES(30) is C1 in component INa_channel_states (dimensionless).
C STATES(31) is C2 in component INa_channel_states (dimensionless).
C STATES(32) is C3 in component INa_channel_states (dimensionless).
C STATES(33) is I1 in component INa_channel_states (dimensionless).
C STATES(34) is I2 in component INa_channel_states (dimensionless).
C ALGBRC(9) is k_I2I1 in component INa_channel_states (per_ms).
C ALGBRC(14) is k_I1O in component INa_channel_states (per_ms).
C ALGBRC(18) is k_OC1 in component INa_channel_states (per_ms).
C ALGBRC(22) is k_C1C2 in component INa_channel_states (per_ms).
C ALGBRC(26) is k_C2C3 in component INa_channel_states (per_ms).
C ALGBRC(30) is k_C3C2 in component INa_channel_states (per_ms).
C ALGBRC(32) is k_C2C1 in component INa_channel_states (per_ms).
C ALGBRC(37) is k_C1O in component INa_channel_states (per_ms).
C ALGBRC(41) is k_OI1 in component INa_channel_states (per_ms).
C ALGBRC(45) is k_I1I2 in component INa_channel_states (per_ms).
C ALGBRC(49) is k_I1C1 in component INa_channel_states (per_ms).
C ALGBRC(53) is k_C1I1 in component INa_channel_states (per_ms).
C CONSTS(153) is phi in component temperature_factor (dimensionless).
C CONSTS(154) is E in component Nernst_potential (mV).
C CONSTS(155) is E in component Nernst_potential (mV).
C CONSTS(90) is R in component constants (J_per_K_mol).
C CONSTS(91) is F in component constants (C_per_mmol).
C CONSTS(92) is z_Ca in component constants (dimensionless).
C CONSTS(93) is z_Na in component constants (dimensionless).
C CONSTS(94) is z_K in component constants (dimensionless).
C CONSTS(95) is R in component constants (J_per_K_mol).
C CONSTS(96) is F in component constants (C_per_mmol).
C CONSTS(97) is z_Ca in component constants (dimensionless).
C CONSTS(98) is z_Na in component constants (dimensionless).
C CONSTS(99) is z_K in component constants (dimensionless).
C CONSTS(100) is R in component constants (J_per_K_mol).
C CONSTS(101) is F in component constants (C_per_mmol).
C CONSTS(102) is z_K in component constants (dimensionless).
C CONSTS(103) is z_Na in component constants (dimensionless).
C CONSTS(104) is z_Ca in component constants (dimensionless).
C CONSTS(105) is R in component constants (J_per_K_mol).
C CONSTS(106) is F in component constants (C_per_mmol).
C CONSTS(107) is z_K in component constants (dimensionless).
C CONSTS(108) is z_Na in component constants (dimensionless).
C CONSTS(109) is z_Ca in component constants (dimensionless).
C CONSTS(110) is R in component constants (J_per_K_mol).
C CONSTS(111) is F in component constants (C_per_mmol).
C CONSTS(112) is z_Na in component constants (dimensionless).
C CONSTS(113) is z_Ca in component constants (dimensionless).
C CONSTS(114) is z_K in component constants (dimensionless).
C CONSTS(115) is R in component constants (J_per_K_mol).
C CONSTS(116) is F in component constants (C_per_mmol).
C CONSTS(117) is z_Na in component constants (dimensionless).
C CONSTS(118) is z_Ca in component constants (dimensionless).
C CONSTS(119) is z_K in component constants (dimensionless).
C CONSTS(120) is R in component constants (J_per_K_mol).
C CONSTS(121) is F in component constants (C_per_mmol).
C CONSTS(122) is z_K in component constants (dimensionless).
C CONSTS(123) is z_Na in component constants (dimensionless).
C CONSTS(124) is z_Ca in component constants (dimensionless).
C RATES(11) is d/dt C3 in component ICaL_channel_states (dimensionless).
C RATES(10) is d/dt C2 in component ICaL_channel_states (dimensionless).
C RATES(9) is d/dt C1 in component ICaL_channel_states (dimensionless).
C RATES(8) is d/dt C0 in component ICaL_channel_states (dimensionless).
C RATES(15) is d/dt C3Ca in component ICaL_channel_states (dimensionless).
C RATES(14) is d/dt C2Ca in component ICaL_channel_states (dimensionless).
C RATES(13) is d/dt C1Ca in component ICaL_channel_states (dimensionless).
C RATES(12) is d/dt C0Ca in component ICaL_channel_states (dimensionless).
C RATES(1) is d/dt O_CaL in component ICaL_channel_states (dimensionless).
C RATES(16) is d/dt ICa in component ICaL_channel_states (dimensionless).
C RATES(17) is d/dt IVs in component ICaL_channel_states (dimensionless).
C RATES(18) is d/dt IVf in component ICaL_channel_states (dimensionless).
C RATES(20) is d/dt IVsCa in component ICaL_channel_states (dimensionless).
C RATES(19) is d/dt IVfCa in component ICaL_channel_states (dimensionless).
C RATES(4) is d/dt X in component gating_kinetics (dimensionless).
C RATES(5) is d/dt X in component gating_kinetics (dimensionless).
C RATES(6) is d/dt X in component gating_kinetics (dimensionless).
C RATES(7) is d/dt X in component gating_kinetics (dimensionless).
C RATES(25) is d/dt C4 in component IBK_channel_states (dimensionless).
C RATES(24) is d/dt C3 in component IBK_channel_states (dimensionless).
C RATES(23) is d/dt C2 in component IBK_channel_states (dimensionless).
C RATES(22) is d/dt C1 in component IBK_channel_states (dimensionless).
C RATES(21) is d/dt C0 in component IBK_channel_states (dimensionless).
C RATES(2) is d/dt O4 in component IBK_channel_states (dimensionless).
C RATES(29) is d/dt O3 in component IBK_channel_states (dimensionless).
C RATES(28) is d/dt O2 in component IBK_channel_states (dimensionless).
C RATES(27) is d/dt O1 in component IBK_channel_states (dimensionless).
C RATES(26) is d/dt O0 in component IBK_channel_states (dimensionless).
C RATES(32) is d/dt C3 in component INa_channel_states (dimensionless).
C RATES(31) is d/dt C2 in component INa_channel_states (dimensionless).
C RATES(30) is d/dt C1 in component INa_channel_states (dimensionless).
C RATES(3) is d/dt O_Na in component INa_channel_states (dimensionless).
C RATES(34) is d/dt I2 in component INa_channel_states (dimensionless).
C RATES(33) is d/dt I1 in component INa_channel_states (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 1000
CONSTS(2) = -100
CONSTS(3) = -90
CONSTS(4) = 5.38843941249284e-5
CONSTS(5) = 153.604280337996
CONSTS(6) = 10.5731241425458
CONSTS(7) = 1.44
CONSTS(8) = 8.314
CONSTS(9) = 310
CONSTS(10) = 96.48534
CONSTS(11) = 310
CONSTS(12) = 297
CONSTS(13) = 2.1
CONSTS(14) = 2
CONSTS(15) = 0.815464741971086
CONSTS(16) = 0.152399266235657
CONSTS(17) = 0.0106805060777161
CONSTS(18) = 0.000332673548872087
CONSTS(19) = 0.0175888495282545
CONSTS(20) = 0.00328711668724504
CONSTS(21) = 0.000230369020877669
CONSTS(22) = 7.1754726923539e-6
CONSTS(23) = 8.38123983500905e-8
CONSTS(24) = 1.1193313274705e-6
CONSTS(25) = 4.0998751301597e-6
CONSTS(26) = 8.84306615061238e-8
CONSTS(27) = 2.41429816075123e-8
CONSTS(28) = 3.88576045134351e-6
CONSTS(29) = 0.0425
CONSTS(30) = 0.0791635737410974
CONSTS(31) = 0.377831534375835
CONSTS(32) = 1.9508
CONSTS(33) = 1.0217
CONSTS(34) = 297
CONSTS(35) = 3.1
CONSTS(36) = 5.4
CONSTS(37) = 0.14714161078933
CONSTS(38) = 0.99994773314105
CONSTS(39) = 4.7803
CONSTS(40) = 763.7564
CONSTS(41) = 80
CONSTS(42) = 0.48379087935899
CONSTS(43) = 0.385183559520031
CONSTS(44) = 0.115002824567753
CONSTS(45) = 0.0152602714149774
CONSTS(46) = 0.000759264410974374
CONSTS(47) = 6.94960798375172e-7
CONSTS(48) = 5.55636826398253e-8
CONSTS(49) = 2.85143702125325e-8
CONSTS(50) = 1.59832806123435e-6
CONSTS(51) = 1.82113764497095e-6
CONSTS(52) = 40633
CONSTS(53) = 11
CONSTS(54) = 1.1
CONSTS(55) = 25.1
CONSTS(56) = 297
CONSTS(57) = 2.45
CONSTS(58) = 140
CONSTS(59) = 0.0119443135223679
CONSTS(60) = 0.0109545368437155
CONSTS(61) = 0.973782548650071
CONSTS(62) = 0.000126882921013389
CONSTS(63) = 0.00318975045717667
CONSTS(64) = 1.96760342050475e-6
CONSTS(65) = 1.38
CONSTS(66) = 87.5
CONSTS(67) = 0.1
CONSTS(68) = 0.35
CONSTS(69) = 1992.1865
CONSTS(70) = 1
CONSTS(71) = 40
CONSTS(72) = 16.197
CONSTS(73) = 0.022488
CONSTS(74) = 0.017512
CONSTS(75) = 50
CONSTS(76) = 3.5e-12
CONSTS(77) = 1
CONSTS(78) = 4
CONSTS(79) = 0.034
CONSTS(80) = 0.012
CONSTS(81) = 0.0009
CONSTS(82) = 0.0001
CONSTS(83) = -73.5049651455872
CONSTS(84) = 0.004914
CONSTS(85) = 5.38843941249284e-5
CONSTS(86) = 153.604280337996
CONSTS(87) = 10.5731241425458
CONSTS(88) = 1
CONSTS(89) = 1
CONSTS(90) = 8.314
CONSTS(91) = 96.48534
CONSTS(92) = 2
CONSTS(93) = 1
CONSTS(94) = 1
CONSTS(95) = 8.314
CONSTS(96) = 96.48534
CONSTS(97) = 2
CONSTS(98) = 1
CONSTS(99) = 1
CONSTS(100) = 8.314
CONSTS(101) = 96.48534
CONSTS(102) = 1
CONSTS(103) = 1
CONSTS(104) = 2
CONSTS(105) = 8.314
CONSTS(106) = 96.48534
CONSTS(107) = 1
CONSTS(108) = 1
CONSTS(109) = 2
CONSTS(110) = 8.314
CONSTS(111) = 96.48534
CONSTS(112) = 1
CONSTS(113) = 2
CONSTS(114) = 1
CONSTS(115) = 8.314
CONSTS(116) = 96.48534
CONSTS(117) = 1
CONSTS(118) = 2
CONSTS(119) = 1
CONSTS(120) = 8.314
CONSTS(121) = 96.48534
CONSTS(122) = 1
CONSTS(123) = 1
CONSTS(124) = 2
CONSTS(125) = CONSTS(13) ** (CONSTS(9) - CONSTS(11))/10.0000
CONSTS(126) = (( CONSTS(90)*CONSTS(9))/( CONSTS(92)*CONSTS(91)))*log(CONSTS(14)/CONSTS(4))
CONSTS(127) = 0.00000
CONSTS(128) = 0.00000
CONSTS(129) = CONSTS(13) ** (CONSTS(9) - CONSTS(12))/10.0000
CONSTS(130) = CONSTS(13) ** (CONSTS(9) - CONSTS(12))/10.0000
CONSTS(131) = CONSTS(35) ** (CONSTS(9) - CONSTS(34))/10.0000
CONSTS(132) = CONSTS(35) ** (CONSTS(9) - CONSTS(34))/10.0000
CONSTS(133) = (( CONSTS(105)*CONSTS(9))/( CONSTS(107)*CONSTS(106)))*log(CONSTS(36)/CONSTS(5))
CONSTS(134) = 4.00000*CONSTS(52)*CONSTS(4)
CONSTS(135) = 3.00000*CONSTS(52)*CONSTS(4)
CONSTS(136) = 2.00000*CONSTS(52)*CONSTS(4)
CONSTS(137) = 1.00000*CONSTS(52)*CONSTS(4)
CONSTS(138) = 4.00000*CONSTS(53)
CONSTS(139) = 3.00000*CONSTS(53)
CONSTS(140) = 2.00000*CONSTS(53)
CONSTS(141) = 1.00000*CONSTS(53)
CONSTS(142) = 4.00000*CONSTS(52)*CONSTS(4)
CONSTS(143) = 3.00000*CONSTS(52)*CONSTS(4)
CONSTS(144) = 2.00000*CONSTS(52)*CONSTS(4)
CONSTS(145) = 1.00000*CONSTS(52)*CONSTS(4)
CONSTS(146) = 4.00000*CONSTS(54)
CONSTS(147) = 3.00000*CONSTS(54)
CONSTS(148) = 2.00000*CONSTS(54)
CONSTS(149) = 1.00000*CONSTS(54)
CONSTS(150) = (( CONSTS(110)*CONSTS(9))/( CONSTS(112)*CONSTS(111)))*log(CONSTS(58)/CONSTS(6))
CONSTS(151) = (( CONSTS(95)*CONSTS(9))/( CONSTS(97)*CONSTS(96)))*log(CONSTS(14)/CONSTS(4))
CONSTS(152) = (( CONSTS(100)*CONSTS(9))/( CONSTS(102)*CONSTS(101)))*log(CONSTS(36)/CONSTS(5))
CONSTS(153) = CONSTS(57) ** (CONSTS(9) - CONSTS(56))/10.0000
CONSTS(154) = (( CONSTS(115)*CONSTS(9))/( CONSTS(117)*CONSTS(116)))*log(CONSTS(58)/CONSTS(6))
CONSTS(155) = (( CONSTS(120)*CONSTS(9))/( CONSTS(122)*CONSTS(121)))*log(CONSTS(36)/CONSTS(5))
STATES(1) = CONSTS(28)
STATES(2) = CONSTS(51)
STATES(3) = CONSTS(64)
STATES(4) = CONSTS(30)
STATES(5) = CONSTS(31)
STATES(6) = CONSTS(37)
STATES(7) = CONSTS(38)
STATES(8) = CONSTS(15)
STATES(9) = CONSTS(16)
STATES(10) = CONSTS(17)
STATES(11) = CONSTS(18)
STATES(12) = CONSTS(19)
STATES(13) = CONSTS(20)
STATES(14) = CONSTS(21)
STATES(15) = CONSTS(22)
STATES(16) = CONSTS(23)
STATES(17) = CONSTS(24)
STATES(18) = CONSTS(25)
STATES(19) = CONSTS(26)
STATES(20) = CONSTS(26)
STATES(21) = CONSTS(42)
STATES(22) = CONSTS(43)
STATES(23) = CONSTS(44)
STATES(24) = CONSTS(45)
STATES(25) = CONSTS(46)
STATES(26) = CONSTS(47)
STATES(27) = CONSTS(48)
STATES(28) = CONSTS(49)
STATES(29) = CONSTS(50)
STATES(30) = CONSTS(59)
STATES(31) = CONSTS(60)
STATES(32) = CONSTS(61)
STATES(33) = CONSTS(62)
STATES(34) = CONSTS(63)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = TERNRY(VOI.LT.CONSTS(1), CONSTS(2), CONSTS(3))
ALGBRC(4) = 1.00000/(1.00000+EXP(- (ALGBRC(1)+60.5000)/5.30000))
RATES(4) = ((ALGBRC(4) - STATES(4))/CONSTS(32))*CONSTS(129)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (ALGBRC(1)+43.0000)/17.3600))
RATES(6) = ((ALGBRC(6) - STATES(6))/CONSTS(39))*CONSTS(131)
ALGBRC(7) = 1.00000/(1.00000+EXP((ALGBRC(1) - 44.9000)/12.0096))
RATES(7) = ((ALGBRC(7) - STATES(7))/CONSTS(40))*CONSTS(132)
ALGBRC(5) = 1.00000/(1.00000+EXP((ALGBRC(1)+75.5000)/4.00000))
ALGBRC(12) = 0.381170*(8.60000+ 14.7000*EXP(- ALGBRC(1)+50.0000 ** 2.00000/( 900.000*1.00000)))
RATES(5) = ((ALGBRC(5) - STATES(5))/ALGBRC(12))*CONSTS(130)
ALGBRC(55) = STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(2)
ALGBRC(8) = 1.00000*EXP(( 8.47188*ALGBRC(1))/( 1.00000*CONSTS(9)))
ALGBRC(34) = 0.0443240*ALGBRC(8)
ALGBRC(13) = 1.00000*EXP(( - 7.77556*ALGBRC(1))/( 1.00000*CONSTS(9)))
ALGBRC(52) = 0.000106000*ALGBRC(13)
RATES(25) = ( - (CONSTS(138)+ALGBRC(34))*STATES(25))/ALGBRC(55)+( CONSTS(137)*STATES(24))/ALGBRC(55)+( ALGBRC(52)*STATES(2))/ALGBRC(55)
ALGBRC(29) = 0.000781000*ALGBRC(8)
ALGBRC(48) = 0.0953120*ALGBRC(13)
RATES(24) = ( - (CONSTS(139)+ALGBRC(29)+CONSTS(137))*STATES(24))/ALGBRC(55)+( CONSTS(136)*STATES(23))/ALGBRC(55)+( ALGBRC(48)*STATES(29))/ALGBRC(55)+( CONSTS(138)*STATES(25))/ALGBRC(55)
ALGBRC(25) = 2.81000e-05*ALGBRC(8)
ALGBRC(44) = 32.6594*ALGBRC(13)
RATES(23) = ( - (CONSTS(140)+ALGBRC(25)+CONSTS(136))*STATES(23))/ALGBRC(55)+( CONSTS(135)*STATES(22))/ALGBRC(55)+( ALGBRC(44)*STATES(28))/ALGBRC(55)+( CONSTS(139)*STATES(24))/ALGBRC(55)
ALGBRC(21) = 0.000869000*ALGBRC(8)
ALGBRC(40) = 144.174*ALGBRC(13)
RATES(22) = ( - (CONSTS(141)+ALGBRC(21)+CONSTS(135))*STATES(22))/ALGBRC(55)+( CONSTS(134)*STATES(21))/ALGBRC(55)+( ALGBRC(40)*STATES(27))/ALGBRC(55)+( CONSTS(140)*STATES(23))/ALGBRC(55)
ALGBRC(17) = 0.0216200*ALGBRC(8)
ALGBRC(36) = 318.108*ALGBRC(13)
RATES(21) = ( - (CONSTS(134)+ALGBRC(17))*STATES(21))/ALGBRC(55)+( CONSTS(141)*STATES(22))/ALGBRC(55)+( ALGBRC(36)*STATES(26))/ALGBRC(55)
RATES(2) = ( - (CONSTS(146)+ALGBRC(52))*STATES(2))/ALGBRC(55)+( CONSTS(145)*STATES(29))/ALGBRC(55)+( ALGBRC(34)*STATES(25))/ALGBRC(55)
RATES(29) = ( - (CONSTS(147)+ALGBRC(48)+CONSTS(145))*STATES(29))/ALGBRC(55)+( CONSTS(144)*STATES(28))/ALGBRC(55)+( ALGBRC(29)*STATES(24))/ALGBRC(55)+( CONSTS(146)*STATES(2))/ALGBRC(55)
RATES(28) = ( - (CONSTS(148)+ALGBRC(44)+CONSTS(144))*STATES(28))/ALGBRC(55)+( CONSTS(143)*STATES(27))/ALGBRC(55)+( ALGBRC(25)*STATES(23))/ALGBRC(55)+( CONSTS(147)*STATES(29))/ALGBRC(55)
RATES(27) = ( - (CONSTS(149)+ALGBRC(40)+CONSTS(143))*STATES(27))/ALGBRC(55)+( CONSTS(142)*STATES(26))/ALGBRC(55)+( ALGBRC(21)*STATES(22))/ALGBRC(55)+( CONSTS(148)*STATES(28))/ALGBRC(55)
RATES(26) = ( - (CONSTS(142)+ALGBRC(36))*STATES(26))/ALGBRC(55)+( CONSTS(149)*STATES(27))/ALGBRC(55)+( ALGBRC(17)*STATES(21))/ALGBRC(55)
ALGBRC(56) = STATES(3)+STATES(30)+STATES(31)+STATES(32)+STATES(33)+STATES(34)
ALGBRC(26) = CONSTS(153)*0.554320*EXP(- 0.0990740+ 0.0364410*ALGBRC(1))
ALGBRC(30) = CONSTS(153)*0.000525480*EXP(- 0.0691020+ 0.00319450*ALGBRC(1))
RATES(32) = ( - ALGBRC(30)*STATES(32))/ALGBRC(56)+( ALGBRC(26)*STATES(31))/ALGBRC(56)
ALGBRC(22) = CONSTS(153)*3.15660*EXP(0.363520+ 0.0771930*ALGBRC(1))
ALGBRC(32) = CONSTS(153)*1.44960*EXP(- 0.156600+ 0.0583530*ALGBRC(1))
RATES(31) = ( - (ALGBRC(32)+ALGBRC(26))*STATES(31))/ALGBRC(56)+( ALGBRC(22)*STATES(30))/ALGBRC(56)+( ALGBRC(30)*STATES(32))/ALGBRC(56)
ALGBRC(18) = CONSTS(153)*2.39100*EXP(- 13.3350 - 0.252890*ALGBRC(1))
ALGBRC(37) = CONSTS(153)*1.53290*EXP(0.00931930+ 0.0410750*ALGBRC(1))
ALGBRC(49) = CONSTS(153)*1.90460*EXP(- 2.48400+ 0.0204060*ALGBRC(1))
ALGBRC(53) = CONSTS(153)*0.000216880*EXP(- 0.0634380+ 0.00466830*ALGBRC(1))
RATES(30) = ( - (ALGBRC(22)+ALGBRC(37)+ALGBRC(53))*STATES(30))/ALGBRC(56)+( ALGBRC(18)*STATES(3))/ALGBRC(56)+( ALGBRC(32)*STATES(31))/ALGBRC(56)+( ALGBRC(49)*STATES(33))/ALGBRC(56)
ALGBRC(14) = CONSTS(153)*0.120520*EXP(- 9.60280+ 0.0830250*ALGBRC(1))
ALGBRC(41) = CONSTS(153)*1.61640*EXP(0.307630+ 0.00605350*ALGBRC(1))
RATES(3) = ( - (ALGBRC(18)+ALGBRC(41))*STATES(3))/ALGBRC(56)+( ALGBRC(37)*STATES(30))/ALGBRC(56)+( ALGBRC(14)*STATES(33))/ALGBRC(56)
ALGBRC(9) = CONSTS(153)*0.00392390*EXP(2.67930+ 0.00614680*ALGBRC(1))
ALGBRC(45) = CONSTS(153)*0.0277350*EXP(0.0514900 - 0.0468650*ALGBRC(1))
RATES(34) = ( - ALGBRC(9)*STATES(34))/ALGBRC(56)+( ALGBRC(45)*STATES(33))/ALGBRC(56)
RATES(33) = ( - (ALGBRC(14)+ALGBRC(45)+ALGBRC(49))*STATES(33))/ALGBRC(56)+( ALGBRC(9)*STATES(34))/ALGBRC(56)+( ALGBRC(53)*STATES(30))/ALGBRC(56)+( ALGBRC(41)*STATES(3))/ALGBRC(56)
ALGBRC(3) = CONSTS(125)*0.731000*EXP(ALGBRC(1)/30.0000)
ALGBRC(24) = 2.00000*ALGBRC(3)
ALGBRC(28) = 1.00000*ALGBRC(3)
ALGBRC(11) = CONSTS(125)*0.214900*EXP(- ALGBRC(1)/40.0000)
ALGBRC(42) = 3.00000*ALGBRC(11)
ALGBRC(46) = 4.00000*ALGBRC(11)
ALGBRC(54) = CONSTS(125)*0.00595600*EXP(- ALGBRC(1)/40.0000)
ALGBRC(50) = CONSTS(125)*0.474200*EXP(ALGBRC(1)/10.0000)
ALGBRC(60) = CONSTS(125)*0.00232000*EXP(- ALGBRC(1)/280.000)
ALGBRC(58) = CONSTS(125)*0.0121300*EXP(ALGBRC(1)/500.000)
ALGBRC(62) = ( ALGBRC(46)*ALGBRC(58)*ALGBRC(54))/( ALGBRC(28)*ALGBRC(60))
ALGBRC(59) = CONSTS(125)*0.0219700*EXP(ALGBRC(1)/500.000)
ALGBRC(57) = CONSTS(125)*0.0140700*EXP(- ALGBRC(1)/300.000)
ALGBRC(61) = ( ALGBRC(46)*ALGBRC(57)*ALGBRC(50))/( ALGBRC(28)*ALGBRC(59))
ALGBRC(65) = STATES(11)+STATES(10)+STATES(9)+STATES(8)+STATES(15)+STATES(14)+STATES(13)+STATES(12)+STATES(1)+STATES(16)+STATES(17)+STATES(18)+STATES(20)+STATES(19)
RATES(11) = (( ALGBRC(24)*STATES(10))/ALGBRC(65) - ( (ALGBRC(28)+ALGBRC(42)+ALGBRC(54)+ALGBRC(50)+CONSTS(128))*STATES(11))/ALGBRC(65))+( ALGBRC(46)*STATES(1))/ALGBRC(65)+( ALGBRC(61)*STATES(18))/ALGBRC(65)+( ALGBRC(62)*STATES(17))/ALGBRC(65)+( CONSTS(127)*STATES(15))/ALGBRC(65)
ALGBRC(20) = 3.00000*ALGBRC(3)
ALGBRC(38) = 2.00000*ALGBRC(11)
RATES(10) = (( ALGBRC(20)*STATES(9))/ALGBRC(65) - ( (ALGBRC(24)+ALGBRC(38)+CONSTS(128))*STATES(10))/ALGBRC(65))+( ALGBRC(42)*STATES(11))/ALGBRC(65)+( CONSTS(127)*STATES(14))/ALGBRC(65)
ALGBRC(16) = 4.00000*ALGBRC(3)
ALGBRC(33) = 1.00000*ALGBRC(11)
RATES(9) = (( ALGBRC(16)*STATES(8))/ALGBRC(65) - ( (ALGBRC(20)+ALGBRC(33)+CONSTS(128))*STATES(9))/ALGBRC(65))+( ALGBRC(38)*STATES(10))/ALGBRC(65)+( CONSTS(127)*STATES(13))/ALGBRC(65)
RATES(8) = ( - (ALGBRC(16)+CONSTS(128))*STATES(8))/ALGBRC(65)+( ALGBRC(33)*STATES(9))/ALGBRC(65)+( CONSTS(127)*STATES(12))/ALGBRC(65)
RATES(15) = ((( CONSTS(128)*STATES(11))/ALGBRC(65)+( ALGBRC(24)*STATES(14))/ALGBRC(65)) - ( (ALGBRC(28)+ALGBRC(42)+ALGBRC(50)+ALGBRC(54)+CONSTS(127))*STATES(15))/ALGBRC(65))+( ALGBRC(46)*STATES(16))/ALGBRC(65)+( ALGBRC(61)*STATES(19))/ALGBRC(65)+( ALGBRC(62)*STATES(20))/ALGBRC(65)
RATES(14) = ((( CONSTS(128)*STATES(10))/ALGBRC(65)+( ALGBRC(20)*STATES(13))/ALGBRC(65)) - ( (ALGBRC(24)+ALGBRC(38)+CONSTS(127))*STATES(14))/ALGBRC(65))+( ALGBRC(42)*STATES(15))/ALGBRC(65)
RATES(13) = ((( CONSTS(128)*STATES(9))/ALGBRC(65)+( ALGBRC(16)*STATES(12))/ALGBRC(65)) - ( (ALGBRC(20)+ALGBRC(33)+CONSTS(127))*STATES(13))/ALGBRC(65))+( ALGBRC(38)*STATES(14))/ALGBRC(65)
RATES(12) = (( CONSTS(128)*STATES(8))/ALGBRC(65) - ( (ALGBRC(16)+CONSTS(127))*STATES(12))/ALGBRC(65))+( ALGBRC(33)*STATES(13))/ALGBRC(65)
RATES(1) = (( ALGBRC(28)*STATES(11))/ALGBRC(65) - ( (ALGBRC(46)+ALGBRC(59)+ALGBRC(60)+CONSTS(128))*STATES(1))/ALGBRC(65))+( ALGBRC(57)*STATES(18))/ALGBRC(65)+( ALGBRC(58)*STATES(17))/ALGBRC(65)+( CONSTS(127)*STATES(16))/ALGBRC(65)
RATES(16) = ((( CONSTS(128)*STATES(1))/ALGBRC(65)+( ALGBRC(28)*STATES(15))/ALGBRC(65)) - ( (ALGBRC(46)+ALGBRC(59)+ALGBRC(60)+CONSTS(127))*STATES(16))/ALGBRC(65))+( ALGBRC(57)*STATES(19))/ALGBRC(65)+( ALGBRC(58)*STATES(20))/ALGBRC(65)
ALGBRC(64) = ALGBRC(60)
ALGBRC(63) = ( ALGBRC(58)*ALGBRC(59))/ALGBRC(57)
RATES(17) = ((( ALGBRC(54)*STATES(11))/ALGBRC(65)+( ALGBRC(60)*STATES(1))/ALGBRC(65)+( ALGBRC(64)*STATES(18))/ALGBRC(65)) - ( (ALGBRC(62)+ALGBRC(58)+ALGBRC(63)+CONSTS(128))*STATES(17))/ALGBRC(65))+( CONSTS(127)*STATES(20))/ALGBRC(65)
RATES(18) = ((( ALGBRC(50)*STATES(11))/ALGBRC(65)+( ALGBRC(59)*STATES(1))/ALGBRC(65)) - ( (ALGBRC(61)+ALGBRC(57)+ALGBRC(64)+CONSTS(128))*STATES(18))/ALGBRC(65))+( ALGBRC(63)*STATES(17))/ALGBRC(65)+( CONSTS(127)*STATES(19))/ALGBRC(65)
RATES(20) = (( CONSTS(128)*STATES(17))/ALGBRC(65)+( ALGBRC(54)*STATES(15))/ALGBRC(65)+( ALGBRC(60)*STATES(16))/ALGBRC(65)+( ALGBRC(64)*STATES(19))/ALGBRC(65)) - ( (ALGBRC(62)+ALGBRC(58)+ALGBRC(63)+CONSTS(127))*STATES(20))/ALGBRC(65)
RATES(19) = ((( CONSTS(128)*STATES(18))/ALGBRC(65)+( ALGBRC(50)*STATES(15))/ALGBRC(65)+( ALGBRC(59)*STATES(16))/ALGBRC(65)) - ( (ALGBRC(61)+ALGBRC(57)+ALGBRC(64)+CONSTS(127))*STATES(19))/ALGBRC(65))+( ALGBRC(63)*STATES(20))/ALGBRC(65)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = TERNRY(VOI.LT.CONSTS(1), CONSTS(2), CONSTS(3))
ALGBRC(4) = 1.00000/(1.00000+EXP(- (ALGBRC(1)+60.5000)/5.30000))
ALGBRC(6) = 1.00000/(1.00000+EXP(- (ALGBRC(1)+43.0000)/17.3600))
ALGBRC(7) = 1.00000/(1.00000+EXP((ALGBRC(1) - 44.9000)/12.0096))
ALGBRC(5) = 1.00000/(1.00000+EXP((ALGBRC(1)+75.5000)/4.00000))
ALGBRC(12) = 0.381170*(8.60000+ 14.7000*EXP(- ALGBRC(1)+50.0000 ** 2.00000/( 900.000*1.00000)))
ALGBRC(55) = STATES(21)+STATES(22)+STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(2)
ALGBRC(8) = 1.00000*EXP(( 8.47188*ALGBRC(1))/( 1.00000*CONSTS(9)))
ALGBRC(34) = 0.0443240*ALGBRC(8)
ALGBRC(13) = 1.00000*EXP(( - 7.77556*ALGBRC(1))/( 1.00000*CONSTS(9)))
ALGBRC(52) = 0.000106000*ALGBRC(13)
ALGBRC(29) = 0.000781000*ALGBRC(8)
ALGBRC(48) = 0.0953120*ALGBRC(13)
ALGBRC(25) = 2.81000e-05*ALGBRC(8)
ALGBRC(44) = 32.6594*ALGBRC(13)
ALGBRC(21) = 0.000869000*ALGBRC(8)
ALGBRC(40) = 144.174*ALGBRC(13)
ALGBRC(17) = 0.0216200*ALGBRC(8)
ALGBRC(36) = 318.108*ALGBRC(13)
ALGBRC(56) = STATES(3)+STATES(30)+STATES(31)+STATES(32)+STATES(33)+STATES(34)
ALGBRC(26) = CONSTS(153)*0.554320*EXP(- 0.0990740+ 0.0364410*ALGBRC(1))
ALGBRC(30) = CONSTS(153)*0.000525480*EXP(- 0.0691020+ 0.00319450*ALGBRC(1))
ALGBRC(22) = CONSTS(153)*3.15660*EXP(0.363520+ 0.0771930*ALGBRC(1))
ALGBRC(32) = CONSTS(153)*1.44960*EXP(- 0.156600+ 0.0583530*ALGBRC(1))
ALGBRC(18) = CONSTS(153)*2.39100*EXP(- 13.3350 - 0.252890*ALGBRC(1))
ALGBRC(37) = CONSTS(153)*1.53290*EXP(0.00931930+ 0.0410750*ALGBRC(1))
ALGBRC(49) = CONSTS(153)*1.90460*EXP(- 2.48400+ 0.0204060*ALGBRC(1))
ALGBRC(53) = CONSTS(153)*0.000216880*EXP(- 0.0634380+ 0.00466830*ALGBRC(1))
ALGBRC(14) = CONSTS(153)*0.120520*EXP(- 9.60280+ 0.0830250*ALGBRC(1))
ALGBRC(41) = CONSTS(153)*1.61640*EXP(0.307630+ 0.00605350*ALGBRC(1))
ALGBRC(9) = CONSTS(153)*0.00392390*EXP(2.67930+ 0.00614680*ALGBRC(1))
ALGBRC(45) = CONSTS(153)*0.0277350*EXP(0.0514900 - 0.0468650*ALGBRC(1))
ALGBRC(3) = CONSTS(125)*0.731000*EXP(ALGBRC(1)/30.0000)
ALGBRC(24) = 2.00000*ALGBRC(3)
ALGBRC(28) = 1.00000*ALGBRC(3)
ALGBRC(11) = CONSTS(125)*0.214900*EXP(- ALGBRC(1)/40.0000)
ALGBRC(42) = 3.00000*ALGBRC(11)
ALGBRC(46) = 4.00000*ALGBRC(11)
ALGBRC(54) = CONSTS(125)*0.00595600*EXP(- ALGBRC(1)/40.0000)
ALGBRC(50) = CONSTS(125)*0.474200*EXP(ALGBRC(1)/10.0000)
ALGBRC(60) = CONSTS(125)*0.00232000*EXP(- ALGBRC(1)/280.000)
ALGBRC(58) = CONSTS(125)*0.0121300*EXP(ALGBRC(1)/500.000)
ALGBRC(62) = ( ALGBRC(46)*ALGBRC(58)*ALGBRC(54))/( ALGBRC(28)*ALGBRC(60))
ALGBRC(59) = CONSTS(125)*0.0219700*EXP(ALGBRC(1)/500.000)
ALGBRC(57) = CONSTS(125)*0.0140700*EXP(- ALGBRC(1)/300.000)
ALGBRC(61) = ( ALGBRC(46)*ALGBRC(57)*ALGBRC(50))/( ALGBRC(28)*ALGBRC(59))
ALGBRC(65) = STATES(11)+STATES(10)+STATES(9)+STATES(8)+STATES(15)+STATES(14)+STATES(13)+STATES(12)+STATES(1)+STATES(16)+STATES(17)+STATES(18)+STATES(20)+STATES(19)
ALGBRC(20) = 3.00000*ALGBRC(3)
ALGBRC(38) = 2.00000*ALGBRC(11)
ALGBRC(16) = 4.00000*ALGBRC(3)
ALGBRC(33) = 1.00000*ALGBRC(11)
ALGBRC(64) = ALGBRC(60)
ALGBRC(63) = ( ALGBRC(58)*ALGBRC(59))/ALGBRC(57)
ALGBRC(2) = STATES(4)*STATES(5)
ALGBRC(10) = STATES(6)*STATES(7)
ALGBRC(15) = ( CONSTS(69)*( EXP(( CONSTS(68)*ALGBRC(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))*CONSTS(6) ** 3.00000*CONSTS(14) - 2.50000*EXP(( (CONSTS(68) - 1.00000)*ALGBRC(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))*CONSTS(58) ** 3.00000*CONSTS(4)))/( (CONSTS(66) ** 3.00000+CONSTS(58) ** 3.00000)*(CONSTS(65)+CONSTS(14))*(1.00000+ CONSTS(67)*EXP(( (CONSTS(68) - 1.00000)*ALGBRC(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))))
ALGBRC(19) = ( CONSTS(72)*CONSTS(36)*CONSTS(6))/( (CONSTS(36)+CONSTS(70))*(CONSTS(71)+CONSTS(6))*(1.00000+ 0.124500*EXP(( - 0.100000*ALGBRC(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))+ 0.0353000*EXP(( - ALGBRC(1)*CONSTS(10))/( CONSTS(8)*CONSTS(9)))))
ALGBRC(23) = CONSTS(7)*STATES(1)*(ALGBRC(1) - CONSTS(126))
ALGBRC(27) = CONSTS(29)*ALGBRC(2)*(ALGBRC(1) - CONSTS(151))
ALGBRC(31) = CONSTS(33)*ALGBRC(10)*(ALGBRC(1) - CONSTS(152))
ALGBRC(35) = CONSTS(41)*STATES(2)*(ALGBRC(1) - CONSTS(133))
ALGBRC(39) = CONSTS(55)*STATES(3)*(ALGBRC(1) - CONSTS(150))
ALGBRC(43) = CONSTS(73)*CONSTS(88)*(ALGBRC(1) - CONSTS(154))
ALGBRC(47) = CONSTS(74)*CONSTS(89)*(ALGBRC(1) - CONSTS(155))
ALGBRC(51) = ALGBRC(23)+ALGBRC(27)+ALGBRC(31)+ALGBRC(35)+ALGBRC(39)+ALGBRC(15)+ALGBRC(19)+ALGBRC(43)+ALGBRC(47)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END