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 146 entries in the algebraic variable array. C There are a total of 41 entries in each of the rate and state variable arrays. C There are a total of 99 entries in the constant variable array. C C C VOI is time in component environment (millisecond). C STATES(1) is Vm in component membrane (millivolt). C CONSTS(1) is R in component membrane (coulomb_millivolt_per_kelvin_millimole). C CONSTS(2) is T in component membrane (kelvin). C CONSTS(3) is F in component membrane (coulomb_per_millimole). C CONSTS(4) is Cm in component membrane (picoF). C ALGBRC(111) is i_tot in component membrane (picoA). C ALGBRC(95) is i_I in component membrane (picoA). C ALGBRC(50) is i_Na in component sodium_current (picoA). C ALGBRC(55) is i_Ca_L in component L_type_Ca_channel (picoA). C ALGBRC(59) is i_Ca_T in component T_type_Ca_channel (picoA). C ALGBRC(81) is i_K1 in component time_independent_potassium_current (picoA). C ALGBRC(82) is i_Kr in component rapid_time_dependent_potassium_current (picoA). C ALGBRC(66) is i_st in component st_channel (picoA). C ALGBRC(70) is i_ha in component ha_channel (picoA). C ALGBRC(83) is i_KACh in component ACh_dependent_potassium_current (picoA). C ALGBRC(109) is i_NaK in component sodium_potassium_pump (picoA). C ALGBRC(100) is i_NaCa in component sodium_calcium_exchanger (picoA). C ALGBRC(86) is i_bNSC in component background_NSC_current (picoA). C ALGBRC(94) is i_Cab in component background_Cab_current (picoA). C ALGBRC(87) is i_Kpl in component background_Kpl_current (picoA). C ALGBRC(91) is i_lCa in component background_lCa_current (picoA). C ALGBRC(93) is i_KATP in component background_KATP_current (picoA). C CONSTS(5) is Nao in component external_ion_concentrations (millimolar). C CONSTS(6) is Cao in component external_ion_concentrations (millimolar). C CONSTS(7) is Ko in component external_ion_concentrations (millimolar). C STATES(2) is Nai in component internal_ion_concentrations (millimolar). C ALGBRC(28) is Cai in component internal_ion_concentrations (millimolar). C STATES(3) is Ki in component internal_ion_concentrations (millimolar). C CONSTS(8) is Vi in component internal_ion_concentrations (micrometre3). C ALGBRC(112) is i_net_Na in component internal_ion_concentrations (picoA). C ALGBRC(113) is i_net_K in component internal_ion_concentrations (picoA). C ALGBRC(102) is i_net_Ca in component internal_ion_concentrations (picoA). C ALGBRC(45) is i_Na_Na in component sodium_current (picoA). C ALGBRC(53) is i_CaL_Na in component L_type_Ca_channel (picoA). C ALGBRC(61) is i_st_Na in component st_channel (picoA). C ALGBRC(68) is i_ha_Na in component ha_channel (picoA). C ALGBRC(85) is i_bNSC_Na in component background_NSC_current (picoA). C ALGBRC(90) is i_lCa_Na in component background_lCa_current (picoA). C ALGBRC(48) is i_Na_K in component sodium_current (picoA). C ALGBRC(54) is i_CaL_K in component L_type_Ca_channel (picoA). C ALGBRC(64) is i_st_K in component st_channel (picoA). C ALGBRC(69) is i_ha_K in component ha_channel (picoA). C ALGBRC(84) is i_bNSC_K in component background_NSC_current (picoA). C ALGBRC(89) is i_lCa_K in component background_lCa_current (picoA). C ALGBRC(52) is i_CaL_Ca in component L_type_Ca_channel (picoA). C ALGBRC(132) is i_RyR in component RyR_channel (picoA). C ALGBRC(121) is i_SR_U in component SR_calcium_pump (picoA). C ALGBRC(126) is i_SR_L in component SR_L_current (picoA). C ALGBRC(146) is dCaidt in component NL_model (millimolar_per_millisecond). C CONSTS(9) is CMDN_max in component internal_ion_concentrations (millimolar). C CONSTS(10) is K_mCMDN in component internal_ion_concentrations (millimolar). C STATES(4) is Ca_Total in component internal_ion_concentrations (millimolar). C ALGBRC(1) is b1 in component internal_ion_concentrations (millimolar). C ALGBRC(15) is c1 in component internal_ion_concentrations (millimolar2). C ALGBRC(32) is CF_Na in component constant_field_equations (millimolar). C ALGBRC(37) is CF_Ca in component constant_field_equations (millimolar). C ALGBRC(42) is CF_K in component constant_field_equations (millimolar). C STATES(5) is ATPi in component ATP_production (millimolar). C ALGBRC(123) is dATPdt in component NL_model (millimolar_per_millisecond). C CONSTS(11) is ProducingRate_Max in component ATP_production (per_millisecond). C CONSTS(12) is Adenosine_Total in component ATP_production (millimolar). C CONSTS(13) is P_Na in component sodium_current (picoA_per_millimolar). C STATES(6) is p_AP_Na in component sodium_current_voltage_dependent_gate (dimensionless). C STATES(7) is y in component sodium_current_ultra_slow_gate (dimensionless). C ALGBRC(2) is p_RI_Na in component sodium_current_voltage_dependent_gate (dimensionless). C STATES(8) is p_RP_Na in component sodium_current_voltage_dependent_gate (dimensionless). C STATES(9) is p_AI_Na in component sodium_current_voltage_dependent_gate (dimensionless). C ALGBRC(16) is k_RP_AP in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(29) is k_AP_RP in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(43) is k_RI_AI in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(38) is k_AI_RI in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(33) is k_AP_AI in component sodium_current_voltage_dependent_gate (per_millisecond). C CONSTS(14) is k_AI_AP in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(46) is k_RP_RI in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(49) is k_RI_RP in component sodium_current_voltage_dependent_gate (per_millisecond). C ALGBRC(3) is alpha_y in component sodium_current_ultra_slow_gate (per_millisecond). C ALGBRC(17) is beta_y in component sodium_current_ultra_slow_gate (per_millisecond). C ALGBRC(51) is p_open_CaL in component L_type_Ca_channel (dimensionless). C ALGBRC(57) is CaDiadic in component L_type_Ca_channel_Ca_dependent_gate (picoA). C CONSTS(15) is P_CaL in component L_type_Ca_channel (picoA_per_millimolar). C STATES(10) is p_AP_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C STATES(11) is p_U in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C STATES(12) is p_UCa in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C STATES(13) is y in component L_type_Ca_channel_ultra_slow_gate (dimensionless). C ALGBRC(4) is p_RI_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C STATES(14) is p_RP_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C STATES(15) is p_AI_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C ALGBRC(18) is k_RP_AP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(30) is k_AP_RP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(34) is k_RI_AI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(39) is k_AI_RI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C CONSTS(16) is k_AP_AI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C CONSTS(17) is k_AI_AP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(44) is k_RP_RI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(47) is k_RI_RP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(56) is iCaL in component L_type_Ca_channel_Ca_dependent_gate (picoA). C ALGBRC(58) is Cacm in component L_type_Ca_channel_Ca_dependent_gate (millimolar). C ALGBRC(67) is p_CCa in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C STATES(16) is p_C in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C CONSTS(18) is k_CCa_UCa in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(19) is k_UCa_CCa in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(20) is k_C_U in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(21) is k_U_C in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(90) is k_UCa_U in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(22) is k_U_UCa in component L_type_Ca_channel_Ca_dependent_gate (per_millimolar_millisecond). C CONSTS(23) is k_CCa_C in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C CONSTS(24) is k_C_CCa in component L_type_Ca_channel_Ca_dependent_gate (per_millimolar_millisecond). C ALGBRC(60) is CaEffC in component L_type_Ca_channel_Ca_dependent_gate (millimolar). C ALGBRC(62) is CaEffU in component L_type_Ca_channel_Ca_dependent_gate (millimolar). C ALGBRC(65) is k_UUCa_Ca in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C ALGBRC(63) is k_CCCa_Ca in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond). C ALGBRC(5) is alpha_y in component L_type_Ca_channel_ultra_slow_gate (per_millisecond). C ALGBRC(19) is beta_y in component L_type_Ca_channel_ultra_slow_gate (per_millisecond). C CONSTS(25) is P_CaT in component T_type_Ca_channel (picoA_per_millimolar). C STATES(17) is y1 in component T_type_Ca_channel_y1_gate (dimensionless). C STATES(18) is y2 in component T_type_Ca_channel_y2_gate (dimensionless). C ALGBRC(6) is alpha_y1 in component T_type_Ca_channel_y1_gate (per_millisecond). C ALGBRC(20) is beta_y1 in component T_type_Ca_channel_y1_gate (per_millisecond). C ALGBRC(7) is alpha_y2 in component T_type_Ca_channel_y2_gate (per_millisecond). C ALGBRC(21) is beta_y2 in component T_type_Ca_channel_y2_gate (per_millisecond). C STATES(19) is y1 in component st_channel_y1_gate (dimensionless). C STATES(20) is y2 in component st_channel_y2_gate (dimensionless). C STATES(21) is y3 in component st_channel_y3_gate (dimensionless). C CONSTS(26) is P_st_Na in component st_channel (picoA_per_millimolar). C CONSTS(27) is P_st_K in component st_channel (picoA_per_millimolar). C ALGBRC(8) is alpha_y1 in component st_channel_y1_gate (per_millisecond). C ALGBRC(22) is beta_y1 in component st_channel_y1_gate (per_millisecond). C ALGBRC(35) is alpha_y2 in component st_channel_y2_gate (per_millisecond). C ALGBRC(40) is beta_y2 in component st_channel_y2_gate (per_millisecond). C ALGBRC(9) is alpha_y3 in component st_channel_y3_gate (per_millisecond). C ALGBRC(23) is beta_y3 in component st_channel_y3_gate (per_millisecond). C STATES(22) is C1 in component ha_channel_voltage_dependent_gate (dimensionless). C STATES(23) is C2 in component ha_channel_voltage_dependent_gate (dimensionless). C CONSTS(28) is P_ha_Na in component ha_channel (picoA_per_millimolar). C CONSTS(29) is P_ha_K in component ha_channel (picoA_per_millimolar). C STATES(24) is O1 in component ha_channel_voltage_dependent_gate (dimensionless). C STATES(25) is O2 in component ha_channel_voltage_dependent_gate (dimensionless). C ALGBRC(10) is O3 in component ha_channel_voltage_dependent_gate (dimensionless). C ALGBRC(24) is alpha in component ha_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(31) is beta in component ha_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(36) is mu in component ha_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(41) is lambda in component ha_channel_voltage_dependent_gate (per_millisecond). C ALGBRC(71) is E_K in component time_independent_potassium_current (millivolt). C CONSTS(91) is g_K1 in component time_independent_potassium_current (nanoS). C CONSTS(30) is P_K1_0 in component time_independent_potassium_current (nanoS_per_picoF). C ALGBRC(75) is fO in component time_independent_potassium_current (dimensionless). C ALGBRC(76) is fO2 in component time_independent_potassium_current (dimensionless). C ALGBRC(78) is fO3 in component time_independent_potassium_current (dimensionless). C ALGBRC(80) is fO4 in component time_independent_potassium_current (dimensionless). C ALGBRC(74) is fB in component time_independent_potassium_current (dimensionless). C ALGBRC(72) is mu in component time_independent_potassium_current (per_millisecond). C ALGBRC(73) is lambda in component time_independent_potassium_current (per_millisecond). C STATES(26) is y in component time_independent_potassium_current_y_gate (dimensionless). C ALGBRC(77) is alpha_y in component time_independent_potassium_current_y_gate (per_millisecond). C ALGBRC(79) is beta_y in component time_independent_potassium_current_y_gate (per_millisecond). C CONSTS(92) is g_Kr in component rapid_time_dependent_potassium_current (nanoS). C CONSTS(31) is P_Kr in component rapid_time_dependent_potassium_current (nanoS_per_picoF). C STATES(27) is y1 in component rapid_time_dependent_potassium_current_y1_gate (dimensionless). C STATES(28) is y2 in component rapid_time_dependent_potassium_current_y2_gate (dimensionless). C STATES(29) is y3 in component rapid_time_dependent_potassium_current_y3_gate (dimensionless). C ALGBRC(11) is alpha_y1 in component rapid_time_dependent_potassium_current_y1_gate (per_millisecond). C ALGBRC(25) is beta_y1 in component rapid_time_dependent_potassium_current_y1_gate (per_millisecond). C ALGBRC(12) is alpha_y2 in component rapid_time_dependent_potassium_current_y2_gate (per_millisecond). C ALGBRC(26) is beta_y2 in component rapid_time_dependent_potassium_current_y2_gate (per_millisecond). C ALGBRC(13) is alpha_y3 in component rapid_time_dependent_potassium_current_y3_gate (per_millisecond). C ALGBRC(27) is beta_y3 in component rapid_time_dependent_potassium_current_y3_gate (per_millisecond). C STATES(30) is y in component ACh_dependent_potassium_current_y_gate (dimensionless). C CONSTS(32) is g_ACh in component ACh_dependent_potassium_current (nanoS). C CONSTS(93) is alpha_y in component ACh_dependent_potassium_current_y_gate (per_millisecond). C ALGBRC(14) is beta_y in component ACh_dependent_potassium_current_y_gate (per_millisecond). C CONSTS(33) is ACh in component ACh_dependent_potassium_current_y_gate (millimolar). C CONSTS(34) is Km_ACh in component ACh_dependent_potassium_current_y_gate (millimolar). C CONSTS(35) is P_bNSC in component background_NSC_current (picoA_per_millimolar). C CONSTS(94) is P_Kpl in component background_Kpl_current (nanoS_per_millimolar). C CONSTS(36) is P_lCa in component background_lCa_current (picoA_per_millimolar). C ALGBRC(88) is p_open in component background_lCa_current (dimensionless). C ALGBRC(92) is p_open in component background_KATP_current (dimensionless). C CONSTS(95) is gamma in component background_KATP_current (nanoS). C CONSTS(37) is P_KATP in component background_KATP_current (nanoS_per_picoF). C CONSTS(38) is N in component background_KATP_current (picoF). C CONSTS(39) is P_Cab in component background_Cab_current (picoA_per_millimolar). C CONSTS(96) is p_E2Na in component sodium_calcium_exchanger (dimensionless). C ALGBRC(96) is p_E1Na in component sodium_calcium_exchanger (dimensionless). C ALGBRC(97) is p_E1Ca in component sodium_calcium_exchanger (dimensionless). C CONSTS(99) is p_E2Ca in component sodium_calcium_exchanger (dimensionless). C ALGBRC(98) is k1 in component sodium_calcium_exchanger (per_millisecond). C ALGBRC(99) is k2 in component sodium_calcium_exchanger (per_millisecond). C CONSTS(40) is k3 in component sodium_calcium_exchanger (per_millisecond). C CONSTS(41) is k4 in component sodium_calcium_exchanger (per_millisecond). C CONSTS(42) is Km_Nai in component sodium_calcium_exchanger (millimolar). C CONSTS(43) is Km_Nao in component sodium_calcium_exchanger (millimolar). C CONSTS(44) is Km_Cai in component sodium_calcium_exchanger (millimolar). C CONSTS(45) is Km_Cao in component sodium_calcium_exchanger (millimolar). C STATES(31) is y in component sodium_calcium_exchanger_y_gate (dimensionless). C CONSTS(46) is P_NaCa in component sodium_calcium_exchanger (picoA_per_picoF). C CONSTS(47) is Partition in component sodium_calcium_exchanger (dimensionless). C ALGBRC(101) is alpha_y in component sodium_calcium_exchanger_y_gate (per_millisecond). C ALGBRC(103) is beta_y in component sodium_calcium_exchanger_y_gate (per_millisecond). C ALGBRC(108) is p_E2Na in component sodium_potassium_pump (dimensionless). C ALGBRC(104) is p_E1Na in component sodium_potassium_pump (dimensionless). C ALGBRC(105) is p_E1K in component sodium_potassium_pump (dimensionless). C ALGBRC(110) is p_E2K in component sodium_potassium_pump (dimensionless). C ALGBRC(106) is k1 in component sodium_potassium_pump (per_millisecond). C CONSTS(48) is k2 in component sodium_potassium_pump (per_millisecond). C CONSTS(49) is k3 in component sodium_potassium_pump (per_millisecond). C CONSTS(50) is k4 in component sodium_potassium_pump (per_millisecond). C CONSTS(51) is Km_Nai in component sodium_potassium_pump (millimolar). C CONSTS(52) is Km_Nao in component sodium_potassium_pump (millimolar). C CONSTS(53) is Km_Ki in component sodium_potassium_pump (millimolar). C CONSTS(54) is Km_Ko in component sodium_potassium_pump (millimolar). C CONSTS(55) is Km_ATP in component sodium_potassium_pump (millimolar). C ALGBRC(107) is Nao_Eff in component sodium_potassium_pump (millimolar). C STATES(32) is y in component sodium_potassium_pump_y_gate (dimensionless). C CONSTS(56) is P_NaK in component sodium_potassium_pump (picoA_per_picoF). C ALGBRC(115) is alpha_y in component sodium_potassium_pump_y_gate (per_millisecond). C ALGBRC(117) is beta_y in component sodium_potassium_pump_y_gate (per_millisecond). C ALGBRC(116) is p_E2Ca in component SR_calcium_pump (dimensionless). C ALGBRC(114) is p_E1Ca in component SR_calcium_pump (dimensionless). C ALGBRC(118) is p_E1 in component SR_calcium_pump (dimensionless). C ALGBRC(119) is p_E2 in component SR_calcium_pump (dimensionless). C CONSTS(57) is k1 in component SR_calcium_pump (per_millisecond). C ALGBRC(120) is k2 in component SR_calcium_pump (per_millisecond). C CONSTS(58) is k3 in component SR_calcium_pump (per_millisecond). C CONSTS(59) is k4 in component SR_calcium_pump (per_millisecond). C CONSTS(60) is Km_CaSR in component SR_calcium_pump (millimolar). C CONSTS(61) is Km_CaCyto in component SR_calcium_pump (millimolar). C CONSTS(62) is Km_ATP in component SR_calcium_pump (millimolar). C CONSTS(63) is i_max in component SR_calcium_pump (picoA). C STATES(33) is Caup in component Ca_concentrations_in_SR (millimolar). C STATES(34) is y in component SR_calcium_pump_y_gate (dimensionless). C ALGBRC(122) is alpha_y in component SR_calcium_pump_y_gate (per_millisecond). C ALGBRC(124) is beta_y in component SR_calcium_pump_y_gate (per_millisecond). C CONSTS(64) is P_RyR in component RyR_channel (picoA_per_millimolar). C ALGBRC(125) is k1 in component RyR_channel (per_millisecond). C ALGBRC(130) is k2 in component RyR_channel (per_millisecond). C ALGBRC(131) is k3 in component RyR_channel (per_millisecond). C CONSTS(65) is k4 in component RyR_channel (per_millisecond). C STATES(35) is p_open_RyR in component RyR_channel (dimensionless). C STATES(36) is p_close_RyR in component RyR_channel (dimensionless). C ALGBRC(129) is Carel in component Ca_concentrations_in_SR (millimolar). C CONSTS(66) is Diadid_Factor in component RyR_channel (per_picoA_millisecond). C ALGBRC(133) is i_SR_T in component SR_T_current (picoA). C CONSTS(67) is P_SR_T in component SR_T_current (picoA_per_millimolar). C CONSTS(68) is P_SR_L in component SR_L_current (picoA_per_millimolar). C STATES(37) is Ca_Total in component Ca_concentrations_in_SR (millimolar). C CONSTS(69) is V_rel in component Ca_concentrations_in_SR (micrometre3). C CONSTS(70) is V_up in component Ca_concentrations_in_SR (micrometre3). C CONSTS(71) is CSQN_max in component Ca_concentrations_in_SR (millimolar). C CONSTS(72) is K_mCSQN in component Ca_concentrations_in_SR (millimolar). C ALGBRC(127) is b1 in component Ca_concentrations_in_SR (millimolar). C ALGBRC(128) is c1 in component Ca_concentrations_in_SR (millimolar2). C CONSTS(97) is EffFraction in component NL_model (dimensionless). C STATES(38) is pCa in component NL_model (dimensionless). C STATES(39) is pCaCB in component NL_model (dimensionless). C STATES(40) is pCB in component NL_model (dimensionless). C ALGBRC(136) is p in component NL_model (dimensionless). C CONSTS(73) is T_t in component NL_model (millimolar). C ALGBRC(140) is Q_a in component NL_model (per_millisecond). C ALGBRC(138) is Q_b in component NL_model (per_millisecond). C ALGBRC(142) is Q_r in component NL_model (per_millisecond). C ALGBRC(143) is Q_d in component NL_model (per_millisecond). C ALGBRC(144) is Q_d1 in component NL_model (per_millisecond). C ALGBRC(145) is Q_d2 in component NL_model (per_millisecond). C CONSTS(74) is Y_1 in component NL_model (per_millimolar_millisecond). C CONSTS(75) is Y_2 in component NL_model (per_millisecond). C CONSTS(76) is Y_3 in component NL_model (per_millisecond). C CONSTS(77) is Y_4 in component NL_model (per_millisecond). C CONSTS(78) is Y_d in component NL_model (millisecond_per_micrometre2). C CONSTS(79) is Z_1 in component NL_model (per_millisecond). C CONSTS(80) is Z_2 in component NL_model (per_millisecond). C CONSTS(81) is Z_3 in component NL_model (per_millimolar_millisecond). C ALGBRC(134) is h in component NL_model (micrometre). C CONSTS(82) is L_a in component NL_model (micrometre). C CONSTS(83) is L in component NL_model (micrometre). C ALGBRC(139) is ForceCB in component NL_model (mN_per_mm2). C STATES(41) is X in component NL_model (micrometre). C ALGBRC(137) is NewCBF in component NL_model (mN_per_mm2_micrometre). C ALGBRC(135) is CBBound in component NL_model (millimolar). C CONSTS(84) is KForceEC in component NL_model (mN_per_mm2_micrometre5). C CONSTS(85) is ZeroForceEL in component NL_model (micrometre). C CONSTS(86) is KForceLinearEc in component NL_model (mN_per_mm2_micrometre). C CONSTS(87) is ForceFactor in component NL_model (mN_per_mm2_micrometre_millimolar). C CONSTS(98) is ForceEcomp in component NL_model (mN_per_mm2). C CONSTS(88) is B in component NL_model (per_millisecond). C CONSTS(89) is h_c in component NL_model (micrometre). C ALGBRC(141) is ForceExt in component NL_model (mN_per_mm2). C RATES(1) is d/dt Vm in component membrane (millivolt). C RATES(2) is d/dt Nai in component internal_ion_concentrations (millimolar). C RATES(3) is d/dt Ki in component internal_ion_concentrations (millimolar). C RATES(4) is d/dt Ca_Total in component internal_ion_concentrations (millimolar). C RATES(5) is d/dt ATPi in component ATP_production (millimolar). C RATES(8) is d/dt p_RP_Na in component sodium_current_voltage_dependent_gate (dimensionless). C RATES(6) is d/dt p_AP_Na in component sodium_current_voltage_dependent_gate (dimensionless). C RATES(9) is d/dt p_AI_Na in component sodium_current_voltage_dependent_gate (dimensionless). C RATES(7) is d/dt y in component sodium_current_ultra_slow_gate (dimensionless). C RATES(14) is d/dt p_RP_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C RATES(10) is d/dt p_AP_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C RATES(15) is d/dt p_AI_CaL in component L_type_Ca_channel_voltage_dependent_gate (dimensionless). C RATES(11) is d/dt p_U in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C RATES(12) is d/dt p_UCa in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C RATES(16) is d/dt p_C in component L_type_Ca_channel_Ca_dependent_gate (dimensionless). C RATES(13) is d/dt y in component L_type_Ca_channel_ultra_slow_gate (dimensionless). C RATES(17) is d/dt y1 in component T_type_Ca_channel_y1_gate (dimensionless). C RATES(18) is d/dt y2 in component T_type_Ca_channel_y2_gate (dimensionless). C RATES(19) is d/dt y1 in component st_channel_y1_gate (dimensionless). C RATES(20) is d/dt y2 in component st_channel_y2_gate (dimensionless). C RATES(21) is d/dt y3 in component st_channel_y3_gate (dimensionless). C RATES(24) is d/dt O1 in component ha_channel_voltage_dependent_gate (dimensionless). C RATES(25) is d/dt O2 in component ha_channel_voltage_dependent_gate (dimensionless). C RATES(22) is d/dt C1 in component ha_channel_voltage_dependent_gate (dimensionless). C RATES(23) is d/dt C2 in component ha_channel_voltage_dependent_gate (dimensionless). C RATES(26) is d/dt y in component time_independent_potassium_current_y_gate (dimensionless). C RATES(27) is d/dt y1 in component rapid_time_dependent_potassium_current_y1_gate (dimensionless). C RATES(28) is d/dt y2 in component rapid_time_dependent_potassium_current_y2_gate (dimensionless). C RATES(29) is d/dt y3 in component rapid_time_dependent_potassium_current_y3_gate (dimensionless). C RATES(30) is d/dt y in component ACh_dependent_potassium_current_y_gate (dimensionless). C RATES(31) is d/dt y in component sodium_calcium_exchanger_y_gate (dimensionless). C RATES(32) is d/dt y in component sodium_potassium_pump_y_gate (dimensionless). C RATES(34) is d/dt y in component SR_calcium_pump_y_gate (dimensionless). C RATES(35) is d/dt p_open_RyR in component RyR_channel (dimensionless). C RATES(36) is d/dt p_close_RyR in component RyR_channel (dimensionless). C RATES(37) is d/dt Ca_Total in component Ca_concentrations_in_SR (millimolar). C RATES(33) is d/dt Caup in component Ca_concentrations_in_SR (millimolar). C RATES(41) is d/dt X in component NL_model (micrometre). C RATES(38) is d/dt pCa in component NL_model (dimensionless). C RATES(39) is d/dt pCaCB in component NL_model (dimensionless). C RATES(40) is d/dt pCB in component NL_model (dimensionless). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = -42.987671350068005 CONSTS(1) = 8.3143 CONSTS(2) = 310 CONSTS(3) = 96.4867 CONSTS(4) = 32 CONSTS(5) = 139.99977217549068 CONSTS(6) = 1.7999668787634697 CONSTS(7) = 5.399964171929831 STATES(2) = 4.6808027298969925 STATES(3) = 139.22099524072777 CONSTS(8) = 2513.25 CONSTS(9) = 0.05 CONSTS(10) = 0.00238 STATES(4) = 0.003834100516470642 STATES(5) = 4.929308214702568 CONSTS(11) = 0.003 CONSTS(12) = 5 CONSTS(13) = 12 STATES(6) = 5.9898633533987405e-5 STATES(7) = 0.02691630436856424 STATES(8) = 2.2262336890623873e-4 STATES(9) = 0.9715915109181413 CONSTS(14) = 8.75e-5 CONSTS(15) = 2112 STATES(10) = 0.010773205742759775 STATES(11) = 0.018412749759150325 STATES(12) = 2.0870506374005444e-4 STATES(13) = 0.8256168329801671 STATES(14) = 0.42833044151885075 STATES(15) = 0.5607106503080197 CONSTS(16) = 0.004 CONSTS(17) = 0.001 STATES(16) = 0.05040310692927668 CONSTS(18) = 0.0003 CONSTS(19) = 0.35 CONSTS(20) = 0.143 CONSTS(21) = 0.35 CONSTS(22) = 6.954 CONSTS(23) = 0.0042 CONSTS(24) = 6.954 CONSTS(25) = 612 STATES(17) = 0.04679297803939428 STATES(18) = 0.004220636800632867 STATES(19) = 0.5375825176280309 STATES(20) = 0.3404904569858008 STATES(21) = 0.5596860236041188 CONSTS(26) = 0.236 CONSTS(27) = 0.138 STATES(22) = 0.9999616586460943 STATES(23) = 3.757612529902994e-5 CONSTS(28) = 1.821 CONSTS(29) = 7.7286 STATES(24) = 7.531035735927369e-7 STATES(25) = 1.1957124194086763e-8 CONSTS(30) = 0.003125 STATES(26) = 0.54434573606543 CONSTS(31) = 0.1045 STATES(27) = 0.7520936568166604 STATES(28) = 0.6128486629679218 STATES(29) = 0.6786243213825371 STATES(30) = 0.018249236601344903 CONSTS(32) = 4.32 CONSTS(33) = 0 CONSTS(34) = 0.0042 CONSTS(35) = 0.152 CONSTS(36) = 0.44 CONSTS(37) = 0.0236 CONSTS(38) = 1000 CONSTS(39) = 0.008 CONSTS(40) = 1 CONSTS(41) = 1 CONSTS(42) = 8.75 CONSTS(43) = 87.5 CONSTS(44) = 0.00138 CONSTS(45) = 1.38 STATES(31) = 0.9218642754061908 CONSTS(46) = 6.81 CONSTS(47) = 0.32 CONSTS(48) = 0.04 CONSTS(49) = 0.01 CONSTS(50) = 0.165 CONSTS(51) = 4.05 CONSTS(52) = 69.8 CONSTS(53) = 32.88 CONSTS(54) = 0.258 CONSTS(55) = 0.094 STATES(32) = 0.6755031143352256 CONSTS(56) = 21 CONSTS(57) = 0.01 CONSTS(58) = 1 CONSTS(59) = 0.01 CONSTS(60) = 0.08 CONSTS(61) = 0.0008 CONSTS(62) = 0.1 CONSTS(63) = 4875 STATES(33) = 3.539652240305415 STATES(34) = 0.8709732600283433 CONSTS(64) = 1860 CONSTS(65) = 0.000849 STATES(35) = 0.011424172905335395 STATES(36) = 0.007354039651246729 CONSTS(66) = -3 CONSTS(67) = 11.58 CONSTS(68) = 13.77 STATES(37) = 5.185011336959393 CONSTS(69) = 1.50795 CONSTS(70) = 3.769875 CONSTS(71) = 10 CONSTS(72) = 0.8 STATES(38) = 0.24298053499618683 STATES(39) = 0.02896443650702234 STATES(40) = 0.0020275353474589783 CONSTS(73) = 0.007 CONSTS(74) = 39 CONSTS(75) = 0.0039 CONSTS(76) = 0.03 CONSTS(77) = 0.12 CONSTS(78) = 0.027 CONSTS(79) = 0.03 CONSTS(80) = 0.0039 CONSTS(81) = 1560 CONSTS(82) = 1.17 CONSTS(83) = 0.9602399999150041 STATES(41) = 0.9552449999150041 CONSTS(84) = 140000 CONSTS(85) = 0.97 CONSTS(86) = 200 CONSTS(87) = 1800000 CONSTS(88) = 1.2 CONSTS(89) = 0.005 CONSTS(90) = ( CONSTS(23)*CONSTS(20)*CONSTS(22)*CONSTS(19))/( CONSTS(21)*CONSTS(24)*CONSTS(18)) CONSTS(91) = CONSTS(30)*CONSTS(4)*CONSTS(7)/5.40000 ** 0.400000 CONSTS(92) = CONSTS(31)*CONSTS(4)*CONSTS(7)/5.40000 ** 0.200000 CONSTS(93) = TERNRY(CONSTS(33).EQ.0.00000, 0.000247500, 0.0123200/(1.00000+CONSTS(34)/CONSTS(33))+0.000247500) CONSTS(94) = 0.00700000*CONSTS(7)/5.40000 ** 0.160000 CONSTS(95) = CONSTS(37)*CONSTS(38)*CONSTS(7)/1.00000 ** 0.240000 CONSTS(96) = 1.00000/(1.00000+ CONSTS(43)/CONSTS(5) ** 3.00000*(1.00000+CONSTS(6)/CONSTS(45))) CONSTS(97) = EXP( - 20.0000*CONSTS(83) - CONSTS(82) ** 2.00000) CONSTS(98) = CONSTS(84)*CONSTS(85) - CONSTS(83) ** 5.00000+ CONSTS(86)*(CONSTS(85) - CONSTS(83)) CONSTS(99) = 1.00000/(1.00000+ (CONSTS(45)/CONSTS(6))*(1.00000+CONSTS(5)/CONSTS(43) ** 3.00000)) RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(14) = 0.0100000*EXP( 0.0133000*(STATES(1)+40.0000)) RATES(30) = CONSTS(93)*(1.00000 - STATES(30)) - ALGBRC(14)*STATES(30) ALGBRC(3) = 1.00000/( 9.00000e+09*EXP(STATES(1)/5.00000)+ 8000.00*EXP(STATES(1)/100.000)) ALGBRC(17) = 1.00000/( 0.0140000*EXP(- STATES(1)/5.00000)+ 4000.00*EXP(- STATES(1)/100.000)) RATES(7) = ALGBRC(3)*(1.00000 - STATES(7)) - ALGBRC(17)*STATES(7) ALGBRC(5) = 1.00000/( 250000.*EXP(STATES(1)/9.00000)+ 58.0000*EXP(STATES(1)/65.0000)) ALGBRC(19) = 1.00000/( 1800.00*EXP(- STATES(1)/14.0000)+ 66.0000*EXP(- STATES(1)/65.0000)) RATES(13) = ALGBRC(5)*(1.00000 - STATES(13)) - ALGBRC(19)*STATES(13) ALGBRC(6) = 1.00000/( 0.0190000*EXP(- STATES(1)/5.60000)+ 0.820000*EXP(- STATES(1)/250.000)) ALGBRC(20) = 1.00000/( 40.0000*EXP(STATES(1)/6.30000)+ 1.50000*EXP(STATES(1)/10000.0)) RATES(17) = ALGBRC(6)*(1.00000 - STATES(17)) - ALGBRC(20)*STATES(17) ALGBRC(7) = 1.00000/( 62000.0*EXP(STATES(1)/10.1000)+ 30.0000*EXP(STATES(1)/3000.00)) ALGBRC(21) = 1.00000/( 0.000600000*EXP(- STATES(1)/6.70000)+ 1.20000*EXP(- STATES(1)/25.0000)) RATES(18) = ALGBRC(7)*(1.00000 - STATES(18)) - ALGBRC(21)*STATES(18) ALGBRC(8) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000)) ALGBRC(22) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000)) RATES(19) = ALGBRC(8)*(1.00000 - STATES(19)) - ALGBRC(22)*STATES(19) ALGBRC(9) = 1.00000/( 400000.*EXP(STATES(1)/9.00000)+ 60.0000*EXP(STATES(1)/65.0000)) ALGBRC(23) = 1.00000/( 700.000*EXP(- STATES(1)/14.0000)+ 60.0000*EXP(- STATES(1)/65.0000)) RATES(21) = ALGBRC(9)*(1.00000 - STATES(21)) - ALGBRC(23)*STATES(21) ALGBRC(11) = 1.00000/( 20.0000*EXP(- STATES(1)/11.5000)+ 5.00000*EXP(- STATES(1)/300.000)) ALGBRC(25) = 1.00000/( 160.000*EXP(STATES(1)/28.0000)+ 200.000*EXP(STATES(1)/1000.00))+1.00000/( 2500.00*EXP(STATES(1)/20.0000)) RATES(27) = ALGBRC(11)*(1.00000 - STATES(27)) - ALGBRC(25)*STATES(27) ALGBRC(12) = 1.00000/( 200.000*EXP(- STATES(1)/13.0000)+ 20.0000*EXP(- STATES(1)/300.000)) ALGBRC(26) = 1.00000/( 1600.00*EXP(STATES(1)/28.0000)+ 2000.00*EXP(STATES(1)/1000.00))+1.00000/( 10000.0*EXP(STATES(1)/20.0000)) RATES(28) = ALGBRC(12)*(1.00000 - STATES(28)) - ALGBRC(26)*STATES(28) ALGBRC(13) = 1.00000/( 10.0000*EXP(STATES(1)/17.0000)+ 2.50000*EXP(STATES(1)/300.000)) ALGBRC(27) = 1.00000/( 0.350000*EXP(- STATES(1)/17.0000)+ 2.00000*EXP(- STATES(1)/150.000)) RATES(29) = ALGBRC(13)*(1.00000 - STATES(29)) - ALGBRC(27)*STATES(29) ALGBRC(18) = 1.00000/( 0.270000*EXP(- STATES(1)/5.90000)+ 1.50000*EXP(- STATES(1)/65.0000)) ALGBRC(30) = 1.00000/( 480.000*EXP(STATES(1)/7.00000)+ 2.20000*EXP(STATES(1)/65.0000)) RATES(10) = ( STATES(14)*ALGBRC(18)+ STATES(15)*CONSTS(17)) - STATES(10)*(ALGBRC(30)+CONSTS(16)) ALGBRC(24) = 1.00000/( 3500.00*EXP(STATES(1)/16.8000)+ 0.300000*EXP(STATES(1)/400.000)) ALGBRC(31) = 1.00000/( 4.00000*EXP(- STATES(1)/14.0000)+ 2.00000*EXP(- STATES(1)/400.000)) RATES(24) = ( ALGBRC(24)*STATES(23)+ ALGBRC(31)*STATES(25)) - (ALGBRC(24)+ALGBRC(31))*STATES(24) ALGBRC(10) = (((1.00000 - STATES(22)) - STATES(23)) - STATES(24)) - STATES(25) RATES(25) = ( ALGBRC(24)*STATES(24)+ ALGBRC(31)*ALGBRC(10)) - (ALGBRC(24)+ALGBRC(31))*STATES(25) ALGBRC(16) = 1.00000/( 0.102700*EXP(- STATES(1)/8.00000)+ 0.250000*EXP(- STATES(1)/50.0000)) ALGBRC(29) = 1.00000/( 26.0000*EXP(STATES(1)/17.0000)+ 0.0200000*EXP(STATES(1)/800.000)) ALGBRC(33) = 1.00000/( 0.800000*EXP(- STATES(1)/400.000)) RATES(6) = ( STATES(8)*ALGBRC(16)+ STATES(9)*CONSTS(14)) - STATES(6)*(ALGBRC(29)+ALGBRC(33)) ALGBRC(4) = ((1.00000 - STATES(10)) - STATES(14)) - STATES(15) ALGBRC(34) = 1.00000/( 0.00180000*EXP(- STATES(1)/7.40000)+ 2.00000*EXP(- STATES(1)/100.000)) ALGBRC(39) = 1.00000/( 2.20000e+06*EXP(STATES(1)/7.40000)+ 11.0000*EXP(STATES(1)/100.000)) RATES(15) = ( ALGBRC(4)*ALGBRC(34)+ STATES(10)*CONSTS(16)) - STATES(15)*(ALGBRC(39)+CONSTS(17)) ALGBRC(35) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000)) ALGBRC(1) = (CONSTS(9) - STATES(4))+CONSTS(10) ALGBRC(15) = CONSTS(10)*STATES(4) ALGBRC(28) = ( (ALGBRC(1) ** 2.00000+ 4.00000*ALGBRC(15)) ** (1.0 / 2) - ALGBRC(1))/2.00000 ALGBRC(40) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+( ALGBRC(28)*2.50000)/(1.00000+EXP(- STATES(1)/5.00000)) RATES(20) = ALGBRC(35)*(1.00000 - STATES(20)) - ALGBRC(40)*STATES(20) ALGBRC(36) = 1.00000/( 4.50000e+07*EXP(STATES(1)/8.00000)+ 500.000*EXP(STATES(1)/200.000)) ALGBRC(41) = 1.00000/( 10.5000*EXP(- STATES(1)/16.4000)+ 0.400000*EXP(- STATES(1)/400.000)) RATES(22) = ALGBRC(41)*STATES(23) - ALGBRC(36)*STATES(22) RATES(23) = ( ALGBRC(36)*STATES(22)+ ALGBRC(31)*STATES(24)) - (ALGBRC(24)+ALGBRC(41))*STATES(23) ALGBRC(2) = ((1.00000 - STATES(8)) - STATES(6)) - STATES(9) ALGBRC(43) = 1.00000/( 0.000102700*EXP(- STATES(1)/8.00000)+ 5.00000*EXP(- STATES(1)/400.000)) ALGBRC(38) = 1.00000/( 1300.00*EXP(STATES(1)/20.0000)+ 0.0400000*EXP(STATES(1)/800.000)) RATES(9) = ( ALGBRC(2)*ALGBRC(43)+ STATES(6)*ALGBRC(33)) - STATES(9)*(ALGBRC(38)+CONSTS(14)) ALGBRC(44) = 0.0400000/(1.00000+( CONSTS(17)*ALGBRC(30)*ALGBRC(34))/( CONSTS(16)*ALGBRC(18)*ALGBRC(39))) ALGBRC(47) = 0.0400000 - ALGBRC(44) RATES(14) = ( STATES(10)*ALGBRC(30)+ ALGBRC(4)*ALGBRC(47)) - STATES(14)*(ALGBRC(44)+ALGBRC(18)) ALGBRC(46) = 0.0100000/(1.00000+( CONSTS(14)*ALGBRC(29)*ALGBRC(43))/( ALGBRC(33)*ALGBRC(16)*ALGBRC(38))) ALGBRC(49) = 0.0100000 - ALGBRC(46) RATES(8) = ( STATES(6)*ALGBRC(29)+ ALGBRC(2)*ALGBRC(49)) - STATES(8)*(ALGBRC(46)+ALGBRC(16)) ALGBRC(37) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(6), ( (( 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(ALGBRC(28) - CONSTS(6)*EXP(( - 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(56) = 0.0676000*ALGBRC(37) ALGBRC(58) = ALGBRC(28) - 0.300000*ALGBRC(56) ALGBRC(60) = ALGBRC(58)*STATES(10) ALGBRC(62) = ALGBRC(60)+ ALGBRC(28)*(1.00000 - STATES(10)) ALGBRC(65) = CONSTS(22)*ALGBRC(62) RATES(11) = ( STATES(16)*CONSTS(20)+ STATES(12)*CONSTS(90)) - STATES(11)*(ALGBRC(65)+CONSTS(21)) ALGBRC(67) = ((1.00000 - STATES(16)) - STATES(11)) - STATES(12) RATES(12) = ( STATES(11)*ALGBRC(65)+ ALGBRC(67)*CONSTS(18)) - STATES(12)*(CONSTS(19)+CONSTS(90)) RATES(16) = ( ALGBRC(67)*CONSTS(23)+ STATES(11)*CONSTS(21)) - STATES(16)*(CONSTS(20)+ CONSTS(24)*ALGBRC(58)*STATES(10)) ALGBRC(71) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(7)/STATES(3)) ALGBRC(77) = 1.00000/( 8000.00*EXP(((STATES(1) - ALGBRC(71)) - 97.0000)/8.50000)+ 7.00000*EXP(((STATES(1) - ALGBRC(71)) - 97.0000)/300.000)) ALGBRC(72) = ( 0.750000*EXP( 0.0350000*((STATES(1) - ALGBRC(71)) - 10.0000)))/(1.00000+EXP( 0.0150000*((STATES(1) - ALGBRC(71)) - 140.000))) ALGBRC(73) = ( 3.00000*EXP( - 0.0480000*((STATES(1) - ALGBRC(71)) - 10.0000))*(1.00000+EXP( 0.0640000*((STATES(1) - ALGBRC(71)) - 38.0000))))/(1.00000+EXP( 0.0300000*((STATES(1) - ALGBRC(71)) - 70.0000))) ALGBRC(75) = ALGBRC(73)/(ALGBRC(72)+ALGBRC(73)) ALGBRC(79) = ( ALGBRC(75) ** 4.00000*1.00000)/( 0.000140000*EXP(- ((STATES(1) - ALGBRC(71)) - 97.0000)/9.10000)+ 0.200000*EXP(- ((STATES(1) - ALGBRC(71)) - 97.0000)/500.000)) RATES(26) = ALGBRC(77)*(1.00000 - STATES(26)) - ALGBRC(79)*STATES(26) ALGBRC(99) = 1.00000*EXP(( (CONSTS(47) - 1.00000)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(101) = ALGBRC(99)*CONSTS(96)+ CONSTS(41)*CONSTS(99) ALGBRC(96) = 1.00000/(1.00000+ CONSTS(42)/STATES(2) ** 3.00000*(1.00000+ALGBRC(28)/CONSTS(44))) ALGBRC(97) = 1.00000/(1.00000+ (CONSTS(44)/ALGBRC(28))*(1.00000+STATES(2)/CONSTS(42) ** 3.00000)) ALGBRC(98) = 1.00000*EXP(( CONSTS(47)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(103) = ALGBRC(98)*ALGBRC(96)+ CONSTS(40)*ALGBRC(97) RATES(31) = ALGBRC(101)*(1.00000 - STATES(31)) - ALGBRC(103)*STATES(31) ALGBRC(32) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(5), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(2) - CONSTS(5)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(85) = CONSTS(35)*ALGBRC(32) ALGBRC(42) = TERNRY(STATES(1).EQ.0.00000, STATES(3), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(3) - CONSTS(7)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(84) = 0.400000*CONSTS(35)*ALGBRC(42) ALGBRC(86) = ALGBRC(84)+ALGBRC(85) ALGBRC(94) = CONSTS(39)*ALGBRC(37) ALGBRC(87) = TERNRY(STATES(1).EQ.- 3.00000, CONSTS(94)*ALGBRC(42)*13.0077, ( CONSTS(94)*ALGBRC(42)*(STATES(1)+3.00000))/(1.00000 - EXP(- (STATES(1)+3.00000)/13.0000))) ALGBRC(88) = 1.00000/(1.00000+0.00120000/ALGBRC(28) ** 3.00000) ALGBRC(90) = CONSTS(36)*ALGBRC(32)*ALGBRC(88) ALGBRC(89) = CONSTS(36)*ALGBRC(42)*ALGBRC(88) ALGBRC(91) = ALGBRC(89)+ALGBRC(90) ALGBRC(92) = 0.800000/(1.00000+STATES(5)/0.100000 ** 2.00000) ALGBRC(93) = CONSTS(95)*(STATES(1) - ALGBRC(71))*ALGBRC(92) ALGBRC(95) = ALGBRC(86)+ALGBRC(94)+ALGBRC(87)+ALGBRC(91)+ALGBRC(93) ALGBRC(45) = CONSTS(13)*ALGBRC(32)*STATES(6)*STATES(7) ALGBRC(48) = 0.100000*CONSTS(13)*ALGBRC(42)*STATES(6)*STATES(7) ALGBRC(50) = ALGBRC(45)+ALGBRC(48) ALGBRC(51) = ( STATES(10)*(STATES(11)+STATES(12))*STATES(13))/(1.00000+1.40000/STATES(5) ** 3.00000) ALGBRC(53) = 1.85000e-05*CONSTS(15)*ALGBRC(32)*ALGBRC(51) ALGBRC(54) = 0.000365000*CONSTS(15)*ALGBRC(42)*ALGBRC(51) ALGBRC(52) = CONSTS(15)*ALGBRC(37)*ALGBRC(51) ALGBRC(55) = ALGBRC(53)+ALGBRC(52)+ALGBRC(54) ALGBRC(59) = CONSTS(25)*ALGBRC(37)*STATES(17)*STATES(18) ALGBRC(74) = ALGBRC(72)/(ALGBRC(72)+ALGBRC(73)) ALGBRC(76) = 2.00000*ALGBRC(75) ** 2.00000*ALGBRC(74) ** 2.00000 ALGBRC(78) = (8.00000/3.00000)*ALGBRC(75) ** 3.00000*ALGBRC(74) ALGBRC(80) = ALGBRC(75) ** 4.00000 ALGBRC(81) = CONSTS(91)*(STATES(1) - ALGBRC(71))*(ALGBRC(80)+ALGBRC(78)+ALGBRC(76))*STATES(26) ALGBRC(82) = CONSTS(92)*(STATES(1) - ALGBRC(71))*( 0.600000*STATES(27)+ 0.400000*STATES(28))*STATES(29) ALGBRC(61) = CONSTS(26)*ALGBRC(32)*STATES(19)*STATES(20)*STATES(21) ALGBRC(64) = CONSTS(27)*ALGBRC(42)*STATES(19)*STATES(20)*STATES(21) ALGBRC(66) = ALGBRC(61)+ALGBRC(64) ALGBRC(68) = CONSTS(28)*ALGBRC(32)*((1.00000 - STATES(22)) - STATES(23)) ALGBRC(69) = CONSTS(29)*ALGBRC(42)*((1.00000 - STATES(22)) - STATES(23)) ALGBRC(70) = ALGBRC(68)+ALGBRC(69) ALGBRC(83) = (( CONSTS(32)*(STATES(1) - ALGBRC(71)))/(1.00000+EXP((STATES(1)+20.0000)/20.0000)))*STATES(30) ALGBRC(107) = CONSTS(5)*EXP(( - 0.820000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(108) = 1.00000/(1.00000+ CONSTS(52)/ALGBRC(107) ** 1.06000*(1.00000+CONSTS(7)/CONSTS(54) ** 1.12000)) ALGBRC(104) = 1.00000/(1.00000+ CONSTS(51)/STATES(2) ** 1.06000*(1.00000+STATES(3)/CONSTS(53) ** 1.12000)) ALGBRC(106) = 0.370000/(1.00000+CONSTS(55)/STATES(5)) ALGBRC(109) = CONSTS(56)*CONSTS(4)*1.00000*( ALGBRC(106)*ALGBRC(104)*STATES(32) - CONSTS(48)*ALGBRC(108)*(1.00000 - STATES(32))) ALGBRC(100) = CONSTS(46)*CONSTS(4)*1.00000*( ALGBRC(98)*ALGBRC(96)*STATES(31) - ALGBRC(99)*CONSTS(96)*(1.00000 - STATES(31))) ALGBRC(111) = ALGBRC(50)+ALGBRC(55)+ALGBRC(59)+ALGBRC(81)+ALGBRC(82)+ALGBRC(95)+ALGBRC(109)+ALGBRC(100)+ALGBRC(66)+ALGBRC(70)+ALGBRC(83) RATES(1) = - ALGBRC(111)/CONSTS(4) ALGBRC(112) = ALGBRC(45)+ALGBRC(53)+ALGBRC(61)+ALGBRC(68)+ALGBRC(85)+ALGBRC(90)+ 3.00000*ALGBRC(109)+ 3.00000*ALGBRC(100) RATES(2) = - ALGBRC(112)/( CONSTS(3)*CONSTS(8)) ALGBRC(113) = (ALGBRC(81)+ALGBRC(82)+ALGBRC(93)+ALGBRC(83)+ALGBRC(48)+ALGBRC(54)+ALGBRC(64)+ALGBRC(69)+ALGBRC(84)+ALGBRC(89)+ALGBRC(87)) - 2.00000*ALGBRC(109) RATES(3) = - ALGBRC(113)/( CONSTS(3)*CONSTS(8)) ALGBRC(110) = 1.00000/(1.00000+ CONSTS(54)/CONSTS(7) ** 1.12000*(1.00000+ALGBRC(107)/CONSTS(52) ** 1.06000)) ALGBRC(115) = CONSTS(48)*ALGBRC(108)+ CONSTS(50)*ALGBRC(110) ALGBRC(105) = 1.00000/(1.00000+ CONSTS(53)/STATES(3) ** 1.12000*(1.00000+STATES(2)/CONSTS(51) ** 1.06000)) ALGBRC(117) = ALGBRC(106)*ALGBRC(104)+ CONSTS(49)*ALGBRC(105) RATES(32) = ALGBRC(115)*(1.00000 - STATES(32)) - ALGBRC(117)*STATES(32) ALGBRC(116) = 1.00000/(1.00000+CONSTS(61)/ALGBRC(28)) ALGBRC(114) = 1.00000/(1.00000+CONSTS(60)/STATES(33)) ALGBRC(120) = 1.00000/(1.00000+CONSTS(62)/STATES(5)) ALGBRC(121) = CONSTS(63)*1.00000*( CONSTS(57)*ALGBRC(114)*STATES(34) - ALGBRC(120)*ALGBRC(116)*(1.00000 - STATES(34))) ALGBRC(123) = - 0.400000*STATES(39)*CONSTS(73) RATES(5) = (( CONSTS(11)*(CONSTS(12) - STATES(5))+ALGBRC(123)) - ALGBRC(109)/( CONSTS(3)*CONSTS(8)))+ALGBRC(121)/( 4.00000*CONSTS(3)*CONSTS(8)) ALGBRC(119) = 1.00000 - ALGBRC(116) ALGBRC(122) = ALGBRC(120)*ALGBRC(116)+ CONSTS(59)*ALGBRC(119) ALGBRC(118) = 1.00000 - ALGBRC(114) ALGBRC(124) = CONSTS(57)*ALGBRC(114)+ CONSTS(58)*ALGBRC(118) RATES(34) = ALGBRC(122)*(1.00000 - STATES(34)) - ALGBRC(124)*STATES(34) ALGBRC(57) = ALGBRC(56)*ALGBRC(51) ALGBRC(125) = 280000.*ALGBRC(28)/1.00000 ** 2.00000+ CONSTS(66)*ALGBRC(57) ALGBRC(127) = (CONSTS(71) - STATES(37))+CONSTS(72) ALGBRC(128) = CONSTS(72)*STATES(37) ALGBRC(129) = ( (ALGBRC(127) ** 2.00000+ 4.00000*ALGBRC(128)) ** (1.0 / 2) - ALGBRC(127))/2.00000 ALGBRC(130) = 0.0800000/(1.00000+0.360000/ALGBRC(129)) RATES(35) = STATES(36)*ALGBRC(125) - STATES(35)*ALGBRC(130) ALGBRC(131) = 0.000377000*ALGBRC(129)/1.00000 ** 2.00000 RATES(36) = ALGBRC(131)*(1.00000 - (STATES(35)+STATES(36))) - (ALGBRC(125)+CONSTS(65))*STATES(36) ALGBRC(132) = CONSTS(64)*(ALGBRC(129) - ALGBRC(28))*STATES(35) ALGBRC(133) = CONSTS(67)*(STATES(33) - ALGBRC(129)) RATES(37) = (ALGBRC(133) - ALGBRC(132))/( 2.00000*CONSTS(3)*CONSTS(69)) ALGBRC(126) = CONSTS(68)*(STATES(33) - ALGBRC(28)) RATES(33) = ((- ALGBRC(121) - ALGBRC(133)) - ALGBRC(126))/( 2.00000*CONSTS(3)*CONSTS(70)) ALGBRC(140) = CONSTS(75)*STATES(38)*CONSTS(97) - CONSTS(80)*STATES(39) ALGBRC(136) = ((1.00000 - STATES(38)) - STATES(39)) - STATES(40) ALGBRC(138) = CONSTS(74)*ALGBRC(28)*ALGBRC(136) - CONSTS(79)*STATES(38) RATES(38) = ALGBRC(138) - ALGBRC(140) ALGBRC(134) = CONSTS(83) - STATES(41) RATES(41) = CONSTS(88)*(ALGBRC(134) - CONSTS(89)) ALGBRC(142) = CONSTS(76)*STATES(39) - CONSTS(81)*STATES(40)*ALGBRC(28) ALGBRC(145) = CONSTS(78)*RATES(41) ** 2.00000*STATES(39) RATES(39) = (ALGBRC(140) - ALGBRC(142)) - ALGBRC(145) ALGBRC(143) = CONSTS(77)*STATES(40) ALGBRC(144) = CONSTS(78)*RATES(41) ** 2.00000*STATES(40) RATES(40) = (ALGBRC(142) - ALGBRC(143)) - ALGBRC(144) ALGBRC(102) = (ALGBRC(52)+ALGBRC(59)+ALGBRC(94)) - 2.00000*ALGBRC(100) ALGBRC(146) = CONSTS(73)*((ALGBRC(145)+ALGBRC(142)) - ALGBRC(138)) RATES(4) = - (((ALGBRC(102) - ALGBRC(121)) - ALGBRC(132)) - ALGBRC(126))/( 2.00000*CONSTS(3)*CONSTS(8))+ALGBRC(146) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(14) = 0.0100000*EXP( 0.0133000*(STATES(1)+40.0000)) ALGBRC(3) = 1.00000/( 9.00000e+09*EXP(STATES(1)/5.00000)+ 8000.00*EXP(STATES(1)/100.000)) ALGBRC(17) = 1.00000/( 0.0140000*EXP(- STATES(1)/5.00000)+ 4000.00*EXP(- STATES(1)/100.000)) ALGBRC(5) = 1.00000/( 250000.*EXP(STATES(1)/9.00000)+ 58.0000*EXP(STATES(1)/65.0000)) ALGBRC(19) = 1.00000/( 1800.00*EXP(- STATES(1)/14.0000)+ 66.0000*EXP(- STATES(1)/65.0000)) ALGBRC(6) = 1.00000/( 0.0190000*EXP(- STATES(1)/5.60000)+ 0.820000*EXP(- STATES(1)/250.000)) ALGBRC(20) = 1.00000/( 40.0000*EXP(STATES(1)/6.30000)+ 1.50000*EXP(STATES(1)/10000.0)) ALGBRC(7) = 1.00000/( 62000.0*EXP(STATES(1)/10.1000)+ 30.0000*EXP(STATES(1)/3000.00)) ALGBRC(21) = 1.00000/( 0.000600000*EXP(- STATES(1)/6.70000)+ 1.20000*EXP(- STATES(1)/25.0000)) ALGBRC(8) = 1.00000/( 0.150000*EXP(- STATES(1)/11.0000)+ 0.200000*EXP(- STATES(1)/700.000)) ALGBRC(22) = 1.00000/( 16.0000*EXP(STATES(1)/8.00000)+ 15.0000*EXP(STATES(1)/50.0000)) ALGBRC(9) = 1.00000/( 400000.*EXP(STATES(1)/9.00000)+ 60.0000*EXP(STATES(1)/65.0000)) ALGBRC(23) = 1.00000/( 700.000*EXP(- STATES(1)/14.0000)+ 60.0000*EXP(- STATES(1)/65.0000)) ALGBRC(11) = 1.00000/( 20.0000*EXP(- STATES(1)/11.5000)+ 5.00000*EXP(- STATES(1)/300.000)) ALGBRC(25) = 1.00000/( 160.000*EXP(STATES(1)/28.0000)+ 200.000*EXP(STATES(1)/1000.00))+1.00000/( 2500.00*EXP(STATES(1)/20.0000)) ALGBRC(12) = 1.00000/( 200.000*EXP(- STATES(1)/13.0000)+ 20.0000*EXP(- STATES(1)/300.000)) ALGBRC(26) = 1.00000/( 1600.00*EXP(STATES(1)/28.0000)+ 2000.00*EXP(STATES(1)/1000.00))+1.00000/( 10000.0*EXP(STATES(1)/20.0000)) ALGBRC(13) = 1.00000/( 10.0000*EXP(STATES(1)/17.0000)+ 2.50000*EXP(STATES(1)/300.000)) ALGBRC(27) = 1.00000/( 0.350000*EXP(- STATES(1)/17.0000)+ 2.00000*EXP(- STATES(1)/150.000)) ALGBRC(18) = 1.00000/( 0.270000*EXP(- STATES(1)/5.90000)+ 1.50000*EXP(- STATES(1)/65.0000)) ALGBRC(30) = 1.00000/( 480.000*EXP(STATES(1)/7.00000)+ 2.20000*EXP(STATES(1)/65.0000)) ALGBRC(24) = 1.00000/( 3500.00*EXP(STATES(1)/16.8000)+ 0.300000*EXP(STATES(1)/400.000)) ALGBRC(31) = 1.00000/( 4.00000*EXP(- STATES(1)/14.0000)+ 2.00000*EXP(- STATES(1)/400.000)) ALGBRC(10) = (((1.00000 - STATES(22)) - STATES(23)) - STATES(24)) - STATES(25) ALGBRC(16) = 1.00000/( 0.102700*EXP(- STATES(1)/8.00000)+ 0.250000*EXP(- STATES(1)/50.0000)) ALGBRC(29) = 1.00000/( 26.0000*EXP(STATES(1)/17.0000)+ 0.0200000*EXP(STATES(1)/800.000)) ALGBRC(33) = 1.00000/( 0.800000*EXP(- STATES(1)/400.000)) ALGBRC(4) = ((1.00000 - STATES(10)) - STATES(14)) - STATES(15) ALGBRC(34) = 1.00000/( 0.00180000*EXP(- STATES(1)/7.40000)+ 2.00000*EXP(- STATES(1)/100.000)) ALGBRC(39) = 1.00000/( 2.20000e+06*EXP(STATES(1)/7.40000)+ 11.0000*EXP(STATES(1)/100.000)) ALGBRC(35) = 1.00000/( 3100.00*EXP(STATES(1)/13.0000)+ 700.000*EXP(STATES(1)/70.0000)) ALGBRC(1) = (CONSTS(9) - STATES(4))+CONSTS(10) ALGBRC(15) = CONSTS(10)*STATES(4) ALGBRC(28) = ( (ALGBRC(1) ** 2.00000+ 4.00000*ALGBRC(15)) ** (1.0 / 2) - ALGBRC(1))/2.00000 ALGBRC(40) = 1.00000/( 95.0000*EXP(- STATES(1)/10.0000)+ 50.0000*EXP(- STATES(1)/700.000))+( ALGBRC(28)*2.50000)/(1.00000+EXP(- STATES(1)/5.00000)) ALGBRC(36) = 1.00000/( 4.50000e+07*EXP(STATES(1)/8.00000)+ 500.000*EXP(STATES(1)/200.000)) ALGBRC(41) = 1.00000/( 10.5000*EXP(- STATES(1)/16.4000)+ 0.400000*EXP(- STATES(1)/400.000)) ALGBRC(2) = ((1.00000 - STATES(8)) - STATES(6)) - STATES(9) ALGBRC(43) = 1.00000/( 0.000102700*EXP(- STATES(1)/8.00000)+ 5.00000*EXP(- STATES(1)/400.000)) ALGBRC(38) = 1.00000/( 1300.00*EXP(STATES(1)/20.0000)+ 0.0400000*EXP(STATES(1)/800.000)) ALGBRC(44) = 0.0400000/(1.00000+( CONSTS(17)*ALGBRC(30)*ALGBRC(34))/( CONSTS(16)*ALGBRC(18)*ALGBRC(39))) ALGBRC(47) = 0.0400000 - ALGBRC(44) ALGBRC(46) = 0.0100000/(1.00000+( CONSTS(14)*ALGBRC(29)*ALGBRC(43))/( ALGBRC(33)*ALGBRC(16)*ALGBRC(38))) ALGBRC(49) = 0.0100000 - ALGBRC(46) ALGBRC(37) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(6), ( (( 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(ALGBRC(28) - CONSTS(6)*EXP(( - 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(56) = 0.0676000*ALGBRC(37) ALGBRC(58) = ALGBRC(28) - 0.300000*ALGBRC(56) ALGBRC(60) = ALGBRC(58)*STATES(10) ALGBRC(62) = ALGBRC(60)+ ALGBRC(28)*(1.00000 - STATES(10)) ALGBRC(65) = CONSTS(22)*ALGBRC(62) ALGBRC(67) = ((1.00000 - STATES(16)) - STATES(11)) - STATES(12) ALGBRC(71) = (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(7)/STATES(3)) ALGBRC(77) = 1.00000/( 8000.00*EXP(((STATES(1) - ALGBRC(71)) - 97.0000)/8.50000)+ 7.00000*EXP(((STATES(1) - ALGBRC(71)) - 97.0000)/300.000)) ALGBRC(72) = ( 0.750000*EXP( 0.0350000*((STATES(1) - ALGBRC(71)) - 10.0000)))/(1.00000+EXP( 0.0150000*((STATES(1) - ALGBRC(71)) - 140.000))) ALGBRC(73) = ( 3.00000*EXP( - 0.0480000*((STATES(1) - ALGBRC(71)) - 10.0000))*(1.00000+EXP( 0.0640000*((STATES(1) - ALGBRC(71)) - 38.0000))))/(1.00000+EXP( 0.0300000*((STATES(1) - ALGBRC(71)) - 70.0000))) ALGBRC(75) = ALGBRC(73)/(ALGBRC(72)+ALGBRC(73)) ALGBRC(79) = ( ALGBRC(75) ** 4.00000*1.00000)/( 0.000140000*EXP(- ((STATES(1) - ALGBRC(71)) - 97.0000)/9.10000)+ 0.200000*EXP(- ((STATES(1) - ALGBRC(71)) - 97.0000)/500.000)) ALGBRC(99) = 1.00000*EXP(( (CONSTS(47) - 1.00000)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(101) = ALGBRC(99)*CONSTS(96)+ CONSTS(41)*CONSTS(99) ALGBRC(96) = 1.00000/(1.00000+ CONSTS(42)/STATES(2) ** 3.00000*(1.00000+ALGBRC(28)/CONSTS(44))) ALGBRC(97) = 1.00000/(1.00000+ (CONSTS(44)/ALGBRC(28))*(1.00000+STATES(2)/CONSTS(42) ** 3.00000)) ALGBRC(98) = 1.00000*EXP(( CONSTS(47)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(103) = ALGBRC(98)*ALGBRC(96)+ CONSTS(40)*ALGBRC(97) ALGBRC(32) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(5), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(2) - CONSTS(5)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(85) = CONSTS(35)*ALGBRC(32) ALGBRC(42) = TERNRY(STATES(1).EQ.0.00000, STATES(3), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(3) - CONSTS(7)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))))) ALGBRC(84) = 0.400000*CONSTS(35)*ALGBRC(42) ALGBRC(86) = ALGBRC(84)+ALGBRC(85) ALGBRC(94) = CONSTS(39)*ALGBRC(37) ALGBRC(87) = TERNRY(STATES(1).EQ.- 3.00000, CONSTS(94)*ALGBRC(42)*13.0077, ( CONSTS(94)*ALGBRC(42)*(STATES(1)+3.00000))/(1.00000 - EXP(- (STATES(1)+3.00000)/13.0000))) ALGBRC(88) = 1.00000/(1.00000+0.00120000/ALGBRC(28) ** 3.00000) ALGBRC(90) = CONSTS(36)*ALGBRC(32)*ALGBRC(88) ALGBRC(89) = CONSTS(36)*ALGBRC(42)*ALGBRC(88) ALGBRC(91) = ALGBRC(89)+ALGBRC(90) ALGBRC(92) = 0.800000/(1.00000+STATES(5)/0.100000 ** 2.00000) ALGBRC(93) = CONSTS(95)*(STATES(1) - ALGBRC(71))*ALGBRC(92) ALGBRC(95) = ALGBRC(86)+ALGBRC(94)+ALGBRC(87)+ALGBRC(91)+ALGBRC(93) ALGBRC(45) = CONSTS(13)*ALGBRC(32)*STATES(6)*STATES(7) ALGBRC(48) = 0.100000*CONSTS(13)*ALGBRC(42)*STATES(6)*STATES(7) ALGBRC(50) = ALGBRC(45)+ALGBRC(48) ALGBRC(51) = ( STATES(10)*(STATES(11)+STATES(12))*STATES(13))/(1.00000+1.40000/STATES(5) ** 3.00000) ALGBRC(53) = 1.85000e-05*CONSTS(15)*ALGBRC(32)*ALGBRC(51) ALGBRC(54) = 0.000365000*CONSTS(15)*ALGBRC(42)*ALGBRC(51) ALGBRC(52) = CONSTS(15)*ALGBRC(37)*ALGBRC(51) ALGBRC(55) = ALGBRC(53)+ALGBRC(52)+ALGBRC(54) ALGBRC(59) = CONSTS(25)*ALGBRC(37)*STATES(17)*STATES(18) ALGBRC(74) = ALGBRC(72)/(ALGBRC(72)+ALGBRC(73)) ALGBRC(76) = 2.00000*ALGBRC(75) ** 2.00000*ALGBRC(74) ** 2.00000 ALGBRC(78) = (8.00000/3.00000)*ALGBRC(75) ** 3.00000*ALGBRC(74) ALGBRC(80) = ALGBRC(75) ** 4.00000 ALGBRC(81) = CONSTS(91)*(STATES(1) - ALGBRC(71))*(ALGBRC(80)+ALGBRC(78)+ALGBRC(76))*STATES(26) ALGBRC(82) = CONSTS(92)*(STATES(1) - ALGBRC(71))*( 0.600000*STATES(27)+ 0.400000*STATES(28))*STATES(29) ALGBRC(61) = CONSTS(26)*ALGBRC(32)*STATES(19)*STATES(20)*STATES(21) ALGBRC(64) = CONSTS(27)*ALGBRC(42)*STATES(19)*STATES(20)*STATES(21) ALGBRC(66) = ALGBRC(61)+ALGBRC(64) ALGBRC(68) = CONSTS(28)*ALGBRC(32)*((1.00000 - STATES(22)) - STATES(23)) ALGBRC(69) = CONSTS(29)*ALGBRC(42)*((1.00000 - STATES(22)) - STATES(23)) ALGBRC(70) = ALGBRC(68)+ALGBRC(69) ALGBRC(83) = (( CONSTS(32)*(STATES(1) - ALGBRC(71)))/(1.00000+EXP((STATES(1)+20.0000)/20.0000)))*STATES(30) ALGBRC(107) = CONSTS(5)*EXP(( - 0.820000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2))) ALGBRC(108) = 1.00000/(1.00000+ CONSTS(52)/ALGBRC(107) ** 1.06000*(1.00000+CONSTS(7)/CONSTS(54) ** 1.12000)) ALGBRC(104) = 1.00000/(1.00000+ CONSTS(51)/STATES(2) ** 1.06000*(1.00000+STATES(3)/CONSTS(53) ** 1.12000)) ALGBRC(106) = 0.370000/(1.00000+CONSTS(55)/STATES(5)) ALGBRC(109) = CONSTS(56)*CONSTS(4)*1.00000*( ALGBRC(106)*ALGBRC(104)*STATES(32) - CONSTS(48)*ALGBRC(108)*(1.00000 - STATES(32))) ALGBRC(100) = CONSTS(46)*CONSTS(4)*1.00000*( ALGBRC(98)*ALGBRC(96)*STATES(31) - ALGBRC(99)*CONSTS(96)*(1.00000 - STATES(31))) ALGBRC(111) = ALGBRC(50)+ALGBRC(55)+ALGBRC(59)+ALGBRC(81)+ALGBRC(82)+ALGBRC(95)+ALGBRC(109)+ALGBRC(100)+ALGBRC(66)+ALGBRC(70)+ALGBRC(83) ALGBRC(112) = ALGBRC(45)+ALGBRC(53)+ALGBRC(61)+ALGBRC(68)+ALGBRC(85)+ALGBRC(90)+ 3.00000*ALGBRC(109)+ 3.00000*ALGBRC(100) ALGBRC(113) = (ALGBRC(81)+ALGBRC(82)+ALGBRC(93)+ALGBRC(83)+ALGBRC(48)+ALGBRC(54)+ALGBRC(64)+ALGBRC(69)+ALGBRC(84)+ALGBRC(89)+ALGBRC(87)) - 2.00000*ALGBRC(109) ALGBRC(110) = 1.00000/(1.00000+ CONSTS(54)/CONSTS(7) ** 1.12000*(1.00000+ALGBRC(107)/CONSTS(52) ** 1.06000)) ALGBRC(115) = CONSTS(48)*ALGBRC(108)+ CONSTS(50)*ALGBRC(110) ALGBRC(105) = 1.00000/(1.00000+ CONSTS(53)/STATES(3) ** 1.12000*(1.00000+STATES(2)/CONSTS(51) ** 1.06000)) ALGBRC(117) = ALGBRC(106)*ALGBRC(104)+ CONSTS(49)*ALGBRC(105) ALGBRC(116) = 1.00000/(1.00000+CONSTS(61)/ALGBRC(28)) ALGBRC(114) = 1.00000/(1.00000+CONSTS(60)/STATES(33)) ALGBRC(120) = 1.00000/(1.00000+CONSTS(62)/STATES(5)) ALGBRC(121) = CONSTS(63)*1.00000*( CONSTS(57)*ALGBRC(114)*STATES(34) - ALGBRC(120)*ALGBRC(116)*(1.00000 - STATES(34))) ALGBRC(123) = - 0.400000*STATES(39)*CONSTS(73) ALGBRC(119) = 1.00000 - ALGBRC(116) ALGBRC(122) = ALGBRC(120)*ALGBRC(116)+ CONSTS(59)*ALGBRC(119) ALGBRC(118) = 1.00000 - ALGBRC(114) ALGBRC(124) = CONSTS(57)*ALGBRC(114)+ CONSTS(58)*ALGBRC(118) ALGBRC(57) = ALGBRC(56)*ALGBRC(51) ALGBRC(125) = 280000.*ALGBRC(28)/1.00000 ** 2.00000+ CONSTS(66)*ALGBRC(57) ALGBRC(127) = (CONSTS(71) - STATES(37))+CONSTS(72) ALGBRC(128) = CONSTS(72)*STATES(37) ALGBRC(129) = ( (ALGBRC(127) ** 2.00000+ 4.00000*ALGBRC(128)) ** (1.0 / 2) - ALGBRC(127))/2.00000 ALGBRC(130) = 0.0800000/(1.00000+0.360000/ALGBRC(129)) ALGBRC(131) = 0.000377000*ALGBRC(129)/1.00000 ** 2.00000 ALGBRC(132) = CONSTS(64)*(ALGBRC(129) - ALGBRC(28))*STATES(35) ALGBRC(133) = CONSTS(67)*(STATES(33) - ALGBRC(129)) ALGBRC(126) = CONSTS(68)*(STATES(33) - ALGBRC(28)) ALGBRC(140) = CONSTS(75)*STATES(38)*CONSTS(97) - CONSTS(80)*STATES(39) ALGBRC(136) = ((1.00000 - STATES(38)) - STATES(39)) - STATES(40) ALGBRC(138) = CONSTS(74)*ALGBRC(28)*ALGBRC(136) - CONSTS(79)*STATES(38) ALGBRC(134) = CONSTS(83) - STATES(41) ALGBRC(142) = CONSTS(76)*STATES(39) - CONSTS(81)*STATES(40)*ALGBRC(28) ALGBRC(145) = CONSTS(78)*RATES(41) ** 2.00000*STATES(39) ALGBRC(143) = CONSTS(77)*STATES(40) ALGBRC(144) = CONSTS(78)*RATES(41) ** 2.00000*STATES(40) ALGBRC(102) = (ALGBRC(52)+ALGBRC(59)+ALGBRC(94)) - 2.00000*ALGBRC(100) ALGBRC(146) = CONSTS(73)*((ALGBRC(145)+ALGBRC(142)) - ALGBRC(138)) ALGBRC(63) = CONSTS(24)*ALGBRC(60) ALGBRC(135) = CONSTS(73)*(STATES(39)+STATES(40)) ALGBRC(137) = CONSTS(87)*ALGBRC(135) ALGBRC(139) = ALGBRC(137)*ALGBRC(134) ALGBRC(141) = - CONSTS(98)+ALGBRC(139) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END