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