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 89 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 151 entries in the constant variable array.
C
C
C CONSTS(1) is Buf_C in component cai (mM).
C STATES(1) is Cai in component cai (mM).
C ALGBRC(1) is Cai_bufc in component cai (dimensionless).
C CONSTS(2) is Cm in component geom (pF).
C CONSTS(3) is F in component phys (C_per_mmol).
C CONSTS(4) is Kbuf_C in component cai (mM).
C CONSTS(145) is Vc in component geom (um3).
C ALGBRC(66) is i_CaL_Ca in component ical (A_per_F).
C ALGBRC(72) is i_CaT in component icat (A_per_F).
C ALGBRC(84) is i_NaCa in component inaca (A_per_F).
C ALGBRC(86) is i_PCa in component ipca (A_per_F).
C ALGBRC(63) is i_b_Ca in component ibca (A_per_F).
C ALGBRC(82) is i_leak in component ileak (mM_per_ms).
C ALGBRC(87) is i_rel in component irel (mM_per_ms).
C ALGBRC(89) is i_up in component iup (mM_per_ms).
C VOI is time in component engine (ms).
C CONSTS(87) is Buf_SR in component casr (mM).
C STATES(2) is Ca_SR in component casr (mM).
C ALGBRC(17) is Ca_SR_bufSR in component casr (dimensionless).
C CONSTS(5) is Kbuf_SR in component casr (mM).
C CONSTS(127) is V_SR in component geom (um3).
C CONSTS(6) is pace in component engine (dimensionless).
C CONSTS(7) is Cao in component extra (mM).
C ALGBRC(33) is E_Ca in component erev (mV).
C ALGBRC(47) is E_K in component erev (mV).
C ALGBRC(61) is E_Na in component erev (mV).
C STATES(3) is Ki in component ki (mM).
C CONSTS(8) is Ko in component extra (mM).
C STATES(4) is Nai in component nai (mM).
C CONSTS(9) is Nao in component extra (mM).
C CONSTS(88) is RTF in component phys (mV).
C CONSTS(10) is VSR_tenT in component geom (um3).
C CONSTS(11) is V_tot in component geom (um3).
C CONSTS(86) is V_tot_tenT in component geom (um3).
C CONSTS(12) is Vc_tenT in component geom (um3).
C STATES(5) is V in component membrane (mV).
C CONSTS(89) is g_b_Ca in component ibca (mS_per_uF).
C CONSTS(90) is g_b_Na in component ibna (mS_per_uF).
C ALGBRC(64) is i_b_Na in component ibna (A_per_F).
C CONSTS(128) is FFRT in component phys (s4_A2_per_g_per_m2_per_mol_times_1e3).
C CONSTS(91) is FRT in component phys (per_mV).
C ALGBRC(4) is alpha_fCa in component ical (dimensionless).
C ALGBRC(20) is beta_fCa in component ical (dimensionless).
C STATES(6) is d in component ical (dimensionless).
C CONSTS(13) is d1 in component ical (mS_per_uF).
C CONSTS(14) is d2 in component ical (mV).
C CONSTS(92) is d3 in component ical (mS_per_uF).
C CONSTS(129) is d4 in component ical (mV).
C CONSTS(15) is d5 in component ical (dimensionless).
C CONSTS(93) is d6 in component ical (mV).
C STATES(7) is f in component ical (dimensionless).
C CONSTS(16) is f1 in component ical (mS_per_uF).
C CONSTS(94) is f2 in component ical (mV).
C CONSTS(95) is f3 in component ical (mS_per_uF).
C CONSTS(130) is f4 in component ical (mV).
C CONSTS(17) is f5 in component ical (dimensionless).
C CONSTS(18) is f6 in component ical (mV).
C STATES(8) is fCa in component ical (dimensionless).
C ALGBRC(50) is fCa_inf in component ical (dimensionless).
C ALGBRC(36) is gamma_fCa in component ical (dimensionless).
C ALGBRC(71) is i_CaL in component ical (A_per_F).
C ALGBRC(68) is i_CaL_K in component ical (A_per_F).
C ALGBRC(70) is i_CaL_Na in component ical (A_per_F).
C ALGBRC(65) is ibarca in component ical (A_per_F).
C ALGBRC(67) is ibark in component ical (A_per_F).
C ALGBRC(69) is ibarna in component ical (A_per_F).
C ALGBRC(2) is ical_d_a in component ical (mS_per_uF).
C ALGBRC(18) is ical_d_b in component ical (mS_per_uF).
C ALGBRC(34) is ical_d_inf in component ical (dimensionless).
C ALGBRC(48) is ical_d_tau in component ical (ms).
C ALGBRC(3) is ical_f_a in component ical (mS_per_uF).
C ALGBRC(19) is ical_f_b in component ical (mS_per_uF).
C ALGBRC(35) is ical_f_inf in component ical (dimensionless).
C ALGBRC(49) is ical_f_tau in component ical (ms).
C ALGBRC(62) is k_fca in component ical (dimensionless).
C CONSTS(19) is p_CaL in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(146) is p_CaL_Ca in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(147) is p_CaL_K in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(148) is p_CaL_Na in component ical (L_per_F_per_ms_times_1e0).
C CONSTS(20) is p_CaL_shannonCa in component ical (dimensionless).
C CONSTS(131) is p_CaL_shannonCap in component ical (dimensionless).
C CONSTS(21) is p_CaL_shannonK in component ical (dimensionless).
C CONSTS(132) is p_CaL_shannonKp in component ical (dimensionless).
C CONSTS(22) is p_CaL_shannonNa in component ical (dimensionless).
C CONSTS(133) is p_CaL_shannonNap in component ical (dimensionless).
C CONSTS(96) is p_CaL_shannonTot in component ical (dimensionless).
C CONSTS(23) is scale in component ical (dimensionless).
C CONSTS(24) is tau_fCa in component ical (ms).
C CONSTS(25) is taud_const in component ical (ms).
C CONSTS(26) is tauf_const in component ical (ms).
C STATES(9) is d in component icat (dimensionless).
C STATES(10) is f in component icat (dimensionless).
C CONSTS(27) is g_CaT in component icat (mS_per_uF).
C ALGBRC(5) is icat_d_inf in component icat (dimensionless).
C ALGBRC(21) is icat_d_tau in component icat (ms).
C ALGBRC(6) is icat_f_inf in component icat (dimensionless).
C ALGBRC(22) is icat_f_tau in component icat (ms).
C CONSTS(97) is Na_frac in component ifunny (dimensionless).
C CONSTS(28) is NatoK_ratio in component ifunny (dimensionless).
C STATES(11) is Xf in component ifunny (dimensionless).
C CONSTS(29) is g_f in component ifunny (mS_per_uF).
C ALGBRC(75) is i_f in component ifunny (A_per_F).
C ALGBRC(73) is i_fK in component ifunny (A_per_F).
C ALGBRC(74) is i_fNa in component ifunny (A_per_F).
C ALGBRC(7) is ifunny_Xf_a in component ifunny (mS_per_uF).
C ALGBRC(23) is ifunny_Xf_b in component ifunny (mS_per_uF).
C ALGBRC(37) is ifunny_Xf_inf in component ifunny (dimensionless).
C ALGBRC(51) is ifunny_Xf_tau in component ifunny (ms).
C CONSTS(30) is xF1 in component ifunny (mS_per_uF).
C CONSTS(98) is xF2 in component ifunny (mV).
C CONSTS(99) is xF3 in component ifunny (mS_per_uF).
C CONSTS(134) is xF4 in component ifunny (mV).
C CONSTS(31) is xF5 in component ifunny (dimensionless).
C CONSTS(32) is xF6 in component ifunny (mV).
C CONSTS(33) is xF_const in component ifunny (ms).
C CONSTS(34) is g_K1 in component ik1 (mS_per_uF).
C ALGBRC(79) is i_K1 in component ik1 (A_per_F).
C ALGBRC(76) is ik1_inf_a in component ik1 (mS_per_uF).
C ALGBRC(77) is ik1_inf_b in component ik1 (mS_per_uF).
C ALGBRC(78) is inf in component ik1 (dimensionless).
C CONSTS(35) is xK11 in component ik1 (mS_per_uF).
C CONSTS(36) is xK12 in component ik1 (mV).
C CONSTS(37) is xK13 in component ik1 (mV).
C CONSTS(38) is xK14 in component ik1 (mV).
C CONSTS(39) is xK15 in component ik1 (mV).
C STATES(12) is Xr1 in component ikr (dimensionless).
C CONSTS(40) is Xr1_1 in component ikr (mS_per_uF).
C CONSTS(41) is Xr1_2 in component ikr (mV).
C CONSTS(100) is Xr1_3 in component ikr (mS_per_uF).
C CONSTS(135) is Xr1_4 in component ikr (mV).
C CONSTS(42) is Xr1_5 in component ikr (dimensionless).
C CONSTS(101) is Xr1_6 in component ikr (mV).
C STATES(13) is Xr2 in component ikr (dimensionless).
C CONSTS(43) is Xr2_1 in component ikr (mS_per_uF).
C CONSTS(102) is Xr2_2 in component ikr (mV).
C CONSTS(103) is Xr2_3 in component ikr (mS_per_uF).
C CONSTS(136) is Xr2_4 in component ikr (mV).
C CONSTS(44) is Xr2_5 in component ikr (dimensionless).
C CONSTS(45) is Xr2_6 in component ikr (mV).
C CONSTS(46) is g_Kr in component ikr (mS_per_uF).
C ALGBRC(80) is i_Kr in component ikr (A_per_F).
C ALGBRC(8) is ikr_Xr1_a in component ikr (mS_per_uF).
C ALGBRC(24) is ikr_Xr1_b in component ikr (mS_per_uF).
C ALGBRC(38) is ikr_Xr1_inf in component ikr (dimensionless).
C ALGBRC(52) is ikr_Xr1_tau in component ikr (ms).
C ALGBRC(9) is ikr_Xr2_a in component ikr (mS_per_uF).
C ALGBRC(25) is ikr_Xr2_b in component ikr (mS_per_uF).
C ALGBRC(39) is ikr_Xr2_inf in component ikr (dimensionless).
C ALGBRC(53) is ikr_Xr2_tau in component ikr (ms).
C CONSTS(47) is tau_1_offset in component ikr (ms).
C CONSTS(48) is tau_2_offset in component ikr (ms).
C STATES(14) is Xs in component iks (dimensionless).
C CONSTS(49) is g_Ks in component iks (mS_per_uF).
C ALGBRC(81) is i_Ks in component iks (A_per_F).
C ALGBRC(10) is iks_Xs_a in component iks (mS_per_uF).
C ALGBRC(26) is iks_Xs_b in component iks (mS_per_uF).
C ALGBRC(40) is iks_Xs_inf in component iks (dimensionless).
C ALGBRC(54) is iks_Xs_tau in component iks (ms).
C CONSTS(50) is ks1 in component iks (mS_per_uF).
C CONSTS(51) is ks2 in component iks (mV).
C CONSTS(104) is ks3 in component iks (mS_per_uF).
C CONSTS(137) is ks4 in component iks (mV).
C CONSTS(52) is ks5 in component iks (dimensionless).
C CONSTS(105) is ks6 in component iks (mV).
C CONSTS(53) is tauks_const in component iks (ms).
C CONSTS(106) is V_leak in component ileak (mS_per_uF).
C CONSTS(54) is g_Na in component ina (mS_per_uF).
C STATES(15) is h in component ina (dimensionless).
C CONSTS(55) is h1 in component ina (mS_per_uF).
C CONSTS(107) is h2 in component ina (mV).
C CONSTS(108) is h3 in component ina (mS_per_uF).
C CONSTS(138) is h4 in component ina (mV).
C CONSTS(56) is h5 in component ina (dimensionless).
C CONSTS(57) is h6 in component ina (mV).
C ALGBRC(83) is i_Na in component ina (A_per_F).
C ALGBRC(11) is ina_h_a in component ina (mS_per_uF).
C ALGBRC(27) is ina_h_b in component ina (mS_per_uF).
C ALGBRC(41) is ina_h_inf in component ina (dimensionless).
C ALGBRC(55) is ina_h_tau in component ina (ms).
C ALGBRC(12) is ina_j_a in component ina (mS_per_uF).
C ALGBRC(28) is ina_j_b in component ina (mS_per_uF).
C ALGBRC(42) is ina_j_inf in component ina (dimensionless).
C ALGBRC(56) is ina_j_tau in component ina (ms).
C ALGBRC(13) is ina_m_a in component ina (mS_per_uF).
C ALGBRC(29) is ina_m_b in component ina (mS_per_uF).
C ALGBRC(43) is ina_m_inf in component ina (dimensionless).
C ALGBRC(57) is ina_m_tau in component ina (ms).
C STATES(16) is j in component ina (dimensionless).
C CONSTS(58) is j1 in component ina (mS_per_uF).
C CONSTS(109) is j2 in component ina (mV).
C CONSTS(139) is j3 in component ina (mS_per_uF).
C CONSTS(149) is j4 in component ina (mV).
C CONSTS(110) is j5 in component ina (dimensionless).
C CONSTS(140) is j6 in component ina (mV).
C STATES(17) is m in component ina (dimensionless).
C CONSTS(59) is m1 in component ina (mS_per_uF).
C CONSTS(60) is m2 in component ina (mV).
C CONSTS(111) is m3 in component ina (mS_per_uF).
C CONSTS(141) is m4 in component ina (mV).
C CONSTS(61) is m5 in component ina (dimensionless).
C CONSTS(112) is m6 in component ina (mV).
C CONSTS(62) is tau_h_const in component ina (ms).
C CONSTS(63) is tau_j_const in component ina (ms).
C CONSTS(64) is tau_m_const in component ina (ms).
C CONSTS(65) is KmCa in component inaca (mM).
C CONSTS(66) is KmNai in component inaca (mM).
C CONSTS(67) is Ksat in component inaca (dimensionless).
C CONSTS(113) is alpha in component inaca (dimensionless).
C CONSTS(114) is gamma in component inaca (dimensionless).
C CONSTS(115) is kNaCa in component inaca (A_per_F).
C CONSTS(68) is Km_K in component inak (mM).
C CONSTS(69) is Km_Na in component inak (mM).
C CONSTS(116) is PNaK in component inak (A_per_F).
C ALGBRC(85) is i_NaK in component inak (A_per_F).
C CONSTS(70) is KPCa in component ipca (mM).
C CONSTS(117) is g_PCa in component ipca (A_per_F).
C STATES(18) is I in component irel (dimensionless).
C CONSTS(71) is MaxSR in component irel (dimensionless).
C CONSTS(72) is MinSR in component irel (dimensionless).
C STATES(19) is O in component irel (dimensionless).
C STATES(20) is R in component irel (dimensionless).
C ALGBRC(14) is RI in component irel (dimensionless).
C CONSTS(73) is ec50SR in component irel (mM).
C ALGBRC(30) is kCaSR in component irel (dimensionless).
C CONSTS(118) is kiCa in component irel (per_mM_per_ms).
C ALGBRC(44) is kiSRCa in component irel (per_mM_per_ms).
C CONSTS(119) is kim in component irel (mS_per_uF).
C CONSTS(142) is koCa in component irel (per_mM2_per_ms).
C ALGBRC(58) is koSRCa in component irel (per_mM2_per_ms).
C CONSTS(120) is kom in component irel (mS_per_uF).
C CONSTS(74) is ks in component irel (mS_per_uF).
C CONSTS(75) is g_to in component ito (mS_per_uF).
C ALGBRC(88) is i_to in component ito (A_per_F).
C ALGBRC(15) is ito_r_a in component ito (mS_per_uF).
C ALGBRC(31) is ito_r_b in component ito (mS_per_uF).
C ALGBRC(45) is ito_r_inf in component ito (dimensionless).
C ALGBRC(59) is ito_r_tau in component ito (ms).
C ALGBRC(16) is ito_s_a in component ito (mS_per_uF).
C ALGBRC(32) is ito_s_b in component ito (mS_per_uF).
C ALGBRC(46) is ito_s_inf in component ito (dimensionless).
C ALGBRC(60) is ito_s_tau in component ito (ms).
C STATES(21) is r in component ito (dimensionless).
C CONSTS(76) is r1 in component ito (mS_per_uF).
C CONSTS(77) is r2 in component ito (mV).
C CONSTS(121) is r3 in component ito (mS_per_uF).
C CONSTS(143) is r4 in component ito (mV).
C CONSTS(78) is r5 in component ito (dimensionless).
C CONSTS(122) is r6 in component ito (mV).
C STATES(22) is s in component ito (dimensionless).
C CONSTS(79) is s1 in component ito (mS_per_uF).
C CONSTS(123) is s2 in component ito (mV).
C CONSTS(124) is s3 in component ito (mS_per_uF).
C CONSTS(144) is s4 in component ito (mV).
C CONSTS(80) is s5 in component ito (dimensionless).
C CONSTS(81) is s6 in component ito (mV).
C CONSTS(82) is tau_r_const in component ito (ms).
C CONSTS(83) is tau_s_const in component ito (ms).
C CONSTS(125) is Kup in component iup (mM).
C CONSTS(126) is VmaxUp in component iup (mM_per_ms).
C CONSTS(151) is i_stim in component stimulus (A_per_F).
C CONSTS(84) is R in component phys (J_per_mol_per_K).
C CONSTS(85) is T in component phys (kelvin).
C CONSTS(150) is amplitude in component stimulus (A_per_F).
C RATES(1) is d/dt Cai in component cai (mM).
C RATES(2) is d/dt Ca_SR in component casr (mM).
C RATES(6) is d/dt d in component ical (dimensionless).
C RATES(7) is d/dt f in component ical (dimensionless).
C RATES(8) is d/dt fCa in component ical (dimensionless).
C RATES(9) is d/dt d in component icat (dimensionless).
C RATES(10) is d/dt f in component icat (dimensionless).
C RATES(11) is d/dt Xf in component ifunny (dimensionless).
C RATES(12) is d/dt Xr1 in component ikr (dimensionless).
C RATES(13) is d/dt Xr2 in component ikr (dimensionless).
C RATES(14) is d/dt Xs in component iks (dimensionless).
C RATES(15) is d/dt h in component ina (dimensionless).
C RATES(16) is d/dt j in component ina (dimensionless).
C RATES(17) is d/dt m in component ina (dimensionless).
C RATES(18) is d/dt I in component irel (dimensionless).
C RATES(19) is d/dt O in component irel (dimensionless).
C RATES(20) is d/dt R in component irel (dimensionless).
C RATES(21) is d/dt r in component ito (dimensionless).
C RATES(22) is d/dt s in component ito (dimensionless).
C RATES(3) is d/dt Ki in component ki (mM).
C RATES(5) is d/dt V in component membrane (mV).
C RATES(4) is d/dt Nai in component nai (mM).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 0.06
STATES(1) = 2.19191642424964390e-04
CONSTS(2) = 60.0
CONSTS(3) = 9.64853415000000041e+01
CONSTS(4) = 0.0006
STATES(2) = 3.35086796732326109e-01
CONSTS(5) = 0.3
CONSTS(6) = 0.0
CONSTS(7) = 1.8
STATES(3) = 1.04748824394112106e+02
CONSTS(8) = 5.4
STATES(4) = 7.16928091250999167
CONSTS(9) = 140.0
CONSTS(10) = 1094.0
CONSTS(11) = 3960.0
CONSTS(12) = 16404.0
STATES(5) = -7.55966016388546791e+01
STATES(6) = 3.94925342652924281e-04
CONSTS(13) = 1.29662941897219994e+01
CONSTS(14) = 7.07914596471100044
CONSTS(15) = 4.49094155069999987e-02
STATES(7) = 1.70990105585540286e-01
CONSTS(16) = 5.12589825999999987e-04
CONSTS(17) = 1.93121122351431995e+03
CONSTS(18) = 5.73002749969900016
STATES(8) = 8.77798946134088598e-01
CONSTS(19) = 3.08027691378999990e-01
CONSTS(20) = 0.00054
CONSTS(21) = 2.7e-07
CONSTS(22) = 1.5e-08
CONSTS(23) = 1.2
CONSTS(24) = 2.0
CONSTS(25) = 1.65824694683000007
CONSTS(26) = 1.00462559171102995e+02
STATES(9) = 2.70195573471577175e-04
STATES(10) = 7.56032904368393432e-01
CONSTS(27) = 0.185
CONSTS(28) = 0.491
STATES(11) = 6.40338504912615469e-03
CONSTS(29) = 0.0435
CONSTS(30) = 5.78970000000000002e-07
CONSTS(31) = 2.00866502378844016e+04
CONSTS(32) = 1.02023528452800001e+01
CONSTS(33) = 2.39452913465299986e+01
CONSTS(34) = 1.33785777797606004e-01
CONSTS(35) = 4.77994972217041014e-01
CONSTS(36) = 2.72427558793486995e+01
CONSTS(37) = 4.92502331781412028
CONSTS(38) = 8.72223760006881932
CONSTS(39) = 5.66361974998243980e+01
STATES(12) = 3.09767485715433222e-01
CONSTS(40) = 5.74885237435000026e-03
CONSTS(41) = 1.36234926362576001e+01
CONSTS(42) = 4.76305711818360011e-02
STATES(13) = 4.50577185148518577e-01
CONSTS(43) = 1.24566405268270002e-02
CONSTS(44) = 3.73426331501040991e+01
CONSTS(45) = 2.20919642353902006e+01
CONSTS(46) = 0.218025
CONSTS(47) = 50.0
CONSTS(48) = 0.0
STATES(14) = 1.53788281650948710e-01
CONSTS(49) = 0.0077
CONSTS(50) = 1.16558447999999992e-03
CONSTS(51) = 6.67268386758935958e+04
CONSTS(52) = 2.80458908250000027e-01
CONSTS(53) = 4.74115000000000034e-06
CONSTS(54) = 9.72061340924100037
STATES(15) = 7.39543607812429227e-01
CONSTS(55) = 3.62659886399999999e-03
CONSTS(56) = 9.66329497711473959e+03
CONSTS(57) = 7.39550356461299963
STATES(16) = 1.24515982574504899e-01
CONSTS(58) = 5.12257182000000044e-04
STATES(17) = 2.97549962926413614e-02
CONSTS(59) = 1.08045846384818006e+02
CONSTS(60) = 1.31070157339409992e+01
CONSTS(61) = 2.32691436700000007e-03
CONSTS(62) = 1.67331502516000014e-01
CONSTS(63) = 9.51088724962000032e-01
CONSTS(64) = 3.19775803839999970e-02
CONSTS(65) = 1.38
CONSTS(66) = 87.5
CONSTS(67) = 0.1
CONSTS(68) = 1.0
CONSTS(69) = 40.0
CONSTS(70) = 0.0005
STATES(18) = 1.42153622323011597e-02
CONSTS(71) = 15.0
CONSTS(72) = 1.0
STATES(19) = 1.65045105312396393e-04
STATES(20) = 1.13120363433751106e-02
CONSTS(73) = 0.45
CONSTS(74) = 12.5
CONSTS(75) = 1.17833333333299997e-01
STATES(21) = 2.67597833344160611e-04
CONSTS(76) = 5.53614181712999975e-02
CONSTS(77) = 1.16842023429669002e+01
CONSTS(78) = 3.98918108037750008
STATES(22) = 7.46802810614006107e-01
CONSTS(79) = 3.44230944300000013e-04
CONSTS(80) = 1.86760536909694991e+02
CONSTS(81) = 8.18093387332270083
CONSTS(82) = 6.96758421171499998e-01
CONSTS(83) = 1.12244577239468999e+01
CONSTS(84) = 8.314472
CONSTS(85) = 310.0
CONSTS(86) = CONSTS(12)+CONSTS(10)
CONSTS(87) = 10.0000*1.20000
CONSTS(88) = ( CONSTS(84)*CONSTS(85))/CONSTS(3)
CONSTS(89) = 0.000592000*0.620000
CONSTS(90) = 0.000290000*1.50000
CONSTS(91) = CONSTS(3)/( CONSTS(84)*CONSTS(85))
CONSTS(92) = CONSTS(15)*CONSTS(13)
CONSTS(93) = - 6.90988
CONSTS(94) = - 49.5057
CONSTS(95) = CONSTS(17)*CONSTS(16)
CONSTS(96) = (CONSTS(20)+CONSTS(22))+CONSTS(21)
CONSTS(97) = CONSTS(28)/(CONSTS(28)+1.00000)
CONSTS(98) = - 14.5897
CONSTS(99) = CONSTS(31)*CONSTS(30)
CONSTS(100) = CONSTS(42)*CONSTS(40)
CONSTS(101) = - 7.06809
CONSTS(102) = - 25.9945
CONSTS(103) = CONSTS(44)*CONSTS(43)
CONSTS(104) = CONSTS(52)*CONSTS(50)
CONSTS(105) = - 18.8670
CONSTS(106) = 8.00000e-05*0.0200000
CONSTS(107) = - 19.8394
CONSTS(108) = CONSTS(56)*CONSTS(55)
CONSTS(109) = - 66.5838
CONSTS(110) = CONSTS(56)
CONSTS(111) = CONSTS(61)*CONSTS(59)
CONSTS(112) = - 7.91773
CONSTS(113) = 2.50000*1.10000
CONSTS(114) = 0.350000*2.00000
CONSTS(115) = 1000.00*1.10000
CONSTS(116) = 1.36200*1.81800
CONSTS(117) = 0.0250000*10.5000
CONSTS(118) = 54.0000*0.342500
CONSTS(119) = 0.00100000*0.557100
CONSTS(120) = 1.50000*0.142900
CONSTS(121) = CONSTS(78)*CONSTS(76)
CONSTS(122) = - 11.0471
CONSTS(123) = - 17.6345
CONSTS(124) = CONSTS(80)*CONSTS(79)
CONSTS(125) = 0.000250000*0.702000
CONSTS(126) = 0.000425000*0.260000
CONSTS(127) = CONSTS(11)*(CONSTS(10)/CONSTS(86))
CONSTS(128) = CONSTS(3)*CONSTS(91)
CONSTS(129) = 1.00000/(1.00000/CONSTS(14)+1.00000/CONSTS(93))
CONSTS(130) = 1.00000/(1.00000/CONSTS(94)+1.00000/CONSTS(18))
CONSTS(131) = CONSTS(20)/CONSTS(96)
CONSTS(132) = CONSTS(21)/CONSTS(96)
CONSTS(133) = CONSTS(22)/CONSTS(96)
CONSTS(134) = 1.00000/(1.00000/CONSTS(98)+1.00000/CONSTS(32))
CONSTS(135) = 1.00000/(1.00000/CONSTS(41)+1.00000/CONSTS(101))
CONSTS(136) = 1.00000/(1.00000/CONSTS(102)+1.00000/CONSTS(45))
CONSTS(137) = 1.00000/(1.00000/CONSTS(51)+1.00000/CONSTS(105))
CONSTS(138) = 1.00000/(1.00000/CONSTS(107)+1.00000/CONSTS(57))
CONSTS(139) = CONSTS(110)*CONSTS(58)
CONSTS(140) = CONSTS(57)
CONSTS(141) = 1.00000/(1.00000/CONSTS(60)+1.00000/CONSTS(112))
CONSTS(142) = 56320.0*11.4302
CONSTS(143) = 1.00000/(1.00000/CONSTS(77)+1.00000/CONSTS(122))
CONSTS(144) = 1.00000/(1.00000/CONSTS(123)+1.00000/CONSTS(81))
CONSTS(145) = CONSTS(11)*(CONSTS(12)/CONSTS(86))
CONSTS(146) = CONSTS(131)*CONSTS(19)
CONSTS(147) = CONSTS(132)*CONSTS(19)
CONSTS(148) = CONSTS(133)*CONSTS(19)
CONSTS(149) = 1.00000/(1.00000/CONSTS(109)+1.00000/CONSTS(140))
CONSTS(150) = - 3.00000
CONSTS(151) = CONSTS(6)*CONSTS(150)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(5)+26.3000)/- 6.00000))
ALGBRC(21) = 1.00000/( 1.06800*EXP((STATES(5)+26.3000)/30.0000)+ 1.06800*EXP((STATES(5)+26.3000)/- 30.0000))
RATES(9) = (ALGBRC(5) - STATES(9))/ALGBRC(21)
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(5)+61.7000)/5.60000))
ALGBRC(22) = 1.00000/( 0.0153000*EXP(- (STATES(5)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(5)+61.7000)/15.3800))
RATES(10) = (ALGBRC(6) - STATES(10))/ALGBRC(22)
ALGBRC(2) = CONSTS(13)*EXP(STATES(5)/CONSTS(14))
ALGBRC(18) = CONSTS(92)*EXP(STATES(5)/CONSTS(129))
ALGBRC(34) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(18))
ALGBRC(48) = 1.00000/(ALGBRC(2)+ALGBRC(18))+CONSTS(25)
RATES(6) = (ALGBRC(34) - STATES(6))/ALGBRC(48)
ALGBRC(3) = CONSTS(16)*EXP(STATES(5)/CONSTS(94))
ALGBRC(19) = CONSTS(95)*EXP(STATES(5)/CONSTS(130))
ALGBRC(35) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(19))
ALGBRC(49) = 1.00000/(ALGBRC(3)+ALGBRC(19))+CONSTS(26)
RATES(7) = (ALGBRC(35) - STATES(7))/ALGBRC(49)
ALGBRC(7) = CONSTS(30)*EXP(STATES(5)/CONSTS(98))
ALGBRC(23) = CONSTS(99)*EXP(STATES(5)/CONSTS(134))
ALGBRC(37) = ALGBRC(7)/(ALGBRC(7)+ALGBRC(23))
ALGBRC(51) = 1.00000/(ALGBRC(7)+ALGBRC(23))+CONSTS(33)
RATES(11) = (ALGBRC(37) - STATES(11))/ALGBRC(51)
ALGBRC(8) = CONSTS(40)*EXP(STATES(5)/CONSTS(41))
ALGBRC(24) = CONSTS(100)*EXP(STATES(5)/CONSTS(135))
ALGBRC(38) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(24))
ALGBRC(52) = 1.00000/(ALGBRC(8)+ALGBRC(24))+CONSTS(47)
RATES(12) = (ALGBRC(38) - STATES(12))/ALGBRC(52)
ALGBRC(9) = CONSTS(43)*EXP(STATES(5)/CONSTS(102))
ALGBRC(25) = CONSTS(103)*EXP(STATES(5)/CONSTS(136))
ALGBRC(39) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(25))
ALGBRC(53) = 1.00000/(ALGBRC(9)+ALGBRC(25))+CONSTS(48)
RATES(13) = (ALGBRC(39) - STATES(13))/ALGBRC(53)
ALGBRC(10) = CONSTS(50)*EXP(STATES(5)/CONSTS(51))
ALGBRC(26) = CONSTS(104)*EXP(STATES(5)/CONSTS(137))
ALGBRC(40) = ALGBRC(10)/(ALGBRC(10)+ALGBRC(26))
ALGBRC(54) = 1.00000/(ALGBRC(10)+ALGBRC(26))+CONSTS(53)
RATES(14) = (ALGBRC(40) - STATES(14))/ALGBRC(54)
ALGBRC(11) = CONSTS(55)*EXP(STATES(5)/CONSTS(107))
ALGBRC(27) = CONSTS(108)*EXP(STATES(5)/CONSTS(138))
ALGBRC(41) = ALGBRC(11)/(ALGBRC(11)+ALGBRC(27))
ALGBRC(55) = 1.00000/(ALGBRC(11)+ALGBRC(27))+CONSTS(62)
RATES(15) = (ALGBRC(41) - STATES(15))/ALGBRC(55)
ALGBRC(12) = CONSTS(58)*EXP(STATES(5)/CONSTS(109))
ALGBRC(28) = CONSTS(139)*EXP(STATES(5)/CONSTS(149))
ALGBRC(42) = ALGBRC(12)/(ALGBRC(12)+ALGBRC(28))
ALGBRC(56) = 1.00000/(ALGBRC(12)+ALGBRC(28))+CONSTS(63)
RATES(16) = (ALGBRC(42) - STATES(16))/ALGBRC(56)
ALGBRC(13) = CONSTS(59)*EXP(STATES(5)/CONSTS(60))
ALGBRC(29) = CONSTS(111)*EXP(STATES(5)/CONSTS(141))
ALGBRC(43) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(29))
ALGBRC(57) = 1.00000/(ALGBRC(13)+ALGBRC(29))+CONSTS(64)
RATES(17) = (ALGBRC(43) - STATES(17))/ALGBRC(57)
ALGBRC(14) = ((1.00000 - STATES(20)) - STATES(19)) - STATES(18)
ALGBRC(30) = CONSTS(71) - (CONSTS(71) - CONSTS(72))/(1.00000+CONSTS(73)/STATES(2) ** 2.50000)
ALGBRC(44) = CONSTS(118)*ALGBRC(30)
ALGBRC(58) = CONSTS(142)/ALGBRC(30)
RATES(18) = (( ( ALGBRC(44)*STATES(1))*STATES(19) - CONSTS(119)*STATES(18)) - CONSTS(120)*STATES(18))+ ( ALGBRC(58)*STATES(1) ** 2.00000)*ALGBRC(14)
RATES(19) = (( ( ALGBRC(58)*STATES(1) ** 2.00000)*STATES(20) - CONSTS(120)*STATES(19)) - ( ALGBRC(44)*STATES(1))*STATES(19))+ CONSTS(119)*STATES(18)
RATES(20) = (( CONSTS(119)*ALGBRC(14) - ( ALGBRC(44)*STATES(1))*STATES(20)) - ( ALGBRC(58)*STATES(1) ** 2.00000)*STATES(20))+ CONSTS(120)*STATES(19)
ALGBRC(15) = CONSTS(76)*EXP(STATES(5)/CONSTS(77))
ALGBRC(31) = CONSTS(121)*EXP(STATES(5)/CONSTS(143))
ALGBRC(45) = ALGBRC(15)/(ALGBRC(15)+ALGBRC(31))
ALGBRC(59) = 1.00000/(ALGBRC(15)+ALGBRC(31))+CONSTS(82)
RATES(21) = (ALGBRC(45) - STATES(21))/ALGBRC(59)
ALGBRC(16) = CONSTS(79)*EXP(STATES(5)/CONSTS(123))
ALGBRC(32) = CONSTS(124)*EXP(STATES(5)/CONSTS(144))
ALGBRC(46) = ALGBRC(16)/(ALGBRC(16)+ALGBRC(32))
ALGBRC(60) = 1.00000/(ALGBRC(16)+ALGBRC(32))+CONSTS(83)
RATES(22) = (ALGBRC(46) - STATES(22))/ALGBRC(60)
ALGBRC(4) = 1.00000/(1.00000+( CONSTS(23)*STATES(1))/0.000325000 ** 8.00000)
ALGBRC(20) = 0.100000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000500000)/0.000100000))
ALGBRC(36) = 0.200000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000750000)/0.000800000))
ALGBRC(50) = (((ALGBRC(4)+ALGBRC(20))+ALGBRC(36))+0.230000)/1.46000
ALGBRC(62) = TERNRY(ALGBRC(50).GT.STATES(8).AND.STATES(5).GT.- 60.0000, 0.00000, 1.00000)
RATES(8) = ( ALGBRC(62)*(ALGBRC(50) - STATES(8)))/CONSTS(24)
ALGBRC(84) = ( CONSTS(115)*( ( EXP( ( CONSTS(114)*STATES(5))*CONSTS(91))*STATES(4) ** 3.00000)*CONSTS(7) - ( ( EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))*CONSTS(9) ** 3.00000)*STATES(1))*CONSTS(113)))/( ( (CONSTS(66) ** 3.00000+CONSTS(9) ** 3.00000)*(CONSTS(65)+CONSTS(7)))*(1.00000+ CONSTS(67)*EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))))
ALGBRC(61) = CONSTS(88)*log(CONSTS(9)/STATES(4))
ALGBRC(64) = CONSTS(90)*(STATES(5) - ALGBRC(61))
ALGBRC(69) = ( ( ( CONSTS(148)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(4))*EXP( STATES(5)*CONSTS(91)) - 0.750000*CONSTS(9)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
ALGBRC(70) = ( ( ALGBRC(69)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(74) = ( ( CONSTS(97)*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(61))
ALGBRC(83) = ( ( ( CONSTS(54)*STATES(17) ** 3.00000)*STATES(15))*STATES(16))*(STATES(5) - ALGBRC(61))
ALGBRC(85) = ( ( CONSTS(116)*CONSTS(8))*STATES(4))/( ( (CONSTS(8)+CONSTS(68))*(STATES(4)+CONSTS(69)))*((1.00000+ 0.124500*EXP( ( - 0.100000*STATES(5))*CONSTS(91)))+ 0.0353000*EXP( - STATES(5)*CONSTS(91))))
RATES(4) = (- CONSTS(2)/( CONSTS(3)*CONSTS(145)))*((((((ALGBRC(83))+ALGBRC(64))+ALGBRC(74))+ 3.00000*ALGBRC(85))+ 3.00000*ALGBRC(84))+ALGBRC(70))
ALGBRC(67) = ( ( ( CONSTS(147)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(3))*EXP( STATES(5)*CONSTS(91)) - 0.750000*CONSTS(8)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
ALGBRC(68) = ( ( ALGBRC(67)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(47) = CONSTS(88)*log(CONSTS(8)/STATES(3))
ALGBRC(73) = ( ( (1.00000 - CONSTS(97))*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(47))
ALGBRC(76) = CONSTS(35)*EXP((STATES(5)+CONSTS(37))/CONSTS(36))
ALGBRC(77) = 1.00000*EXP((STATES(5)+CONSTS(39))/CONSTS(38))
ALGBRC(78) = ALGBRC(76)/(ALGBRC(76)+ALGBRC(77))
ALGBRC(79) = ( ( CONSTS(34)* (CONSTS(8)/5.40000) ** (1.0 / 2))*ALGBRC(78))*(STATES(5) - ALGBRC(47))
ALGBRC(80) = ( ( ( CONSTS(46)* (CONSTS(8)/5.40000) ** (1.0 / 2))*STATES(12))*STATES(13))*(STATES(5) - ALGBRC(47))
ALGBRC(81) = ( CONSTS(49)*STATES(14) ** 2.00000)*(STATES(5) - ALGBRC(47))
ALGBRC(88) = ( ( CONSTS(75)*STATES(21))*STATES(22))*(STATES(5) - ALGBRC(47))
RATES(3) = (- CONSTS(2)/( CONSTS(3)*CONSTS(145)))*(((((((ALGBRC(79))+ALGBRC(88))+ALGBRC(80))+ALGBRC(81))+ALGBRC(73)) - 2.00000*ALGBRC(85))+ALGBRC(68))
ALGBRC(33) = ( 0.500000*CONSTS(88))*log(CONSTS(7)/STATES(1))
ALGBRC(72) = ( ( CONSTS(27)*STATES(9))*STATES(10))*(STATES(5) - ALGBRC(33))
ALGBRC(86) = ( CONSTS(117)*STATES(1))/(STATES(1)+CONSTS(70))
ALGBRC(63) = CONSTS(89)*(STATES(5) - ALGBRC(33))
ALGBRC(65) = ( ( ( ( CONSTS(146)*4.00000)*STATES(5))*CONSTS(128))*( ( 0.341000*STATES(1))*EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 0.341000*CONSTS(7)))/(EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 1.00000)
ALGBRC(66) = ( ( ALGBRC(65)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(71) = (ALGBRC(66)+ALGBRC(70))+ALGBRC(68)
ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
RATES(5) = - ((((((((((((((ALGBRC(79))+ALGBRC(88))+ALGBRC(80))+ALGBRC(81))+ALGBRC(71))+ALGBRC(72))+ALGBRC(85))+ALGBRC(83))+ALGBRC(84))+ALGBRC(86))+ALGBRC(75))+ALGBRC(64))+ALGBRC(63))+CONSTS(151))
ALGBRC(1) = 1.00000/(1.00000+( CONSTS(1)*CONSTS(4))/STATES(1)+CONSTS(4) ** 2.00000)
ALGBRC(82) = (STATES(2) - STATES(1))*CONSTS(106)
ALGBRC(87) = ( ( CONSTS(74)*STATES(19))*(STATES(2) - STATES(1)))*(CONSTS(127)/CONSTS(145))
ALGBRC(89) = CONSTS(126)/(1.00000+CONSTS(125) ** 2.00000/STATES(1) ** 2.00000)
RATES(1) = ALGBRC(1)*(((- ALGBRC(89)+ALGBRC(82))+ALGBRC(87)) - (CONSTS(2)/( ( 2.00000*CONSTS(145))*CONSTS(3)))*(((((ALGBRC(66))+ALGBRC(72))+ALGBRC(63))+ALGBRC(86)) - 2.00000*ALGBRC(84)))
ALGBRC(17) = 1.00000/(1.00000+( CONSTS(87)*CONSTS(5))/STATES(2)+CONSTS(5) ** 2.00000)
RATES(2) = (( ALGBRC(17)*CONSTS(145))/CONSTS(127))*((ALGBRC(89) - ALGBRC(87)) - ALGBRC(82))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(5) = 1.00000/(1.00000+EXP((STATES(5)+26.3000)/- 6.00000))
ALGBRC(21) = 1.00000/( 1.06800*EXP((STATES(5)+26.3000)/30.0000)+ 1.06800*EXP((STATES(5)+26.3000)/- 30.0000))
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(5)+61.7000)/5.60000))
ALGBRC(22) = 1.00000/( 0.0153000*EXP(- (STATES(5)+61.7000)/83.3000)+ 0.0150000*EXP((STATES(5)+61.7000)/15.3800))
ALGBRC(2) = CONSTS(13)*EXP(STATES(5)/CONSTS(14))
ALGBRC(18) = CONSTS(92)*EXP(STATES(5)/CONSTS(129))
ALGBRC(34) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(18))
ALGBRC(48) = 1.00000/(ALGBRC(2)+ALGBRC(18))+CONSTS(25)
ALGBRC(3) = CONSTS(16)*EXP(STATES(5)/CONSTS(94))
ALGBRC(19) = CONSTS(95)*EXP(STATES(5)/CONSTS(130))
ALGBRC(35) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(19))
ALGBRC(49) = 1.00000/(ALGBRC(3)+ALGBRC(19))+CONSTS(26)
ALGBRC(7) = CONSTS(30)*EXP(STATES(5)/CONSTS(98))
ALGBRC(23) = CONSTS(99)*EXP(STATES(5)/CONSTS(134))
ALGBRC(37) = ALGBRC(7)/(ALGBRC(7)+ALGBRC(23))
ALGBRC(51) = 1.00000/(ALGBRC(7)+ALGBRC(23))+CONSTS(33)
ALGBRC(8) = CONSTS(40)*EXP(STATES(5)/CONSTS(41))
ALGBRC(24) = CONSTS(100)*EXP(STATES(5)/CONSTS(135))
ALGBRC(38) = ALGBRC(8)/(ALGBRC(8)+ALGBRC(24))
ALGBRC(52) = 1.00000/(ALGBRC(8)+ALGBRC(24))+CONSTS(47)
ALGBRC(9) = CONSTS(43)*EXP(STATES(5)/CONSTS(102))
ALGBRC(25) = CONSTS(103)*EXP(STATES(5)/CONSTS(136))
ALGBRC(39) = ALGBRC(9)/(ALGBRC(9)+ALGBRC(25))
ALGBRC(53) = 1.00000/(ALGBRC(9)+ALGBRC(25))+CONSTS(48)
ALGBRC(10) = CONSTS(50)*EXP(STATES(5)/CONSTS(51))
ALGBRC(26) = CONSTS(104)*EXP(STATES(5)/CONSTS(137))
ALGBRC(40) = ALGBRC(10)/(ALGBRC(10)+ALGBRC(26))
ALGBRC(54) = 1.00000/(ALGBRC(10)+ALGBRC(26))+CONSTS(53)
ALGBRC(11) = CONSTS(55)*EXP(STATES(5)/CONSTS(107))
ALGBRC(27) = CONSTS(108)*EXP(STATES(5)/CONSTS(138))
ALGBRC(41) = ALGBRC(11)/(ALGBRC(11)+ALGBRC(27))
ALGBRC(55) = 1.00000/(ALGBRC(11)+ALGBRC(27))+CONSTS(62)
ALGBRC(12) = CONSTS(58)*EXP(STATES(5)/CONSTS(109))
ALGBRC(28) = CONSTS(139)*EXP(STATES(5)/CONSTS(149))
ALGBRC(42) = ALGBRC(12)/(ALGBRC(12)+ALGBRC(28))
ALGBRC(56) = 1.00000/(ALGBRC(12)+ALGBRC(28))+CONSTS(63)
ALGBRC(13) = CONSTS(59)*EXP(STATES(5)/CONSTS(60))
ALGBRC(29) = CONSTS(111)*EXP(STATES(5)/CONSTS(141))
ALGBRC(43) = ALGBRC(13)/(ALGBRC(13)+ALGBRC(29))
ALGBRC(57) = 1.00000/(ALGBRC(13)+ALGBRC(29))+CONSTS(64)
ALGBRC(14) = ((1.00000 - STATES(20)) - STATES(19)) - STATES(18)
ALGBRC(30) = CONSTS(71) - (CONSTS(71) - CONSTS(72))/(1.00000+CONSTS(73)/STATES(2) ** 2.50000)
ALGBRC(44) = CONSTS(118)*ALGBRC(30)
ALGBRC(58) = CONSTS(142)/ALGBRC(30)
ALGBRC(15) = CONSTS(76)*EXP(STATES(5)/CONSTS(77))
ALGBRC(31) = CONSTS(121)*EXP(STATES(5)/CONSTS(143))
ALGBRC(45) = ALGBRC(15)/(ALGBRC(15)+ALGBRC(31))
ALGBRC(59) = 1.00000/(ALGBRC(15)+ALGBRC(31))+CONSTS(82)
ALGBRC(16) = CONSTS(79)*EXP(STATES(5)/CONSTS(123))
ALGBRC(32) = CONSTS(124)*EXP(STATES(5)/CONSTS(144))
ALGBRC(46) = ALGBRC(16)/(ALGBRC(16)+ALGBRC(32))
ALGBRC(60) = 1.00000/(ALGBRC(16)+ALGBRC(32))+CONSTS(83)
ALGBRC(4) = 1.00000/(1.00000+( CONSTS(23)*STATES(1))/0.000325000 ** 8.00000)
ALGBRC(20) = 0.100000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000500000)/0.000100000))
ALGBRC(36) = 0.200000/(1.00000+EXP(( CONSTS(23)*STATES(1) - 0.000750000)/0.000800000))
ALGBRC(50) = (((ALGBRC(4)+ALGBRC(20))+ALGBRC(36))+0.230000)/1.46000
ALGBRC(62) = TERNRY(ALGBRC(50).GT.STATES(8).AND.STATES(5).GT.- 60.0000, 0.00000, 1.00000)
ALGBRC(84) = ( CONSTS(115)*( ( EXP( ( CONSTS(114)*STATES(5))*CONSTS(91))*STATES(4) ** 3.00000)*CONSTS(7) - ( ( EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))*CONSTS(9) ** 3.00000)*STATES(1))*CONSTS(113)))/( ( (CONSTS(66) ** 3.00000+CONSTS(9) ** 3.00000)*(CONSTS(65)+CONSTS(7)))*(1.00000+ CONSTS(67)*EXP( ( (CONSTS(114) - 1.00000)*STATES(5))*CONSTS(91))))
ALGBRC(61) = CONSTS(88)*log(CONSTS(9)/STATES(4))
ALGBRC(64) = CONSTS(90)*(STATES(5) - ALGBRC(61))
ALGBRC(69) = ( ( ( CONSTS(148)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(4))*EXP( STATES(5)*CONSTS(91)) - 0.750000*CONSTS(9)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
ALGBRC(70) = ( ( ALGBRC(69)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(74) = ( ( CONSTS(97)*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(61))
ALGBRC(83) = ( ( ( CONSTS(54)*STATES(17) ** 3.00000)*STATES(15))*STATES(16))*(STATES(5) - ALGBRC(61))
ALGBRC(85) = ( ( CONSTS(116)*CONSTS(8))*STATES(4))/( ( (CONSTS(8)+CONSTS(68))*(STATES(4)+CONSTS(69)))*((1.00000+ 0.124500*EXP( ( - 0.100000*STATES(5))*CONSTS(91)))+ 0.0353000*EXP( - STATES(5)*CONSTS(91))))
ALGBRC(67) = ( ( ( CONSTS(147)*STATES(5))*CONSTS(128))*( ( 0.750000*STATES(3))*EXP( STATES(5)*CONSTS(91)) - 0.750000*CONSTS(8)))/(EXP( STATES(5)*CONSTS(91)) - 1.00000)
ALGBRC(68) = ( ( ALGBRC(67)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(47) = CONSTS(88)*log(CONSTS(8)/STATES(3))
ALGBRC(73) = ( ( (1.00000 - CONSTS(97))*CONSTS(29))*STATES(11))*(STATES(5) - ALGBRC(47))
ALGBRC(76) = CONSTS(35)*EXP((STATES(5)+CONSTS(37))/CONSTS(36))
ALGBRC(77) = 1.00000*EXP((STATES(5)+CONSTS(39))/CONSTS(38))
ALGBRC(78) = ALGBRC(76)/(ALGBRC(76)+ALGBRC(77))
ALGBRC(79) = ( ( CONSTS(34)* (CONSTS(8)/5.40000) ** (1.0 / 2))*ALGBRC(78))*(STATES(5) - ALGBRC(47))
ALGBRC(80) = ( ( ( CONSTS(46)* (CONSTS(8)/5.40000) ** (1.0 / 2))*STATES(12))*STATES(13))*(STATES(5) - ALGBRC(47))
ALGBRC(81) = ( CONSTS(49)*STATES(14) ** 2.00000)*(STATES(5) - ALGBRC(47))
ALGBRC(88) = ( ( CONSTS(75)*STATES(21))*STATES(22))*(STATES(5) - ALGBRC(47))
ALGBRC(33) = ( 0.500000*CONSTS(88))*log(CONSTS(7)/STATES(1))
ALGBRC(72) = ( ( CONSTS(27)*STATES(9))*STATES(10))*(STATES(5) - ALGBRC(33))
ALGBRC(86) = ( CONSTS(117)*STATES(1))/(STATES(1)+CONSTS(70))
ALGBRC(63) = CONSTS(89)*(STATES(5) - ALGBRC(33))
ALGBRC(65) = ( ( ( ( CONSTS(146)*4.00000)*STATES(5))*CONSTS(128))*( ( 0.341000*STATES(1))*EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 0.341000*CONSTS(7)))/(EXP( ( 2.00000*STATES(5))*CONSTS(91)) - 1.00000)
ALGBRC(66) = ( ( ALGBRC(65)*STATES(6))*STATES(7))*STATES(8)
ALGBRC(71) = (ALGBRC(66)+ALGBRC(70))+ALGBRC(68)
ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
ALGBRC(1) = 1.00000/(1.00000+( CONSTS(1)*CONSTS(4))/STATES(1)+CONSTS(4) ** 2.00000)
ALGBRC(82) = (STATES(2) - STATES(1))*CONSTS(106)
ALGBRC(87) = ( ( CONSTS(74)*STATES(19))*(STATES(2) - STATES(1)))*(CONSTS(127)/CONSTS(145))
ALGBRC(89) = CONSTS(126)/(1.00000+CONSTS(125) ** 2.00000/STATES(1) ** 2.00000)
ALGBRC(17) = 1.00000/(1.00000+( CONSTS(87)*CONSTS(5))/STATES(2)+CONSTS(5) ** 2.00000)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END