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