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 118 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 111 entries in the constant variable array.
C
C
C VOI is t in component environment (millisecond).
C STATES(1) is V in component membrane (millivolt).
C CONSTS(1) is R in component membrane (joule_per_kmole_kelvin).
C CONSTS(2) is T in component membrane (kelvin).
C CONSTS(3) is Frdy in component membrane (coulomb_per_mole).
C CONSTS(96) is FoRT in component membrane (per_milliV).
C CONSTS(4) is C_mem in component membrane (farad).
C ALGBRC(97) is I_Na_tot_junc in component Sodium_Concentrations (picoA_per_picoF).
C ALGBRC(103) is I_Na_tot_sl in component Sodium_Concentrations (picoA_per_picoF).
C ALGBRC(58) is I_ClCa in component Ca_Activated_Cl_Current (picoA_per_picoF).
C ALGBRC(59) is I_Cl_Bk in component Background_Cl_Current (picoA_per_picoF).
C ALGBRC(112) is I_Ca_tot_junc in component Calcium_Concentrations (picoA_per_picoF).
C ALGBRC(115) is I_Ca_tot_sl in component Calcium_Concentrations (picoA_per_picoF).
C ALGBRC(109) is I_K_tot in component Potassium_Concentrations (picoA_per_picoF).
C ALGBRC(104) is I_Na_tot in component membrane (picoA_per_picoF).
C ALGBRC(60) is I_Cl_tot in component membrane (picoA_per_picoF).
C ALGBRC(117) is I_Ca_tot in component membrane (picoA_per_picoF).
C ALGBRC(118) is I_tot in component membrane (picoA_per_picoF).
C ALGBRC(11) is I_Stim in component membrane (picoA_per_picoF).
C CONSTS(5) is I_Stim_Start in component membrane (millisecond).
C CONSTS(6) is I_Stim_End in component membrane (millisecond).
C CONSTS(7) is I_Stim_Amplitude in component membrane (picoA_per_picoF).
C CONSTS(8) is I_Stim_Period in component membrane (millisecond).
C CONSTS(9) is I_Stim_PulseDuration in component membrane (millisecond).
C CONSTS(10) is F_junc in component membrane (dimensionless).
C CONSTS(97) is F_sl in component membrane (dimensionless).
C CONSTS(11) is F_junc_Ca_L in component membrane (dimensionless).
C CONSTS(98) is F_sl_Ca_L in component membrane (dimensionless).
C CONSTS(99) is V_cell in component membrane (litre).
C CONSTS(12) is Length_cell in component membrane (decimetre).
C CONSTS(13) is Radius_cell in component membrane (decimetre).
C CONSTS(101) is V_myo in component membrane (litre).
C CONSTS(104) is V_sl in component membrane (litre).
C CONSTS(105) is V_sr in component membrane (litre).
C CONSTS(106) is V_junc in component membrane (litre).
C CONSTS(14) is J_Ca_juncsl in component membrane (L_per_millis).
C CONSTS(15) is J_Ca_slmyo in component membrane (L_per_millis).
C CONSTS(16) is J_Na_juncsl in component membrane (L_per_millis).
C CONSTS(17) is J_Na_slmyo in component membrane (L_per_millis).
C ALGBRC(95) is I_Na_junc in component Fast_Na_Current (picoA_per_picoF).
C ALGBRC(99) is I_Na_sl in component Fast_Na_Current (picoA_per_picoF).
C ALGBRC(100) is I_Na in component Fast_Na_Current (picoA_per_picoF).
C ALGBRC(98) is E_Na_sl in component Nerst_Potentials (millivolt).
C ALGBRC(94) is E_Na_junc in component Nerst_Potentials (millivolt).
C CONSTS(18) is G_Na in component Fast_Na_Current (milliS_per_microF).
C STATES(2) is m in component Fast_Na_Current (dimensionless).
C ALGBRC(1) is m_ss in component Fast_Na_Current (dimensionless).
C ALGBRC(13) is tau_m in component Fast_Na_Current (millisecond).
C ALGBRC(2) is a_h in component Fast_Na_Current (dimensionless).
C ALGBRC(14) is b_h in component Fast_Na_Current (dimensionless).
C ALGBRC(27) is tau_h in component Fast_Na_Current (millisecond).
C STATES(3) is h in component Fast_Na_Current (dimensionless).
C ALGBRC(34) is h_ss in component Fast_Na_Current (dimensionless).
C ALGBRC(3) is a_j in component Fast_Na_Current (dimensionless).
C ALGBRC(15) is b_j in component Fast_Na_Current (dimensionless).
C ALGBRC(28) is tau_j in component Fast_Na_Current (millisecond).
C STATES(4) is j in component Fast_Na_Current (dimensionless).
C ALGBRC(35) is j_ss in component Fast_Na_Current (dimensionless).
C ALGBRC(96) is I_Na_Bk_junc in component Background_Na_Current (picoA_per_picoF).
C ALGBRC(101) is I_Na_Bk_sl in component Background_Na_Current (picoA_per_picoF).
C ALGBRC(102) is I_Na_Bk in component Background_Na_Current (picoA_per_picoF).
C CONSTS(19) is G_Na_B in component Background_Na_Current (milliS_per_microF).
C CONSTS(20) is Na_o in component Sodium_Concentrations (millimolar).
C CONSTS(100) is sigma in component Na_K_Pump_Current (dimensionless).
C ALGBRC(25) is f_NaK in component Na_K_Pump_Current (dimensionless).
C ALGBRC(32) is I_NaK_junc in component Na_K_Pump_Current (picoA_per_picoF).
C ALGBRC(39) is I_NaK_sl in component Na_K_Pump_Current (picoA_per_picoF).
C ALGBRC(44) is I_NaK in component Na_K_Pump_Current (picoA_per_picoF).
C CONSTS(21) is Ibar_NaK in component Na_K_Pump_Current (picoA_per_picoF).
C CONSTS(22) is K_o in component Potassium_Concentrations (millimolar).
C STATES(5) is Na_j in component Sodium_Concentrations (millimolar).
C STATES(6) is Na_sl in component Sodium_Concentrations (millimolar).
C CONSTS(23) is Km_Ko in component Na_K_Pump_Current (millimolar).
C CONSTS(24) is Km_Naip in component Na_K_Pump_Current (millimolar).
C ALGBRC(4) is x_r_ss in component Rapidly_Activating_K_Current (dimensionless).
C STATES(7) is x_Kr in component Rapidly_Activating_K_Current (dimensionless).
C ALGBRC(45) is r_Kr in component Rapidly_Activating_K_Current (dimensionless).
C CONSTS(102) is E_K in component Nerst_Potentials (millivolt).
C CONSTS(25) is G_Kr in component Rapidly_Activating_K_Current (milliS_per_microF).
C ALGBRC(16) is tau_xr in component Rapidly_Activating_K_Current (millisecond).
C ALGBRC(46) is I_Kr in component Rapidly_Activating_K_Current (picoA_per_picoF).
C ALGBRC(5) is x_s_ss in component Slowly_Activating_K_Current (dimensionless).
C ALGBRC(17) is tau_xs in component Slowly_Activating_K_Current (millisecond).
C STATES(8) is x_Ks in component Slowly_Activating_K_Current (dimensionless).
C ALGBRC(105) is E_Ks in component Nerst_Potentials (millivolt).
C CONSTS(26) is G_Ks_junc in component Slowly_Activating_K_Current (milliS_per_microF).
C CONSTS(27) is G_Ks_sl in component Slowly_Activating_K_Current (milliS_per_microF).
C ALGBRC(106) is I_Ks_junc in component Slowly_Activating_K_Current (picoA_per_picoF).
C ALGBRC(107) is I_Ks_sl in component Slowly_Activating_K_Current (picoA_per_picoF).
C ALGBRC(108) is I_Ks in component Slowly_Activating_K_Current (picoA_per_picoF).
C ALGBRC(47) is kp_Kp in component Plateau_K_Current (dimensionless).
C CONSTS(28) is G_Kp in component Plateau_K_Current (milliS_per_microF).
C ALGBRC(48) is I_Kp in component Plateau_K_Current (picoA_per_picoF).
C ALGBRC(6) is x_to_ss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(7) is y_to_ss in component Transient_Outward_K_Current (dimensionless).
C ALGBRC(18) is tau_x_tos in component Transient_Outward_K_Current (millisecond).
C ALGBRC(19) is tau_y_tos in component Transient_Outward_K_Current (millisecond).
C ALGBRC(20) is tau_x_tof in component Transient_Outward_K_Current (millisecond).
C ALGBRC(21) is tau_y_tof in component Transient_Outward_K_Current (millisecond).
C STATES(9) is x_to_s in component Transient_Outward_K_Current (dimensionless).
C STATES(10) is y_to_s in component Transient_Outward_K_Current (dimensionless).
C STATES(11) is x_to_f in component Transient_Outward_K_Current (dimensionless).
C STATES(12) is y_to_f in component Transient_Outward_K_Current (dimensionless).
C CONSTS(29) is G_to_s in component Transient_Outward_K_Current (milliS_per_microF).
C CONSTS(30) is G_to_f in component Transient_Outward_K_Current (milliS_per_microF).
C ALGBRC(49) is I_to_s in component Transient_Outward_K_Current (picoA_per_picoF).
C ALGBRC(50) is I_to_f in component Transient_Outward_K_Current (picoA_per_picoF).
C ALGBRC(51) is I_to in component Transient_Outward_K_Current (picoA_per_picoF).
C ALGBRC(52) is a_K1 in component Inward_Rectifier_K_Current (dimensionless).
C ALGBRC(53) is b_K1 in component Inward_Rectifier_K_Current (dimensionless).
C ALGBRC(54) is K1_ss in component Inward_Rectifier_K_Current (dimensionless).
C ALGBRC(55) is I_K1 in component Inward_Rectifier_K_Current (picoA_per_picoF).
C CONSTS(31) is G_K1 in component Inward_Rectifier_K_Current (milliS_per_microF).
C ALGBRC(56) is I_ClCa_junc in component Ca_Activated_Cl_Current (picoA_per_picoF).
C ALGBRC(57) is I_ClCa_sl in component Ca_Activated_Cl_Current (picoA_per_picoF).
C STATES(13) is Ca_j in component Calcium_Concentrations (millimolar).
C STATES(14) is Ca_sl in component Calcium_Concentrations (millimolar).
C CONSTS(103) is E_Cl in component Nerst_Potentials (millivolt).
C CONSTS(32) is G_ClCa in component Ca_Activated_Cl_Current (milliS_per_microF).
C CONSTS(33) is Kd_ClCa in component Ca_Activated_Cl_Current (millimolar).
C CONSTS(34) is G_Cl_B in component Background_Cl_Current (milliS_per_microF).
C ALGBRC(8) is d_ss in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(22) is alpha_d in component L_Type_Calcium_Current (millisecond).
C ALGBRC(29) is beta_d in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(36) is gamma_d in component L_Type_Calcium_Current (millisecond).
C ALGBRC(41) is tau_d in component L_Type_Calcium_Current (millisecond).
C ALGBRC(9) is f_ss in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(23) is alpha_f in component L_Type_Calcium_Current (millisecond).
C ALGBRC(30) is beta_f in component L_Type_Calcium_Current (millisecond).
C ALGBRC(37) is gamma_f in component L_Type_Calcium_Current (millisecond).
C ALGBRC(42) is tau_f in component L_Type_Calcium_Current (millisecond).
C ALGBRC(10) is f_2_ss in component L_Type_Calcium_Current (dimensionless).
C ALGBRC(24) is alpha_f_2 in component L_Type_Calcium_Current (millisecond).
C ALGBRC(31) is beta_f_2 in component L_Type_Calcium_Current (millisecond).
C ALGBRC(38) is gamma_f_2 in component L_Type_Calcium_Current (millisecond).
C ALGBRC(43) is tau_f_2 in component L_Type_Calcium_Current (millisecond).
C STATES(15) is d in component L_Type_Calcium_Current (dimensionless).
C STATES(16) is f in component L_Type_Calcium_Current (dimensionless).
C STATES(17) is f_2 in component L_Type_Calcium_Current (dimensionless).
C STATES(18) is f_Ca_B_j in component L_Type_Calcium_Current (dimensionless).
C STATES(19) is f_Ca_B_sl in component L_Type_Calcium_Current (dimensionless).
C CONSTS(35) is K_i in component Potassium_Concentrations (millimolar).
C CONSTS(36) is Ca_o in component Calcium_Concentrations (millimolar).
C ALGBRC(61) is Ibar_Ca_j in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(62) is Ibar_Ca_sl in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(63) is Ibar_K in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(64) is Ibar_Na_j in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(65) is Ibar_Na_sl in component L_Type_Calcium_Current (picoA_per_picoF).
C CONSTS(37) is p_Ca in component L_Type_Calcium_Current (litre_per_farad_millisecond).
C CONSTS(38) is p_Na in component L_Type_Calcium_Current (litre_per_farad_millisecond).
C CONSTS(39) is p_K in component L_Type_Calcium_Current (litre_per_farad_millisecond).
C ALGBRC(66) is I_Ca_junc in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(67) is I_Ca_sl in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(68) is I_Ca in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(70) is I_Ca_Na_junc in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(71) is I_Ca_Na_sl in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(73) is I_Ca_Na in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(69) is I_Ca_K in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(75) is I_Ca_L in component L_Type_Calcium_Current (picoA_per_picoF).
C ALGBRC(72) is Ka_junc in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(74) is Ka_sl in component Na_Ca_Exchanger_Current (dimensionless).
C CONSTS(40) is Kd_act in component Na_Ca_Exchanger_Current (millimolar).
C CONSTS(41) is Km_Ca_i in component Na_Ca_Exchanger_Current (millimolar).
C CONSTS(42) is Km_Ca_o in component Na_Ca_Exchanger_Current (millimolar).
C CONSTS(43) is Km_Na_i in component Na_Ca_Exchanger_Current (millimolar).
C CONSTS(44) is Km_Na_o in component Na_Ca_Exchanger_Current (millimolar).
C CONSTS(45) is k_sat in component Na_Ca_Exchanger_Current (dimensionless).
C CONSTS(46) is nu in component Na_Ca_Exchanger_Current (dimensionless).
C ALGBRC(76) is s1_junc in component Na_Ca_Exchanger_Current (millimolar4).
C ALGBRC(78) is s2_junc in component Na_Ca_Exchanger_Current (millimolar4).
C ALGBRC(80) is s3_junc in component Na_Ca_Exchanger_Current (millimolar4).
C ALGBRC(77) is s1_sl in component Na_Ca_Exchanger_Current (millimolar4).
C ALGBRC(79) is s2_sl in component Na_Ca_Exchanger_Current (millimolar4).
C ALGBRC(81) is s3_sl in component Na_Ca_Exchanger_Current (millimolar4).
C CONSTS(47) is Ibar_NCX in component Na_Ca_Exchanger_Current (picoA_per_picoF).
C ALGBRC(83) is I_ncx_sl in component Na_Ca_Exchanger_Current (picoA_per_picoF).
C ALGBRC(82) is I_ncx_junc in component Na_Ca_Exchanger_Current (picoA_per_picoF).
C ALGBRC(85) is I_ncx in component Na_Ca_Exchanger_Current (picoA_per_picoF).
C CONSTS(48) is Km_P_Ca in component Sarcolemmal_Ca_Pump_Current (millimolar).
C CONSTS(49) is Ibar_PMCA in component Sarcolemmal_Ca_Pump_Current (picoA_per_picoF).
C ALGBRC(84) is I_pCa_junc in component Sarcolemmal_Ca_Pump_Current (picoA_per_picoF).
C ALGBRC(86) is I_pCa_sl in component Sarcolemmal_Ca_Pump_Current (picoA_per_picoF).
C ALGBRC(88) is I_pCa in component Sarcolemmal_Ca_Pump_Current (picoA_per_picoF).
C ALGBRC(113) is E_Ca_sl in component Nerst_Potentials (millivolt).
C ALGBRC(110) is E_Ca_junc in component Nerst_Potentials (millivolt).
C CONSTS(50) is G_Ca_B in component Background_Ca_Current (milliS_per_microF).
C ALGBRC(114) is I_Ca_Bk_sl in component Background_Ca_Current (picoA_per_picoF).
C ALGBRC(111) is I_Ca_Bk_junc in component Background_Ca_Current (picoA_per_picoF).
C ALGBRC(116) is I_Ca_Bk in component Background_Ca_Current (picoA_per_picoF).
C CONSTS(51) is Max_SR in component SR_Fluxes (dimensionless).
C CONSTS(52) is Min_SR in component SR_Fluxes (dimensionless).
C STATES(20) is Ca_SR in component Calcium_Concentrations (millimolar).
C CONSTS(53) is ec50_SR in component SR_Fluxes (millimolar).
C ALGBRC(12) is k_Ca_SR in component SR_Fluxes (dimensionless).
C CONSTS(54) is ko_Ca in component SR_Fluxes (per_milliM2_per_millis).
C ALGBRC(26) is ko_SR_Ca in component SR_Fluxes (per_milliM2_per_millis).
C CONSTS(55) is ki_Ca in component SR_Fluxes (per_milliM_per_millis).
C ALGBRC(33) is ki_SR_Ca in component SR_Fluxes (per_milliM_per_millis).
C CONSTS(56) is ki_m in component SR_Fluxes (per_millis).
C CONSTS(57) is ko_m in component SR_Fluxes (per_millis).
C CONSTS(58) is ks in component SR_Fluxes (per_millis).
C CONSTS(59) is Km_f in component SR_Fluxes (millimolar).
C CONSTS(60) is Km_r in component SR_Fluxes (millimolar).
C ALGBRC(87) is J_SR_Ca_rel in component SR_Fluxes (milliM_per_millis).
C ALGBRC(90) is J_SR_leak in component SR_Fluxes (milliM_per_millis).
C ALGBRC(89) is J_ser_Ca in component SR_Fluxes (milliM_per_millis).
C ALGBRC(40) is RI in component SR_Fluxes (dimensionless).
C STATES(21) is Ry_Rr in component SR_Fluxes (dimensionless).
C STATES(22) is Ry_Ro in component SR_Fluxes (dimensionless).
C STATES(23) is Ry_Ri in component SR_Fluxes (dimensionless).
C STATES(24) is Ca_i in component Calcium_Concentrations (millimolar).
C CONSTS(61) is V_max_SR_CaP in component SR_Fluxes (milliM_per_millis).
C CONSTS(62) is hill_SR_CaP in component SR_Fluxes (dimensionless).
C CONSTS(63) is Mg_i in component Magnesium_Concentrations (millimolar).
C STATES(25) is TnC_l in component Cytosolic_Ca_Buffers (millimolar).
C STATES(26) is TnC_h_c in component Cytosolic_Ca_Buffers (millimolar).
C STATES(27) is TnC_h_m in component Cytosolic_Ca_Buffers (millimolar).
C STATES(28) is CaM in component Cytosolic_Ca_Buffers (millimolar).
C STATES(29) is Myo_c in component Cytosolic_Ca_Buffers (millimolar).
C STATES(30) is Myo_m in component Cytosolic_Ca_Buffers (millimolar).
C STATES(31) is SRB in component Cytosolic_Ca_Buffers (millimolar).
C CONSTS(64) is k_off_TnC_l in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(65) is k_off_TnC_h_Ca in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(66) is k_off_TnC_h_Mg in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(67) is k_off_CaM in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(68) is k_off_myo_Ca in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(69) is k_off_myo_Mg in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(70) is k_off_SR in component Cytosolic_Ca_Buffers (per_millis).
C CONSTS(71) is k_on_TnC_l in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(72) is k_on_TnC_h_Ca in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(73) is k_on_TnC_h_Mg in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(74) is k_on_CaM in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(75) is k_on_myo_Ca in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(76) is k_on_myo_Mg in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(77) is k_on_SR in component Cytosolic_Ca_Buffers (per_milliM_per_millis).
C CONSTS(78) is B_max_TnC_low in component Cytosolic_Ca_Buffers (millimolar).
C CONSTS(79) is B_max_TnC_high in component Cytosolic_Ca_Buffers (millimolar).
C CONSTS(80) is B_max_CaM in component Cytosolic_Ca_Buffers (millimolar).
C CONSTS(81) is B_max_myosin in component Cytosolic_Ca_Buffers (millimolar).
C CONSTS(82) is B_max_SR in component Cytosolic_Ca_Buffers (millimolar).
C ALGBRC(93) is J_Ca_B_cytosol in component Cytosolic_Ca_Buffers (milliM_per_millis).
C CONSTS(107) is B_max_SL_low_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C CONSTS(108) is B_max_SL_low_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C CONSTS(109) is B_max_SL_high_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C CONSTS(110) is B_max_SL_high_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C CONSTS(83) is k_off_sl_l in component Junctional_and_SL_Ca_Buffers (per_millis).
C CONSTS(84) is k_off_sl_h in component Junctional_and_SL_Ca_Buffers (per_millis).
C CONSTS(85) is k_on_sl_l in component Junctional_and_SL_Ca_Buffers (per_milliM_per_millis).
C CONSTS(86) is k_on_sl_h in component Junctional_and_SL_Ca_Buffers (per_milliM_per_millis).
C STATES(32) is SLL_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C STATES(33) is SLL_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C STATES(34) is SLH_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C STATES(35) is SLH_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C ALGBRC(91) is J_Ca_B_junction in component Junctional_and_SL_Ca_Buffers (milliM_per_millis).
C ALGBRC(92) is J_Ca_B_sl in component Junctional_and_SL_Ca_Buffers (milliM_per_millis).
C CONSTS(87) is B_max_Na_j in component Sodium_Concentrations (millimolar).
C CONSTS(88) is B_max_Na_sl in component Sodium_Concentrations (millimolar).
C CONSTS(89) is k_off_Na in component Sodium_Concentrations (per_millis).
C CONSTS(90) is k_on_Na in component Sodium_Concentrations (per_milliM_per_millis).
C STATES(36) is Na_i in component Sodium_Concentrations (millimolar).
C STATES(37) is Na_B_j in component Sodium_Concentrations (millimolar).
C STATES(38) is Na_B_sl in component Sodium_Concentrations (millimolar).
C CONSTS(111) is B_max_csqn in component Calcium_Concentrations (millimolar).
C STATES(39) is Csqn_b in component Calcium_Concentrations (millimolar).
C CONSTS(91) is k_off_csqn in component Calcium_Concentrations (per_millis).
C CONSTS(92) is k_on_csqn in component Calcium_Concentrations (per_milliM_per_millis).
C CONSTS(93) is Cl_i in component Chlorine_Concentrations (millimolar).
C CONSTS(94) is Cl_o in component Chlorine_Concentrations (millimolar).
C CONSTS(95) is p_Na_K in component Nerst_Potentials (dimensionless).
C RATES(1) is d/dt V in component membrane (millivolt).
C RATES(2) is d/dt m in component Fast_Na_Current (dimensionless).
C RATES(3) is d/dt h in component Fast_Na_Current (dimensionless).
C RATES(4) is d/dt j in component Fast_Na_Current (dimensionless).
C RATES(7) is d/dt x_Kr in component Rapidly_Activating_K_Current (dimensionless).
C RATES(8) is d/dt x_Ks in component Slowly_Activating_K_Current (dimensionless).
C RATES(9) is d/dt x_to_s in component Transient_Outward_K_Current (dimensionless).
C RATES(10) is d/dt y_to_s in component Transient_Outward_K_Current (dimensionless).
C RATES(11) is d/dt x_to_f in component Transient_Outward_K_Current (dimensionless).
C RATES(12) is d/dt y_to_f in component Transient_Outward_K_Current (dimensionless).
C RATES(15) is d/dt d in component L_Type_Calcium_Current (dimensionless).
C RATES(16) is d/dt f in component L_Type_Calcium_Current (dimensionless).
C RATES(17) is d/dt f_2 in component L_Type_Calcium_Current (dimensionless).
C RATES(18) is d/dt f_Ca_B_j in component L_Type_Calcium_Current (dimensionless).
C RATES(19) is d/dt f_Ca_B_sl in component L_Type_Calcium_Current (dimensionless).
C RATES(21) is d/dt Ry_Rr in component SR_Fluxes (dimensionless).
C RATES(22) is d/dt Ry_Ro in component SR_Fluxes (dimensionless).
C RATES(23) is d/dt Ry_Ri in component SR_Fluxes (dimensionless).
C RATES(25) is d/dt TnC_l in component Cytosolic_Ca_Buffers (millimolar).
C RATES(26) is d/dt TnC_h_c in component Cytosolic_Ca_Buffers (millimolar).
C RATES(27) is d/dt TnC_h_m in component Cytosolic_Ca_Buffers (millimolar).
C RATES(28) is d/dt CaM in component Cytosolic_Ca_Buffers (millimolar).
C RATES(29) is d/dt Myo_c in component Cytosolic_Ca_Buffers (millimolar).
C RATES(30) is d/dt Myo_m in component Cytosolic_Ca_Buffers (millimolar).
C RATES(31) is d/dt SRB in component Cytosolic_Ca_Buffers (millimolar).
C RATES(32) is d/dt SLL_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C RATES(33) is d/dt SLL_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C RATES(34) is d/dt SLH_j in component Junctional_and_SL_Ca_Buffers (millimolar).
C RATES(35) is d/dt SLH_sl in component Junctional_and_SL_Ca_Buffers (millimolar).
C RATES(5) is d/dt Na_j in component Sodium_Concentrations (millimolar).
C RATES(37) is d/dt Na_B_j in component Sodium_Concentrations (millimolar).
C RATES(6) is d/dt Na_sl in component Sodium_Concentrations (millimolar).
C RATES(38) is d/dt Na_B_sl in component Sodium_Concentrations (millimolar).
C RATES(36) is d/dt Na_i in component Sodium_Concentrations (millimolar).
C RATES(39) is d/dt Csqn_b in component Calcium_Concentrations (millimolar).
C RATES(13) is d/dt Ca_j in component Calcium_Concentrations (millimolar).
C RATES(14) is d/dt Ca_sl in component Calcium_Concentrations (millimolar).
C RATES(24) is d/dt Ca_i in component Calcium_Concentrations (millimolar).
C RATES(20) is d/dt Ca_SR in component Calcium_Concentrations (millimolar).
C
SUBROUTINE initConsts(CONSTS, RATES, STATES)
REAL CONSTS(*), RATES(*), STATES(*)
STATES(1) = -8.410546e1
CONSTS(1) = 8314
CONSTS(2) = 310
CONSTS(3) = 96485
CONSTS(4) = 1.381e-10
CONSTS(5) = 0
CONSTS(6) = 100000
CONSTS(7) = 40
CONSTS(8) = 1000
CONSTS(9) = 1
CONSTS(10) = 0.11
CONSTS(11) = 0.9
CONSTS(12) = 100e-5
CONSTS(13) = 10.25e-5
CONSTS(14) = 8.2413e-13
CONSTS(15) = 3.7243e-12
CONSTS(16) = 1.8313e-14
CONSTS(17) = 1.6386e-12
CONSTS(18) = 18.86
STATES(2) = 2.176608e-3
STATES(3) = 7.126555e-1
STATES(4) = 7.119893e-1
CONSTS(19) = 0.597e-3
CONSTS(20) = 140
CONSTS(21) = 0.99
CONSTS(22) = 5.4
STATES(5) = 1.00211e1
STATES(6) = 1.001974e1
CONSTS(23) = 1.5
CONSTS(24) = 11
STATES(7) = 1.896559e-2
CONSTS(25) = 0.035
STATES(8) = 3.55636e-3
CONSTS(26) = 0.0035
CONSTS(27) = 0.0035
CONSTS(28) = 0.002
STATES(9) = 3.592503e-4
STATES(10) = 8.161309e-1
STATES(11) = 3.592405e-4
STATES(12) = 9.999976e-1
CONSTS(29) = 3.7596e-2
CONSTS(30) = 1.404e-3
CONSTS(31) = 0.57153
STATES(13) = 2.048633e-4
STATES(14) = 1.188246e-4
CONSTS(32) = 0.054813
CONSTS(33) = 100e-3
CONSTS(34) = 9e-3
STATES(15) = 1.879996e-6
STATES(16) = 9.789409e-1
STATES(17) = 9.993986e-1
STATES(18) = 2.861794e-2
STATES(19) = 1.69833e-2
CONSTS(35) = 138
CONSTS(36) = 1.8
CONSTS(37) = 1.9887e-4
CONSTS(38) = 3.0375e-9
CONSTS(39) = 5.4675e-8
CONSTS(40) = 0.15e-3
CONSTS(41) = 3.59e-3
CONSTS(42) = 1.3
CONSTS(43) = 12.29
CONSTS(44) = 87.5
CONSTS(45) = 0.32
CONSTS(46) = 0.27
CONSTS(47) = 4.5
CONSTS(48) = 0.5e-3
CONSTS(49) = 0.0673
CONSTS(50) = 5.513e-4
CONSTS(51) = 15
CONSTS(52) = 1
STATES(20) = 6.138856e-1
CONSTS(53) = 0.45
CONSTS(54) = 10
CONSTS(55) = 0.5
CONSTS(56) = 0.005
CONSTS(57) = 0.06
CONSTS(58) = 25
CONSTS(59) = 0.246e-3
CONSTS(60) = 1.7
STATES(21) = 8.888214e-1
STATES(22) = 1.149876e-6
STATES(23) = 1.43831e-7
STATES(24) = 9.719632e-5
CONSTS(61) = 5.3114e-3
CONSTS(62) = 1.787
CONSTS(63) = 1
STATES(25) = 9.811535e-3
STATES(26) = 1.225802e-1
STATES(27) = 8.128604e-3
STATES(28) = 3.288063e-4
STATES(29) = 2.522168e-3
STATES(30) = 1.369514e-1
STATES(31) = 2.38683e-3
CONSTS(64) = 19.6e-3
CONSTS(65) = 0.032e-3
CONSTS(66) = 3.33e-3
CONSTS(67) = 238e-3
CONSTS(68) = 0.46e-3
CONSTS(69) = 0.057e-3
CONSTS(70) = 60e-3
CONSTS(71) = 32.7
CONSTS(72) = 2.37
CONSTS(73) = 3e-3
CONSTS(74) = 34
CONSTS(75) = 13.8
CONSTS(76) = 0.0157
CONSTS(77) = 100
CONSTS(78) = 70e-3
CONSTS(79) = 140e-3
CONSTS(80) = 24e-3
CONSTS(81) = 140e-3
CONSTS(82) = 17.1e-3
CONSTS(83) = 1.3
CONSTS(84) = 30e-3
CONSTS(85) = 100
CONSTS(86) = 100
STATES(32) = 8.606485e-3
STATES(33) = 1.101044e-2
STATES(34) = 8.078504e-2
STATES(35) = 1.238366e-1
CONSTS(87) = 7.561
CONSTS(88) = 1.65
CONSTS(89) = 1e-3
CONSTS(90) = 0.1e-3
STATES(36) = 1.001989e1
STATES(37) = 3.785209e0
STATES(38) = 8.259271e-1
STATES(39) = 1.262853
CONSTS(91) = 65
CONSTS(92) = 100
CONSTS(93) = 15
CONSTS(94) = 150
CONSTS(95) = 0.01833
CONSTS(96) = CONSTS(3)/( CONSTS(1)*CONSTS(2))
CONSTS(97) = 1.00000 - CONSTS(10)
CONSTS(98) = 1.00000 - CONSTS(11)
CONSTS(99) = 3.14159265358979*CONSTS(13) ** 2.00000*CONSTS(12)
CONSTS(100) = (EXP(CONSTS(20)/67.3000) - 1.00000)/7.00000
CONSTS(101) = 0.650000*CONSTS(99)
CONSTS(102) = (1.00000/CONSTS(96))*log(CONSTS(22)/CONSTS(35))
CONSTS(103) = (1.00000/CONSTS(96))*log(CONSTS(93)/CONSTS(94))
CONSTS(104) = 0.0200000*CONSTS(99)
CONSTS(105) = 0.0350000*CONSTS(99)
CONSTS(106) = 0.000539000*CONSTS(99)
CONSTS(107) = ( 0.000460000*CONSTS(101))/CONSTS(106)
CONSTS(108) = ( 0.0374000*CONSTS(101))/CONSTS(104)
CONSTS(109) = ( 0.000165000*CONSTS(101))/CONSTS(106)
CONSTS(110) = ( 0.0134000*CONSTS(101))/CONSTS(104)
CONSTS(111) = ( 0.140000*CONSTS(101))/CONSTS(105)
RETURN
END
SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
RATES(18) = 1.70000*STATES(13)*(1.00000 - STATES(18)) - 0.0119000*STATES(18)
RATES(19) = 1.70000*STATES(14)*(1.00000 - STATES(19)) - 0.0119000*STATES(19)
RATES(36) = (CONSTS(17)/CONSTS(101))*(STATES(6) - STATES(36))
ALGBRC(1) = 1.00000/1.00000+EXP(- (56.8600+STATES(1))/9.03000) ** 2.00000
ALGBRC(13) = 0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
RATES(2) = (ALGBRC(1) - STATES(2))/ALGBRC(13)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0000)/5.00000))
ALGBRC(16) = 3300.00/( (1.00000+EXP((- 22.0000 - STATES(1))/9.00000))*(1.00000+EXP((STATES(1)+11.0000)/9.00000)))+230.000/(1.00000+EXP((STATES(1)+40.0000)/20.0000))
RATES(7) = (ALGBRC(4) - STATES(7))/ALGBRC(16)
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+3.80000)/14.2500))
ALGBRC(17) = 990.100/(1.00000+EXP(- (STATES(1)+2.43600)/14.1200))
RATES(8) = (ALGBRC(5) - STATES(8))/ALGBRC(17)
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.0000)/13.0000))
ALGBRC(18) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
RATES(9) = (ALGBRC(6) - STATES(9))/ALGBRC(18)
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+19.5000)/5.00000))
ALGBRC(19) = 800.000/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
RATES(10) = (ALGBRC(7) - STATES(10))/ALGBRC(19)
ALGBRC(20) = 8.50000*EXP(- (STATES(1)+45.0000)/50.0000 ** 2.00000)+0.500000
RATES(11) = (ALGBRC(6) - STATES(11))/ALGBRC(20)
ALGBRC(21) = 85.0000*EXP(- STATES(1)+40.0000 ** 2.00000/220.000)+7.00000
RATES(12) = (ALGBRC(7) - STATES(12))/ALGBRC(21)
ALGBRC(12) = CONSTS(51) - (CONSTS(51) - CONSTS(52))/(1.00000+CONSTS(53)/STATES(20) ** 2.50000)
ALGBRC(26) = CONSTS(54)/ALGBRC(12)
ALGBRC(33) = CONSTS(55)*ALGBRC(12)
RATES(22) = ( ALGBRC(26)*STATES(13) ** 2.00000*STATES(21) - CONSTS(57)*STATES(22)) - ( ALGBRC(33)*STATES(13)*STATES(22) - CONSTS(56)*STATES(23))
ALGBRC(2) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, 0.0570000*EXP(- (STATES(1)+80.0000)/6.80000))
ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 5.92310/(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 2.70000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.348500*STATES(1)))
ALGBRC(27) = 1.00000/(ALGBRC(2)+ALGBRC(14))
ALGBRC(34) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
RATES(3) = (ALGBRC(34) - STATES(3))/ALGBRC(27)
ALGBRC(3) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 25428.0*EXP( 0.244400*STATES(1)) - 6.94800e-06*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 0.600000*EXP( 0.0570000*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.0242400*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
ALGBRC(28) = 1.00000/(ALGBRC(3)+ALGBRC(15))
ALGBRC(35) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
RATES(4) = (ALGBRC(35) - STATES(4))/ALGBRC(28)
ALGBRC(40) = ((1.00000 - STATES(21)) - STATES(22)) - STATES(23)
RATES(21) = ( CONSTS(56)*ALGBRC(40) - ALGBRC(33)*STATES(13)*STATES(21)) - ( ALGBRC(26)*STATES(13) ** 2.00000*STATES(21) - CONSTS(57)*STATES(22))
RATES(23) = ( ALGBRC(33)*STATES(13)*STATES(22) - CONSTS(56)*STATES(23)) - ( CONSTS(57)*STATES(23) - ALGBRC(26)*STATES(13) ** 2.00000*ALGBRC(40))
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+5.00000)/6.00000))
ALGBRC(22) = 1.40000/(1.00000+EXP((- 35.0000 - STATES(1))/13.0000))+0.250000
ALGBRC(29) = 1.40000/(1.00000+EXP((STATES(1)+5.00000)/5.00000))
ALGBRC(36) = 1.00000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(41) = ALGBRC(22)*ALGBRC(29)+ALGBRC(36)
RATES(15) = (ALGBRC(8) - STATES(15))/ALGBRC(41)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+20.0000)/7.00000))
ALGBRC(23) = 1102.50*EXP(- (STATES(1)+27.0000)/15.0000 ** 2.00000)
ALGBRC(30) = 200.000/(1.00000+EXP((13.0000 - STATES(1))/10.0000))
ALGBRC(37) = 180.000/(1.00000+EXP((STATES(1)+30.0000)/10.0000))+20.0000
ALGBRC(42) = ALGBRC(23)+ALGBRC(30)+ALGBRC(37)
RATES(16) = (ALGBRC(9) - STATES(16))/ALGBRC(42)
ALGBRC(10) = 0.670000/(1.00000+EXP((STATES(1)+35.0000)/7.00000))+0.330000
ALGBRC(24) = 300.000*EXP(- STATES(1)+25.0000 ** 2.00000/170.000)
ALGBRC(31) = 31.0000/(1.00000+EXP((25.0000 - STATES(1))/10.0000))
ALGBRC(38) = 16.0000/(1.00000+EXP((STATES(1)+30.0000)/10.0000))
ALGBRC(43) = ALGBRC(24)+ALGBRC(31)+ALGBRC(38)
RATES(17) = (ALGBRC(10) - STATES(17))/ALGBRC(43)
RATES(25) = CONSTS(71)*STATES(24)*(CONSTS(78) - STATES(25)) - CONSTS(64)*STATES(25)
RATES(32) = CONSTS(85)*STATES(13)*(CONSTS(107) - STATES(32)) - CONSTS(83)*STATES(32)
RATES(26) = CONSTS(72)*STATES(24)*((CONSTS(79) - STATES(26)) - STATES(27)) - CONSTS(65)*STATES(26)
RATES(39) = CONSTS(92)*STATES(20)*(CONSTS(111) - STATES(39)) - CONSTS(91)*STATES(39)
RATES(33) = CONSTS(85)*STATES(14)*(CONSTS(108) - STATES(33)) - CONSTS(83)*STATES(33)
RATES(27) = CONSTS(73)*CONSTS(63)*((CONSTS(79) - STATES(26)) - STATES(27)) - CONSTS(66)*STATES(27)
ALGBRC(87) = CONSTS(58)*STATES(22)*(STATES(20) - STATES(13))
ALGBRC(90) = 5.34800e-06*(STATES(20) - STATES(13))
ALGBRC(89) = ( CONSTS(61)*(STATES(24)/CONSTS(59) ** CONSTS(62) - STATES(20)/CONSTS(60) ** CONSTS(62)))/(1.00000+STATES(24)/CONSTS(59) ** CONSTS(62)+STATES(20)/CONSTS(60) ** CONSTS(62))
RATES(20) = (ALGBRC(89) - (( ALGBRC(90)*CONSTS(101))/CONSTS(105)+ALGBRC(87))) - RATES(39)
RATES(34) = CONSTS(86)*STATES(13)*(CONSTS(109) - STATES(34)) - CONSTS(84)*STATES(34)
RATES(28) = CONSTS(74)*STATES(24)*(CONSTS(80) - STATES(28)) - CONSTS(67)*STATES(28)
RATES(35) = CONSTS(86)*STATES(14)*(CONSTS(110) - STATES(35)) - CONSTS(84)*STATES(35)
RATES(29) = CONSTS(75)*STATES(24)*((CONSTS(81) - STATES(29)) - STATES(30)) - CONSTS(68)*STATES(29)
RATES(30) = CONSTS(76)*CONSTS(63)*((CONSTS(81) - STATES(29)) - STATES(30)) - CONSTS(69)*STATES(30)
RATES(31) = CONSTS(77)*STATES(24)*(CONSTS(82) - STATES(31)) - CONSTS(70)*STATES(31)
RATES(37) = CONSTS(90)*STATES(5)*(CONSTS(87) - STATES(37)) - CONSTS(89)*STATES(37)
RATES(38) = CONSTS(90)*STATES(6)*(CONSTS(88) - STATES(38)) - CONSTS(89)*STATES(38)
ALGBRC(93) = RATES(25)+RATES(26)+RATES(27)+RATES(28)+RATES(29)+RATES(30)+RATES(31)
RATES(24) = (( - ALGBRC(89)*CONSTS(105))/CONSTS(101) - ALGBRC(93))+ (CONSTS(15)/CONSTS(101))*(STATES(14) - STATES(24))
ALGBRC(94) = (1.00000/CONSTS(96))*log(CONSTS(20)/STATES(5))
ALGBRC(95) = CONSTS(10)*CONSTS(18)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(94))
ALGBRC(96) = CONSTS(10)*CONSTS(19)*(STATES(1) - ALGBRC(94))
ALGBRC(25) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(96))+ 0.0365000*CONSTS(100)*EXP( - STATES(1)*CONSTS(96)))
ALGBRC(32) = ( CONSTS(10)*CONSTS(21)*ALGBRC(25)*CONSTS(22))/( (1.00000+CONSTS(24)/STATES(5) ** 4.00000)*(CONSTS(22)+CONSTS(23)))
ALGBRC(64) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(5)*EXP( STATES(1)*CONSTS(96)) - CONSTS(20)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(70) = CONSTS(11)*ALGBRC(64)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18))
ALGBRC(72) = 1.00000/(1.00000+CONSTS(40)/STATES(13) ** 2.00000)
ALGBRC(76) = EXP( CONSTS(46)*STATES(1)*CONSTS(96))*STATES(5) ** 3.00000*CONSTS(36)
ALGBRC(78) = EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))*CONSTS(20) ** 3.00000*STATES(13)
ALGBRC(80) = CONSTS(41)*CONSTS(20) ** 3.00000*(1.00000+STATES(5)/CONSTS(43) ** 3.00000)+ CONSTS(44) ** 3.00000*STATES(13)*(1.00000+STATES(13)/CONSTS(41))+ CONSTS(42)*STATES(5) ** 3.00000+ STATES(5) ** 3.00000*CONSTS(36)+ CONSTS(20) ** 3.00000*STATES(13)
ALGBRC(82) = ( CONSTS(10)*CONSTS(47)*ALGBRC(72)*(ALGBRC(76) - ALGBRC(78)))/( ALGBRC(80)*(1.00000+ CONSTS(45)*EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))))
ALGBRC(97) = ALGBRC(95)+ALGBRC(96)+ 3.00000*ALGBRC(82)+ 3.00000*ALGBRC(32)+ALGBRC(70)
RATES(5) = (( - ALGBRC(97)*CONSTS(4))/( CONSTS(106)*CONSTS(3))+ (CONSTS(16)/CONSTS(106))*(STATES(6) - STATES(5))) - RATES(37)
ALGBRC(98) = (1.00000/CONSTS(96))*log(CONSTS(20)/STATES(6))
ALGBRC(99) = CONSTS(97)*CONSTS(18)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(98))
ALGBRC(101) = CONSTS(97)*CONSTS(19)*(STATES(1) - ALGBRC(98))
ALGBRC(39) = ( CONSTS(97)*CONSTS(21)*ALGBRC(25)*CONSTS(22))/( (1.00000+CONSTS(24)/STATES(6) ** 4.00000)*(CONSTS(22)+CONSTS(23)))
ALGBRC(65) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(6)*EXP( STATES(1)*CONSTS(96)) - CONSTS(20)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(71) = CONSTS(98)*ALGBRC(65)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19))
ALGBRC(74) = 1.00000/(1.00000+CONSTS(40)/STATES(14) ** 2.00000)
ALGBRC(77) = EXP( CONSTS(46)*STATES(1)*CONSTS(96))*STATES(6) ** 3.00000*CONSTS(36)
ALGBRC(79) = EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))*CONSTS(20) ** 3.00000*STATES(14)
ALGBRC(81) = CONSTS(41)*CONSTS(20) ** 3.00000*(1.00000+STATES(6)/CONSTS(43) ** 3.00000)+ CONSTS(44) ** 3.00000*STATES(14)*(1.00000+STATES(14)/CONSTS(41))+ CONSTS(42)*STATES(6) ** 3.00000+ STATES(6) ** 3.00000*CONSTS(36)+ CONSTS(20) ** 3.00000*STATES(14)
ALGBRC(83) = ( CONSTS(97)*CONSTS(47)*ALGBRC(74)*(ALGBRC(77) - ALGBRC(79)))/( ALGBRC(81)*(1.00000+ CONSTS(45)*EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))))
ALGBRC(103) = ALGBRC(99)+ALGBRC(101)+ 3.00000*ALGBRC(83)+ 3.00000*ALGBRC(39)+ALGBRC(71)
RATES(6) = (( - ALGBRC(103)*CONSTS(4))/( CONSTS(104)*CONSTS(3))+ (CONSTS(16)/CONSTS(104))*(STATES(5) - STATES(6))+ (CONSTS(17)/CONSTS(104))*(STATES(36) - STATES(6))) - RATES(38)
ALGBRC(61) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(13)*EXP( 2.00000*STATES(1)*CONSTS(96)) - CONSTS(36)))/(EXP( 2.00000*STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(66) = CONSTS(11)*ALGBRC(61)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18))
ALGBRC(84) = ( CONSTS(10)*CONSTS(49)*STATES(13) ** 1.60000)/(CONSTS(48) ** 1.60000+STATES(13) ** 1.60000)
ALGBRC(110) = (1.00000/( 2.00000*CONSTS(96)))*log(CONSTS(36)/STATES(13))
ALGBRC(111) = CONSTS(10)*CONSTS(50)*(STATES(1) - ALGBRC(110))
ALGBRC(112) = (ALGBRC(66)+ALGBRC(111)+ALGBRC(84)) - 2.00000*ALGBRC(82)
ALGBRC(91) = RATES(32)+RATES(34)
RATES(13) = ((( - ALGBRC(112)*CONSTS(4))/( CONSTS(106)*2.00000*CONSTS(3))+ (CONSTS(14)/CONSTS(106))*(STATES(14) - STATES(13))) - ALGBRC(91))+( ALGBRC(87)*CONSTS(105))/CONSTS(106)+( ALGBRC(90)*CONSTS(101))/CONSTS(106)
ALGBRC(62) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(14)*EXP( 2.00000*STATES(1)*CONSTS(96)) - CONSTS(36)))/(EXP( 2.00000*STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(67) = CONSTS(98)*ALGBRC(62)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19))
ALGBRC(86) = ( CONSTS(97)*CONSTS(49)*STATES(14) ** 1.60000)/(CONSTS(48) ** 1.60000+STATES(14) ** 1.60000)
ALGBRC(113) = (1.00000/( 2.00000*CONSTS(96)))*log(CONSTS(36)/STATES(14))
ALGBRC(114) = CONSTS(97)*CONSTS(50)*(STATES(1) - ALGBRC(113))
ALGBRC(115) = (ALGBRC(67)+ALGBRC(114)+ALGBRC(86)) - 2.00000*ALGBRC(83)
ALGBRC(92) = RATES(33)+RATES(35)
RATES(14) = (( - ALGBRC(115)*CONSTS(4))/( CONSTS(104)*2.00000*CONSTS(3))+ (CONSTS(14)/CONSTS(104))*(STATES(13) - STATES(14))+ (CONSTS(15)/CONSTS(104))*(STATES(24) - STATES(14))) - ALGBRC(92)
ALGBRC(44) = ALGBRC(32)+ALGBRC(39)
ALGBRC(45) = 1.00000/(1.00000+EXP((STATES(1)+74.0000)/24.0000))
ALGBRC(46) = CONSTS(25)* (CONSTS(22)/5.40000) ** (1.0 / 2)*STATES(7)*ALGBRC(45)*(STATES(1) - CONSTS(102))
ALGBRC(105) = (1.00000/CONSTS(96))*log((CONSTS(22)+ CONSTS(95)*CONSTS(20))/(CONSTS(35)+ CONSTS(95)*STATES(36)))
ALGBRC(106) = CONSTS(10)*CONSTS(26)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105))
ALGBRC(107) = CONSTS(97)*CONSTS(27)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105))
ALGBRC(108) = ALGBRC(106)+ALGBRC(107)
ALGBRC(47) = 1.00000/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(48) = CONSTS(28)*ALGBRC(47)*(STATES(1) - CONSTS(102))
ALGBRC(49) = CONSTS(29)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(102))
ALGBRC(50) = CONSTS(30)*STATES(11)*STATES(12)*(STATES(1) - CONSTS(102))
ALGBRC(51) = ALGBRC(49)+ALGBRC(50)
ALGBRC(52) = 4.09380/(1.00000+EXP( 0.121650*((STATES(1) - CONSTS(102)) - 49.9344)))
ALGBRC(53) = ( 15.7197*EXP( 0.0673900*((STATES(1) - CONSTS(102)) - 3.25710))+EXP( 0.0617500*((STATES(1) - CONSTS(102)) - 594.310)))/(1.00000+EXP( - 0.162850*((STATES(1) - CONSTS(102))+14.2067)))
ALGBRC(54) = ALGBRC(52)/(ALGBRC(52)+ALGBRC(53))
ALGBRC(55) = CONSTS(31)* (CONSTS(22)/5.40000) ** (1.0 / 2)*ALGBRC(54)*(STATES(1) - CONSTS(102))
ALGBRC(63) = ( CONSTS(39)*STATES(1)*CONSTS(3)*CONSTS(96)*( CONSTS(35)*EXP( STATES(1)*CONSTS(96)) - CONSTS(22)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(69) = ALGBRC(63)*STATES(15)*STATES(16)*STATES(17)*( CONSTS(11)*(1.00000 - STATES(18))+ CONSTS(98)*(1.00000 - STATES(19)))
ALGBRC(109) = ((ALGBRC(51)+ALGBRC(46)+ALGBRC(108)+ALGBRC(55)) - 2.00000*ALGBRC(44))+ALGBRC(69)+ALGBRC(48)
ALGBRC(104) = ALGBRC(97)+ALGBRC(103)
ALGBRC(56) = (( CONSTS(10)*CONSTS(32))/(1.00000+CONSTS(33)/STATES(13)))*(STATES(1) - CONSTS(103))
ALGBRC(57) = (( CONSTS(97)*CONSTS(32))/(1.00000+CONSTS(33)/STATES(14)))*(STATES(1) - CONSTS(103))
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(59) = CONSTS(34)*(STATES(1) - CONSTS(103))
ALGBRC(60) = ALGBRC(58)+ALGBRC(59)
ALGBRC(117) = ALGBRC(112)+ALGBRC(115)
ALGBRC(118) = ALGBRC(104)+ALGBRC(60)+ALGBRC(117)+ALGBRC(109)
ALGBRC(11) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(8))*CONSTS(8).LE.CONSTS(9), CONSTS(7), 0.00000)
RATES(1) = - (ALGBRC(118) - ALGBRC(11))
RETURN
END
SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
ALGBRC(1) = 1.00000/1.00000+EXP(- (56.8600+STATES(1))/9.03000) ** 2.00000
ALGBRC(13) = 0.129200*EXP(- (STATES(1)+45.7900)/15.5400 ** 2.00000)+ 0.0648700*EXP(- (STATES(1) - 4.82300)/51.1200 ** 2.00000)
ALGBRC(4) = 1.00000/(1.00000+EXP(- (STATES(1)+10.0000)/5.00000))
ALGBRC(16) = 3300.00/( (1.00000+EXP((- 22.0000 - STATES(1))/9.00000))*(1.00000+EXP((STATES(1)+11.0000)/9.00000)))+230.000/(1.00000+EXP((STATES(1)+40.0000)/20.0000))
ALGBRC(5) = 1.00000/(1.00000+EXP(- (STATES(1)+3.80000)/14.2500))
ALGBRC(17) = 990.100/(1.00000+EXP(- (STATES(1)+2.43600)/14.1200))
ALGBRC(6) = 1.00000/(1.00000+EXP(- (STATES(1) - 19.0000)/13.0000))
ALGBRC(18) = 9.00000/(1.00000+EXP((STATES(1)+3.00000)/15.0000))+0.500000
ALGBRC(7) = 1.00000/(1.00000+EXP((STATES(1)+19.5000)/5.00000))
ALGBRC(19) = 800.000/(1.00000+EXP((STATES(1)+60.0000)/10.0000))+30.0000
ALGBRC(20) = 8.50000*EXP(- (STATES(1)+45.0000)/50.0000 ** 2.00000)+0.500000
ALGBRC(21) = 85.0000*EXP(- STATES(1)+40.0000 ** 2.00000/220.000)+7.00000
ALGBRC(12) = CONSTS(51) - (CONSTS(51) - CONSTS(52))/(1.00000+CONSTS(53)/STATES(20) ** 2.50000)
ALGBRC(26) = CONSTS(54)/ALGBRC(12)
ALGBRC(33) = CONSTS(55)*ALGBRC(12)
ALGBRC(2) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, 0.0570000*EXP(- (STATES(1)+80.0000)/6.80000))
ALGBRC(14) = TERNRY(STATES(1).GE.- 40.0000, 5.92310/(1.00000+EXP(- (STATES(1)+10.6600)/11.1000)), 2.70000*EXP( 0.0790000*STATES(1))+ 310000.*EXP( 0.348500*STATES(1)))
ALGBRC(27) = 1.00000/(ALGBRC(2)+ALGBRC(14))
ALGBRC(34) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
ALGBRC(3) = TERNRY(STATES(1).GE.- 40.0000, 0.00000, ( ( - 25428.0*EXP( 0.244400*STATES(1)) - 6.94800e-06*EXP( - 0.0439100*STATES(1)))*(STATES(1)+37.7800))/(1.00000+EXP( 0.311000*(STATES(1)+79.2300))))
ALGBRC(15) = TERNRY(STATES(1).GE.- 40.0000, ( 0.600000*EXP( 0.0570000*STATES(1)))/(1.00000+EXP( - 0.100000*(STATES(1)+32.0000))), ( 0.0242400*EXP( - 0.0105200*STATES(1)))/(1.00000+EXP( - 0.137800*(STATES(1)+40.1400))))
ALGBRC(28) = 1.00000/(ALGBRC(3)+ALGBRC(15))
ALGBRC(35) = 1.00000/1.00000+EXP((STATES(1)+71.5500)/7.43000) ** 2.00000
ALGBRC(40) = ((1.00000 - STATES(21)) - STATES(22)) - STATES(23)
ALGBRC(8) = 1.00000/(1.00000+EXP(- (STATES(1)+5.00000)/6.00000))
ALGBRC(22) = 1.40000/(1.00000+EXP((- 35.0000 - STATES(1))/13.0000))+0.250000
ALGBRC(29) = 1.40000/(1.00000+EXP((STATES(1)+5.00000)/5.00000))
ALGBRC(36) = 1.00000/(1.00000+EXP((50.0000 - STATES(1))/20.0000))
ALGBRC(41) = ALGBRC(22)*ALGBRC(29)+ALGBRC(36)
ALGBRC(9) = 1.00000/(1.00000+EXP((STATES(1)+20.0000)/7.00000))
ALGBRC(23) = 1102.50*EXP(- (STATES(1)+27.0000)/15.0000 ** 2.00000)
ALGBRC(30) = 200.000/(1.00000+EXP((13.0000 - STATES(1))/10.0000))
ALGBRC(37) = 180.000/(1.00000+EXP((STATES(1)+30.0000)/10.0000))+20.0000
ALGBRC(42) = ALGBRC(23)+ALGBRC(30)+ALGBRC(37)
ALGBRC(10) = 0.670000/(1.00000+EXP((STATES(1)+35.0000)/7.00000))+0.330000
ALGBRC(24) = 300.000*EXP(- STATES(1)+25.0000 ** 2.00000/170.000)
ALGBRC(31) = 31.0000/(1.00000+EXP((25.0000 - STATES(1))/10.0000))
ALGBRC(38) = 16.0000/(1.00000+EXP((STATES(1)+30.0000)/10.0000))
ALGBRC(43) = ALGBRC(24)+ALGBRC(31)+ALGBRC(38)
ALGBRC(87) = CONSTS(58)*STATES(22)*(STATES(20) - STATES(13))
ALGBRC(90) = 5.34800e-06*(STATES(20) - STATES(13))
ALGBRC(89) = ( CONSTS(61)*(STATES(24)/CONSTS(59) ** CONSTS(62) - STATES(20)/CONSTS(60) ** CONSTS(62)))/(1.00000+STATES(24)/CONSTS(59) ** CONSTS(62)+STATES(20)/CONSTS(60) ** CONSTS(62))
ALGBRC(93) = RATES(25)+RATES(26)+RATES(27)+RATES(28)+RATES(29)+RATES(30)+RATES(31)
ALGBRC(94) = (1.00000/CONSTS(96))*log(CONSTS(20)/STATES(5))
ALGBRC(95) = CONSTS(10)*CONSTS(18)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(94))
ALGBRC(96) = CONSTS(10)*CONSTS(19)*(STATES(1) - ALGBRC(94))
ALGBRC(25) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(96))+ 0.0365000*CONSTS(100)*EXP( - STATES(1)*CONSTS(96)))
ALGBRC(32) = ( CONSTS(10)*CONSTS(21)*ALGBRC(25)*CONSTS(22))/( (1.00000+CONSTS(24)/STATES(5) ** 4.00000)*(CONSTS(22)+CONSTS(23)))
ALGBRC(64) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(5)*EXP( STATES(1)*CONSTS(96)) - CONSTS(20)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(70) = CONSTS(11)*ALGBRC(64)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18))
ALGBRC(72) = 1.00000/(1.00000+CONSTS(40)/STATES(13) ** 2.00000)
ALGBRC(76) = EXP( CONSTS(46)*STATES(1)*CONSTS(96))*STATES(5) ** 3.00000*CONSTS(36)
ALGBRC(78) = EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))*CONSTS(20) ** 3.00000*STATES(13)
ALGBRC(80) = CONSTS(41)*CONSTS(20) ** 3.00000*(1.00000+STATES(5)/CONSTS(43) ** 3.00000)+ CONSTS(44) ** 3.00000*STATES(13)*(1.00000+STATES(13)/CONSTS(41))+ CONSTS(42)*STATES(5) ** 3.00000+ STATES(5) ** 3.00000*CONSTS(36)+ CONSTS(20) ** 3.00000*STATES(13)
ALGBRC(82) = ( CONSTS(10)*CONSTS(47)*ALGBRC(72)*(ALGBRC(76) - ALGBRC(78)))/( ALGBRC(80)*(1.00000+ CONSTS(45)*EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))))
ALGBRC(97) = ALGBRC(95)+ALGBRC(96)+ 3.00000*ALGBRC(82)+ 3.00000*ALGBRC(32)+ALGBRC(70)
ALGBRC(98) = (1.00000/CONSTS(96))*log(CONSTS(20)/STATES(6))
ALGBRC(99) = CONSTS(97)*CONSTS(18)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(98))
ALGBRC(101) = CONSTS(97)*CONSTS(19)*(STATES(1) - ALGBRC(98))
ALGBRC(39) = ( CONSTS(97)*CONSTS(21)*ALGBRC(25)*CONSTS(22))/( (1.00000+CONSTS(24)/STATES(6) ** 4.00000)*(CONSTS(22)+CONSTS(23)))
ALGBRC(65) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(6)*EXP( STATES(1)*CONSTS(96)) - CONSTS(20)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(71) = CONSTS(98)*ALGBRC(65)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19))
ALGBRC(74) = 1.00000/(1.00000+CONSTS(40)/STATES(14) ** 2.00000)
ALGBRC(77) = EXP( CONSTS(46)*STATES(1)*CONSTS(96))*STATES(6) ** 3.00000*CONSTS(36)
ALGBRC(79) = EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))*CONSTS(20) ** 3.00000*STATES(14)
ALGBRC(81) = CONSTS(41)*CONSTS(20) ** 3.00000*(1.00000+STATES(6)/CONSTS(43) ** 3.00000)+ CONSTS(44) ** 3.00000*STATES(14)*(1.00000+STATES(14)/CONSTS(41))+ CONSTS(42)*STATES(6) ** 3.00000+ STATES(6) ** 3.00000*CONSTS(36)+ CONSTS(20) ** 3.00000*STATES(14)
ALGBRC(83) = ( CONSTS(97)*CONSTS(47)*ALGBRC(74)*(ALGBRC(77) - ALGBRC(79)))/( ALGBRC(81)*(1.00000+ CONSTS(45)*EXP( (CONSTS(46) - 1.00000)*STATES(1)*CONSTS(96))))
ALGBRC(103) = ALGBRC(99)+ALGBRC(101)+ 3.00000*ALGBRC(83)+ 3.00000*ALGBRC(39)+ALGBRC(71)
ALGBRC(61) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(13)*EXP( 2.00000*STATES(1)*CONSTS(96)) - CONSTS(36)))/(EXP( 2.00000*STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(66) = CONSTS(11)*ALGBRC(61)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18))
ALGBRC(84) = ( CONSTS(10)*CONSTS(49)*STATES(13) ** 1.60000)/(CONSTS(48) ** 1.60000+STATES(13) ** 1.60000)
ALGBRC(110) = (1.00000/( 2.00000*CONSTS(96)))*log(CONSTS(36)/STATES(13))
ALGBRC(111) = CONSTS(10)*CONSTS(50)*(STATES(1) - ALGBRC(110))
ALGBRC(112) = (ALGBRC(66)+ALGBRC(111)+ALGBRC(84)) - 2.00000*ALGBRC(82)
ALGBRC(91) = RATES(32)+RATES(34)
ALGBRC(62) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(96)*( STATES(14)*EXP( 2.00000*STATES(1)*CONSTS(96)) - CONSTS(36)))/(EXP( 2.00000*STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(67) = CONSTS(98)*ALGBRC(62)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19))
ALGBRC(86) = ( CONSTS(97)*CONSTS(49)*STATES(14) ** 1.60000)/(CONSTS(48) ** 1.60000+STATES(14) ** 1.60000)
ALGBRC(113) = (1.00000/( 2.00000*CONSTS(96)))*log(CONSTS(36)/STATES(14))
ALGBRC(114) = CONSTS(97)*CONSTS(50)*(STATES(1) - ALGBRC(113))
ALGBRC(115) = (ALGBRC(67)+ALGBRC(114)+ALGBRC(86)) - 2.00000*ALGBRC(83)
ALGBRC(92) = RATES(33)+RATES(35)
ALGBRC(44) = ALGBRC(32)+ALGBRC(39)
ALGBRC(45) = 1.00000/(1.00000+EXP((STATES(1)+74.0000)/24.0000))
ALGBRC(46) = CONSTS(25)* (CONSTS(22)/5.40000) ** (1.0 / 2)*STATES(7)*ALGBRC(45)*(STATES(1) - CONSTS(102))
ALGBRC(105) = (1.00000/CONSTS(96))*log((CONSTS(22)+ CONSTS(95)*CONSTS(20))/(CONSTS(35)+ CONSTS(95)*STATES(36)))
ALGBRC(106) = CONSTS(10)*CONSTS(26)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105))
ALGBRC(107) = CONSTS(97)*CONSTS(27)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105))
ALGBRC(108) = ALGBRC(106)+ALGBRC(107)
ALGBRC(47) = 1.00000/(1.00000+EXP(7.48800 - STATES(1)/5.98000))
ALGBRC(48) = CONSTS(28)*ALGBRC(47)*(STATES(1) - CONSTS(102))
ALGBRC(49) = CONSTS(29)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(102))
ALGBRC(50) = CONSTS(30)*STATES(11)*STATES(12)*(STATES(1) - CONSTS(102))
ALGBRC(51) = ALGBRC(49)+ALGBRC(50)
ALGBRC(52) = 4.09380/(1.00000+EXP( 0.121650*((STATES(1) - CONSTS(102)) - 49.9344)))
ALGBRC(53) = ( 15.7197*EXP( 0.0673900*((STATES(1) - CONSTS(102)) - 3.25710))+EXP( 0.0617500*((STATES(1) - CONSTS(102)) - 594.310)))/(1.00000+EXP( - 0.162850*((STATES(1) - CONSTS(102))+14.2067)))
ALGBRC(54) = ALGBRC(52)/(ALGBRC(52)+ALGBRC(53))
ALGBRC(55) = CONSTS(31)* (CONSTS(22)/5.40000) ** (1.0 / 2)*ALGBRC(54)*(STATES(1) - CONSTS(102))
ALGBRC(63) = ( CONSTS(39)*STATES(1)*CONSTS(3)*CONSTS(96)*( CONSTS(35)*EXP( STATES(1)*CONSTS(96)) - CONSTS(22)))/(EXP( STATES(1)*CONSTS(96)) - 1.00000)
ALGBRC(69) = ALGBRC(63)*STATES(15)*STATES(16)*STATES(17)*( CONSTS(11)*(1.00000 - STATES(18))+ CONSTS(98)*(1.00000 - STATES(19)))
ALGBRC(109) = ((ALGBRC(51)+ALGBRC(46)+ALGBRC(108)+ALGBRC(55)) - 2.00000*ALGBRC(44))+ALGBRC(69)+ALGBRC(48)
ALGBRC(104) = ALGBRC(97)+ALGBRC(103)
ALGBRC(56) = (( CONSTS(10)*CONSTS(32))/(1.00000+CONSTS(33)/STATES(13)))*(STATES(1) - CONSTS(103))
ALGBRC(57) = (( CONSTS(97)*CONSTS(32))/(1.00000+CONSTS(33)/STATES(14)))*(STATES(1) - CONSTS(103))
ALGBRC(58) = ALGBRC(56)+ALGBRC(57)
ALGBRC(59) = CONSTS(34)*(STATES(1) - CONSTS(103))
ALGBRC(60) = ALGBRC(58)+ALGBRC(59)
ALGBRC(117) = ALGBRC(112)+ALGBRC(115)
ALGBRC(118) = ALGBRC(104)+ALGBRC(60)+ALGBRC(117)+ALGBRC(109)
ALGBRC(11) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(8))*CONSTS(8).LE.CONSTS(9), CONSTS(7), 0.00000)
ALGBRC(68) = ALGBRC(66)+ALGBRC(67)
ALGBRC(73) = ALGBRC(70)+ALGBRC(71)
ALGBRC(75) = ALGBRC(68)+ALGBRC(69)+ALGBRC(73)
ALGBRC(85) = ALGBRC(82)+ALGBRC(83)
ALGBRC(88) = ALGBRC(84)+ALGBRC(86)
ALGBRC(100) = ALGBRC(95)+ALGBRC(99)
ALGBRC(102) = ALGBRC(96)+ALGBRC(101)
ALGBRC(116) = ALGBRC(111)+ALGBRC(114)
RETURN
END
REAL FUNCTION TERNRY(TEST, VALA, VALB)
LOGICAL TEST
REAL VALA, VALB
IF (TEST) THEN
TERNRY = VALA
ELSE
TERNRY = VALB
ENDIF
RETURN
END