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 128 entries in the algebraic variable array.
C There are a total of 45 entries in each of the rate and state variable arrays.
C There are a total of 177 entries in the constant variable array.
C
C
C VOI is time in component environment (millisecond).
C STATES(1) is V in component cell (millivolt).
C ALGBRC(51) is i_Na in component INa (microA_per_microF).
C ALGBRC(54) is i_Nab in component INab (microA_per_microF).
C ALGBRC(58) is i_NaK in component INaK (microA_per_microF).
C ALGBRC(60) is i_Kr in component IKr (microA_per_microF).
C ALGBRC(68) is i_Ks in component IKs (microA_per_microF).
C ALGBRC(72) is i_tos in component Itos (microA_per_microF).
C ALGBRC(73) is i_tof in component Itof (microA_per_microF).
C ALGBRC(77) is i_K1 in component IK1 (microA_per_microF).
C ALGBRC(78) is i_Kp in component IKp (microA_per_microF).
C ALGBRC(97) is i_NaCa in component INaCa (microA_per_microF).
C ALGBRC(79) is i_Cl_Ca in component ICl_Ca (microA_per_microF).
C ALGBRC(80) is i_Clb in component IClb (microA_per_microF).
C ALGBRC(89) is i_CaL in component ICaL (microA_per_microF).
C ALGBRC(106) is i_Cab in component ICab (microA_per_microF).
C ALGBRC(103) is i_Cap in component ICap (microA_per_microF).
C ALGBRC(7) is i_Stim in component cell (microA_per_microF).
C CONSTS(1) is stim_start in component cell (millisecond).
C CONSTS(2) is stim_period in component cell (millisecond).
C CONSTS(3) is stim_duration in component cell (millisecond).
C CONSTS(4) is stim_amplitude in component cell (microA_per_microF).
C CONSTS(5) is Ko in component model_parameters (millimolar).
C CONSTS(6) is Nao in component model_parameters (millimolar).
C CONSTS(7) is Cao in component model_parameters (millimolar).
C CONSTS(8) is Clo in component model_parameters (millimolar).
C CONSTS(9) is Ki in component model_parameters (millimolar).
C CONSTS(10) is Mgi in component model_parameters (millimolar).
C CONSTS(11) is Mgo in component model_parameters (millimolar).
C CONSTS(12) is Cli in component model_parameters (millimolar).
C CONSTS(13) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(14) is T in component model_parameters (kelvin).
C CONSTS(15) is F in component model_parameters (coulomb_per_mole).
C CONSTS(170) is Cm in component model_parameters (farad).
C CONSTS(16) is Cm_per_area in component model_parameters (farad_per_cm2).
C CONSTS(17) is cell_length in component model_parameters (micrometre).
C CONSTS(18) is cell_radius in component model_parameters (micrometre).
C CONSTS(19) is jct_radius in component model_parameters (micrometre).
C CONSTS(20) is jct_depth in component model_parameters (micrometre).
C CONSTS(21) is SL_depth in component model_parameters (micrometre).
C CONSTS(171) is Vol_Cell in component model_parameters (litre).
C CONSTS(172) is Vol_SR in component model_parameters (litre).
C CONSTS(173) is Vol_SL in component model_parameters (litre).
C CONSTS(174) is Vol_jct in component model_parameters (litre).
C CONSTS(175) is Vol_cytosol in component model_parameters (litre).
C CONSTS(176) is Vol_mito in component model_parameters (litre).
C ALGBRC(22) is E_Na_jct in component reversal_potentials (millivolt).
C ALGBRC(40) is E_Na_SL in component reversal_potentials (millivolt).
C ALGBRC(45) is E_Ca_jct in component reversal_potentials (millivolt).
C ALGBRC(47) is E_Ca_SL in component reversal_potentials (millivolt).
C CONSTS(151) is E_K in component reversal_potentials (millivolt).
C CONSTS(152) is E_Cl in component reversal_potentials (millivolt).
C STATES(2) is Nai in component Na_buffer (millimolar).
C STATES(3) is Na_jct in component Na_buffer (millimolar).
C STATES(4) is Na_SL in component Na_buffer (millimolar).
C STATES(5) is Cai in component Ca_buffer (millimolar).
C STATES(6) is Ca_jct in component Ca_buffer (millimolar).
C STATES(7) is Ca_SL in component Ca_buffer (millimolar).
C ALGBRC(49) is i_Na_jct in component INa (microA_per_microF).
C ALGBRC(50) is i_Na_SL in component INa (microA_per_microF).
C CONSTS(22) is G_INa in component INa (milliS_per_microF).
C CONSTS(23) is Fx_Na_jct in component INa (dimensionless).
C CONSTS(24) is Fx_Na_SL in component INa (dimensionless).
C STATES(8) is m in component INa_m_gate (dimensionless).
C STATES(9) is h in component INa_h_gate (dimensionless).
C STATES(10) is j in component INa_j_gate (dimensionless).
C ALGBRC(48) is openProb in component INa (dimensionless).
C ALGBRC(1) is alpha_h in component INa_h_gate (per_millisecond).
C ALGBRC(18) is beta_h in component INa_h_gate (per_millisecond).
C ALGBRC(37) is tau_h in component INa_h_gate (millisecond).
C ALGBRC(42) is h_infinity in component INa_h_gate (dimensionless).
C ALGBRC(2) is alpha_j in component INa_j_gate (per_millisecond).
C ALGBRC(19) is beta_j in component INa_j_gate (per_millisecond).
C ALGBRC(38) is tau_j in component INa_j_gate (millisecond).
C ALGBRC(43) is j_infinity in component INa_j_gate (dimensionless).
C ALGBRC(3) is alpha_m in component INa_m_gate (per_millisecond).
C ALGBRC(20) is beta_m in component INa_m_gate (per_millisecond).
C ALGBRC(39) is tau_m in component INa_m_gate (millisecond).
C ALGBRC(44) is m_infinity in component INa_m_gate (dimensionless).
C ALGBRC(52) is i_Nab_jct in component INab (microA_per_microF).
C ALGBRC(53) is i_Nab_SL in component INab (microA_per_microF).
C CONSTS(25) is G_NaBk in component INab (milliS_per_microF).
C CONSTS(26) is Fx_NaBk_jct in component INab (dimensionless).
C CONSTS(27) is Fx_NaBk_SL in component INab (dimensionless).
C ALGBRC(56) is i_NaK_jct in component INaK (microA_per_microF).
C ALGBRC(57) is i_NaK_SL in component INaK (microA_per_microF).
C ALGBRC(55) is f_NaK in component INaK (dimensionless).
C CONSTS(28) is Q10_NaK in component INaK (dimensionless).
C CONSTS(153) is Q_NaK in component INaK (dimensionless).
C CONSTS(29) is H_NaK in component INaK (dimensionless).
C CONSTS(30) is Q10_Km_Nai in component INaK (dimensionless).
C CONSTS(154) is Q_Km_Nai in component INaK (dimensionless).
C CONSTS(31) is Km_Nai in component INaK (millimolar).
C CONSTS(32) is Km_Ko in component INaK (millimolar).
C CONSTS(33) is I_NaK_max in component INaK (microA_per_microF).
C CONSTS(155) is sigma in component INaK (dimensionless).
C CONSTS(34) is Fx_NaK_jct in component INaK (dimensionless).
C CONSTS(35) is Fx_NaK_SL in component INaK (dimensionless).
C CONSTS(156) is G_IKr in component IKr (milliS_per_microF).
C STATES(11) is Xr in component IKr_Xr_gate (dimensionless).
C ALGBRC(59) is Rr in component IKr_Rr_gate (dimensionless).
C ALGBRC(4) is Xr_infinity in component IKr_Xr_gate (dimensionless).
C ALGBRC(21) is tau_Xr in component IKr_Xr_gate (millisecond).
C ALGBRC(66) is i_Ks_jct in component IKs (microA_per_microF).
C ALGBRC(67) is i_Ks_SL in component IKs (microA_per_microF).
C ALGBRC(64) is G_Ks_SL in component IKs (milliS_per_microF).
C ALGBRC(63) is G_Ks_jct in component IKs (milliS_per_microF).
C CONSTS(36) is Fx_Ks_jct in component IKs (dimensionless).
C CONSTS(37) is Fx_Ks_SL in component IKs (dimensionless).
C ALGBRC(5) is E_Ks_jct in component IKs (millivolt).
C ALGBRC(6) is E_Ks_SL in component IKs (millivolt).
C ALGBRC(65) is E_Ks in component IKs (millivolt).
C CONSTS(38) is pKNa in component IKs (dimensionless).
C ALGBRC(61) is pCa_jct in component IKs (dimensionless).
C ALGBRC(62) is pCa_SL in component IKs (dimensionless).
C STATES(12) is Xs in component IKs_Xs_gate (dimensionless).
C ALGBRC(8) is Xs_infinity in component IKs_Xs_gate (dimensionless).
C ALGBRC(23) is tau_Xs in component IKs_Xs_gate (millisecond).
C CONSTS(39) is G_tos in component Itos (milliS_per_microF).
C STATES(13) is Y_tos in component Itos_Y_gate (dimensionless).
C STATES(14) is X_tos in component Itos_X_gate (dimensionless).
C ALGBRC(71) is R_tos in component Itos_R_gate (dimensionless).
C ALGBRC(9) is X_tos_infinity in component Itos_X_gate (dimensionless).
C ALGBRC(24) is tau_X_tos in component Itos_X_gate (millisecond).
C ALGBRC(10) is Y_tos_infinity in component Itos_Y_gate (dimensionless).
C ALGBRC(25) is tau_Y_tos in component Itos_Y_gate (millisecond).
C STATES(15) is R_tos_other in component Itos_R_gate (dimensionless).
C ALGBRC(69) is R_tos_infinity in component Itos_R_gate (dimensionless).
C ALGBRC(70) is tau_R_tos in component Itos_R_gate (millisecond).
C CONSTS(40) is G_tof in component Itof (milliS_per_microF).
C STATES(16) is Y_tof in component Itof_Y_gate (dimensionless).
C STATES(17) is X_tof in component Itof_X_gate (dimensionless).
C ALGBRC(11) is X_tof_infinity in component Itof_X_gate (dimensionless).
C ALGBRC(26) is tau_X_tof in component Itof_X_gate (millisecond).
C ALGBRC(12) is Y_tof_infinity in component Itof_Y_gate (dimensionless).
C ALGBRC(27) is tau_Y_tof in component Itof_Y_gate (millisecond).
C CONSTS(157) is G_K1 in component IK1 (milliS_per_microF).
C ALGBRC(76) is K1_infinity in component IK1_K1_gate (dimensionless).
C ALGBRC(74) is alpha_K1 in component IK1_K1_gate (per_millisecond).
C ALGBRC(75) is beta_K1 in component IK1_K1_gate (per_millisecond).
C CONSTS(41) is g_Kp in component IKp (milliS_per_microF).
C CONSTS(42) is G_Cl in component ICl_Ca (milliS_per_microF).
C CONSTS(43) is Kd_ClCa in component ICl_Ca (millimolar).
C CONSTS(44) is Fx_Cl_jct in component ICl_Ca (dimensionless).
C CONSTS(45) is Fx_Cl_SL in component ICl_Ca (dimensionless).
C CONSTS(46) is G_ClBk in component IClb (milliS_per_microF).
C ALGBRC(83) is i_CaL_Ca_SL in component ICaL (microA_per_microF).
C ALGBRC(86) is i_CaL_Ca_jct in component ICaL (microA_per_microF).
C ALGBRC(84) is i_CaL_Na_SL in component ICaL (microA_per_microF).
C ALGBRC(87) is i_CaL_Na_jct in component ICaL (microA_per_microF).
C ALGBRC(88) is i_CaL_K in component ICaL (microA_per_microF).
C CONSTS(47) is PCa in component ICaL (litre_per_farad_millisecond).
C CONSTS(48) is PNa in component ICaL (litre_per_farad_millisecond).
C CONSTS(49) is PK in component ICaL (litre_per_farad_millisecond).
C CONSTS(50) is Fx_ICaL_jct in component ICaL (dimensionless).
C CONSTS(51) is Fx_ICaL_SL in component ICaL (dimensionless).
C CONSTS(52) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(53) is gamma_Cao in component ICaL (dimensionless).
C CONSTS(54) is gamma_Nai in component ICaL (dimensionless).
C CONSTS(55) is gamma_Nao in component ICaL (dimensionless).
C CONSTS(56) is gamma_Ki in component ICaL (dimensionless).
C CONSTS(57) is gamma_Ko in component ICaL (dimensionless).
C CONSTS(58) is Q10_CaL in component ICaL (dimensionless).
C CONSTS(158) is Q_CaL in component ICaL (dimensionless).
C STATES(18) is d in component ICaL_d_gate (dimensionless).
C STATES(19) is f in component ICaL_f_gate (dimensionless).
C ALGBRC(82) is fCa_SL in component ICaL_fCa_gate (dimensionless).
C ALGBRC(85) is fCa_jct in component ICaL_fCa_gate (dimensionless).
C ALGBRC(81) is temp in component ICaL (coulomb_per_mole).
C ALGBRC(13) is d_infinity in component ICaL_d_gate (dimensionless).
C ALGBRC(28) is tau_d in component ICaL_d_gate (millisecond).
C ALGBRC(14) is f_infinity in component ICaL_f_gate (dimensionless).
C ALGBRC(29) is tau_f in component ICaL_f_gate (millisecond).
C STATES(20) is fCaB_SL in component ICaL_fCa_gate (dimensionless).
C STATES(21) is fCaB_jct in component ICaL_fCa_gate (dimensionless).
C ALGBRC(95) is i_NaCa_SL in component INaCa (microA_per_microF).
C ALGBRC(94) is i_NaCa_jct in component INaCa (microA_per_microF).
C CONSTS(59) is V_max in component INaCa (microA_per_microF).
C CONSTS(60) is Fx_NCX_jct in component INaCa (dimensionless).
C CONSTS(61) is Fx_NCX_SL in component INaCa (dimensionless).
C CONSTS(62) is Q10_NCX in component INaCa (dimensionless).
C CONSTS(159) is Q_NCX in component INaCa (dimensionless).
C CONSTS(63) is K_mNai in component INaCa (millimolar).
C CONSTS(64) is K_mCao in component INaCa (millimolar).
C CONSTS(65) is K_mNao in component INaCa (millimolar).
C CONSTS(66) is K_mCai in component INaCa (millimolar).
C ALGBRC(92) is Ka_SL in component INaCa (dimensionless).
C ALGBRC(93) is Ka_jct in component INaCa (dimensionless).
C CONSTS(67) is Kd_act in component INaCa (millimolar).
C CONSTS(68) is ksat in component INaCa (dimensionless).
C CONSTS(69) is eta in component INaCa (dimensionless).
C CONSTS(70) is HNa in component INaCa (dimensionless).
C ALGBRC(90) is temp_jct in component INaCa (millimolar4).
C ALGBRC(91) is temp_SL in component INaCa (millimolar4).
C ALGBRC(102) is i_Cap_SL in component ICap (microA_per_microF).
C ALGBRC(100) is i_Cap_jct in component ICap (microA_per_microF).
C CONSTS(71) is V_max in component ICap (millimolar_per_millisecond).
C CONSTS(72) is Fx_SLCaP_jct in component ICap (dimensionless).
C CONSTS(73) is Fx_SLCaP_SL in component ICap (dimensionless).
C CONSTS(74) is Q10_SLCaP in component ICap (dimensionless).
C CONSTS(160) is Q_SLCaP in component ICap (dimensionless).
C CONSTS(75) is Km in component ICap (millimolar).
C CONSTS(76) is H in component ICap (dimensionless).
C CONSTS(177) is V_maxAF_old in component ICap (microA_per_microF).
C CONSTS(77) is V_maxAF in component ICap (microA_per_microF).
C ALGBRC(104) is i_Cab_jct in component ICab (microA_per_microF).
C ALGBRC(105) is i_Cab_SL in component ICab (microA_per_microF).
C CONSTS(78) is G_CaBk in component ICab (milliS_per_microF).
C CONSTS(79) is Fx_CaBk_jct in component ICab (dimensionless).
C CONSTS(80) is Fx_CaBk_SL in component ICab (dimensionless).
C ALGBRC(107) is j_rel_SR in component Jrel_SR (millimolar_per_millisecond).
C CONSTS(81) is Max_SR in component Jrel_SR (dimensionless).
C CONSTS(82) is Min_SR in component Jrel_SR (dimensionless).
C CONSTS(83) is EC50_SR in component Jrel_SR (millimolar).
C STATES(22) is R in component Jrel_SR (dimensionless).
C STATES(23) is I in component Jrel_SR (dimensionless).
C ALGBRC(46) is RI in component Jrel_SR (dimensionless).
C STATES(24) is O in component Jrel_SR (dimensionless).
C CONSTS(84) is ks in component Jrel_SR (per_millisecond).
C CONSTS(85) is koCa in component Jrel_SR (per_millimolar2_per_millisecond).
C CONSTS(86) is kom in component Jrel_SR (per_millisecond).
C CONSTS(87) is kiCa in component Jrel_SR (per_millimolar_per_millisecond).
C CONSTS(88) is kim in component Jrel_SR (per_millisecond).
C CONSTS(89) is HSR in component Jrel_SR (dimensionless).
C STATES(25) is Ca_SR in component Ca_buffer (millimolar).
C ALGBRC(15) is kCaSR in component Jrel_SR (dimensionless).
C ALGBRC(30) is koSRCa in component Jrel_SR (per_millimolar2_per_millisecond).
C ALGBRC(41) is kiSRCa in component Jrel_SR (per_millimolar_per_millisecond).
C ALGBRC(108) is j_leak_SR in component Jleak_SR (millimolar_per_millisecond).
C CONSTS(90) is KSRleak in component Jleak_SR (per_millisecond).
C ALGBRC(109) is j_pump_SR in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(91) is V_max in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(92) is Q10_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(161) is Q_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(93) is Kmf in component Jpump_SR (millimolar).
C CONSTS(94) is Kmr in component Jpump_SR (millimolar).
C CONSTS(95) is H in component Jpump_SR (dimensionless).
C ALGBRC(96) is J_Na_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(99) is J_Na_SL_cytosol in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(111) is J_Ca_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(112) is J_Ca_SL_cytosol in component ion_diffusion (millimole_per_millisecond).
C CONSTS(96) is D_Ca_jct_SL in component ion_diffusion (dm2_per_second).
C CONSTS(97) is D_Ca_SL_cytosol in component ion_diffusion (dm2_per_second).
C CONSTS(98) is D_Na_jct_SL in component ion_diffusion (dm2_per_second).
C CONSTS(99) is D_Na_SL_cytosol in component ion_diffusion (dm2_per_second).
C CONSTS(100) is x_jct_SL in component ion_diffusion (micrometre).
C CONSTS(101) is x_SL_cytosol in component ion_diffusion (micrometre).
C CONSTS(102) is A_jct_SL in component ion_diffusion (cm2).
C CONSTS(103) is A_SL_cytosol in component ion_diffusion (cm2).
C CONSTS(162) is j_Na_jct_SL in component ion_diffusion (litre_per_millisecond).
C CONSTS(163) is j_Na_SL_cytosol in component ion_diffusion (litre_per_millisecond).
C CONSTS(164) is j_Ca_jct_SL in component ion_diffusion (litre_per_millisecond).
C CONSTS(165) is j_Ca_SL_cytosol in component ion_diffusion (litre_per_millisecond).
C CONSTS(166) is j_Na_jct_SL_new in component ion_diffusion (litre_per_millisecond).
C CONSTS(167) is j_Na_SL_cytosol_new in component ion_diffusion (litre_per_millisecond).
C CONSTS(168) is j_Ca_jct_SL_new in component ion_diffusion (litre_per_millisecond).
C CONSTS(169) is j_Ca_SL_cytosol_new in component ion_diffusion (litre_per_millisecond).
C STATES(26) is Na_SL_buf in component Na_buffer (millimolar).
C STATES(27) is Na_jct_buf in component Na_buffer (millimolar).
C CONSTS(104) is Bmax_SL in component Na_buffer (millimolar).
C CONSTS(105) is Bmax_jct in component Na_buffer (millimolar).
C CONSTS(106) is kon in component Na_buffer (per_millimolar_per_millisecond).
C CONSTS(107) is koff in component Na_buffer (per_millisecond).
C ALGBRC(98) is dNa_jct_buf in component Na_buffer (millimolar_per_millisecond).
C ALGBRC(101) is dNa_SL_buf in component Na_buffer (millimolar_per_millisecond).
C STATES(28) is Ca_SLB_SL in component Ca_buffer (millimolar).
C STATES(29) is Ca_SLB_jct in component Ca_buffer (millimolar).
C STATES(30) is Ca_SLHigh_SL in component Ca_buffer (millimolar).
C STATES(31) is Ca_SLHigh_jct in component Ca_buffer (millimolar).
C STATES(32) is Ca_Calsequestrin in component Ca_buffer (millimolar).
C CONSTS(108) is Bmax_SLB_SL in component Ca_buffer (millimolar).
C CONSTS(109) is Bmax_SLB_jct in component Ca_buffer (millimolar).
C CONSTS(110) is Bmax_SLHigh_SL in component Ca_buffer (millimolar).
C CONSTS(111) is Bmax_SLHigh_jct in component Ca_buffer (millimolar).
C CONSTS(112) is Bmax_Calsequestrin in component Ca_buffer (millimolar).
C CONSTS(113) is kon_SL in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(114) is kon_Calsequestrin in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(115) is koff_SLB in component Ca_buffer (per_millisecond).
C CONSTS(116) is koff_SLHigh in component Ca_buffer (per_millisecond).
C CONSTS(117) is koff_Calsequestrin in component Ca_buffer (per_millisecond).
C ALGBRC(110) is dCalsequestrin in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(128) is dCa_cytosol_tot_bound in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(115) is dCa_SLB_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(113) is dCa_SLB_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(118) is dCa_SLHigh_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(116) is dCa_SLHigh_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(119) is dCa_jct_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(122) is dCa_SL_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(121) is i_Ca_jct_tot in component Ca_buffer (microA_per_microF).
C ALGBRC(124) is i_Ca_SL_tot in component Ca_buffer (microA_per_microF).
C STATES(33) is Ca_TroponinC in component cytosolic_Ca_buffer (millimolar).
C STATES(34) is Ca_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar).
C STATES(35) is Mg_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar).
C STATES(36) is Ca_Calmodulin in component cytosolic_Ca_buffer (millimolar).
C STATES(37) is Ca_Myosin in component cytosolic_Ca_buffer (millimolar).
C STATES(38) is Mg_Myosin in component cytosolic_Ca_buffer (millimolar).
C STATES(39) is Ca_SRB in component cytosolic_Ca_buffer (millimolar).
C CONSTS(118) is Bmax_TroponinC in component cytosolic_Ca_buffer (millimolar).
C CONSTS(119) is Bmax_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(120) is Bmax_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(121) is Bmax_Calmodulin in component cytosolic_Ca_buffer (millimolar).
C CONSTS(122) is Bmax_Myosin_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(123) is Bmax_Myosin_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(124) is Bmax_SRB in component cytosolic_Ca_buffer (millimolar).
C CONSTS(125) is kon_TroponinC in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(126) is kon_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(127) is kon_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(128) is kon_Calmodulin in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(129) is kon_Myosin_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(130) is kon_Myosin_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(131) is kon_SRB in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(132) is koff_TroponinC in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(133) is koff_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(134) is koff_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(135) is koff_Calmodulin in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(136) is koff_Myosin_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(137) is koff_Myosin_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(138) is koff_SRB in component cytosolic_Ca_buffer (per_millisecond).
C ALGBRC(114) is dCa_TroponinC in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(117) is dCa_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(120) is dMg_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(123) is dCa_Calmodulin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(125) is dCa_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(126) is dMg_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(127) is dCa_SRB in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C CONSTS(139) is Indo1 in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(140) is Fluo3 in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(40) is Ca_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(41) is Ca_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(42) is Ca_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(43) is Ca_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(44) is Ca_Fluo3_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C STATES(45) is Ca_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(141) is Bmax_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(142) is Bmax_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(143) is Bmax_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(144) is Bmax_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(145) is Bmax_Fluo3_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(146) is Bmax_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(147) is kon_Indo1 in component indo_fluo_Ca_buffer_not_connected (per_millimolar2_per_millisecond).
C CONSTS(148) is kon_Fluo3 in component indo_fluo_Ca_buffer_not_connected (per_millimolar2_per_millisecond).
C CONSTS(149) is koff_Indo1 in component indo_fluo_Ca_buffer_not_connected (per_millisecond).
C CONSTS(150) is koff_Fluo3 in component indo_fluo_Ca_buffer_not_connected (per_millisecond).
C ALGBRC(16) is Indo1Bound in component indo_fluo_Ca_buffer_not_connected (millimolar).
C ALGBRC(17) is Fluo3Bound in component indo_fluo_Ca_buffer_not_connected (millimolar).
C ALGBRC(31) is dCa_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(32) is dCa_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(33) is dCa_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(34) is dCa_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(35) is dCa_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(36) is dCa_Fluo3_SL in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C RATES(1) is d/dt V in component cell (millivolt).
C RATES(9) is d/dt h in component INa_h_gate (dimensionless).
C RATES(10) is d/dt j in component INa_j_gate (dimensionless).
C RATES(8) is d/dt m in component INa_m_gate (dimensionless).
C RATES(11) is d/dt Xr in component IKr_Xr_gate (dimensionless).
C RATES(12) is d/dt Xs in component IKs_Xs_gate (dimensionless).
C RATES(14) is d/dt X_tos in component Itos_X_gate (dimensionless).
C RATES(13) is d/dt Y_tos in component Itos_Y_gate (dimensionless).
C RATES(15) is d/dt R_tos_other in component Itos_R_gate (dimensionless).
C RATES(17) is d/dt X_tof in component Itof_X_gate (dimensionless).
C RATES(16) is d/dt Y_tof in component Itof_Y_gate (dimensionless).
C RATES(18) is d/dt d in component ICaL_d_gate (dimensionless).
C RATES(19) is d/dt f in component ICaL_f_gate (dimensionless).
C RATES(20) is d/dt fCaB_SL in component ICaL_fCa_gate (dimensionless).
C RATES(21) is d/dt fCaB_jct in component ICaL_fCa_gate (dimensionless).
C RATES(22) is d/dt R in component Jrel_SR (dimensionless).
C RATES(24) is d/dt O in component Jrel_SR (dimensionless).
C RATES(23) is d/dt I in component Jrel_SR (dimensionless).
C RATES(27) is d/dt Na_jct_buf in component Na_buffer (millimolar).
C RATES(26) is d/dt Na_SL_buf in component Na_buffer (millimolar).
C RATES(3) is d/dt Na_jct in component Na_buffer (millimolar).
C RATES(4) is d/dt Na_SL in component Na_buffer (millimolar).
C RATES(2) is d/dt Nai in component Na_buffer (millimolar).
C RATES(32) is d/dt Ca_Calsequestrin in component Ca_buffer (millimolar).
C RATES(28) is d/dt Ca_SLB_SL in component Ca_buffer (millimolar).
C RATES(29) is d/dt Ca_SLB_jct in component Ca_buffer (millimolar).
C RATES(30) is d/dt Ca_SLHigh_SL in component Ca_buffer (millimolar).
C RATES(31) is d/dt Ca_SLHigh_jct in component Ca_buffer (millimolar).
C RATES(25) is d/dt Ca_SR in component Ca_buffer (millimolar).
C RATES(6) is d/dt Ca_jct in component Ca_buffer (millimolar).
C RATES(7) is d/dt Ca_SL in component Ca_buffer (millimolar).
C RATES(5) is d/dt Cai in component Ca_buffer (millimolar).
C RATES(33) is d/dt Ca_TroponinC in component cytosolic_Ca_buffer (millimolar).
C RATES(34) is d/dt Ca_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar).
C RATES(35) is d/dt Mg_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar).
C RATES(36) is d/dt Ca_Calmodulin in component cytosolic_Ca_buffer (millimolar).
C RATES(37) is d/dt Ca_Myosin in component cytosolic_Ca_buffer (millimolar).
C RATES(38) is d/dt Mg_Myosin in component cytosolic_Ca_buffer (millimolar).
C RATES(39) is d/dt Ca_SRB in component cytosolic_Ca_buffer (millimolar).
C RATES(40) is d/dt Ca_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C RATES(41) is d/dt Ca_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C RATES(42) is d/dt Ca_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C RATES(43) is d/dt Ca_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C RATES(44) is d/dt Ca_Fluo3_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C RATES(45) is d/dt Ca_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -85.719687955637
CONSTS(1) = 100
CONSTS(2) = 1000
CONSTS(3) = 5
CONSTS(4) = 9.5
CONSTS(5) = 5.4
CONSTS(6) = 140
CONSTS(7) = 1.8
CONSTS(8) = 150
CONSTS(9) = 135
CONSTS(10) = 1
CONSTS(11) = 1
CONSTS(12) = 15
CONSTS(13) = 8314.3
CONSTS(14) = 310
CONSTS(15) = 96486.7
CONSTS(16) = 2e-6
CONSTS(17) = 100
CONSTS(18) = 10.25
CONSTS(19) = 0.16
CONSTS(20) = 0.015
CONSTS(21) = 0.045
STATES(2) = 8.874461106492
STATES(3) = 8.872823559072
STATES(4) = 8.874077316753
STATES(5) = 0.000087350002
STATES(6) = 0.000174843061
STATES(7) = 0.000106395937
CONSTS(22) = 16
CONSTS(23) = 0.11
CONSTS(24) = 0.89
STATES(8) = 0.001370685156
STATES(9) = 0.987140350343
STATES(10) = 0.991822731369
CONSTS(25) = 0.297e-3
CONSTS(26) = 0.11
CONSTS(27) = 0.89
CONSTS(28) = 1.63
CONSTS(29) = 4
CONSTS(30) = 1.49
CONSTS(31) = 11
CONSTS(32) = 1.5
CONSTS(33) = 1.91
CONSTS(34) = 0.11
CONSTS(35) = 0.89
STATES(11) = 0.008471550841
CONSTS(36) = 0.11
CONSTS(37) = 0.89
CONSTS(38) = 0.01833
STATES(12) = 0.00687399199
CONSTS(39) = 0.02
STATES(13) = 0.293519921626
STATES(14) = 0.004011272375
STATES(15) = 0.383430556383
CONSTS(40) = 0.06
STATES(16) = 0.9946314893
STATES(17) = 0.00401120993
CONSTS(41) = 0.001
CONSTS(42) = 0.109625
CONSTS(43) = 0.1
CONSTS(44) = 0.11
CONSTS(45) = 0.89
CONSTS(46) = 0.009
CONSTS(47) = 5.4e-4
CONSTS(48) = 1.5e-8
CONSTS(49) = 2.7e-7
CONSTS(50) = 0.9
CONSTS(51) = 0.1
CONSTS(52) = 0.341
CONSTS(53) = 0.341
CONSTS(54) = 0.75
CONSTS(55) = 0.75
CONSTS(56) = 0.75
CONSTS(57) = 0.75
CONSTS(58) = 1.8
STATES(18) = 0.000006997531
STATES(19) = 1.000675515962
STATES(20) = 0.015352888928
STATES(21) = 0.024609183734
CONSTS(59) = 9
CONSTS(60) = 0.11
CONSTS(61) = 0.89
CONSTS(62) = 1.57
CONSTS(63) = 12.29
CONSTS(64) = 1.3
CONSTS(65) = 87.5
CONSTS(66) = 0.00359
CONSTS(67) = 0.000256
CONSTS(68) = 0.27
CONSTS(69) = 0.35
CONSTS(70) = 3
CONSTS(71) = 2.2e-6
CONSTS(72) = 0.11
CONSTS(73) = 0.89
CONSTS(74) = 2.35
CONSTS(75) = 0.0005
CONSTS(76) = 1.6
CONSTS(77) = 0.0673
CONSTS(78) = 0.0002513
CONSTS(79) = 0.11
CONSTS(80) = 0.89
CONSTS(81) = 15
CONSTS(82) = 1
CONSTS(83) = 0.45
STATES(22) = 0.884673513138
STATES(23) = 0.00000009272
STATES(24) = 0.000000711264
CONSTS(84) = 25
CONSTS(85) = 10
CONSTS(86) = 0.06
CONSTS(87) = 0.5
CONSTS(88) = 0.005
CONSTS(89) = 2.5
STATES(25) = 0.545611267699
CONSTS(90) = 5.348e-6
CONSTS(91) = 286e-6
CONSTS(92) = 2.6
CONSTS(93) = 0.000246
CONSTS(94) = 1.7
CONSTS(95) = 1.787
CONSTS(96) = 1.64e-6
CONSTS(97) = 1.22e-6
CONSTS(98) = 1.09e-5
CONSTS(99) = 1.79e-5
CONSTS(100) = 0.5
CONSTS(101) = 0.45
CONSTS(102) = 3.01e-6
CONSTS(103) = 1.3e-4
STATES(26) = 0.776121392467
STATES(27) = 3.557055389701
CONSTS(104) = 1.65
CONSTS(105) = 3.7
CONSTS(106) = 0.0001
CONSTS(107) = 1e-3
STATES(28) = 0.009868629147
STATES(29) = 0.007780801995
STATES(30) = 0.114438990328
STATES(31) = 0.077503874257
STATES(32) = 1.186496899338
CONSTS(108) = 0.0374
CONSTS(109) = 0.0046
CONSTS(110) = 0.0134
CONSTS(111) = 0.00165
CONSTS(112) = 0.14
CONSTS(113) = 100
CONSTS(114) = 100
CONSTS(115) = 1.3
CONSTS(116) = 30e-3
CONSTS(117) = 65
STATES(33) = 0.008963736337
STATES(34) = 0.117995194438
STATES(35) = 0.010337654274
STATES(36) = 0.000295961245
STATES(37) = 0.001984672275
STATES(38) = 0.137497736234
STATES(39) = 0.002177112381
CONSTS(118) = 0.07
CONSTS(119) = 0.14
CONSTS(120) = 0.14
CONSTS(121) = 0.024
CONSTS(122) = 0.14
CONSTS(123) = 0.14
CONSTS(124) = 0.019
CONSTS(125) = 32.7
CONSTS(126) = 2.37
CONSTS(127) = 3e-3
CONSTS(128) = 34
CONSTS(129) = 13.8
CONSTS(130) = 15.7e-3
CONSTS(131) = 100
CONSTS(132) = 19.6e-3
CONSTS(133) = 0.032e-3
CONSTS(134) = 3.33e-3
CONSTS(135) = 238e-3
CONSTS(136) = 0.46e-3
CONSTS(137) = 0.057e-3
CONSTS(138) = 60e-3
CONSTS(139) = 0
CONSTS(140) = 0
STATES(40) = 0
STATES(41) = 0
STATES(42) = 0
STATES(43) = 0
STATES(44) = 0
STATES(45) = 0
CONSTS(141) = 0.025
CONSTS(142) = 0.00077
CONSTS(143) = 0.00002
CONSTS(144) = 0.025
CONSTS(145) = 0.00077
CONSTS(146) = 0.00002
CONSTS(147) = 100
CONSTS(148) = 100
CONSTS(149) = 60e-3
CONSTS(150) = 110e-3
CONSTS(151) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(5)/CONSTS(9))
CONSTS(152) = (( - CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(8)/CONSTS(12))
CONSTS(153) = CONSTS(28) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(154) = CONSTS(30) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(155) = (EXP(CONSTS(6)/67.3000) - 1.00000)/7.00000
CONSTS(156) = 0.0300000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(157) = 0.900000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(158) = CONSTS(58) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(159) = CONSTS(62) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(160) = CONSTS(74) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(161) = CONSTS(92) ** (CONSTS(14) - 310.000)/10.0000
CONSTS(162) = ( CONSTS(98)*CONSTS(102))/CONSTS(100)
CONSTS(163) = ( CONSTS(99)*CONSTS(103))/CONSTS(101)
CONSTS(164) = ( CONSTS(96)*CONSTS(102))/CONSTS(100)
CONSTS(165) = ( CONSTS(97)*CONSTS(103))/CONSTS(101)
CONSTS(166) = 1.83130e-14
CONSTS(167) = 1.63860e-12
CONSTS(168) = 8.24130e-13
CONSTS(169) = 3.72430e-12
CONSTS(170) = ( (( CONSTS(16)*2.00000*CONSTS(18))/10000.0)* 3.14159265358979*CONSTS(17))/10000.0
CONSTS(171) = ( 3.14159*CONSTS(18)/1000.00 ** 2.00000*CONSTS(17))/1000.00 ** 3.00000
CONSTS(172) = 0.0350000*CONSTS(171)
CONSTS(173) = 0.0200000*CONSTS(171)
CONSTS(174) = 0.000510000*CONSTS(171)
CONSTS(175) = 0.650000*CONSTS(171)
CONSTS(176) = CONSTS(171) - (CONSTS(172)+CONSTS(173)+CONSTS(174)+CONSTS(175))
CONSTS(177) = ( CONSTS(71)*CONSTS(175)*2.00000*CONSTS(15))/CONSTS(170)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(20) = 1.70000*STATES(7)*(1.00000 - STATES(20)) - 11.9000*STATES(20)
RATES(21) = 1.70000*STATES(6)*(1.00000 - STATES(21)) - 11.9000*STATES(21)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (50.0000+STATES(1))/7.50000))
ALGBRC(21) = 1.00000/(( 0.00138000*(STATES(1)+7.00000))/(1.00000 - EXP( - 0.123000*(STATES(1)+7.00000)))+( 0.000610000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000))
RATES(11) = (ALGBRC(4) - STATES(11))/ALGBRC(21)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(23) = 1.00000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(- 1.00000+EXP( 0.0687000*(STATES(1)+30.0000))))
RATES(12) = (ALGBRC(8) - STATES(12))/ALGBRC(23)
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(24) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
RATES(14) = (ALGBRC(9) - STATES(14))/ALGBRC(24)
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(25) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
RATES(13) = (ALGBRC(10) - STATES(13))/ALGBRC(25)
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(26) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
RATES(17) = (ALGBRC(11) - STATES(17))/ALGBRC(26)
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(27) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
RATES(16) = (ALGBRC(12) - STATES(16))/ALGBRC(27)
ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(28) = ( 1.00000*ALGBRC(13)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
RATES(18) = (ALGBRC(13) - STATES(18))/ALGBRC(28)
ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(29) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
RATES(19) = (ALGBRC(14) - STATES(19))/ALGBRC(29)
ALGBRC(16) = STATES(40)+( STATES(42)*CONSTS(174))/CONSTS(175)+( STATES(41)*CONSTS(173))/CONSTS(175)
ALGBRC(31) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(5)*(CONSTS(141) - STATES(40)) - CONSTS(149)*STATES(40)
RATES(40) = ALGBRC(31)
ALGBRC(33) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(7)*(( CONSTS(142)*CONSTS(175))/CONSTS(173) - STATES(41)) - CONSTS(149)*STATES(41)
RATES(41) = ALGBRC(33)
ALGBRC(32) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(6)*(( CONSTS(143)*CONSTS(175))/CONSTS(174) - STATES(42)) - CONSTS(149)*STATES(42)
RATES(42) = ALGBRC(32)
ALGBRC(17) = STATES(43)+( STATES(45)*CONSTS(174))/CONSTS(175)+( STATES(44)*CONSTS(173))/CONSTS(175)
ALGBRC(34) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(5)*(CONSTS(144) - STATES(43)) - CONSTS(150)*STATES(43)
RATES(43) = ALGBRC(34)
ALGBRC(36) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(7)*(( CONSTS(145)*CONSTS(175))/CONSTS(173) - STATES(44)) - CONSTS(150)*STATES(44)
RATES(44) = ALGBRC(36)
ALGBRC(35) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(6)*(( CONSTS(146)*CONSTS(175))/CONSTS(174) - STATES(45)) - CONSTS(150)*STATES(45)
RATES(45) = ALGBRC(35)
ALGBRC(15) = CONSTS(81) - (CONSTS(81) - CONSTS(82))/(1.00000+CONSTS(83)/STATES(25) ** CONSTS(89))
ALGBRC(30) = CONSTS(85)/ALGBRC(15)
ALGBRC(41) = CONSTS(87)*ALGBRC(15)
RATES(24) = ( ALGBRC(30)*STATES(6) ** 2.00000*STATES(22) - CONSTS(86)*STATES(24)) - ( ALGBRC(41)*STATES(6)*STATES(24) - CONSTS(88)*STATES(23))
ALGBRC(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(18) = TERNRY(STATES(1).LT.- 40.0000, 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)), 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))))
ALGBRC(37) = 1.00000/(ALGBRC(1)+ALGBRC(18))
ALGBRC(42) = ALGBRC(1)/(ALGBRC(1)+ALGBRC(18))
RATES(9) = (ALGBRC(42) - STATES(9))/ALGBRC(37)
ALGBRC(2) = TERNRY(STATES(1).LT.- 40.0000, (( ( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/1.00000)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))), 0.00000)
ALGBRC(19) = TERNRY(STATES(1).LT.- 40.0000, ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))), ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))
ALGBRC(38) = 1.00000/(ALGBRC(2)+ALGBRC(19))
ALGBRC(43) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(19))
RATES(10) = (ALGBRC(43) - STATES(10))/ALGBRC(38)
ALGBRC(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(20) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(39) = 1.00000/(ALGBRC(3)+ALGBRC(20))
ALGBRC(44) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(20))
RATES(8) = (ALGBRC(44) - STATES(8))/ALGBRC(39)
ALGBRC(46) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
RATES(22) = ( CONSTS(88)*ALGBRC(46) - ALGBRC(41)*STATES(6)*STATES(22)) - ( ALGBRC(30)*STATES(6) ** 2.00000*STATES(22) - CONSTS(86)*STATES(24))
RATES(23) = ( ALGBRC(41)*STATES(6)*STATES(24) - CONSTS(88)*STATES(23)) - ( CONSTS(86)*STATES(23) - ALGBRC(30)*STATES(6) ** 2.00000*ALGBRC(46))
ALGBRC(69) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(71) = ALGBRC(69)
ALGBRC(70) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
RATES(15) = (ALGBRC(69) - ALGBRC(71))/ALGBRC(70)
ALGBRC(98) = CONSTS(106)*STATES(3)*(CONSTS(105) - STATES(27)) - CONSTS(107)*STATES(27)
RATES(27) = ALGBRC(98)
ALGBRC(22) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(6)/STATES(3))
ALGBRC(48) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(49) = CONSTS(23)*CONSTS(22)*ALGBRC(48)*(STATES(1) - ALGBRC(22))
ALGBRC(52) = CONSTS(26)*CONSTS(25)*(STATES(1) - ALGBRC(22))
ALGBRC(55) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))+ 0.0365000*CONSTS(155)*EXP(( - STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(56) = ( (( CONSTS(34)*CONSTS(153)*CONSTS(33)*ALGBRC(55))/(1.00000+( CONSTS(154)*CONSTS(31))/STATES(3) ** CONSTS(29)))*CONSTS(5))/(CONSTS(5)+CONSTS(32))
ALGBRC(85) = 1.00000 - STATES(21)
ALGBRC(81) = ( STATES(18)*STATES(19)*CONSTS(158)*STATES(1)*CONSTS(15) ** 2.00000)/( CONSTS(13)*CONSTS(14))
ALGBRC(87) = ( ALGBRC(81)*ALGBRC(85)*CONSTS(50)*CONSTS(48)*( CONSTS(54)*STATES(3)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(55)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(93) = 1.00000/(1.00000+CONSTS(67)/STATES(6) ** 3.00000)
ALGBRC(90) = ( EXP(( CONSTS(69)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*STATES(3) ** CONSTS(70)*CONSTS(7) - EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*CONSTS(6) ** CONSTS(70)*STATES(6))/(1.00000+ CONSTS(68)*EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(94) = ( CONSTS(60)*CONSTS(59)*ALGBRC(93)*CONSTS(159)*ALGBRC(90))/( CONSTS(66)*CONSTS(6) ** CONSTS(70)*(1.00000+STATES(3)/CONSTS(63) ** CONSTS(70))+ CONSTS(65) ** CONSTS(70)*STATES(6)*(1.00000+STATES(6)/CONSTS(66))+ CONSTS(64)*STATES(3) ** CONSTS(70)+ STATES(3) ** CONSTS(70)*CONSTS(7)+ CONSTS(6) ** CONSTS(70)*STATES(6))
ALGBRC(96) = (STATES(3) - STATES(4))*CONSTS(162)
RATES(3) = (( - CONSTS(170)*(ALGBRC(49)+ 3.00000*ALGBRC(94)+ALGBRC(52)+ 3.00000*ALGBRC(56)+ALGBRC(87)))/( CONSTS(174)*CONSTS(15)) - ALGBRC(96)/CONSTS(174)) - ALGBRC(98)
ALGBRC(99) = (STATES(4) - STATES(2))*CONSTS(163)
RATES(2) = ALGBRC(99)/CONSTS(175)
ALGBRC(101) = CONSTS(106)*STATES(4)*(CONSTS(104) - STATES(26)) - CONSTS(107)*STATES(26)
RATES(26) = ALGBRC(101)
ALGBRC(40) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(6)/STATES(4))
ALGBRC(50) = CONSTS(24)*CONSTS(22)*ALGBRC(48)*(STATES(1) - ALGBRC(40))
ALGBRC(53) = CONSTS(27)*CONSTS(25)*(STATES(1) - ALGBRC(40))
ALGBRC(57) = ( (( CONSTS(35)*CONSTS(153)*CONSTS(33)*ALGBRC(55))/(1.00000+( CONSTS(154)*CONSTS(31))/STATES(4) ** CONSTS(29)))*CONSTS(5))/(CONSTS(5)+CONSTS(32))
ALGBRC(82) = 1.00000 - STATES(20)
ALGBRC(84) = ( ALGBRC(81)*ALGBRC(82)*CONSTS(51)*CONSTS(48)*( CONSTS(54)*STATES(4)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(55)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(92) = 1.00000/(1.00000+CONSTS(67)/STATES(7) ** 3.00000)
ALGBRC(91) = ( EXP(( CONSTS(69)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*STATES(4) ** CONSTS(70)*CONSTS(7) - EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*CONSTS(6) ** CONSTS(70)*STATES(7))/(1.00000+ CONSTS(68)*EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(95) = ( CONSTS(61)*CONSTS(59)*ALGBRC(92)*CONSTS(159)*ALGBRC(91))/( CONSTS(66)*CONSTS(6) ** CONSTS(70)*(1.00000+STATES(4)/CONSTS(63) ** CONSTS(70))+ CONSTS(65) ** CONSTS(70)*STATES(7)*(1.00000+STATES(7)/CONSTS(66))+ CONSTS(64)*STATES(4) ** CONSTS(70)+ STATES(4) ** CONSTS(70)*CONSTS(7)+ CONSTS(6) ** CONSTS(70)*STATES(7))
RATES(4) = (( - CONSTS(170)*(ALGBRC(50)+ 3.00000*ALGBRC(95)+ALGBRC(53)+ 3.00000*ALGBRC(57)+ALGBRC(84)))/( CONSTS(173)*CONSTS(15))+(ALGBRC(96) - ALGBRC(99))/CONSTS(173)) - ALGBRC(101)
ALGBRC(51) = ALGBRC(49)+ALGBRC(50)
ALGBRC(54) = ALGBRC(52)+ALGBRC(53)
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(59) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(60) = CONSTS(156)*STATES(11)*ALGBRC(59)*(STATES(1) - CONSTS(151))
ALGBRC(61) = - arbitrary_log(STATES(6)/1.00000+3.00000, 10)
ALGBRC(63) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP(- 7.20000+ALGBRC(61)/0.600000)))
ALGBRC(65) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log((CONSTS(5)+ CONSTS(38)*CONSTS(6))/(CONSTS(9)+ CONSTS(38)*STATES(2)))
ALGBRC(66) = CONSTS(36)*ALGBRC(63)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(65))
ALGBRC(62) = - arbitrary_log(STATES(7)/1.00000+3.00000, 10)
ALGBRC(64) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP(- 7.20000+ALGBRC(62)/0.600000)))
ALGBRC(67) = CONSTS(37)*ALGBRC(64)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(65))
ALGBRC(68) = ALGBRC(66)+ALGBRC(67)
ALGBRC(72) = CONSTS(39)*STATES(14)*(STATES(13)+ 0.500000*ALGBRC(71))*(STATES(1) - CONSTS(151))
ALGBRC(73) = CONSTS(40)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(151))
ALGBRC(74) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(151)+59.2150))))
ALGBRC(75) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(151))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(151)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(151))+4.75300)))
ALGBRC(76) = ALGBRC(74)/(ALGBRC(74)+ALGBRC(75))
ALGBRC(77) = CONSTS(157)*ALGBRC(76)*(STATES(1) - CONSTS(151))
ALGBRC(78) = ( CONSTS(41)*(STATES(1) - CONSTS(151)))/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(97) = ALGBRC(94)+ALGBRC(95)
ALGBRC(79) = CONSTS(42)*(STATES(1) - CONSTS(152))*(CONSTS(44)/(1.00000+CONSTS(43)/STATES(6))+CONSTS(45)/(1.00000+CONSTS(43)/STATES(7)))
ALGBRC(80) = CONSTS(46)*(STATES(1) - CONSTS(152))
ALGBRC(83) = ( ALGBRC(81)*ALGBRC(82)*CONSTS(51)*CONSTS(47)*4.00000*( CONSTS(52)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(53)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(86) = ( ALGBRC(81)*ALGBRC(85)*CONSTS(50)*CONSTS(47)*4.00000*( CONSTS(52)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(53)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(88) = ( ALGBRC(81)*( ALGBRC(82)*CONSTS(51)+ ALGBRC(85)*CONSTS(50))*CONSTS(49)*( CONSTS(56)*CONSTS(9)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(57)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(89) = ALGBRC(83)+ALGBRC(86)+ALGBRC(84)+ALGBRC(87)+ALGBRC(88)
ALGBRC(45) = (( CONSTS(13)*CONSTS(14))/( 2.00000*CONSTS(15)))*log(CONSTS(7)/STATES(6))
ALGBRC(104) = CONSTS(78)*CONSTS(79)*(STATES(1) - ALGBRC(45))
ALGBRC(47) = (( CONSTS(13)*CONSTS(14))/( 2.00000*CONSTS(15)))*log(CONSTS(7)/STATES(7))
ALGBRC(105) = CONSTS(78)*CONSTS(80)*(STATES(1) - ALGBRC(47))
ALGBRC(106) = ALGBRC(105)+ALGBRC(104)
ALGBRC(102) = ( CONSTS(160)*CONSTS(77)*CONSTS(73))/(1.00000+CONSTS(75)/STATES(7) ** CONSTS(76))
ALGBRC(100) = ( CONSTS(160)*CONSTS(77)*CONSTS(72))/(1.00000+CONSTS(75)/STATES(6) ** CONSTS(76))
ALGBRC(103) = ALGBRC(100)+ALGBRC(102)
ALGBRC(7) = TERNRY(VOI - INT(VOI/CONSTS(2))*CONSTS(2).GE.CONSTS(1).AND.VOI - INT(VOI/CONSTS(2))*CONSTS(2).LE.CONSTS(1)+CONSTS(3), - CONSTS(4), 0.00000)
RATES(1) = - (ALGBRC(51)+ALGBRC(54)+ALGBRC(58)+ALGBRC(60)+ALGBRC(68)+ALGBRC(72)+ALGBRC(73)+ALGBRC(77)+ALGBRC(78)+ALGBRC(97)+ALGBRC(79)+ALGBRC(80)+ALGBRC(89)+ALGBRC(106)+ALGBRC(103)+ALGBRC(7))
ALGBRC(110) = CONSTS(114)*STATES(25)*(( CONSTS(112)*CONSTS(175))/CONSTS(172) - STATES(32)) - CONSTS(117)*STATES(32)
RATES(32) = ALGBRC(110)
ALGBRC(107) = CONSTS(84)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(108) = CONSTS(90)*(STATES(25) - STATES(6))
ALGBRC(109) = ( (( CONSTS(161)*CONSTS(91)*CONSTS(175))/CONSTS(172))*(STATES(5)/CONSTS(93) ** CONSTS(95) - STATES(25)/CONSTS(94) ** CONSTS(95)))/(1.00000+STATES(5)/CONSTS(93) ** CONSTS(95)+STATES(25)/CONSTS(94) ** CONSTS(95))
RATES(25) = (ALGBRC(109) - (( ALGBRC(108)*CONSTS(175))/CONSTS(172)+ALGBRC(107))) - ALGBRC(110)
ALGBRC(113) = CONSTS(113)*STATES(6)*(( CONSTS(109)*CONSTS(175))/CONSTS(174) - STATES(29)) - CONSTS(115)*STATES(29)
RATES(29) = ALGBRC(113)
ALGBRC(115) = CONSTS(113)*STATES(7)*(( CONSTS(108)*CONSTS(175))/CONSTS(173) - STATES(28)) - CONSTS(115)*STATES(28)
RATES(28) = ALGBRC(115)
ALGBRC(116) = CONSTS(113)*STATES(6)*(( CONSTS(111)*CONSTS(175))/CONSTS(174) - STATES(31)) - CONSTS(116)*STATES(31)
RATES(31) = ALGBRC(116)
ALGBRC(114) = CONSTS(125)*STATES(5)*(CONSTS(118) - STATES(33)) - CONSTS(132)*STATES(33)
RATES(33) = ALGBRC(114)
ALGBRC(118) = CONSTS(113)*STATES(7)*(( CONSTS(110)*CONSTS(175))/CONSTS(173) - STATES(30)) - CONSTS(116)*STATES(30)
RATES(30) = ALGBRC(118)
ALGBRC(117) = CONSTS(126)*STATES(5)*(CONSTS(119) - (STATES(34)+STATES(35))) - CONSTS(133)*STATES(34)
RATES(34) = ALGBRC(117)
ALGBRC(111) = (STATES(6) - STATES(7))*CONSTS(164)
ALGBRC(119) = ALGBRC(113)+ALGBRC(116)
ALGBRC(121) = (ALGBRC(86) - 2.00000*ALGBRC(94))+ALGBRC(104)+ALGBRC(100)
RATES(6) = ((( - ALGBRC(121)*CONSTS(170))/( CONSTS(174)*2.00000*CONSTS(15)) - ALGBRC(111)/CONSTS(174))+( ALGBRC(107)*CONSTS(172))/CONSTS(174)+( ALGBRC(108)*CONSTS(175))/CONSTS(174)) - 1.00000*ALGBRC(119)
ALGBRC(120) = CONSTS(127)*CONSTS(10)*(CONSTS(120) - (STATES(34)+STATES(35))) - CONSTS(134)*STATES(35)
RATES(35) = ALGBRC(120)
ALGBRC(112) = (STATES(7) - STATES(5))*CONSTS(165)
ALGBRC(122) = ALGBRC(115)+ALGBRC(118)
ALGBRC(124) = (ALGBRC(83) - 2.00000*ALGBRC(95))+ALGBRC(105)+ALGBRC(102)
RATES(7) = (( - ALGBRC(124)*CONSTS(170))/( CONSTS(173)*2.00000*CONSTS(15))+(ALGBRC(111) - ALGBRC(112))/CONSTS(173)) - 1.00000*ALGBRC(122)
ALGBRC(123) = CONSTS(128)*STATES(5)*(CONSTS(121) - STATES(36)) - CONSTS(135)*STATES(36)
RATES(36) = ALGBRC(123)
ALGBRC(125) = CONSTS(129)*STATES(5)*(CONSTS(122) - (STATES(37)+STATES(38))) - CONSTS(136)*STATES(37)
RATES(37) = ALGBRC(125)
ALGBRC(126) = CONSTS(130)*CONSTS(10)*(CONSTS(123) - (STATES(37)+STATES(38))) - CONSTS(137)*STATES(38)
RATES(38) = ALGBRC(126)
ALGBRC(127) = CONSTS(131)*STATES(5)*(CONSTS(124) - STATES(39)) - CONSTS(138)*STATES(39)
RATES(39) = ALGBRC(127)
ALGBRC(128) = ALGBRC(114)+ALGBRC(117)+ALGBRC(120)+ALGBRC(123)+ALGBRC(125)+ALGBRC(126)+ALGBRC(127)
RATES(5) = (( - ALGBRC(109)*CONSTS(172))/CONSTS(175)+ALGBRC(112)/CONSTS(175)) - 1.00000*ALGBRC(128)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (50.0000+STATES(1))/7.50000))
ALGBRC(21) = 1.00000/(( 0.00138000*(STATES(1)+7.00000))/(1.00000 - EXP( - 0.123000*(STATES(1)+7.00000)))+( 0.000610000*(STATES(1)+10.0000))/(EXP( 0.145000*(STATES(1)+10.0000)) - 1.00000))
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(23) = 1.00000/(( 7.19000e-05*(STATES(1)+30.0000))/(1.00000 - EXP( - 0.148000*(STATES(1)+30.0000)))+( 0.000131000*(STATES(1)+30.0000))/(- 1.00000+EXP( 0.0687000*(STATES(1)+30.0000))))
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(24) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(25) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
ALGBRC(11) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(26) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
ALGBRC(12) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(27) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
ALGBRC(13) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(28) = ( 1.00000*ALGBRC(13)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
ALGBRC(14) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(29) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
ALGBRC(16) = STATES(40)+( STATES(42)*CONSTS(174))/CONSTS(175)+( STATES(41)*CONSTS(173))/CONSTS(175)
ALGBRC(31) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(5)*(CONSTS(141) - STATES(40)) - CONSTS(149)*STATES(40)
ALGBRC(33) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(7)*(( CONSTS(142)*CONSTS(175))/CONSTS(173) - STATES(41)) - CONSTS(149)*STATES(41)
ALGBRC(32) = CONSTS(147)*(CONSTS(139) - ALGBRC(16))*STATES(6)*(( CONSTS(143)*CONSTS(175))/CONSTS(174) - STATES(42)) - CONSTS(149)*STATES(42)
ALGBRC(17) = STATES(43)+( STATES(45)*CONSTS(174))/CONSTS(175)+( STATES(44)*CONSTS(173))/CONSTS(175)
ALGBRC(34) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(5)*(CONSTS(144) - STATES(43)) - CONSTS(150)*STATES(43)
ALGBRC(36) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(7)*(( CONSTS(145)*CONSTS(175))/CONSTS(173) - STATES(44)) - CONSTS(150)*STATES(44)
ALGBRC(35) = CONSTS(148)*(CONSTS(140) - ALGBRC(17))*STATES(6)*(( CONSTS(146)*CONSTS(175))/CONSTS(174) - STATES(45)) - CONSTS(150)*STATES(45)
ALGBRC(15) = CONSTS(81) - (CONSTS(81) - CONSTS(82))/(1.00000+CONSTS(83)/STATES(25) ** CONSTS(89))
ALGBRC(30) = CONSTS(85)/ALGBRC(15)
ALGBRC(41) = CONSTS(87)*ALGBRC(15)
ALGBRC(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(18) = TERNRY(STATES(1).LT.- 40.0000, 3.56000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.350000*STATES(1)), 1.00000/( 0.130000*(1.00000+EXP((STATES(1)+10.6600)/- 11.1000))))
ALGBRC(37) = 1.00000/(ALGBRC(1)+ALGBRC(18))
ALGBRC(42) = ALGBRC(1)/(ALGBRC(1)+ALGBRC(18))
ALGBRC(2) = TERNRY(STATES(1).LT.- 40.0000, (( ( - 127140.*EXP( 0.244400*STATES(1)) - 3.47400e-05*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/1.00000)/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))), 0.00000)
ALGBRC(19) = TERNRY(STATES(1).LT.- 40.0000, ( 0.121200*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))), ( 0.300000*EXP( - 2.53500e-07*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))))
ALGBRC(38) = 1.00000/(ALGBRC(2)+ALGBRC(19))
ALGBRC(43) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(19))
ALGBRC(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(20) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(39) = 1.00000/(ALGBRC(3)+ALGBRC(20))
ALGBRC(44) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(20))
ALGBRC(46) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
ALGBRC(69) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(71) = ALGBRC(69)
ALGBRC(70) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
ALGBRC(98) = CONSTS(106)*STATES(3)*(CONSTS(105) - STATES(27)) - CONSTS(107)*STATES(27)
ALGBRC(22) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(6)/STATES(3))
ALGBRC(48) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(49) = CONSTS(23)*CONSTS(22)*ALGBRC(48)*(STATES(1) - ALGBRC(22))
ALGBRC(52) = CONSTS(26)*CONSTS(25)*(STATES(1) - ALGBRC(22))
ALGBRC(55) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))+ 0.0365000*CONSTS(155)*EXP(( - STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(56) = ( (( CONSTS(34)*CONSTS(153)*CONSTS(33)*ALGBRC(55))/(1.00000+( CONSTS(154)*CONSTS(31))/STATES(3) ** CONSTS(29)))*CONSTS(5))/(CONSTS(5)+CONSTS(32))
ALGBRC(85) = 1.00000 - STATES(21)
ALGBRC(81) = ( STATES(18)*STATES(19)*CONSTS(158)*STATES(1)*CONSTS(15) ** 2.00000)/( CONSTS(13)*CONSTS(14))
ALGBRC(87) = ( ALGBRC(81)*ALGBRC(85)*CONSTS(50)*CONSTS(48)*( CONSTS(54)*STATES(3)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(55)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(93) = 1.00000/(1.00000+CONSTS(67)/STATES(6) ** 3.00000)
ALGBRC(90) = ( EXP(( CONSTS(69)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*STATES(3) ** CONSTS(70)*CONSTS(7) - EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*CONSTS(6) ** CONSTS(70)*STATES(6))/(1.00000+ CONSTS(68)*EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(94) = ( CONSTS(60)*CONSTS(59)*ALGBRC(93)*CONSTS(159)*ALGBRC(90))/( CONSTS(66)*CONSTS(6) ** CONSTS(70)*(1.00000+STATES(3)/CONSTS(63) ** CONSTS(70))+ CONSTS(65) ** CONSTS(70)*STATES(6)*(1.00000+STATES(6)/CONSTS(66))+ CONSTS(64)*STATES(3) ** CONSTS(70)+ STATES(3) ** CONSTS(70)*CONSTS(7)+ CONSTS(6) ** CONSTS(70)*STATES(6))
ALGBRC(96) = (STATES(3) - STATES(4))*CONSTS(162)
ALGBRC(99) = (STATES(4) - STATES(2))*CONSTS(163)
ALGBRC(101) = CONSTS(106)*STATES(4)*(CONSTS(104) - STATES(26)) - CONSTS(107)*STATES(26)
ALGBRC(40) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log(CONSTS(6)/STATES(4))
ALGBRC(50) = CONSTS(24)*CONSTS(22)*ALGBRC(48)*(STATES(1) - ALGBRC(40))
ALGBRC(53) = CONSTS(27)*CONSTS(25)*(STATES(1) - ALGBRC(40))
ALGBRC(57) = ( (( CONSTS(35)*CONSTS(153)*CONSTS(33)*ALGBRC(55))/(1.00000+( CONSTS(154)*CONSTS(31))/STATES(4) ** CONSTS(29)))*CONSTS(5))/(CONSTS(5)+CONSTS(32))
ALGBRC(82) = 1.00000 - STATES(20)
ALGBRC(84) = ( ALGBRC(81)*ALGBRC(82)*CONSTS(51)*CONSTS(48)*( CONSTS(54)*STATES(4)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(55)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(92) = 1.00000/(1.00000+CONSTS(67)/STATES(7) ** 3.00000)
ALGBRC(91) = ( EXP(( CONSTS(69)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*STATES(4) ** CONSTS(70)*CONSTS(7) - EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14)))*CONSTS(6) ** CONSTS(70)*STATES(7))/(1.00000+ CONSTS(68)*EXP(( (CONSTS(69) - 1.00000)*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))))
ALGBRC(95) = ( CONSTS(61)*CONSTS(59)*ALGBRC(92)*CONSTS(159)*ALGBRC(91))/( CONSTS(66)*CONSTS(6) ** CONSTS(70)*(1.00000+STATES(4)/CONSTS(63) ** CONSTS(70))+ CONSTS(65) ** CONSTS(70)*STATES(7)*(1.00000+STATES(7)/CONSTS(66))+ CONSTS(64)*STATES(4) ** CONSTS(70)+ STATES(4) ** CONSTS(70)*CONSTS(7)+ CONSTS(6) ** CONSTS(70)*STATES(7))
ALGBRC(51) = ALGBRC(49)+ALGBRC(50)
ALGBRC(54) = ALGBRC(52)+ALGBRC(53)
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(59) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(60) = CONSTS(156)*STATES(11)*ALGBRC(59)*(STATES(1) - CONSTS(151))
ALGBRC(61) = - arbitrary_log(STATES(6)/1.00000+3.00000, 10)
ALGBRC(63) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP(- 7.20000+ALGBRC(61)/0.600000)))
ALGBRC(65) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log((CONSTS(5)+ CONSTS(38)*CONSTS(6))/(CONSTS(9)+ CONSTS(38)*STATES(2)))
ALGBRC(66) = CONSTS(36)*ALGBRC(63)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(65))
ALGBRC(62) = - arbitrary_log(STATES(7)/1.00000+3.00000, 10)
ALGBRC(64) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP(- 7.20000+ALGBRC(62)/0.600000)))
ALGBRC(67) = CONSTS(37)*ALGBRC(64)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(65))
ALGBRC(68) = ALGBRC(66)+ALGBRC(67)
ALGBRC(72) = CONSTS(39)*STATES(14)*(STATES(13)+ 0.500000*ALGBRC(71))*(STATES(1) - CONSTS(151))
ALGBRC(73) = CONSTS(40)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(151))
ALGBRC(74) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(151)+59.2150))))
ALGBRC(75) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(151))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(151)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(151))+4.75300)))
ALGBRC(76) = ALGBRC(74)/(ALGBRC(74)+ALGBRC(75))
ALGBRC(77) = CONSTS(157)*ALGBRC(76)*(STATES(1) - CONSTS(151))
ALGBRC(78) = ( CONSTS(41)*(STATES(1) - CONSTS(151)))/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(97) = ALGBRC(94)+ALGBRC(95)
ALGBRC(79) = CONSTS(42)*(STATES(1) - CONSTS(152))*(CONSTS(44)/(1.00000+CONSTS(43)/STATES(6))+CONSTS(45)/(1.00000+CONSTS(43)/STATES(7)))
ALGBRC(80) = CONSTS(46)*(STATES(1) - CONSTS(152))
ALGBRC(83) = ( ALGBRC(81)*ALGBRC(82)*CONSTS(51)*CONSTS(47)*4.00000*( CONSTS(52)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(53)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(86) = ( ALGBRC(81)*ALGBRC(85)*CONSTS(50)*CONSTS(47)*4.00000*( CONSTS(52)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(53)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(88) = ( ALGBRC(81)*( ALGBRC(82)*CONSTS(51)+ ALGBRC(85)*CONSTS(50))*CONSTS(49)*( CONSTS(56)*CONSTS(9)*EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - CONSTS(57)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(15))/( CONSTS(13)*CONSTS(14))) - 1.00000)
ALGBRC(89) = ALGBRC(83)+ALGBRC(86)+ALGBRC(84)+ALGBRC(87)+ALGBRC(88)
ALGBRC(45) = (( CONSTS(13)*CONSTS(14))/( 2.00000*CONSTS(15)))*log(CONSTS(7)/STATES(6))
ALGBRC(104) = CONSTS(78)*CONSTS(79)*(STATES(1) - ALGBRC(45))
ALGBRC(47) = (( CONSTS(13)*CONSTS(14))/( 2.00000*CONSTS(15)))*log(CONSTS(7)/STATES(7))
ALGBRC(105) = CONSTS(78)*CONSTS(80)*(STATES(1) - ALGBRC(47))
ALGBRC(106) = ALGBRC(105)+ALGBRC(104)
ALGBRC(102) = ( CONSTS(160)*CONSTS(77)*CONSTS(73))/(1.00000+CONSTS(75)/STATES(7) ** CONSTS(76))
ALGBRC(100) = ( CONSTS(160)*CONSTS(77)*CONSTS(72))/(1.00000+CONSTS(75)/STATES(6) ** CONSTS(76))
ALGBRC(103) = ALGBRC(100)+ALGBRC(102)
ALGBRC(7) = TERNRY(VOI - INT(VOI/CONSTS(2))*CONSTS(2).GE.CONSTS(1).AND.VOI - INT(VOI/CONSTS(2))*CONSTS(2).LE.CONSTS(1)+CONSTS(3), - CONSTS(4), 0.00000)
ALGBRC(110) = CONSTS(114)*STATES(25)*(( CONSTS(112)*CONSTS(175))/CONSTS(172) - STATES(32)) - CONSTS(117)*STATES(32)
ALGBRC(107) = CONSTS(84)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(108) = CONSTS(90)*(STATES(25) - STATES(6))
ALGBRC(109) = ( (( CONSTS(161)*CONSTS(91)*CONSTS(175))/CONSTS(172))*(STATES(5)/CONSTS(93) ** CONSTS(95) - STATES(25)/CONSTS(94) ** CONSTS(95)))/(1.00000+STATES(5)/CONSTS(93) ** CONSTS(95)+STATES(25)/CONSTS(94) ** CONSTS(95))
ALGBRC(113) = CONSTS(113)*STATES(6)*(( CONSTS(109)*CONSTS(175))/CONSTS(174) - STATES(29)) - CONSTS(115)*STATES(29)
ALGBRC(115) = CONSTS(113)*STATES(7)*(( CONSTS(108)*CONSTS(175))/CONSTS(173) - STATES(28)) - CONSTS(115)*STATES(28)
ALGBRC(116) = CONSTS(113)*STATES(6)*(( CONSTS(111)*CONSTS(175))/CONSTS(174) - STATES(31)) - CONSTS(116)*STATES(31)
ALGBRC(114) = CONSTS(125)*STATES(5)*(CONSTS(118) - STATES(33)) - CONSTS(132)*STATES(33)
ALGBRC(118) = CONSTS(113)*STATES(7)*(( CONSTS(110)*CONSTS(175))/CONSTS(173) - STATES(30)) - CONSTS(116)*STATES(30)
ALGBRC(117) = CONSTS(126)*STATES(5)*(CONSTS(119) - (STATES(34)+STATES(35))) - CONSTS(133)*STATES(34)
ALGBRC(111) = (STATES(6) - STATES(7))*CONSTS(164)
ALGBRC(119) = ALGBRC(113)+ALGBRC(116)
ALGBRC(121) = (ALGBRC(86) - 2.00000*ALGBRC(94))+ALGBRC(104)+ALGBRC(100)
ALGBRC(120) = CONSTS(127)*CONSTS(10)*(CONSTS(120) - (STATES(34)+STATES(35))) - CONSTS(134)*STATES(35)
ALGBRC(112) = (STATES(7) - STATES(5))*CONSTS(165)
ALGBRC(122) = ALGBRC(115)+ALGBRC(118)
ALGBRC(124) = (ALGBRC(83) - 2.00000*ALGBRC(95))+ALGBRC(105)+ALGBRC(102)
ALGBRC(123) = CONSTS(128)*STATES(5)*(CONSTS(121) - STATES(36)) - CONSTS(135)*STATES(36)
ALGBRC(125) = CONSTS(129)*STATES(5)*(CONSTS(122) - (STATES(37)+STATES(38))) - CONSTS(136)*STATES(37)
ALGBRC(126) = CONSTS(130)*CONSTS(10)*(CONSTS(123) - (STATES(37)+STATES(38))) - CONSTS(137)*STATES(38)
ALGBRC(127) = CONSTS(131)*STATES(5)*(CONSTS(124) - STATES(39)) - CONSTS(138)*STATES(39)
ALGBRC(128) = ALGBRC(114)+ALGBRC(117)+ALGBRC(120)+ALGBRC(123)+ALGBRC(125)+ALGBRC(126)+ALGBRC(127)
ALGBRC(5) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log((CONSTS(5)+ CONSTS(38)*CONSTS(6))/(CONSTS(9)+ CONSTS(38)*STATES(3)))
ALGBRC(6) = (( CONSTS(13)*CONSTS(14))/CONSTS(15))*log((CONSTS(5)+ CONSTS(38)*CONSTS(6))/(CONSTS(9)+ CONSTS(38)*STATES(4)))
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END