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 106 entries in the algebraic variable array.
C There are a total of 22 entries in each of the rate and state variable arrays.
C There are a total of 93 entries in the constant variable array.
C
C
C VOI is time in component environment (ms).
C STATES(1) is V in component membrane (mV).
C CONSTS(1) is R in component membrane (mJ_per_mole_K).
C CONSTS(2) is T in component membrane (kelvin).
C CONSTS(3) is F in component membrane (C_per_mole).
C CONSTS(4) is Cm in component membrane (uF).
C ALGBRC(1) is I_Stim in component membrane (uA).
C CONSTS(5) is stim_period in component membrane (ms).
C CONSTS(6) is stim_duration in component membrane (ms).
C CONSTS(7) is stim_amplitude in component membrane (uA).
C ALGBRC(26) is i_Na in component sodium_current (nanoA).
C ALGBRC(28) is i_t in component Ca_independent_transient_outward_K_current (nanoA).
C ALGBRC(29) is i_ss in component steady_state_outward_K_current (nanoA).
C ALGBRC(33) is i_f in component hyperpolarisation_activated_current (nanoA).
C ALGBRC(30) is i_K1 in component inward_rectifier (nanoA).
C ALGBRC(34) is i_B_Na in component background_currents (nanoA).
C ALGBRC(36) is i_B_K in component background_currents (nanoA).
C ALGBRC(22) is i_NaK in component sodium_potassium_pump (uA).
C ALGBRC(106) is I_LCC in component convert_hinch (uA).
C ALGBRC(39) is I_NaCa in component convert_hinch (uA).
C ALGBRC(42) is I_pCa in component convert_hinch (uA).
C ALGBRC(45) is I_CaB in component convert_hinch (uA).
C ALGBRC(12) is FVRT in component membrane (dimensionless).
C ALGBRC(21) is FVRT_Ca in component membrane (dimensionless).
C CONSTS(8) is V_myo in component cell_geometry (um3).
C CONSTS(9) is V_SR in component cell_geometry (um3).
C CONSTS(10) is V_myo_uL in component cell_geometry (uL).
C CONSTS(11) is V_SR_uL in component cell_geometry (uL).
C ALGBRC(105) is I_LCC in component LCC_current (mM_per_ms).
C ALGBRC(38) is I_NaCa in component Na_Ca_Exchanger (mM_per_ms).
C ALGBRC(41) is I_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
C ALGBRC(44) is I_CaB in component Background_Ca_current (mM_per_ms).
C ALGBRC(101) is I_RyR in component RyR_current (mM_per_ms).
C ALGBRC(102) is I_RyR in component convert_hinch (mM_per_ms).
C CONSTS(12) is i_NaK_max in component sodium_potassium_pump (uA).
C CONSTS(13) is K_m_K in component sodium_potassium_pump (mM).
C CONSTS(14) is K_m_Na in component sodium_potassium_pump (mM).
C CONSTS(83) is sigma in component sodium_potassium_pump (dimensionless).
C CONSTS(15) is K_o in component standard_ionic_concentrations (millimolar).
C CONSTS(16) is Na_o in component standard_ionic_concentrations (millimolar).
C STATES(2) is Na_i in component intracellular_ion_concentrations (mM).
C STATES(3) is Ca_i in component intracellular_ion_concentrations (mM).
C STATES(4) is Ca_SR in component intracellular_ion_concentrations (mM).
C STATES(5) is K_i in component intracellular_ion_concentrations (mM).
C STATES(6) is TRPN in component intracellular_ion_concentrations (mM).
C ALGBRC(31) is i_f_Na in component hyperpolarisation_activated_current (nanoA).
C ALGBRC(32) is i_f_K in component hyperpolarisation_activated_current (nanoA).
C ALGBRC(40) is I_SERCA in component SERCA (mM_per_ms).
C ALGBRC(46) is I_SR in component SR_Ca_leak_current (mM_per_ms).
C ALGBRC(24) is Ca_b in component intracellular_ion_concentrations (mM).
C CONSTS(17) is Ca_TRPN_Max in component troponin (mM).
C ALGBRC(62) is J_TRPN in component troponin (mM_per_ms).
C ALGBRC(47) is beta_CMDN in component calmodulin_Ca_buffer (dimensionless).
C ALGBRC(25) is E_Na in component sodium_current (millivolt).
C CONSTS(18) is g_Na in component sodium_current (microS).
C CONSTS(84) is g_Na_endo in component sodium_current (microS).
C STATES(7) is m in component sodium_current_m_gate (dimensionless).
C STATES(8) is h in component sodium_current_h_gate (dimensionless).
C STATES(9) is j in component sodium_current_j_gate (dimensionless).
C ALGBRC(27) is E_K in component Ca_independent_transient_outward_K_current (millivolt).
C CONSTS(19) is g_t in component Ca_independent_transient_outward_K_current (microS).
C CONSTS(85) is g_t_endo in component Ca_independent_transient_outward_K_current (microS).
C CONSTS(20) is a_endo in component Ca_independent_transient_outward_K_current (dimensionless).
C CONSTS(21) is b_endo in component Ca_independent_transient_outward_K_current (dimensionless).
C STATES(10) is r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C STATES(11) is s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C STATES(12) is s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C CONSTS(22) is g_ss in component steady_state_outward_K_current (microS).
C STATES(13) is r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C STATES(14) is s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C CONSTS(23) is g_K1 in component inward_rectifier (microS).
C CONSTS(24) is g_f in component hyperpolarisation_activated_current (microS).
C CONSTS(25) is f_Na in component hyperpolarisation_activated_current (dimensionless).
C CONSTS(86) is f_K in component hyperpolarisation_activated_current (dimensionless).
C STATES(15) is y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(37) is i_B in component background_currents (nanoA).
C ALGBRC(35) is i_B_Ca in component background_currents (nanoA).
C CONSTS(26) is g_B_Na in component background_currents (microS).
C CONSTS(27) is g_B_Ca in component background_currents (microS).
C CONSTS(28) is g_B_K in component background_currents (microS).
C CONSTS(29) is E_Ca in component background_currents (millivolt).
C CONSTS(30) is Ca_o in component standard_ionic_concentrations (millimolar).
C CONSTS(31) is g_D in component CaRU (um3_per_ms).
C CONSTS(32) is J_R in component CaRU (um3_per_ms).
C CONSTS(33) is J_L in component CaRU (um3_per_ms).
C CONSTS(34) is N in component CaRU (dimensionless).
C CONSTS(35) is K_mNa in component Na_Ca_Exchanger (mM).
C CONSTS(36) is K_mCa in component Na_Ca_Exchanger (mM).
C CONSTS(37) is eta in component Na_Ca_Exchanger (dimensionless).
C CONSTS(38) is k_sat in component Na_Ca_Exchanger (dimensionless).
C CONSTS(39) is g_NCX in component Na_Ca_Exchanger (mM_per_ms).
C CONSTS(40) is g_SERCA in component SERCA (mM_per_ms).
C CONSTS(41) is K_SERCA in component SERCA (mM).
C CONSTS(42) is g_pCa in component Sarcolemmal_Ca_pump (mM_per_ms).
C CONSTS(43) is K_mpCa in component Sarcolemmal_Ca_pump (mM).
C ALGBRC(43) is E_Ca in component Background_Ca_current (mV).
C CONSTS(44) is g_CaB in component Background_Ca_current (mM_per_mV_ms).
C CONSTS(45) is g_SRl in component SR_Ca_leak_current (per_ms).
C CONSTS(46) is k_CMDN in component calmodulin_Ca_buffer (mM).
C CONSTS(47) is B_CMDN in component calmodulin_Ca_buffer (mM).
C STATES(16) is z in component tropomyosin (dimensionless).
C ALGBRC(55) is z_max in component tropomyosin (dimensionless).
C CONSTS(48) is k_on in component troponin (per_mM_per_ms).
C CONSTS(49) is k_Ref_off in component troponin (per_ms).
C CONSTS(50) is gamma_trpn in component troponin (dimensionless).
C CONSTS(51) is alpha_0 in component tropomyosin (per_ms).
C CONSTS(52) is alpha_r1 in component tropomyosin (per_ms).
C CONSTS(53) is alpha_r2 in component tropomyosin (per_ms).
C CONSTS(54) is n_Rel in component tropomyosin (dimensionless).
C CONSTS(55) is K_z in component tropomyosin (dimensionless).
C CONSTS(56) is n_Hill in component tropomyosin (dimensionless).
C CONSTS(57) is Ca_50ref in component tropomyosin (mM).
C CONSTS(58) is z_p in component tropomyosin (dimensionless).
C CONSTS(59) is beta_1 in component tropomyosin (dimensionless).
C ALGBRC(52) is Ca_50 in component tropomyosin (mM).
C ALGBRC(53) is Ca_TRPN_50 in component tropomyosin (mM).
C CONSTS(87) is K_2 in component tropomyosin (per_ms).
C CONSTS(91) is K_1 in component tropomyosin (per_ms).
C ALGBRC(54) is alpha_Tm in component tropomyosin (per_ms).
C ALGBRC(48) is beta_Tm in component tropomyosin (per_ms).
C CONSTS(60) is beta_0 in component filament_overlap (dimensionless).
C ALGBRC(50) is lambda in component Myofilaments (dimensionless).
C ALGBRC(61) is k_off in component troponin (per_ms).
C ALGBRC(60) is Tension in component Cross_Bridges (N_per_mm2).
C CONSTS(61) is T_ref in component length_independent_tension (N_per_mm2).
C ALGBRC(49) is ExtensionRatio in component Myofilaments (dimensionless).
C CONSTS(88) is dExtensionRatiodt in component Myofilaments (per_ms).
C ALGBRC(51) is lambda_prev in component Myofilaments (dimensionless).
C ALGBRC(56) is overlap in component filament_overlap (dimensionless).
C ALGBRC(57) is T_Base in component length_independent_tension (N_per_mm2).
C ALGBRC(58) is T_0 in component isometric_tension (N_per_mm2).
C ALGBRC(59) is Q in component Cross_Bridges (dimensionless).
C CONSTS(62) is a in component Cross_Bridges (dimensionless).
C STATES(17) is Q_1 in component Cross_Bridges (dimensionless).
C STATES(18) is Q_2 in component Cross_Bridges (dimensionless).
C STATES(19) is Q_3 in component Cross_Bridges (dimensionless).
C CONSTS(63) is A_1 in component Cross_Bridges (dimensionless).
C CONSTS(64) is A_2 in component Cross_Bridges (dimensionless).
C CONSTS(65) is A_3 in component Cross_Bridges (dimensionless).
C CONSTS(66) is alpha_1 in component Cross_Bridges (per_ms).
C CONSTS(67) is alpha_2 in component Cross_Bridges (per_ms).
C CONSTS(68) is alpha_3 in component Cross_Bridges (per_ms).
C ALGBRC(2) is m_infinity in component sodium_current_m_gate (dimensionless).
C ALGBRC(13) is tau_m in component sodium_current_m_gate (second).
C ALGBRC(3) is h_infinity in component sodium_current_h_gate (dimensionless).
C ALGBRC(14) is tau_h in component sodium_current_h_gate (second).
C ALGBRC(4) is j_infinity in component sodium_current_j_gate (dimensionless).
C ALGBRC(15) is tau_j in component sodium_current_j_gate (second).
C ALGBRC(16) is tau_r in component Ca_independent_transient_outward_K_current_r_gate (second).
C ALGBRC(5) is r_infinity in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C ALGBRC(17) is tau_s_endo in component Ca_independent_transient_outward_K_current_s_gate (second).
C ALGBRC(6) is s_infinity in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C ALGBRC(18) is tau_s_slow_endo in component Ca_independent_transient_outward_K_current_s_slow_gate (second).
C ALGBRC(7) is s_slow_infinity in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C ALGBRC(19) is tau_r_ss in component steady_state_outward_K_current_r_ss_gate (second).
C ALGBRC(8) is r_ss_infinity in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C CONSTS(89) is tau_s_ss in component steady_state_outward_K_current_s_ss_gate (second).
C ALGBRC(9) is s_ss_infinity in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C ALGBRC(20) is tau_y in component hyperpolarisation_activated_current_y_gate (second).
C ALGBRC(10) is y_infinity in component hyperpolarisation_activated_current_y_gate (dimensionless).
C ALGBRC(82) is C_oc in component DS_Calcium_Concentrations (mM).
C ALGBRC(80) is C_co in component DS_Calcium_Concentrations (mM).
C CONSTS(69) is V_L in component CaRU_Transitions (mV).
C CONSTS(70) is del_VL in component CaRU_Transitions (mV).
C CONSTS(71) is phi_L in component CaRU_Transitions (dimensionless).
C CONSTS(72) is t_L in component CaRU_Transitions (ms).
C CONSTS(73) is tau_L in component CaRU_Transitions (ms).
C CONSTS(90) is t_R in component CaRU_Transitions (ms).
C CONSTS(74) is tau_R in component CaRU_Transitions (ms).
C CONSTS(75) is phi_R in component CaRU_Transitions (dimensionless).
C CONSTS(76) is theta_R in component CaRU_Transitions (dimensionless).
C CONSTS(77) is K_RyR in component CaRU_Transitions (mM).
C CONSTS(78) is K_L in component CaRU_Transitions (mM).
C CONSTS(79) is a in component CaRU_Transitions (dimensionless).
C CONSTS(80) is b in component CaRU_Transitions (dimensionless).
C CONSTS(81) is c in component CaRU_Transitions (dimensionless).
C CONSTS(82) is d in component CaRU_Transitions (dimensionless).
C ALGBRC(63) is expVL in component CaRU_Transitions (dimensionless).
C ALGBRC(64) is alpha_p in component CaRU_Transitions (per_ms).
C CONSTS(92) is alpha_m in component CaRU_Transitions (per_ms).
C ALGBRC(83) is beta_poc in component CaRU_Transitions (per_ms).
C ALGBRC(65) is beta_pcc in component CaRU_Transitions (per_ms).
C CONSTS(93) is beta_m in component CaRU_Transitions (per_ms).
C ALGBRC(81) is epsilon_pco in component CaRU_Transitions (per_ms).
C ALGBRC(66) is epsilon_pcc in component CaRU_Transitions (per_ms).
C ALGBRC(67) is epsilon_m in component CaRU_Transitions (per_ms).
C ALGBRC(84) is mu_poc in component CaRU_Transitions (per_ms).
C ALGBRC(68) is mu_pcc in component CaRU_Transitions (per_ms).
C ALGBRC(85) is mu_moc in component CaRU_Transitions (per_ms).
C ALGBRC(69) is mu_mcc in component CaRU_Transitions (per_ms).
C ALGBRC(86) is denom in component CaRU_states (per_ms3).
C ALGBRC(87) is y_oc in component CaRU_states (dimensionless).
C ALGBRC(88) is y_co in component CaRU_states (dimensionless).
C ALGBRC(89) is y_oo in component CaRU_states (dimensionless).
C ALGBRC(91) is y_cc in component CaRU_states (dimensionless).
C ALGBRC(70) is y_ci in component CaRU_states (dimensionless).
C ALGBRC(71) is y_oi in component CaRU_states (dimensionless).
C ALGBRC(72) is y_ic in component CaRU_states (dimensionless).
C ALGBRC(73) is y_io in component CaRU_states (dimensionless).
C ALGBRC(93) is y_ii in component CaRU_states (dimensionless).
C ALGBRC(94) is r_1 in component CaRU_reduced_states (per_ms).
C ALGBRC(96) is r_2 in component CaRU_reduced_states (per_ms).
C ALGBRC(74) is r_3 in component CaRU_reduced_states (per_ms).
C ALGBRC(75) is r_4 in component CaRU_reduced_states (per_ms).
C ALGBRC(98) is r_5 in component CaRU_reduced_states (per_ms).
C ALGBRC(76) is r_6 in component CaRU_reduced_states (per_ms).
C ALGBRC(77) is r_7 in component CaRU_reduced_states (per_ms).
C ALGBRC(78) is r_8 in component CaRU_reduced_states (per_ms).
C STATES(20) is z_1 in component CaRU_reduced_states (dimensionless).
C STATES(21) is z_2 in component CaRU_reduced_states (dimensionless).
C STATES(22) is z_3 in component CaRU_reduced_states (dimensionless).
C ALGBRC(79) is z_4 in component CaRU_reduced_states (dimensionless).
C ALGBRC(11) is C_cc in component DS_Calcium_Concentrations (mM).
C ALGBRC(23) is C_oo in component DS_Calcium_Concentrations (mM).
C ALGBRC(97) is J_Loo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(95) is J_Loc in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(90) is J_Rco in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(92) is J_Roo in component LCC_and_RyR_fluxes (um3_mM_per_ms).
C ALGBRC(99) is J_R1 in component RyR_current (um3_mM_per_ms).
C ALGBRC(100) is J_R3 in component RyR_current (um3_mM_per_ms).
C ALGBRC(103) is J_L1 in component LCC_current (um3_mM_per_ms).
C ALGBRC(104) is J_L2 in component LCC_current (um3_mM_per_ms).
C RATES(1) is d/dt V in component membrane (mV).
C RATES(2) is d/dt Na_i in component intracellular_ion_concentrations (mM).
C RATES(5) is d/dt K_i in component intracellular_ion_concentrations (mM).
C RATES(6) is d/dt TRPN in component intracellular_ion_concentrations (mM).
C RATES(3) is d/dt Ca_i in component intracellular_ion_concentrations (mM).
C RATES(4) is d/dt Ca_SR in component intracellular_ion_concentrations (mM).
C RATES(16) is d/dt z in component tropomyosin (dimensionless).
C RATES(17) is d/dt Q_1 in component Cross_Bridges (dimensionless).
C RATES(18) is d/dt Q_2 in component Cross_Bridges (dimensionless).
C RATES(19) is d/dt Q_3 in component Cross_Bridges (dimensionless).
C RATES(7) is d/dt m in component sodium_current_m_gate (dimensionless).
C RATES(8) is d/dt h in component sodium_current_h_gate (dimensionless).
C RATES(9) is d/dt j in component sodium_current_j_gate (dimensionless).
C RATES(10) is d/dt r in component Ca_independent_transient_outward_K_current_r_gate (dimensionless).
C RATES(11) is d/dt s in component Ca_independent_transient_outward_K_current_s_gate (dimensionless).
C RATES(12) is d/dt s_slow in component Ca_independent_transient_outward_K_current_s_slow_gate (dimensionless).
C RATES(13) is d/dt r_ss in component steady_state_outward_K_current_r_ss_gate (dimensionless).
C RATES(14) is d/dt s_ss in component steady_state_outward_K_current_s_ss_gate (dimensionless).
C RATES(15) is d/dt y in component hyperpolarisation_activated_current_y_gate (dimensionless).
C RATES(20) is d/dt z_1 in component CaRU_reduced_states (dimensionless).
C RATES(21) is d/dt z_2 in component CaRU_reduced_states (dimensionless).
C RATES(22) is d/dt z_3 in component CaRU_reduced_states (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -80.50146
CONSTS(1) = 8314.5
CONSTS(2) = 295
CONSTS(3) = 96487
CONSTS(4) = 0.0001
CONSTS(5) = 1e3
CONSTS(6) = 10
CONSTS(7) = -0.6e-3
CONSTS(8) = 25.85e3
CONSTS(9) = 2.098e3
CONSTS(10) = 25.85e-6
CONSTS(11) = 2.098e-6
CONSTS(12) = 0.95e-4
CONSTS(13) = 1.5
CONSTS(14) = 10
CONSTS(15) = 5.4
CONSTS(16) = 140
STATES(2) = 10.73519
STATES(3) = 0.00007901351
STATES(4) = 700.0e-3
STATES(5) = 139.2751
STATES(6) = 0.067593139865
CONSTS(17) = 70e-3
CONSTS(18) = 0.8
STATES(7) = 0.004164108
STATES(8) = 0.6735613
STATES(9) = 0.6729362
CONSTS(19) = 0.035
CONSTS(20) = 0.583
CONSTS(21) = 0.417
STATES(10) = 0.002191519
STATES(11) = 0.9842542
STATES(12) = 0.6421196
CONSTS(22) = 0.007
STATES(13) = 0.002907171
STATES(14) = 0.3142767
CONSTS(23) = 0.024
CONSTS(24) = 0.00145
CONSTS(25) = 0.2
STATES(15) = 0.003578708
CONSTS(26) = 0.00008015
CONSTS(27) = 0.0000324
CONSTS(28) = 0.000138
CONSTS(29) = 65
CONSTS(30) = 1.2
CONSTS(31) = 0.065
CONSTS(32) = 0.02
CONSTS(33) = 9.13e-4
CONSTS(34) = 50000
CONSTS(35) = 87.5
CONSTS(36) = 1.38
CONSTS(37) = 0.35
CONSTS(38) = 0.1
CONSTS(39) = 38.5e-3
CONSTS(40) = 0.45e-3
CONSTS(41) = 0.5e-3
CONSTS(42) = 0.0035e-3
CONSTS(43) = 0.5e-3
CONSTS(44) = 2.6875e-8
CONSTS(45) = 1.8951e-5
CONSTS(46) = 2.382e-3
CONSTS(47) = 50e-3
STATES(16) = 0.014417937837
CONSTS(48) = 100
CONSTS(49) = 0.2
CONSTS(50) = 2
CONSTS(51) = 8e-3
CONSTS(52) = 2e-3
CONSTS(53) = 1.75e-3
CONSTS(54) = 3
CONSTS(55) = 0.15
CONSTS(56) = 3
CONSTS(57) = 1.05e-3
CONSTS(58) = 0.85
CONSTS(59) = -4
CONSTS(60) = 4.9
CONSTS(61) = 56.2
CONSTS(62) = 0.35
STATES(17) = 0
STATES(18) = 0
STATES(19) = 0
CONSTS(63) = -29
CONSTS(64) = 138
CONSTS(65) = 129
CONSTS(66) = 0.03
CONSTS(67) = 0.13
CONSTS(68) = 0.625
CONSTS(69) = -2
CONSTS(70) = 7
CONSTS(71) = 2.35
CONSTS(72) = 1
CONSTS(73) = 650
CONSTS(74) = 2.43
CONSTS(75) = 0.05
CONSTS(76) = 0.012
CONSTS(77) = 41e-3
CONSTS(78) = 0.22e-3
CONSTS(79) = 0.0625
CONSTS(80) = 14
CONSTS(81) = 0.01
CONSTS(82) = 100
STATES(20) = 0.98859
STATES(21) = 0.0087302
STATES(22) = 0.0026566
CONSTS(83) = (EXP(CONSTS(16)/67.3000) - 1.00000)/7.00000
CONSTS(84) = 1.33000*CONSTS(18)
CONSTS(85) = 0.464700*CONSTS(19)
CONSTS(86) = 1.00000 - CONSTS(25)
CONSTS(87) = (( CONSTS(53)*CONSTS(58) ** CONSTS(54))/(CONSTS(58) ** CONSTS(54)+CONSTS(55) ** CONSTS(54)))*(1.00000 - ( CONSTS(54)*CONSTS(55) ** CONSTS(54))/(CONSTS(58) ** CONSTS(54)+CONSTS(55) ** CONSTS(54)))
CONSTS(88) = 0.00000
CONSTS(89) = 2.10000
CONSTS(90) = 1.17000*CONSTS(72)
CONSTS(91) = ( CONSTS(53)*CONSTS(58) ** CONSTS(54) - 1.00000*CONSTS(54)*CONSTS(55) ** CONSTS(54))/CONSTS(58) ** CONSTS(54)+CONSTS(55) ** CONSTS(54) ** 2.00000
CONSTS(92) = CONSTS(71)/CONSTS(72)
CONSTS(93) = CONSTS(75)/CONSTS(90)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(17) = CONSTS(63)*CONSTS(88) - CONSTS(66)*STATES(17)
RATES(18) = CONSTS(64)*CONSTS(88) - CONSTS(67)*STATES(18)
RATES(19) = CONSTS(65)*CONSTS(88) - CONSTS(68)*STATES(19)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000))
RATES(14) = ((ALGBRC(9) - STATES(14))/CONSTS(89))*0.00100000
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+45.0000)/- 6.50000))
ALGBRC(13) = 0.00136000/(( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))+ 0.0800000*EXP(- STATES(1)/11.0000))
RATES(7) = ((ALGBRC(2) - STATES(7))/ALGBRC(13))*0.00100000
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000))
ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 0.000453700*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 0.00349000/( 0.135000*EXP(- (STATES(1)+80.0000)/6.80000)+ 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1))))
RATES(8) = ((ALGBRC(3) - STATES(8))/ALGBRC(14))*0.00100000
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000))
ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 0.0116300*(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))/EXP( - 2.53500e-07*STATES(1)), 0.00349000/( ((STATES(1)+37.7800)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))*( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))+( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400)))))
RATES(9) = ((ALGBRC(4) - STATES(9))/ALGBRC(15))*0.00100000
ALGBRC(16) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200))
RATES(10) = ((ALGBRC(5) - STATES(10))/ALGBRC(16))*0.00100000
ALGBRC(17) = 0.550000*EXP(- (STATES(1)+70.0000)/25.0000 ** 2.00000)+0.0490000
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
RATES(11) = ((ALGBRC(6) - STATES(11))/ALGBRC(17))*0.00100000
ALGBRC(18) = 3.30000*EXP(( (- (STATES(1)+70.0000)/30.0000)*(STATES(1)+70.0000))/30.0000)+0.0490000
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
RATES(12) = ((ALGBRC(7) - STATES(12))/ALGBRC(18))*0.00100000
ALGBRC(19) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+11.5000)/- 11.8200))
RATES(13) = ((ALGBRC(8) - STATES(13))/ALGBRC(19))*0.00100000
ALGBRC(20) = 1.00000/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900))
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800))
RATES(15) = ((ALGBRC(10) - STATES(15))/ALGBRC(20))*0.00100000
ALGBRC(1) = TERNRY(VOI - INT(VOI/CONSTS(5))*CONSTS(5).GE.0.00000.AND.VOI - INT(VOI/CONSTS(5))*CONSTS(5).LE.CONSTS(6), CONSTS(7), 0.00000)
ALGBRC(27) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(15)/STATES(5))
ALGBRC(28) = CONSTS(85)*STATES(10)*( CONSTS(20)*STATES(11)+ CONSTS(21)*STATES(12))*(STATES(1) - ALGBRC(27))
ALGBRC(29) = CONSTS(22)*STATES(13)*STATES(14)*(STATES(1) - ALGBRC(27))
ALGBRC(30) = ( (48.0000/(EXP((STATES(1)+37.0000)/25.0000)+EXP((STATES(1)+37.0000)/- 25.0000))+10.0000)*0.00100000)/(1.00000+EXP((STATES(1) - (ALGBRC(27)+76.7700))/- 17.0000))+( CONSTS(23)*(STATES(1) - (ALGBRC(27)+1.73000)))/( (1.00000+EXP(( 1.61300*CONSTS(3)*(STATES(1) - (ALGBRC(27)+1.73000)))/( CONSTS(1)*CONSTS(2))))*(1.00000+EXP((CONSTS(15) - 0.998800)/- 0.124000)))
ALGBRC(36) = CONSTS(28)*(STATES(1) - ALGBRC(27))
ALGBRC(22) = ( (( (( CONSTS(12)*1.00000)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*CONSTS(83)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*CONSTS(15))/(CONSTS(15)+CONSTS(13)))*1.00000)/(1.00000+CONSTS(14)/STATES(2) ** 4.00000)
ALGBRC(32) = CONSTS(24)*STATES(15)*CONSTS(86)*(STATES(1) - ALGBRC(27))
RATES(5) = ( - (ALGBRC(1)+ALGBRC(29)*0.00100000+ALGBRC(36)*0.00100000+ALGBRC(28)*0.00100000+ALGBRC(30)*0.00100000+ALGBRC(32)*0.00100000+ ALGBRC(22)*- 2.00000)*1.00000)/( CONSTS(10)*CONSTS(3))
ALGBRC(25) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(16)/STATES(2))
ALGBRC(26) = CONSTS(84)*STATES(7) ** 3.00000*STATES(8)*STATES(9)*(STATES(1) - ALGBRC(25))
ALGBRC(34) = CONSTS(26)*(STATES(1) - ALGBRC(25))
ALGBRC(12) = ( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))
ALGBRC(38) = ( CONSTS(39)*( EXP( CONSTS(37)*ALGBRC(12))*STATES(2) ** 3.00000*CONSTS(30) - EXP( (CONSTS(37) - 1.00000)*ALGBRC(12))*CONSTS(16) ** 3.00000*STATES(3)))/( (CONSTS(16) ** 3.00000+CONSTS(35) ** 3.00000)*(CONSTS(30)+CONSTS(36))*(1.00000+ CONSTS(38)*EXP( (CONSTS(37) - 1.00000)*ALGBRC(12))))
ALGBRC(39) = ALGBRC(38)*CONSTS(10)*CONSTS(3)
ALGBRC(31) = CONSTS(24)*STATES(15)*CONSTS(25)*(STATES(1) - ALGBRC(25))
RATES(2) = ( - (ALGBRC(26)*0.00100000+ALGBRC(34)*0.00100000+ ALGBRC(39)*3.00000+ ALGBRC(22)*3.00000+ALGBRC(31)*0.00100000)*1.00000)/( CONSTS(10)*CONSTS(3))
ALGBRC(24) = CONSTS(17) - STATES(6)
ALGBRC(49) = TERNRY(VOI.GT.300000., 1.00000, 1.00000)
ALGBRC(50) = TERNRY(ALGBRC(49).GT.0.800000.AND.ALGBRC(49).LE.1.15000, ALGBRC(49), TERNRY(ALGBRC(49).GT.1.15000, 1.15000, 0.800000)
ALGBRC(52) = CONSTS(57)*(1.00000+ CONSTS(59)*(ALGBRC(50) - 1.00000))
ALGBRC(53) = ( ALGBRC(52)*CONSTS(17))/(ALGBRC(52)+ (CONSTS(49)/CONSTS(48))*(1.00000 - ( (1.00000+ CONSTS(60)*(ALGBRC(50) - 1.00000))*0.500000)/CONSTS(50)))
ALGBRC(54) = CONSTS(51)*ALGBRC(24)/ALGBRC(53) ** CONSTS(56)
ALGBRC(48) = CONSTS(52)+( CONSTS(53)*STATES(16) ** CONSTS(54) - 1.00000)/(STATES(16) ** CONSTS(54)+CONSTS(55) ** CONSTS(54))
RATES(16) = ALGBRC(54)*(1.00000 - STATES(16)) - ALGBRC(48)*STATES(16)
ALGBRC(56) = 1.00000+ CONSTS(60)*(ALGBRC(50) - 1.00000)
ALGBRC(55) = (CONSTS(51)/ALGBRC(53)/CONSTS(17) ** CONSTS(56) - CONSTS(87))/(CONSTS(52)+CONSTS(91)+CONSTS(51)/ALGBRC(53)/CONSTS(17) ** CONSTS(56))
ALGBRC(57) = ( CONSTS(61)*STATES(16))/ALGBRC(55)
ALGBRC(58) = ALGBRC(57)*ALGBRC(56)
ALGBRC(59) = STATES(17)+STATES(18)+STATES(19)
ALGBRC(60) = TERNRY(ALGBRC(59).LT.0.00000, ( ALGBRC(58)*( CONSTS(62)*ALGBRC(59)+1.00000))/(1.00000 - ALGBRC(59)), ( ALGBRC(58)*(1.00000+ (CONSTS(62)+2.00000)*ALGBRC(59)))/(1.00000+ALGBRC(59)))
ALGBRC(61) = TERNRY(1.00000 - ALGBRC(60)/( CONSTS(50)*CONSTS(61)).GT.0.100000, CONSTS(49)*(1.00000 - ALGBRC(60)/( CONSTS(50)*CONSTS(61))), CONSTS(49)*0.100000)
ALGBRC(62) = (CONSTS(17) - STATES(6))*ALGBRC(61) - STATES(3)*STATES(6)*CONSTS(48)
RATES(6) = ALGBRC(62)
ALGBRC(21) = 2.00000*ALGBRC(12)
ALGBRC(82) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-09, (STATES(3)+( (CONSTS(33)/CONSTS(31))*CONSTS(30)*ALGBRC(21)*EXP(- ALGBRC(21)))/(1.00000 - EXP(- ALGBRC(21))))/(1.00000+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), (STATES(3)+ (CONSTS(33)/CONSTS(31))*CONSTS(30))/(1.00000+CONSTS(33)/CONSTS(31)))
ALGBRC(84) = (ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000)/( CONSTS(74)*(ALGBRC(82) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(68) = (STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000)/( CONSTS(74)*(STATES(3) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(63) = EXP((STATES(1) - CONSTS(69))/CONSTS(70))
ALGBRC(64) = ALGBRC(63)/( CONSTS(72)*(ALGBRC(63)+1.00000))
ALGBRC(65) = STATES(3) ** 2.00000/( CONSTS(90)*(STATES(3) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(83) = ALGBRC(82) ** 2.00000/( CONSTS(90)*(ALGBRC(82) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(86) = (ALGBRC(64)+CONSTS(92))*( (CONSTS(92)+CONSTS(93)+ALGBRC(83))*(CONSTS(93)+ALGBRC(65))+ ALGBRC(64)*(CONSTS(93)+ALGBRC(83)))
ALGBRC(87) = ( ALGBRC(64)*CONSTS(93)*(ALGBRC(64)+CONSTS(92)+CONSTS(93)+ALGBRC(65)))/ALGBRC(86)
ALGBRC(91) = ( CONSTS(92)*CONSTS(93)*(CONSTS(92)+ALGBRC(64)+CONSTS(93)+ALGBRC(83)))/ALGBRC(86)
ALGBRC(94) = ALGBRC(87)*ALGBRC(84)+ ALGBRC(91)*ALGBRC(68)
ALGBRC(85) = ( CONSTS(76)*CONSTS(82)*(ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))/( CONSTS(74)*( CONSTS(82)*ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))
ALGBRC(69) = ( CONSTS(76)*CONSTS(82)*(STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))/( CONSTS(74)*( CONSTS(82)*STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))
ALGBRC(96) = ( ALGBRC(64)*ALGBRC(85)+ CONSTS(92)*ALGBRC(69))/(ALGBRC(64)+CONSTS(92))
ALGBRC(66) = ( STATES(3)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*CONSTS(78)*(ALGBRC(63)+1.00000))
ALGBRC(77) = ( CONSTS(92)*ALGBRC(66))/(ALGBRC(64)+CONSTS(92))
ALGBRC(67) = ( CONSTS(80)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*( CONSTS(80)*ALGBRC(63)+CONSTS(79)))
ALGBRC(78) = ALGBRC(67)
ALGBRC(79) = ((1.00000 - STATES(20)) - STATES(21)) - STATES(22)
RATES(21) = ( ALGBRC(94)*STATES(20) - (ALGBRC(96)+ALGBRC(77))*STATES(21))+ ALGBRC(78)*ALGBRC(79)
ALGBRC(80) = (STATES(3)+ (CONSTS(32)/CONSTS(31))*STATES(4))/(1.00000+CONSTS(32)/CONSTS(31))
ALGBRC(81) = ( ALGBRC(80)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*CONSTS(78)*(ALGBRC(63)+1.00000))
ALGBRC(88) = ( CONSTS(92)*( ALGBRC(65)*(CONSTS(92)+CONSTS(93)+ALGBRC(83))+ ALGBRC(83)*ALGBRC(64)))/ALGBRC(86)
ALGBRC(98) = ALGBRC(88)*ALGBRC(81)+ ALGBRC(91)*ALGBRC(66)
ALGBRC(76) = ALGBRC(67)
RATES(20) = - (ALGBRC(94)+ALGBRC(98))*STATES(20)+ ALGBRC(96)*STATES(21)+ ALGBRC(76)*STATES(22)
ALGBRC(74) = ( CONSTS(93)*ALGBRC(68))/(CONSTS(93)+ALGBRC(65))
ALGBRC(75) = ALGBRC(69)
RATES(22) = ( ALGBRC(98)*STATES(20) - (ALGBRC(76)+ALGBRC(74))*STATES(22))+ ALGBRC(75)*ALGBRC(79)
ALGBRC(89) = ( ALGBRC(64)*( ALGBRC(83)*(ALGBRC(64)+CONSTS(93)+ALGBRC(65))+ ALGBRC(65)*CONSTS(92)))/ALGBRC(86)
ALGBRC(90) = ( CONSTS(32)*(STATES(4) - STATES(3)))/(1.00000+CONSTS(32)/CONSTS(31))
ALGBRC(92) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( CONSTS(32)*((STATES(4) - STATES(3))+ (( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*(STATES(4) - CONSTS(30)*EXP(- ALGBRC(21)))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), ( CONSTS(32)*((STATES(4) - STATES(3))+ (( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(4) - CONSTS(30)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(99) = ALGBRC(89)*ALGBRC(92)+ ALGBRC(90)*ALGBRC(88)
ALGBRC(100) = ( ALGBRC(90)*ALGBRC(65))/(CONSTS(93)+ALGBRC(65))
ALGBRC(101) = ( ( STATES(20)*ALGBRC(99)+ STATES(22)*ALGBRC(100))*CONSTS(34))/CONSTS(8)
ALGBRC(102) = 1.50000*ALGBRC(101)
ALGBRC(40) = ( CONSTS(40)*STATES(3) ** 2.00000)/(CONSTS(41) ** 2.00000+STATES(3) ** 2.00000)
ALGBRC(46) = CONSTS(45)*(STATES(4) - STATES(3))
RATES(4) = (CONSTS(10)/CONSTS(11))*((- ALGBRC(102)+ALGBRC(40)) - ALGBRC(46))
ALGBRC(33) = ALGBRC(31)+ALGBRC(32)
ALGBRC(97) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( (( CONSTS(33)*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*(( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(3))+ (CONSTS(32)/CONSTS(31))*( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(4))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(ALGBRC(21)))), ( (( CONSTS(33)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(30)*EXP(- 1.00000e-05) - STATES(3))+ (CONSTS(32)/CONSTS(31))*( CONSTS(30)*EXP(- 1.00000e-05) - STATES(4))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(95) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( (( CONSTS(33)*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(3)))/(1.00000+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), ( (( CONSTS(33)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(30)*EXP(- 1.00000e-05) - STATES(3)))/(1.00000+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(103) = ALGBRC(97)*ALGBRC(89)+ ALGBRC(95)*ALGBRC(87)
ALGBRC(104) = ( ALGBRC(95)*ALGBRC(64))/(ALGBRC(64)+CONSTS(92))
ALGBRC(105) = ( ( STATES(20)*ALGBRC(103)+ STATES(21)*ALGBRC(104))*CONSTS(34))/CONSTS(8)
ALGBRC(106) = - 1.50000*ALGBRC(105)*2.00000*CONSTS(10)*CONSTS(3)
ALGBRC(41) = ( CONSTS(42)*STATES(3))/(CONSTS(43)+STATES(3))
ALGBRC(42) = ALGBRC(41)*2.00000*CONSTS(10)*CONSTS(3)
ALGBRC(43) = (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(30)/STATES(3))
ALGBRC(44) = CONSTS(44)*(ALGBRC(43) - STATES(1))
ALGBRC(45) = - ALGBRC(44)*2.00000*CONSTS(10)*CONSTS(3)
RATES(1) = - (ALGBRC(26)*0.00100000+ALGBRC(28)*0.00100000+ALGBRC(29)*0.00100000+ALGBRC(33)*0.00100000+ALGBRC(30)*0.00100000+ALGBRC(34)*0.00100000+ALGBRC(36)*0.00100000+ALGBRC(22)+ALGBRC(1)+ALGBRC(45)+ALGBRC(39)+ALGBRC(42)+ALGBRC(106))/CONSTS(4)
ALGBRC(47) = 1.00000+( CONSTS(46)*CONSTS(47))/CONSTS(46)+STATES(3) ** 2.00000 ** - 1.00000
RATES(3) = ALGBRC(47)*(((ALGBRC(102) - ALGBRC(40))+ALGBRC(46)+ALGBRC(62)) - ( - 2.00000*ALGBRC(39)+ALGBRC(42)+ALGBRC(45)+ALGBRC(106))/( 2.00000*CONSTS(10)*CONSTS(3)))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+87.5000)/10.3000))
ALGBRC(2) = 1.00000/(1.00000+EXP((STATES(1)+45.0000)/- 6.50000))
ALGBRC(13) = 0.00136000/(( 0.320000*(STATES(1)+47.1300))/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))+ 0.0800000*EXP(- STATES(1)/11.0000))
ALGBRC(3) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000))
ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 0.000453700*(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 0.00349000/( 0.135000*EXP(- (STATES(1)+80.0000)/6.80000)+ 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1))))
ALGBRC(4) = 1.00000/(1.00000+EXP((STATES(1)+76.1000)/6.07000))
ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 0.0116300*(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))/EXP( - 2.53500e-07*STATES(1)), 0.00349000/( ((STATES(1)+37.7800)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))*( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))+( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400)))))
ALGBRC(16) = 1.00000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(1)+10.6000)/- 11.4200))
ALGBRC(17) = 0.550000*EXP(- (STATES(1)+70.0000)/25.0000 ** 2.00000)+0.0490000
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
ALGBRC(18) = 3.30000*EXP(( (- (STATES(1)+70.0000)/30.0000)*(STATES(1)+70.0000))/30.0000)+0.0490000
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+45.3000)/6.88410))
ALGBRC(19) = 10.0000/( 45.1600*EXP( 0.0357700*(STATES(1)+50.0000))+ 98.9000*EXP( - 0.100000*(STATES(1)+38.0000)))
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+11.5000)/- 11.8200))
ALGBRC(20) = 1.00000/( 0.118850*EXP((STATES(1)+80.0000)/28.3700)+ 0.562300*EXP((STATES(1)+80.0000)/- 14.1900))
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+138.600)/10.4800))
ALGBRC(1) = TERNRY(VOI - INT(VOI/CONSTS(5))*CONSTS(5).GE.0.00000.AND.VOI - INT(VOI/CONSTS(5))*CONSTS(5).LE.CONSTS(6), CONSTS(7), 0.00000)
ALGBRC(27) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(15)/STATES(5))
ALGBRC(28) = CONSTS(85)*STATES(10)*( CONSTS(20)*STATES(11)+ CONSTS(21)*STATES(12))*(STATES(1) - ALGBRC(27))
ALGBRC(29) = CONSTS(22)*STATES(13)*STATES(14)*(STATES(1) - ALGBRC(27))
ALGBRC(30) = ( (48.0000/(EXP((STATES(1)+37.0000)/25.0000)+EXP((STATES(1)+37.0000)/- 25.0000))+10.0000)*0.00100000)/(1.00000+EXP((STATES(1) - (ALGBRC(27)+76.7700))/- 17.0000))+( CONSTS(23)*(STATES(1) - (ALGBRC(27)+1.73000)))/( (1.00000+EXP(( 1.61300*CONSTS(3)*(STATES(1) - (ALGBRC(27)+1.73000)))/( CONSTS(1)*CONSTS(2))))*(1.00000+EXP((CONSTS(15) - 0.998800)/- 0.124000)))
ALGBRC(36) = CONSTS(28)*(STATES(1) - ALGBRC(27))
ALGBRC(22) = ( (( (( CONSTS(12)*1.00000)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ 0.0365000*CONSTS(83)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*CONSTS(15))/(CONSTS(15)+CONSTS(13)))*1.00000)/(1.00000+CONSTS(14)/STATES(2) ** 4.00000)
ALGBRC(32) = CONSTS(24)*STATES(15)*CONSTS(86)*(STATES(1) - ALGBRC(27))
ALGBRC(25) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(16)/STATES(2))
ALGBRC(26) = CONSTS(84)*STATES(7) ** 3.00000*STATES(8)*STATES(9)*(STATES(1) - ALGBRC(25))
ALGBRC(34) = CONSTS(26)*(STATES(1) - ALGBRC(25))
ALGBRC(12) = ( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))
ALGBRC(38) = ( CONSTS(39)*( EXP( CONSTS(37)*ALGBRC(12))*STATES(2) ** 3.00000*CONSTS(30) - EXP( (CONSTS(37) - 1.00000)*ALGBRC(12))*CONSTS(16) ** 3.00000*STATES(3)))/( (CONSTS(16) ** 3.00000+CONSTS(35) ** 3.00000)*(CONSTS(30)+CONSTS(36))*(1.00000+ CONSTS(38)*EXP( (CONSTS(37) - 1.00000)*ALGBRC(12))))
ALGBRC(39) = ALGBRC(38)*CONSTS(10)*CONSTS(3)
ALGBRC(31) = CONSTS(24)*STATES(15)*CONSTS(25)*(STATES(1) - ALGBRC(25))
ALGBRC(24) = CONSTS(17) - STATES(6)
ALGBRC(49) = TERNRY(VOI.GT.300000., 1.00000, 1.00000)
ALGBRC(50) = TERNRY(ALGBRC(49).GT.0.800000.AND.ALGBRC(49).LE.1.15000, ALGBRC(49), TERNRY(ALGBRC(49).GT.1.15000, 1.15000, 0.800000)
ALGBRC(52) = CONSTS(57)*(1.00000+ CONSTS(59)*(ALGBRC(50) - 1.00000))
ALGBRC(53) = ( ALGBRC(52)*CONSTS(17))/(ALGBRC(52)+ (CONSTS(49)/CONSTS(48))*(1.00000 - ( (1.00000+ CONSTS(60)*(ALGBRC(50) - 1.00000))*0.500000)/CONSTS(50)))
ALGBRC(54) = CONSTS(51)*ALGBRC(24)/ALGBRC(53) ** CONSTS(56)
ALGBRC(48) = CONSTS(52)+( CONSTS(53)*STATES(16) ** CONSTS(54) - 1.00000)/(STATES(16) ** CONSTS(54)+CONSTS(55) ** CONSTS(54))
ALGBRC(56) = 1.00000+ CONSTS(60)*(ALGBRC(50) - 1.00000)
ALGBRC(55) = (CONSTS(51)/ALGBRC(53)/CONSTS(17) ** CONSTS(56) - CONSTS(87))/(CONSTS(52)+CONSTS(91)+CONSTS(51)/ALGBRC(53)/CONSTS(17) ** CONSTS(56))
ALGBRC(57) = ( CONSTS(61)*STATES(16))/ALGBRC(55)
ALGBRC(58) = ALGBRC(57)*ALGBRC(56)
ALGBRC(59) = STATES(17)+STATES(18)+STATES(19)
ALGBRC(60) = TERNRY(ALGBRC(59).LT.0.00000, ( ALGBRC(58)*( CONSTS(62)*ALGBRC(59)+1.00000))/(1.00000 - ALGBRC(59)), ( ALGBRC(58)*(1.00000+ (CONSTS(62)+2.00000)*ALGBRC(59)))/(1.00000+ALGBRC(59)))
ALGBRC(61) = TERNRY(1.00000 - ALGBRC(60)/( CONSTS(50)*CONSTS(61)).GT.0.100000, CONSTS(49)*(1.00000 - ALGBRC(60)/( CONSTS(50)*CONSTS(61))), CONSTS(49)*0.100000)
ALGBRC(62) = (CONSTS(17) - STATES(6))*ALGBRC(61) - STATES(3)*STATES(6)*CONSTS(48)
ALGBRC(21) = 2.00000*ALGBRC(12)
ALGBRC(82) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-09, (STATES(3)+( (CONSTS(33)/CONSTS(31))*CONSTS(30)*ALGBRC(21)*EXP(- ALGBRC(21)))/(1.00000 - EXP(- ALGBRC(21))))/(1.00000+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), (STATES(3)+ (CONSTS(33)/CONSTS(31))*CONSTS(30))/(1.00000+CONSTS(33)/CONSTS(31)))
ALGBRC(84) = (ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000)/( CONSTS(74)*(ALGBRC(82) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(68) = (STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000)/( CONSTS(74)*(STATES(3) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(63) = EXP((STATES(1) - CONSTS(69))/CONSTS(70))
ALGBRC(64) = ALGBRC(63)/( CONSTS(72)*(ALGBRC(63)+1.00000))
ALGBRC(65) = STATES(3) ** 2.00000/( CONSTS(90)*(STATES(3) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(83) = ALGBRC(82) ** 2.00000/( CONSTS(90)*(ALGBRC(82) ** 2.00000+CONSTS(77) ** 2.00000))
ALGBRC(86) = (ALGBRC(64)+CONSTS(92))*( (CONSTS(92)+CONSTS(93)+ALGBRC(83))*(CONSTS(93)+ALGBRC(65))+ ALGBRC(64)*(CONSTS(93)+ALGBRC(83)))
ALGBRC(87) = ( ALGBRC(64)*CONSTS(93)*(ALGBRC(64)+CONSTS(92)+CONSTS(93)+ALGBRC(65)))/ALGBRC(86)
ALGBRC(91) = ( CONSTS(92)*CONSTS(93)*(CONSTS(92)+ALGBRC(64)+CONSTS(93)+ALGBRC(83)))/ALGBRC(86)
ALGBRC(94) = ALGBRC(87)*ALGBRC(84)+ ALGBRC(91)*ALGBRC(68)
ALGBRC(85) = ( CONSTS(76)*CONSTS(82)*(ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))/( CONSTS(74)*( CONSTS(82)*ALGBRC(82) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))
ALGBRC(69) = ( CONSTS(76)*CONSTS(82)*(STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))/( CONSTS(74)*( CONSTS(82)*STATES(3) ** 2.00000+ CONSTS(81)*CONSTS(77) ** 2.00000))
ALGBRC(96) = ( ALGBRC(64)*ALGBRC(85)+ CONSTS(92)*ALGBRC(69))/(ALGBRC(64)+CONSTS(92))
ALGBRC(66) = ( STATES(3)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*CONSTS(78)*(ALGBRC(63)+1.00000))
ALGBRC(77) = ( CONSTS(92)*ALGBRC(66))/(ALGBRC(64)+CONSTS(92))
ALGBRC(67) = ( CONSTS(80)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*( CONSTS(80)*ALGBRC(63)+CONSTS(79)))
ALGBRC(78) = ALGBRC(67)
ALGBRC(79) = ((1.00000 - STATES(20)) - STATES(21)) - STATES(22)
ALGBRC(80) = (STATES(3)+ (CONSTS(32)/CONSTS(31))*STATES(4))/(1.00000+CONSTS(32)/CONSTS(31))
ALGBRC(81) = ( ALGBRC(80)*(ALGBRC(63)+CONSTS(79)))/( CONSTS(73)*CONSTS(78)*(ALGBRC(63)+1.00000))
ALGBRC(88) = ( CONSTS(92)*( ALGBRC(65)*(CONSTS(92)+CONSTS(93)+ALGBRC(83))+ ALGBRC(83)*ALGBRC(64)))/ALGBRC(86)
ALGBRC(98) = ALGBRC(88)*ALGBRC(81)+ ALGBRC(91)*ALGBRC(66)
ALGBRC(76) = ALGBRC(67)
ALGBRC(74) = ( CONSTS(93)*ALGBRC(68))/(CONSTS(93)+ALGBRC(65))
ALGBRC(75) = ALGBRC(69)
ALGBRC(89) = ( ALGBRC(64)*( ALGBRC(83)*(ALGBRC(64)+CONSTS(93)+ALGBRC(65))+ ALGBRC(65)*CONSTS(92)))/ALGBRC(86)
ALGBRC(90) = ( CONSTS(32)*(STATES(4) - STATES(3)))/(1.00000+CONSTS(32)/CONSTS(31))
ALGBRC(92) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( CONSTS(32)*((STATES(4) - STATES(3))+ (( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*(STATES(4) - CONSTS(30)*EXP(- ALGBRC(21)))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), ( CONSTS(32)*((STATES(4) - STATES(3))+ (( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(STATES(4) - CONSTS(30)*EXP(- 1.00000e-05))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(99) = ALGBRC(89)*ALGBRC(92)+ ALGBRC(90)*ALGBRC(88)
ALGBRC(100) = ( ALGBRC(90)*ALGBRC(65))/(CONSTS(93)+ALGBRC(65))
ALGBRC(101) = ( ( STATES(20)*ALGBRC(99)+ STATES(22)*ALGBRC(100))*CONSTS(34))/CONSTS(8)
ALGBRC(102) = 1.50000*ALGBRC(101)
ALGBRC(40) = ( CONSTS(40)*STATES(3) ** 2.00000)/(CONSTS(41) ** 2.00000+STATES(3) ** 2.00000)
ALGBRC(46) = CONSTS(45)*(STATES(4) - STATES(3))
ALGBRC(33) = ALGBRC(31)+ALGBRC(32)
ALGBRC(97) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( (( CONSTS(33)*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*(( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(3))+ (CONSTS(32)/CONSTS(31))*( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(4))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(ALGBRC(21)))), ( (( CONSTS(33)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*(( CONSTS(30)*EXP(- 1.00000e-05) - STATES(3))+ (CONSTS(32)/CONSTS(31))*( CONSTS(30)*EXP(- 1.00000e-05) - STATES(4))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(95) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-05, ( (( CONSTS(33)*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21))))*( CONSTS(30)*EXP(- ALGBRC(21)) - STATES(3)))/(1.00000+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), ( (( CONSTS(33)*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05)))*( CONSTS(30)*EXP(- 1.00000e-05) - STATES(3)))/(1.00000+( (CONSTS(33)/CONSTS(31))*1.00000e-05)/(1.00000 - EXP(- 1.00000e-05))))
ALGBRC(103) = ALGBRC(97)*ALGBRC(89)+ ALGBRC(95)*ALGBRC(87)
ALGBRC(104) = ( ALGBRC(95)*ALGBRC(64))/(ALGBRC(64)+CONSTS(92))
ALGBRC(105) = ( ( STATES(20)*ALGBRC(103)+ STATES(21)*ALGBRC(104))*CONSTS(34))/CONSTS(8)
ALGBRC(106) = - 1.50000*ALGBRC(105)*2.00000*CONSTS(10)*CONSTS(3)
ALGBRC(41) = ( CONSTS(42)*STATES(3))/(CONSTS(43)+STATES(3))
ALGBRC(42) = ALGBRC(41)*2.00000*CONSTS(10)*CONSTS(3)
ALGBRC(43) = (( CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3)))*log(CONSTS(30)/STATES(3))
ALGBRC(44) = CONSTS(44)*(ALGBRC(43) - STATES(1))
ALGBRC(45) = - ALGBRC(44)*2.00000*CONSTS(10)*CONSTS(3)
ALGBRC(47) = 1.00000+( CONSTS(46)*CONSTS(47))/CONSTS(46)+STATES(3) ** 2.00000 ** - 1.00000
ALGBRC(11) = STATES(3)
ALGBRC(23) = TERNRY(ABS(ALGBRC(21)).GT.1.00000e-09, (STATES(3)+ (CONSTS(32)/CONSTS(31))*STATES(4)+( (CONSTS(33)/CONSTS(31))*CONSTS(30)*ALGBRC(21)*EXP(- ALGBRC(21)))/(1.00000 - EXP(- ALGBRC(21))))/(1.00000+CONSTS(32)/CONSTS(31)+( (CONSTS(33)/CONSTS(31))*ALGBRC(21))/(1.00000 - EXP(- ALGBRC(21)))), (STATES(3)+ (CONSTS(32)/CONSTS(31))*STATES(4)+ (CONSTS(33)/CONSTS(31))*CONSTS(30))/(1.00000+CONSTS(32)/CONSTS(31)+CONSTS(33)/CONSTS(31)))
ALGBRC(35) = CONSTS(27)*(STATES(1) - CONSTS(29))
ALGBRC(37) = ALGBRC(34)+ALGBRC(35)+ALGBRC(36)
ALGBRC(51) = ALGBRC(49)
ALGBRC(70) = CONSTS(92)/(ALGBRC(64)+CONSTS(92))
ALGBRC(71) = ALGBRC(64)/(ALGBRC(64)+CONSTS(92))
ALGBRC(72) = CONSTS(93)/(ALGBRC(65)+CONSTS(93))
ALGBRC(73) = ALGBRC(65)/(ALGBRC(65)+CONSTS(93))
ALGBRC(93) = (((((((1.00000 - ALGBRC(87)) - ALGBRC(88)) - ALGBRC(89)) - ALGBRC(91)) - ALGBRC(70)) - ALGBRC(72)) - ALGBRC(71)) - ALGBRC(73)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END