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 90 entries in the algebraic variable array.
C There are a total of 33 entries in each of the rate and state variable arrays.
C There are a total of 104 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is R in component Membrane (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component Membrane (kelvin).
C CONSTS(3) is F in component Membrane (coulomb_per_mole).
C CONSTS(4) is C in component Membrane (microF).
C CONSTS(86) is RTONF in component Membrane (millivolt).
C ALGBRC(52) is i_f in component i_f (nanoA).
C ALGBRC(53) is i_NaK in component i_NaK (nanoA).
C ALGBRC(67) is i_NaCa in component i_NaCa (nanoA).
C ALGBRC(69) is i_Na in component i_Na (nanoA).
C ALGBRC(77) is i_Kr in component i_Kr (nanoA).
C ALGBRC(81) is i_Ks in component i_Ks (nanoA).
C ALGBRC(75) is i_to in component i_to (nanoA).
C ALGBRC(73) is i_CaL in component i_CaL (nanoA).
C ALGBRC(74) is i_CaT in component i_CaT (nanoA).
C ALGBRC(84) is i_KACh in component i_KACh (nanoA).
C ALGBRC(8) is V in component Membrane (millivolt).
C CONSTS(5) is clamp_mode in component Membrane (dimensionless).
C ALGBRC(4) is V_clamp in component Voltage_clamp (millivolt).
C STATES(1) is V_ode in component Membrane (millivolt).
C ALGBRC(87) is i_tot in component Membrane (nanoA).
C CONSTS(6) is t_holding in component Voltage_clamp (second).
C CONSTS(7) is t_test in component Voltage_clamp (second).
C CONSTS(8) is V_test in component Voltage_clamp (millivolt).
C CONSTS(9) is V_holding in component Voltage_clamp (millivolt).
C CONSTS(10) is Iva_3_uM in component Rate_modulation_experiments (dimensionless).
C CONSTS(11) is Cs_5_mM in component Rate_modulation_experiments (dimensionless).
C CONSTS(12) is ACh in component Rate_modulation_experiments (millimolar).
C CONSTS(13) is Iso_1_uM in component Rate_modulation_experiments (dimensionless).
C CONSTS(14) is BAPTA_10_mM in component Rate_modulation_experiments (dimensionless).
C ALGBRC(17) is Nai in component Nai_concentration (millimolar).
C CONSTS(15) is Nao in component Ionic_values (millimolar).
C CONSTS(16) is Ki in component Ionic_values (millimolar).
C CONSTS(17) is Ko in component Ionic_values (millimolar).
C STATES(2) is Ca_sub in component Ca_dynamics (millimolar).
C CONSTS(18) is Cao in component Ionic_values (millimolar).
C ALGBRC(33) is E_Na in component Ionic_values (millivolt).
C CONSTS(88) is E_K in component Ionic_values (millivolt).
C ALGBRC(1) is E_Ca in component Ionic_values (millivolt).
C CONSTS(99) is V_sub in component Cell_parameters (millimetre3).
C CONSTS(101) is V_i in component Cell_parameters (millimetre3).
C ALGBRC(47) is i_fNa in component i_f (nanoA).
C ALGBRC(72) is i_siNa in component i_CaL (nanoA).
C STATES(3) is Nai_ in component Nai_concentration (millimolar).
C ALGBRC(43) is ICs_on_Icontrol in component i_f (dimensionless).
C ALGBRC(50) is i_fK in component i_f (nanoA).
C CONSTS(83) is g_f_Na in component i_f (microS).
C CONSTS(84) is g_f_K in component i_f (microS).
C CONSTS(19) is Km_f in component i_f (millimolar).
C STATES(4) is y in component i_f_y_gate (dimensionless).
C ALGBRC(9) is tau_y in component i_f_y_gate (second).
C ALGBRC(26) is y_infinity in component i_f_y_gate (dimensionless).
C CONSTS(90) is ACh_shift in component i_f_y_gate (millivolt).
C CONSTS(91) is Iso_shift in component i_f_y_gate (millivolt).
C CONSTS(20) is Km_Kp in component i_NaK (millimolar).
C CONSTS(21) is Km_Nap in component i_NaK (millimolar).
C CONSTS(22) is i_NaK_max in component i_NaK (nanoA).
C CONSTS(92) is Iso_increase in component i_NaK (dimensionless).
C CONSTS(23) is K_NaCa in component i_NaCa (nanoA).
C ALGBRC(64) is x1 in component i_NaCa (dimensionless).
C ALGBRC(60) is x2 in component i_NaCa (dimensionless).
C ALGBRC(65) is x3 in component i_NaCa (dimensionless).
C ALGBRC(66) is x4 in component i_NaCa (dimensionless).
C ALGBRC(55) is k41 in component i_NaCa (dimensionless).
C CONSTS(93) is k34 in component i_NaCa (dimensionless).
C ALGBRC(63) is k23 in component i_NaCa (dimensionless).
C ALGBRC(62) is k21 in component i_NaCa (dimensionless).
C ALGBRC(59) is k32 in component i_NaCa (dimensionless).
C ALGBRC(54) is k43 in component i_NaCa (dimensionless).
C ALGBRC(57) is k12 in component i_NaCa (dimensionless).
C ALGBRC(58) is k14 in component i_NaCa (dimensionless).
C CONSTS(24) is Qci in component i_NaCa (dimensionless).
C CONSTS(25) is Qn in component i_NaCa (dimensionless).
C CONSTS(26) is Qco in component i_NaCa (dimensionless).
C CONSTS(27) is K3ni in component i_NaCa (millimolar).
C CONSTS(28) is Kci in component i_NaCa (millimolar).
C CONSTS(29) is K1ni in component i_NaCa (millimolar).
C CONSTS(30) is K2ni in component i_NaCa (millimolar).
C CONSTS(31) is Kcni in component i_NaCa (millimolar).
C CONSTS(32) is K3no in component i_NaCa (millimolar).
C CONSTS(33) is K1no in component i_NaCa (millimolar).
C CONSTS(34) is K2no in component i_NaCa (millimolar).
C CONSTS(35) is Kco in component i_NaCa (millimolar).
C ALGBRC(61) is do in component i_NaCa (dimensionless).
C ALGBRC(56) is di in component i_NaCa (dimensionless).
C CONSTS(36) is g_Na in component i_Na (microS).
C ALGBRC(68) is E_mh in component i_Na (millivolt).
C STATES(5) is m in component i_Na_m_gate (dimensionless).
C STATES(6) is h in component i_Na_h_gate (dimensionless).
C ALGBRC(27) is alpha_m in component i_Na_m_gate (per_second).
C ALGBRC(41) is beta_m in component i_Na_m_gate (per_second).
C CONSTS(37) is delta_m in component i_Na_m_gate (millivolt).
C ALGBRC(10) is E0_m in component i_Na_m_gate (millivolt).
C ALGBRC(11) is alpha_h in component i_Na_h_gate (per_second).
C ALGBRC(28) is beta_h in component i_Na_h_gate (per_second).
C ALGBRC(70) is i_siCa in component i_CaL (nanoA).
C ALGBRC(71) is i_siK in component i_CaL (nanoA).
C CONSTS(95) is ACh_block in component i_CaL (dimensionless).
C CONSTS(38) is P_CaL in component i_CaL (nanoA_per_millimolar).
C STATES(7) is dL in component i_CaL_dL_gate (dimensionless).
C STATES(8) is fL in component i_CaL_fL_gate (dimensionless).
C STATES(9) is fCa in component i_CaL_fCa_gate (dimensionless).
C CONSTS(94) is Iso_increase in component i_CaL (dimensionless).
C ALGBRC(12) is dL_infinity in component i_CaL_dL_gate (dimensionless).
C ALGBRC(51) is tau_dL in component i_CaL_dL_gate (second).
C ALGBRC(42) is alpha_dL in component i_CaL_dL_gate (per_second).
C ALGBRC(49) is beta_dL in component i_CaL_dL_gate (per_second).
C ALGBRC(29) is adVm in component i_CaL_dL_gate (millivolt).
C ALGBRC(46) is bdVm in component i_CaL_dL_gate (millivolt).
C CONSTS(96) is Iso_shift in component i_CaL_dL_gate (millivolt).
C CONSTS(97) is Iso_slope in component i_CaL_dL_gate (dimensionless).
C ALGBRC(13) is fL_infinity in component i_CaL_fL_gate (dimensionless).
C ALGBRC(30) is tau_fL in component i_CaL_fL_gate (second).
C CONSTS(39) is alpha_fCa in component i_CaL_fCa_gate (per_second).
C ALGBRC(2) is fCa_infinity in component i_CaL_fCa_gate (dimensionless).
C ALGBRC(6) is tau_fCa in component i_CaL_fCa_gate (second).
C CONSTS(40) is Km_fCa in component i_CaL_fCa_gate (millimolar).
C CONSTS(41) is P_CaT in component i_CaT (nanoA_per_millimolar).
C STATES(10) is dT in component i_CaT_dT_gate (dimensionless).
C STATES(11) is fT in component i_CaT_fT_gate (dimensionless).
C ALGBRC(14) is dT_infinity in component i_CaT_dT_gate (dimensionless).
C ALGBRC(31) is tau_dT in component i_CaT_dT_gate (second).
C ALGBRC(15) is fT_infinity in component i_CaT_fT_gate (dimensionless).
C ALGBRC(32) is tau_fT in component i_CaT_fT_gate (second).
C ALGBRC(76) is j_SRCarel in component Ca_SR_release (millimolar_per_second).
C STATES(12) is R in component Ca_SR_release (dimensionless).
C STATES(13) is O in component Ca_SR_release (dimensionless).
C STATES(14) is I in component Ca_SR_release (dimensionless).
C STATES(15) is RI in component Ca_SR_release (dimensionless).
C CONSTS(42) is ks in component Ca_SR_release (per_second).
C CONSTS(43) is MaxSR in component Ca_SR_release (dimensionless).
C CONSTS(44) is MinSR in component Ca_SR_release (dimensionless).
C CONSTS(45) is EC50_SR in component Ca_SR_release (millimolar).
C CONSTS(46) is HSR in component Ca_SR_release (dimensionless).
C ALGBRC(7) is koSRCa in component Ca_SR_release (per_millimolar2_second).
C ALGBRC(16) is kiSRCa in component Ca_SR_release (per_millimolar_second).
C CONSTS(47) is koCa in component Ca_SR_release (per_millimolar2_second).
C CONSTS(48) is kiCa in component Ca_SR_release (per_millimolar_second).
C ALGBRC(3) is kCaSR in component Ca_SR_release (dimensionless).
C CONSTS(49) is kim in component Ca_SR_release (per_second).
C CONSTS(50) is kom in component Ca_SR_release (per_second).
C STATES(16) is Ca_jsr in component Ca_dynamics (millimolar).
C ALGBRC(78) is j_Ca_dif in component Ca_intracellular_fluxes (millimolar_per_second).
C ALGBRC(80) is j_up in component Ca_intracellular_fluxes (millimolar_per_second).
C ALGBRC(83) is j_tr in component Ca_intracellular_fluxes (millimolar_per_second).
C CONSTS(51) is tau_dif_Ca in component Ca_intracellular_fluxes (second).
C CONSTS(52) is tau_tr in component Ca_intracellular_fluxes (second).
C CONSTS(89) is P_up in component Ca_intracellular_fluxes (millimolar_per_second).
C CONSTS(53) is P_up_basal in component Ca_intracellular_fluxes (millimolar_per_second).
C CONSTS(85) is b_up in component Ca_intracellular_fluxes (dimensionless).
C CONSTS(54) is K_up in component Ca_intracellular_fluxes (millimolar).
C STATES(17) is Ca_nsr in component Ca_dynamics (millimolar).
C STATES(18) is Cai in component Ca_dynamics (millimolar).
C CONSTS(55) is TC_tot in component Ca_buffering (millimolar).
C CONSTS(56) is TMC_tot in component Ca_buffering (millimolar).
C CONSTS(57) is CM_tot in component Ca_buffering (millimolar).
C CONSTS(58) is CQ_tot in component Ca_buffering (millimolar).
C ALGBRC(82) is delta_fTC in component Ca_buffering (per_second).
C ALGBRC(85) is delta_fTMC in component Ca_buffering (per_second).
C ALGBRC(89) is delta_fCMs in component Ca_buffering (per_second).
C ALGBRC(88) is delta_fCMi in component Ca_buffering (per_second).
C ALGBRC(86) is delta_fCQ in component Ca_buffering (per_second).
C ALGBRC(5) is delta_fTMM in component Ca_buffering (per_second).
C STATES(19) is fTMM in component Ca_buffering (dimensionless).
C STATES(20) is fCMi in component Ca_buffering (dimensionless).
C STATES(21) is fCMs in component Ca_buffering (dimensionless).
C STATES(22) is fTC in component Ca_buffering (dimensionless).
C STATES(23) is fTMC in component Ca_buffering (dimensionless).
C STATES(24) is fCQ in component Ca_buffering (dimensionless).
C CONSTS(59) is kf_TC in component Ca_buffering (per_millimolar_second).
C CONSTS(60) is kf_TMM in component Ca_buffering (per_millimolar_second).
C CONSTS(61) is kf_TMC in component Ca_buffering (per_millimolar_second).
C CONSTS(62) is kf_CM in component Ca_buffering (per_millimolar_second).
C CONSTS(63) is kf_CQ in component Ca_buffering (per_millimolar_second).
C CONSTS(64) is kb_TC in component Ca_buffering (per_second).
C CONSTS(65) is kb_TMC in component Ca_buffering (per_second).
C CONSTS(66) is kb_TMM in component Ca_buffering (per_second).
C CONSTS(67) is kb_CM in component Ca_buffering (per_second).
C CONSTS(68) is kb_CQ in component Ca_buffering (per_second).
C CONSTS(69) is Mgi in component Ca_buffering (millimolar).
C CONSTS(100) is V_jsr in component Cell_parameters (millimetre3).
C CONSTS(102) is V_nsr in component Cell_parameters (millimetre3).
C ALGBRC(90) is BAPTA in component Ca_dynamics (millimolar).
C CONSTS(70) is kfBAPTA in component Ca_dynamics (per_millimolar_second).
C CONSTS(71) is kbBAPTA in component Ca_dynamics (per_second).
C STATES(25) is fBAPTA in component Ca_dynamics (millimolar).
C STATES(26) is fBAPTA_sub in component Ca_dynamics (millimolar).
C CONSTS(72) is T in component Ca_dynamics (second).
C CONSTS(98) is V_cell in component Cell_parameters (millimetre3).
C CONSTS(73) is V_jsr_part in component Cell_parameters (dimensionless).
C CONSTS(74) is V_i_part in component Cell_parameters (dimensionless).
C CONSTS(75) is V_nsr_part in component Cell_parameters (dimensionless).
C CONSTS(76) is R_cell in component Cell_parameters (micrometre).
C CONSTS(77) is L_cell in component Cell_parameters (micrometre).
C CONSTS(78) is L_sub in component Cell_parameters (micrometre).
C CONSTS(79) is g_to in component i_to (microS).
C STATES(27) is q in component i_to_q_gate (dimensionless).
C STATES(28) is r in component i_to_r_gate (dimensionless).
C ALGBRC(18) is q_infinity in component i_to_q_gate (dimensionless).
C ALGBRC(34) is tau_q in component i_to_q_gate (second).
C ALGBRC(19) is r_infinity in component i_to_r_gate (dimensionless).
C ALGBRC(35) is tau_r in component i_to_r_gate (second).
C CONSTS(80) is g_Kr in component i_Kr (microS).
C STATES(29) is paS in component i_Kr_pa_gate (dimensionless).
C STATES(30) is paF in component i_Kr_pa_gate (dimensionless).
C STATES(31) is piy in component i_Kr_pi_gate (dimensionless).
C ALGBRC(20) is pa_infinity in component i_Kr_pa_gate (dimensionless).
C ALGBRC(21) is alfapaF in component i_Kr_pa_gate (per_second).
C ALGBRC(22) is betapaF in component i_Kr_pa_gate (per_second).
C ALGBRC(36) is tau_paS in component i_Kr_pa_gate (second).
C ALGBRC(37) is tau_paF in component i_Kr_pa_gate (second).
C ALGBRC(38) is pi_infinity in component i_Kr_pi_gate (dimensionless).
C ALGBRC(23) is tau_pi in component i_Kr_pi_gate (second).
C CONSTS(87) is g_Ks in component i_Ks (microS).
C ALGBRC(79) is E_Ks in component i_Ks (millivolt).
C STATES(32) is n in component i_Ks_n_gate (dimensionless).
C ALGBRC(24) is n_infinity in component i_Ks_n_gate (dimensionless).
C ALGBRC(48) is tau_n in component i_Ks_n_gate (second).
C CONSTS(103) is Iso_shift in component i_Ks_n_gate (millivolt).
C CONSTS(81) is shift in component i_Ks_n_gate (millivolt).
C ALGBRC(39) is alpha_n in component i_Ks_n_gate (per_second).
C ALGBRC(44) is beta_n in component i_Ks_n_gate (per_second).
C CONSTS(82) is g_KACh in component i_KACh (microS).
C STATES(33) is a in component i_KACh_a_gate (dimensionless).
C CONSTS(104) is alpha_a in component i_KACh_a_gate (per_second).
C ALGBRC(25) is beta_a in component i_KACh_a_gate (per_second).
C ALGBRC(40) is a_infinity in component i_KACh_a_gate (dimensionless).
C ALGBRC(45) is tau_a in component i_KACh_a_gate (second).
C RATES(1) is d/dt V_ode in component Membrane (millivolt).
C RATES(3) is d/dt Nai_ in component Nai_concentration (millimolar).
C RATES(4) is d/dt y in component i_f_y_gate (dimensionless).
C RATES(5) is d/dt m in component i_Na_m_gate (dimensionless).
C RATES(6) is d/dt h in component i_Na_h_gate (dimensionless).
C RATES(7) is d/dt dL in component i_CaL_dL_gate (dimensionless).
C RATES(8) is d/dt fL in component i_CaL_fL_gate (dimensionless).
C RATES(9) is d/dt fCa in component i_CaL_fCa_gate (dimensionless).
C RATES(10) is d/dt dT in component i_CaT_dT_gate (dimensionless).
C RATES(11) is d/dt fT in component i_CaT_fT_gate (dimensionless).
C RATES(12) is d/dt R in component Ca_SR_release (dimensionless).
C RATES(13) is d/dt O in component Ca_SR_release (dimensionless).
C RATES(14) is d/dt I in component Ca_SR_release (dimensionless).
C RATES(15) is d/dt RI in component Ca_SR_release (dimensionless).
C RATES(22) is d/dt fTC in component Ca_buffering (dimensionless).
C RATES(23) is d/dt fTMC in component Ca_buffering (dimensionless).
C RATES(19) is d/dt fTMM in component Ca_buffering (dimensionless).
C RATES(20) is d/dt fCMi in component Ca_buffering (dimensionless).
C RATES(21) is d/dt fCMs in component Ca_buffering (dimensionless).
C RATES(24) is d/dt fCQ in component Ca_buffering (dimensionless).
C RATES(18) is d/dt Cai in component Ca_dynamics (millimolar).
C RATES(25) is d/dt fBAPTA in component Ca_dynamics (millimolar).
C RATES(2) is d/dt Ca_sub in component Ca_dynamics (millimolar).
C RATES(26) is d/dt fBAPTA_sub in component Ca_dynamics (millimolar).
C RATES(17) is d/dt Ca_nsr in component Ca_dynamics (millimolar).
C RATES(16) is d/dt Ca_jsr in component Ca_dynamics (millimolar).
C RATES(27) is d/dt q in component i_to_q_gate (dimensionless).
C RATES(28) is d/dt r in component i_to_r_gate (dimensionless).
C RATES(29) is d/dt paS in component i_Kr_pa_gate (dimensionless).
C RATES(30) is d/dt paF in component i_Kr_pa_gate (dimensionless).
C RATES(31) is d/dt piy in component i_Kr_pi_gate (dimensionless).
C RATES(32) is d/dt n in component i_Ks_n_gate (dimensionless).
C RATES(33) is d/dt a in component i_KACh_a_gate (dimensionless).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 8314.472
CONSTS(2) = 310
CONSTS(3) = 96485.3415
CONSTS(4) = 3.2e-5
CONSTS(5) = 0
STATES(1) = -52
CONSTS(6) = 0.5
CONSTS(7) = 0.5
CONSTS(8) = -35
CONSTS(9) = -45
CONSTS(10) = 0
CONSTS(11) = 0
CONSTS(12) = 0
CONSTS(13) = 0
CONSTS(14) = 0
CONSTS(15) = 140
CONSTS(16) = 140
CONSTS(17) = 5.4
STATES(2) = 1e-5
CONSTS(18) = 1.8
STATES(3) = 7.5
CONSTS(19) = 45
STATES(4) = 0.181334538702451
CONSTS(20) = 1.4
CONSTS(21) = 14
CONSTS(22) = 0.063
CONSTS(23) = 4
CONSTS(24) = 0.1369
CONSTS(25) = 0.4315
CONSTS(26) = 0
CONSTS(27) = 26.44
CONSTS(28) = 0.0207
CONSTS(29) = 395.3
CONSTS(30) = 2.289
CONSTS(31) = 26.44
CONSTS(32) = 4.663
CONSTS(33) = 1628
CONSTS(34) = 561.4
CONSTS(35) = 3.663
CONSTS(36) = 0.0125
STATES(5) = 0.440131579215766
STATES(6) = 1.3676940140066e-5
CONSTS(37) = 1e-5
CONSTS(38) = 0.2
STATES(7) = 0
STATES(8) = 0.497133507285601
STATES(9) = 0.697998543259722
CONSTS(39) = 0.01
CONSTS(40) = 0.00035
CONSTS(41) = 0.02
STATES(10) = 0
STATES(11) = 0
STATES(12) = 0.912317231017262
STATES(13) = 1.7340201253e-7
STATES(14) = 7.86181717518e-8
STATES(15) = 0.211148145512825
CONSTS(42) = 250000000
CONSTS(43) = 15
CONSTS(44) = 1
CONSTS(45) = 0.45
CONSTS(46) = 2.5
CONSTS(47) = 10000
CONSTS(48) = 500
CONSTS(49) = 5
CONSTS(50) = 60
STATES(16) = 0.316762674605
CONSTS(51) = 4e-5
CONSTS(52) = 0.04
CONSTS(53) = 12
CONSTS(54) = 0.0006
STATES(17) = 1.05386465080816
STATES(18) = 1e-5
CONSTS(55) = 0.031
CONSTS(56) = 0.062
CONSTS(57) = 0.045
CONSTS(58) = 10
STATES(19) = 0.501049376634
STATES(20) = 0.0373817991524254
STATES(21) = 0.054381370046
STATES(22) = 0.0180519400676086
STATES(23) = 0.281244308217086
STATES(24) = 0.299624275428735
CONSTS(59) = 88800
CONSTS(60) = 2277
CONSTS(61) = 227700
CONSTS(62) = 227700
CONSTS(63) = 534
CONSTS(64) = 446
CONSTS(65) = 7.51
CONSTS(66) = 751
CONSTS(67) = 542
CONSTS(68) = 445
CONSTS(69) = 2.5
CONSTS(70) = 940000
CONSTS(71) = 119.38
STATES(25) = 0
STATES(26) = 0
CONSTS(72) = 6.928
CONSTS(73) = 0.0012
CONSTS(74) = 0.46
CONSTS(75) = 0.0116
CONSTS(76) = 4
CONSTS(77) = 70
CONSTS(78) = 0.02
CONSTS(79) = 0.002
STATES(27) = 0.506139850982478
STATES(28) = 0.0144605370597924
CONSTS(80) = 0.0021637
STATES(29) = 0.322999177802891
STATES(30) = 0.0990510403258968
STATES(31) = 0.705410877258545
STATES(32) = 0
CONSTS(81) = 0
CONSTS(82) = 0.00864
STATES(33) = 0
CONSTS(83) = TERNRY(CONSTS(10).GE.1.00000, 0.0300000*(1.00000 - 0.660000), 0.0300000)
CONSTS(84) = TERNRY(CONSTS(10).GE.1.00000, 0.0300000*(1.00000 - 0.660000), 0.0300000)
CONSTS(85) = TERNRY(CONSTS(13).GT.0.00000, - 0.250000, TERNRY(CONSTS(12).GT.0.00000, ( 0.700000*CONSTS(12))/(9.00000e-05+CONSTS(12)), 0.00000)
CONSTS(86) = ( CONSTS(1)*CONSTS(2))/CONSTS(3)
CONSTS(87) = TERNRY(CONSTS(13).GT.0.00000, 1.20000*0.00165760, 0.00165760)
CONSTS(88) = CONSTS(86)*log(CONSTS(17)/CONSTS(16))
CONSTS(89) = CONSTS(53)*(1.00000 - CONSTS(85))
CONSTS(90) = TERNRY(CONSTS(12).GT.0.00000, - 1.00000 - ( 9.89800* 1.00000*CONSTS(12) ** 0.618000)/( 1.00000*CONSTS(12) ** 0.618000+0.00122423), 0.00000)
CONSTS(91) = TERNRY(CONSTS(13).GT.0.00000, 7.50000, 0.00000)
CONSTS(92) = TERNRY(CONSTS(13).GT.0.00000, 1.20000, 1.00000)
CONSTS(93) = CONSTS(15)/(CONSTS(32)+CONSTS(15))
CONSTS(94) = TERNRY(CONSTS(13).GT.0.00000, 1.23000, 1.00000)
CONSTS(95) = ( 0.310000*CONSTS(12))/(CONSTS(12)+9.00000e-05)
CONSTS(96) = TERNRY(CONSTS(13).GT.0.00000, - 8.00000, 0.00000)
CONSTS(97) = TERNRY(CONSTS(13).GT.0.00000, 0.690000, 1.00000)
CONSTS(98) = 1.00000e-09* 3.14159265358979*CONSTS(76) ** 2.00000*CONSTS(77)
CONSTS(99) = 1.00000e-09*2.00000* 3.14159265358979*CONSTS(78)*(CONSTS(76) - CONSTS(78)/2.00000)*CONSTS(77)
CONSTS(100) = CONSTS(73)*CONSTS(98)
CONSTS(101) = CONSTS(74)*CONSTS(98) - CONSTS(99)
CONSTS(102) = CONSTS(75)*CONSTS(98)
CONSTS(103) = TERNRY(CONSTS(13).GT.0.00000, - 14.0000, 0.00000)
CONSTS(104) = (3.59880 - 0.0256410)/(1.00000+1.21550e-06/ 1.00000*CONSTS(12) ** 1.69510)+0.0256410
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(23)+STATES(19))) - CONSTS(66)*STATES(19)
RATES(19) = ALGBRC(5)
ALGBRC(2) = CONSTS(40)/(CONSTS(40)+STATES(2))
ALGBRC(6) = ( 0.00100000*ALGBRC(2))/CONSTS(39)
RATES(9) = (ALGBRC(2) - STATES(9))/ALGBRC(6)
ALGBRC(3) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(16) ** CONSTS(46))
ALGBRC(7) = CONSTS(47)/ALGBRC(3)
ALGBRC(16) = CONSTS(48)*ALGBRC(3)
RATES(12) = ( CONSTS(49)*STATES(15) - ALGBRC(16)*STATES(2)*STATES(12)) - ( ALGBRC(7)*STATES(2) ** 2.00000*STATES(12) - CONSTS(50)*STATES(13))
RATES(13) = ( ALGBRC(7)*STATES(2) ** 2.00000*STATES(12) - CONSTS(50)*STATES(13)) - ( ALGBRC(16)*STATES(2)*STATES(13) - CONSTS(49)*STATES(14))
RATES(14) = ( ALGBRC(16)*STATES(2)*STATES(13) - CONSTS(49)*STATES(14)) - ( CONSTS(50)*STATES(14) - ALGBRC(7)*STATES(2) ** 2.00000*STATES(15))
RATES(15) = ( CONSTS(50)*STATES(14) - ALGBRC(7)*STATES(2) ** 2.00000*STATES(15)) - ( CONSTS(49)*STATES(15) - ALGBRC(16)*STATES(2)*STATES(12))
ALGBRC(4) = TERNRY(VOI.GT.CONSTS(6).AND.VOI.LT.CONSTS(6)+CONSTS(7), CONSTS(8), CONSTS(9))
ALGBRC(8) = TERNRY(CONSTS(5).GE.1.00000, ALGBRC(4), STATES(1))
ALGBRC(9) = 0.716653/( 0.0708000*EXP(- (((ALGBRC(8)+5.00000) - CONSTS(90)) - CONSTS(91))/20.2791)+ 10.6000*EXP(((ALGBRC(8) - CONSTS(90)) - CONSTS(91))/18.0000))
ALGBRC(26) = 1.00000/(1.00000+EXP((((ALGBRC(8)+52.5000) - CONSTS(90)) - CONSTS(91))/9.00000))
RATES(4) = (ALGBRC(26) - STATES(4))/ALGBRC(9)
ALGBRC(11) = 20.0000*EXP( - 0.125000*(ALGBRC(8)+75.0000))
ALGBRC(28) = 2000.00/( 320.000*EXP( - 0.100000*(ALGBRC(8)+75.0000))+1.00000)
RATES(6) = ALGBRC(11)*(1.00000 - STATES(6)) - ALGBRC(28)*STATES(6)
ALGBRC(13) = 1.00000/(1.00000+EXP((ALGBRC(8)+37.4000)/5.30000))
ALGBRC(30) = 0.00100000*(44.3000+ 230.000*EXP(- (ALGBRC(8)+36.0000)/10.0000 ** 2.00000))
RATES(8) = (ALGBRC(13) - STATES(8))/ALGBRC(30)
ALGBRC(14) = 1.00000/(1.00000+EXP(- (ALGBRC(8)+38.3000)/5.50000))
ALGBRC(31) = 0.00100000/( 1.06800*EXP((ALGBRC(8)+38.3000)/30.0000)+ 1.06800*EXP(- (ALGBRC(8)+38.3000)/30.0000))
RATES(10) = (ALGBRC(14) - STATES(10))/ALGBRC(31)
ALGBRC(15) = 1.00000/(1.00000+EXP((ALGBRC(8)+58.7000)/3.80000))
ALGBRC(32) = 1.00000/( 16.6700*EXP(- (ALGBRC(8)+75.0000)/83.3000)+ 16.6700*EXP((ALGBRC(8)+75.0000)/15.3800))
RATES(11) = (ALGBRC(15) - STATES(11))/ALGBRC(32)
ALGBRC(18) = 1.00000/(1.00000+EXP((ALGBRC(8)+49.0000)/13.0000))
ALGBRC(34) = 0.00100000*0.600000*(65.1700/( 0.570000*EXP( - 0.0800000*(ALGBRC(8)+44.0000))+ 0.0650000*EXP( 0.100000*(ALGBRC(8)+45.9300)))+10.1000)
RATES(27) = (ALGBRC(18) - STATES(27))/ALGBRC(34)
ALGBRC(19) = 1.00000/(1.00000+EXP(- (ALGBRC(8) - 19.3000)/15.0000))
ALGBRC(35) = 0.00100000*0.660000*1.40000*(15.5900/( 1.03700*EXP( 0.0900000*(ALGBRC(8)+30.6100))+ 0.369000*EXP( - 0.120000*(ALGBRC(8)+23.8400)))+2.98000)
RATES(28) = (ALGBRC(19) - STATES(28))/ALGBRC(35)
ALGBRC(20) = 1.00000/(1.00000+EXP(- (ALGBRC(8)+14.8000)/8.50000))
ALGBRC(36) = 0.846554/( 4.20000*EXP(ALGBRC(8)/17.0000)+ 0.150000*EXP(- ALGBRC(8)/21.6000))
RATES(29) = (ALGBRC(20) - STATES(29))/ALGBRC(36)
ALGBRC(37) = 1.00000/( 30.0000*EXP(ALGBRC(8)/10.0000)+EXP(- ALGBRC(8)/12.0000))
RATES(30) = (ALGBRC(20) - STATES(30))/ALGBRC(37)
ALGBRC(38) = 1.00000/(1.00000+EXP((ALGBRC(8)+28.6000)/17.1000))
ALGBRC(23) = 1.00000/( 100.000*EXP(- ALGBRC(8)/54.6450)+ 656.000*EXP(ALGBRC(8)/106.157))
RATES(31) = (ALGBRC(38) - STATES(31))/ALGBRC(23)
ALGBRC(10) = ALGBRC(8)+41.0000
ALGBRC(27) = TERNRY(ABS(ALGBRC(10)).LT.CONSTS(37), 2000.00, ( 200.000*ALGBRC(10))/(1.00000 - EXP( - 0.100000*ALGBRC(10))))
ALGBRC(41) = 8000.00*EXP( - 0.0560000*(ALGBRC(8)+66.0000))
RATES(5) = ALGBRC(27)*(1.00000 - STATES(5)) - ALGBRC(41)*STATES(5)
ALGBRC(25) = 10.0000*EXP( 0.0133000*(ALGBRC(8)+40.0000))
ALGBRC(40) = CONSTS(104)/(CONSTS(104)+ALGBRC(25))
ALGBRC(45) = 1.00000/(CONSTS(104)+ALGBRC(25))
RATES(33) = (ALGBRC(40) - STATES(33))/ALGBRC(45)
ALGBRC(24) = (14.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/12.0000)))/(14.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/12.0000))+ 1.00000*EXP(- (ALGBRC(8) - CONSTS(103))/45.0000))
ALGBRC(39) = 28.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/3.00000))
ALGBRC(44) = 1.00000*EXP(- (((ALGBRC(8) - CONSTS(103)) - CONSTS(81)) - 5.00000)/25.0000)
ALGBRC(48) = 1.00000/(ALGBRC(39)+ALGBRC(44))
RATES(32) = (ALGBRC(24) - STATES(32))/ALGBRC(48)
ALGBRC(12) = 1.00000/(1.00000+EXP(- ((ALGBRC(8)+20.3000) - CONSTS(96))/( CONSTS(97)*4.20000)))
ALGBRC(29) = TERNRY(ALGBRC(8).EQ.- 41.8000, - 41.8000, TERNRY(ALGBRC(8).EQ.0.00000, 0.00000, TERNRY(ALGBRC(8).EQ.- 6.80000, - 6.80001, ALGBRC(8))
ALGBRC(42) = ( - 0.0283900*((ALGBRC(29)+41.8000) - CONSTS(96)))/(EXP(- ((ALGBRC(29)+41.8000) - CONSTS(96))/2.50000) - 1.00000) - ( 0.0849000*((ALGBRC(29)+6.80000) - CONSTS(96)))/(EXP(- ((ALGBRC(29)+6.80000) - CONSTS(96))/4.80000) - 1.00000)
ALGBRC(46) = TERNRY(ALGBRC(8).EQ.- 1.80000, - 1.80001, ALGBRC(8))
ALGBRC(49) = ( 0.0114300*((ALGBRC(46)+1.80000) - CONSTS(96)))/(EXP(((ALGBRC(46)+1.80000) - CONSTS(96))/2.50000) - 1.00000)
ALGBRC(51) = 0.00100000/(ALGBRC(42)+ALGBRC(49))
RATES(7) = (ALGBRC(12) - STATES(7))/ALGBRC(51)
ALGBRC(17) = TERNRY(CONSTS(14).GT.0.00000, 7.50000, STATES(3))
ALGBRC(33) = CONSTS(86)*log(CONSTS(15)/ALGBRC(17))
ALGBRC(53) = CONSTS(92)*CONSTS(22)*1.00000+CONSTS(20)/CONSTS(17) ** 1.20000 ** - 1.00000*1.00000+CONSTS(21)/ALGBRC(17) ** 1.30000 ** - 1.00000*1.00000+EXP(- ((ALGBRC(8) - ALGBRC(33))+110.000)/20.0000) ** - 1.00000
ALGBRC(55) = EXP(( - CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86)))
ALGBRC(61) = 1.00000+ (CONSTS(18)/CONSTS(35))*(1.00000+EXP(( CONSTS(26)*ALGBRC(8))/CONSTS(86)))+ (CONSTS(15)/CONSTS(33))*(1.00000+ (CONSTS(15)/CONSTS(34))*(1.00000+CONSTS(15)/CONSTS(32)))
ALGBRC(63) = ( (( (CONSTS(15)/CONSTS(33))*CONSTS(15))/CONSTS(34))*(1.00000+CONSTS(15)/CONSTS(32))*EXP(( - CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86))))/ALGBRC(61)
ALGBRC(62) = ( (CONSTS(18)/CONSTS(35))*EXP(( CONSTS(26)*ALGBRC(8))/CONSTS(86)))/ALGBRC(61)
ALGBRC(59) = EXP(( CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86)))
ALGBRC(54) = ALGBRC(17)/(CONSTS(27)+ALGBRC(17))
ALGBRC(64) = ALGBRC(55)*CONSTS(93)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(ALGBRC(54)+ALGBRC(55))
ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(28))*(1.00000+EXP(( - CONSTS(24)*ALGBRC(8))/CONSTS(86))+ALGBRC(17)/CONSTS(31))+ (ALGBRC(17)/CONSTS(29))*(1.00000+ (ALGBRC(17)/CONSTS(30))*(1.00000+ALGBRC(17)/CONSTS(27)))
ALGBRC(57) = ( (STATES(2)/CONSTS(28))*EXP(( - CONSTS(24)*ALGBRC(8))/CONSTS(86)))/ALGBRC(56)
ALGBRC(58) = ( (( (ALGBRC(17)/CONSTS(29))*ALGBRC(17))/CONSTS(30))*(1.00000+ALGBRC(17)/CONSTS(27))*EXP(( CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86))))/ALGBRC(56)
ALGBRC(60) = ALGBRC(59)*ALGBRC(54)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(93)+ALGBRC(59))
ALGBRC(65) = ALGBRC(58)*ALGBRC(54)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(ALGBRC(54)+ALGBRC(55))
ALGBRC(66) = ALGBRC(63)*CONSTS(93)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(93)+ALGBRC(59))
ALGBRC(67) = ( CONSTS(23)*( ALGBRC(60)*ALGBRC(62) - ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
ALGBRC(68) = CONSTS(86)*log((CONSTS(15)+ 0.120000*CONSTS(17))/(ALGBRC(17)+ 0.120000*CONSTS(16)))
ALGBRC(69) = CONSTS(36)*STATES(5) ** 3.00000*STATES(6)*(ALGBRC(8) - ALGBRC(68))
ALGBRC(43) = TERNRY(CONSTS(11).GE.1.00000, (10.6015/5.00000)/(10.6015/5.00000+EXP(( - 0.710000*ALGBRC(8))/25.0000)), 1.00000)
ALGBRC(47) = (( STATES(4) ** 2.00000*CONSTS(17))/(CONSTS(17)+CONSTS(19)))*CONSTS(83)*(ALGBRC(8) - ALGBRC(33))*ALGBRC(43)
ALGBRC(72) = (( 1.85000e-05*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))))*(ALGBRC(17) - CONSTS(15)*EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
RATES(3) = ( - 1.00000*(ALGBRC(69)+ALGBRC(47)+ALGBRC(72)+ 3.00000*ALGBRC(53)+ 3.00000*ALGBRC(67)))/( 1.00000*(CONSTS(101)+CONSTS(99))*CONSTS(3))
ALGBRC(82) = CONSTS(59)*STATES(18)*(1.00000 - STATES(22)) - CONSTS(64)*STATES(22)
RATES(22) = ALGBRC(82)
ALGBRC(80) = CONSTS(89)/(1.00000+CONSTS(54)/STATES(18))
ALGBRC(83) = (STATES(17) - STATES(16))/CONSTS(52)
RATES(17) = ALGBRC(80) - ( ALGBRC(83)*CONSTS(100))/CONSTS(102)
ALGBRC(85) = CONSTS(61)*STATES(18)*(1.00000 - (STATES(23)+STATES(19))) - CONSTS(65)*STATES(23)
RATES(23) = ALGBRC(85)
ALGBRC(86) = CONSTS(63)*STATES(16)*(1.00000 - STATES(24)) - CONSTS(68)*STATES(24)
RATES(24) = ALGBRC(86)
ALGBRC(76) = CONSTS(42)*STATES(13)*(STATES(16) - STATES(2))
RATES(16) = ALGBRC(83) - (ALGBRC(76)+ CONSTS(58)*ALGBRC(86))
ALGBRC(50) = (( STATES(4) ** 2.00000*CONSTS(17))/(CONSTS(17)+CONSTS(19)))*CONSTS(84)*(ALGBRC(8) - CONSTS(88))*ALGBRC(43)
ALGBRC(52) = ALGBRC(47)+ALGBRC(50)
ALGBRC(77) = CONSTS(80)*(ALGBRC(8) - CONSTS(88))*( 0.900000*STATES(30)+ 0.100000*STATES(29))*STATES(31)
ALGBRC(79) = CONSTS(86)*log((CONSTS(17)+ 0.00000*CONSTS(15))/(CONSTS(16)+ 0.00000*ALGBRC(17)))
ALGBRC(81) = CONSTS(87)*(ALGBRC(8) - ALGBRC(79))*STATES(32) ** 2.00000
ALGBRC(75) = CONSTS(79)*(ALGBRC(8) - CONSTS(88))*STATES(27)*STATES(28)
ALGBRC(70) = (( 2.00000*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000)*2.00000)/CONSTS(86)))))*(STATES(2) - CONSTS(18)*EXP(( - 2.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
ALGBRC(71) = (( 0.000365000*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))))*(CONSTS(16) - CONSTS(17)*EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
ALGBRC(73) = (ALGBRC(70)+ALGBRC(71)+ALGBRC(72))*(1.00000 - CONSTS(95))*1.00000*CONSTS(94)
ALGBRC(74) = (( 2.00000*CONSTS(41)*ALGBRC(8))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*ALGBRC(8)*2.00000)/CONSTS(86)))))*(STATES(2) - CONSTS(18)*EXP(( - 2.00000*ALGBRC(8))/CONSTS(86)))*STATES(10)*STATES(11)
ALGBRC(84) = TERNRY(CONSTS(12).GT.0.00000, CONSTS(82)*(ALGBRC(8) - CONSTS(88))*(1.00000+EXP((ALGBRC(8)+20.0000)/20.0000))*STATES(33), 0.00000)
ALGBRC(87) = ALGBRC(52)+ALGBRC(77)+ALGBRC(81)+ALGBRC(75)+ALGBRC(53)+ALGBRC(67)+ALGBRC(69)+ALGBRC(73)+ALGBRC(74)+ALGBRC(84)
RATES(1) = - ALGBRC(87)/CONSTS(4)
ALGBRC(88) = CONSTS(62)*STATES(18)*(1.00000 - STATES(20)) - CONSTS(67)*STATES(20)
RATES(20) = ALGBRC(88)
ALGBRC(89) = CONSTS(62)*STATES(2)*(1.00000 - STATES(21)) - CONSTS(67)*STATES(21)
RATES(21) = ALGBRC(89)
ALGBRC(78) = (STATES(2) - STATES(18))/CONSTS(51)
ALGBRC(90) = TERNRY(CONSTS(14).GT.0.00000.AND.VOI.GT.CONSTS(72), 10.0000, 0.00000)
RATES(18) = (( 1.00000*( ALGBRC(78)*CONSTS(99) - ALGBRC(80)*CONSTS(102)))/CONSTS(101) - ( CONSTS(57)*ALGBRC(88)+ CONSTS(55)*ALGBRC(82)+ CONSTS(56)*ALGBRC(85))) - ( CONSTS(70)*STATES(18)*(ALGBRC(90) - STATES(25)) - CONSTS(71)*STATES(25))
RATES(25) = CONSTS(70)*STATES(18)*(ALGBRC(90) - STATES(25)) - CONSTS(71)*STATES(25)
RATES(2) = (( ALGBRC(76)*CONSTS(100))/CONSTS(99) - (((ALGBRC(70)+ALGBRC(74)) - 2.00000*ALGBRC(67))/( 2.00000*CONSTS(3)*CONSTS(99))+ALGBRC(78)+ CONSTS(57)*ALGBRC(89))) - ( CONSTS(70)*STATES(2)*(ALGBRC(90) - STATES(26)) - CONSTS(71)*STATES(26))
RATES(26) = CONSTS(70)*STATES(2)*(ALGBRC(90) - STATES(26)) - CONSTS(71)*STATES(26)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = CONSTS(60)*CONSTS(69)*(1.00000 - (STATES(23)+STATES(19))) - CONSTS(66)*STATES(19)
ALGBRC(2) = CONSTS(40)/(CONSTS(40)+STATES(2))
ALGBRC(6) = ( 0.00100000*ALGBRC(2))/CONSTS(39)
ALGBRC(3) = CONSTS(43) - (CONSTS(43) - CONSTS(44))/(1.00000+CONSTS(45)/STATES(16) ** CONSTS(46))
ALGBRC(7) = CONSTS(47)/ALGBRC(3)
ALGBRC(16) = CONSTS(48)*ALGBRC(3)
ALGBRC(4) = TERNRY(VOI.GT.CONSTS(6).AND.VOI.LT.CONSTS(6)+CONSTS(7), CONSTS(8), CONSTS(9))
ALGBRC(8) = TERNRY(CONSTS(5).GE.1.00000, ALGBRC(4), STATES(1))
ALGBRC(9) = 0.716653/( 0.0708000*EXP(- (((ALGBRC(8)+5.00000) - CONSTS(90)) - CONSTS(91))/20.2791)+ 10.6000*EXP(((ALGBRC(8) - CONSTS(90)) - CONSTS(91))/18.0000))
ALGBRC(26) = 1.00000/(1.00000+EXP((((ALGBRC(8)+52.5000) - CONSTS(90)) - CONSTS(91))/9.00000))
ALGBRC(11) = 20.0000*EXP( - 0.125000*(ALGBRC(8)+75.0000))
ALGBRC(28) = 2000.00/( 320.000*EXP( - 0.100000*(ALGBRC(8)+75.0000))+1.00000)
ALGBRC(13) = 1.00000/(1.00000+EXP((ALGBRC(8)+37.4000)/5.30000))
ALGBRC(30) = 0.00100000*(44.3000+ 230.000*EXP(- (ALGBRC(8)+36.0000)/10.0000 ** 2.00000))
ALGBRC(14) = 1.00000/(1.00000+EXP(- (ALGBRC(8)+38.3000)/5.50000))
ALGBRC(31) = 0.00100000/( 1.06800*EXP((ALGBRC(8)+38.3000)/30.0000)+ 1.06800*EXP(- (ALGBRC(8)+38.3000)/30.0000))
ALGBRC(15) = 1.00000/(1.00000+EXP((ALGBRC(8)+58.7000)/3.80000))
ALGBRC(32) = 1.00000/( 16.6700*EXP(- (ALGBRC(8)+75.0000)/83.3000)+ 16.6700*EXP((ALGBRC(8)+75.0000)/15.3800))
ALGBRC(18) = 1.00000/(1.00000+EXP((ALGBRC(8)+49.0000)/13.0000))
ALGBRC(34) = 0.00100000*0.600000*(65.1700/( 0.570000*EXP( - 0.0800000*(ALGBRC(8)+44.0000))+ 0.0650000*EXP( 0.100000*(ALGBRC(8)+45.9300)))+10.1000)
ALGBRC(19) = 1.00000/(1.00000+EXP(- (ALGBRC(8) - 19.3000)/15.0000))
ALGBRC(35) = 0.00100000*0.660000*1.40000*(15.5900/( 1.03700*EXP( 0.0900000*(ALGBRC(8)+30.6100))+ 0.369000*EXP( - 0.120000*(ALGBRC(8)+23.8400)))+2.98000)
ALGBRC(20) = 1.00000/(1.00000+EXP(- (ALGBRC(8)+14.8000)/8.50000))
ALGBRC(36) = 0.846554/( 4.20000*EXP(ALGBRC(8)/17.0000)+ 0.150000*EXP(- ALGBRC(8)/21.6000))
ALGBRC(37) = 1.00000/( 30.0000*EXP(ALGBRC(8)/10.0000)+EXP(- ALGBRC(8)/12.0000))
ALGBRC(38) = 1.00000/(1.00000+EXP((ALGBRC(8)+28.6000)/17.1000))
ALGBRC(23) = 1.00000/( 100.000*EXP(- ALGBRC(8)/54.6450)+ 656.000*EXP(ALGBRC(8)/106.157))
ALGBRC(10) = ALGBRC(8)+41.0000
ALGBRC(27) = TERNRY(ABS(ALGBRC(10)).LT.CONSTS(37), 2000.00, ( 200.000*ALGBRC(10))/(1.00000 - EXP( - 0.100000*ALGBRC(10))))
ALGBRC(41) = 8000.00*EXP( - 0.0560000*(ALGBRC(8)+66.0000))
ALGBRC(25) = 10.0000*EXP( 0.0133000*(ALGBRC(8)+40.0000))
ALGBRC(40) = CONSTS(104)/(CONSTS(104)+ALGBRC(25))
ALGBRC(45) = 1.00000/(CONSTS(104)+ALGBRC(25))
ALGBRC(24) = (14.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/12.0000)))/(14.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/12.0000))+ 1.00000*EXP(- (ALGBRC(8) - CONSTS(103))/45.0000))
ALGBRC(39) = 28.0000/(1.00000+EXP(- ((ALGBRC(8) - 40.0000) - CONSTS(103))/3.00000))
ALGBRC(44) = 1.00000*EXP(- (((ALGBRC(8) - CONSTS(103)) - CONSTS(81)) - 5.00000)/25.0000)
ALGBRC(48) = 1.00000/(ALGBRC(39)+ALGBRC(44))
ALGBRC(12) = 1.00000/(1.00000+EXP(- ((ALGBRC(8)+20.3000) - CONSTS(96))/( CONSTS(97)*4.20000)))
ALGBRC(29) = TERNRY(ALGBRC(8).EQ.- 41.8000, - 41.8000, TERNRY(ALGBRC(8).EQ.0.00000, 0.00000, TERNRY(ALGBRC(8).EQ.- 6.80000, - 6.80001, ALGBRC(8))
ALGBRC(42) = ( - 0.0283900*((ALGBRC(29)+41.8000) - CONSTS(96)))/(EXP(- ((ALGBRC(29)+41.8000) - CONSTS(96))/2.50000) - 1.00000) - ( 0.0849000*((ALGBRC(29)+6.80000) - CONSTS(96)))/(EXP(- ((ALGBRC(29)+6.80000) - CONSTS(96))/4.80000) - 1.00000)
ALGBRC(46) = TERNRY(ALGBRC(8).EQ.- 1.80000, - 1.80001, ALGBRC(8))
ALGBRC(49) = ( 0.0114300*((ALGBRC(46)+1.80000) - CONSTS(96)))/(EXP(((ALGBRC(46)+1.80000) - CONSTS(96))/2.50000) - 1.00000)
ALGBRC(51) = 0.00100000/(ALGBRC(42)+ALGBRC(49))
ALGBRC(17) = TERNRY(CONSTS(14).GT.0.00000, 7.50000, STATES(3))
ALGBRC(33) = CONSTS(86)*log(CONSTS(15)/ALGBRC(17))
ALGBRC(53) = CONSTS(92)*CONSTS(22)*1.00000+CONSTS(20)/CONSTS(17) ** 1.20000 ** - 1.00000*1.00000+CONSTS(21)/ALGBRC(17) ** 1.30000 ** - 1.00000*1.00000+EXP(- ((ALGBRC(8) - ALGBRC(33))+110.000)/20.0000) ** - 1.00000
ALGBRC(55) = EXP(( - CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86)))
ALGBRC(61) = 1.00000+ (CONSTS(18)/CONSTS(35))*(1.00000+EXP(( CONSTS(26)*ALGBRC(8))/CONSTS(86)))+ (CONSTS(15)/CONSTS(33))*(1.00000+ (CONSTS(15)/CONSTS(34))*(1.00000+CONSTS(15)/CONSTS(32)))
ALGBRC(63) = ( (( (CONSTS(15)/CONSTS(33))*CONSTS(15))/CONSTS(34))*(1.00000+CONSTS(15)/CONSTS(32))*EXP(( - CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86))))/ALGBRC(61)
ALGBRC(62) = ( (CONSTS(18)/CONSTS(35))*EXP(( CONSTS(26)*ALGBRC(8))/CONSTS(86)))/ALGBRC(61)
ALGBRC(59) = EXP(( CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86)))
ALGBRC(54) = ALGBRC(17)/(CONSTS(27)+ALGBRC(17))
ALGBRC(64) = ALGBRC(55)*CONSTS(93)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(62)*ALGBRC(59)*(ALGBRC(54)+ALGBRC(55))
ALGBRC(56) = 1.00000+ (STATES(2)/CONSTS(28))*(1.00000+EXP(( - CONSTS(24)*ALGBRC(8))/CONSTS(86))+ALGBRC(17)/CONSTS(31))+ (ALGBRC(17)/CONSTS(29))*(1.00000+ (ALGBRC(17)/CONSTS(30))*(1.00000+ALGBRC(17)/CONSTS(27)))
ALGBRC(57) = ( (STATES(2)/CONSTS(28))*EXP(( - CONSTS(24)*ALGBRC(8))/CONSTS(86)))/ALGBRC(56)
ALGBRC(58) = ( (( (ALGBRC(17)/CONSTS(29))*ALGBRC(17))/CONSTS(30))*(1.00000+ALGBRC(17)/CONSTS(27))*EXP(( CONSTS(25)*ALGBRC(8))/( 2.00000*CONSTS(86))))/ALGBRC(56)
ALGBRC(60) = ALGBRC(59)*ALGBRC(54)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(55)*ALGBRC(57)*(CONSTS(93)+ALGBRC(59))
ALGBRC(65) = ALGBRC(58)*ALGBRC(54)*(ALGBRC(63)+ALGBRC(62))+ ALGBRC(57)*ALGBRC(63)*(ALGBRC(54)+ALGBRC(55))
ALGBRC(66) = ALGBRC(63)*CONSTS(93)*(ALGBRC(58)+ALGBRC(57))+ ALGBRC(58)*ALGBRC(62)*(CONSTS(93)+ALGBRC(59))
ALGBRC(67) = ( CONSTS(23)*( ALGBRC(60)*ALGBRC(62) - ALGBRC(64)*ALGBRC(57)))/(ALGBRC(64)+ALGBRC(60)+ALGBRC(65)+ALGBRC(66))
ALGBRC(68) = CONSTS(86)*log((CONSTS(15)+ 0.120000*CONSTS(17))/(ALGBRC(17)+ 0.120000*CONSTS(16)))
ALGBRC(69) = CONSTS(36)*STATES(5) ** 3.00000*STATES(6)*(ALGBRC(8) - ALGBRC(68))
ALGBRC(43) = TERNRY(CONSTS(11).GE.1.00000, (10.6015/5.00000)/(10.6015/5.00000+EXP(( - 0.710000*ALGBRC(8))/25.0000)), 1.00000)
ALGBRC(47) = (( STATES(4) ** 2.00000*CONSTS(17))/(CONSTS(17)+CONSTS(19)))*CONSTS(83)*(ALGBRC(8) - ALGBRC(33))*ALGBRC(43)
ALGBRC(72) = (( 1.85000e-05*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))))*(ALGBRC(17) - CONSTS(15)*EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
ALGBRC(82) = CONSTS(59)*STATES(18)*(1.00000 - STATES(22)) - CONSTS(64)*STATES(22)
ALGBRC(80) = CONSTS(89)/(1.00000+CONSTS(54)/STATES(18))
ALGBRC(83) = (STATES(17) - STATES(16))/CONSTS(52)
ALGBRC(85) = CONSTS(61)*STATES(18)*(1.00000 - (STATES(23)+STATES(19))) - CONSTS(65)*STATES(23)
ALGBRC(86) = CONSTS(63)*STATES(16)*(1.00000 - STATES(24)) - CONSTS(68)*STATES(24)
ALGBRC(76) = CONSTS(42)*STATES(13)*(STATES(16) - STATES(2))
ALGBRC(50) = (( STATES(4) ** 2.00000*CONSTS(17))/(CONSTS(17)+CONSTS(19)))*CONSTS(84)*(ALGBRC(8) - CONSTS(88))*ALGBRC(43)
ALGBRC(52) = ALGBRC(47)+ALGBRC(50)
ALGBRC(77) = CONSTS(80)*(ALGBRC(8) - CONSTS(88))*( 0.900000*STATES(30)+ 0.100000*STATES(29))*STATES(31)
ALGBRC(79) = CONSTS(86)*log((CONSTS(17)+ 0.00000*CONSTS(15))/(CONSTS(16)+ 0.00000*ALGBRC(17)))
ALGBRC(81) = CONSTS(87)*(ALGBRC(8) - ALGBRC(79))*STATES(32) ** 2.00000
ALGBRC(75) = CONSTS(79)*(ALGBRC(8) - CONSTS(88))*STATES(27)*STATES(28)
ALGBRC(70) = (( 2.00000*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000)*2.00000)/CONSTS(86)))))*(STATES(2) - CONSTS(18)*EXP(( - 2.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
ALGBRC(71) = (( 0.000365000*CONSTS(38)*(ALGBRC(8) - 0.00000))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))))*(CONSTS(16) - CONSTS(17)*EXP(( - 1.00000*(ALGBRC(8) - 0.00000))/CONSTS(86)))*STATES(7)*STATES(8)*STATES(9)
ALGBRC(73) = (ALGBRC(70)+ALGBRC(71)+ALGBRC(72))*(1.00000 - CONSTS(95))*1.00000*CONSTS(94)
ALGBRC(74) = (( 2.00000*CONSTS(41)*ALGBRC(8))/( CONSTS(86)*(1.00000 - EXP(( - 1.00000*ALGBRC(8)*2.00000)/CONSTS(86)))))*(STATES(2) - CONSTS(18)*EXP(( - 2.00000*ALGBRC(8))/CONSTS(86)))*STATES(10)*STATES(11)
ALGBRC(84) = TERNRY(CONSTS(12).GT.0.00000, CONSTS(82)*(ALGBRC(8) - CONSTS(88))*(1.00000+EXP((ALGBRC(8)+20.0000)/20.0000))*STATES(33), 0.00000)
ALGBRC(87) = ALGBRC(52)+ALGBRC(77)+ALGBRC(81)+ALGBRC(75)+ALGBRC(53)+ALGBRC(67)+ALGBRC(69)+ALGBRC(73)+ALGBRC(74)+ALGBRC(84)
ALGBRC(88) = CONSTS(62)*STATES(18)*(1.00000 - STATES(20)) - CONSTS(67)*STATES(20)
ALGBRC(89) = CONSTS(62)*STATES(2)*(1.00000 - STATES(21)) - CONSTS(67)*STATES(21)
ALGBRC(78) = (STATES(2) - STATES(18))/CONSTS(51)
ALGBRC(90) = TERNRY(CONSTS(14).GT.0.00000.AND.VOI.GT.CONSTS(72), 10.0000, 0.00000)
ALGBRC(1) = 0.500000*CONSTS(86)*log(CONSTS(18)/STATES(2))
ALGBRC(21) = (1.00000/(1.00000+EXP(- (ALGBRC(8)+23.2000)/6.60000)))/(0.846554/( 37.2000*EXP(ALGBRC(8)/11.9000)+ 0.960000*EXP(- ALGBRC(8)/18.5000)))
ALGBRC(22) = 4.00000*(( 37.2000*EXP(ALGBRC(8)/15.9000)+ 0.960000*EXP(- ALGBRC(8)/22.5000))/0.846554 - (1.00000/(1.00000+EXP(- (ALGBRC(8)+23.2000)/10.6000)))/(0.846554/( 37.2000*EXP(ALGBRC(8)/15.9000)+ 0.960000*EXP(- ALGBRC(8)/22.5000))))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END