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 108 entries in the algebraic variable array.
C There are a total of 55 entries in each of the rate and state variable arrays.
C There are a total of 130 entries in the constant variable array.
C
C
C VOI is time in component environment (second).
C CONSTS(1) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(2) is T in component model_parameters (kelvin).
C CONSTS(3) is F in component model_parameters (coulomb_per_mole).
C CONSTS(4) is Na_e in component model_parameters (millimolar).
C CONSTS(5) is Ca_e in component model_parameters (millimolar).
C CONSTS(6) is K_e in component model_parameters (millimolar).
C CONSTS(7) is ATP_i in component model_parameters (millimolar).
C CONSTS(128) is Vt in component model_parameters (cm3).
C CONSTS(85) is Vd in component model_parameters (cm3).
C CONSTS(87) is Vmyo in component model_parameters (cm3).
C CONSTS(81) is Sms in component model_parameters (cm2).
C CONSTS(91) is Smt in component model_parameters (cm2).
C CONSTS(130) is Cms in component model_parameters (microF).
C CONSTS(129) is Cmt in component model_parameters (microF).
C CONSTS(90) is Rst in component model_parameters (ohm).
C CONSTS(88) is VSRup in component model_parameters (cm3).
C CONSTS(86) is VSRrel in component model_parameters (cm3).
C CONSTS(84) is Vc in component model_parameters (cm3).
C CONSTS(89) is pt in component model_parameters (dimensionless).
C CONSTS(8) is Rot in component model_parameters (ohm_cm).
C CONSTS(9) is rt in component model_parameters (cm).
C CONSTS(10) is Lt in component model_parameters (cm).
C CONSTS(11) is rc in component model_parameters (cm).
C CONSTS(12) is Lc in component model_parameters (cm).
C CONSTS(13) is ptcm in component model_parameters (per_cm2).
C ALGBRC(15) is i_circ in component i_circ (microA).
C STATES(1) is Vm_s in component Vm_s (millivolt).
C STATES(2) is Vm_t in component Vm_t (millivolt).
C CONSTS(92) is g_Na_s in component membrane_permeabilities (milliS).
C CONSTS(93) is g_Naps_s in component membrane_permeabilities (milliS).
C CONSTS(94) is q_Kr_s in component membrane_permeabilities (milliS).
C CONSTS(95) is q_Ks_s in component membrane_permeabilities (milliS).
C CONSTS(96) is g_Kp_s in component membrane_permeabilities (milliS).
C CONSTS(97) is g_Kto_s in component membrane_permeabilities (milliS).
C CONSTS(98) is g_K1_s in component membrane_permeabilities (milliS).
C CONSTS(99) is g_KNa_s in component membrane_permeabilities (milliS).
C CONSTS(100) is g_KATP_s in component membrane_permeabilities (milliS).
C CONSTS(101) is g_Nab_s in component membrane_permeabilities (milliS).
C CONSTS(102) is g_Cab_s in component membrane_permeabilities (milliS).
C CONSTS(103) is i_pCa_max_s in component membrane_permeabilities (microA).
C CONSTS(104) is i_NaK_max_s in component membrane_permeabilities (microA).
C CONSTS(105) is i_NaCa_max_s in component membrane_permeabilities (microA).
C CONSTS(106) is P_CaL_s in component membrane_permeabilities (litre_per_second).
C CONSTS(107) is P_KL_s in component membrane_permeabilities (litre_per_second).
C CONSTS(108) is P_nsNa_s in component membrane_permeabilities (litre_per_second).
C CONSTS(109) is P_nsK_s in component membrane_permeabilities (litre_per_second).
C CONSTS(110) is g_Na_t in component membrane_permeabilities (milliS).
C CONSTS(111) is g_Naps_t in component membrane_permeabilities (milliS).
C CONSTS(112) is q_Kr_t in component membrane_permeabilities (milliS).
C CONSTS(113) is q_Ks_t in component membrane_permeabilities (milliS).
C CONSTS(114) is g_Kp_t in component membrane_permeabilities (milliS).
C CONSTS(115) is g_Kto_t in component membrane_permeabilities (milliS).
C CONSTS(116) is g_K1_t in component membrane_permeabilities (milliS).
C CONSTS(117) is g_KNa_t in component membrane_permeabilities (milliS).
C CONSTS(118) is g_KATP_t in component membrane_permeabilities (milliS).
C CONSTS(119) is g_Nab_t in component membrane_permeabilities (milliS).
C CONSTS(120) is g_Cab_t in component membrane_permeabilities (milliS).
C CONSTS(121) is i_pCa_max_t in component membrane_permeabilities (microA).
C CONSTS(122) is i_NaK_max_t in component membrane_permeabilities (microA).
C CONSTS(123) is i_NaCa_max_t in component membrane_permeabilities (microA).
C CONSTS(124) is P_CaL_t in component membrane_permeabilities (litre_per_second).
C CONSTS(125) is P_KL_t in component membrane_permeabilities (litre_per_second).
C CONSTS(126) is P_nsNa_t in component membrane_permeabilities (litre_per_second).
C CONSTS(127) is P_nsK_t in component membrane_permeabilities (litre_per_second).
C CONSTS(14) is g_Na in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(15) is fNat in component membrane_permeabilities (dimensionless).
C CONSTS(16) is g_Naps in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(17) is fNapst in component membrane_permeabilities (dimensionless).
C CONSTS(18) is q_Kr in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(19) is fKrt in component membrane_permeabilities (dimensionless).
C CONSTS(20) is q_Ks in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(21) is fKst in component membrane_permeabilities (dimensionless).
C CONSTS(22) is g_Kp in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(23) is fKpt in component membrane_permeabilities (dimensionless).
C CONSTS(24) is g_Kto in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(25) is fKtot in component membrane_permeabilities (dimensionless).
C CONSTS(26) is g_K1 in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(27) is fK1t in component membrane_permeabilities (dimensionless).
C CONSTS(28) is g_KNa in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(29) is fKNat in component membrane_permeabilities (dimensionless).
C CONSTS(30) is g_KATP in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(31) is fKATPt in component membrane_permeabilities (dimensionless).
C CONSTS(32) is g_Nab in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(33) is fNabt in component membrane_permeabilities (dimensionless).
C CONSTS(34) is g_Cab in component membrane_permeabilities (milliS_per_cm2).
C CONSTS(35) is fCabt in component membrane_permeabilities (dimensionless).
C CONSTS(36) is P_CaL in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(37) is fCaLt in component membrane_permeabilities (dimensionless).
C CONSTS(38) is P_KL in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(39) is P_nsNa in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(40) is fnsNat in component membrane_permeabilities (dimensionless).
C CONSTS(41) is P_nsK in component membrane_permeabilities (litre_per_second_cm2).
C CONSTS(42) is fnsKt in component membrane_permeabilities (dimensionless).
C CONSTS(43) is i_NaCa_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(44) is fNaCat in component membrane_permeabilities (dimensionless).
C CONSTS(45) is i_NaK_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(46) is fNaKt in component membrane_permeabilities (dimensionless).
C CONSTS(47) is i_pCa_max in component membrane_permeabilities (microA_per_cm2).
C CONSTS(48) is fpCat in component membrane_permeabilities (dimensionless).
C ALGBRC(30) is i_Kext in component i_Kext (microA).
C CONSTS(49) is stim_Period in component i_Kext (second).
C ALGBRC(36) is i_Na_s in component i_Na_s (microA).
C ALGBRC(33) is E_Na_s in component i_Na_s (millivolt).
C STATES(3) is Na_i in component ion_concentrations (millimolar).
C STATES(4) is K_i in component ion_concentrations (millimolar).
C STATES(5) is m in component i_Na_s_m_gate (dimensionless).
C STATES(6) is h in component i_Na_s_h_gate (dimensionless).
C ALGBRC(1) is m_infinity in component i_Na_s_m_gate (dimensionless).
C ALGBRC(16) is tau_m in component i_Na_s_m_gate (second).
C ALGBRC(2) is h_infinity in component i_Na_s_h_gate (dimensionless).
C ALGBRC(17) is tau_h in component i_Na_s_h_gate (second).
C ALGBRC(40) is i_Na_t in component i_Na_t (microA).
C ALGBRC(39) is E_Na_t in component i_Na_t (millivolt).
C STATES(7) is Na_t in component ion_concentrations (millimolar).
C STATES(8) is K_t in component ion_concentrations (millimolar).
C STATES(9) is m in component i_Na_t_m_gate (dimensionless).
C STATES(10) is h in component i_Na_t_h_gate (dimensionless).
C ALGBRC(3) is m_infinity in component i_Na_t_m_gate (dimensionless).
C ALGBRC(18) is tau_m in component i_Na_t_m_gate (second).
C ALGBRC(4) is h_infinity in component i_Na_t_h_gate (dimensionless).
C ALGBRC(19) is tau_h in component i_Na_t_h_gate (second).
C ALGBRC(41) is i_Naps_s in component i_Naps_s (microA).
C ALGBRC(42) is i_Naps_t in component i_Naps_t (microA).
C ALGBRC(43) is i_CaL_s in component i_CaL_s (microA).
C ALGBRC(45) is i_KL_s in component i_CaL_s (microA).
C ALGBRC(5) is alfas in component i_CaL_s (per_second).
C ALGBRC(20) is betas in component i_CaL_s (per_second).
C ALGBRC(31) is gama in component i_CaL_s (per_second).
C CONSTS(50) is omega in component i_CaL_s (per_second).
C CONSTS(51) is a in component i_CaL_s (dimensionless).
C CONSTS(52) is b in component i_CaL_s (dimensionless).
C CONSTS(53) is f in component i_CaL_s (per_second).
C CONSTS(54) is g in component i_CaL_s (per_second).
C CONSTS(55) is f2 in component i_CaL_s (per_second).
C CONSTS(56) is g2 in component i_CaL_s (per_second).
C ALGBRC(34) is alfa2s in component i_CaL_s (per_second).
C ALGBRC(37) is beta2s in component i_CaL_s (per_second).
C STATES(11) is Cst in component i_CaL_s (dimensionless).
C STATES(12) is C1 in component i_CaL_s (dimensionless).
C STATES(13) is C2 in component i_CaL_s (dimensionless).
C STATES(14) is C3 in component i_CaL_s (dimensionless).
C STATES(15) is C4 in component i_CaL_s (dimensionless).
C STATES(16) is Co in component i_CaL_s (dimensionless).
C STATES(17) is Ccast in component i_CaL_s (dimensionless).
C STATES(18) is Cca1 in component i_CaL_s (dimensionless).
C STATES(19) is Cca2 in component i_CaL_s (dimensionless).
C STATES(20) is Cca3 in component i_CaL_s (dimensionless).
C STATES(21) is Cca4 in component i_CaL_s (dimensionless).
C STATES(22) is Ccao in component i_CaL_s (dimensionless).
C STATES(23) is Ca_ss in component ion_concentrations (millimolar).
C ALGBRC(44) is i_CaL_t in component i_CaL_t (microA).
C STATES(24) is y in component i_CaL_s_y_gate (dimensionless).
C ALGBRC(6) is y_infinity in component i_CaL_s_y_gate (dimensionless).
C ALGBRC(21) is tau_y in component i_CaL_s_y_gate (second).
C ALGBRC(46) is i_KL_t in component i_CaL_t (microA).
C ALGBRC(22) is alfat in component i_CaL_t (per_second).
C ALGBRC(32) is betat in component i_CaL_t (per_second).
C ALGBRC(7) is gama in component i_CaL_t (per_second).
C CONSTS(57) is omega in component i_CaL_t (per_second).
C CONSTS(58) is a in component i_CaL_t (dimensionless).
C CONSTS(59) is b in component i_CaL_t (dimensionless).
C CONSTS(60) is f in component i_CaL_t (per_second).
C CONSTS(61) is g in component i_CaL_t (per_second).
C CONSTS(62) is f2 in component i_CaL_t (per_second).
C CONSTS(63) is g2 in component i_CaL_t (per_second).
C ALGBRC(35) is alfa2t in component i_CaL_t (per_second).
C ALGBRC(38) is beta2t in component i_CaL_t (per_second).
C STATES(25) is TCst in component i_CaL_t (dimensionless).
C STATES(26) is TC1 in component i_CaL_t (dimensionless).
C STATES(27) is TC2 in component i_CaL_t (dimensionless).
C STATES(28) is TC3 in component i_CaL_t (dimensionless).
C STATES(29) is TC4 in component i_CaL_t (dimensionless).
C STATES(30) is TCo in component i_CaL_t (dimensionless).
C STATES(31) is TCcast in component i_CaL_t (dimensionless).
C STATES(32) is TCca1 in component i_CaL_t (dimensionless).
C STATES(33) is TCca2 in component i_CaL_t (dimensionless).
C STATES(34) is TCca3 in component i_CaL_t (dimensionless).
C STATES(35) is TCca4 in component i_CaL_t (dimensionless).
C STATES(36) is TCcao in component i_CaL_t (dimensionless).
C STATES(37) is Ca_t in component ion_concentrations (millimolar).
C STATES(38) is y in component i_CaL_t_y_gate (dimensionless).
C ALGBRC(8) is y_infinity in component i_CaL_t_y_gate (dimensionless).
C ALGBRC(23) is tau_y in component i_CaL_t_y_gate (second).
C ALGBRC(49) is i_Kr_s in component i_Kr_s (microA).
C ALGBRC(47) is E_Kr_s in component i_Kr_s (millivolt).
C STATES(39) is xr in component i_Kr_s_xr_gate (dimensionless).
C ALGBRC(48) is xri in component i_Kr_s_xri_gate (dimensionless).
C ALGBRC(9) is xr_infinity in component i_Kr_s_xr_gate (dimensionless).
C ALGBRC(24) is tau_xr in component i_Kr_s_xr_gate (second).
C ALGBRC(52) is i_Kr_t in component i_Kr_t (microA).
C ALGBRC(50) is E_Kr_t in component i_Kr_t (millivolt).
C STATES(40) is xr in component i_Kr_t_xr_gate (dimensionless).
C ALGBRC(51) is xri in component i_Kr_t_xri_gate (dimensionless).
C ALGBRC(10) is xr_infinity in component i_Kr_t_xr_gate (dimensionless).
C ALGBRC(25) is tau_xr in component i_Kr_t_xr_gate (second).
C ALGBRC(54) is i_Ks_s in component i_Ks_s (microA).
C ALGBRC(53) is E_Ks_s in component i_Ks_s (millivolt).
C CONSTS(64) is PRNaK in component i_Ks_s (dimensionless).
C STATES(41) is Ca_i in component ion_concentrations (millimolar).
C STATES(42) is xs in component i_Ks_s_xs_gate (dimensionless).
C ALGBRC(11) is xs_infinity in component i_Ks_s_xs_gate (dimensionless).
C ALGBRC(26) is tau_xs in component i_Ks_s_xs_gate (second).
C ALGBRC(56) is i_Ks_t in component i_Ks_t (microA).
C ALGBRC(55) is E_Ks_t in component i_Ks_t (millivolt).
C CONSTS(65) is PRNaK in component i_Ks_t (dimensionless).
C STATES(43) is xs in component i_Ks_t_xs_gate (dimensionless).
C ALGBRC(12) is xs_infinity in component i_Ks_t_xs_gate (dimensionless).
C ALGBRC(27) is tau_xs in component i_Ks_t_xs_gate (second).
C ALGBRC(60) is i_K1_s in component i_K1_s (microA).
C ALGBRC(57) is E_K1_s in component i_K1_s (millivolt).
C ALGBRC(58) is aK1s in component i_K1_s (dimensionless).
C ALGBRC(59) is bK1s in component i_K1_s (dimensionless).
C ALGBRC(64) is i_K1_t in component i_K1_t (microA).
C ALGBRC(61) is E_K1_t in component i_K1_t (millivolt).
C ALGBRC(62) is aK1t in component i_K1_t (dimensionless).
C ALGBRC(63) is bK1t in component i_K1_t (dimensionless).
C ALGBRC(67) is i_Kp_s in component i_Kp_s (microA).
C ALGBRC(66) is kps in component i_Kp_s (dimensionless).
C ALGBRC(65) is E_Kp_s in component i_Kp_s (millivolt).
C ALGBRC(70) is i_Kp_t in component i_Kp_t (microA).
C ALGBRC(69) is kpt in component i_Kp_t (dimensionless).
C ALGBRC(68) is E_Kp_t in component i_Kp_t (millivolt).
C ALGBRC(72) is i_Kto_s in component i_Kto_s (microA).
C CONSTS(66) is fr2 in component i_Kto_s (dimensionless).
C ALGBRC(71) is E_Kto_s in component i_Kto_s (millivolt).
C STATES(44) is r2 in component i_Kto_s_r2_gate (dimensionless).
C STATES(45) is r3 in component i_Kto_s_r3_gate (dimensionless).
C CONSTS(67) is tau_r2 in component i_Kto_s_r2_gate (second).
C ALGBRC(13) is r2_infinity in component i_Kto_s_r2_gate (dimensionless).
C CONSTS(68) is tau_r3 in component i_Kto_s_r3_gate (second).
C ALGBRC(28) is r3_infinity in component i_Kto_s_r3_gate (dimensionless).
C ALGBRC(74) is i_Kto_t in component i_Kto_t (microA).
C CONSTS(69) is fr2 in component i_Kto_t (dimensionless).
C ALGBRC(73) is E_Kto_t in component i_Kto_t (millivolt).
C STATES(46) is r2 in component i_Kto_t_r2_gate (dimensionless).
C STATES(47) is r3 in component i_Kto_t_r3_gate (dimensionless).
C CONSTS(70) is tau_r2 in component i_Kto_t_r2_gate (second).
C ALGBRC(14) is r2_infinity in component i_Kto_t_r2_gate (dimensionless).
C ALGBRC(29) is r3_infinity in component i_Kto_t_r3_gate (dimensionless).
C CONSTS(71) is tau_r3 in component i_Kto_t_r3_gate (second).
C ALGBRC(76) is i_KNa_s in component i_KNa_s (microA).
C ALGBRC(75) is E_KNa_s in component i_KNa_s (millivolt).
C ALGBRC(78) is i_KNa_t in component i_KNa_t (microA).
C ALGBRC(77) is E_KNa_t in component i_KNa_t (millivolt).
C ALGBRC(79) is i_nsNa_s in component i_nsNa_s (microA).
C ALGBRC(80) is i_nsNa_t in component i_nsNa_t (microA).
C ALGBRC(81) is i_nsK_s in component i_nsK_s (microA).
C ALGBRC(82) is i_nsK_t in component i_nsK_t (microA).
C ALGBRC(83) is i_Nab_s in component i_Nab_s (microA).
C ALGBRC(84) is i_Nab_t in component i_Nab_t (microA).
C ALGBRC(86) is i_Cab_s in component i_Cab_s (microA).
C ALGBRC(85) is E_Ca_s in component i_Cab_s (millivolt).
C ALGBRC(88) is i_Cab_t in component i_Cab_t (microA).
C ALGBRC(87) is E_Ca_t in component i_Cab_t (millivolt).
C ALGBRC(89) is i_NaCa_s in component i_NaCa_s (microA).
C ALGBRC(90) is i_NaCa_t in component i_NaCa_t (microA).
C ALGBRC(91) is i_NaK_s in component i_NaK_s (microA).
C ALGBRC(92) is i_NaK_t in component i_NaK_t (microA).
C ALGBRC(94) is i_pCa_s in component i_pCa_s (microA).
C ALGBRC(95) is i_pCa_t in component i_pCa_t (microA).
C ALGBRC(99) is i_KATP_s in component i_KATP_s (microA).
C ALGBRC(96) is E_KATP_s in component i_KATP_s (millivolt).
C ALGBRC(103) is i_KATP_t in component i_KATP_t (microA).
C ALGBRC(101) is E_KATP_t in component i_KATP_t (millivolt).
C ALGBRC(93) is JteNa in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(97) is JteCa in component t_tubular_ion_fluxes (millimolar_per_second).
C ALGBRC(105) is JteK in component t_tubular_ion_fluxes (millimolar_per_second).
C CONSTS(72) is tau_Na in component t_tubular_ion_fluxes (second).
C CONSTS(73) is tau_Ca in component t_tubular_ion_fluxes (second).
C CONSTS(74) is tau_K in component t_tubular_ion_fluxes (second).
C ALGBRC(98) is JCaSRup in component JCaSRup (millimolar_per_second).
C ALGBRC(100) is JCaSRleak in component JCaSRleak (millimolar_per_second).
C STATES(48) is CaSRup in component CaSRup (millimolar).
C ALGBRC(102) is Jtr in component Jtr (millimolar_per_second).
C CONSTS(75) is tau_tr in component Jtr (second).
C STATES(49) is CaSRrel in component CaSRrel (millimolar).
C ALGBRC(104) is JCaSRrel in component JCaSRrel (millimolar_per_second).
C CONSTS(82) is kap in component JCaSRrel (per_millimolar4_per_second).
C CONSTS(76) is kam in component JCaSRrel (per_second).
C CONSTS(83) is kbp in component JCaSRrel (per_millimolar3_per_second).
C CONSTS(77) is kbm in component JCaSRrel (per_second).
C CONSTS(78) is kcp in component JCaSRrel (per_second).
C CONSTS(79) is kcm in component JCaSRrel (per_second).
C STATES(50) is F1 in component JCaSRrel (dimensionless).
C STATES(51) is F2 in component JCaSRrel (dimensionless).
C STATES(52) is F3 in component JCaSRrel (dimensionless).
C STATES(53) is F4 in component JCaSRrel (dimensionless).
C ALGBRC(106) is JCad in component JCad (millimolar_per_second).
C CONSTS(80) is tau_d in component JCad (second).
C STATES(54) is BTRH in component ion_concentrations (millimolar).
C STATES(55) is BTRL in component ion_concentrations (millimolar).
C ALGBRC(107) is dBTRH in component ion_concentrations (millimolar_per_second).
C ALGBRC(108) is dBTRL in component ion_concentrations (millimolar_per_second).
C RATES(5) is d/dt m in component i_Na_s_m_gate (dimensionless).
C RATES(6) is d/dt h in component i_Na_s_h_gate (dimensionless).
C RATES(9) is d/dt m in component i_Na_t_m_gate (dimensionless).
C RATES(10) is d/dt h in component i_Na_t_h_gate (dimensionless).
C RATES(11) is d/dt Cst in component i_CaL_s (dimensionless).
C RATES(12) is d/dt C1 in component i_CaL_s (dimensionless).
C RATES(13) is d/dt C2 in component i_CaL_s (dimensionless).
C RATES(14) is d/dt C3 in component i_CaL_s (dimensionless).
C RATES(15) is d/dt C4 in component i_CaL_s (dimensionless).
C RATES(16) is d/dt Co in component i_CaL_s (dimensionless).
C RATES(17) is d/dt Ccast in component i_CaL_s (dimensionless).
C RATES(18) is d/dt Cca1 in component i_CaL_s (dimensionless).
C RATES(19) is d/dt Cca2 in component i_CaL_s (dimensionless).
C RATES(20) is d/dt Cca3 in component i_CaL_s (dimensionless).
C RATES(21) is d/dt Cca4 in component i_CaL_s (dimensionless).
C RATES(22) is d/dt Ccao in component i_CaL_s (dimensionless).
C RATES(24) is d/dt y in component i_CaL_s_y_gate (dimensionless).
C RATES(25) is d/dt TCst in component i_CaL_t (dimensionless).
C RATES(26) is d/dt TC1 in component i_CaL_t (dimensionless).
C RATES(27) is d/dt TC2 in component i_CaL_t (dimensionless).
C RATES(28) is d/dt TC3 in component i_CaL_t (dimensionless).
C RATES(29) is d/dt TC4 in component i_CaL_t (dimensionless).
C RATES(30) is d/dt TCo in component i_CaL_t (dimensionless).
C RATES(31) is d/dt TCcast in component i_CaL_t (dimensionless).
C RATES(32) is d/dt TCca1 in component i_CaL_t (dimensionless).
C RATES(33) is d/dt TCca2 in component i_CaL_t (dimensionless).
C RATES(34) is d/dt TCca3 in component i_CaL_t (dimensionless).
C RATES(35) is d/dt TCca4 in component i_CaL_t (dimensionless).
C RATES(36) is d/dt TCcao in component i_CaL_t (dimensionless).
C RATES(38) is d/dt y in component i_CaL_t_y_gate (dimensionless).
C RATES(39) is d/dt xr in component i_Kr_s_xr_gate (dimensionless).
C RATES(40) is d/dt xr in component i_Kr_t_xr_gate (dimensionless).
C RATES(42) is d/dt xs in component i_Ks_s_xs_gate (dimensionless).
C RATES(43) is d/dt xs in component i_Ks_t_xs_gate (dimensionless).
C RATES(44) is d/dt r2 in component i_Kto_s_r2_gate (dimensionless).
C RATES(45) is d/dt r3 in component i_Kto_s_r3_gate (dimensionless).
C RATES(46) is d/dt r2 in component i_Kto_t_r2_gate (dimensionless).
C RATES(47) is d/dt r3 in component i_Kto_t_r3_gate (dimensionless).
C RATES(50) is d/dt F1 in component JCaSRrel (dimensionless).
C RATES(51) is d/dt F2 in component JCaSRrel (dimensionless).
C RATES(52) is d/dt F3 in component JCaSRrel (dimensionless).
C RATES(53) is d/dt F4 in component JCaSRrel (dimensionless).
C RATES(49) is d/dt CaSRrel in component CaSRrel (millimolar).
C RATES(48) is d/dt CaSRup in component CaSRup (millimolar).
C RATES(23) is d/dt Ca_ss in component ion_concentrations (millimolar).
C RATES(41) is d/dt Ca_i in component ion_concentrations (millimolar).
C RATES(54) is d/dt BTRH in component ion_concentrations (millimolar).
C RATES(55) is d/dt BTRL in component ion_concentrations (millimolar).
C RATES(3) is d/dt Na_i in component ion_concentrations (millimolar).
C RATES(4) is d/dt K_i in component ion_concentrations (millimolar).
C RATES(7) is d/dt Na_t in component ion_concentrations (millimolar).
C RATES(37) is d/dt Ca_t in component ion_concentrations (millimolar).
C RATES(8) is d/dt K_t in component ion_concentrations (millimolar).
C RATES(1) is d/dt Vm_s in component Vm_s (millivolt).
C RATES(2) is d/dt Vm_t in component Vm_t (millivolt).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
CONSTS(1) = 8310
CONSTS(2) = 310
CONSTS(3) = 96500
CONSTS(4) = 140
CONSTS(5) = 1.8
CONSTS(6) = 5.4
CONSTS(7) = 6.8
CONSTS(8) = 83.33
CONSTS(9) = 148.15e-7
CONSTS(10) = 5.93e-4
CONSTS(11) = 5.93e-4
CONSTS(12) = 0.013
CONSTS(13) = 21028875
STATES(1) = -8.5207812e1
STATES(2) = -8.5208034e1
CONSTS(14) = 30
CONSTS(15) = 0.64
CONSTS(16) = 0.0053
CONSTS(17) = 0.526
CONSTS(18) = 0.8
CONSTS(19) = 0.526
CONSTS(20) = 0.1975
CONSTS(21) = 0.526
CONSTS(22) = 0.006
CONSTS(23) = 0.526
CONSTS(24) = 0
CONSTS(25) = 0.526
CONSTS(26) = 0.75
CONSTS(27) = 0.8
CONSTS(28) = 0.12848
CONSTS(29) = 0.526
CONSTS(30) = 1
CONSTS(31) = 0.526
CONSTS(32) = 0.00075
CONSTS(33) = 0.526
CONSTS(34) = 0.0021
CONSTS(35) = 0.526
CONSTS(36) = 48e-4
CONSTS(37) = 0.64
CONSTS(38) = 48e-7
CONSTS(39) = 1.75e-7
CONSTS(40) = 0.526
CONSTS(41) = 1.75e-7
CONSTS(42) = 0.526
CONSTS(43) = 0.00025
CONSTS(44) = 0.526
CONSTS(45) = 1.5
CONSTS(46) = 0.526
CONSTS(47) = 1.15
CONSTS(48) = 0.2
CONSTS(49) = 1
STATES(3) = 1.1120279e1
STATES(4) = 1.3678926e2
STATES(5) = 3.4130493e-3
STATES(6) = 8.2699973e-1
STATES(7) = 1.3993239e2
STATES(8) = 5.4140321e0
STATES(9) = 3.4129472e-3
STATES(10) = 8.2700551e-1
CONSTS(50) = 10
CONSTS(51) = 2
CONSTS(52) = 2
CONSTS(53) = 300
CONSTS(54) = 2000
CONSTS(55) = 5
CONSTS(56) = 7000
STATES(11) = 9.9797984e-1
STATES(12) = 7.5722514e-5
STATES(13) = 2.1545646e-9
STATES(14) = 2.7463178e-14
STATES(15) = -5.7569102e-17
STATES(16) = -1.3496934e-17
STATES(17) = 1.9155065e-3
STATES(18) = 5.8135383e-7
STATES(19) = 6.615513e-11
STATES(20) = 3.3426555e-15
STATES(21) = 4.6640361e-19
STATES(22) = 3.7986638e-22
STATES(23) = 9.5977033e-5
STATES(24) = 7.8134852e-1
CONSTS(57) = 10
CONSTS(58) = 2
CONSTS(59) = 2
CONSTS(60) = 300
CONSTS(61) = 2000
CONSTS(62) = 5
CONSTS(63) = 7000
STATES(25) = 9.9798934e-1
STATES(26) = 7.5720257e-5
STATES(27) = 2.1544156e-9
STATES(28) = 2.7459475e-14
STATES(29) = -5.7367189e-17
STATES(30) = -1.3449454e-17
STATES(31) = 1.9155245e-3
STATES(32) = 5.8133643e-7
STATES(33) = 6.6150547e-11
STATES(34) = 3.3423167e-15
STATES(35) = 4.6165955e-19
STATES(36) = 3.75926e-22
STATES(37) = 1.8476402e0
STATES(38) = 7.813477e-1
STATES(39) = 2.0469344e-4
STATES(40) = 2.0469344e-4
CONSTS(64) = 0.01833
STATES(41) = 8.8787034e-5
STATES(42) = 6.1359896e-3
CONSTS(65) = 0.01833
STATES(43) = 6.1359896e-3
CONSTS(66) = 1
STATES(44) = 9.625701e-1
STATES(45) = 9.6250049e-1
CONSTS(67) = 0.02
CONSTS(68) = 0.08
CONSTS(69) = 1
STATES(46) = 9.625701e-1
STATES(47) = 9.6250049e-1
CONSTS(70) = 0.02
CONSTS(71) = 0.08
CONSTS(72) = 0.2
CONSTS(73) = 0.24
CONSTS(74) = 0.2
STATES(48) = 1.0200296e0
CONSTS(75) = 0.03448
STATES(49) = 1.0326252e0
CONSTS(76) = 576
CONSTS(77) = 1930
CONSTS(78) = 18
CONSTS(79) = 0.8
STATES(50) = 8.053082e-1
STATES(51) = 2.6639195e-4
STATES(52) = 4.9455459e-10
STATES(53) = 1.9442578e-1
CONSTS(80) = 0.003125
STATES(54) = 9.6700747e-1
STATES(55) = 8.1740868e-2
CONSTS(81) = 2.00000*3.14159*CONSTS(11) ** 2.00000+ 2.00000*3.14159*CONSTS(11)*CONSTS(12)
CONSTS(82) = 3.00000*1.21500e+13
CONSTS(83) = 4.05000e+09
CONSTS(84) = 3.14159*CONSTS(11) ** 2.00000*CONSTS(12)
CONSTS(85) = CONSTS(84)*0.680000*6.00000e-05
CONSTS(86) = CONSTS(84)*0.00316000
CONSTS(87) = CONSTS(84)*0.680000
CONSTS(88) = CONSTS(84)*0.0550000
CONSTS(89) = CONSTS(13)*2.00000*3.14159*CONSTS(11)*CONSTS(12)
CONSTS(90) = (( CONSTS(8)*CONSTS(10))/2.00000)/( 3.14159*CONSTS(9) ** 2.00000*CONSTS(89))
CONSTS(91) = 2.00000*3.14159*CONSTS(9)*CONSTS(10)*CONSTS(89)
CONSTS(92) = CONSTS(14)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(15))
CONSTS(93) = CONSTS(16)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(17))
CONSTS(94) = CONSTS(18)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(19))
CONSTS(95) = CONSTS(20)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(21))
CONSTS(96) = CONSTS(22)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(23))
CONSTS(97) = CONSTS(24)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(25))
CONSTS(98) = CONSTS(26)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(27))
CONSTS(99) = CONSTS(28)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(29))
CONSTS(100) = CONSTS(30)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(31))
CONSTS(101) = CONSTS(32)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(33))
CONSTS(102) = CONSTS(34)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(35))
CONSTS(103) = CONSTS(47)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(48))
CONSTS(104) = CONSTS(45)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(46))
CONSTS(105) = CONSTS(43)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(44))
CONSTS(106) = CONSTS(36)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(37))
CONSTS(107) = CONSTS(38)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(37))
CONSTS(108) = CONSTS(39)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(40))
CONSTS(109) = CONSTS(41)*(CONSTS(81)+CONSTS(91))*(1.00000 - CONSTS(42))
CONSTS(110) = CONSTS(14)*(CONSTS(81)+CONSTS(91))*CONSTS(15)
CONSTS(111) = CONSTS(16)*(CONSTS(81)+CONSTS(91))*CONSTS(17)
CONSTS(112) = CONSTS(18)*(CONSTS(81)+CONSTS(91))*CONSTS(19)
CONSTS(113) = CONSTS(20)*(CONSTS(81)+CONSTS(91))*CONSTS(21)
CONSTS(114) = CONSTS(22)*(CONSTS(81)+CONSTS(91))*CONSTS(23)
CONSTS(115) = CONSTS(24)*(CONSTS(81)+CONSTS(91))*CONSTS(25)
CONSTS(116) = CONSTS(26)*(CONSTS(81)+CONSTS(91))*CONSTS(27)
CONSTS(117) = CONSTS(28)*(CONSTS(81)+CONSTS(91))*CONSTS(29)
CONSTS(118) = CONSTS(30)*(CONSTS(81)+CONSTS(91))*CONSTS(31)
CONSTS(119) = CONSTS(32)*(CONSTS(81)+CONSTS(91))*CONSTS(33)
CONSTS(120) = CONSTS(34)*(CONSTS(81)+CONSTS(91))*CONSTS(35)
CONSTS(121) = CONSTS(47)*(CONSTS(81)+CONSTS(91))*CONSTS(48)
CONSTS(122) = CONSTS(45)*(CONSTS(81)+CONSTS(91))*CONSTS(46)
CONSTS(123) = CONSTS(43)*(CONSTS(81)+CONSTS(91))*CONSTS(44)
CONSTS(124) = CONSTS(36)*(CONSTS(81)+CONSTS(91))*CONSTS(37)
CONSTS(125) = CONSTS(38)*(CONSTS(81)+CONSTS(91))*CONSTS(37)
CONSTS(126) = CONSTS(39)*(CONSTS(81)+CONSTS(91))*CONSTS(40)
CONSTS(127) = CONSTS(41)*(CONSTS(81)+CONSTS(91))*CONSTS(42)
CONSTS(128) = 3.14159*CONSTS(9) ** 2.00000*CONSTS(10)*CONSTS(89)
CONSTS(129) = CONSTS(91)*1.00000
CONSTS(130) = CONSTS(81)*1.00000
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(16) = CONSTS(53)*STATES(15) - CONSTS(54)*STATES(16)
RATES(22) = CONSTS(55)*STATES(21) - CONSTS(56)*STATES(22)
RATES(30) = CONSTS(60)*STATES(29) - CONSTS(61)*STATES(30)
RATES(36) = CONSTS(62)*STATES(35) - CONSTS(63)*STATES(36)
RATES(50) = STATES(51)*CONSTS(76) - STATES(50)*CONSTS(82)*STATES(23) ** 4.00000
RATES(51) = ( STATES(50)*CONSTS(82)*STATES(23) ** 4.00000+ STATES(52)*CONSTS(77)+ STATES(53)*CONSTS(79)) - STATES(51)*(CONSTS(76)+ CONSTS(83)*STATES(23) ** 3.00000+CONSTS(78))
RATES(52) = STATES(51)*CONSTS(83)*STATES(23) ** 3.00000 - STATES(52)*CONSTS(77)
RATES(53) = STATES(51)*CONSTS(78) - STATES(53)*CONSTS(79)
ALGBRC(13) = 1.00600/(1.00000+EXP((STATES(1)+36.6900+10.0000)/12.4300))
RATES(44) = (ALGBRC(13) - STATES(44))/CONSTS(67)
ALGBRC(14) = 1.00600/(1.00000+EXP((STATES(2)+36.6900+10.0000)/12.4300))
RATES(46) = (ALGBRC(14) - STATES(46))/CONSTS(70)
ALGBRC(1) = 1.00000/(1.00000+EXP(((STATES(1)+52.2000) - 9.00000)/- 7.40000))
ALGBRC(16) = (0.00100000/( 101.600*EXP( STATES(1)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(1)))+0.000100000)/4.83000
RATES(5) = (ALGBRC(1) - STATES(5))/ALGBRC(16)
ALGBRC(2) = 1.00000/(1.00000+EXP(((STATES(1)+85.6000) - 9.00000)/5.50000))
ALGBRC(17) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(1))+ 6.53700*EXP( STATES(1)*0.0801600))+0.000500000)/4.83000
RATES(6) = (ALGBRC(2) - STATES(6))/ALGBRC(17)
ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(2)+52.2000) - 9.00000)/- 7.40000))
ALGBRC(18) = (0.00100000/( 101.600*EXP( STATES(2)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(2)))+0.000100000)/4.83000
RATES(9) = (ALGBRC(3) - STATES(9))/ALGBRC(18)
ALGBRC(4) = 1.00000/(1.00000+EXP(((STATES(2)+85.6000) - 9.00000)/5.50000))
ALGBRC(19) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(2))+ 6.53700*EXP( STATES(2)*0.0801600))+0.000500000)/4.83000
RATES(10) = (ALGBRC(4) - STATES(10))/ALGBRC(19)
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/6.00000))
ALGBRC(21) = 0.00100000/(0.0200000+ 0.0197000*EXP(- (STATES(1)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(1)+40.0000)/9.50000)*4.00000))
RATES(24) = (ALGBRC(6) - STATES(24))/ALGBRC(21)
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(2)+35.0000)/6.00000))
ALGBRC(23) = 0.00100000/(0.0200000+ 0.0197000*EXP(- (STATES(2)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(2)+40.0000)/9.50000)*4.00000))
RATES(38) = (ALGBRC(8) - STATES(38))/ALGBRC(23)
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+21.5000)/7.50000))
ALGBRC(24) = 0.00100000/(( 0.00138000*(STATES(1)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(1)+14.2000)))+( 0.000610000*(STATES(1)+38.9000))/(EXP( 0.145000*(STATES(1)+38.9000)) - 1.00000))
RATES(39) = (ALGBRC(9) - STATES(39))/ALGBRC(24)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(2)+21.5000)/7.50000))
ALGBRC(25) = 0.00100000/(( 0.00138000*(STATES(2)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(2)+14.2000)))+( 0.000610000*(STATES(2)+38.9000))/(EXP( 0.145000*(STATES(2)+38.9000)) - 1.00000))
RATES(40) = (ALGBRC(10) - STATES(40))/ALGBRC(25)
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(26) = 0.00100000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(EXP( 0.0687000*(STATES(1)+30.0000)) - 1.00000))
RATES(42) = (ALGBRC(11) - STATES(42))/ALGBRC(26)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(2) - 1.50000)/16.7000))
ALGBRC(27) = 0.00100000/(( 7.19000e-05*(STATES(2)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(2)+30.0000)))+( 0.000131000*(STATES(2)+30.0000))/(EXP( 0.0687000*(STATES(2)+30.0000)) - 1.00000))
RATES(43) = (ALGBRC(12) - STATES(43))/ALGBRC(27)
ALGBRC(28) = ALGBRC(13)
RATES(45) = (ALGBRC(28) - STATES(45))/CONSTS(68)
ALGBRC(29) = ALGBRC(14)
RATES(47) = (ALGBRC(29) - STATES(47))/CONSTS(71)
ALGBRC(5) = 400.000*EXP((STATES(1)+12.0000)/10.0000)
ALGBRC(20) = 50.0000*EXP(- (STATES(1)+12.0000)/13.0000)
ALGBRC(31) = 187.500*STATES(23)
RATES(11) = ( ALGBRC(20)*STATES(12)+ CONSTS(50)*STATES(17)) - ( 4.00000*ALGBRC(5)+ALGBRC(31))*STATES(11)
RATES(12) = ( 4.00000*ALGBRC(5)*STATES(11)+ 2.00000*ALGBRC(20)*STATES(13)+ (CONSTS(50)/CONSTS(52))*STATES(18)) - (ALGBRC(20)+ 3.00000*ALGBRC(5)+ ALGBRC(31)*CONSTS(51))*STATES(12)
RATES(13) = ( 3.00000*ALGBRC(5)*STATES(12)+ 3.00000*ALGBRC(20)*STATES(14)+ (CONSTS(50)/CONSTS(52) ** 2.00000)*STATES(19)) - ( 2.00000*ALGBRC(20)+ 2.00000*ALGBRC(5)+ ALGBRC(31)*CONSTS(51) ** 2.00000)*STATES(13)
RATES(14) = ( 2.00000*ALGBRC(5)*STATES(13)+ 4.00000*ALGBRC(20)*STATES(15)+ (CONSTS(50)/CONSTS(52) ** 3.00000)*STATES(20)) - ( 3.00000*ALGBRC(20)+ALGBRC(5)+ ALGBRC(31)*CONSTS(51) ** 3.00000)*STATES(14)
RATES(15) = ( ALGBRC(5)*STATES(14)+ CONSTS(54)*STATES(16)+ (CONSTS(50)/CONSTS(52) ** 4.00000)*STATES(21)) - ( 4.00000*ALGBRC(20)+CONSTS(53)+ ALGBRC(31)*CONSTS(51) ** 4.00000)*STATES(15)
ALGBRC(22) = 400.000*EXP((STATES(2)+12.0000)/10.0000)
ALGBRC(32) = 50.0000*EXP(- (STATES(2)+12.0000)/13.0000)
ALGBRC(7) = 187.500*STATES(23)
RATES(25) = ( ALGBRC(32)*STATES(26)+ CONSTS(57)*STATES(31)) - ( 4.00000*ALGBRC(22)+ALGBRC(7))*STATES(25)
RATES(26) = ( 4.00000*ALGBRC(22)*STATES(25)+ 2.00000*ALGBRC(32)*STATES(27)+ (CONSTS(57)/CONSTS(59))*STATES(32)) - (ALGBRC(32)+ 3.00000*ALGBRC(22)+ ALGBRC(7)*CONSTS(58))*STATES(26)
RATES(27) = ( 3.00000*ALGBRC(22)*STATES(26)+ 3.00000*ALGBRC(32)*STATES(28)+ (CONSTS(57)/CONSTS(59) ** 2.00000)*STATES(33)) - ( 2.00000*ALGBRC(32)+ 2.00000*ALGBRC(22)+ ALGBRC(7)*CONSTS(58) ** 2.00000)*STATES(27)
RATES(28) = ( 2.00000*ALGBRC(22)*STATES(27)+ 4.00000*ALGBRC(32)*STATES(29)+ (CONSTS(57)/CONSTS(59) ** 3.00000)*STATES(34)) - ( 3.00000*ALGBRC(32)+ALGBRC(22)+ ALGBRC(7)*CONSTS(58) ** 3.00000)*STATES(28)
RATES(29) = ( ALGBRC(22)*STATES(28)+ CONSTS(61)*STATES(30)+ (CONSTS(57)/CONSTS(59) ** 4.00000)*STATES(35)) - ( 4.00000*ALGBRC(32)+CONSTS(60)+ ALGBRC(7)*CONSTS(58) ** 4.00000)*STATES(29)
ALGBRC(34) = ALGBRC(5)*CONSTS(51)
ALGBRC(37) = ALGBRC(20)/CONSTS(52)
RATES(17) = ( ALGBRC(37)*STATES(18)+ ALGBRC(31)*STATES(11)) - ( 4.00000*ALGBRC(34)+CONSTS(50))*STATES(17)
RATES(18) = ( 4.00000*ALGBRC(34)*STATES(17)+ 2.00000*ALGBRC(37)*STATES(19)+ ALGBRC(31)*CONSTS(51)*STATES(12)) - (ALGBRC(37)+ 3.00000*ALGBRC(34)+CONSTS(50)/CONSTS(52))*STATES(18)
RATES(19) = ( 3.00000*ALGBRC(34)*STATES(18)+ 3.00000*ALGBRC(37)*STATES(20)+ ALGBRC(31)*CONSTS(51) ** 2.00000*STATES(13)) - ( 2.00000*ALGBRC(37)+ 2.00000*ALGBRC(34)+CONSTS(50)/CONSTS(52) ** 2.00000)*STATES(19)
RATES(20) = ( 2.00000*ALGBRC(34)*STATES(19)+ 4.00000*ALGBRC(37)*STATES(21)+ ALGBRC(31)*CONSTS(51) ** 3.00000*STATES(14)) - ( 3.00000*ALGBRC(37)+ALGBRC(34)+CONSTS(50)/CONSTS(52) ** 3.00000)*STATES(20)
RATES(21) = ( ALGBRC(34)*STATES(20)+ CONSTS(56)*STATES(22)+ ALGBRC(31)*CONSTS(51) ** 4.00000*STATES(15)) - ( 4.00000*ALGBRC(37)+CONSTS(55)+CONSTS(50)/CONSTS(52) ** 4.00000)*STATES(21)
ALGBRC(35) = ALGBRC(22)*CONSTS(58)
ALGBRC(38) = ALGBRC(32)/CONSTS(59)
RATES(31) = ( ALGBRC(38)*STATES(32)+ ALGBRC(7)*STATES(25)) - ( 4.00000*ALGBRC(35)+CONSTS(57))*STATES(31)
RATES(32) = ( 4.00000*ALGBRC(35)*STATES(31)+ 2.00000*ALGBRC(38)*STATES(33)+ ALGBRC(7)*CONSTS(58)*STATES(26)) - (ALGBRC(38)+ 3.00000*ALGBRC(35)+CONSTS(57)/CONSTS(59))*STATES(32)
RATES(33) = ( 3.00000*ALGBRC(35)*STATES(32)+ 3.00000*ALGBRC(38)*STATES(34)+ ALGBRC(7)*CONSTS(58) ** 2.00000*STATES(27)) - ( 2.00000*ALGBRC(38)+ 2.00000*ALGBRC(35)+CONSTS(57)/CONSTS(59) ** 2.00000)*STATES(33)
RATES(34) = ( 2.00000*ALGBRC(35)*STATES(33)+ 4.00000*ALGBRC(38)*STATES(35)+ ALGBRC(7)*CONSTS(58) ** 3.00000*STATES(28)) - ( 3.00000*ALGBRC(38)+ALGBRC(35)+CONSTS(57)/CONSTS(59) ** 3.00000)*STATES(34)
RATES(35) = ( ALGBRC(35)*STATES(34)+ CONSTS(63)*STATES(36)+ ALGBRC(7)*CONSTS(58) ** 4.00000*STATES(29)) - ( 4.00000*ALGBRC(38)+CONSTS(62)+CONSTS(57)/CONSTS(59) ** 4.00000)*STATES(35)
ALGBRC(36) = CONSTS(92)*STATES(5) ** 3.00000*STATES(6)*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(40) = CONSTS(110)*STATES(9) ** 3.00000*STATES(10)*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(41) = (CONSTS(93)/(1.00000+EXP((- 54.0000 - STATES(1))/8.00000)))*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(42) = (CONSTS(111)/(1.00000+EXP((- 54.0000 - STATES(2))/8.00000)))*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(79) = (( (( 1.00000*CONSTS(108)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*CONSTS(4)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(80) = (( (( 1.00000*CONSTS(126)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*STATES(7)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(33) = ( log(CONSTS(4)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(83) = CONSTS(101)*(STATES(1) - ALGBRC(33))
ALGBRC(39) = ( log(STATES(7)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(84) = CONSTS(119)*(STATES(2) - ALGBRC(39))
ALGBRC(89) = ( CONSTS(105)*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5) - CONSTS(4) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5)+ CONSTS(4) ** 3.00000*STATES(41)))
ALGBRC(90) = ( CONSTS(123)*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37) - STATES(7) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37)+ STATES(7) ** 3.00000*STATES(41)))
ALGBRC(91) = ( (( (CONSTS(104)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(CONSTS(4)/67.3000) - 1.00000)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*CONSTS(6))/(CONSTS(6)+1.50000)
ALGBRC(92) = ( (( (CONSTS(122)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(STATES(7)/67.3000) - 1.00000)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*STATES(8))/(STATES(8)+1.50000)
RATES(3) = - (ALGBRC(36)+ALGBRC(40)+ALGBRC(41)+ALGBRC(42)+ALGBRC(79)+ALGBRC(80)+ALGBRC(83)+ALGBRC(84)+ 3.00000*ALGBRC(89)+ 3.00000*ALGBRC(90)+ 3.00000*ALGBRC(91)+ 3.00000*ALGBRC(92))/( CONSTS(3)*CONSTS(87))
ALGBRC(93) = (( 1.00000*CONSTS(128))/CONSTS(72))*(STATES(7) - CONSTS(4))
RATES(7) = ((ALGBRC(40)+ALGBRC(42)+ALGBRC(80)+ALGBRC(84)+ 3.00000*ALGBRC(90)+ 3.00000*ALGBRC(92))/CONSTS(3) - 1.00000*ALGBRC(93))/CONSTS(128)
ALGBRC(44) = ( (( 1.00000*CONSTS(124)*4.00000*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.341000*STATES(37)))/(EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(87) = ( log(STATES(37)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
ALGBRC(88) = CONSTS(120)*(STATES(2) - ALGBRC(87))
ALGBRC(95) = ( CONSTS(121)*STATES(41))/(0.000500000+STATES(41))
ALGBRC(97) = (( 1.00000*CONSTS(128))/CONSTS(73))*(STATES(37) - CONSTS(5))
RATES(37) = (( - 2.00000*ALGBRC(90)+ALGBRC(44)+ALGBRC(88)+ALGBRC(95))/( 2.00000*CONSTS(3)) - 1.00000*ALGBRC(97))/CONSTS(128)
ALGBRC(15) = 1000.00*(STATES(2)/CONSTS(90) - STATES(1)/CONSTS(90))
ALGBRC(30) = TERNRY(VOI - INT(VOI/CONSTS(49))*CONSTS(49).GE.0.00000.AND.VOI - INT(VOI/CONSTS(49))*CONSTS(49).LE.0.00100000, 45.0000*(CONSTS(81)+CONSTS(91)), 0.00000)
ALGBRC(43) = ( (( 1.00000*CONSTS(106)*4.00000*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.341000*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(45) = ( (( (( 1.00000*CONSTS(107))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(47) = ( log(CONSTS(6)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(1)+9.00000)/22.4000))
ALGBRC(49) = CONSTS(94)*0.0261400* (CONSTS(6)/5.40000) ** (1.0 / 2)*STATES(39)*ALGBRC(48)*(STATES(1) - ALGBRC(47))
ALGBRC(53) = ( log((CONSTS(6)+ CONSTS(64)*CONSTS(4))/(STATES(4)+ CONSTS(64)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(54) = CONSTS(95)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(42) ** 2.00000*(STATES(1) - ALGBRC(53))
ALGBRC(57) = ALGBRC(47)
ALGBRC(58) = 1020.00/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(57)) - 59.2150)))
ALGBRC(59) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(57))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(57)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(57))+4.75300)))
ALGBRC(60) = (( CONSTS(98)* (CONSTS(6)/5.40000) ** (1.0 / 2)*ALGBRC(58))/(ALGBRC(58)+ALGBRC(59)))*(STATES(1) - ALGBRC(57))
ALGBRC(66) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/5.00000))
ALGBRC(65) = ALGBRC(47)
ALGBRC(67) = CONSTS(96)*ALGBRC(66)*(STATES(1) - ALGBRC(65))
ALGBRC(71) = ALGBRC(47)
ALGBRC(72) = (CONSTS(97)/( (1.00000+EXP((STATES(1)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(1) - 45.8000)/25.8700))))*( STATES(44)*CONSTS(66)+ STATES(45)*(1.00000 - CONSTS(66)))*(STATES(1) - ALGBRC(71))
ALGBRC(75) = ALGBRC(47)
ALGBRC(76) = (( CONSTS(99)*(0.800000 - 0.650000/(1.00000+EXP((STATES(1)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(1) - ALGBRC(75))
ALGBRC(81) = (( (( 1.00000*CONSTS(109)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(85) = ( log(CONSTS(5)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
ALGBRC(86) = CONSTS(102)*(STATES(1) - ALGBRC(85))
ALGBRC(94) = ( CONSTS(103)*STATES(41))/(0.000500000+STATES(41))
ALGBRC(96) = ALGBRC(47)
ALGBRC(99) = (( CONSTS(100)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*CONSTS(6)/4.00000 ** 0.240000*(STATES(1) - ALGBRC(96))
RATES(1) = ( 1000.00*(((((((((((((((((((ALGBRC(30)+ALGBRC(15)) - ALGBRC(36)) - ALGBRC(41)) - ALGBRC(43)) - ALGBRC(45)) - ALGBRC(54)) - ALGBRC(49)) - ALGBRC(60)) - ALGBRC(67)) - ALGBRC(76)) - ALGBRC(99)) - ALGBRC(79)) - ALGBRC(81)) - ALGBRC(83)) - ALGBRC(86)) - ALGBRC(89)) - ALGBRC(91)) - ALGBRC(94)) - ALGBRC(72)))/CONSTS(130)
ALGBRC(98) = ( 1.00000*1000.00*CONSTS(87)*0.00180000*STATES(41) ** 2.00000)/(STATES(41) ** 2.00000+0.000500000 ** 2.00000)
ALGBRC(100) = 1250.00*CONSTS(87)*5.80000e-05*(STATES(48) - STATES(41))
ALGBRC(102) = (( 1.00000*CONSTS(86))/CONSTS(75))*(STATES(48) - STATES(49))
RATES(48) = ((ALGBRC(98) - ALGBRC(100)) - ALGBRC(102))/( 1.00000*CONSTS(88))
ALGBRC(104) = 1800.00*CONSTS(86)*(STATES(51)+STATES(52))*(STATES(49) - STATES(23))
RATES(49) = ( (1.00000/(1.00000+( 15.0000*0.800000)/0.800000+STATES(49) ** 2.00000))*(ALGBRC(102) - ALGBRC(104)))/CONSTS(86)
ALGBRC(46) = ( (( (( 1.00000*CONSTS(125))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(50) = ( log(STATES(8)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(2)+9.00000)/22.4000))
ALGBRC(52) = CONSTS(112)*0.0261400* (STATES(8)/5.40000) ** (1.0 / 2)*STATES(40)*ALGBRC(51)*(STATES(2) - ALGBRC(50))
ALGBRC(55) = ( log((STATES(8)+ CONSTS(65)*STATES(7))/(STATES(4)+ CONSTS(65)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(56) = CONSTS(113)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(43) ** 2.00000*(STATES(2) - ALGBRC(55))
ALGBRC(61) = ALGBRC(50)
ALGBRC(62) = 1020.00/(1.00000+EXP( 0.238500*((STATES(2) - ALGBRC(61)) - 59.2150)))
ALGBRC(63) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(2) - ALGBRC(61))+5.47600))+EXP( 0.0617500*((STATES(2) - ALGBRC(61)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(2) - ALGBRC(61))+4.75300)))
ALGBRC(64) = (( CONSTS(116)* (STATES(8)/5.40000) ** (1.0 / 2)*ALGBRC(62))/(ALGBRC(62)+ALGBRC(63)))*(STATES(2) - ALGBRC(61))
ALGBRC(69) = 1.00000/(1.00000+EXP((20.0000 - STATES(2))/5.00000))
ALGBRC(68) = ALGBRC(50)
ALGBRC(70) = CONSTS(114)*ALGBRC(69)*(STATES(2) - ALGBRC(68))
ALGBRC(73) = ALGBRC(50)
ALGBRC(74) = (CONSTS(115)/( (1.00000+EXP((STATES(2)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(2) - 45.8000)/25.8700))))*( STATES(46)*CONSTS(69)+ STATES(47)*(1.00000 - CONSTS(69)))*(STATES(2) - ALGBRC(73))
ALGBRC(77) = ALGBRC(50)
ALGBRC(78) = (( CONSTS(117)*(0.800000 - 0.650000/(1.00000+EXP((STATES(2)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(2) - ALGBRC(77))
ALGBRC(82) = (( (( 1.00000*CONSTS(127)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(101) = ALGBRC(50)
ALGBRC(103) = (( CONSTS(118)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*STATES(8)/4.00000 ** 0.240000*(STATES(2) - ALGBRC(101))
RATES(4) = - (((( - 0.00000*ALGBRC(30)+ALGBRC(54)+ALGBRC(56)+ALGBRC(49)+ALGBRC(52)+ALGBRC(60)+ALGBRC(64)+ALGBRC(67)+ALGBRC(70)+ALGBRC(45)+ALGBRC(46)+ALGBRC(76)+ALGBRC(78)+ALGBRC(81)+ALGBRC(82)) - 2.00000*ALGBRC(91)) - 2.00000*ALGBRC(92))+ALGBRC(99)+ALGBRC(103)+ALGBRC(72)+ALGBRC(74))/( CONSTS(3)*CONSTS(87))
RATES(2) = ( 1000.00*((((((((((((((((((- ALGBRC(15) - ALGBRC(40)) - ALGBRC(42)) - ALGBRC(44)) - ALGBRC(46)) - ALGBRC(56)) - ALGBRC(52)) - ALGBRC(64)) - ALGBRC(70)) - ALGBRC(78)) - ALGBRC(103)) - ALGBRC(80)) - ALGBRC(82)) - ALGBRC(84)) - ALGBRC(88)) - ALGBRC(90)) - ALGBRC(92)) - ALGBRC(95)) - ALGBRC(74)))/CONSTS(129)
ALGBRC(106) = ( 1.00000*(STATES(23) - STATES(41))*CONSTS(87))/CONSTS(80)
RATES(23) = (1.00000/(1.00000+( 0.0500000*0.00238000)/0.00238000+STATES(23) ** 2.00000))*((- (ALGBRC(43)+ALGBRC(44))/( 2.00000*CONSTS(3)*CONSTS(85))+( 1.00000*ALGBRC(104))/CONSTS(85)) - ( 1.00000*ALGBRC(106))/CONSTS(85))
ALGBRC(105) = (( 1.00000*CONSTS(128))/CONSTS(74))*(STATES(8) - CONSTS(6))
RATES(8) = ((((ALGBRC(56)+ALGBRC(52)+ALGBRC(64)+ALGBRC(70)+ALGBRC(78)+ALGBRC(46)+ALGBRC(82)) - 2.00000*ALGBRC(92))+ALGBRC(103)+ALGBRC(74))/CONSTS(3) - 1.00000*ALGBRC(105))/CONSTS(128)
ALGBRC(107) = 20000.0*STATES(41)*(1.00000 - STATES(54)) - 0.0700000*STATES(54)
RATES(54) = ALGBRC(107)
ALGBRC(108) = 40000.0*STATES(41)*(1.00000 - STATES(55)) - STATES(55)*40.0000
RATES(41) = (1.00000/(1.00000+( 0.0500000*0.00238000)/0.00238000+STATES(41) ** 2.00000))*(((( - 2.00000*(ALGBRC(89)+ALGBRC(90))+ALGBRC(86)+ALGBRC(88)+ALGBRC(94)+ALGBRC(95))/( - 2.00000*CONSTS(3)*CONSTS(87))+( 1.00000*((ALGBRC(106)+ALGBRC(100)) - ALGBRC(98)))/CONSTS(87)) - ALGBRC(107)*0.140000) - ALGBRC(108)*0.0700000)
RATES(55) = ALGBRC(108)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(13) = 1.00600/(1.00000+EXP((STATES(1)+36.6900+10.0000)/12.4300))
ALGBRC(14) = 1.00600/(1.00000+EXP((STATES(2)+36.6900+10.0000)/12.4300))
ALGBRC(1) = 1.00000/(1.00000+EXP(((STATES(1)+52.2000) - 9.00000)/- 7.40000))
ALGBRC(16) = (0.00100000/( 101.600*EXP( STATES(1)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(1)))+0.000100000)/4.83000
ALGBRC(2) = 1.00000/(1.00000+EXP(((STATES(1)+85.6000) - 9.00000)/5.50000))
ALGBRC(17) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(1))+ 6.53700*EXP( STATES(1)*0.0801600))+0.000500000)/4.83000
ALGBRC(3) = 1.00000/(1.00000+EXP(((STATES(2)+52.2000) - 9.00000)/- 7.40000))
ALGBRC(18) = (0.00100000/( 101.600*EXP( STATES(2)*0.113500)+ 0.0226800*EXP( - 0.0717000*STATES(2)))+0.000100000)/4.83000
ALGBRC(4) = 1.00000/(1.00000+EXP(((STATES(2)+85.6000) - 9.00000)/5.50000))
ALGBRC(19) = (0.00100000/( 1.13810e-06*EXP( - 0.101700*STATES(2))+ 6.53700*EXP( STATES(2)*0.0801600))+0.000500000)/4.83000
ALGBRC(6) = 1.00000/(1.00000+EXP((STATES(1)+35.0000)/6.00000))
ALGBRC(21) = 0.00100000/(0.0200000+ 0.0197000*EXP(- (STATES(1)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(1)+40.0000)/9.50000)*4.00000))
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(2)+35.0000)/6.00000))
ALGBRC(23) = 0.00100000/(0.0200000+ 0.0197000*EXP(- (STATES(2)+10.0000)*0.0337000 ** 2.00000))+0.550000/(1.00000+EXP( ((STATES(2)+40.0000)/9.50000)*4.00000))
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+21.5000)/7.50000))
ALGBRC(24) = 0.00100000/(( 0.00138000*(STATES(1)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(1)+14.2000)))+( 0.000610000*(STATES(1)+38.9000))/(EXP( 0.145000*(STATES(1)+38.9000)) - 1.00000))
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(2)+21.5000)/7.50000))
ALGBRC(25) = 0.00100000/(( 0.00138000*(STATES(2)+14.2000))/(1.00000 - EXP( - 0.123000*(STATES(2)+14.2000)))+( 0.000610000*(STATES(2)+38.9000))/(EXP( 0.145000*(STATES(2)+38.9000)) - 1.00000))
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(26) = 0.00100000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(EXP( 0.0687000*(STATES(1)+30.0000)) - 1.00000))
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(2) - 1.50000)/16.7000))
ALGBRC(27) = 0.00100000/(( 7.19000e-05*(STATES(2)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(2)+30.0000)))+( 0.000131000*(STATES(2)+30.0000))/(EXP( 0.0687000*(STATES(2)+30.0000)) - 1.00000))
ALGBRC(28) = ALGBRC(13)
ALGBRC(29) = ALGBRC(14)
ALGBRC(5) = 400.000*EXP((STATES(1)+12.0000)/10.0000)
ALGBRC(20) = 50.0000*EXP(- (STATES(1)+12.0000)/13.0000)
ALGBRC(31) = 187.500*STATES(23)
ALGBRC(22) = 400.000*EXP((STATES(2)+12.0000)/10.0000)
ALGBRC(32) = 50.0000*EXP(- (STATES(2)+12.0000)/13.0000)
ALGBRC(7) = 187.500*STATES(23)
ALGBRC(34) = ALGBRC(5)*CONSTS(51)
ALGBRC(37) = ALGBRC(20)/CONSTS(52)
ALGBRC(35) = ALGBRC(22)*CONSTS(58)
ALGBRC(38) = ALGBRC(32)/CONSTS(59)
ALGBRC(36) = CONSTS(92)*STATES(5) ** 3.00000*STATES(6)*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(40) = CONSTS(110)*STATES(9) ** 3.00000*STATES(10)*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(41) = (CONSTS(93)/(1.00000+EXP((- 54.0000 - STATES(1))/8.00000)))*(STATES(1) - ( log((CONSTS(4)+ 0.120000*CONSTS(6))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(42) = (CONSTS(111)/(1.00000+EXP((- 54.0000 - STATES(2))/8.00000)))*(STATES(2) - ( log((STATES(7)+ 0.120000*STATES(8))/(STATES(3)+ 0.120000*STATES(4)))*CONSTS(1)*CONSTS(2))/CONSTS(3))
ALGBRC(79) = (( (( 1.00000*CONSTS(108)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*CONSTS(4)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(80) = (( (( 1.00000*CONSTS(126)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(3)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*STATES(7)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(33) = ( log(CONSTS(4)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(83) = CONSTS(101)*(STATES(1) - ALGBRC(33))
ALGBRC(39) = ( log(STATES(7)/STATES(3))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(84) = CONSTS(119)*(STATES(2) - ALGBRC(39))
ALGBRC(89) = ( CONSTS(105)*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5) - CONSTS(4) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*CONSTS(5)+ CONSTS(4) ** 3.00000*STATES(41)))
ALGBRC(90) = ( CONSTS(123)*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37) - STATES(7) ** 3.00000*STATES(41)))/(1.00000+ 0.000100000*EXP(( - 0.850000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*( EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))*STATES(3) ** 3.00000*STATES(37)+ STATES(7) ** 3.00000*STATES(41)))
ALGBRC(91) = ( (( (CONSTS(104)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(CONSTS(4)/67.3000) - 1.00000)*EXP(( - STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*CONSTS(6))/(CONSTS(6)+1.50000)
ALGBRC(92) = ( (( (CONSTS(122)/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))+ (( 0.0365000*1.00000)/7.00000)*(EXP(STATES(7)/67.3000) - 1.00000)*EXP(( - STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2)))))*1.00000)/(1.00000+10.0000/STATES(3) ** 1.50000))*STATES(8))/(STATES(8)+1.50000)
ALGBRC(93) = (( 1.00000*CONSTS(128))/CONSTS(72))*(STATES(7) - CONSTS(4))
ALGBRC(44) = ( (( 1.00000*CONSTS(124)*4.00000*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.341000*STATES(37)))/(EXP(( 2.00000*STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(87) = ( log(STATES(37)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
ALGBRC(88) = CONSTS(120)*(STATES(2) - ALGBRC(87))
ALGBRC(95) = ( CONSTS(121)*STATES(41))/(0.000500000+STATES(41))
ALGBRC(97) = (( 1.00000*CONSTS(128))/CONSTS(73))*(STATES(37) - CONSTS(5))
ALGBRC(15) = 1000.00*(STATES(2)/CONSTS(90) - STATES(1)/CONSTS(90))
ALGBRC(30) = TERNRY(VOI - INT(VOI/CONSTS(49))*CONSTS(49).GE.0.00000.AND.VOI - INT(VOI/CONSTS(49))*CONSTS(49).LE.0.00100000, 45.0000*(CONSTS(81)+CONSTS(91)), 0.00000)
ALGBRC(43) = ( (( 1.00000*CONSTS(106)*4.00000*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.00100000*EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.341000*CONSTS(5)))/(EXP(( 2.00000*STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(45) = ( (( (( 1.00000*CONSTS(107))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(16)+STATES(22))*STATES(24)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(47) = ( log(CONSTS(6)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(48) = 1.00000/(1.00000+EXP((STATES(1)+9.00000)/22.4000))
ALGBRC(49) = CONSTS(94)*0.0261400* (CONSTS(6)/5.40000) ** (1.0 / 2)*STATES(39)*ALGBRC(48)*(STATES(1) - ALGBRC(47))
ALGBRC(53) = ( log((CONSTS(6)+ CONSTS(64)*CONSTS(4))/(STATES(4)+ CONSTS(64)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(54) = CONSTS(95)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(42) ** 2.00000*(STATES(1) - ALGBRC(53))
ALGBRC(57) = ALGBRC(47)
ALGBRC(58) = 1020.00/(1.00000+EXP( 0.238500*((STATES(1) - ALGBRC(57)) - 59.2150)))
ALGBRC(59) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(1) - ALGBRC(57))+5.47600))+EXP( 0.0617500*((STATES(1) - ALGBRC(57)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - ALGBRC(57))+4.75300)))
ALGBRC(60) = (( CONSTS(98)* (CONSTS(6)/5.40000) ** (1.0 / 2)*ALGBRC(58))/(ALGBRC(58)+ALGBRC(59)))*(STATES(1) - ALGBRC(57))
ALGBRC(66) = 1.00000/(1.00000+EXP((20.0000 - STATES(1))/5.00000))
ALGBRC(65) = ALGBRC(47)
ALGBRC(67) = CONSTS(96)*ALGBRC(66)*(STATES(1) - ALGBRC(65))
ALGBRC(71) = ALGBRC(47)
ALGBRC(72) = (CONSTS(97)/( (1.00000+EXP((STATES(1)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(1) - 45.8000)/25.8700))))*( STATES(44)*CONSTS(66)+ STATES(45)*(1.00000 - CONSTS(66)))*(STATES(1) - ALGBRC(71))
ALGBRC(75) = ALGBRC(47)
ALGBRC(76) = (( CONSTS(99)*(0.800000 - 0.650000/(1.00000+EXP((STATES(1)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(1) - ALGBRC(75))
ALGBRC(81) = (( (( 1.00000*CONSTS(109)*STATES(1)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(85) = ( log(CONSTS(5)/STATES(41))*CONSTS(1)*CONSTS(2))/( 2.00000*CONSTS(3))
ALGBRC(86) = CONSTS(102)*(STATES(1) - ALGBRC(85))
ALGBRC(94) = ( CONSTS(103)*STATES(41))/(0.000500000+STATES(41))
ALGBRC(96) = ALGBRC(47)
ALGBRC(99) = (( CONSTS(100)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*CONSTS(6)/4.00000 ** 0.240000*(STATES(1) - ALGBRC(96))
ALGBRC(98) = ( 1.00000*1000.00*CONSTS(87)*0.00180000*STATES(41) ** 2.00000)/(STATES(41) ** 2.00000+0.000500000 ** 2.00000)
ALGBRC(100) = 1250.00*CONSTS(87)*5.80000e-05*(STATES(48) - STATES(41))
ALGBRC(102) = (( 1.00000*CONSTS(86))/CONSTS(75))*(STATES(48) - STATES(49))
ALGBRC(104) = 1800.00*CONSTS(86)*(STATES(51)+STATES(52))*(STATES(49) - STATES(23))
ALGBRC(46) = ( (( (( 1.00000*CONSTS(125))/(1.00000 - (ALGBRC(43)+ALGBRC(44))/( 0.458000*(CONSTS(81)+CONSTS(91)))))*(STATES(30)+STATES(36))*STATES(38)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000)
ALGBRC(50) = ( log(STATES(8)/STATES(4))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(51) = 1.00000/(1.00000+EXP((STATES(2)+9.00000)/22.4000))
ALGBRC(52) = CONSTS(112)*0.0261400* (STATES(8)/5.40000) ** (1.0 / 2)*STATES(40)*ALGBRC(51)*(STATES(2) - ALGBRC(50))
ALGBRC(55) = ( log((STATES(8)+ CONSTS(65)*STATES(7))/(STATES(4)+ CONSTS(65)*STATES(3)))*CONSTS(1)*CONSTS(2))/CONSTS(3)
ALGBRC(56) = CONSTS(113)*(0.230800+0.769200/(1.00000+EXP((- arbitrary_log( 1.00000*STATES(41), 10) - 4.20000)/0.600000)))*STATES(43) ** 2.00000*(STATES(2) - ALGBRC(55))
ALGBRC(61) = ALGBRC(50)
ALGBRC(62) = 1020.00/(1.00000+EXP( 0.238500*((STATES(2) - ALGBRC(61)) - 59.2150)))
ALGBRC(63) = ( 1000.00*( 0.491240*EXP( 0.0803200*((STATES(2) - ALGBRC(61))+5.47600))+EXP( 0.0617500*((STATES(2) - ALGBRC(61)) - 594.310))))/(1.00000+EXP( - 0.514300*((STATES(2) - ALGBRC(61))+4.75300)))
ALGBRC(64) = (( CONSTS(116)* (STATES(8)/5.40000) ** (1.0 / 2)*ALGBRC(62))/(ALGBRC(62)+ALGBRC(63)))*(STATES(2) - ALGBRC(61))
ALGBRC(69) = 1.00000/(1.00000+EXP((20.0000 - STATES(2))/5.00000))
ALGBRC(68) = ALGBRC(50)
ALGBRC(70) = CONSTS(114)*ALGBRC(69)*(STATES(2) - ALGBRC(68))
ALGBRC(73) = ALGBRC(50)
ALGBRC(74) = (CONSTS(115)/( (1.00000+EXP((STATES(2)+57.5300)/- 5.86300))*(1.00000+EXP((STATES(2) - 45.8000)/25.8700))))*( STATES(46)*CONSTS(69)+ STATES(47)*(1.00000 - CONSTS(69)))*(STATES(2) - ALGBRC(73))
ALGBRC(77) = ALGBRC(50)
ALGBRC(78) = (( CONSTS(117)*(0.800000 - 0.650000/(1.00000+EXP((STATES(2)+125.000)/15.0000)))*0.850000)/(1.00000+66.0000/STATES(3) ** 2.80000))*(STATES(2) - ALGBRC(77))
ALGBRC(82) = (( (( 1.00000*CONSTS(127)*STATES(2)*CONSTS(3) ** 2.00000)/( CONSTS(1)*CONSTS(2)))*( 0.750000*STATES(4)*EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 0.750000*STATES(8)))/(EXP(( STATES(2)*CONSTS(3))/( CONSTS(1)*CONSTS(2))) - 1.00000))/(1.00000+0.00250000/STATES(41) ** 3.00000)
ALGBRC(101) = ALGBRC(50)
ALGBRC(103) = (( CONSTS(118)*1.00000)/(1.00000+CONSTS(7)/0.114000 ** 2.00000))*STATES(8)/4.00000 ** 0.240000*(STATES(2) - ALGBRC(101))
ALGBRC(106) = ( 1.00000*(STATES(23) - STATES(41))*CONSTS(87))/CONSTS(80)
ALGBRC(105) = (( 1.00000*CONSTS(128))/CONSTS(74))*(STATES(8) - CONSTS(6))
ALGBRC(107) = 20000.0*STATES(41)*(1.00000 - STATES(54)) - 0.0700000*STATES(54)
ALGBRC(108) = 40000.0*STATES(41)*(1.00000 - STATES(55)) - STATES(55)*40.0000
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END