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 131 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 164 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(58) is i_Na in component INa (microA_per_microF).
C ALGBRC(61) is i_Nab in component INab (microA_per_microF).
C ALGBRC(65) is i_NaK in component INaK (microA_per_microF).
C ALGBRC(67) is i_Kr in component IKr (microA_per_microF).
C ALGBRC(75) is i_Ks in component IKs (microA_per_microF).
C ALGBRC(76) is i_tos in component Itos (microA_per_microF).
C ALGBRC(77) is i_tof in component Itof (microA_per_microF).
C ALGBRC(81) is i_K1 in component IK1 (microA_per_microF).
C ALGBRC(100) is i_NaCa in component INaCa (microA_per_microF).
C ALGBRC(82) is i_Cl_Ca in component ICl_Ca (microA_per_microF).
C ALGBRC(83) is i_Clb in component IClb (microA_per_microF).
C ALGBRC(92) is i_CaL in component ICaL (microA_per_microF).
C ALGBRC(109) is i_Cab in component ICab (microA_per_microF).
C ALGBRC(106) 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 Cli in component model_parameters (millimolar).
C CONSTS(12) is R in component model_parameters (joule_per_kilomole_kelvin).
C CONSTS(13) is T in component model_parameters (kelvin).
C CONSTS(14) is F in component model_parameters (coulomb_per_mole).
C CONSTS(157) is Cm in component model_parameters (farad).
C CONSTS(15) is Cm_per_area in component model_parameters (farad_per_cm2).
C CONSTS(16) is cell_length in component model_parameters (micrometre).
C CONSTS(17) is cell_radius in component model_parameters (micrometre).
C CONSTS(158) is Vol_Cell in component model_parameters (litre).
C CONSTS(159) is Vol_SR in component model_parameters (litre).
C CONSTS(160) is Vol_SL in component model_parameters (litre).
C CONSTS(161) is Vol_jct in component model_parameters (litre).
C CONSTS(162) is Vol_cytosol in component model_parameters (litre).
C CONSTS(163) is Vol_mito in component model_parameters (litre).
C ALGBRC(27) is E_Na_jct in component reversal_potentials (millivolt).
C ALGBRC(47) is E_Na_SL in component reversal_potentials (millivolt).
C ALGBRC(52) is E_Ca_jct in component reversal_potentials (millivolt).
C ALGBRC(54) is E_Ca_SL in component reversal_potentials (millivolt).
C CONSTS(146) is E_K in component reversal_potentials (millivolt).
C CONSTS(147) 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(56) is i_Na_jct in component INa (microA_per_microF).
C ALGBRC(57) is i_Na_SL in component INa (microA_per_microF).
C CONSTS(18) is G_INa in component INa (milliS_per_microF).
C CONSTS(19) is Fx_Na_jct in component INa (dimensionless).
C CONSTS(20) 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(55) is openProb in component INa (dimensionless).
C ALGBRC(1) is alpha_h in component INa_h_gate (per_millisecond).
C ALGBRC(23) is beta_h in component INa_h_gate (per_millisecond).
C ALGBRC(44) is tau_h in component INa_h_gate (millisecond).
C ALGBRC(49) is h_infinity in component INa_h_gate (dimensionless).
C ALGBRC(2) is alpha_j in component INa_j_gate (per_millisecond).
C ALGBRC(24) is beta_j in component INa_j_gate (per_millisecond).
C ALGBRC(45) is tau_j in component INa_j_gate (millisecond).
C ALGBRC(50) is j_infinity in component INa_j_gate (dimensionless).
C ALGBRC(3) is alpha_m in component INa_m_gate (per_millisecond).
C ALGBRC(25) is beta_m in component INa_m_gate (per_millisecond).
C ALGBRC(46) is tau_m in component INa_m_gate (millisecond).
C ALGBRC(51) is m_infinity in component INa_m_gate (dimensionless).
C ALGBRC(59) is i_Nab_jct in component INab (microA_per_microF).
C ALGBRC(60) is i_Nab_SL in component INab (microA_per_microF).
C CONSTS(21) is G_NaBk in component INab (milliS_per_microF).
C CONSTS(22) is Fx_NaBk_jct in component INab (dimensionless).
C CONSTS(23) is Fx_NaBk_SL in component INab (dimensionless).
C ALGBRC(63) is i_NaK_jct in component INaK (microA_per_microF).
C ALGBRC(64) is i_NaK_SL in component INaK (microA_per_microF).
C ALGBRC(62) is f_NaK in component INaK (dimensionless).
C CONSTS(24) is Q10_NaK in component INaK (dimensionless).
C CONSTS(148) is Q_NaK in component INaK (dimensionless).
C CONSTS(25) is H_NaK in component INaK (dimensionless).
C CONSTS(26) is Q10_Km_Nai in component INaK (dimensionless).
C CONSTS(149) is Q_Km_Nai in component INaK (dimensionless).
C CONSTS(27) is Km_Nai in component INaK (millimolar).
C CONSTS(28) is Km_Ko in component INaK (millimolar).
C CONSTS(29) is I_NaK_max in component INaK (microA_per_microF).
C CONSTS(150) is sigma in component INaK (dimensionless).
C CONSTS(30) is Fx_NaK_jct in component INaK (dimensionless).
C CONSTS(31) is Fx_NaK_SL in component INaK (dimensionless).
C CONSTS(151) is G_IKr in component IKr (milliS_per_microF).
C STATES(11) is Xr in component IKr_Xr_gate (dimensionless).
C ALGBRC(66) is Rr in component IKr_Rr_gate (dimensionless).
C ALGBRC(4) is Xr_infinity in component IKr_Xr_gate (dimensionless).
C ALGBRC(26) is tau_Xr in component IKr_Xr_gate (millisecond).
C ALGBRC(73) is i_Ks_jct in component IKs (microA_per_microF).
C ALGBRC(74) is i_Ks_SL in component IKs (microA_per_microF).
C ALGBRC(71) is G_Ks_SL in component IKs (milliS_per_microF).
C ALGBRC(70) is G_Ks_jct in component IKs (milliS_per_microF).
C CONSTS(32) is Fx_Ks_jct in component IKs (dimensionless).
C CONSTS(33) 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(72) is E_Ks in component IKs (millivolt).
C CONSTS(34) is pKNa in component IKs (dimensionless).
C ALGBRC(68) is pCa_jct in component IKs (dimensionless).
C ALGBRC(69) 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(28) is tau_Xs in component IKs_Xs_gate (millisecond).
C CONSTS(35) 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 STATES(15) is R_tos in component Itos_R_gate (dimensionless).
C ALGBRC(9) is X_tos_infinity in component Itos_X_gate (dimensionless).
C ALGBRC(29) 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(30) is tau_Y_tos in component Itos_Y_gate (millisecond).
C ALGBRC(11) is R_tos_infinity in component Itos_R_gate (dimensionless).
C ALGBRC(31) is R_tos_other in component Itos_R_gate (dimensionless).
C ALGBRC(32) is tau_R_tos in component Itos_R_gate (millisecond).
C CONSTS(36) 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(12) is X_tof_infinity in component Itof_X_gate (dimensionless).
C ALGBRC(33) is tau_X_tof in component Itof_X_gate (millisecond).
C ALGBRC(13) is Y_tof_infinity in component Itof_Y_gate (dimensionless).
C ALGBRC(34) is tau_Y_tof in component Itof_Y_gate (millisecond).
C CONSTS(152) is G_K1 in component IK1 (milliS_per_microF).
C ALGBRC(80) is K1_infinity in component IK1_K1_gate (dimensionless).
C ALGBRC(78) is alpha_K1 in component IK1_K1_gate (per_millisecond).
C ALGBRC(79) is beta_K1 in component IK1_K1_gate (per_millisecond).
C CONSTS(37) is G_Cl in component ICl_Ca (milliS_per_microF).
C CONSTS(38) is Kd_ClCa in component ICl_Ca (millimolar).
C CONSTS(39) is Fx_Cl_jct in component ICl_Ca (dimensionless).
C CONSTS(40) is Fx_Cl_SL in component ICl_Ca (dimensionless).
C CONSTS(41) is G_ClBk in component IClb (milliS_per_microF).
C ALGBRC(86) is i_CaL_Ca_SL in component ICaL (microA_per_microF).
C ALGBRC(89) is i_CaL_Ca_jct in component ICaL (microA_per_microF).
C ALGBRC(87) is i_CaL_Na_SL in component ICaL (microA_per_microF).
C ALGBRC(90) is i_CaL_Na_jct in component ICaL (microA_per_microF).
C ALGBRC(91) is i_CaL_K in component ICaL (microA_per_microF).
C CONSTS(42) is PCa in component ICaL (litre_per_farad_millisecond).
C CONSTS(43) is PNa in component ICaL (litre_per_farad_millisecond).
C CONSTS(44) is PK in component ICaL (litre_per_farad_millisecond).
C CONSTS(45) is Fx_ICaL_jct in component ICaL (dimensionless).
C CONSTS(46) is Fx_ICaL_SL in component ICaL (dimensionless).
C CONSTS(47) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(48) is gamma_Cao in component ICaL (dimensionless).
C CONSTS(49) is gamma_Nai in component ICaL (dimensionless).
C CONSTS(50) is gamma_Nao in component ICaL (dimensionless).
C CONSTS(51) is gamma_Ki in component ICaL (dimensionless).
C CONSTS(52) is gamma_Ko in component ICaL (dimensionless).
C CONSTS(53) is Q10_CaL in component ICaL (dimensionless).
C CONSTS(153) 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(85) is fCa_SL in component ICaL_fCa_gate (dimensionless).
C ALGBRC(88) is fCa_jct in component ICaL_fCa_gate (dimensionless).
C ALGBRC(84) is temp in component ICaL (coulomb_per_mole).
C ALGBRC(14) is d_infinity in component ICaL_d_gate (dimensionless).
C ALGBRC(35) is tau_d in component ICaL_d_gate (millisecond).
C ALGBRC(15) is f_infinity in component ICaL_f_gate (dimensionless).
C ALGBRC(36) 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(98) is i_NaCa_SL in component INaCa (microA_per_microF).
C ALGBRC(97) is i_NaCa_jct in component INaCa (microA_per_microF).
C CONSTS(54) is V_max in component INaCa (microA_per_microF).
C CONSTS(55) is Fx_NCX_jct in component INaCa (dimensionless).
C CONSTS(56) is Fx_NCX_SL in component INaCa (dimensionless).
C CONSTS(57) is Q10_NCX in component INaCa (dimensionless).
C CONSTS(154) is Q_NCX in component INaCa (dimensionless).
C CONSTS(58) is K_mNai in component INaCa (millimolar).
C CONSTS(59) is K_mCao in component INaCa (millimolar).
C CONSTS(60) is K_mNao in component INaCa (millimolar).
C CONSTS(61) is K_mCai in component INaCa (millimolar).
C ALGBRC(95) is Ka_SL in component INaCa (dimensionless).
C ALGBRC(96) is Ka_jct in component INaCa (dimensionless).
C CONSTS(62) is Kd_act in component INaCa (millimolar).
C CONSTS(63) is ksat in component INaCa (dimensionless).
C CONSTS(64) is eta in component INaCa (dimensionless).
C CONSTS(65) is HNa in component INaCa (dimensionless).
C ALGBRC(93) is temp_jct in component INaCa (millimolar4).
C ALGBRC(94) is temp_SL in component INaCa (millimolar4).
C ALGBRC(105) is i_Cap_SL in component ICap (microA_per_microF).
C ALGBRC(103) is i_Cap_jct in component ICap (microA_per_microF).
C CONSTS(66) is V_max in component ICap (millimolar_per_millisecond).
C CONSTS(67) is Fx_SLCaP_jct in component ICap (dimensionless).
C CONSTS(68) is Fx_SLCaP_SL in component ICap (dimensionless).
C CONSTS(69) is Q10_SLCaP in component ICap (dimensionless).
C CONSTS(155) is Q_SLCaP in component ICap (dimensionless).
C CONSTS(70) is Km in component ICap (millimolar).
C CONSTS(71) is H in component ICap (dimensionless).
C CONSTS(164) is V_maxAF_old in component ICap (microA_per_microF).
C CONSTS(72) is V_maxAF in component ICap (microA_per_microF).
C ALGBRC(107) is i_Cab_jct in component ICab (microA_per_microF).
C ALGBRC(108) is i_Cab_SL in component ICab (microA_per_microF).
C CONSTS(73) is G_CaBk in component ICab (milliS_per_microF).
C CONSTS(74) is Fx_CaBk_jct in component ICab (dimensionless).
C CONSTS(75) is Fx_CaBk_SL in component ICab (dimensionless).
C ALGBRC(110) is j_rel_SR in component Jrel_SR (millimolar_per_millisecond).
C CONSTS(76) is Max_SR in component Jrel_SR (dimensionless).
C CONSTS(77) is Min_SR in component Jrel_SR (dimensionless).
C CONSTS(78) 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(53) is RI in component Jrel_SR (dimensionless).
C STATES(24) is O in component Jrel_SR (dimensionless).
C CONSTS(79) is ks in component Jrel_SR (per_millisecond).
C CONSTS(80) is koCa in component Jrel_SR (per_millimolar2_per_millisecond).
C CONSTS(81) is kom in component Jrel_SR (per_millisecond).
C CONSTS(82) is kiCa in component Jrel_SR (per_millimolar_per_millisecond).
C CONSTS(83) is kim in component Jrel_SR (per_millisecond).
C CONSTS(84) is HSR in component Jrel_SR (dimensionless).
C STATES(25) is Ca_SR in component Ca_buffer (millimolar).
C ALGBRC(16) is kCaSR in component Jrel_SR (dimensionless).
C ALGBRC(37) is koSRCa in component Jrel_SR (per_millimolar2_per_millisecond).
C ALGBRC(48) is kiSRCa in component Jrel_SR (per_millimolar_per_millisecond).
C ALGBRC(111) is j_leak_SR in component Jleak_SR (millimolar_per_millisecond).
C CONSTS(85) is KSRleak in component Jleak_SR (per_millisecond).
C ALGBRC(112) is j_pump_SR in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(86) is V_max in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(87) is Q10_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(156) is Q_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(88) is Kmf in component Jpump_SR (millimolar).
C CONSTS(89) is Kmr in component Jpump_SR (millimolar).
C CONSTS(90) is H in component Jpump_SR (dimensionless).
C ALGBRC(99) is J_Na_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(102) is J_Na_SL_cytosol in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(114) is J_Ca_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(115) is J_Ca_SL_cytosol in component ion_diffusion (millimole_per_millisecond).
C CONSTS(91) is D_Ca_jct_SL in component ion_diffusion (dm2_per_second).
C CONSTS(92) is D_Ca_SL_cytosol in component ion_diffusion (dm2_per_second).
C CONSTS(93) is D_Na_jct_SL in component ion_diffusion (dm2_per_second).
C CONSTS(94) is D_Na_SL_cytosol in component ion_diffusion (dm2_per_second).
C CONSTS(95) is x_jct_SL in component ion_diffusion (micrometre).
C CONSTS(96) is x_SL_cytosol in component ion_diffusion (micrometre).
C CONSTS(97) is A_jct_SL in component ion_diffusion (cm2).
C CONSTS(98) is A_SL_cytosol in component ion_diffusion (cm2).
C ALGBRC(17) is J_Na_jct_SL_old in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(18) is J_Na_SL_cytosol_old in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(19) is J_Ca_jct_SL_old in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(20) is J_Ca_SL_cytosol_old in component ion_diffusion (millimole_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(99) is Bmax_SL in component Na_buffer (millimolar).
C CONSTS(100) is Bmax_jct in component Na_buffer (millimolar).
C CONSTS(101) is kon in component Na_buffer (per_millimolar_per_millisecond).
C CONSTS(102) is koff in component Na_buffer (per_millisecond).
C ALGBRC(101) is dNa_jct_buf in component Na_buffer (millimolar_per_millisecond).
C ALGBRC(104) 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(103) is Bmax_SLB_SL in component Ca_buffer (millimolar).
C CONSTS(104) is Bmax_SLB_jct in component Ca_buffer (millimolar).
C CONSTS(105) is Bmax_SLHigh_SL in component Ca_buffer (millimolar).
C CONSTS(106) is Bmax_SLHigh_jct in component Ca_buffer (millimolar).
C CONSTS(107) is Bmax_Calsequestrin in component Ca_buffer (millimolar).
C CONSTS(108) is kon_SL in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(109) is kon_Calsequestrin in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(110) is koff_SLB in component Ca_buffer (per_millisecond).
C CONSTS(111) is koff_SLHigh in component Ca_buffer (per_millisecond).
C CONSTS(112) is koff_Calsequestrin in component Ca_buffer (per_millisecond).
C ALGBRC(113) is dCalsequestrin in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(131) is dCa_cytosol_tot_bound in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(118) is dCa_SLB_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(116) is dCa_SLB_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(121) is dCa_SLHigh_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(119) is dCa_SLHigh_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(122) is dCa_jct_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(125) is dCa_SL_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(124) is i_Ca_jct_tot in component Ca_buffer (microA_per_microF).
C ALGBRC(127) 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(113) is Bmax_TroponinC in component cytosolic_Ca_buffer (millimolar).
C CONSTS(114) is Bmax_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(115) is Bmax_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(116) is Bmax_Calmodulin in component cytosolic_Ca_buffer (millimolar).
C CONSTS(117) is Bmax_Myosin_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(118) is Bmax_Myosin_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(119) is Bmax_SRB in component cytosolic_Ca_buffer (millimolar).
C CONSTS(120) is kon_TroponinC in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(121) is kon_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(122) is kon_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(123) is kon_Calmodulin in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(124) is kon_Myosin_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(125) is kon_Myosin_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(126) is kon_SRB in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(127) is koff_TroponinC in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(128) is koff_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(129) is koff_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(130) is koff_Calmodulin in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(131) is koff_Myosin_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(132) is koff_Myosin_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(133) is koff_SRB in component cytosolic_Ca_buffer (per_millisecond).
C ALGBRC(117) is dCa_TroponinC in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(120) is dCa_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(123) is dMg_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(126) is dCa_Calmodulin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(128) is dCa_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(129) is dMg_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(130) is dCa_SRB in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C CONSTS(134) is Indo1 in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(135) 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(136) is Bmax_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(137) is Bmax_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(138) is Bmax_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(139) is Bmax_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(140) is Bmax_Fluo3_SL in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(141) is Bmax_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar).
C CONSTS(142) is kon_Indo1 in component indo_fluo_Ca_buffer_not_connected (per_millimolar2_per_millisecond).
C CONSTS(143) is kon_Fluo3 in component indo_fluo_Ca_buffer_not_connected (per_millimolar2_per_millisecond).
C CONSTS(144) is koff_Indo1 in component indo_fluo_Ca_buffer_not_connected (per_millisecond).
C CONSTS(145) is koff_Fluo3 in component indo_fluo_Ca_buffer_not_connected (per_millisecond).
C ALGBRC(21) is Indo1Bound in component indo_fluo_Ca_buffer_not_connected (millimolar).
C ALGBRC(22) is Fluo3Bound in component indo_fluo_Ca_buffer_not_connected (millimolar).
C ALGBRC(38) is dCa_Indo1_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(39) is dCa_Indo1_jct in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(40) is dCa_Indo1_SL in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(41) is dCa_Fluo3_Cytosol in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(42) is dCa_Fluo3_jct in component indo_fluo_Ca_buffer_not_connected (millimolar_per_millisecond).
C ALGBRC(43) 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 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) = 15
CONSTS(12) = 8314.3
CONSTS(13) = 310
CONSTS(14) = 96486.7
CONSTS(15) = 2e-6
CONSTS(16) = 100
CONSTS(17) = 10.25
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(18) = 16
CONSTS(19) = 0.11
CONSTS(20) = 0.89
STATES(8) = 0.001370685156
STATES(9) = 0.987140350343
STATES(10) = 0.991822731369
CONSTS(21) = 0.297e-3
CONSTS(22) = 0.11
CONSTS(23) = 0.89
CONSTS(24) = 1.63
CONSTS(25) = 4
CONSTS(26) = 1.49
CONSTS(27) = 11
CONSTS(28) = 1.5
CONSTS(29) = 1.91
CONSTS(30) = 0.11
CONSTS(31) = 0.89
STATES(11) = 0.008471550841
CONSTS(32) = 0.11
CONSTS(33) = 0.89
CONSTS(34) = 0.01833
STATES(12) = 0.00687399199
CONSTS(35) = 0.06
STATES(13) = 0.293519921626
STATES(14) = 0.004011272375
STATES(15) = 0.383430556383
CONSTS(36) = 0.02
STATES(16) = 0.9946314893
STATES(17) = 0.00401120993
CONSTS(37) = 0.109625
CONSTS(38) = 0.1
CONSTS(39) = 0.11
CONSTS(40) = 0.89
CONSTS(41) = 0.009
CONSTS(42) = 5.4e-4
CONSTS(43) = 1.5e-8
CONSTS(44) = 2.7e-7
CONSTS(45) = 0.9
CONSTS(46) = 0.1
CONSTS(47) = 0.341
CONSTS(48) = 0.341
CONSTS(49) = 0.75
CONSTS(50) = 0.75
CONSTS(51) = 0.75
CONSTS(52) = 0.75
CONSTS(53) = 1.8
STATES(18) = 0.000006997531
STATES(19) = 1.000675515962
STATES(20) = 0.015352888928
STATES(21) = 0.024609183734
CONSTS(54) = 9
CONSTS(55) = 0.11
CONSTS(56) = 0.89
CONSTS(57) = 1.57
CONSTS(58) = 12.29
CONSTS(59) = 1.3
CONSTS(60) = 87.5
CONSTS(61) = 0.00359
CONSTS(62) = 0.000256
CONSTS(63) = 0.27
CONSTS(64) = 0.35
CONSTS(65) = 3
CONSTS(66) = 2.2e-6
CONSTS(67) = 0.11
CONSTS(68) = 0.89
CONSTS(69) = 2.35
CONSTS(70) = 0.0005
CONSTS(71) = 1.6
CONSTS(72) = 0.0673
CONSTS(73) = 0.0002513
CONSTS(74) = 0.11
CONSTS(75) = 0.89
CONSTS(76) = 15
CONSTS(77) = 1
CONSTS(78) = 0.45
STATES(22) = 0.884673513138
STATES(23) = 0.00000009272
STATES(24) = 0.000000711264
CONSTS(79) = 25
CONSTS(80) = 10
CONSTS(81) = 0.06
CONSTS(82) = 0.5
CONSTS(83) = 0.005
CONSTS(84) = 2.5
STATES(25) = 0.545611267699
CONSTS(85) = 5.348e-6
CONSTS(86) = 286e-6
CONSTS(87) = 2.6
CONSTS(88) = 0.000246
CONSTS(89) = 1.7
CONSTS(90) = 1.787
CONSTS(91) = 1.64e-6
CONSTS(92) = 1.22e-6
CONSTS(93) = 1.09e-5
CONSTS(94) = 1.79e-5
CONSTS(95) = 0.5
CONSTS(96) = 0.45
CONSTS(97) = 3.01e-6
CONSTS(98) = 1.3e-4
STATES(26) = 0.776121392467
STATES(27) = 3.557055389701
CONSTS(99) = 1.65
CONSTS(100) = 7.561
CONSTS(101) = 0.0001
CONSTS(102) = 1e-3
STATES(28) = 0.009868629147
STATES(29) = 0.007780801995
STATES(30) = 0.114438990328
STATES(31) = 0.077503874257
STATES(32) = 1.186496899338
CONSTS(103) = 0.0374
CONSTS(104) = 0.0046
CONSTS(105) = 0.0134
CONSTS(106) = 0.00165
CONSTS(107) = 0.14
CONSTS(108) = 100
CONSTS(109) = 100
CONSTS(110) = 1.3
CONSTS(111) = 30e-3
CONSTS(112) = 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(113) = 0.07
CONSTS(114) = 0.14
CONSTS(115) = 0.14
CONSTS(116) = 0.024
CONSTS(117) = 0.14
CONSTS(118) = 0.14
CONSTS(119) = 0.0171
CONSTS(120) = 32.7
CONSTS(121) = 2.37
CONSTS(122) = 3e-3
CONSTS(123) = 34
CONSTS(124) = 13.8
CONSTS(125) = 15.7e-3
CONSTS(126) = 100
CONSTS(127) = 19.6e-3
CONSTS(128) = 0.032e-3
CONSTS(129) = 3.33e-3
CONSTS(130) = 238e-3
CONSTS(131) = 0.46e-3
CONSTS(132) = 0.057e-3
CONSTS(133) = 60e-3
CONSTS(134) = 0
CONSTS(135) = 0
STATES(40) = 0
STATES(41) = 0
STATES(42) = 0
STATES(43) = 0
STATES(44) = 0
STATES(45) = 0
CONSTS(136) = 0.025
CONSTS(137) = 0.00077
CONSTS(138) = 0.00002
CONSTS(139) = 0.025
CONSTS(140) = 0.00077
CONSTS(141) = 0.00002
CONSTS(142) = 100
CONSTS(143) = 100
CONSTS(144) = 60e-3
CONSTS(145) = 110e-3
CONSTS(146) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(5)/CONSTS(9))
CONSTS(147) = (( - CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(8)/CONSTS(11))
CONSTS(148) = CONSTS(24) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(149) = CONSTS(26) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(150) = (EXP(CONSTS(6)/67.3000) - 1.00000)/7.00000
CONSTS(151) = 0.0300000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(152) = 0.900000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(153) = CONSTS(53) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(154) = CONSTS(57) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(155) = CONSTS(69) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(156) = CONSTS(87) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(157) = ( (( CONSTS(15)*2.00000*CONSTS(17))/10000.0)* 3.14159265358979*CONSTS(16))/10000.0
CONSTS(158) = ( 3.14159*CONSTS(17)/1000.00 ** 2.00000*CONSTS(16))/1000.00 ** 3.00000
CONSTS(159) = 0.0350000*CONSTS(158)
CONSTS(160) = 0.0200000*CONSTS(158)
CONSTS(161) = 0.000510000*CONSTS(158)
CONSTS(162) = 0.650000*CONSTS(158)
CONSTS(163) = CONSTS(158) - (CONSTS(159)+CONSTS(160)+CONSTS(161)+CONSTS(162))
CONSTS(164) = ( CONSTS(66)*CONSTS(162)*2.00000*CONSTS(14))/CONSTS(157)
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)) - 0.0119000*STATES(20)
RATES(21) = 1.70000*STATES(6)*(1.00000 - STATES(21)) - 0.0119000*STATES(21)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (50.0000+STATES(1))/7.50000))
ALGBRC(26) = 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(26)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(28) = 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(28)
ALGBRC(9) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(29) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
RATES(14) = (ALGBRC(9) - STATES(14))/ALGBRC(29)
ALGBRC(10) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(30) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
RATES(13) = (ALGBRC(10) - STATES(13))/ALGBRC(30)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(32) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
RATES(15) = (ALGBRC(11) - STATES(15))/ALGBRC(32)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(33) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
RATES(17) = (ALGBRC(12) - STATES(17))/ALGBRC(33)
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(34) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
RATES(16) = (ALGBRC(13) - STATES(16))/ALGBRC(34)
ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(35) = ( 1.00000*ALGBRC(14)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
RATES(18) = (ALGBRC(14) - STATES(18))/ALGBRC(35)
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(36) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
RATES(19) = (ALGBRC(15) - STATES(19))/ALGBRC(36)
ALGBRC(21) = STATES(40)+( STATES(42)*CONSTS(161))/CONSTS(162)+( STATES(41)*CONSTS(160))/CONSTS(162)
ALGBRC(38) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(5)*(CONSTS(136) - STATES(40)) - CONSTS(144)*STATES(40)
RATES(40) = ALGBRC(38)
ALGBRC(40) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(7)*(( CONSTS(137)*CONSTS(162))/CONSTS(160) - STATES(41)) - CONSTS(144)*STATES(41)
RATES(41) = ALGBRC(40)
ALGBRC(39) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(6)*(( CONSTS(138)*CONSTS(162))/CONSTS(161) - STATES(42)) - CONSTS(144)*STATES(42)
RATES(42) = ALGBRC(39)
ALGBRC(22) = STATES(43)+( STATES(45)*CONSTS(161))/CONSTS(162)+( STATES(44)*CONSTS(160))/CONSTS(162)
ALGBRC(41) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(5)*(CONSTS(139) - STATES(43)) - CONSTS(145)*STATES(43)
RATES(43) = ALGBRC(41)
ALGBRC(43) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(7)*(( CONSTS(140)*CONSTS(162))/CONSTS(160) - STATES(44)) - CONSTS(145)*STATES(44)
RATES(44) = ALGBRC(43)
ALGBRC(42) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(6)*(( CONSTS(141)*CONSTS(162))/CONSTS(161) - STATES(45)) - CONSTS(145)*STATES(45)
RATES(45) = ALGBRC(42)
ALGBRC(16) = CONSTS(76) - (CONSTS(76) - CONSTS(77))/(1.00000+CONSTS(78)/STATES(25) ** CONSTS(84))
ALGBRC(37) = CONSTS(80)/ALGBRC(16)
ALGBRC(48) = CONSTS(82)*ALGBRC(16)
RATES(24) = ( ALGBRC(37)*STATES(6) ** 2.00000*STATES(22) - CONSTS(81)*STATES(24)) - ( ALGBRC(48)*STATES(6)*STATES(24) - CONSTS(83)*STATES(23))
ALGBRC(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(23) = 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(44) = 1.00000/(ALGBRC(1)+ALGBRC(23))
ALGBRC(49) = ALGBRC(1)/(ALGBRC(1)+ALGBRC(23))
RATES(9) = (ALGBRC(49) - STATES(9))/ALGBRC(44)
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(24) = 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(45) = 1.00000/(ALGBRC(2)+ALGBRC(24))
ALGBRC(50) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(24))
RATES(10) = (ALGBRC(50) - STATES(10))/ALGBRC(45)
ALGBRC(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(25) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(46) = 1.00000/(ALGBRC(3)+ALGBRC(25))
ALGBRC(51) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(25))
RATES(8) = (ALGBRC(51) - STATES(8))/ALGBRC(46)
ALGBRC(53) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
RATES(22) = ( CONSTS(83)*ALGBRC(53) - ALGBRC(48)*STATES(6)*STATES(22)) - ( ALGBRC(37)*STATES(6) ** 2.00000*STATES(22) - CONSTS(81)*STATES(24))
RATES(23) = ( ALGBRC(48)*STATES(6)*STATES(24) - CONSTS(83)*STATES(23)) - ( CONSTS(81)*STATES(23) - ALGBRC(37)*STATES(6) ** 2.00000*ALGBRC(53))
ALGBRC(101) = CONSTS(101)*STATES(3)*(CONSTS(100) - STATES(27)) - CONSTS(102)*STATES(27)
RATES(27) = ALGBRC(101)
ALGBRC(27) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(3))
ALGBRC(55) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(56) = CONSTS(19)*CONSTS(18)*ALGBRC(55)*(STATES(1) - ALGBRC(27))
ALGBRC(59) = CONSTS(22)*CONSTS(21)*(STATES(1) - ALGBRC(27))
ALGBRC(62) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))+ 0.0365000*CONSTS(150)*EXP(( - STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(63) = ( (( CONSTS(30)*CONSTS(148)*CONSTS(29)*ALGBRC(62))/(1.00000+( CONSTS(149)*CONSTS(27))/STATES(3) ** CONSTS(25)))*CONSTS(5))/(CONSTS(5)+CONSTS(28))
ALGBRC(88) = 1.00000 - STATES(21)
ALGBRC(84) = ( 0.450000*STATES(18)*STATES(19)*CONSTS(153)*STATES(1)*CONSTS(14) ** 2.00000)/( CONSTS(12)*CONSTS(13))
ALGBRC(90) = ( ALGBRC(84)*ALGBRC(88)*CONSTS(45)*CONSTS(43)*( CONSTS(49)*STATES(3)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(50)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(96) = 1.00000/(1.00000+CONSTS(62)/STATES(6) ** 3.00000)
ALGBRC(93) = ( EXP(( CONSTS(64)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(3) ** CONSTS(65)*CONSTS(7) - EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(65)*STATES(6))/(1.00000+ CONSTS(63)*EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(97) = ( CONSTS(55)*CONSTS(54)*ALGBRC(96)*CONSTS(154)*ALGBRC(93))/( CONSTS(61)*CONSTS(6) ** CONSTS(65)*(1.00000+STATES(3)/CONSTS(58) ** CONSTS(65))+ CONSTS(60) ** CONSTS(65)*STATES(6)*(1.00000+STATES(6)/CONSTS(61))+ CONSTS(59)*STATES(3) ** CONSTS(65)+ STATES(3) ** CONSTS(65)*CONSTS(7)+ CONSTS(6) ** CONSTS(65)*STATES(6))
ALGBRC(99) = (STATES(3) - STATES(4))*1.83130e-14
RATES(3) = (( - CONSTS(157)*(ALGBRC(56)+ 3.00000*ALGBRC(97)+ALGBRC(59)+ 3.00000*ALGBRC(63)+ALGBRC(90)))/( CONSTS(161)*CONSTS(14)) - ALGBRC(99)/CONSTS(161)) - ALGBRC(101)
ALGBRC(102) = (STATES(4) - STATES(2))*1.63860e-12
RATES(2) = ALGBRC(102)/CONSTS(162)
ALGBRC(104) = CONSTS(101)*STATES(4)*(CONSTS(99) - STATES(26)) - CONSTS(102)*STATES(26)
RATES(26) = ALGBRC(104)
ALGBRC(47) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(4))
ALGBRC(57) = CONSTS(20)*CONSTS(18)*ALGBRC(55)*(STATES(1) - ALGBRC(47))
ALGBRC(60) = CONSTS(23)*CONSTS(21)*(STATES(1) - ALGBRC(47))
ALGBRC(64) = ( (( CONSTS(31)*CONSTS(148)*CONSTS(29)*ALGBRC(62))/(1.00000+( CONSTS(149)*CONSTS(27))/STATES(4) ** CONSTS(25)))*CONSTS(5))/(CONSTS(5)+CONSTS(28))
ALGBRC(85) = 1.00000 - STATES(20)
ALGBRC(87) = ( ALGBRC(84)*ALGBRC(85)*CONSTS(46)*CONSTS(43)*( CONSTS(49)*STATES(4)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(50)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(95) = 1.00000/(1.00000+CONSTS(62)/STATES(7) ** 3.00000)
ALGBRC(94) = ( EXP(( CONSTS(64)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(4) ** CONSTS(65)*CONSTS(7) - EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(65)*STATES(7))/(1.00000+ CONSTS(63)*EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(98) = ( CONSTS(56)*CONSTS(54)*ALGBRC(95)*CONSTS(154)*ALGBRC(94))/( CONSTS(61)*CONSTS(6) ** CONSTS(65)*(1.00000+STATES(4)/CONSTS(58) ** CONSTS(65))+ CONSTS(60) ** CONSTS(65)*STATES(7)*(1.00000+STATES(7)/CONSTS(61))+ CONSTS(59)*STATES(4) ** CONSTS(65)+ STATES(4) ** CONSTS(65)*CONSTS(7)+ CONSTS(6) ** CONSTS(65)*STATES(7))
RATES(4) = (( - CONSTS(157)*(ALGBRC(57)+ 3.00000*ALGBRC(98)+ALGBRC(60)+ 3.00000*ALGBRC(64)+ALGBRC(87)))/( CONSTS(160)*CONSTS(14))+(ALGBRC(99) - ALGBRC(102))/CONSTS(160)) - ALGBRC(104)
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(61) = ALGBRC(59)+ALGBRC(60)
ALGBRC(65) = ALGBRC(63)+ALGBRC(64)
ALGBRC(66) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(67) = CONSTS(151)*STATES(11)*ALGBRC(66)*(STATES(1) - CONSTS(146))
ALGBRC(68) = - arbitrary_log(STATES(6)/1.00000, 10)+3.00000
ALGBRC(70) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(68))/0.600000)))
ALGBRC(72) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(34)*CONSTS(6))/(CONSTS(9)+ CONSTS(34)*STATES(2)))
ALGBRC(73) = CONSTS(32)*ALGBRC(70)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(72))
ALGBRC(69) = - arbitrary_log(STATES(7)/1.00000, 10)+3.00000
ALGBRC(71) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(69))/0.600000)))
ALGBRC(74) = CONSTS(33)*ALGBRC(71)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(72))
ALGBRC(75) = ALGBRC(73)+ALGBRC(74)
ALGBRC(76) = CONSTS(35)*STATES(14)*(STATES(13)+ 0.500000*STATES(15))*(STATES(1) - CONSTS(146))
ALGBRC(77) = CONSTS(36)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(146))
ALGBRC(78) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(146)+59.2150))))
ALGBRC(79) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(146))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(146)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(146))+4.75300)))
ALGBRC(80) = ALGBRC(78)/(ALGBRC(78)+ALGBRC(79))
ALGBRC(81) = CONSTS(152)*ALGBRC(80)*(STATES(1) - CONSTS(146))
ALGBRC(100) = ALGBRC(97)+ALGBRC(98)
ALGBRC(82) = CONSTS(37)*(STATES(1) - CONSTS(147))*(CONSTS(39)/(1.00000+CONSTS(38)/STATES(6))+CONSTS(40)/(1.00000+CONSTS(38)/STATES(7)))
ALGBRC(83) = CONSTS(41)*(STATES(1) - CONSTS(147))
ALGBRC(86) = ( ALGBRC(84)*ALGBRC(85)*CONSTS(46)*CONSTS(42)*4.00000*( CONSTS(47)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(48)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(89) = ( ALGBRC(84)*ALGBRC(88)*CONSTS(45)*CONSTS(42)*4.00000*( CONSTS(47)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(48)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(91) = ( ALGBRC(84)*( ALGBRC(85)*CONSTS(46)+ ALGBRC(88)*CONSTS(45))*CONSTS(44)*( CONSTS(51)*CONSTS(9)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(52)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(92) = ALGBRC(86)+ALGBRC(89)+ALGBRC(87)+ALGBRC(90)+ALGBRC(91)
ALGBRC(52) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(6))
ALGBRC(107) = CONSTS(73)*CONSTS(74)*(STATES(1) - ALGBRC(52))
ALGBRC(54) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(7))
ALGBRC(108) = CONSTS(73)*CONSTS(75)*(STATES(1) - ALGBRC(54))
ALGBRC(109) = ALGBRC(108)+ALGBRC(107)
ALGBRC(105) = ( CONSTS(155)*CONSTS(72)*CONSTS(68))/(1.00000+CONSTS(70)/STATES(7) ** CONSTS(71))
ALGBRC(103) = ( CONSTS(155)*CONSTS(72)*CONSTS(67))/(1.00000+CONSTS(70)/STATES(6) ** CONSTS(71))
ALGBRC(106) = ALGBRC(103)+ALGBRC(105)
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(58)+ALGBRC(61)+ALGBRC(65)+ALGBRC(67)+ALGBRC(75)+ALGBRC(76)+ALGBRC(77)+ALGBRC(81)+ALGBRC(100)+ALGBRC(82)+ALGBRC(83)+ALGBRC(92)+ALGBRC(109)+ALGBRC(106)+ALGBRC(7))
ALGBRC(113) = CONSTS(109)*STATES(25)*(( CONSTS(107)*CONSTS(162))/CONSTS(159) - STATES(32)) - CONSTS(112)*STATES(32)
RATES(32) = ALGBRC(113)
ALGBRC(110) = CONSTS(79)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(111) = CONSTS(85)*(STATES(25) - STATES(6))
ALGBRC(112) = ( (( CONSTS(156)*CONSTS(86)*CONSTS(162))/CONSTS(159))*(STATES(5)/CONSTS(88) ** CONSTS(90) - STATES(25)/CONSTS(89) ** CONSTS(90)))/(1.00000+STATES(5)/CONSTS(88) ** CONSTS(90)+STATES(25)/CONSTS(89) ** CONSTS(90))
RATES(25) = (ALGBRC(112) - (( ALGBRC(111)*CONSTS(162))/CONSTS(159)+ALGBRC(110))) - ALGBRC(113)
ALGBRC(116) = CONSTS(108)*STATES(6)*(( CONSTS(104)*0.100000*CONSTS(162))/CONSTS(161) - STATES(29)) - CONSTS(110)*STATES(29)
RATES(29) = ALGBRC(116)
ALGBRC(118) = CONSTS(108)*STATES(7)*(( CONSTS(103)*CONSTS(162))/CONSTS(160) - STATES(28)) - CONSTS(110)*STATES(28)
RATES(28) = ALGBRC(118)
ALGBRC(119) = CONSTS(108)*STATES(6)*(( CONSTS(106)*0.100000*CONSTS(162))/CONSTS(161) - STATES(31)) - CONSTS(111)*STATES(31)
RATES(31) = ALGBRC(119)
ALGBRC(117) = CONSTS(120)*STATES(5)*(CONSTS(113) - STATES(33)) - CONSTS(127)*STATES(33)
RATES(33) = ALGBRC(117)
ALGBRC(121) = CONSTS(108)*STATES(7)*(( CONSTS(105)*CONSTS(162))/CONSTS(160) - STATES(30)) - CONSTS(111)*STATES(30)
RATES(30) = ALGBRC(121)
ALGBRC(120) = CONSTS(121)*STATES(5)*(CONSTS(114) - (STATES(34)+STATES(35))) - CONSTS(128)*STATES(34)
RATES(34) = ALGBRC(120)
ALGBRC(114) = (STATES(6) - STATES(7))*8.24130e-13
ALGBRC(122) = ALGBRC(116)+ALGBRC(119)
ALGBRC(124) = (ALGBRC(89) - 2.00000*ALGBRC(97))+ALGBRC(107)+ALGBRC(103)
RATES(6) = ((( - ALGBRC(124)*CONSTS(157))/( CONSTS(161)*2.00000*CONSTS(14)) - ALGBRC(114)/CONSTS(161))+( ALGBRC(110)*CONSTS(159))/CONSTS(161)+( ALGBRC(111)*CONSTS(162))/CONSTS(161)) - 1.00000*ALGBRC(122)
ALGBRC(123) = CONSTS(122)*CONSTS(10)*(CONSTS(115) - (STATES(34)+STATES(35))) - CONSTS(129)*STATES(35)
RATES(35) = ALGBRC(123)
ALGBRC(115) = (STATES(7) - STATES(5))*3.72430e-12
ALGBRC(125) = ALGBRC(118)+ALGBRC(121)
ALGBRC(127) = (ALGBRC(86) - 2.00000*ALGBRC(98))+ALGBRC(108)+ALGBRC(105)
RATES(7) = (( - ALGBRC(127)*CONSTS(157))/( CONSTS(160)*2.00000*CONSTS(14))+(ALGBRC(114) - ALGBRC(115))/CONSTS(160)) - 1.00000*ALGBRC(125)
ALGBRC(126) = CONSTS(123)*STATES(5)*(CONSTS(116) - STATES(36)) - CONSTS(130)*STATES(36)
RATES(36) = ALGBRC(126)
ALGBRC(128) = CONSTS(124)*STATES(5)*(CONSTS(117) - (STATES(37)+STATES(38))) - CONSTS(131)*STATES(37)
RATES(37) = ALGBRC(128)
ALGBRC(129) = CONSTS(125)*CONSTS(10)*(CONSTS(118) - (STATES(37)+STATES(38))) - CONSTS(132)*STATES(38)
RATES(38) = ALGBRC(129)
ALGBRC(130) = CONSTS(126)*STATES(5)*(CONSTS(119) - STATES(39)) - CONSTS(133)*STATES(39)
RATES(39) = ALGBRC(130)
ALGBRC(131) = ALGBRC(117)+ALGBRC(120)+ALGBRC(123)+ALGBRC(126)+ALGBRC(128)+ALGBRC(129)+ALGBRC(130)
RATES(5) = (( - ALGBRC(112)*CONSTS(159))/CONSTS(162)+ALGBRC(115)/CONSTS(162)) - 1.00000*ALGBRC(131)
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(26) = 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(28) = 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(29) = 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(30) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(32) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(33) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(34) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
ALGBRC(14) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(35) = ( 1.00000*ALGBRC(14)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
ALGBRC(15) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(36) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
ALGBRC(21) = STATES(40)+( STATES(42)*CONSTS(161))/CONSTS(162)+( STATES(41)*CONSTS(160))/CONSTS(162)
ALGBRC(38) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(5)*(CONSTS(136) - STATES(40)) - CONSTS(144)*STATES(40)
ALGBRC(40) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(7)*(( CONSTS(137)*CONSTS(162))/CONSTS(160) - STATES(41)) - CONSTS(144)*STATES(41)
ALGBRC(39) = CONSTS(142)*(CONSTS(134) - ALGBRC(21))*STATES(6)*(( CONSTS(138)*CONSTS(162))/CONSTS(161) - STATES(42)) - CONSTS(144)*STATES(42)
ALGBRC(22) = STATES(43)+( STATES(45)*CONSTS(161))/CONSTS(162)+( STATES(44)*CONSTS(160))/CONSTS(162)
ALGBRC(41) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(5)*(CONSTS(139) - STATES(43)) - CONSTS(145)*STATES(43)
ALGBRC(43) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(7)*(( CONSTS(140)*CONSTS(162))/CONSTS(160) - STATES(44)) - CONSTS(145)*STATES(44)
ALGBRC(42) = CONSTS(143)*(CONSTS(135) - ALGBRC(22))*STATES(6)*(( CONSTS(141)*CONSTS(162))/CONSTS(161) - STATES(45)) - CONSTS(145)*STATES(45)
ALGBRC(16) = CONSTS(76) - (CONSTS(76) - CONSTS(77))/(1.00000+CONSTS(78)/STATES(25) ** CONSTS(84))
ALGBRC(37) = CONSTS(80)/ALGBRC(16)
ALGBRC(48) = CONSTS(82)*ALGBRC(16)
ALGBRC(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(23) = 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(44) = 1.00000/(ALGBRC(1)+ALGBRC(23))
ALGBRC(49) = ALGBRC(1)/(ALGBRC(1)+ALGBRC(23))
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(24) = 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(45) = 1.00000/(ALGBRC(2)+ALGBRC(24))
ALGBRC(50) = ALGBRC(2)/(ALGBRC(2)+ALGBRC(24))
ALGBRC(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(25) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(46) = 1.00000/(ALGBRC(3)+ALGBRC(25))
ALGBRC(51) = ALGBRC(3)/(ALGBRC(3)+ALGBRC(25))
ALGBRC(53) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
ALGBRC(101) = CONSTS(101)*STATES(3)*(CONSTS(100) - STATES(27)) - CONSTS(102)*STATES(27)
ALGBRC(27) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(3))
ALGBRC(55) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(56) = CONSTS(19)*CONSTS(18)*ALGBRC(55)*(STATES(1) - ALGBRC(27))
ALGBRC(59) = CONSTS(22)*CONSTS(21)*(STATES(1) - ALGBRC(27))
ALGBRC(62) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))+ 0.0365000*CONSTS(150)*EXP(( - STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(63) = ( (( CONSTS(30)*CONSTS(148)*CONSTS(29)*ALGBRC(62))/(1.00000+( CONSTS(149)*CONSTS(27))/STATES(3) ** CONSTS(25)))*CONSTS(5))/(CONSTS(5)+CONSTS(28))
ALGBRC(88) = 1.00000 - STATES(21)
ALGBRC(84) = ( 0.450000*STATES(18)*STATES(19)*CONSTS(153)*STATES(1)*CONSTS(14) ** 2.00000)/( CONSTS(12)*CONSTS(13))
ALGBRC(90) = ( ALGBRC(84)*ALGBRC(88)*CONSTS(45)*CONSTS(43)*( CONSTS(49)*STATES(3)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(50)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(96) = 1.00000/(1.00000+CONSTS(62)/STATES(6) ** 3.00000)
ALGBRC(93) = ( EXP(( CONSTS(64)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(3) ** CONSTS(65)*CONSTS(7) - EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(65)*STATES(6))/(1.00000+ CONSTS(63)*EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(97) = ( CONSTS(55)*CONSTS(54)*ALGBRC(96)*CONSTS(154)*ALGBRC(93))/( CONSTS(61)*CONSTS(6) ** CONSTS(65)*(1.00000+STATES(3)/CONSTS(58) ** CONSTS(65))+ CONSTS(60) ** CONSTS(65)*STATES(6)*(1.00000+STATES(6)/CONSTS(61))+ CONSTS(59)*STATES(3) ** CONSTS(65)+ STATES(3) ** CONSTS(65)*CONSTS(7)+ CONSTS(6) ** CONSTS(65)*STATES(6))
ALGBRC(99) = (STATES(3) - STATES(4))*1.83130e-14
ALGBRC(102) = (STATES(4) - STATES(2))*1.63860e-12
ALGBRC(104) = CONSTS(101)*STATES(4)*(CONSTS(99) - STATES(26)) - CONSTS(102)*STATES(26)
ALGBRC(47) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(4))
ALGBRC(57) = CONSTS(20)*CONSTS(18)*ALGBRC(55)*(STATES(1) - ALGBRC(47))
ALGBRC(60) = CONSTS(23)*CONSTS(21)*(STATES(1) - ALGBRC(47))
ALGBRC(64) = ( (( CONSTS(31)*CONSTS(148)*CONSTS(29)*ALGBRC(62))/(1.00000+( CONSTS(149)*CONSTS(27))/STATES(4) ** CONSTS(25)))*CONSTS(5))/(CONSTS(5)+CONSTS(28))
ALGBRC(85) = 1.00000 - STATES(20)
ALGBRC(87) = ( ALGBRC(84)*ALGBRC(85)*CONSTS(46)*CONSTS(43)*( CONSTS(49)*STATES(4)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(50)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(95) = 1.00000/(1.00000+CONSTS(62)/STATES(7) ** 3.00000)
ALGBRC(94) = ( EXP(( CONSTS(64)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(4) ** CONSTS(65)*CONSTS(7) - EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(65)*STATES(7))/(1.00000+ CONSTS(63)*EXP(( (CONSTS(64) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(98) = ( CONSTS(56)*CONSTS(54)*ALGBRC(95)*CONSTS(154)*ALGBRC(94))/( CONSTS(61)*CONSTS(6) ** CONSTS(65)*(1.00000+STATES(4)/CONSTS(58) ** CONSTS(65))+ CONSTS(60) ** CONSTS(65)*STATES(7)*(1.00000+STATES(7)/CONSTS(61))+ CONSTS(59)*STATES(4) ** CONSTS(65)+ STATES(4) ** CONSTS(65)*CONSTS(7)+ CONSTS(6) ** CONSTS(65)*STATES(7))
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(61) = ALGBRC(59)+ALGBRC(60)
ALGBRC(65) = ALGBRC(63)+ALGBRC(64)
ALGBRC(66) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(67) = CONSTS(151)*STATES(11)*ALGBRC(66)*(STATES(1) - CONSTS(146))
ALGBRC(68) = - arbitrary_log(STATES(6)/1.00000, 10)+3.00000
ALGBRC(70) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(68))/0.600000)))
ALGBRC(72) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(34)*CONSTS(6))/(CONSTS(9)+ CONSTS(34)*STATES(2)))
ALGBRC(73) = CONSTS(32)*ALGBRC(70)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(72))
ALGBRC(69) = - arbitrary_log(STATES(7)/1.00000, 10)+3.00000
ALGBRC(71) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(69))/0.600000)))
ALGBRC(74) = CONSTS(33)*ALGBRC(71)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(72))
ALGBRC(75) = ALGBRC(73)+ALGBRC(74)
ALGBRC(76) = CONSTS(35)*STATES(14)*(STATES(13)+ 0.500000*STATES(15))*(STATES(1) - CONSTS(146))
ALGBRC(77) = CONSTS(36)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(146))
ALGBRC(78) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(146)+59.2150))))
ALGBRC(79) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(146))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(146)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(146))+4.75300)))
ALGBRC(80) = ALGBRC(78)/(ALGBRC(78)+ALGBRC(79))
ALGBRC(81) = CONSTS(152)*ALGBRC(80)*(STATES(1) - CONSTS(146))
ALGBRC(100) = ALGBRC(97)+ALGBRC(98)
ALGBRC(82) = CONSTS(37)*(STATES(1) - CONSTS(147))*(CONSTS(39)/(1.00000+CONSTS(38)/STATES(6))+CONSTS(40)/(1.00000+CONSTS(38)/STATES(7)))
ALGBRC(83) = CONSTS(41)*(STATES(1) - CONSTS(147))
ALGBRC(86) = ( ALGBRC(84)*ALGBRC(85)*CONSTS(46)*CONSTS(42)*4.00000*( CONSTS(47)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(48)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(89) = ( ALGBRC(84)*ALGBRC(88)*CONSTS(45)*CONSTS(42)*4.00000*( CONSTS(47)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(48)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(91) = ( ALGBRC(84)*( ALGBRC(85)*CONSTS(46)+ ALGBRC(88)*CONSTS(45))*CONSTS(44)*( CONSTS(51)*CONSTS(9)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(52)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(92) = ALGBRC(86)+ALGBRC(89)+ALGBRC(87)+ALGBRC(90)+ALGBRC(91)
ALGBRC(52) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(6))
ALGBRC(107) = CONSTS(73)*CONSTS(74)*(STATES(1) - ALGBRC(52))
ALGBRC(54) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(7))
ALGBRC(108) = CONSTS(73)*CONSTS(75)*(STATES(1) - ALGBRC(54))
ALGBRC(109) = ALGBRC(108)+ALGBRC(107)
ALGBRC(105) = ( CONSTS(155)*CONSTS(72)*CONSTS(68))/(1.00000+CONSTS(70)/STATES(7) ** CONSTS(71))
ALGBRC(103) = ( CONSTS(155)*CONSTS(72)*CONSTS(67))/(1.00000+CONSTS(70)/STATES(6) ** CONSTS(71))
ALGBRC(106) = ALGBRC(103)+ALGBRC(105)
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(113) = CONSTS(109)*STATES(25)*(( CONSTS(107)*CONSTS(162))/CONSTS(159) - STATES(32)) - CONSTS(112)*STATES(32)
ALGBRC(110) = CONSTS(79)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(111) = CONSTS(85)*(STATES(25) - STATES(6))
ALGBRC(112) = ( (( CONSTS(156)*CONSTS(86)*CONSTS(162))/CONSTS(159))*(STATES(5)/CONSTS(88) ** CONSTS(90) - STATES(25)/CONSTS(89) ** CONSTS(90)))/(1.00000+STATES(5)/CONSTS(88) ** CONSTS(90)+STATES(25)/CONSTS(89) ** CONSTS(90))
ALGBRC(116) = CONSTS(108)*STATES(6)*(( CONSTS(104)*0.100000*CONSTS(162))/CONSTS(161) - STATES(29)) - CONSTS(110)*STATES(29)
ALGBRC(118) = CONSTS(108)*STATES(7)*(( CONSTS(103)*CONSTS(162))/CONSTS(160) - STATES(28)) - CONSTS(110)*STATES(28)
ALGBRC(119) = CONSTS(108)*STATES(6)*(( CONSTS(106)*0.100000*CONSTS(162))/CONSTS(161) - STATES(31)) - CONSTS(111)*STATES(31)
ALGBRC(117) = CONSTS(120)*STATES(5)*(CONSTS(113) - STATES(33)) - CONSTS(127)*STATES(33)
ALGBRC(121) = CONSTS(108)*STATES(7)*(( CONSTS(105)*CONSTS(162))/CONSTS(160) - STATES(30)) - CONSTS(111)*STATES(30)
ALGBRC(120) = CONSTS(121)*STATES(5)*(CONSTS(114) - (STATES(34)+STATES(35))) - CONSTS(128)*STATES(34)
ALGBRC(114) = (STATES(6) - STATES(7))*8.24130e-13
ALGBRC(122) = ALGBRC(116)+ALGBRC(119)
ALGBRC(124) = (ALGBRC(89) - 2.00000*ALGBRC(97))+ALGBRC(107)+ALGBRC(103)
ALGBRC(123) = CONSTS(122)*CONSTS(10)*(CONSTS(115) - (STATES(34)+STATES(35))) - CONSTS(129)*STATES(35)
ALGBRC(115) = (STATES(7) - STATES(5))*3.72430e-12
ALGBRC(125) = ALGBRC(118)+ALGBRC(121)
ALGBRC(127) = (ALGBRC(86) - 2.00000*ALGBRC(98))+ALGBRC(108)+ALGBRC(105)
ALGBRC(126) = CONSTS(123)*STATES(5)*(CONSTS(116) - STATES(36)) - CONSTS(130)*STATES(36)
ALGBRC(128) = CONSTS(124)*STATES(5)*(CONSTS(117) - (STATES(37)+STATES(38))) - CONSTS(131)*STATES(37)
ALGBRC(129) = CONSTS(125)*CONSTS(10)*(CONSTS(118) - (STATES(37)+STATES(38))) - CONSTS(132)*STATES(38)
ALGBRC(130) = CONSTS(126)*STATES(5)*(CONSTS(119) - STATES(39)) - CONSTS(133)*STATES(39)
ALGBRC(131) = ALGBRC(117)+ALGBRC(120)+ALGBRC(123)+ALGBRC(126)+ALGBRC(128)+ALGBRC(129)+ALGBRC(130)
ALGBRC(5) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(34)*CONSTS(6))/(CONSTS(9)+ CONSTS(34)*STATES(3)))
ALGBRC(6) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(34)*CONSTS(6))/(CONSTS(9)+ CONSTS(34)*STATES(4)))
ALGBRC(17) = ( CONSTS(93)*CONSTS(97)*(STATES(3) - STATES(4)))/CONSTS(95)
ALGBRC(18) = ( CONSTS(94)*CONSTS(98)*(STATES(4) - STATES(2)))/CONSTS(96)
ALGBRC(19) = ( CONSTS(91)*CONSTS(97)*(STATES(6) - STATES(7)))/CONSTS(95)
ALGBRC(20) = ( CONSTS(92)*CONSTS(98)*(STATES(7) - STATES(5)))/CONSTS(96)
ALGBRC(31) = ALGBRC(11)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END