C
C There are a total of 80 entries in the algebraic variable array.
C There are a total of 38 entries in each of the rate and state variable arrays.
C There are a total of 65 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is V in component membrane (millivolt).
C CONSTS(1) is R in component membrane (joule_per_mole_kelvin).
C CONSTS(2) is T in component membrane (kelvin).
C CONSTS(3) is F in component membrane (coulomb_per_millimole).
C CONSTS(4) is Cm in component membrane (picoF).
C ALGBRC(14) is I_stim in component I_stim (picoA).
C ALGBRC(54) is i_Na in component sodium_current (picoA).
C ALGBRC(22) is i_Ca_L in component L_type_Ca_channel (picoA).
C ALGBRC(30) is i_Ca_T in component T_type_Ca_channel (picoA).
C ALGBRC(38) is i_Kv in component voltage_dep_K_channel (picoA).
C ALGBRC(46) is i_BK in component Calcium_voltage_activated_potassium_channel (picoA).
C ALGBRC(58) is i_NCX in component Na_Ca_exchanger (picoA).
C ALGBRC(61) is i_NaK in component sodium_potassium_pump (picoA).
C ALGBRC(73) is i_NsK in component non_specific_K_current (picoA).
C ALGBRC(67) is i_NsNa in component non_specific_Na_current (picoA).
C ALGBRC(9) is local_time in component I_stim (millisecond).
C CONSTS(5) is period in component I_stim (millisecond).
C ALGBRC(1) is stim_start in component I_stim (millisecond).
C CONSTS(6) is Gcouple in component I_stim (nanoS).
C CONSTS(7) is t_ICCpeak in component I_stim (millisecond).
C CONSTS(8) is t_ICCplateau in component I_stim (millisecond).
C CONSTS(9) is t_slope in component I_stim (millisecond).
C CONSTS(10) is f_1 in component I_stim (millisecond).
C CONSTS(11) is f_2 in component I_stim (millisecond).
C CONSTS(12) is V_ICCrest in component I_stim (millivolt).
C CONSTS(13) is V_ICCamp in component I_stim (millivolt).
C ALGBRC(18) is E_Ca in component L_type_Ca_channel (millivolt).
C CONSTS(14) is g_CaL in component L_type_Ca_channel (nanoS).
C CONSTS(15) is Cao in component ionic_concentrations (millimolar).
C STATES(2) is Cai in component ionic_concentrations (millimolar).
C STATES(3) is O_CaL in component L_type_Ca_channel_states (dimensionless).
C STATES(4) is C0 in component L_type_Ca_channel_states (dimensionless).
C STATES(5) is C1 in component L_type_Ca_channel_states (dimensionless).
C STATES(6) is C2 in component L_type_Ca_channel_states (dimensionless).
C STATES(7) is C3 in component L_type_Ca_channel_states (dimensionless).
C STATES(8) is C0Ca in component L_type_Ca_channel_states (dimensionless).
C STATES(9) is C1Ca in component L_type_Ca_channel_states (dimensionless).
C STATES(10) is C2Ca in component L_type_Ca_channel_states (dimensionless).
C STATES(11) is C3Ca in component L_type_Ca_channel_states (dimensionless).
C STATES(12) is ICa in component L_type_Ca_channel_states (dimensionless).
C STATES(13) is IVs in component L_type_Ca_channel_states (dimensionless).
C STATES(14) is IVf in component L_type_Ca_channel_states (dimensionless).
C STATES(15) is IVfCa in component L_type_Ca_channel_states (dimensionless).
C STATES(16) is IVsCa in component L_type_Ca_channel_states (dimensionless).
C ALGBRC(2) is alpha in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(10) is beta in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(15) is alpha_0 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(19) is alpha_1 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(23) is alpha_2 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(27) is alpha_3 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(31) is beta_0 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(35) is beta_1 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(39) is beta_2 in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(43) is beta_3 in component L_type_Ca_channel_states (per_millisecond).
C CONSTS(54) is delta in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(78) is theta in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(51) is phi_s in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(47) is phi_f in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(71) is omega_s in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(68) is omega_f in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(76) is omega_fs in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(74) is omega_sf in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(65) is psi_s in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(62) is psi_f in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(59) is xi_s in component L_type_Ca_channel_states (per_millisecond).
C ALGBRC(55) is xi_f in component L_type_Ca_channel_states (per_millisecond).
C CONSTS(45) is T_correction_Ca in component L_type_Ca_channel_states (dimensionless).
C CONSTS(16) is Q10Ca in component L_type_Ca_channel_states (dimensionless).
C ALGBRC(80) is norm in component L_type_Ca_channel_states (dimensionless).
C CONSTS(17) is g_CaT in component T_type_Ca_channel (nanoS).
C ALGBRC(26) is E_Ca in component T_type_Ca_channel (millivolt).
C STATES(17) is d_CaT in component T_type_Ca_channel_d_gate (dimensionless).
C STATES(18) is f_CaT in component T_type_Ca_channel_f_gate (dimensionless).
C ALGBRC(3) is d_CaT_inf in component T_type_Ca_channel_d_gate (dimensionless).
C CONSTS(55) is tau_d_CaT in component T_type_Ca_channel_d_gate (millisecond).
C CONSTS(46) is T_correction_Ca in component T_type_Ca_channel_d_gate (dimensionless).
C CONSTS(18) is Q10Ca in component T_type_Ca_channel_d_gate (dimensionless).
C ALGBRC(4) is f_CaT_inf in component T_type_Ca_channel_f_gate (dimensionless).
C ALGBRC(11) is tau_f_CaT in component T_type_Ca_channel_f_gate (millisecond).
C CONSTS(47) is T_correction_Ca in component T_type_Ca_channel_f_gate (dimensionless).
C CONSTS(19) is Q10Ca in component T_type_Ca_channel_f_gate (dimensionless).
C CONSTS(20) is g_Kv in component voltage_dep_K_channel (nanoS).
C ALGBRC(34) is E_K in component voltage_dep_K_channel (millivolt).
C STATES(19) is x_Kv in component voltage_dep_K_channel_x_gate (dimensionless).
C STATES(20) is y_Kv in component voltage_dep_K_channel_y_gate (dimensionless).
C CONSTS(21) is Ko in component ionic_concentrations (millimolar).
C STATES(21) is Ki in component ionic_concentrations (millimolar).
C ALGBRC(5) is x_Kv_inf in component voltage_dep_K_channel_x_gate (dimensionless).
C CONSTS(56) is tau_x_Kv in component voltage_dep_K_channel_x_gate (millisecond).
C CONSTS(48) is T_correction_K in component voltage_dep_K_channel_x_gate (dimensionless).
C CONSTS(22) is Q10K in component voltage_dep_K_channel_x_gate (dimensionless).
C ALGBRC(6) is y_Kv_inf in component voltage_dep_K_channel_y_gate (dimensionless).
C CONSTS(57) is tau_y_Kv in component voltage_dep_K_channel_y_gate (millisecond).
C CONSTS(49) is T_correction_K in component voltage_dep_K_channel_y_gate (dimensionless).
C CONSTS(23) is Q10K in component voltage_dep_K_channel_y_gate (dimensionless).
C ALGBRC(42) is E_K in component Calcium_voltage_activated_potassium_channel (millivolt).
C CONSTS(24) is g_BK in component Calcium_voltage_activated_potassium_channel (nanoS).
C STATES(22) is O4 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(23) is C0 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(24) is C1 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(25) is C2 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(26) is C3 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(27) is C4 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(28) is O0 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(29) is O1 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(30) is O2 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C STATES(31) is O3 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C ALGBRC(7) is alpha in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(12) is beta in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(50) is k_on in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(51) is k_off_C in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(52) is k_off_O in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(16) is k_C0O0 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(20) is k_C1O1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(24) is k_C2O2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(28) is k_C3O3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(32) is k_C4O4 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(36) is k_O0C0 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(40) is k_O1C1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(44) is k_O2C2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(48) is k_O3C3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(52) is k_O4C4 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(56) is k_C0C1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(60) is k_C1C2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(63) is k_C2C3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(66) is k_C3C4 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(58) is k_C4C3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(59) is k_C3C2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(60) is k_C2C1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(61) is k_C1C0 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(69) is k_O0O1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(72) is k_O1O2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(75) is k_O2O3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(77) is k_O3O4 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(62) is k_O4O3 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(63) is k_O3O2 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(64) is k_O2O1 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C CONSTS(65) is k_O1O0 in component Calcium_voltage_activated_potassium_channel_states (per_millisecond).
C ALGBRC(79) is norm in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C ALGBRC(50) is E_Na in component sodium_current (millivolt).
C CONSTS(25) is g_Na in component sodium_current (nanoS).
C CONSTS(26) is Nao in component ionic_concentrations (millimolar).
C STATES(32) is Nai in component ionic_concentrations (millimolar).
C STATES(33) is O_Na in component Na_channel_states (dimensionless).
C STATES(34) is C1 in component Na_channel_states (dimensionless).
C STATES(35) is C2 in component Na_channel_states (dimensionless).
C STATES(36) is C3 in component Na_channel_states (dimensionless).
C STATES(37) is I1 in component Na_channel_states (dimensionless).
C STATES(38) is I2 in component Na_channel_states (dimensionless).
C ALGBRC(8) is k_I2I1 in component Na_channel_states (per_millisecond).
C ALGBRC(13) is k_I1O in component Na_channel_states (per_millisecond).
C ALGBRC(17) is k_OC1 in component Na_channel_states (per_millisecond).
C ALGBRC(21) is k_C1C2 in component Na_channel_states (per_millisecond).
C ALGBRC(25) is k_C2C3 in component Na_channel_states (per_millisecond).
C ALGBRC(29) is k_C3C2 in component Na_channel_states (per_millisecond).
C ALGBRC(33) is k_C2C1 in component Na_channel_states (per_millisecond).
C ALGBRC(37) is k_C1O in component Na_channel_states (per_millisecond).
C ALGBRC(41) is k_OI1 in component Na_channel_states (per_millisecond).
C ALGBRC(45) is k_I1I2 in component Na_channel_states (per_millisecond).
C ALGBRC(49) is k_I1C1 in component Na_channel_states (per_millisecond).
C ALGBRC(53) is k_C1I1 in component Na_channel_states (per_millisecond).
C ALGBRC(57) is norm in component Na_channel_states (dimensionless).
C CONSTS(53) is T_correction_Na in component Na_channel_states (dimensionless).
C CONSTS(27) is Q10Na in component Na_channel_states (dimensionless).
C CONSTS(28) is P_NCX in component Na_Ca_exchanger (picoA).
C CONSTS(29) is K_mCa in component Na_Ca_exchanger (millimolar).
C CONSTS(30) is K_mNai in component Na_Ca_exchanger (millimolar).
C CONSTS(31) is k_sat in component Na_Ca_exchanger (dimensionless).
C CONSTS(32) is gamma in component Na_Ca_exchanger (dimensionless).
C CONSTS(33) is P_NaK in component sodium_potassium_pump (picoA).
C CONSTS(34) is K_mK in component sodium_potassium_pump (millimolar).
C CONSTS(35) is K_mNa in component sodium_potassium_pump (millimolar).
C CONSTS(36) is g_NsNa in component non_specific_Na_current (nanoS).
C ALGBRC(64) is E_Na in component non_specific_Na_current (millivolt).
C CONSTS(37) is g_NsK in component non_specific_K_current (nanoS).
C ALGBRC(70) is E_K in component non_specific_K_current (millivolt).
C CONSTS(38) is CRT_total in component ionic_concentrations (millimolar).
C CONSTS(39) is CaM_total in component ionic_concentrations (millimolar).
C CONSTS(40) is K_D_CRT in component ionic_concentrations (millimolar).
C CONSTS(41) is K_D_CaM in component ionic_concentrations (millimolar4).
C CONSTS(42) is n_CRT in component ionic_concentrations (dimensionless).
C CONSTS(43) is n_CaM in component ionic_concentrations (dimensionless).
C CONSTS(44) is V_myo in component ionic_concentrations (litre).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(7) is d/dt C3 in component L_type_Ca_channel_states (dimensionless).
C RATES(6) is d/dt C2 in component L_type_Ca_channel_states (dimensionless).
C RATES(5) is d/dt C1 in component L_type_Ca_channel_states (dimensionless).
C RATES(4) is d/dt C0 in component L_type_Ca_channel_states (dimensionless).
C RATES(11) is d/dt C3Ca in component L_type_Ca_channel_states (dimensionless).
C RATES(10) is d/dt C2Ca in component L_type_Ca_channel_states (dimensionless).
C RATES(9) is d/dt C1Ca in component L_type_Ca_channel_states (dimensionless).
C RATES(8) is d/dt C0Ca in component L_type_Ca_channel_states (dimensionless).
C RATES(3) is d/dt O_CaL in component L_type_Ca_channel_states (dimensionless).
C RATES(12) is d/dt ICa in component L_type_Ca_channel_states (dimensionless).
C RATES(13) is d/dt IVs in component L_type_Ca_channel_states (dimensionless).
C RATES(14) is d/dt IVf in component L_type_Ca_channel_states (dimensionless).
C RATES(16) is d/dt IVsCa in component L_type_Ca_channel_states (dimensionless).
C RATES(15) is d/dt IVfCa in component L_type_Ca_channel_states (dimensionless).
C RATES(17) is d/dt d_CaT in component T_type_Ca_channel_d_gate (dimensionless).
C RATES(18) is d/dt f_CaT in component T_type_Ca_channel_f_gate (dimensionless).
C RATES(19) is d/dt x_Kv in component voltage_dep_K_channel_x_gate (dimensionless).
C RATES(20) is d/dt y_Kv in component voltage_dep_K_channel_y_gate (dimensionless).
C RATES(27) is d/dt C4 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(26) is d/dt C3 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(25) is d/dt C2 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(24) is d/dt C1 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(23) is d/dt C0 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(22) is d/dt O4 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(31) is d/dt O3 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(30) is d/dt O2 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(29) is d/dt O1 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(28) is d/dt O0 in component Calcium_voltage_activated_potassium_channel_states (dimensionless).
C RATES(36) is d/dt C3 in component Na_channel_states (dimensionless).
C RATES(35) is d/dt C2 in component Na_channel_states (dimensionless).
C RATES(34) is d/dt C1 in component Na_channel_states (dimensionless).
C RATES(33) is d/dt O_Na in component Na_channel_states (dimensionless).
C RATES(38) is d/dt I2 in component Na_channel_states (dimensionless).
C RATES(37) is d/dt I1 in component Na_channel_states (dimensionless).
C RATES(32) is d/dt Nai in component ionic_concentrations (millimolar).
C RATES(21) is d/dt Ki in component ionic_concentrations (millimolar).
C RATES(2) is d/dt Cai in component ionic_concentrations (millimolar).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -73.5049651455872
      CONSTS(1) = 8.314
      CONSTS(2) = 310
      CONSTS(3) = 96.48534
      CONSTS(4) = 50
      CONSTS(5) = 10000
      CONSTS(6) = 2.6
      CONSTS(7) = 300
      CONSTS(8) = 9700
      CONSTS(9) = 600
      CONSTS(10) = 12000
      CONSTS(11) = 300
      CONSTS(12) = -57
      CONSTS(13) = 33.5
      CONSTS(14) = 1.44
      CONSTS(15) = 2
      STATES(2) = 5.38843941249284e-5
      STATES(3) = 3.88576045134351e-6
      STATES(4) = 0.815464741971086
      STATES(5) = 0.152399266235657
      STATES(6) = 0.0106805060777161
      STATES(7) = 0.000332673548872087
      STATES(8) = 0.0175888495282545
      STATES(9) = 0.00328711668724504
      STATES(10) = 0.000230369020877669
      STATES(11) = 7.1754726923539e-6
      STATES(12) = 8.38123983500905e-8
      STATES(13) = 1.1193313274705e-6
      STATES(14) = 4.0998751301597e-6
      STATES(15) = 8.84306615061238e-8
      STATES(16) = 2.41429816075123e-8
      CONSTS(16) = 2.1
      CONSTS(17) = 0.0425
      STATES(17) = 0.0791635737410974
      STATES(18) = 0.377831534375835
      CONSTS(18) = 2.1
      CONSTS(19) = 2.1
      CONSTS(20) = 1.0217
      STATES(19) = 0.14714161078933
      STATES(20) = 0.99994773314105
      CONSTS(21) = 5.4
      STATES(21) = 153.604280337996
      CONSTS(22) = 3.1
      CONSTS(23) = 3.1
      CONSTS(24) = 80
      STATES(22) = 1.82113764497095e-6
      STATES(23) = 0.48379087935899
      STATES(24) = 0.385183559520031
      STATES(25) = 0.115002824567753
      STATES(26) = 0.0152602714149774
      STATES(27) = 0.000759264410974374
      STATES(28) = 6.94960798375172e-7
      STATES(29) = 5.55636826398253e-8
      STATES(30) = 2.85143702125325e-8
      STATES(31) = 1.59832806123435e-6
      CONSTS(25) = 25.1
      CONSTS(26) = 140
      STATES(32) = 10.5731241425458
      STATES(33) = 1.96760342050475e-6
      STATES(34) = 0.0119443135223679
      STATES(35) = 0.0109545368437155
      STATES(36) = 0.973782548650071
      STATES(37) = 0.000126882921013389
      STATES(38) = 0.00318975045717667
      CONSTS(27) = 2.45
      CONSTS(28) = 1992.335
      CONSTS(29) = 1.38
      CONSTS(30) = 87.5
      CONSTS(31) = 0.1
      CONSTS(32) = 0.35
      CONSTS(33) = 9.26
      CONSTS(34) = 1
      CONSTS(35) = 40
      CONSTS(36) = 0.022488
      CONSTS(37) = 0.017512
      CONSTS(38) = 0.034
      CONSTS(39) = 0.012
      CONSTS(40) = 0.0009
      CONSTS(41) = 0.0001
      CONSTS(42) = 1
      CONSTS(43) = 4
      CONSTS(44) = 3.50000e-12
      CONSTS(45) = CONSTS(16) ** (CONSTS(2) - 310.000)/10.0000
      CONSTS(46) = CONSTS(18) ** (CONSTS(2) - 297.000)/10.0000
      CONSTS(47) = CONSTS(19) ** (CONSTS(2) - 297.000)/10.0000
      CONSTS(48) = CONSTS(22) ** (CONSTS(2) - 297.000)/10.0000
      CONSTS(49) = CONSTS(23) ** (CONSTS(2) - 297.000)/10.0000
      CONSTS(50) = 40633.0
      CONSTS(51) = 11.0000
      CONSTS(52) = 1.10000
      CONSTS(53) =  1.00000*CONSTS(27) ** (CONSTS(2) - 297.000)/10.0000
      CONSTS(54) =  CONSTS(45)*0.0100000
      CONSTS(55) = 1.90580/CONSTS(46)
      CONSTS(56) = 4.78030/CONSTS(48)
      CONSTS(57) = 763.756/CONSTS(49)
      CONSTS(58) =  4.00000*CONSTS(51)
      CONSTS(59) =  3.00000*CONSTS(51)
      CONSTS(60) =  2.00000*CONSTS(51)
      CONSTS(61) =  1.00000*CONSTS(51)
      CONSTS(62) =  4.00000*CONSTS(52)
      CONSTS(63) =  3.00000*CONSTS(52)
      CONSTS(64) =  2.00000*CONSTS(52)
      CONSTS(65) =  1.00000*CONSTS(52)
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+60.5000)/5.30000))
      RATES(17) = (ALGBRC(3) - STATES(17))/CONSTS(55)
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+43.0000)/17.3600))
      RATES(19) = (ALGBRC(5) - STATES(19))/CONSTS(56)
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1) - 44.9000)/12.0096))
      RATES(20) = (ALGBRC(6) - STATES(20))/CONSTS(57)
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.5000)/4.00000))
      ALGBRC(11) = ( 0.381170*(8.60000+ 14.7000*EXP(( - (STATES(1)+50.0000)*(STATES(1)+50.0000))/900.000)))/CONSTS(47)
      RATES(18) = (ALGBRC(4) - STATES(18))/ALGBRC(11)
      ALGBRC(25) =  CONSTS(53)*0.554320*EXP(- 0.0990740+ 0.0364410*STATES(1))
      ALGBRC(29) =  CONSTS(53)*0.000525480*EXP(- 0.0691020+ 0.00319450*STATES(1))
      ALGBRC(57) = STATES(33)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(38)
      RATES(36) = ( - ALGBRC(29)*STATES(36))/ALGBRC(57)+( ALGBRC(25)*STATES(35))/ALGBRC(57)
      ALGBRC(21) =  CONSTS(53)*3.15660*EXP(0.363520+ 0.0771930*STATES(1))
      ALGBRC(33) =  CONSTS(53)*1.44960*EXP(- 0.156600+ 0.0583530*STATES(1))
      RATES(35) = ( - (ALGBRC(33)+ALGBRC(25))*STATES(35))/ALGBRC(57)+( ALGBRC(21)*STATES(34))/ALGBRC(57)+( ALGBRC(29)*STATES(36))/ALGBRC(57)
      ALGBRC(17) =  CONSTS(53)*2.39100*EXP(- 13.3350 -  0.252890*STATES(1))
      ALGBRC(37) =  CONSTS(53)*1.53290*EXP(0.00931930+ 0.0410750*STATES(1))
      ALGBRC(49) =  CONSTS(53)*1.90460*EXP(- 2.48400+ 0.0204060*STATES(1))
      ALGBRC(53) =  CONSTS(53)*0.000216880*EXP(- 0.0634380+ 0.00466830*STATES(1))
      RATES(34) = ( - (ALGBRC(21)+ALGBRC(37)+ALGBRC(53))*STATES(34))/ALGBRC(57)+( ALGBRC(17)*STATES(33))/ALGBRC(57)+( ALGBRC(33)*STATES(35))/ALGBRC(57)+( ALGBRC(49)*STATES(37))/ALGBRC(57)
      ALGBRC(13) =  CONSTS(53)*0.120520*EXP(- 9.60280+ 0.0830250*STATES(1))
      ALGBRC(41) =  CONSTS(53)*1.61640*EXP(0.307630+ 0.00605350*STATES(1))
      RATES(33) = ( - (ALGBRC(17)+ALGBRC(41))*STATES(33))/ALGBRC(57)+( ALGBRC(37)*STATES(34))/ALGBRC(57)+( ALGBRC(13)*STATES(37))/ALGBRC(57)
      ALGBRC(8) =  CONSTS(53)*0.00392390*EXP(2.67930+ 0.00614680*STATES(1))
      ALGBRC(45) =  CONSTS(53)*0.0277350*EXP(0.0514900 -  0.0468650*STATES(1))
      RATES(38) = ( - ALGBRC(8)*STATES(38))/ALGBRC(57)+( ALGBRC(45)*STATES(37))/ALGBRC(57)
      RATES(37) = ( - (ALGBRC(13)+ALGBRC(45)+ALGBRC(49))*STATES(37))/ALGBRC(57)+( ALGBRC(8)*STATES(38))/ALGBRC(57)+( ALGBRC(53)*STATES(34))/ALGBRC(57)+( ALGBRC(41)*STATES(33))/ALGBRC(57)
      ALGBRC(18) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(15)/STATES(2))
      ALGBRC(22) =  CONSTS(14)*STATES(3)*(STATES(1) - ALGBRC(18))
      ALGBRC(26) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(15)/STATES(2))
      ALGBRC(30) =  CONSTS(17)*STATES(17)*STATES(18)*(STATES(1) - ALGBRC(26))
      ALGBRC(58) = ( CONSTS(28)*( EXP(( CONSTS(32)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(32) ** 3.00000*CONSTS(15) -  2.50000*EXP(( (CONSTS(32) - 1.00000)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*CONSTS(26) ** 3.00000*STATES(2)))/( (CONSTS(30) ** 3.00000+CONSTS(26) ** 3.00000)*(CONSTS(29)+CONSTS(15))*(1.00000+ CONSTS(31)*EXP(( (CONSTS(32) - 1.00000)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))
      RATES(2) = (( - 1.00000e-15*((ALGBRC(22)+ALGBRC(30)) -  ALGBRC(58)*2.00000))/( 2.00000*CONSTS(44)*CONSTS(3)))/(1.00000+( CONSTS(42)*CONSTS(38)*CONSTS(40)*STATES(2) ** CONSTS(42) - 1.00000)/STATES(2) ** CONSTS(42)+CONSTS(40) ** 2.00000+( CONSTS(43)*CONSTS(39)*CONSTS(41)*STATES(2) ** CONSTS(43) - 1.00000)/STATES(2) ** CONSTS(43)+CONSTS(41) ** 2.00000)
      ALGBRC(50) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(26)/STATES(32))
      ALGBRC(54) =  CONSTS(25)*STATES(33)*(STATES(1) - ALGBRC(50))
      ALGBRC(61) = ( CONSTS(33)*CONSTS(21)*STATES(32))/( (CONSTS(21)+CONSTS(34))*(CONSTS(35)+STATES(32))*(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0353000*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(64) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(26)/STATES(32))
      ALGBRC(67) =  CONSTS(36)*(STATES(1) - ALGBRC(64))
      RATES(32) = ( (- 1.00000/1.00000e+15)*(ALGBRC(54)+ALGBRC(67)+ ALGBRC(58)*3.00000+ ALGBRC(61)*3.00000))/( CONSTS(44)*CONSTS(3))
      ALGBRC(1) = TERNRY(VOI.GE. CONSTS(5)*0.00000.AND.VOI.LE. CONSTS(5)*1.00000,  CONSTS(5)*0.00000, TERNRY(VOI.GE. CONSTS(5)*1.00000.AND.VOI.LE. CONSTS(5)*2.00000,  CONSTS(5)*1.00000, TERNRY(VOI.GE. CONSTS(5)*2.00000.AND.VOI.LT. CONSTS(5)*3.00000,  CONSTS(5)*2.00000, TERNRY(VOI.GE. CONSTS(5)*3.00000.AND.VOI.LT. CONSTS(5)*4.00000,  CONSTS(5)*3.00000, 0.00000)
      ALGBRC(9) = VOI - ALGBRC(1)
      ALGBRC(14) = TERNRY(ALGBRC(9).LT.CONSTS(7),  CONSTS(6)*(STATES(1) - (CONSTS(12)+( CONSTS(13)*ALGBRC(9))/CONSTS(11))),  CONSTS(6)*(STATES(1) - (CONSTS(12)+( CONSTS(13)*(1.00000+EXP(- CONSTS(10)/( 2.00000*CONSTS(9))))*1.00000)/(1.00000+EXP(((ALGBRC(9) - CONSTS(11)) -  0.500000*CONSTS(10))/CONSTS(9))))))
      ALGBRC(34) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(38) =  CONSTS(20)*STATES(19)*STATES(20)*(STATES(1) - ALGBRC(34))
      ALGBRC(42) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(46) =  CONSTS(24)*STATES(22)*(STATES(1) - ALGBRC(42))
      ALGBRC(70) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(73) =  CONSTS(37)*(STATES(1) - ALGBRC(70))
      RATES(1) =  (- 1.00000/CONSTS(4))*(ALGBRC(54)+ALGBRC(22)+ALGBRC(30)+ALGBRC(38)+ALGBRC(46)+ALGBRC(58)+ALGBRC(61)+ALGBRC(73)+ALGBRC(67)+ALGBRC(14))
      RATES(21) = ( - 1.00000e-15*((ALGBRC(38)+ALGBRC(46)+ALGBRC(73)+ALGBRC(14)) -  ALGBRC(61)*2.00000))/( CONSTS(44)*CONSTS(3))
      ALGBRC(7) =  1.00000*EXP(( 8.47188*STATES(1))/( 1.00000*CONSTS(2)))
      ALGBRC(32) =  0.0443240*ALGBRC(7)
      ALGBRC(12) =  1.00000*EXP(( - 7.77556*STATES(1))/( 1.00000*CONSTS(2)))
      ALGBRC(52) =  0.000106000*ALGBRC(12)
      ALGBRC(66) =  1.00000*CONSTS(50)*STATES(2)
      ALGBRC(79) = STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)+STATES(31)+STATES(22)
      RATES(27) = ( - (CONSTS(58)+ALGBRC(32))*STATES(27))/ALGBRC(79)+( ALGBRC(66)*STATES(26))/ALGBRC(79)+( ALGBRC(52)*STATES(22))/ALGBRC(79)
      ALGBRC(28) =  0.000781000*ALGBRC(7)
      ALGBRC(48) =  0.0953120*ALGBRC(12)
      ALGBRC(63) =  2.00000*CONSTS(50)*STATES(2)
      RATES(26) = ( - (CONSTS(59)+ALGBRC(28)+ALGBRC(66))*STATES(26))/ALGBRC(79)+( ALGBRC(63)*STATES(25))/ALGBRC(79)+( ALGBRC(48)*STATES(31))/ALGBRC(79)+( CONSTS(58)*STATES(27))/ALGBRC(79)
      ALGBRC(24) =  2.81000e-05*ALGBRC(7)
      ALGBRC(44) =  32.6594*ALGBRC(12)
      ALGBRC(60) =  3.00000*CONSTS(50)*STATES(2)
      RATES(25) = ( - (CONSTS(60)+ALGBRC(24)+ALGBRC(63))*STATES(25))/ALGBRC(79)+( ALGBRC(60)*STATES(24))/ALGBRC(79)+( ALGBRC(44)*STATES(30))/ALGBRC(79)+( CONSTS(59)*STATES(26))/ALGBRC(79)
      ALGBRC(20) =  0.000869000*ALGBRC(7)
      ALGBRC(40) =  144.174*ALGBRC(12)
      ALGBRC(56) =  4.00000*CONSTS(50)*STATES(2)
      RATES(24) = ( - (CONSTS(61)+ALGBRC(20)+ALGBRC(60))*STATES(24))/ALGBRC(79)+( ALGBRC(56)*STATES(23))/ALGBRC(79)+( ALGBRC(40)*STATES(29))/ALGBRC(79)+( CONSTS(60)*STATES(25))/ALGBRC(79)
      ALGBRC(16) =  0.0216200*ALGBRC(7)
      ALGBRC(36) =  318.108*ALGBRC(12)
      RATES(23) = ( - (ALGBRC(56)+ALGBRC(16))*STATES(23))/ALGBRC(79)+( CONSTS(61)*STATES(24))/ALGBRC(79)+( ALGBRC(36)*STATES(28))/ALGBRC(79)
      ALGBRC(77) =  1.00000*CONSTS(50)*STATES(2)
      RATES(22) = ( - (CONSTS(62)+ALGBRC(52))*STATES(22))/ALGBRC(79)+( ALGBRC(77)*STATES(31))/ALGBRC(79)+( ALGBRC(32)*STATES(27))/ALGBRC(79)
      ALGBRC(75) =  2.00000*CONSTS(50)*STATES(2)
      RATES(31) = ( - (CONSTS(63)+ALGBRC(48)+ALGBRC(77))*STATES(31))/ALGBRC(79)+( ALGBRC(75)*STATES(30))/ALGBRC(79)+( ALGBRC(28)*STATES(26))/ALGBRC(79)+( CONSTS(62)*STATES(22))/ALGBRC(79)
      ALGBRC(72) =  3.00000*CONSTS(50)*STATES(2)
      RATES(30) = ( - (CONSTS(64)+ALGBRC(44)+ALGBRC(75))*STATES(30))/ALGBRC(79)+( ALGBRC(72)*STATES(29))/ALGBRC(79)+( ALGBRC(24)*STATES(25))/ALGBRC(79)+( CONSTS(63)*STATES(31))/ALGBRC(79)
      ALGBRC(69) =  4.00000*CONSTS(50)*STATES(2)
      RATES(29) = ( - (CONSTS(65)+ALGBRC(40)+ALGBRC(72))*STATES(29))/ALGBRC(79)+( ALGBRC(69)*STATES(28))/ALGBRC(79)+( ALGBRC(20)*STATES(24))/ALGBRC(79)+( CONSTS(64)*STATES(30))/ALGBRC(79)
      RATES(28) = ( - (ALGBRC(69)+ALGBRC(36))*STATES(28))/ALGBRC(79)+( CONSTS(65)*STATES(29))/ALGBRC(79)+( ALGBRC(16)*STATES(23))/ALGBRC(79)
      ALGBRC(2) =  CONSTS(45)*0.731000*EXP(STATES(1)/30.0000)
      ALGBRC(23) =  2.00000*ALGBRC(2)
      ALGBRC(27) =  1.00000*ALGBRC(2)
      ALGBRC(10) =  CONSTS(45)*0.214900*EXP(- STATES(1)/40.0000)
      ALGBRC(39) =  3.00000*ALGBRC(10)
      ALGBRC(43) =  4.00000*ALGBRC(10)
      ALGBRC(78) = ( CONSTS(45)*4.00000)/(1.00000+1.00000/STATES(2))
      ALGBRC(51) =  CONSTS(45)*0.0595600*EXP(- STATES(1)/40.0000)
      ALGBRC(47) =  CONSTS(45)*0.474200*EXP(STATES(1)/10.0000)
      ALGBRC(65) =  CONSTS(45)*0.00232000*EXP(- STATES(1)/280.000)
      ALGBRC(59) =  CONSTS(45)*0.0121300*EXP(STATES(1)/500.000)
      ALGBRC(71) = ( ALGBRC(43)*ALGBRC(59)*ALGBRC(51))/( ALGBRC(27)*ALGBRC(65))
      ALGBRC(62) =  CONSTS(45)*0.0219700*EXP(STATES(1)/500.000)
      ALGBRC(55) =  CONSTS(45)*0.0140700*EXP(- STATES(1)/300.000)
      ALGBRC(68) = ( ALGBRC(43)*ALGBRC(55)*ALGBRC(47))/( ALGBRC(27)*ALGBRC(62))
      ALGBRC(80) = STATES(7)+STATES(6)+STATES(5)+STATES(4)+STATES(11)+STATES(10)+STATES(9)+STATES(8)+STATES(3)+STATES(12)+STATES(13)+STATES(14)+STATES(16)+STATES(15)
      RATES(7) = (( ALGBRC(23)*STATES(6))/ALGBRC(80) - ( (ALGBRC(27)+ALGBRC(39)+ALGBRC(51)+ALGBRC(47)+ALGBRC(78))*STATES(7))/ALGBRC(80))+( ALGBRC(43)*STATES(3))/ALGBRC(80)+( ALGBRC(68)*STATES(14))/ALGBRC(80)+( ALGBRC(71)*STATES(13))/ALGBRC(80)+( CONSTS(54)*STATES(11))/ALGBRC(80)
      ALGBRC(19) =  3.00000*ALGBRC(2)
      ALGBRC(35) =  2.00000*ALGBRC(10)
      RATES(6) = (( ALGBRC(19)*STATES(5))/ALGBRC(80) - ( (ALGBRC(23)+ALGBRC(35)+ALGBRC(78))*STATES(6))/ALGBRC(80))+( ALGBRC(39)*STATES(7))/ALGBRC(80)+( CONSTS(54)*STATES(10))/ALGBRC(80)
      ALGBRC(15) =  4.00000*ALGBRC(2)
      ALGBRC(31) =  1.00000*ALGBRC(10)
      RATES(5) = (( ALGBRC(15)*STATES(4))/ALGBRC(80) - ( (ALGBRC(19)+ALGBRC(31)+ALGBRC(78))*STATES(5))/ALGBRC(80))+( ALGBRC(35)*STATES(6))/ALGBRC(80)+( CONSTS(54)*STATES(9))/ALGBRC(80)
      RATES(4) = ( - (ALGBRC(15)+ALGBRC(78))*STATES(4))/ALGBRC(80)+( ALGBRC(31)*STATES(5))/ALGBRC(80)+( CONSTS(54)*STATES(8))/ALGBRC(80)
      RATES(11) = ((( ALGBRC(78)*STATES(7))/ALGBRC(80)+( ALGBRC(23)*STATES(10))/ALGBRC(80)) - ( (ALGBRC(27)+ALGBRC(39)+ALGBRC(47)+ALGBRC(51)+CONSTS(54))*STATES(11))/ALGBRC(80))+( ALGBRC(43)*STATES(12))/ALGBRC(80)+( ALGBRC(68)*STATES(15))/ALGBRC(80)+( ALGBRC(71)*STATES(16))/ALGBRC(80)
      RATES(10) = ((( ALGBRC(78)*STATES(6))/ALGBRC(80)+( ALGBRC(19)*STATES(9))/ALGBRC(80)) - ( (ALGBRC(23)+ALGBRC(35)+CONSTS(54))*STATES(10))/ALGBRC(80))+( ALGBRC(39)*STATES(11))/ALGBRC(80)
      RATES(9) = ((( ALGBRC(78)*STATES(5))/ALGBRC(80)+( ALGBRC(15)*STATES(8))/ALGBRC(80)) - ( (ALGBRC(19)+ALGBRC(31)+CONSTS(54))*STATES(9))/ALGBRC(80))+( ALGBRC(35)*STATES(10))/ALGBRC(80)
      RATES(8) = (( ALGBRC(78)*STATES(4))/ALGBRC(80) - ( (ALGBRC(15)+CONSTS(54))*STATES(8))/ALGBRC(80))+( ALGBRC(31)*STATES(9))/ALGBRC(80)
      RATES(3) = (( ALGBRC(27)*STATES(7))/ALGBRC(80) - ( (ALGBRC(43)+ALGBRC(62)+ALGBRC(65)+ALGBRC(78))*STATES(3))/ALGBRC(80))+( ALGBRC(55)*STATES(14))/ALGBRC(80)+( ALGBRC(59)*STATES(13))/ALGBRC(80)+( CONSTS(54)*STATES(12))/ALGBRC(80)
      RATES(12) = ((( ALGBRC(78)*STATES(3))/ALGBRC(80)+( ALGBRC(27)*STATES(11))/ALGBRC(80)) - ( (ALGBRC(43)+ALGBRC(62)+ALGBRC(65)+CONSTS(54))*STATES(12))/ALGBRC(80))+( ALGBRC(55)*STATES(15))/ALGBRC(80)+( ALGBRC(59)*STATES(16))/ALGBRC(80)
      ALGBRC(76) = ALGBRC(65)
      ALGBRC(74) = ( ALGBRC(59)*ALGBRC(62))/ALGBRC(55)
      RATES(13) = ((( ALGBRC(51)*STATES(7))/ALGBRC(80)+( ALGBRC(65)*STATES(3))/ALGBRC(80)+( ALGBRC(76)*STATES(14))/ALGBRC(80)) - ( (ALGBRC(71)+ALGBRC(59)+ALGBRC(74)+ALGBRC(78))*STATES(13))/ALGBRC(80))+( CONSTS(54)*STATES(16))/ALGBRC(80)
      RATES(14) = ((( ALGBRC(47)*STATES(7))/ALGBRC(80)+( ALGBRC(62)*STATES(3))/ALGBRC(80)) - ( (ALGBRC(68)+ALGBRC(55)+ALGBRC(76)+ALGBRC(78))*STATES(14))/ALGBRC(80))+( ALGBRC(74)*STATES(13))/ALGBRC(80)+( CONSTS(54)*STATES(15))/ALGBRC(80)
      RATES(16) = (( ALGBRC(78)*STATES(13))/ALGBRC(80)+( ALGBRC(51)*STATES(11))/ALGBRC(80)+( ALGBRC(65)*STATES(12))/ALGBRC(80)+( ALGBRC(76)*STATES(15))/ALGBRC(80)) - ( (ALGBRC(71)+ALGBRC(59)+ALGBRC(74)+CONSTS(54))*STATES(16))/ALGBRC(80)
      RATES(15) = ((( ALGBRC(78)*STATES(14))/ALGBRC(80)+( ALGBRC(47)*STATES(11))/ALGBRC(80)+( ALGBRC(62)*STATES(12))/ALGBRC(80)) - ( (ALGBRC(68)+ALGBRC(55)+ALGBRC(76)+CONSTS(54))*STATES(15))/ALGBRC(80))+( ALGBRC(74)*STATES(16))/ALGBRC(80)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/(1.00000+EXP(- (STATES(1)+60.5000)/5.30000))
      ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+43.0000)/17.3600))
      ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1) - 44.9000)/12.0096))
      ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+75.5000)/4.00000))
      ALGBRC(11) = ( 0.381170*(8.60000+ 14.7000*EXP(( - (STATES(1)+50.0000)*(STATES(1)+50.0000))/900.000)))/CONSTS(47)
      ALGBRC(25) =  CONSTS(53)*0.554320*EXP(- 0.0990740+ 0.0364410*STATES(1))
      ALGBRC(29) =  CONSTS(53)*0.000525480*EXP(- 0.0691020+ 0.00319450*STATES(1))
      ALGBRC(57) = STATES(33)+STATES(34)+STATES(35)+STATES(36)+STATES(37)+STATES(38)
      ALGBRC(21) =  CONSTS(53)*3.15660*EXP(0.363520+ 0.0771930*STATES(1))
      ALGBRC(33) =  CONSTS(53)*1.44960*EXP(- 0.156600+ 0.0583530*STATES(1))
      ALGBRC(17) =  CONSTS(53)*2.39100*EXP(- 13.3350 -  0.252890*STATES(1))
      ALGBRC(37) =  CONSTS(53)*1.53290*EXP(0.00931930+ 0.0410750*STATES(1))
      ALGBRC(49) =  CONSTS(53)*1.90460*EXP(- 2.48400+ 0.0204060*STATES(1))
      ALGBRC(53) =  CONSTS(53)*0.000216880*EXP(- 0.0634380+ 0.00466830*STATES(1))
      ALGBRC(13) =  CONSTS(53)*0.120520*EXP(- 9.60280+ 0.0830250*STATES(1))
      ALGBRC(41) =  CONSTS(53)*1.61640*EXP(0.307630+ 0.00605350*STATES(1))
      ALGBRC(8) =  CONSTS(53)*0.00392390*EXP(2.67930+ 0.00614680*STATES(1))
      ALGBRC(45) =  CONSTS(53)*0.0277350*EXP(0.0514900 -  0.0468650*STATES(1))
      ALGBRC(18) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(15)/STATES(2))
      ALGBRC(22) =  CONSTS(14)*STATES(3)*(STATES(1) - ALGBRC(18))
      ALGBRC(26) =  (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(15)/STATES(2))
      ALGBRC(30) =  CONSTS(17)*STATES(17)*STATES(18)*(STATES(1) - ALGBRC(26))
      ALGBRC(58) = ( CONSTS(28)*( EXP(( CONSTS(32)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(32) ** 3.00000*CONSTS(15) -  2.50000*EXP(( (CONSTS(32) - 1.00000)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*CONSTS(26) ** 3.00000*STATES(2)))/( (CONSTS(30) ** 3.00000+CONSTS(26) ** 3.00000)*(CONSTS(29)+CONSTS(15))*(1.00000+ CONSTS(31)*EXP(( (CONSTS(32) - 1.00000)*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(50) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(26)/STATES(32))
      ALGBRC(54) =  CONSTS(25)*STATES(33)*(STATES(1) - ALGBRC(50))
      ALGBRC(61) = ( CONSTS(33)*CONSTS(21)*STATES(32))/( (CONSTS(21)+CONSTS(34))*(CONSTS(35)+STATES(32))*(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0353000*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(64) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(26)/STATES(32))
      ALGBRC(67) =  CONSTS(36)*(STATES(1) - ALGBRC(64))
      ALGBRC(1) = TERNRY(VOI.GE. CONSTS(5)*0.00000.AND.VOI.LE. CONSTS(5)*1.00000,  CONSTS(5)*0.00000, TERNRY(VOI.GE. CONSTS(5)*1.00000.AND.VOI.LE. CONSTS(5)*2.00000,  CONSTS(5)*1.00000, TERNRY(VOI.GE. CONSTS(5)*2.00000.AND.VOI.LT. CONSTS(5)*3.00000,  CONSTS(5)*2.00000, TERNRY(VOI.GE. CONSTS(5)*3.00000.AND.VOI.LT. CONSTS(5)*4.00000,  CONSTS(5)*3.00000, 0.00000)
      ALGBRC(9) = VOI - ALGBRC(1)
      ALGBRC(14) = TERNRY(ALGBRC(9).LT.CONSTS(7),  CONSTS(6)*(STATES(1) - (CONSTS(12)+( CONSTS(13)*ALGBRC(9))/CONSTS(11))),  CONSTS(6)*(STATES(1) - (CONSTS(12)+( CONSTS(13)*(1.00000+EXP(- CONSTS(10)/( 2.00000*CONSTS(9))))*1.00000)/(1.00000+EXP(((ALGBRC(9) - CONSTS(11)) -  0.500000*CONSTS(10))/CONSTS(9))))))
      ALGBRC(34) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(38) =  CONSTS(20)*STATES(19)*STATES(20)*(STATES(1) - ALGBRC(34))
      ALGBRC(42) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(46) =  CONSTS(24)*STATES(22)*(STATES(1) - ALGBRC(42))
      ALGBRC(70) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(21)/STATES(21))
      ALGBRC(73) =  CONSTS(37)*(STATES(1) - ALGBRC(70))
      ALGBRC(7) =  1.00000*EXP(( 8.47188*STATES(1))/( 1.00000*CONSTS(2)))
      ALGBRC(32) =  0.0443240*ALGBRC(7)
      ALGBRC(12) =  1.00000*EXP(( - 7.77556*STATES(1))/( 1.00000*CONSTS(2)))
      ALGBRC(52) =  0.000106000*ALGBRC(12)
      ALGBRC(66) =  1.00000*CONSTS(50)*STATES(2)
      ALGBRC(79) = STATES(23)+STATES(24)+STATES(25)+STATES(26)+STATES(27)+STATES(28)+STATES(29)+STATES(30)+STATES(31)+STATES(22)
      ALGBRC(28) =  0.000781000*ALGBRC(7)
      ALGBRC(48) =  0.0953120*ALGBRC(12)
      ALGBRC(63) =  2.00000*CONSTS(50)*STATES(2)
      ALGBRC(24) =  2.81000e-05*ALGBRC(7)
      ALGBRC(44) =  32.6594*ALGBRC(12)
      ALGBRC(60) =  3.00000*CONSTS(50)*STATES(2)
      ALGBRC(20) =  0.000869000*ALGBRC(7)
      ALGBRC(40) =  144.174*ALGBRC(12)
      ALGBRC(56) =  4.00000*CONSTS(50)*STATES(2)
      ALGBRC(16) =  0.0216200*ALGBRC(7)
      ALGBRC(36) =  318.108*ALGBRC(12)
      ALGBRC(77) =  1.00000*CONSTS(50)*STATES(2)
      ALGBRC(75) =  2.00000*CONSTS(50)*STATES(2)
      ALGBRC(72) =  3.00000*CONSTS(50)*STATES(2)
      ALGBRC(69) =  4.00000*CONSTS(50)*STATES(2)
      ALGBRC(2) =  CONSTS(45)*0.731000*EXP(STATES(1)/30.0000)
      ALGBRC(23) =  2.00000*ALGBRC(2)
      ALGBRC(27) =  1.00000*ALGBRC(2)
      ALGBRC(10) =  CONSTS(45)*0.214900*EXP(- STATES(1)/40.0000)
      ALGBRC(39) =  3.00000*ALGBRC(10)
      ALGBRC(43) =  4.00000*ALGBRC(10)
      ALGBRC(78) = ( CONSTS(45)*4.00000)/(1.00000+1.00000/STATES(2))
      ALGBRC(51) =  CONSTS(45)*0.0595600*EXP(- STATES(1)/40.0000)
      ALGBRC(47) =  CONSTS(45)*0.474200*EXP(STATES(1)/10.0000)
      ALGBRC(65) =  CONSTS(45)*0.00232000*EXP(- STATES(1)/280.000)
      ALGBRC(59) =  CONSTS(45)*0.0121300*EXP(STATES(1)/500.000)
      ALGBRC(71) = ( ALGBRC(43)*ALGBRC(59)*ALGBRC(51))/( ALGBRC(27)*ALGBRC(65))
      ALGBRC(62) =  CONSTS(45)*0.0219700*EXP(STATES(1)/500.000)
      ALGBRC(55) =  CONSTS(45)*0.0140700*EXP(- STATES(1)/300.000)
      ALGBRC(68) = ( ALGBRC(43)*ALGBRC(55)*ALGBRC(47))/( ALGBRC(27)*ALGBRC(62))
      ALGBRC(80) = STATES(7)+STATES(6)+STATES(5)+STATES(4)+STATES(11)+STATES(10)+STATES(9)+STATES(8)+STATES(3)+STATES(12)+STATES(13)+STATES(14)+STATES(16)+STATES(15)
      ALGBRC(19) =  3.00000*ALGBRC(2)
      ALGBRC(35) =  2.00000*ALGBRC(10)
      ALGBRC(15) =  4.00000*ALGBRC(2)
      ALGBRC(31) =  1.00000*ALGBRC(10)
      ALGBRC(76) = ALGBRC(65)
      ALGBRC(74) = ( ALGBRC(59)*ALGBRC(62))/ALGBRC(55)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END