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 110 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(95) 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_Amplitude in component membrane (picoA_per_picoF). C CONSTS(7) is I_Stim_Period in component membrane (millisecond). C CONSTS(8) is I_Stim_PulseDuration in component membrane (millisecond). C CONSTS(9) is F_junc in component membrane (dimensionless). C CONSTS(96) is F_sl in component membrane (dimensionless). C CONSTS(10) is F_junc_Ca_L in component membrane (dimensionless). C CONSTS(97) is F_sl_Ca_L in component membrane (dimensionless). C CONSTS(98) is V_cell in component membrane (litre). C CONSTS(11) is Length_cell in component membrane (decimetre). C CONSTS(12) is Radius_cell in component membrane (decimetre). C CONSTS(100) is V_myo in component membrane (litre). C CONSTS(103) is V_sl in component membrane (litre). C CONSTS(104) is V_sr in component membrane (litre). C CONSTS(105) is V_junc in component membrane (litre). C CONSTS(13) is J_Ca_juncsl in component membrane (L_per_millis). C CONSTS(14) is J_Ca_slmyo in component membrane (L_per_millis). C CONSTS(15) is J_Na_juncsl in component membrane (L_per_millis). C CONSTS(16) 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(17) 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(18) is G_Na_B in component Background_Na_Current (milliS_per_microF). C CONSTS(19) is Na_o in component Sodium_Concentrations (millimolar). C CONSTS(99) 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(20) is Ibar_NaK in component Na_K_Pump_Current (picoA_per_picoF). C CONSTS(21) 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(22) is Km_Ko in component Na_K_Pump_Current (millimolar). C CONSTS(23) 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(101) is E_K in component Nerst_Potentials (millivolt). C CONSTS(24) 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(25) is G_Ks_junc in component Slowly_Activating_K_Current (milliS_per_microF). C CONSTS(26) 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(27) 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(28) is G_to_s in component Transient_Outward_K_Current (milliS_per_microF). C CONSTS(29) 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(30) 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(102) is E_Cl in component Nerst_Potentials (millivolt). C CONSTS(31) is G_ClCa in component Ca_Activated_Cl_Current (milliS_per_microF). C CONSTS(32) is Kd_ClCa in component Ca_Activated_Cl_Current (millimolar). C CONSTS(33) 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(34) is K_i in component Potassium_Concentrations (millimolar). C CONSTS(35) 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(36) is p_Ca in component L_Type_Calcium_Current (litre_per_farad_millisecond). C CONSTS(37) is p_Na in component L_Type_Calcium_Current (litre_per_farad_millisecond). C CONSTS(38) 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(39) is Kd_act in component Na_Ca_Exchanger_Current (millimolar). C CONSTS(40) is Km_Ca_i in component Na_Ca_Exchanger_Current (millimolar). C CONSTS(41) is Km_Ca_o in component Na_Ca_Exchanger_Current (millimolar). C CONSTS(42) is Km_Na_i in component Na_Ca_Exchanger_Current (millimolar). C CONSTS(43) is Km_Na_o in component Na_Ca_Exchanger_Current (millimolar). C CONSTS(44) is k_sat in component Na_Ca_Exchanger_Current (dimensionless). C CONSTS(45) 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(46) 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(47) is Km_P_Ca in component Sarcolemmal_Ca_Pump_Current (millimolar). C CONSTS(48) 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(49) 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(50) is Max_SR in component SR_Fluxes (dimensionless). C CONSTS(51) is Min_SR in component SR_Fluxes (dimensionless). C STATES(20) is Ca_SR in component Calcium_Concentrations (millimolar). C CONSTS(52) is ec50_SR in component SR_Fluxes (millimolar). C ALGBRC(12) is k_Ca_SR in component SR_Fluxes (dimensionless). C CONSTS(53) 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(54) 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(55) is ki_m in component SR_Fluxes (per_millis). C CONSTS(56) is ko_m in component SR_Fluxes (per_millis). C CONSTS(57) is ks in component SR_Fluxes (per_millis). C CONSTS(58) is Km_f in component SR_Fluxes (millimolar). C CONSTS(59) 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(60) is V_max_SR_CaP in component SR_Fluxes (milliM_per_millis). C CONSTS(61) is hill_SR_CaP in component SR_Fluxes (dimensionless). C CONSTS(62) 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(63) is k_off_TnC_l in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(64) is k_off_TnC_h_Ca in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(65) is k_off_TnC_h_Mg in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(66) is k_off_CaM in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(67) is k_off_myo_Ca in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(68) is k_off_myo_Mg in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(69) is k_off_SR in component Cytosolic_Ca_Buffers (per_millis). C CONSTS(70) is k_on_TnC_l in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(71) is k_on_TnC_h_Ca in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(72) is k_on_TnC_h_Mg in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(73) is k_on_CaM in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(74) is k_on_myo_Ca in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(75) is k_on_myo_Mg in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(76) is k_on_SR in component Cytosolic_Ca_Buffers (per_milliM_per_millis). C CONSTS(77) is B_max_TnC_low in component Cytosolic_Ca_Buffers (millimolar). C CONSTS(78) is B_max_TnC_high in component Cytosolic_Ca_Buffers (millimolar). C CONSTS(79) is B_max_CaM in component Cytosolic_Ca_Buffers (millimolar). C CONSTS(80) is B_max_myosin in component Cytosolic_Ca_Buffers (millimolar). C CONSTS(81) 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(106) is B_max_SL_low_j in component Junctional_and_SL_Ca_Buffers (millimolar). C CONSTS(107) is B_max_SL_low_sl in component Junctional_and_SL_Ca_Buffers (millimolar). C CONSTS(108) is B_max_SL_high_j in component Junctional_and_SL_Ca_Buffers (millimolar). C CONSTS(109) is B_max_SL_high_sl in component Junctional_and_SL_Ca_Buffers (millimolar). C CONSTS(82) is k_off_sl_l in component Junctional_and_SL_Ca_Buffers (per_millis). C CONSTS(83) is k_off_sl_h in component Junctional_and_SL_Ca_Buffers (per_millis). C CONSTS(84) is k_on_sl_l in component Junctional_and_SL_Ca_Buffers (per_milliM_per_millis). C CONSTS(85) 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(86) is B_max_Na_j in component Sodium_Concentrations (millimolar). C CONSTS(87) is B_max_Na_sl in component Sodium_Concentrations (millimolar). C CONSTS(88) is k_off_Na in component Sodium_Concentrations (per_millis). C CONSTS(89) 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(110) is B_max_csqn in component Calcium_Concentrations (millimolar). C STATES(39) is Csqn_b in component Calcium_Concentrations (millimolar). C CONSTS(90) is k_off_csqn in component Calcium_Concentrations (per_millis). C CONSTS(91) is k_on_csqn in component Calcium_Concentrations (per_milliM_per_millis). C CONSTS(92) is Cl_i in component Chlorine_Concentrations (millimolar). C CONSTS(93) is Cl_o in component Chlorine_Concentrations (millimolar). C CONSTS(94) 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.413368e1 CONSTS(1) = 8314 CONSTS(2) = 310 CONSTS(3) = 96485 CONSTS(4) = 1.381e-10 CONSTS(5) = 0 CONSTS(6) = 40 CONSTS(7) = 1000 CONSTS(8) = 1 CONSTS(9) = 0.11 CONSTS(10) = 0.9 CONSTS(11) = 100e-5 CONSTS(12) = 10.25e-5 CONSTS(13) = 8.2413e-13 CONSTS(14) = 3.7243e-12 CONSTS(15) = 1.8313e-14 CONSTS(16) = 1.6386e-12 CONSTS(17) = 18.86 STATES(2) = 2.163678e-3 STATES(3) = 7.13497e-1 STATES(4) = 7.128671e-1 CONSTS(18) = 0.597e-3 CONSTS(19) = 140 CONSTS(20) = 0.99 CONSTS(21) = 5.4 STATES(5) = 1.007931e1 STATES(6) = 1.00781e1 CONSTS(22) = 1.5 CONSTS(23) = 11 STATES(7) = 1.516232e-2 CONSTS(24) = 0.035 STATES(8) = 3.549354e-3 CONSTS(25) = 0.0035 CONSTS(26) = 0.0035 CONSTS(27) = 0.002 STATES(9) = 3.584727e-4 STATES(10) = 8.087629e-1 STATES(11) = 3.584625e-4 STATES(12) = 9.999976e-1 CONSTS(28) = 0.0156 CONSTS(29) = 0.1144 CONSTS(30) = 0.57153 STATES(13) = 2.038197e-4 STATES(14) = 1.184305e-4 CONSTS(31) = 0.054813 CONSTS(32) = 100e-3 CONSTS(33) = 9e-3 STATES(15) = 1.871177e-6 STATES(16) = 9.804391e-1 STATES(17) = 9.99401e-1 STATES(18) = 2.847118e-2 STATES(19) = 1.692189e-2 CONSTS(34) = 138 CONSTS(35) = 1.8 CONSTS(36) = 1.9887e-4 CONSTS(37) = 3.0375e-9 CONSTS(38) = 5.4675e-8 CONSTS(39) = 0.15e-3 CONSTS(40) = 3.59e-3 CONSTS(41) = 1.3 CONSTS(42) = 12.29 CONSTS(43) = 87.5 CONSTS(44) = 0.32 CONSTS(45) = 0.27 CONSTS(46) = 4.5 CONSTS(47) = 0.5e-3 CONSTS(48) = 0.0673 CONSTS(49) = 5.513e-4 CONSTS(50) = 15 CONSTS(51) = 1 STATES(20) = 6.093596e-1 CONSTS(52) = 0.45 CONSTS(53) = 10 CONSTS(54) = 0.5 CONSTS(55) = 0.005 CONSTS(56) = 0.06 CONSTS(57) = 25 CONSTS(58) = 0.246e-3 CONSTS(59) = 1.7 STATES(21) = 8.886338e-1 STATES(22) = 1.126209e-6 STATES(23) = 1.411382e-7 STATES(24) = 9.658067e-5 CONSTS(60) = 5.3114e-3 CONSTS(61) = 1.787 CONSTS(62) = 1 STATES(25) = 9.757237e-3 STATES(26) = 1.225914e-1 STATES(27) = 8.12201e-3 STATES(28) = 3.267494e-4 STATES(29) = 2.520383e-3 STATES(30) = 1.369529e-1 STATES(31) = 2.373753e-3 CONSTS(63) = 19.6e-3 CONSTS(64) = 0.032e-3 CONSTS(65) = 3.33e-3 CONSTS(66) = 238e-3 CONSTS(67) = 0.46e-3 CONSTS(68) = 0.057e-3 CONSTS(69) = 60e-3 CONSTS(70) = 32.7 CONSTS(71) = 2.37 CONSTS(72) = 3e-3 CONSTS(73) = 34 CONSTS(74) = 13.8 CONSTS(75) = 0.0157 CONSTS(76) = 100 CONSTS(77) = 70e-3 CONSTS(78) = 140e-3 CONSTS(79) = 24e-3 CONSTS(80) = 140e-3 CONSTS(81) = 17.1e-3 CONSTS(82) = 1.3 CONSTS(83) = 30e-3 CONSTS(84) = 100 CONSTS(85) = 100 STATES(32) = 8.563314e-3 STATES(33) = 1.097424e-2 STATES(34) = 8.053908e-2 STATES(35) = 1.235381e-1 CONSTS(86) = 7.561 CONSTS(87) = 1.65 CONSTS(88) = 1e-3 CONSTS(89) = 0.1e-3 STATES(36) = 1.007825e1 STATES(37) = 3.796195e0 STATES(38) = 8.283308e-1 STATES(39) = 1.258048e0 CONSTS(90) = 65 CONSTS(91) = 100 CONSTS(92) = 15 CONSTS(93) = 150 CONSTS(94) = 0.01833 CONSTS(95) = CONSTS(3)/( CONSTS(1)*CONSTS(2)) CONSTS(96) = 1.00000 - CONSTS(9) CONSTS(97) = 1.00000 - CONSTS(10) CONSTS(98) = 3.14159265358979*CONSTS(12) ** 2.00000*CONSTS(11) CONSTS(99) = (EXP(CONSTS(19)/67.3000) - 1.00000)/7.00000 CONSTS(100) = 0.650000*CONSTS(98) CONSTS(101) = (1.00000/CONSTS(95))*log(CONSTS(21)/CONSTS(34)) CONSTS(102) = (1.00000/CONSTS(95))*log(CONSTS(92)/CONSTS(93)) CONSTS(103) = 0.0200000*CONSTS(98) CONSTS(104) = 0.0350000*CONSTS(98) CONSTS(105) = 0.000539000*CONSTS(98) CONSTS(106) = ( 0.000460000*CONSTS(100))/CONSTS(105) CONSTS(107) = ( 0.0374000*CONSTS(100))/CONSTS(103) CONSTS(108) = ( 0.000165000*CONSTS(100))/CONSTS(105) CONSTS(109) = ( 0.0134000*CONSTS(100))/CONSTS(103) CONSTS(110) = ( 0.140000*CONSTS(100))/CONSTS(104) 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(16)/CONSTS(100))*(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(50) - (CONSTS(50) - CONSTS(51))/(1.00000+CONSTS(52)/STATES(20) ** 2.50000) ALGBRC(26) = CONSTS(53)/ALGBRC(12) ALGBRC(33) = CONSTS(54)*ALGBRC(12) RATES(22) = ( ALGBRC(26)*STATES(13) ** 2.00000*STATES(21) - CONSTS(56)*STATES(22)) - ( ALGBRC(33)*STATES(13)*STATES(22) - CONSTS(55)*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(55)*ALGBRC(40) - ALGBRC(33)*STATES(13)*STATES(21)) - ( ALGBRC(26)*STATES(13) ** 2.00000*STATES(21) - CONSTS(56)*STATES(22)) RATES(23) = ( ALGBRC(33)*STATES(13)*STATES(22) - CONSTS(55)*STATES(23)) - ( CONSTS(56)*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(70)*STATES(24)*(CONSTS(77) - STATES(25)) - CONSTS(63)*STATES(25) RATES(32) = CONSTS(84)*STATES(13)*(CONSTS(106) - STATES(32)) - CONSTS(82)*STATES(32) RATES(26) = CONSTS(71)*STATES(24)*((CONSTS(78) - STATES(26)) - STATES(27)) - CONSTS(64)*STATES(26) RATES(39) = CONSTS(91)*STATES(20)*(CONSTS(110) - STATES(39)) - CONSTS(90)*STATES(39) RATES(33) = CONSTS(84)*STATES(14)*(CONSTS(107) - STATES(33)) - CONSTS(82)*STATES(33) RATES(27) = CONSTS(72)*CONSTS(62)*((CONSTS(78) - STATES(26)) - STATES(27)) - CONSTS(65)*STATES(27) ALGBRC(87) = CONSTS(57)*STATES(22)*(STATES(20) - STATES(13)) ALGBRC(90) = 5.34800e-06*(STATES(20) - STATES(13)) ALGBRC(89) = ( CONSTS(60)*(STATES(24)/CONSTS(58) ** CONSTS(61) - STATES(20)/CONSTS(59) ** CONSTS(61)))/(1.00000+STATES(24)/CONSTS(58) ** CONSTS(61)+STATES(20)/CONSTS(59) ** CONSTS(61)) RATES(20) = (ALGBRC(89) - (( ALGBRC(90)*CONSTS(100))/CONSTS(104)+ALGBRC(87))) - RATES(39) RATES(34) = CONSTS(85)*STATES(13)*(CONSTS(108) - STATES(34)) - CONSTS(83)*STATES(34) RATES(28) = CONSTS(73)*STATES(24)*(CONSTS(79) - STATES(28)) - CONSTS(66)*STATES(28) RATES(35) = CONSTS(85)*STATES(14)*(CONSTS(109) - STATES(35)) - CONSTS(83)*STATES(35) RATES(29) = CONSTS(74)*STATES(24)*((CONSTS(80) - STATES(29)) - STATES(30)) - CONSTS(67)*STATES(29) RATES(30) = CONSTS(75)*CONSTS(62)*((CONSTS(80) - STATES(29)) - STATES(30)) - CONSTS(68)*STATES(30) RATES(31) = CONSTS(76)*STATES(24)*(CONSTS(81) - STATES(31)) - CONSTS(69)*STATES(31) RATES(37) = CONSTS(89)*STATES(5)*(CONSTS(86) - STATES(37)) - CONSTS(88)*STATES(37) RATES(38) = CONSTS(89)*STATES(6)*(CONSTS(87) - STATES(38)) - CONSTS(88)*STATES(38) ALGBRC(93) = RATES(25)+RATES(26)+RATES(27)+RATES(28)+RATES(29)+RATES(30)+RATES(31) RATES(24) = (( - ALGBRC(89)*CONSTS(104))/CONSTS(100) - ALGBRC(93))+ (CONSTS(14)/CONSTS(100))*(STATES(14) - STATES(24)) ALGBRC(94) = (1.00000/CONSTS(95))*log(CONSTS(19)/STATES(5)) ALGBRC(95) = CONSTS(9)*CONSTS(17)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(94)) ALGBRC(96) = CONSTS(9)*CONSTS(18)*(STATES(1) - ALGBRC(94)) ALGBRC(25) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(95))+ 0.0365000*CONSTS(99)*EXP( - STATES(1)*CONSTS(95))) ALGBRC(32) = ( CONSTS(9)*CONSTS(20)*ALGBRC(25)*CONSTS(21))/( (1.00000+CONSTS(23)/STATES(5) ** 4.00000)*(CONSTS(21)+CONSTS(22))) ALGBRC(64) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(5)*EXP( STATES(1)*CONSTS(95)) - CONSTS(19)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(70) = CONSTS(10)*ALGBRC(64)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18)) ALGBRC(72) = 1.00000/(1.00000+CONSTS(39)/STATES(13) ** 2.00000) ALGBRC(76) = EXP( CONSTS(45)*STATES(1)*CONSTS(95))*STATES(5) ** 3.00000*CONSTS(35) ALGBRC(78) = EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95))*CONSTS(19) ** 3.00000*STATES(13) ALGBRC(80) = CONSTS(40)*CONSTS(19) ** 3.00000*(1.00000+STATES(5)/CONSTS(42) ** 3.00000)+ CONSTS(43) ** 3.00000*STATES(13)*(1.00000+STATES(13)/CONSTS(40))+ CONSTS(41)*STATES(5) ** 3.00000+ STATES(5) ** 3.00000*CONSTS(35)+ CONSTS(19) ** 3.00000*STATES(13) ALGBRC(82) = ( CONSTS(9)*CONSTS(46)*ALGBRC(72)*(ALGBRC(76) - ALGBRC(78)))/( ALGBRC(80)*(1.00000+ CONSTS(44)*EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95)))) ALGBRC(97) = ALGBRC(95)+ALGBRC(96)+ 3.00000*ALGBRC(82)+ 3.00000*ALGBRC(32)+ALGBRC(70) RATES(5) = (( - ALGBRC(97)*CONSTS(4))/( CONSTS(105)*CONSTS(3))+ (CONSTS(15)/CONSTS(105))*(STATES(6) - STATES(5))) - RATES(37) ALGBRC(98) = (1.00000/CONSTS(95))*log(CONSTS(19)/STATES(6)) ALGBRC(99) = CONSTS(96)*CONSTS(17)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(98)) ALGBRC(101) = CONSTS(96)*CONSTS(18)*(STATES(1) - ALGBRC(98)) ALGBRC(39) = ( CONSTS(96)*CONSTS(20)*ALGBRC(25)*CONSTS(21))/( (1.00000+CONSTS(23)/STATES(6) ** 4.00000)*(CONSTS(21)+CONSTS(22))) ALGBRC(65) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(6)*EXP( STATES(1)*CONSTS(95)) - CONSTS(19)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(71) = CONSTS(97)*ALGBRC(65)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19)) ALGBRC(74) = 1.00000/(1.00000+CONSTS(39)/STATES(14) ** 2.00000) ALGBRC(77) = EXP( CONSTS(45)*STATES(1)*CONSTS(95))*STATES(6) ** 3.00000*CONSTS(35) ALGBRC(79) = EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95))*CONSTS(19) ** 3.00000*STATES(14) ALGBRC(81) = CONSTS(40)*CONSTS(19) ** 3.00000*(1.00000+STATES(6)/CONSTS(42) ** 3.00000)+ CONSTS(43) ** 3.00000*STATES(14)*(1.00000+STATES(14)/CONSTS(40))+ CONSTS(41)*STATES(6) ** 3.00000+ STATES(6) ** 3.00000*CONSTS(35)+ CONSTS(19) ** 3.00000*STATES(14) ALGBRC(83) = ( CONSTS(96)*CONSTS(46)*ALGBRC(74)*(ALGBRC(77) - ALGBRC(79)))/( ALGBRC(81)*(1.00000+ CONSTS(44)*EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95)))) ALGBRC(103) = ALGBRC(99)+ALGBRC(101)+ 3.00000*ALGBRC(83)+ 3.00000*ALGBRC(39)+ALGBRC(71) RATES(6) = (( - ALGBRC(103)*CONSTS(4))/( CONSTS(103)*CONSTS(3))+ (CONSTS(15)/CONSTS(103))*(STATES(5) - STATES(6))+ (CONSTS(16)/CONSTS(103))*(STATES(36) - STATES(6))) - RATES(38) ALGBRC(61) = ( CONSTS(36)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(13)*EXP( 2.00000*STATES(1)*CONSTS(95)) - CONSTS(35)))/(EXP( 2.00000*STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(66) = CONSTS(10)*ALGBRC(61)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18)) ALGBRC(84) = ( CONSTS(9)*CONSTS(48)*STATES(13) ** 1.60000)/(CONSTS(47) ** 1.60000+STATES(13) ** 1.60000) ALGBRC(110) = (1.00000/( 2.00000*CONSTS(95)))*log(CONSTS(35)/STATES(13)) ALGBRC(111) = CONSTS(9)*CONSTS(49)*(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(105)*2.00000*CONSTS(3))+ (CONSTS(13)/CONSTS(105))*(STATES(14) - STATES(13))) - ALGBRC(91))+( ALGBRC(87)*CONSTS(104))/CONSTS(105)+( ALGBRC(90)*CONSTS(100))/CONSTS(105) ALGBRC(62) = ( CONSTS(36)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(14)*EXP( 2.00000*STATES(1)*CONSTS(95)) - CONSTS(35)))/(EXP( 2.00000*STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(67) = CONSTS(97)*ALGBRC(62)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19)) ALGBRC(86) = ( CONSTS(96)*CONSTS(48)*STATES(14) ** 1.60000)/(CONSTS(47) ** 1.60000+STATES(14) ** 1.60000) ALGBRC(113) = (1.00000/( 2.00000*CONSTS(95)))*log(CONSTS(35)/STATES(14)) ALGBRC(114) = CONSTS(96)*CONSTS(49)*(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(103)*2.00000*CONSTS(3))+ (CONSTS(13)/CONSTS(103))*(STATES(13) - STATES(14))+ (CONSTS(14)/CONSTS(103))*(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(24)* (CONSTS(21)/5.40000) ** (1.0 / 2)*STATES(7)*ALGBRC(45)*(STATES(1) - CONSTS(101)) ALGBRC(105) = (1.00000/CONSTS(95))*log((CONSTS(21)+ CONSTS(94)*CONSTS(19))/(CONSTS(34)+ CONSTS(94)*STATES(36))) ALGBRC(106) = CONSTS(9)*CONSTS(25)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105)) ALGBRC(107) = CONSTS(96)*CONSTS(26)*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(27)*ALGBRC(47)*(STATES(1) - CONSTS(101)) ALGBRC(49) = CONSTS(28)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(101)) ALGBRC(50) = CONSTS(29)*STATES(11)*STATES(12)*(STATES(1) - CONSTS(101)) ALGBRC(51) = ALGBRC(49)+ALGBRC(50) ALGBRC(52) = 4.09380/(1.00000+EXP( 0.121650*((STATES(1) - CONSTS(101)) - 49.9344))) ALGBRC(53) = ( 15.7197*EXP( 0.0673900*((STATES(1) - CONSTS(101)) - 3.25710))+EXP( 0.0617500*((STATES(1) - CONSTS(101)) - 594.310)))/(1.00000+EXP( - 0.162850*((STATES(1) - CONSTS(101))+14.2067))) ALGBRC(54) = ALGBRC(52)/(ALGBRC(52)+ALGBRC(53)) ALGBRC(55) = CONSTS(30)* (CONSTS(21)/5.40000) ** (1.0 / 2)*ALGBRC(54)*(STATES(1) - CONSTS(101)) ALGBRC(63) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(95)*( CONSTS(34)*EXP( STATES(1)*CONSTS(95)) - CONSTS(21)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(69) = ALGBRC(63)*STATES(15)*STATES(16)*STATES(17)*( CONSTS(10)*(1.00000 - STATES(18))+ CONSTS(97)*(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(9)*CONSTS(31))/(1.00000+CONSTS(32)/STATES(13)))*(STATES(1) - CONSTS(102)) ALGBRC(57) = (( CONSTS(96)*CONSTS(31))/(1.00000+CONSTS(32)/STATES(14)))*(STATES(1) - CONSTS(102)) ALGBRC(58) = ALGBRC(56)+ALGBRC(57) ALGBRC(59) = CONSTS(33)*(STATES(1) - CONSTS(102)) 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 - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(6), 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(50) - (CONSTS(50) - CONSTS(51))/(1.00000+CONSTS(52)/STATES(20) ** 2.50000) ALGBRC(26) = CONSTS(53)/ALGBRC(12) ALGBRC(33) = CONSTS(54)*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(57)*STATES(22)*(STATES(20) - STATES(13)) ALGBRC(90) = 5.34800e-06*(STATES(20) - STATES(13)) ALGBRC(89) = ( CONSTS(60)*(STATES(24)/CONSTS(58) ** CONSTS(61) - STATES(20)/CONSTS(59) ** CONSTS(61)))/(1.00000+STATES(24)/CONSTS(58) ** CONSTS(61)+STATES(20)/CONSTS(59) ** CONSTS(61)) ALGBRC(93) = RATES(25)+RATES(26)+RATES(27)+RATES(28)+RATES(29)+RATES(30)+RATES(31) ALGBRC(94) = (1.00000/CONSTS(95))*log(CONSTS(19)/STATES(5)) ALGBRC(95) = CONSTS(9)*CONSTS(17)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(94)) ALGBRC(96) = CONSTS(9)*CONSTS(18)*(STATES(1) - ALGBRC(94)) ALGBRC(25) = 1.00000/(1.00000+ 0.124500*EXP( - 0.100000*STATES(1)*CONSTS(95))+ 0.0365000*CONSTS(99)*EXP( - STATES(1)*CONSTS(95))) ALGBRC(32) = ( CONSTS(9)*CONSTS(20)*ALGBRC(25)*CONSTS(21))/( (1.00000+CONSTS(23)/STATES(5) ** 4.00000)*(CONSTS(21)+CONSTS(22))) ALGBRC(64) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(5)*EXP( STATES(1)*CONSTS(95)) - CONSTS(19)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(70) = CONSTS(10)*ALGBRC(64)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18)) ALGBRC(72) = 1.00000/(1.00000+CONSTS(39)/STATES(13) ** 2.00000) ALGBRC(76) = EXP( CONSTS(45)*STATES(1)*CONSTS(95))*STATES(5) ** 3.00000*CONSTS(35) ALGBRC(78) = EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95))*CONSTS(19) ** 3.00000*STATES(13) ALGBRC(80) = CONSTS(40)*CONSTS(19) ** 3.00000*(1.00000+STATES(5)/CONSTS(42) ** 3.00000)+ CONSTS(43) ** 3.00000*STATES(13)*(1.00000+STATES(13)/CONSTS(40))+ CONSTS(41)*STATES(5) ** 3.00000+ STATES(5) ** 3.00000*CONSTS(35)+ CONSTS(19) ** 3.00000*STATES(13) ALGBRC(82) = ( CONSTS(9)*CONSTS(46)*ALGBRC(72)*(ALGBRC(76) - ALGBRC(78)))/( ALGBRC(80)*(1.00000+ CONSTS(44)*EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95)))) ALGBRC(97) = ALGBRC(95)+ALGBRC(96)+ 3.00000*ALGBRC(82)+ 3.00000*ALGBRC(32)+ALGBRC(70) ALGBRC(98) = (1.00000/CONSTS(95))*log(CONSTS(19)/STATES(6)) ALGBRC(99) = CONSTS(96)*CONSTS(17)*STATES(2) ** 3.00000*STATES(3)*STATES(4)*(STATES(1) - ALGBRC(98)) ALGBRC(101) = CONSTS(96)*CONSTS(18)*(STATES(1) - ALGBRC(98)) ALGBRC(39) = ( CONSTS(96)*CONSTS(20)*ALGBRC(25)*CONSTS(21))/( (1.00000+CONSTS(23)/STATES(6) ** 4.00000)*(CONSTS(21)+CONSTS(22))) ALGBRC(65) = ( CONSTS(37)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(6)*EXP( STATES(1)*CONSTS(95)) - CONSTS(19)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(71) = CONSTS(97)*ALGBRC(65)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19)) ALGBRC(74) = 1.00000/(1.00000+CONSTS(39)/STATES(14) ** 2.00000) ALGBRC(77) = EXP( CONSTS(45)*STATES(1)*CONSTS(95))*STATES(6) ** 3.00000*CONSTS(35) ALGBRC(79) = EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95))*CONSTS(19) ** 3.00000*STATES(14) ALGBRC(81) = CONSTS(40)*CONSTS(19) ** 3.00000*(1.00000+STATES(6)/CONSTS(42) ** 3.00000)+ CONSTS(43) ** 3.00000*STATES(14)*(1.00000+STATES(14)/CONSTS(40))+ CONSTS(41)*STATES(6) ** 3.00000+ STATES(6) ** 3.00000*CONSTS(35)+ CONSTS(19) ** 3.00000*STATES(14) ALGBRC(83) = ( CONSTS(96)*CONSTS(46)*ALGBRC(74)*(ALGBRC(77) - ALGBRC(79)))/( ALGBRC(81)*(1.00000+ CONSTS(44)*EXP( (CONSTS(45) - 1.00000)*STATES(1)*CONSTS(95)))) ALGBRC(103) = ALGBRC(99)+ALGBRC(101)+ 3.00000*ALGBRC(83)+ 3.00000*ALGBRC(39)+ALGBRC(71) ALGBRC(61) = ( CONSTS(36)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(13)*EXP( 2.00000*STATES(1)*CONSTS(95)) - CONSTS(35)))/(EXP( 2.00000*STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(66) = CONSTS(10)*ALGBRC(61)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(18)) ALGBRC(84) = ( CONSTS(9)*CONSTS(48)*STATES(13) ** 1.60000)/(CONSTS(47) ** 1.60000+STATES(13) ** 1.60000) ALGBRC(110) = (1.00000/( 2.00000*CONSTS(95)))*log(CONSTS(35)/STATES(13)) ALGBRC(111) = CONSTS(9)*CONSTS(49)*(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(36)*STATES(1)*CONSTS(3)*CONSTS(95)*( STATES(14)*EXP( 2.00000*STATES(1)*CONSTS(95)) - CONSTS(35)))/(EXP( 2.00000*STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(67) = CONSTS(97)*ALGBRC(62)*STATES(15)*STATES(16)*STATES(17)*(1.00000 - STATES(19)) ALGBRC(86) = ( CONSTS(96)*CONSTS(48)*STATES(14) ** 1.60000)/(CONSTS(47) ** 1.60000+STATES(14) ** 1.60000) ALGBRC(113) = (1.00000/( 2.00000*CONSTS(95)))*log(CONSTS(35)/STATES(14)) ALGBRC(114) = CONSTS(96)*CONSTS(49)*(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(24)* (CONSTS(21)/5.40000) ** (1.0 / 2)*STATES(7)*ALGBRC(45)*(STATES(1) - CONSTS(101)) ALGBRC(105) = (1.00000/CONSTS(95))*log((CONSTS(21)+ CONSTS(94)*CONSTS(19))/(CONSTS(34)+ CONSTS(94)*STATES(36))) ALGBRC(106) = CONSTS(9)*CONSTS(25)*STATES(8) ** 2.00000*(STATES(1) - ALGBRC(105)) ALGBRC(107) = CONSTS(96)*CONSTS(26)*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(27)*ALGBRC(47)*(STATES(1) - CONSTS(101)) ALGBRC(49) = CONSTS(28)*STATES(9)*STATES(10)*(STATES(1) - CONSTS(101)) ALGBRC(50) = CONSTS(29)*STATES(11)*STATES(12)*(STATES(1) - CONSTS(101)) ALGBRC(51) = ALGBRC(49)+ALGBRC(50) ALGBRC(52) = 4.09380/(1.00000+EXP( 0.121650*((STATES(1) - CONSTS(101)) - 49.9344))) ALGBRC(53) = ( 15.7197*EXP( 0.0673900*((STATES(1) - CONSTS(101)) - 3.25710))+EXP( 0.0617500*((STATES(1) - CONSTS(101)) - 594.310)))/(1.00000+EXP( - 0.162850*((STATES(1) - CONSTS(101))+14.2067))) ALGBRC(54) = ALGBRC(52)/(ALGBRC(52)+ALGBRC(53)) ALGBRC(55) = CONSTS(30)* (CONSTS(21)/5.40000) ** (1.0 / 2)*ALGBRC(54)*(STATES(1) - CONSTS(101)) ALGBRC(63) = ( CONSTS(38)*STATES(1)*CONSTS(3)*CONSTS(95)*( CONSTS(34)*EXP( STATES(1)*CONSTS(95)) - CONSTS(21)))/(EXP( STATES(1)*CONSTS(95)) - 1.00000) ALGBRC(69) = ALGBRC(63)*STATES(15)*STATES(16)*STATES(17)*( CONSTS(10)*(1.00000 - STATES(18))+ CONSTS(97)*(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(9)*CONSTS(31))/(1.00000+CONSTS(32)/STATES(13)))*(STATES(1) - CONSTS(102)) ALGBRC(57) = (( CONSTS(96)*CONSTS(31))/(1.00000+CONSTS(32)/STATES(14)))*(STATES(1) - CONSTS(102)) ALGBRC(58) = ALGBRC(56)+ALGBRC(57) ALGBRC(59) = CONSTS(33)*(STATES(1) - CONSTS(102)) 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 - CONSTS(5)) - INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(6), 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