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 111 entries in the algebraic variable array.
C There are a total of 39 entries in each of the rate and state variable arrays.
C There are a total of 137 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(37) is i_Na in component INa (microA_per_microF).
C ALGBRC(40) is i_Nab in component INab (microA_per_microF).
C ALGBRC(44) is i_NaK in component INaK (microA_per_microF).
C ALGBRC(46) is i_Kr in component IKr (microA_per_microF).
C ALGBRC(54) is i_Ks in component IKs (microA_per_microF).
C ALGBRC(55) is i_Kp in component IKp (microA_per_microF).
C ALGBRC(56) is i_tos in component Itos (microA_per_microF).
C ALGBRC(57) is i_tof in component Itof (microA_per_microF).
C ALGBRC(61) is i_K1 in component IK1 (microA_per_microF).
C ALGBRC(80) is i_NaCa in component INaCa (microA_per_microF).
C ALGBRC(62) is i_Cl_Ca in component ICl_Ca (microA_per_microF).
C ALGBRC(63) is i_Clb in component IClb (microA_per_microF).
C ALGBRC(72) is i_CaL in component ICaL (microA_per_microF).
C ALGBRC(89) is i_Cab in component ICab (microA_per_microF).
C ALGBRC(86) is i_Cap in component ICap (microA_per_microF).
C ALGBRC(5) 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(15) is Cm in component model_parameters (farad).
C CONSTS(16) is cell_length in component model_parameters (micrometre).
C CONSTS(17) is cell_radius in component model_parameters (micrometre).
C CONSTS(133) is Vol_Cell in component model_parameters (litre).
C CONSTS(134) is Vol_SR in component model_parameters (litre).
C CONSTS(135) is Vol_SL in component model_parameters (litre).
C CONSTS(136) is Vol_jct in component model_parameters (litre).
C CONSTS(137) is Vol_myo in component model_parameters (litre).
C ALGBRC(19) is E_Na_jct in component reversal_potentials (millivolt).
C ALGBRC(29) is E_Na_SL in component reversal_potentials (millivolt).
C ALGBRC(31) is E_Ca_jct in component reversal_potentials (millivolt).
C ALGBRC(33) is E_Ca_SL in component reversal_potentials (millivolt).
C CONSTS(124) is E_K in component reversal_potentials (millivolt).
C CONSTS(125) 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(35) is i_Na_jct in component INa (microA_per_microF).
C ALGBRC(36) 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(34) is openProb in component INa (dimensionless).
C ALGBRC(1) is alpha_h in component INa_h_gate (per_millisecond).
C ALGBRC(15) is beta_h in component INa_h_gate (per_millisecond).
C ALGBRC(2) is alpha_j in component INa_j_gate (per_millisecond).
C ALGBRC(16) is beta_j in component INa_j_gate (per_millisecond).
C ALGBRC(3) is alpha_m in component INa_m_gate (per_millisecond).
C ALGBRC(17) is beta_m in component INa_m_gate (per_millisecond).
C ALGBRC(38) is i_Nab_jct in component INab (microA_per_microF).
C ALGBRC(39) 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(42) is i_NaK_jct in component INaK (microA_per_microF).
C ALGBRC(43) is i_NaK_SL in component INaK (microA_per_microF).
C ALGBRC(41) is f_NaK in component INaK (dimensionless).
C CONSTS(24) is H_NaK in component INaK (dimensionless).
C CONSTS(25) is Km_Nai in component INaK (millimolar).
C CONSTS(26) is Km_Ko in component INaK (millimolar).
C CONSTS(27) is I_NaK_max in component INaK (microA_per_microF).
C CONSTS(126) is sigma in component INaK (dimensionless).
C CONSTS(28) is Fx_NaK_jct in component INaK (dimensionless).
C CONSTS(29) is Fx_NaK_SL in component INaK (dimensionless).
C CONSTS(127) is G_IKr in component IKr (milliS_per_microF).
C STATES(11) is Xr in component IKr_Xr_gate (dimensionless).
C ALGBRC(45) is Rr in component IKr_Rr_gate (dimensionless).
C ALGBRC(4) is Xr_infinity in component IKr_Xr_gate (dimensionless).
C ALGBRC(18) is tau_Xr in component IKr_Xr_gate (millisecond).
C ALGBRC(52) is i_Ks_jct in component IKs (microA_per_microF).
C ALGBRC(53) is i_Ks_SL in component IKs (microA_per_microF).
C ALGBRC(50) is G_Ks_SL in component IKs (milliS_per_microF).
C ALGBRC(49) is G_Ks_jct in component IKs (milliS_per_microF).
C CONSTS(30) is Fx_Ks_jct in component IKs (dimensionless).
C CONSTS(31) is Fx_Ks_SL in component IKs (dimensionless).
C ALGBRC(51) is E_Ks in component IKs (millivolt).
C CONSTS(32) is pKNa in component IKs (dimensionless).
C ALGBRC(47) is pCa_jct in component IKs (dimensionless).
C ALGBRC(48) is pCa_SL in component IKs (dimensionless).
C STATES(12) is Xs in component IKs_Xs_gate (dimensionless).
C ALGBRC(6) is Xs_infinity in component IKs_Xs_gate (dimensionless).
C ALGBRC(20) is tau_Xs in component IKs_Xs_gate (millisecond).
C CONSTS(33) is g_Kp in component IKp (milliS_per_microF).
C CONSTS(34) 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(7) is X_tos_infinity in component Itos_X_gate (dimensionless).
C ALGBRC(21) is tau_X_tos in component Itos_X_gate (millisecond).
C ALGBRC(8) is Y_tos_infinity in component Itos_Y_gate (dimensionless).
C ALGBRC(22) is tau_Y_tos in component Itos_Y_gate (millisecond).
C ALGBRC(9) is R_tos_infinity in component Itos_R_gate (dimensionless).
C ALGBRC(23) is tau_R_tos in component Itos_R_gate (millisecond).
C CONSTS(35) 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(10) is X_tof_infinity in component Itof_X_gate (dimensionless).
C ALGBRC(24) is tau_X_tof in component Itof_X_gate (millisecond).
C ALGBRC(11) is Y_tof_infinity in component Itof_Y_gate (dimensionless).
C ALGBRC(25) is tau_Y_tof in component Itof_Y_gate (millisecond).
C CONSTS(128) is G_K1 in component IK1 (milliS_per_microF).
C ALGBRC(60) is K1_infinity in component IK1_K1_gate (dimensionless).
C ALGBRC(58) is alpha_K1 in component IK1_K1_gate (per_millisecond).
C ALGBRC(59) is beta_K1 in component IK1_K1_gate (per_millisecond).
C CONSTS(36) is G_Cl in component ICl_Ca (milliS_per_microF).
C CONSTS(37) is Kd_ClCa in component ICl_Ca (millimolar).
C CONSTS(38) is Fx_Cl_jct in component ICl_Ca (dimensionless).
C CONSTS(39) is Fx_Cl_SL in component ICl_Ca (dimensionless).
C CONSTS(40) is G_ClBk in component IClb (milliS_per_microF).
C ALGBRC(66) is i_CaL_Ca_SL in component ICaL (microA_per_microF).
C ALGBRC(69) is i_CaL_Ca_jct in component ICaL (microA_per_microF).
C ALGBRC(67) is i_CaL_Na_SL in component ICaL (microA_per_microF).
C ALGBRC(70) is i_CaL_Na_jct in component ICaL (microA_per_microF).
C ALGBRC(71) is i_CaL_K in component ICaL (microA_per_microF).
C CONSTS(41) is PCa in component ICaL (litre_per_farad_millisecond).
C CONSTS(42) is PNa in component ICaL (litre_per_farad_millisecond).
C CONSTS(43) is PK in component ICaL (litre_per_farad_millisecond).
C CONSTS(44) is Fx_ICaL_jct in component ICaL (dimensionless).
C CONSTS(45) is Fx_ICaL_SL in component ICaL (dimensionless).
C CONSTS(46) is gamma_Cai in component ICaL (dimensionless).
C CONSTS(47) is gamma_Cao in component ICaL (dimensionless).
C CONSTS(48) is gamma_Nai in component ICaL (dimensionless).
C CONSTS(49) is gamma_Nao in component ICaL (dimensionless).
C CONSTS(50) is gamma_Ki in component ICaL (dimensionless).
C CONSTS(51) is gamma_Ko in component ICaL (dimensionless).
C CONSTS(52) is Q10_CaL in component ICaL (dimensionless).
C CONSTS(129) 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(65) is fCa_SL in component ICaL_fCa_gate (dimensionless).
C ALGBRC(68) is fCa_jct in component ICaL_fCa_gate (dimensionless).
C ALGBRC(64) is temp in component ICaL (coulomb_per_mole).
C ALGBRC(12) is d_infinity in component ICaL_d_gate (dimensionless).
C ALGBRC(26) is tau_d in component ICaL_d_gate (millisecond).
C ALGBRC(13) is f_infinity in component ICaL_f_gate (dimensionless).
C ALGBRC(27) 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(78) is i_NaCa_SL in component INaCa (microA_per_microF).
C ALGBRC(77) is i_NaCa_jct in component INaCa (microA_per_microF).
C CONSTS(53) is V_max in component INaCa (microA_per_microF).
C CONSTS(54) is Fx_NCX_jct in component INaCa (dimensionless).
C CONSTS(55) is Fx_NCX_SL in component INaCa (dimensionless).
C CONSTS(56) is Q10_NCX in component INaCa (dimensionless).
C CONSTS(130) is Q_NCX in component INaCa (dimensionless).
C CONSTS(57) is K_mNai in component INaCa (millimolar).
C CONSTS(58) is K_mCao in component INaCa (millimolar).
C CONSTS(59) is K_mNao in component INaCa (millimolar).
C CONSTS(60) is K_mCai in component INaCa (millimolar).
C ALGBRC(75) is Ka_SL in component INaCa (dimensionless).
C ALGBRC(76) is Ka_jct in component INaCa (dimensionless).
C CONSTS(61) is Kd_act in component INaCa (millimolar).
C CONSTS(62) is ksat in component INaCa (dimensionless).
C CONSTS(63) is eta in component INaCa (dimensionless).
C CONSTS(64) is HNa in component INaCa (dimensionless).
C ALGBRC(73) is temp_jct in component INaCa (millimolar4).
C ALGBRC(74) is temp_SL in component INaCa (millimolar4).
C ALGBRC(85) is i_Cap_SL in component ICap (microA_per_microF).
C ALGBRC(83) is i_Cap_jct in component ICap (microA_per_microF).
C CONSTS(65) is Fx_SLCaP_jct in component ICap (dimensionless).
C CONSTS(66) is Fx_SLCaP_SL in component ICap (dimensionless).
C CONSTS(67) is Q10_SLCaP in component ICap (dimensionless).
C CONSTS(131) is Q_SLCaP in component ICap (dimensionless).
C CONSTS(68) is Km in component ICap (millimolar).
C CONSTS(69) is H in component ICap (dimensionless).
C CONSTS(70) is V_maxAF in component ICap (microA_per_microF).
C ALGBRC(87) is i_Cab_jct in component ICab (microA_per_microF).
C ALGBRC(88) is i_Cab_SL in component ICab (microA_per_microF).
C CONSTS(71) is G_CaBk in component ICab (milliS_per_microF).
C CONSTS(72) is Fx_CaBk_jct in component ICab (dimensionless).
C CONSTS(73) is Fx_CaBk_SL in component ICab (dimensionless).
C ALGBRC(90) is j_rel_SR in component Jrel_SR (millimolar_per_millisecond).
C CONSTS(74) is Max_SR in component Jrel_SR (dimensionless).
C CONSTS(75) is Min_SR in component Jrel_SR (dimensionless).
C CONSTS(76) 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(32) is RI in component Jrel_SR (dimensionless).
C STATES(24) is O in component Jrel_SR (dimensionless).
C CONSTS(77) is ks in component Jrel_SR (per_millisecond).
C CONSTS(78) is koCa in component Jrel_SR (per_millimolar2_per_millisecond).
C CONSTS(79) is kom in component Jrel_SR (per_millisecond).
C CONSTS(80) is kiCa in component Jrel_SR (per_millimolar_per_millisecond).
C CONSTS(81) is kim in component Jrel_SR (per_millisecond).
C CONSTS(82) is HSR in component Jrel_SR (dimensionless).
C STATES(25) is Ca_SR in component Ca_buffer (millimolar).
C ALGBRC(14) is kCaSR in component Jrel_SR (dimensionless).
C ALGBRC(28) is koSRCa in component Jrel_SR (per_millimolar2_per_millisecond).
C ALGBRC(30) is kiSRCa in component Jrel_SR (per_millimolar_per_millisecond).
C ALGBRC(91) is j_leak_SR in component Jleak_SR (millimolar_per_millisecond).
C CONSTS(83) is KSRleak in component Jleak_SR (per_millisecond).
C ALGBRC(92) is j_pump_SR in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(84) is V_max in component Jpump_SR (millimolar_per_millisecond).
C CONSTS(85) is Q10_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(132) is Q_SRCaP in component Jpump_SR (dimensionless).
C CONSTS(86) is Kmf in component Jpump_SR (millimolar).
C CONSTS(87) is Kmr in component Jpump_SR (millimolar).
C CONSTS(88) is H in component Jpump_SR (dimensionless).
C ALGBRC(79) is J_Na_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(82) is J_Na_SL_myo in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(94) is J_Ca_jct_SL in component ion_diffusion (millimole_per_millisecond).
C ALGBRC(95) is J_Ca_SL_myo 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(89) is Bmax_SL in component Na_buffer (millimolar).
C CONSTS(90) is Bmax_jct in component Na_buffer (millimolar).
C CONSTS(91) is kon in component Na_buffer (per_millimolar_per_millisecond).
C CONSTS(92) is koff in component Na_buffer (per_millisecond).
C ALGBRC(81) is dNa_jct_buf in component Na_buffer (millimolar_per_millisecond).
C ALGBRC(84) 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(93) is Bmax_SLB_SL in component Ca_buffer (millimolar).
C CONSTS(94) is Bmax_SLB_jct in component Ca_buffer (millimolar).
C CONSTS(95) is Bmax_SLHigh_SL in component Ca_buffer (millimolar).
C CONSTS(96) is Bmax_SLHigh_jct in component Ca_buffer (millimolar).
C CONSTS(97) is Bmax_Calsequestrin in component Ca_buffer (millimolar).
C CONSTS(98) is kon_SL in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(99) is kon_Calsequestrin in component Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(100) is koff_SLB in component Ca_buffer (per_millisecond).
C CONSTS(101) is koff_SLHigh in component Ca_buffer (per_millisecond).
C CONSTS(102) is koff_Calsequestrin in component Ca_buffer (per_millisecond).
C ALGBRC(93) is dCalsequestrin in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(111) is dCa_cytosol_tot_bound in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(98) is dCa_SLB_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(96) is dCa_SLB_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(101) is dCa_SLHigh_SL in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(99) is dCa_SLHigh_jct in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(102) is dCa_jct_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(105) is dCa_SL_tot_bound in component Ca_buffer (millimolar_per_millisecond).
C ALGBRC(104) is i_Ca_jct_tot in component Ca_buffer (microA_per_microF).
C ALGBRC(107) 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(103) is Bmax_TroponinC in component cytosolic_Ca_buffer (millimolar).
C CONSTS(104) is Bmax_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(105) is Bmax_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(106) is Bmax_Calmodulin in component cytosolic_Ca_buffer (millimolar).
C CONSTS(107) is Bmax_Myosin_Ca in component cytosolic_Ca_buffer (millimolar).
C CONSTS(108) is Bmax_Myosin_Mg in component cytosolic_Ca_buffer (millimolar).
C CONSTS(109) is Bmax_SRB in component cytosolic_Ca_buffer (millimolar).
C CONSTS(110) is kon_TroponinC in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(111) is kon_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(112) is kon_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(113) is kon_Calmodulin in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(114) is kon_Myosin_Ca in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(115) is kon_Myosin_Mg in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(116) is kon_SRB in component cytosolic_Ca_buffer (per_millimolar_per_millisecond).
C CONSTS(117) is koff_TroponinC in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(118) is koff_TroponinC_Ca_Mg_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(119) is koff_TroponinC_Ca_Mg_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(120) is koff_Calmodulin in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(121) is koff_Myosin_Ca in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(122) is koff_Myosin_Mg in component cytosolic_Ca_buffer (per_millisecond).
C CONSTS(123) is koff_SRB in component cytosolic_Ca_buffer (per_millisecond).
C ALGBRC(97) is dCa_TroponinC in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(100) is dCa_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(103) is dMg_TroponinC_Ca_Mg in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(106) is dCa_Calmodulin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(108) is dCa_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(109) is dMg_Myosin in component cytosolic_Ca_buffer (millimolar_per_millisecond).
C ALGBRC(110) is dCa_SRB in component cytosolic_Ca_buffer (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
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -8.556885e1
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) = 96485
CONSTS(15) = 1.381e-10
CONSTS(16) = 100
CONSTS(17) = 10.25
STATES(2) = 8.80853
STATES(3) = 8.80329
STATES(4) = 8.80733
STATES(5) = 8.597401e-5
STATES(6) = 1.737475e-4
STATES(7) = 1.031812e-4
CONSTS(18) = 16
CONSTS(19) = 0.11
CONSTS(20) = 0.89
STATES(8) = 1.405627e-3
STATES(9) = 9.867005e-1
STATES(10) = 9.91562e-1
CONSTS(21) = 0.297e-3
CONSTS(22) = 0.11
CONSTS(23) = 0.89
CONSTS(24) = 4
CONSTS(25) = 11
CONSTS(26) = 1.5
CONSTS(27) = 1.90719
CONSTS(28) = 0.11
CONSTS(29) = 0.89
STATES(11) = 8.641386e-3
CONSTS(30) = 0.11
CONSTS(31) = 0.89
CONSTS(32) = 0.01833
STATES(12) = 5.412034e-3
CONSTS(33) = 0.001
CONSTS(34) = 0.06
STATES(13) = 9.945511e-1
STATES(14) = 4.051574e-3
STATES(15) = 0.9946
CONSTS(35) = 0.02
STATES(16) = 9.945511e-1
STATES(17) = 4.051574e-3
CONSTS(36) = 0.109625
CONSTS(37) = 0.1
CONSTS(38) = 0.11
CONSTS(39) = 0.89
CONSTS(40) = 0.009
CONSTS(41) = 5.4e-4
CONSTS(42) = 1.5e-8
CONSTS(43) = 2.7e-7
CONSTS(44) = 0.9
CONSTS(45) = 0.1
CONSTS(46) = 0.341
CONSTS(47) = 0.341
CONSTS(48) = 0.75
CONSTS(49) = 0.75
CONSTS(50) = 0.75
CONSTS(51) = 0.75
CONSTS(52) = 1.8
STATES(18) = 7.175662e-6
STATES(19) = 1.000681
STATES(20) = 1.452605e-2
STATES(21) = 2.421991e-2
CONSTS(53) = 9
CONSTS(54) = 0.11
CONSTS(55) = 0.89
CONSTS(56) = 1.57
CONSTS(57) = 12.29
CONSTS(58) = 1.3
CONSTS(59) = 87.5
CONSTS(60) = 0.00359
CONSTS(61) = 0.000256
CONSTS(62) = 0.27
CONSTS(63) = 0.35
CONSTS(64) = 3
CONSTS(65) = 0.11
CONSTS(66) = 0.89
CONSTS(67) = 2.35
CONSTS(68) = 0.0005
CONSTS(69) = 1.6
CONSTS(70) = 0.0673
CONSTS(71) = 0.0002513
CONSTS(72) = 0.11
CONSTS(73) = 0.89
CONSTS(74) = 15
CONSTS(75) = 1
CONSTS(76) = 0.45
STATES(22) = 8.884332e-1
STATES(23) = 1.024274e-7
STATES(24) = 8.156628e-7
CONSTS(77) = 25
CONSTS(78) = 10
CONSTS(79) = 0.06
CONSTS(80) = 0.5
CONSTS(81) = 0.005
CONSTS(82) = 2.5
STATES(25) = 5.545201e-1
CONSTS(83) = 5.348e-6
CONSTS(84) = 5.3114e-3
CONSTS(85) = 2.6
CONSTS(86) = 0.000246
CONSTS(87) = 1.7
CONSTS(88) = 1.787
STATES(26) = 7.720854e-1
STATES(27) = 3.539892
CONSTS(89) = 1.65
CONSTS(90) = 7.561
CONSTS(91) = 0.0001
CONSTS(92) = 1e-3
STATES(28) = 1.110363e-1
STATES(29) = 9.566355e-3
STATES(30) = 7.297378e-2
STATES(31) = 7.347888e-3
STATES(32) = 1.242988
CONSTS(93) = 0.0374
CONSTS(94) = 0.0046
CONSTS(95) = 0.0134
CONSTS(96) = 0.00165
CONSTS(97) = 0.14
CONSTS(98) = 100
CONSTS(99) = 100
CONSTS(100) = 1.3
CONSTS(101) = 30e-3
CONSTS(102) = 65
STATES(33) = 8.773191e-3
STATES(34) = 1.078283e-1
STATES(35) = 1.524002e-2
STATES(36) = 2.911916e-4
STATES(37) = 1.298754e-3
STATES(38) = 1.381982e-1
STATES(39) = 2.143165e-3
CONSTS(103) = 0.07
CONSTS(104) = 0.14
CONSTS(105) = 0.14
CONSTS(106) = 0.024
CONSTS(107) = 0.14
CONSTS(108) = 0.14
CONSTS(109) = 0.0171
CONSTS(110) = 32.7
CONSTS(111) = 2.37
CONSTS(112) = 3e-3
CONSTS(113) = 34
CONSTS(114) = 13.8
CONSTS(115) = 15.7e-3
CONSTS(116) = 100
CONSTS(117) = 19.6e-3
CONSTS(118) = 0.032e-3
CONSTS(119) = 3.33e-3
CONSTS(120) = 238e-3
CONSTS(121) = 0.46e-3
CONSTS(122) = 0.057e-3
CONSTS(123) = 60e-3
CONSTS(124) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(5)/CONSTS(9))
CONSTS(125) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(11)/CONSTS(8))
CONSTS(126) = (EXP(CONSTS(6)/67.3000) - 1.00000)/7.00000
CONSTS(127) = 0.0300000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(128) = 0.900000* (CONSTS(5)/5.40000) ** (1.0 / 2)
CONSTS(129) = CONSTS(52) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(130) = CONSTS(56) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(131) = CONSTS(67) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(132) = CONSTS(85) ** (CONSTS(13) - 310.000)/10.0000
CONSTS(133) = ( 3.14159*CONSTS(17)/1000.00 ** 2.00000*CONSTS(16))/1000.00 ** 3.00000
CONSTS(134) = 0.0350000*CONSTS(133)
CONSTS(135) = 0.0200000*CONSTS(133)
CONSTS(136) = 0.0539000*0.0100000*CONSTS(133)
CONSTS(137) = 0.650000*CONSTS(133)
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(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(15) = 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))))
RATES(9) = ALGBRC(1)*(1.00000 - STATES(9)) - ALGBRC(15)*STATES(9)
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(16) = 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))))
RATES(10) = ALGBRC(2)*(1.00000 - STATES(10)) - ALGBRC(16)*STATES(10)
ALGBRC(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(17) = 0.0800000*EXP(- STATES(1)/11.0000)
RATES(8) = ALGBRC(3)*(1.00000 - STATES(8)) - ALGBRC(17)*STATES(8)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (50.0000+STATES(1))/7.50000))
ALGBRC(18) = 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(18)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(20) = 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(6) - STATES(12))/ALGBRC(20)
ALGBRC(7) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(21) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
RATES(14) = (ALGBRC(7) - STATES(14))/ALGBRC(21)
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(22) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
RATES(13) = (ALGBRC(8) - STATES(13))/ALGBRC(22)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(23) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
RATES(15) = (ALGBRC(9) - STATES(15))/ALGBRC(23)
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(24) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
RATES(17) = (ALGBRC(10) - STATES(17))/ALGBRC(24)
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(25) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
RATES(16) = (ALGBRC(11) - STATES(16))/ALGBRC(25)
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(26) = ( 1.00000*ALGBRC(12)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
RATES(18) = (ALGBRC(12) - STATES(18))/ALGBRC(26)
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(27) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
RATES(19) = (ALGBRC(13) - STATES(19))/ALGBRC(27)
ALGBRC(14) = CONSTS(74) - (CONSTS(74) - CONSTS(75))/(1.00000+CONSTS(76)/STATES(25) ** CONSTS(82))
ALGBRC(28) = CONSTS(78)/ALGBRC(14)
ALGBRC(30) = CONSTS(80)*ALGBRC(14)
RATES(24) = ( ALGBRC(28)*STATES(6) ** 2.00000*STATES(22) - CONSTS(79)*STATES(24)) - ( ALGBRC(30)*STATES(6)*STATES(24) - CONSTS(81)*STATES(23))
ALGBRC(32) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
RATES(22) = ( CONSTS(81)*ALGBRC(32) - ALGBRC(30)*STATES(6)*STATES(22)) - ( ALGBRC(28)*STATES(6) ** 2.00000*STATES(22) - CONSTS(79)*STATES(24))
RATES(23) = ( ALGBRC(30)*STATES(6)*STATES(24) - CONSTS(81)*STATES(23)) - ( CONSTS(79)*STATES(23) - ALGBRC(28)*STATES(6) ** 2.00000*ALGBRC(32))
ALGBRC(81) = CONSTS(91)*STATES(3)*(CONSTS(90) - STATES(27)) - CONSTS(92)*STATES(27)
RATES(27) = ALGBRC(81)
ALGBRC(19) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(3))
ALGBRC(34) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(35) = CONSTS(19)*CONSTS(18)*ALGBRC(34)*(STATES(1) - ALGBRC(19))
ALGBRC(38) = CONSTS(22)*CONSTS(21)*(STATES(1) - ALGBRC(19))
ALGBRC(41) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))+ 0.0365000*CONSTS(126)*EXP(( - STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(42) = ( (( CONSTS(28)*CONSTS(27)*ALGBRC(41))/(1.00000+CONSTS(25)/STATES(3) ** CONSTS(24)))*CONSTS(5))/(CONSTS(5)+CONSTS(26))
ALGBRC(68) = 1.00000 - STATES(21)
ALGBRC(64) = ( 0.450000*STATES(18)*STATES(19)*CONSTS(129)*STATES(1)*CONSTS(14) ** 2.00000)/( CONSTS(12)*CONSTS(13))
ALGBRC(70) = ( ALGBRC(64)*ALGBRC(68)*CONSTS(44)*CONSTS(42)*( CONSTS(48)*STATES(3)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(49)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(76) = 1.00000/(1.00000+CONSTS(61)/STATES(6) ** 3.00000)
ALGBRC(73) = ( EXP(( CONSTS(63)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(3) ** CONSTS(64)*CONSTS(7) - EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(64)*STATES(6))/(1.00000+ CONSTS(62)*EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(77) = ( CONSTS(54)*CONSTS(53)*ALGBRC(76)*CONSTS(130)*ALGBRC(73))/( CONSTS(60)*CONSTS(6) ** CONSTS(64)*(1.00000+STATES(3)/CONSTS(57) ** CONSTS(64))+ CONSTS(59) ** CONSTS(64)*STATES(6)*(1.00000+STATES(6)/CONSTS(60))+ CONSTS(58)*STATES(3) ** CONSTS(64)+ STATES(3) ** CONSTS(64)*CONSTS(7)+ CONSTS(6) ** CONSTS(64)*STATES(6))
ALGBRC(79) = (STATES(3) - STATES(4))*1.83130e-14
RATES(3) = (( - CONSTS(15)*(ALGBRC(35)+ 3.00000*ALGBRC(77)+ALGBRC(38)+ 3.00000*ALGBRC(42)+ALGBRC(70)))/( CONSTS(136)*CONSTS(14)) - ALGBRC(79)/CONSTS(136)) - ALGBRC(81)
ALGBRC(82) = (STATES(4) - STATES(2))*1.63860e-12
RATES(2) = ALGBRC(82)/CONSTS(137)
ALGBRC(84) = CONSTS(91)*STATES(4)*(CONSTS(89) - STATES(26)) - CONSTS(92)*STATES(26)
RATES(26) = ALGBRC(84)
ALGBRC(29) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(4))
ALGBRC(36) = CONSTS(20)*CONSTS(18)*ALGBRC(34)*(STATES(1) - ALGBRC(29))
ALGBRC(39) = CONSTS(23)*CONSTS(21)*(STATES(1) - ALGBRC(29))
ALGBRC(43) = ( (( CONSTS(29)*CONSTS(27)*ALGBRC(41))/(1.00000+CONSTS(25)/STATES(4) ** CONSTS(24)))*CONSTS(5))/(CONSTS(5)+CONSTS(26))
ALGBRC(65) = 1.00000 - STATES(20)
ALGBRC(67) = ( ALGBRC(64)*ALGBRC(65)*CONSTS(45)*CONSTS(42)*( CONSTS(48)*STATES(4)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(49)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(75) = 1.00000/(1.00000+CONSTS(61)/STATES(7) ** 3.00000)
ALGBRC(74) = ( EXP(( CONSTS(63)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(4) ** CONSTS(64)*CONSTS(7) - EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(64)*STATES(7))/(1.00000+ CONSTS(62)*EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(78) = ( CONSTS(55)*CONSTS(53)*ALGBRC(75)*CONSTS(130)*ALGBRC(74))/( CONSTS(60)*CONSTS(6) ** CONSTS(64)*(1.00000+STATES(4)/CONSTS(57) ** CONSTS(64))+ CONSTS(59) ** CONSTS(64)*STATES(7)*(1.00000+STATES(7)/CONSTS(60))+ CONSTS(58)*STATES(4) ** CONSTS(64)+ STATES(4) ** CONSTS(64)*CONSTS(7)+ CONSTS(6) ** CONSTS(64)*STATES(7))
RATES(4) = (( - CONSTS(15)*(ALGBRC(36)+ 3.00000*ALGBRC(78)+ALGBRC(39)+ 3.00000*ALGBRC(43)+ALGBRC(67)))/( CONSTS(135)*CONSTS(14))+(ALGBRC(79) - ALGBRC(82))/CONSTS(135)) - ALGBRC(84)
ALGBRC(37) = ALGBRC(35)+ALGBRC(36)
ALGBRC(40) = ALGBRC(38)+ALGBRC(39)
ALGBRC(44) = ALGBRC(42)+ALGBRC(43)
ALGBRC(45) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(46) = CONSTS(127)*STATES(11)*ALGBRC(45)*(STATES(1) - CONSTS(124))
ALGBRC(47) = - arbitrary_log(STATES(6)/1.00000, 10)+3.00000
ALGBRC(49) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(47))/0.600000)))
ALGBRC(51) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(32)*CONSTS(6))/(CONSTS(9)+ CONSTS(32)*STATES(2)))
ALGBRC(52) = CONSTS(30)*ALGBRC(49)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(51))
ALGBRC(48) = - arbitrary_log(STATES(7)/1.00000, 10)+3.00000
ALGBRC(50) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(48))/0.600000)))
ALGBRC(53) = CONSTS(31)*ALGBRC(50)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(51))
ALGBRC(54) = ALGBRC(52)+ALGBRC(53)
ALGBRC(55) = ( CONSTS(33)*(STATES(1) - CONSTS(124)))/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(56) = CONSTS(34)*STATES(14)*(STATES(13)+ 0.500000*STATES(15))*(STATES(1) - CONSTS(124))
ALGBRC(57) = CONSTS(35)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(124))
ALGBRC(58) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(124)+59.2150))))
ALGBRC(59) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(124))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(124)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(124))+4.75300)))
ALGBRC(60) = ALGBRC(58)/(ALGBRC(58)+ALGBRC(59))
ALGBRC(61) = CONSTS(128)*ALGBRC(60)*(STATES(1) - CONSTS(124))
ALGBRC(80) = ALGBRC(77)+ALGBRC(78)
ALGBRC(62) = CONSTS(36)*(STATES(1) - CONSTS(125))*(CONSTS(38)/(1.00000+CONSTS(37)/STATES(6))+CONSTS(39)/(1.00000+CONSTS(37)/STATES(7)))
ALGBRC(63) = CONSTS(40)*(STATES(1) - CONSTS(125))
ALGBRC(66) = ( ALGBRC(64)*ALGBRC(65)*CONSTS(45)*CONSTS(41)*4.00000*( CONSTS(46)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(47)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(69) = ( ALGBRC(64)*ALGBRC(68)*CONSTS(44)*CONSTS(41)*4.00000*( CONSTS(46)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(47)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(71) = ( ALGBRC(64)*( ALGBRC(65)*CONSTS(45)+ ALGBRC(68)*CONSTS(44))*CONSTS(43)*( CONSTS(50)*CONSTS(9)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(51)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(72) = ALGBRC(66)+ALGBRC(69)+ALGBRC(67)+ALGBRC(70)+ALGBRC(71)
ALGBRC(31) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(6))
ALGBRC(87) = CONSTS(71)*CONSTS(72)*(STATES(1) - ALGBRC(31))
ALGBRC(33) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(7))
ALGBRC(88) = CONSTS(71)*CONSTS(73)*(STATES(1) - ALGBRC(33))
ALGBRC(89) = ALGBRC(88)+ALGBRC(87)
ALGBRC(85) = ( CONSTS(131)*CONSTS(70)*CONSTS(66))/(1.00000+CONSTS(68)/STATES(7) ** CONSTS(69))
ALGBRC(83) = ( CONSTS(131)*CONSTS(70)*CONSTS(65))/(1.00000+CONSTS(68)/STATES(6) ** CONSTS(69))
ALGBRC(86) = ALGBRC(83)+ALGBRC(85)
ALGBRC(5) = 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(37)+ALGBRC(40)+ALGBRC(44)+ALGBRC(46)+ALGBRC(54)+ALGBRC(56)+ALGBRC(57)+ALGBRC(61)+ALGBRC(80)+ALGBRC(62)+ALGBRC(63)+ALGBRC(72)+ALGBRC(89)+ALGBRC(86)+ALGBRC(55)+ALGBRC(5))
ALGBRC(93) = CONSTS(99)*STATES(25)*(( CONSTS(97)*CONSTS(137))/CONSTS(134) - STATES(32)) - CONSTS(102)*STATES(32)
RATES(32) = ALGBRC(93)
ALGBRC(90) = CONSTS(77)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(91) = CONSTS(83)*(STATES(25) - STATES(6))
ALGBRC(92) = ( CONSTS(132)*CONSTS(84)*(STATES(5)/CONSTS(86) ** CONSTS(88) - STATES(25)/CONSTS(87) ** CONSTS(88)))/(1.00000+STATES(5)/CONSTS(86) ** CONSTS(88)+STATES(25)/CONSTS(87) ** CONSTS(88))
RATES(25) = (ALGBRC(92) - (( ALGBRC(91)*CONSTS(137))/CONSTS(134)+ALGBRC(90))) - ALGBRC(93)
ALGBRC(96) = CONSTS(98)*STATES(6)*(( CONSTS(94)*0.100000*CONSTS(137))/CONSTS(136) - STATES(29)) - CONSTS(100)*STATES(29)
RATES(29) = ALGBRC(96)
ALGBRC(98) = CONSTS(98)*STATES(7)*(( CONSTS(93)*CONSTS(137))/CONSTS(135) - STATES(28)) - CONSTS(100)*STATES(28)
RATES(28) = ALGBRC(98)
ALGBRC(99) = CONSTS(98)*STATES(6)*(( CONSTS(96)*0.100000*CONSTS(137))/CONSTS(136) - STATES(31)) - CONSTS(101)*STATES(31)
RATES(31) = ALGBRC(99)
ALGBRC(97) = CONSTS(110)*STATES(5)*(CONSTS(103) - STATES(33)) - CONSTS(117)*STATES(33)
RATES(33) = ALGBRC(97)
ALGBRC(101) = CONSTS(98)*STATES(7)*(( CONSTS(95)*CONSTS(137))/CONSTS(135) - STATES(30)) - CONSTS(101)*STATES(30)
RATES(30) = ALGBRC(101)
ALGBRC(100) = CONSTS(111)*STATES(5)*(CONSTS(104) - (STATES(34)+STATES(35))) - CONSTS(118)*STATES(34)
RATES(34) = ALGBRC(100)
ALGBRC(94) = (STATES(6) - STATES(7))*8.24130e-13
ALGBRC(102) = ALGBRC(96)+ALGBRC(99)
ALGBRC(104) = (ALGBRC(69) - 2.00000*ALGBRC(77))+ALGBRC(87)+ALGBRC(83)
RATES(6) = ((( - ALGBRC(104)*CONSTS(15))/( CONSTS(136)*2.00000*CONSTS(14)) - ALGBRC(94)/CONSTS(136))+( ALGBRC(90)*CONSTS(134))/CONSTS(136)+( ALGBRC(91)*CONSTS(137))/CONSTS(136)) - 1.00000*ALGBRC(102)
ALGBRC(103) = CONSTS(112)*CONSTS(10)*(CONSTS(105) - (STATES(34)+STATES(35))) - CONSTS(119)*STATES(35)
RATES(35) = ALGBRC(103)
ALGBRC(95) = (STATES(7) - STATES(5))*3.72430e-12
ALGBRC(105) = ALGBRC(98)+ALGBRC(101)
ALGBRC(107) = (ALGBRC(66) - 2.00000*ALGBRC(78))+ALGBRC(88)+ALGBRC(85)
RATES(7) = (( - ALGBRC(107)*CONSTS(15))/( CONSTS(135)*2.00000*CONSTS(14))+(ALGBRC(94) - ALGBRC(95))/CONSTS(135)) - 1.00000*ALGBRC(105)
ALGBRC(106) = CONSTS(113)*STATES(5)*(CONSTS(106) - STATES(36)) - CONSTS(120)*STATES(36)
RATES(36) = ALGBRC(106)
ALGBRC(108) = CONSTS(114)*STATES(5)*(CONSTS(107) - (STATES(37)+STATES(38))) - CONSTS(121)*STATES(37)
RATES(37) = ALGBRC(108)
ALGBRC(109) = CONSTS(115)*CONSTS(10)*(CONSTS(108) - (STATES(37)+STATES(38))) - CONSTS(122)*STATES(38)
RATES(38) = ALGBRC(109)
ALGBRC(110) = CONSTS(116)*STATES(5)*(CONSTS(109) - STATES(39)) - CONSTS(123)*STATES(39)
RATES(39) = ALGBRC(110)
ALGBRC(111) = ALGBRC(97)+ALGBRC(100)+ALGBRC(103)+ALGBRC(106)+ALGBRC(108)+ALGBRC(109)+ALGBRC(110)
RATES(5) = (( - ALGBRC(92)*CONSTS(134))/CONSTS(137)+ALGBRC(95)/CONSTS(137)) - 1.00000*ALGBRC(111)
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = TERNRY(STATES(1).LT.- 40.0000, 0.135000*EXP((80.0000+STATES(1))/- 6.80000), 0.00000)
ALGBRC(15) = 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(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(16) = 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(3) = (( 0.320000*(STATES(1)+47.1300))/1.00000)/(1.00000 - EXP( - 0.100000*(STATES(1)+47.1300)))
ALGBRC(17) = 0.0800000*EXP(- STATES(1)/11.0000)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (50.0000+STATES(1))/7.50000))
ALGBRC(18) = 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(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 1.50000)/16.7000))
ALGBRC(20) = 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(7) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(21) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
ALGBRC(8) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(22) = 3000.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(23) = 2800.00/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+220.000
ALGBRC(10) = 1.00000/(1.00000+EXP(- (STATES(1)+3.00000)/15.0000))
ALGBRC(24) = 3.50000*EXP(- STATES(1)/30.0000 ** 2.00000)+1.50000
ALGBRC(11) = 1.00000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))
ALGBRC(25) = 20.0000/(1.00000+EXP((STATES(1)+33.5000)/10.0000))+20.0000
ALGBRC(12) = 1.00000/(1.00000+EXP(- (STATES(1)+14.5000)/6.00000))
ALGBRC(26) = ( 1.00000*ALGBRC(12)*(1.00000 - EXP(- (STATES(1)+14.5000)/6.00000)))/( 0.0350000*(STATES(1)+14.5000))
ALGBRC(13) = 1.00000/(1.00000+EXP((STATES(1)+35.0600)/3.60000))+0.600000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(27) = 1.00000/( 0.0197000*EXP(- 0.0337000*(STATES(1)+14.5000) ** 2.00000)+0.0200000)
ALGBRC(14) = CONSTS(74) - (CONSTS(74) - CONSTS(75))/(1.00000+CONSTS(76)/STATES(25) ** CONSTS(82))
ALGBRC(28) = CONSTS(78)/ALGBRC(14)
ALGBRC(30) = CONSTS(80)*ALGBRC(14)
ALGBRC(32) = ((1.00000 - STATES(22)) - STATES(24)) - STATES(23)
ALGBRC(81) = CONSTS(91)*STATES(3)*(CONSTS(90) - STATES(27)) - CONSTS(92)*STATES(27)
ALGBRC(19) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(3))
ALGBRC(34) = STATES(8) ** 3.00000*STATES(9)*STATES(10)
ALGBRC(35) = CONSTS(19)*CONSTS(18)*ALGBRC(34)*(STATES(1) - ALGBRC(19))
ALGBRC(38) = CONSTS(22)*CONSTS(21)*(STATES(1) - ALGBRC(19))
ALGBRC(41) = 1.00000/(1.00000+ 0.124500*EXP(( - 0.100000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))+ 0.0365000*CONSTS(126)*EXP(( - STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(42) = ( (( CONSTS(28)*CONSTS(27)*ALGBRC(41))/(1.00000+CONSTS(25)/STATES(3) ** CONSTS(24)))*CONSTS(5))/(CONSTS(5)+CONSTS(26))
ALGBRC(68) = 1.00000 - STATES(21)
ALGBRC(64) = ( 0.450000*STATES(18)*STATES(19)*CONSTS(129)*STATES(1)*CONSTS(14) ** 2.00000)/( CONSTS(12)*CONSTS(13))
ALGBRC(70) = ( ALGBRC(64)*ALGBRC(68)*CONSTS(44)*CONSTS(42)*( CONSTS(48)*STATES(3)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(49)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(76) = 1.00000/(1.00000+CONSTS(61)/STATES(6) ** 3.00000)
ALGBRC(73) = ( EXP(( CONSTS(63)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(3) ** CONSTS(64)*CONSTS(7) - EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(64)*STATES(6))/(1.00000+ CONSTS(62)*EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(77) = ( CONSTS(54)*CONSTS(53)*ALGBRC(76)*CONSTS(130)*ALGBRC(73))/( CONSTS(60)*CONSTS(6) ** CONSTS(64)*(1.00000+STATES(3)/CONSTS(57) ** CONSTS(64))+ CONSTS(59) ** CONSTS(64)*STATES(6)*(1.00000+STATES(6)/CONSTS(60))+ CONSTS(58)*STATES(3) ** CONSTS(64)+ STATES(3) ** CONSTS(64)*CONSTS(7)+ CONSTS(6) ** CONSTS(64)*STATES(6))
ALGBRC(79) = (STATES(3) - STATES(4))*1.83130e-14
ALGBRC(82) = (STATES(4) - STATES(2))*1.63860e-12
ALGBRC(84) = CONSTS(91)*STATES(4)*(CONSTS(89) - STATES(26)) - CONSTS(92)*STATES(26)
ALGBRC(29) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log(CONSTS(6)/STATES(4))
ALGBRC(36) = CONSTS(20)*CONSTS(18)*ALGBRC(34)*(STATES(1) - ALGBRC(29))
ALGBRC(39) = CONSTS(23)*CONSTS(21)*(STATES(1) - ALGBRC(29))
ALGBRC(43) = ( (( CONSTS(29)*CONSTS(27)*ALGBRC(41))/(1.00000+CONSTS(25)/STATES(4) ** CONSTS(24)))*CONSTS(5))/(CONSTS(5)+CONSTS(26))
ALGBRC(65) = 1.00000 - STATES(20)
ALGBRC(67) = ( ALGBRC(64)*ALGBRC(65)*CONSTS(45)*CONSTS(42)*( CONSTS(48)*STATES(4)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(49)*CONSTS(6)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(75) = 1.00000/(1.00000+CONSTS(61)/STATES(7) ** 3.00000)
ALGBRC(74) = ( EXP(( CONSTS(63)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*STATES(4) ** CONSTS(64)*CONSTS(7) - EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13)))*CONSTS(6) ** CONSTS(64)*STATES(7))/(1.00000+ CONSTS(62)*EXP(( (CONSTS(63) - 1.00000)*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))))
ALGBRC(78) = ( CONSTS(55)*CONSTS(53)*ALGBRC(75)*CONSTS(130)*ALGBRC(74))/( CONSTS(60)*CONSTS(6) ** CONSTS(64)*(1.00000+STATES(4)/CONSTS(57) ** CONSTS(64))+ CONSTS(59) ** CONSTS(64)*STATES(7)*(1.00000+STATES(7)/CONSTS(60))+ CONSTS(58)*STATES(4) ** CONSTS(64)+ STATES(4) ** CONSTS(64)*CONSTS(7)+ CONSTS(6) ** CONSTS(64)*STATES(7))
ALGBRC(37) = ALGBRC(35)+ALGBRC(36)
ALGBRC(40) = ALGBRC(38)+ALGBRC(39)
ALGBRC(44) = ALGBRC(42)+ALGBRC(43)
ALGBRC(45) = 1.00000/(1.00000+EXP((33.0000+STATES(1))/22.4000))
ALGBRC(46) = CONSTS(127)*STATES(11)*ALGBRC(45)*(STATES(1) - CONSTS(124))
ALGBRC(47) = - arbitrary_log(STATES(6)/1.00000, 10)+3.00000
ALGBRC(49) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(47))/0.600000)))
ALGBRC(51) = (( CONSTS(12)*CONSTS(13))/CONSTS(14))*log((CONSTS(5)+ CONSTS(32)*CONSTS(6))/(CONSTS(9)+ CONSTS(32)*STATES(2)))
ALGBRC(52) = CONSTS(30)*ALGBRC(49)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(51))
ALGBRC(48) = - arbitrary_log(STATES(7)/1.00000, 10)+3.00000
ALGBRC(50) = 0.0700000*(0.0570000+0.190000/(1.00000+EXP((- 7.20000+ALGBRC(48))/0.600000)))
ALGBRC(53) = CONSTS(31)*ALGBRC(50)*STATES(12) ** 2.00000*(STATES(1) - ALGBRC(51))
ALGBRC(54) = ALGBRC(52)+ALGBRC(53)
ALGBRC(55) = ( CONSTS(33)*(STATES(1) - CONSTS(124)))/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(56) = CONSTS(34)*STATES(14)*(STATES(13)+ 0.500000*STATES(15))*(STATES(1) - CONSTS(124))
ALGBRC(57) = CONSTS(35)*STATES(17)*STATES(16)*(STATES(1) - CONSTS(124))
ALGBRC(58) = 1.02000/(1.00000+EXP( 0.238500*(STATES(1) - (CONSTS(124)+59.2150))))
ALGBRC(59) = ( 0.491240*EXP( 0.0803200*((STATES(1) - CONSTS(124))+5.47600))+ 1.00000*EXP( 0.0617500*(STATES(1) - (CONSTS(124)+594.310))))/(1.00000+EXP( - 0.514300*((STATES(1) - CONSTS(124))+4.75300)))
ALGBRC(60) = ALGBRC(58)/(ALGBRC(58)+ALGBRC(59))
ALGBRC(61) = CONSTS(128)*ALGBRC(60)*(STATES(1) - CONSTS(124))
ALGBRC(80) = ALGBRC(77)+ALGBRC(78)
ALGBRC(62) = CONSTS(36)*(STATES(1) - CONSTS(125))*(CONSTS(38)/(1.00000+CONSTS(37)/STATES(6))+CONSTS(39)/(1.00000+CONSTS(37)/STATES(7)))
ALGBRC(63) = CONSTS(40)*(STATES(1) - CONSTS(125))
ALGBRC(66) = ( ALGBRC(64)*ALGBRC(65)*CONSTS(45)*CONSTS(41)*4.00000*( CONSTS(46)*STATES(7)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(47)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(69) = ( ALGBRC(64)*ALGBRC(68)*CONSTS(44)*CONSTS(41)*4.00000*( CONSTS(46)*STATES(6)*EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(47)*CONSTS(7)))/(EXP(( 2.00000*STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(71) = ( ALGBRC(64)*( ALGBRC(65)*CONSTS(45)+ ALGBRC(68)*CONSTS(44))*CONSTS(43)*( CONSTS(50)*CONSTS(9)*EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - CONSTS(51)*CONSTS(5)))/(EXP(( STATES(1)*CONSTS(14))/( CONSTS(12)*CONSTS(13))) - 1.00000)
ALGBRC(72) = ALGBRC(66)+ALGBRC(69)+ALGBRC(67)+ALGBRC(70)+ALGBRC(71)
ALGBRC(31) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(6))
ALGBRC(87) = CONSTS(71)*CONSTS(72)*(STATES(1) - ALGBRC(31))
ALGBRC(33) = (( CONSTS(12)*CONSTS(13))/( 2.00000*CONSTS(14)))*log(CONSTS(7)/STATES(7))
ALGBRC(88) = CONSTS(71)*CONSTS(73)*(STATES(1) - ALGBRC(33))
ALGBRC(89) = ALGBRC(88)+ALGBRC(87)
ALGBRC(85) = ( CONSTS(131)*CONSTS(70)*CONSTS(66))/(1.00000+CONSTS(68)/STATES(7) ** CONSTS(69))
ALGBRC(83) = ( CONSTS(131)*CONSTS(70)*CONSTS(65))/(1.00000+CONSTS(68)/STATES(6) ** CONSTS(69))
ALGBRC(86) = ALGBRC(83)+ALGBRC(85)
ALGBRC(5) = 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(93) = CONSTS(99)*STATES(25)*(( CONSTS(97)*CONSTS(137))/CONSTS(134) - STATES(32)) - CONSTS(102)*STATES(32)
ALGBRC(90) = CONSTS(77)*STATES(24)*(STATES(25) - STATES(6))
ALGBRC(91) = CONSTS(83)*(STATES(25) - STATES(6))
ALGBRC(92) = ( CONSTS(132)*CONSTS(84)*(STATES(5)/CONSTS(86) ** CONSTS(88) - STATES(25)/CONSTS(87) ** CONSTS(88)))/(1.00000+STATES(5)/CONSTS(86) ** CONSTS(88)+STATES(25)/CONSTS(87) ** CONSTS(88))
ALGBRC(96) = CONSTS(98)*STATES(6)*(( CONSTS(94)*0.100000*CONSTS(137))/CONSTS(136) - STATES(29)) - CONSTS(100)*STATES(29)
ALGBRC(98) = CONSTS(98)*STATES(7)*(( CONSTS(93)*CONSTS(137))/CONSTS(135) - STATES(28)) - CONSTS(100)*STATES(28)
ALGBRC(99) = CONSTS(98)*STATES(6)*(( CONSTS(96)*0.100000*CONSTS(137))/CONSTS(136) - STATES(31)) - CONSTS(101)*STATES(31)
ALGBRC(97) = CONSTS(110)*STATES(5)*(CONSTS(103) - STATES(33)) - CONSTS(117)*STATES(33)
ALGBRC(101) = CONSTS(98)*STATES(7)*(( CONSTS(95)*CONSTS(137))/CONSTS(135) - STATES(30)) - CONSTS(101)*STATES(30)
ALGBRC(100) = CONSTS(111)*STATES(5)*(CONSTS(104) - (STATES(34)+STATES(35))) - CONSTS(118)*STATES(34)
ALGBRC(94) = (STATES(6) - STATES(7))*8.24130e-13
ALGBRC(102) = ALGBRC(96)+ALGBRC(99)
ALGBRC(104) = (ALGBRC(69) - 2.00000*ALGBRC(77))+ALGBRC(87)+ALGBRC(83)
ALGBRC(103) = CONSTS(112)*CONSTS(10)*(CONSTS(105) - (STATES(34)+STATES(35))) - CONSTS(119)*STATES(35)
ALGBRC(95) = (STATES(7) - STATES(5))*3.72430e-12
ALGBRC(105) = ALGBRC(98)+ALGBRC(101)
ALGBRC(107) = (ALGBRC(66) - 2.00000*ALGBRC(78))+ALGBRC(88)+ALGBRC(85)
ALGBRC(106) = CONSTS(113)*STATES(5)*(CONSTS(106) - STATES(36)) - CONSTS(120)*STATES(36)
ALGBRC(108) = CONSTS(114)*STATES(5)*(CONSTS(107) - (STATES(37)+STATES(38))) - CONSTS(121)*STATES(37)
ALGBRC(109) = CONSTS(115)*CONSTS(10)*(CONSTS(108) - (STATES(37)+STATES(38))) - CONSTS(122)*STATES(38)
ALGBRC(110) = CONSTS(116)*STATES(5)*(CONSTS(109) - STATES(39)) - CONSTS(123)*STATES(39)
ALGBRC(111) = ALGBRC(97)+ALGBRC(100)+ALGBRC(103)+ALGBRC(106)+ALGBRC(108)+ALGBRC(109)+ALGBRC(110)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END