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