C
C There are a total of 141 entries in the algebraic variable array.
C There are a total of 37 entries in each of the rate and state variable arrays.
C There are a total of 101 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(1) is i_ext in component membrane (picoA).
C ALGBRC(106) is i_tot in component membrane (picoA).
C ALGBRC(90) is i_I in component membrane (picoA).
C ALGBRC(46) is i_Na in component sodium_current (picoA).
C ALGBRC(51) is i_Ca_L in component L_type_Ca_channel (picoA).
C ALGBRC(55) is i_Ca_T in component T_type_Ca_channel (picoA).
C ALGBRC(71) is i_K1 in component time_independent_potassium_current (picoA).
C ALGBRC(72) is i_Kr in component rapid_time_dependent_potassium_current (picoA).
C ALGBRC(75) is i_Ks in component slow_time_dependent_potassium_current (picoA).
C ALGBRC(78) is i_to in component transient_outward_current (picoA).
C ALGBRC(104) is i_NaK in component sodium_potassium_pump (picoA).
C ALGBRC(95) is i_NaCa in component sodium_calcium_exchanger (picoA).
C ALGBRC(81) is i_bNSC in component background_NSC_current (picoA).
C ALGBRC(89) is i_Cab in component background_Cab_current (picoA).
C ALGBRC(82) is i_Kpl in component background_Kpl_current (picoA).
C ALGBRC(86) is i_lCa in component background_lCa_current (picoA).
C ALGBRC(88) is i_KATP in component background_KATP_current (picoA).
C CONSTS(5) is stim_start in component membrane (millisecond).
C CONSTS(6) is stim_end in component membrane (millisecond).
C CONSTS(7) is stim_period in component membrane (millisecond).
C CONSTS(8) is stim_duration in component membrane (millisecond).
C CONSTS(9) is stim_amplitude in component membrane (picoA).
C CONSTS(10) is Nao in component external_ion_concentrations (millimolar).
C CONSTS(11) is Cao in component external_ion_concentrations (millimolar).
C CONSTS(12) is Ko in component external_ion_concentrations (millimolar).
C STATES(2) is Nai in component internal_ion_concentrations (millimolar).
C ALGBRC(30) is Cai in component internal_ion_concentrations (millimolar).
C STATES(3) is Ki in component internal_ion_concentrations (millimolar).
C CONSTS(13) is Vi in component internal_ion_concentrations (micrometre3).
C ALGBRC(107) is i_net_Na in component internal_ion_concentrations (picoA).
C ALGBRC(108) is i_net_K in component internal_ion_concentrations (picoA).
C ALGBRC(97) is i_net_Ca in component internal_ion_concentrations (picoA).
C ALGBRC(43) is i_Na_Na in component sodium_current (picoA).
C ALGBRC(49) is i_CaL_Na in component L_type_Ca_channel (picoA).
C ALGBRC(80) is i_bNSC_Na in component background_NSC_current (picoA).
C ALGBRC(85) is i_lCa_Na in component background_lCa_current (picoA).
C ALGBRC(76) is i_to_K in component transient_outward_current (picoA).
C ALGBRC(77) is i_to_Na in component transient_outward_current (picoA).
C ALGBRC(73) is i_Ks_K in component slow_time_dependent_potassium_current (picoA).
C ALGBRC(74) is i_Ks_Na in component slow_time_dependent_potassium_current (picoA).
C ALGBRC(45) is i_Na_K in component sodium_current (picoA).
C ALGBRC(50) is i_CaL_K in component L_type_Ca_channel (picoA).
C ALGBRC(79) is i_bNSC_K in component background_NSC_current (picoA).
C ALGBRC(84) is i_lCa_K in component background_lCa_current (picoA).
C ALGBRC(48) is i_CaL_Ca in component L_type_Ca_channel (picoA).
C ALGBRC(127) is i_RyR in component RyR_channel (picoA).
C ALGBRC(116) is i_SR_U in component SR_calcium_pump (picoA).
C ALGBRC(121) is i_SR_L in component SR_L_current (picoA).
C ALGBRC(141) is dCaidt in component NL_model (millimolar_per_millisecond).
C CONSTS(14) is CMDN_max in component internal_ion_concentrations (millimolar).
C CONSTS(15) is K_mCMDN in component internal_ion_concentrations (millimolar).
C STATES(4) is Ca_Total in component internal_ion_concentrations (millimolar).
C ALGBRC(14) is b1 in component internal_ion_concentrations (millimolar).
C ALGBRC(27) is c1 in component internal_ion_concentrations (millimolar2).
C ALGBRC(33) is CF_Na in component constant_field_equations (millimolar).
C ALGBRC(37) is CF_Ca in component constant_field_equations (millimolar).
C ALGBRC(40) is CF_K in component constant_field_equations (millimolar).
C STATES(5) is ATPi in component ATP_production (millimolar).
C ALGBRC(118) is dATPdt in component NL_model (millimolar_per_millisecond).
C CONSTS(16) is ProducingRate_Max in component ATP_production (per_millisecond).
C CONSTS(17) is Adenosine_Total in component ATP_production (millimolar).
C CONSTS(18) 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(15) is k_RP_AP in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(28) is k_AP_RP in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(38) is k_RI_AI in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(34) is k_AI_RI in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(31) is k_AP_AI in component sodium_current_voltage_dependent_gate (per_millisecond).
C CONSTS(19) is k_AI_AP in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(41) is k_RP_RI in component sodium_current_voltage_dependent_gate (per_millisecond).
C ALGBRC(44) 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(16) is beta_y in component sodium_current_ultra_slow_gate (per_millisecond).
C ALGBRC(47) is p_open_CaL in component L_type_Ca_channel (dimensionless).
C ALGBRC(53) is CaDiadic in component L_type_Ca_channel_Ca_dependent_gate (picoA).
C CONSTS(20) 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(17) is k_RP_AP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(29) is k_AP_RP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(32) is k_RI_AI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(35) is k_AI_RI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C CONSTS(21) is k_AP_AI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C CONSTS(22) is k_AI_AP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(39) is k_RP_RI in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(42) is k_RI_RP in component L_type_Ca_channel_voltage_dependent_gate (per_millisecond).
C ALGBRC(52) is iCaL in component L_type_Ca_channel_Ca_dependent_gate (picoA).
C ALGBRC(54) is Cacm in component L_type_Ca_channel_Ca_dependent_gate (millimolar).
C ALGBRC(63) 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(23) is k_CCa_UCa in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(24) is k_UCa_CCa in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(25) is k_C_U in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(26) is k_U_C in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(93) is k_UCa_U in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(27) is k_U_UCa in component L_type_Ca_channel_Ca_dependent_gate (per_millimolar_millisecond).
C CONSTS(28) is k_CCa_C in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C CONSTS(29) is k_C_CCa in component L_type_Ca_channel_Ca_dependent_gate (per_millimolar_millisecond).
C ALGBRC(56) is CaEffC in component L_type_Ca_channel_Ca_dependent_gate (millimolar).
C ALGBRC(58) is CaEffU in component L_type_Ca_channel_Ca_dependent_gate (millimolar).
C ALGBRC(61) is k_UUCa_Ca in component L_type_Ca_channel_Ca_dependent_gate (per_millisecond).
C ALGBRC(59) 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(18) is beta_y in component L_type_Ca_channel_ultra_slow_gate (per_millisecond).
C CONSTS(30) 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(19) 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(20) is beta_y2 in component T_type_Ca_channel_y2_gate (per_millisecond).
C ALGBRC(57) is E_K in component time_independent_potassium_current (millivolt).
C CONSTS(94) is g_K1 in component time_independent_potassium_current (nanoS).
C CONSTS(31) is P_K1_0 in component time_independent_potassium_current (nanoS_per_picoF).
C ALGBRC(65) is fO in component time_independent_potassium_current (dimensionless).
C ALGBRC(66) is fO2 in component time_independent_potassium_current (dimensionless).
C ALGBRC(68) is fO3 in component time_independent_potassium_current (dimensionless).
C ALGBRC(70) is fO4 in component time_independent_potassium_current (dimensionless).
C ALGBRC(64) is fB in component time_independent_potassium_current (dimensionless).
C ALGBRC(60) is mu in component time_independent_potassium_current (per_millisecond).
C ALGBRC(62) is lambda in component time_independent_potassium_current (per_millisecond).
C STATES(19) is y in component time_independent_potassium_current_y_gate (dimensionless).
C ALGBRC(67) is alpha_y in component time_independent_potassium_current_y_gate (per_millisecond).
C ALGBRC(69) is beta_y in component time_independent_potassium_current_y_gate (per_millisecond).
C CONSTS(95) is g_Kr in component rapid_time_dependent_potassium_current (nanoS).
C CONSTS(32) is P_Kr in component rapid_time_dependent_potassium_current (nanoS_per_picoF).
C STATES(20) is y1 in component rapid_time_dependent_potassium_current_y1_gate (dimensionless).
C STATES(21) is y2 in component rapid_time_dependent_potassium_current_y2_gate (dimensionless).
C STATES(22) is y3 in component rapid_time_dependent_potassium_current_y3_gate (dimensionless).
C ALGBRC(8) is alpha_y1 in component rapid_time_dependent_potassium_current_y1_gate (per_millisecond).
C ALGBRC(21) is beta_y1 in component rapid_time_dependent_potassium_current_y1_gate (per_millisecond).
C ALGBRC(9) is alpha_y2 in component rapid_time_dependent_potassium_current_y2_gate (per_millisecond).
C ALGBRC(22) is beta_y2 in component rapid_time_dependent_potassium_current_y2_gate (per_millisecond).
C ALGBRC(10) is alpha_y3 in component rapid_time_dependent_potassium_current_y3_gate (per_millisecond).
C ALGBRC(23) is beta_y3 in component rapid_time_dependent_potassium_current_y3_gate (per_millisecond).
C STATES(23) is y1 in component slow_time_dependent_potassium_current_y1_gate (dimensionless).
C STATES(24) is y2 in component slow_time_dependent_potassium_current_y2_gate (dimensionless).
C CONSTS(33) is P_Ks_K in component slow_time_dependent_potassium_current (picoA_per_millimolar).
C CONSTS(34) is P_Ks_Na in component slow_time_dependent_potassium_current (picoA_per_millimolar).
C ALGBRC(11) is alpha_y1 in component slow_time_dependent_potassium_current_y1_gate (per_millisecond).
C ALGBRC(24) is beta_y1 in component slow_time_dependent_potassium_current_y1_gate (per_millisecond).
C ALGBRC(36) is alpha_y2 in component slow_time_dependent_potassium_current_y2_gate (per_millisecond).
C CONSTS(35) is beta_y2 in component slow_time_dependent_potassium_current_y2_gate (per_millisecond).
C STATES(25) is y1 in component transient_outward_current_y1_gate (dimensionless).
C STATES(26) is y2 in component transient_outward_current_y2_gate (dimensionless).
C CONSTS(36) is P_to_K in component transient_outward_current (picoA_per_millimolar).
C CONSTS(37) is P_to_Na in component transient_outward_current (picoA_per_millimolar).
C ALGBRC(12) is alpha_y1 in component transient_outward_current_y1_gate (per_millisecond).
C ALGBRC(25) is beta_y1 in component transient_outward_current_y1_gate (per_millisecond).
C ALGBRC(13) is alpha_y2 in component transient_outward_current_y2_gate (per_millisecond).
C ALGBRC(26) is beta_y2 in component transient_outward_current_y2_gate (per_millisecond).
C CONSTS(38) is P_bNSC in component background_NSC_current (picoA_per_millimolar).
C CONSTS(96) is P_Kpl in component background_Kpl_current (nanoS_per_millimolar).
C CONSTS(39) is P_lCa in component background_lCa_current (picoA_per_millimolar).
C ALGBRC(83) is p_open in component background_lCa_current (dimensionless).
C ALGBRC(87) is p_open in component background_KATP_current (dimensionless).
C CONSTS(97) is gamma in component background_KATP_current (nanoS).
C CONSTS(40) is P_KATP in component background_KATP_current (nanoS_per_picoF).
C CONSTS(41) is N in component background_KATP_current (picoF).
C CONSTS(42) is P_Cab in component background_Cab_current (picoA_per_millimolar).
C CONSTS(98) is p_E2Na in component sodium_calcium_exchanger (dimensionless).
C ALGBRC(91) is p_E1Na in component sodium_calcium_exchanger (dimensionless).
C ALGBRC(92) is p_E1Ca in component sodium_calcium_exchanger (dimensionless).
C CONSTS(101) is p_E2Ca in component sodium_calcium_exchanger (dimensionless).
C ALGBRC(93) is k1 in component sodium_calcium_exchanger (per_millisecond).
C ALGBRC(94) is k2 in component sodium_calcium_exchanger (per_millisecond).
C CONSTS(43) is k3 in component sodium_calcium_exchanger (per_millisecond).
C CONSTS(44) is k4 in component sodium_calcium_exchanger (per_millisecond).
C CONSTS(45) is Km_Nai in component sodium_calcium_exchanger (millimolar).
C CONSTS(46) is Km_Nao in component sodium_calcium_exchanger (millimolar).
C CONSTS(47) is Km_Cai in component sodium_calcium_exchanger (millimolar).
C CONSTS(48) is Km_Cao in component sodium_calcium_exchanger (millimolar).
C STATES(27) is y in component sodium_calcium_exchanger_y_gate (dimensionless).
C CONSTS(49) is P_NaCa in component sodium_calcium_exchanger (picoA_per_picoF).
C CONSTS(50) is Partition in component sodium_calcium_exchanger (dimensionless).
C ALGBRC(96) is alpha_y in component sodium_calcium_exchanger_y_gate (per_millisecond).
C ALGBRC(98) is beta_y in component sodium_calcium_exchanger_y_gate (per_millisecond).
C ALGBRC(103) is p_E2Na in component sodium_potassium_pump (dimensionless).
C ALGBRC(99) is p_E1Na in component sodium_potassium_pump (dimensionless).
C ALGBRC(100) is p_E1K in component sodium_potassium_pump (dimensionless).
C ALGBRC(105) is p_E2K in component sodium_potassium_pump (dimensionless).
C ALGBRC(101) is k1 in component sodium_potassium_pump (per_millisecond).
C CONSTS(51) is k2 in component sodium_potassium_pump (per_millisecond).
C CONSTS(52) is k3 in component sodium_potassium_pump (per_millisecond).
C CONSTS(53) is k4 in component sodium_potassium_pump (per_millisecond).
C CONSTS(54) is Km_Nai in component sodium_potassium_pump (millimolar).
C CONSTS(55) is Km_Nao in component sodium_potassium_pump (millimolar).
C CONSTS(56) is Km_Ki in component sodium_potassium_pump (millimolar).
C CONSTS(57) is Km_Ko in component sodium_potassium_pump (millimolar).
C CONSTS(58) is Km_ATP in component sodium_potassium_pump (millimolar).
C ALGBRC(102) is Nao_Eff in component sodium_potassium_pump (millimolar).
C STATES(28) is y in component sodium_potassium_pump_y_gate (dimensionless).
C CONSTS(59) is P_NaK in component sodium_potassium_pump (picoA_per_picoF).
C ALGBRC(110) is alpha_y in component sodium_potassium_pump_y_gate (per_millisecond).
C ALGBRC(112) is beta_y in component sodium_potassium_pump_y_gate (per_millisecond).
C ALGBRC(111) is p_E2Ca in component SR_calcium_pump (dimensionless).
C ALGBRC(109) is p_E1Ca in component SR_calcium_pump (dimensionless).
C ALGBRC(113) is p_E1 in component SR_calcium_pump (dimensionless).
C ALGBRC(114) is p_E2 in component SR_calcium_pump (dimensionless).
C CONSTS(60) is k1 in component SR_calcium_pump (per_millisecond).
C ALGBRC(115) is k2 in component SR_calcium_pump (per_millisecond).
C CONSTS(61) is k3 in component SR_calcium_pump (per_millisecond).
C CONSTS(62) is k4 in component SR_calcium_pump (per_millisecond).
C CONSTS(63) is Km_CaSR in component SR_calcium_pump (millimolar).
C CONSTS(64) is Km_CaCyto in component SR_calcium_pump (millimolar).
C CONSTS(65) is Km_ATP in component SR_calcium_pump (millimolar).
C CONSTS(66) is i_max in component SR_calcium_pump (picoA).
C STATES(29) is Caup in component Ca_concentrations_in_SR (millimolar).
C STATES(30) is y in component SR_calcium_pump_y_gate (dimensionless).
C ALGBRC(117) is alpha_y in component SR_calcium_pump_y_gate (per_millisecond).
C ALGBRC(119) is beta_y in component SR_calcium_pump_y_gate (per_millisecond).
C CONSTS(67) is P_RyR in component RyR_channel (picoA_per_millimolar).
C ALGBRC(120) is k1 in component RyR_channel (per_millisecond).
C ALGBRC(125) is k2 in component RyR_channel (per_millisecond).
C ALGBRC(126) is k3 in component RyR_channel (per_millisecond).
C CONSTS(68) is k4 in component RyR_channel (per_millisecond).
C STATES(31) is p_open_RyR in component RyR_channel (dimensionless).
C STATES(32) is p_close_RyR in component RyR_channel (dimensionless).
C ALGBRC(124) is Carel in component Ca_concentrations_in_SR (millimolar).
C CONSTS(69) is Diadid_Factor in component RyR_channel (per_picoA_millisecond).
C ALGBRC(128) is i_SR_T in component SR_T_current (picoA).
C CONSTS(70) is P_SR_T in component SR_T_current (picoA_per_millimolar).
C CONSTS(71) is P_SR_L in component SR_L_current (picoA_per_millimolar).
C STATES(33) is Ca_Total in component Ca_concentrations_in_SR (millimolar).
C CONSTS(72) is V_rel in component Ca_concentrations_in_SR (micrometre3).
C CONSTS(73) is V_up in component Ca_concentrations_in_SR (micrometre3).
C CONSTS(74) is CSQN_max in component Ca_concentrations_in_SR (millimolar).
C CONSTS(75) is K_mCSQN in component Ca_concentrations_in_SR (millimolar).
C ALGBRC(122) is b1 in component Ca_concentrations_in_SR (millimolar).
C ALGBRC(123) is c1 in component Ca_concentrations_in_SR (millimolar2).
C CONSTS(99) is EffFraction in component NL_model (dimensionless).
C STATES(34) is pCa in component NL_model (dimensionless).
C STATES(35) is pCaCB in component NL_model (dimensionless).
C STATES(36) is pCB in component NL_model (dimensionless).
C ALGBRC(131) is p in component NL_model (dimensionless).
C CONSTS(76) is T_t in component NL_model (millimolar).
C ALGBRC(135) is Q_a in component NL_model (per_millisecond).
C ALGBRC(133) is Q_b in component NL_model (per_millisecond).
C ALGBRC(137) is Q_r in component NL_model (per_millisecond).
C ALGBRC(138) is Q_d in component NL_model (per_millisecond).
C ALGBRC(139) is Q_d1 in component NL_model (per_millisecond).
C ALGBRC(140) is Q_d2 in component NL_model (per_millisecond).
C CONSTS(77) is Y_1 in component NL_model (per_millimolar_millisecond).
C CONSTS(78) is Y_2 in component NL_model (per_millisecond).
C CONSTS(79) is Y_3 in component NL_model (per_millisecond).
C CONSTS(80) is Y_4 in component NL_model (per_millisecond).
C CONSTS(81) is Y_d in component NL_model (millisecond_per_micrometre2).
C CONSTS(82) is Z_1 in component NL_model (per_millisecond).
C CONSTS(83) is Z_2 in component NL_model (per_millisecond).
C CONSTS(84) is Z_3 in component NL_model (per_millimolar_millisecond).
C ALGBRC(129) is h in component NL_model (micrometre).
C CONSTS(85) is L_a in component NL_model (micrometre).
C CONSTS(86) is L in component NL_model (micrometre).
C ALGBRC(134) is ForceCB in component NL_model (mN_per_mm2).
C STATES(37) is X in component NL_model (micrometre).
C ALGBRC(132) is NewCBF in component NL_model (mN_per_mm2_micrometre).
C ALGBRC(130) is CBBound in component NL_model (millimolar).
C CONSTS(87) is KForceEC in component NL_model (mN_per_mm2_micrometre5).
C CONSTS(88) is ZeroForceEL in component NL_model (micrometre).
C CONSTS(89) is KForceLinearEc in component NL_model (mN_per_mm2_micrometre).
C CONSTS(90) is ForceFactor in component NL_model (mN_per_mm2_micrometre_millimolar).
C CONSTS(100) is ForceEcomp in component NL_model (mN_per_mm2).
C CONSTS(91) is B in component NL_model (per_millisecond).
C CONSTS(92) is h_c in component NL_model (micrometre).
C ALGBRC(136) 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 y in component time_independent_potassium_current_y_gate (dimensionless).
C RATES(20) is d/dt y1 in component rapid_time_dependent_potassium_current_y1_gate (dimensionless).
C RATES(21) is d/dt y2 in component rapid_time_dependent_potassium_current_y2_gate (dimensionless).
C RATES(22) is d/dt y3 in component rapid_time_dependent_potassium_current_y3_gate (dimensionless).
C RATES(23) is d/dt y1 in component slow_time_dependent_potassium_current_y1_gate (dimensionless).
C RATES(24) is d/dt y2 in component slow_time_dependent_potassium_current_y2_gate (dimensionless).
C RATES(25) is d/dt y1 in component transient_outward_current_y1_gate (dimensionless).
C RATES(26) is d/dt y2 in component transient_outward_current_y2_gate (dimensionless).
C RATES(27) is d/dt y in component sodium_calcium_exchanger_y_gate (dimensionless).
C RATES(28) is d/dt y in component sodium_potassium_pump_y_gate (dimensionless).
C RATES(30) is d/dt y in component SR_calcium_pump_y_gate (dimensionless).
C RATES(31) is d/dt p_open_RyR in component RyR_channel (dimensionless).
C RATES(32) is d/dt p_close_RyR in component RyR_channel (dimensionless).
C RATES(33) is d/dt Ca_Total in component Ca_concentrations_in_SR (millimolar).
C RATES(29) is d/dt Caup in component Ca_concentrations_in_SR (millimolar).
C RATES(37) is d/dt X in component NL_model (micrometre).
C RATES(34) is d/dt pCa in component NL_model (dimensionless).
C RATES(35) is d/dt pCaCB in component NL_model (dimensionless).
C RATES(36) is d/dt pCB in component NL_model (dimensionless).
C
      SUBROUTINE initConsts(CONSTS, RATES, STATES)
      REAL CONSTS(*), RATES(*), STATES(*)
      STATES(1) = -85.95752434460744
      CONSTS(1) = 8.3143
      CONSTS(2) = 310
      CONSTS(3) = 96.4867
      CONSTS(4) = 132
      CONSTS(5) = 50
      CONSTS(6) = 1000000
      CONSTS(7) = 400
      CONSTS(8) = 2
      CONSTS(9) = -4000
      CONSTS(10) = 140
      CONSTS(11) = 1.8
      CONSTS(12) = 5.4
      STATES(2) = 4.925761439682025
      STATES(3) = 143.1837333000449
      CONSTS(13) = 8000
      CONSTS(14) = 0.05
      CONSTS(15) = 0.00238
      STATES(4) = 4.0180173572968586e-4
      STATES(5) = 4.657102729020499
      CONSTS(16) = 0.003
      CONSTS(17) = 5
      CONSTS(18) = 2860
      STATES(6) = 1.779648367445368e-5
      STATES(7) = 0.5861887862983165
      STATES(8) = 0.3556412697995689
      STATES(9) = 0.40285968661346977
      CONSTS(19) = 0.0000875
      CONSTS(20) = 8712
      STATES(10) = 1.5445004166497696e-6
      STATES(11) = 0.17246483915629204
      STATES(12) = 6.098246017787626e-5
      STATES(13) = 0.9985266538252986
      STATES(14) = 0.9968480629364956
      STATES(15) = 8.77325391245903e-4
      CONSTS(21) = 0.004
      CONSTS(22) = 0.001
      STATES(16) = 0.4250747299372254
      CONSTS(23) = 0.0003
      CONSTS(24) = 0.35
      CONSTS(25) = 0.143
      CONSTS(26) = 0.35
      CONSTS(27) = 6.954
      CONSTS(28) = 0.0042
      CONSTS(29) = 6.954
      CONSTS(30) = 612
      STATES(17) = 1.6882718240109127e-5
      STATES(18) = 0.8585352091865849
      CONSTS(31) = 1.146
      STATES(19) = 0.6080573900752752
      CONSTS(32) = 0.00864
      STATES(20) = 0.0018339931180983765
      STATES(21) = 0.20443083454225305
      STATES(22) = 0.967887666264921
      STATES(23) = 0.09738789658609195
      STATES(24) = 0.09745345578743213
      CONSTS(33) = 5.04
      CONSTS(34) = 0.2016
      CONSTS(35) = 0.004444
      STATES(25) = 7.956883250874798e-4
      STATES(26) = 0.9999125083105881
      CONSTS(36) = 0.033
      CONSTS(37) = 0.00297
      CONSTS(38) = 0.385
      CONSTS(39) = 0.11
      CONSTS(40) = 0.0236
      CONSTS(41) = 2333
      CONSTS(42) = 0.04
      CONSTS(43) = 1
      CONSTS(44) = 1
      CONSTS(45) = 8.75
      CONSTS(46) = 87.5
      CONSTS(47) = 0.00138
      CONSTS(48) = 1.38
      STATES(27) = 0.9891789193465331
      CONSTS(49) = 6.81
      CONSTS(50) = 0.32
      CONSTS(51) = 0.04
      CONSTS(52) = 0.01
      CONSTS(53) = 0.165
      CONSTS(54) = 4.05
      CONSTS(55) = 69.8
      CONSTS(56) = 32.88
      CONSTS(57) = 0.258
      CONSTS(58) = 0.094
      STATES(28) = 0.5910747147428818
      CONSTS(59) = 21
      CONSTS(60) = 0.01
      CONSTS(61) = 1
      CONSTS(62) = 0.01
      CONSTS(63) = 0.08
      CONSTS(64) = 0.0008
      CONSTS(65) = 0.1
      CONSTS(66) = 162500
      STATES(29) = 2.611712901567567
      STATES(30) = 0.46108441538480216
      CONSTS(67) = 62000
      CONSTS(68) = 0.000849
      STATES(31) = 3.4314360001543243e-4
      STATES(32) = 0.19135178123107768
      CONSTS(69) = -150
      CONSTS(70) = 386
      CONSTS(71) = 459
      STATES(33) = 9.455741736977666
      CONSTS(72) = 160
      CONSTS(73) = 400
      CONSTS(74) = 10
      CONSTS(75) = 0.8
      STATES(34) = 0.02490898775497523
      STATES(35) = 0.001990153835322864
      STATES(36) = 4.2941813853474524e-4
      CONSTS(76) = 0.07
      CONSTS(77) = 39
      CONSTS(78) = 0.0039
      CONSTS(79) = 0.03
      CONSTS(80) = 0.12
      CONSTS(81) = 0.027
      CONSTS(82) = 0.03
      CONSTS(83) = 0.0039
      CONSTS(84) = 1560
      CONSTS(85) = 1.17
      CONSTS(86) = 0.9623799975411884
      STATES(37) = 0.9573749975411884
      CONSTS(87) = 140000
      CONSTS(88) = 0.97
      CONSTS(89) = 200
      CONSTS(90) = 1800000
      CONSTS(91) = 1.2
      CONSTS(92) = 0.005
      CONSTS(93) = ( CONSTS(28)*CONSTS(25)*CONSTS(27)*CONSTS(24))/( CONSTS(26)*CONSTS(29)*CONSTS(23))
      CONSTS(94) =  CONSTS(31)*CONSTS(4)*CONSTS(12)/5.40000 ** 0.400000
      CONSTS(95) =  CONSTS(32)*CONSTS(4)*CONSTS(12)/5.40000 ** 0.200000
      CONSTS(96) =  0.000110000*CONSTS(12)/5.40000 ** 0.160000
      CONSTS(97) =  CONSTS(40)*CONSTS(41)*CONSTS(12)/1.00000 ** 0.240000
      CONSTS(98) = 1.00000/(1.00000+ CONSTS(46)/CONSTS(10) ** 3.00000*(1.00000+CONSTS(11)/CONSTS(48)))
      CONSTS(99) = EXP( - 20.0000*CONSTS(86) - CONSTS(85) ** 2.00000)
      CONSTS(100) =  CONSTS(87)*CONSTS(88) - CONSTS(86) ** 5.00000+ CONSTS(89)*(CONSTS(88) - CONSTS(86))
      CONSTS(101) = 1.00000/(1.00000+ (CONSTS(48)/CONSTS(11))*(1.00000+CONSTS(10)/CONSTS(46) ** 3.00000))
      RETURN
      END
      SUBROUTINE computeRates(VOI, CONSTS,  RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/( 9.00000e+09*EXP(STATES(1)/5.00000)+ 8000.00*EXP(STATES(1)/100.000))
      ALGBRC(16) = 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(16)*STATES(7)
      ALGBRC(5) = 1.00000/( 250000.*EXP(STATES(1)/9.00000)+ 58.0000*EXP(STATES(1)/65.0000))
      ALGBRC(18) = 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(18)*STATES(13)
      ALGBRC(6) = 1.00000/( 0.0190000*EXP(- STATES(1)/5.60000)+ 0.820000*EXP(- STATES(1)/250.000))
      ALGBRC(19) = 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(19)*STATES(17)
      ALGBRC(7) = 1.00000/( 62000.0*EXP(STATES(1)/10.1000)+ 30.0000*EXP(STATES(1)/3000.00))
      ALGBRC(20) = 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(20)*STATES(18)
      ALGBRC(8) = 1.00000/( 20.0000*EXP(- STATES(1)/11.5000)+ 5.00000*EXP(- STATES(1)/300.000))
      ALGBRC(21) = 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(20) =  ALGBRC(8)*(1.00000 - STATES(20)) -  ALGBRC(21)*STATES(20)
      ALGBRC(9) = 1.00000/( 200.000*EXP(- STATES(1)/13.0000)+ 20.0000*EXP(- STATES(1)/300.000))
      ALGBRC(22) = 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(21) =  ALGBRC(9)*(1.00000 - STATES(21)) -  ALGBRC(22)*STATES(21)
      ALGBRC(10) = 1.00000/( 10.0000*EXP(STATES(1)/17.0000)+ 2.50000*EXP(STATES(1)/300.000))
      ALGBRC(23) = 1.00000/( 0.350000*EXP(- STATES(1)/17.0000)+ 2.00000*EXP(- STATES(1)/150.000))
      RATES(22) =  ALGBRC(10)*(1.00000 - STATES(22)) -  ALGBRC(23)*STATES(22)
      ALGBRC(11) = 1.00000/( 85.0000*EXP(- STATES(1)/10.5000)+ 370.000*EXP(- STATES(1)/62.0000))
      ALGBRC(24) = 1.00000/( 1450.00*EXP(STATES(1)/20.0000)+ 260.000*EXP(STATES(1)/100.000))
      RATES(23) =  ALGBRC(11)*(1.00000 - STATES(23)) -  ALGBRC(24)*STATES(23)
      ALGBRC(12) = 1.00000/( 11.0000*EXP(- STATES(1)/28.0000)+ 0.200000*EXP(- STATES(1)/400.000))
      ALGBRC(25) = 1.00000/( 4.40000*EXP(STATES(1)/16.0000)+ 0.200000*EXP(STATES(1)/500.000))
      RATES(25) =  ALGBRC(12)*(1.00000 - STATES(25)) -  ALGBRC(25)*STATES(25)
      ALGBRC(13) = ( 0.00380000*EXP(- (STATES(1)+13.5000)/11.3000))/(1.00000+ 0.0513350*EXP(- (STATES(1)+13.5000)/11.3000))
      ALGBRC(26) = ( 0.00380000*EXP((STATES(1)+13.5000)/11.3000))/(1.00000+ 0.0670830*EXP((STATES(1)+13.5000)/11.3000))
      RATES(26) =  ALGBRC(13)*(1.00000 - STATES(26)) -  ALGBRC(26)*STATES(26)
      ALGBRC(17) = 1.00000/( 0.270000*EXP(- STATES(1)/5.90000)+ 1.50000*EXP(- STATES(1)/65.0000))
      ALGBRC(29) = 1.00000/( 480.000*EXP(STATES(1)/7.00000)+ 2.20000*EXP(STATES(1)/65.0000))
      RATES(10) = ( STATES(14)*ALGBRC(17)+ STATES(15)*CONSTS(22)) -  STATES(10)*(ALGBRC(29)+CONSTS(21))
      ALGBRC(15) = 1.00000/( 0.102700*EXP(- STATES(1)/8.00000)+ 0.250000*EXP(- STATES(1)/50.0000))
      ALGBRC(28) = 1.00000/( 26.0000*EXP(STATES(1)/17.0000)+ 0.0200000*EXP(STATES(1)/800.000))
      ALGBRC(31) = 1.00000/( 0.800000*EXP(- STATES(1)/400.000))
      RATES(6) = ( STATES(8)*ALGBRC(15)+ STATES(9)*CONSTS(19)) -  STATES(6)*(ALGBRC(28)+ALGBRC(31))
      ALGBRC(4) = ((1.00000 - STATES(10)) - STATES(14)) - STATES(15)
      ALGBRC(32) = 1.00000/( 0.00180000*EXP(- STATES(1)/7.40000)+ 2.00000*EXP(- STATES(1)/100.000))
      ALGBRC(35) = 1.00000/( 2.20000e+06*EXP(STATES(1)/7.40000)+ 11.0000*EXP(STATES(1)/100.000))
      RATES(15) = ( ALGBRC(4)*ALGBRC(32)+ STATES(10)*CONSTS(21)) -  STATES(15)*(ALGBRC(35)+CONSTS(22))
      ALGBRC(14) = (CONSTS(14) - STATES(4))+CONSTS(15)
      ALGBRC(27) =  CONSTS(15)*STATES(4)
      ALGBRC(30) = ( (ALGBRC(14) ** 2.00000+ 4.00000*ALGBRC(27)) ** (1.0 / 2) - ALGBRC(14))/2.00000
      ALGBRC(36) =  3.70000*ALGBRC(30)
      RATES(24) =  ALGBRC(36)*(1.00000 - STATES(24)) -  CONSTS(35)*STATES(24)
      ALGBRC(2) = ((1.00000 - STATES(8)) - STATES(6)) - STATES(9)
      ALGBRC(38) = 1.00000/( 0.000102700*EXP(- STATES(1)/8.00000)+ 5.00000*EXP(- STATES(1)/400.000))
      ALGBRC(34) = 1.00000/( 1300.00*EXP(STATES(1)/20.0000)+ 0.0400000*EXP(STATES(1)/800.000))
      RATES(9) = ( ALGBRC(2)*ALGBRC(38)+ STATES(6)*ALGBRC(31)) -  STATES(9)*(ALGBRC(34)+CONSTS(19))
      ALGBRC(39) = 0.0400000/(1.00000+( CONSTS(22)*ALGBRC(29)*ALGBRC(32))/( CONSTS(21)*ALGBRC(17)*ALGBRC(35)))
      ALGBRC(42) = 0.0400000 - ALGBRC(39)
      RATES(14) = ( STATES(10)*ALGBRC(29)+ ALGBRC(4)*ALGBRC(42)) -  STATES(14)*(ALGBRC(39)+ALGBRC(17))
      ALGBRC(41) = 0.0100000/(1.00000+( CONSTS(19)*ALGBRC(28)*ALGBRC(38))/( ALGBRC(31)*ALGBRC(15)*ALGBRC(34)))
      ALGBRC(44) = 0.0100000 - ALGBRC(41)
      RATES(8) = ( STATES(6)*ALGBRC(28)+ ALGBRC(2)*ALGBRC(44)) -  STATES(8)*(ALGBRC(41)+ALGBRC(15))
      ALGBRC(37) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(11), ( (( 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(ALGBRC(30) -  CONSTS(11)*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(52) =  0.0676000*ALGBRC(37)
      ALGBRC(54) = ALGBRC(30) -  0.300000*ALGBRC(52)
      ALGBRC(56) =  ALGBRC(54)*STATES(10)
      ALGBRC(58) = ALGBRC(56)+ ALGBRC(30)*(1.00000 - STATES(10))
      ALGBRC(61) =  CONSTS(27)*ALGBRC(58)
      RATES(11) = ( STATES(16)*CONSTS(25)+ STATES(12)*CONSTS(93)) -  STATES(11)*(ALGBRC(61)+CONSTS(26))
      ALGBRC(63) = ((1.00000 - STATES(16)) - STATES(11)) - STATES(12)
      RATES(12) = ( STATES(11)*ALGBRC(61)+ ALGBRC(63)*CONSTS(23)) -  STATES(12)*(CONSTS(24)+CONSTS(93))
      RATES(16) = ( ALGBRC(63)*CONSTS(28)+ STATES(11)*CONSTS(26)) -  STATES(16)*(CONSTS(25)+ CONSTS(29)*ALGBRC(54)*STATES(10))
      ALGBRC(57) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(12)/STATES(3))
      ALGBRC(67) = 1.00000/( 8000.00*EXP(((STATES(1) - ALGBRC(57)) - 97.0000)/8.50000)+ 7.00000*EXP(((STATES(1) - ALGBRC(57)) - 97.0000)/300.000))
      ALGBRC(60) = ( 0.750000*EXP( 0.0350000*((STATES(1) - ALGBRC(57)) - 10.0000)))/(1.00000+EXP( 0.0150000*((STATES(1) - ALGBRC(57)) - 140.000)))
      ALGBRC(62) = ( 3.00000*EXP( - 0.0480000*((STATES(1) - ALGBRC(57)) - 10.0000))*(1.00000+EXP( 0.0640000*((STATES(1) - ALGBRC(57)) - 38.0000))))/(1.00000+EXP( 0.0300000*((STATES(1) - ALGBRC(57)) - 70.0000)))
      ALGBRC(65) = ALGBRC(62)/(ALGBRC(60)+ALGBRC(62))
      ALGBRC(69) = ( ALGBRC(65) ** 4.00000*1.00000)/( 0.000140000*EXP(- ((STATES(1) - ALGBRC(57)) - 97.0000)/9.10000)+ 0.200000*EXP(- ((STATES(1) - ALGBRC(57)) - 97.0000)/500.000))
      RATES(19) =  ALGBRC(67)*(1.00000 - STATES(19)) -  ALGBRC(69)*STATES(19)
      ALGBRC(94) =  1.00000*EXP(( (CONSTS(50) - 1.00000)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(96) =  ALGBRC(94)*CONSTS(98)+ CONSTS(44)*CONSTS(101)
      ALGBRC(91) = 1.00000/(1.00000+ CONSTS(45)/STATES(2) ** 3.00000*(1.00000+ALGBRC(30)/CONSTS(47)))
      ALGBRC(92) = 1.00000/(1.00000+ (CONSTS(47)/ALGBRC(30))*(1.00000+STATES(2)/CONSTS(45) ** 3.00000))
      ALGBRC(93) =  1.00000*EXP(( CONSTS(50)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(98) =  ALGBRC(93)*ALGBRC(91)+ CONSTS(43)*ALGBRC(92)
      RATES(27) =  ALGBRC(96)*(1.00000 - STATES(27)) -  ALGBRC(98)*STATES(27)
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) -  INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(9), 0.00000)
      ALGBRC(33) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(10), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(2) -  CONSTS(10)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(80) =  CONSTS(38)*ALGBRC(33)
      ALGBRC(40) = TERNRY(STATES(1).EQ.0.00000, STATES(3), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(3) -  CONSTS(12)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(79) =  0.400000*CONSTS(38)*ALGBRC(40)
      ALGBRC(81) = ALGBRC(79)+ALGBRC(80)
      ALGBRC(89) =  CONSTS(42)*ALGBRC(37)
      ALGBRC(82) = TERNRY(STATES(1).EQ.- 3.00000,  CONSTS(96)*ALGBRC(40)*13.0077, ( CONSTS(96)*ALGBRC(40)*(STATES(1)+3.00000))/(1.00000 - EXP(- (STATES(1)+3.00000)/13.0000)))
      ALGBRC(83) = 1.00000/(1.00000+0.00120000/ALGBRC(30) ** 3.00000)
      ALGBRC(85) =  CONSTS(39)*ALGBRC(33)*ALGBRC(83)
      ALGBRC(84) =  CONSTS(39)*ALGBRC(40)*ALGBRC(83)
      ALGBRC(86) = ALGBRC(84)+ALGBRC(85)
      ALGBRC(87) = 0.800000/(1.00000+STATES(5)/0.100000 ** 2.00000)
      ALGBRC(88) =  CONSTS(97)*(STATES(1) - ALGBRC(57))*ALGBRC(87)
      ALGBRC(90) = ALGBRC(81)+ALGBRC(89)+ALGBRC(82)+ALGBRC(86)+ALGBRC(88)
      ALGBRC(43) =  CONSTS(18)*ALGBRC(33)*STATES(6)*STATES(7)
      ALGBRC(45) =  0.100000*CONSTS(18)*ALGBRC(40)*STATES(6)*STATES(7)
      ALGBRC(46) = ALGBRC(43)+ALGBRC(45)
      ALGBRC(47) = ( STATES(10)*(STATES(11)+STATES(12))*STATES(13))/(1.00000+1.40000/STATES(5) ** 3.00000)
      ALGBRC(49) =  1.85000e-05*CONSTS(20)*ALGBRC(33)*ALGBRC(47)
      ALGBRC(50) =  0.000365000*CONSTS(20)*ALGBRC(40)*ALGBRC(47)
      ALGBRC(48) =  CONSTS(20)*ALGBRC(37)*ALGBRC(47)
      ALGBRC(51) = ALGBRC(49)+ALGBRC(48)+ALGBRC(50)
      ALGBRC(55) =  CONSTS(30)*ALGBRC(37)*STATES(17)*STATES(18)
      ALGBRC(64) = ALGBRC(60)/(ALGBRC(60)+ALGBRC(62))
      ALGBRC(66) =  2.00000*ALGBRC(65) ** 2.00000*ALGBRC(64) ** 2.00000
      ALGBRC(68) =  (8.00000/3.00000)*ALGBRC(65) ** 3.00000*ALGBRC(64)
      ALGBRC(70) = ALGBRC(65) ** 4.00000
      ALGBRC(71) =  CONSTS(94)*(STATES(1) - ALGBRC(57))*(ALGBRC(70)+ALGBRC(68)+ALGBRC(66))*STATES(19)
      ALGBRC(72) =  CONSTS(95)*(STATES(1) - ALGBRC(57))*( 0.600000*STATES(20)+ 0.400000*STATES(21))*STATES(22)
      ALGBRC(73) =  CONSTS(33)*ALGBRC(40)*STATES(23) ** 2.00000*( 0.900000*STATES(24)+0.100000)
      ALGBRC(74) =  CONSTS(34)*ALGBRC(33)*STATES(23) ** 2.00000*( 0.900000*STATES(24)+0.100000)
      ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
      ALGBRC(76) =  CONSTS(36)*ALGBRC(40)*STATES(25) ** 3.00000*STATES(26)
      ALGBRC(77) =  CONSTS(37)*ALGBRC(33)*STATES(25) ** 3.00000*STATES(26)
      ALGBRC(78) = ALGBRC(77)+ALGBRC(76)
      ALGBRC(102) =  CONSTS(10)*EXP(( - 0.820000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(103) = 1.00000/(1.00000+ CONSTS(55)/ALGBRC(102) ** 1.06000*(1.00000+CONSTS(12)/CONSTS(57) ** 1.12000))
      ALGBRC(99) = 1.00000/(1.00000+ CONSTS(54)/STATES(2) ** 1.06000*(1.00000+STATES(3)/CONSTS(56) ** 1.12000))
      ALGBRC(101) = 0.370000/(1.00000+CONSTS(58)/STATES(5))
      ALGBRC(104) =  CONSTS(59)*CONSTS(4)*1.00000*( ALGBRC(101)*ALGBRC(99)*STATES(28) -  CONSTS(51)*ALGBRC(103)*(1.00000 - STATES(28)))
      ALGBRC(95) =  CONSTS(49)*CONSTS(4)*1.00000*( ALGBRC(93)*ALGBRC(91)*STATES(27) -  ALGBRC(94)*CONSTS(98)*(1.00000 - STATES(27)))
      ALGBRC(106) = ALGBRC(46)+ALGBRC(51)+ALGBRC(55)+ALGBRC(71)+ALGBRC(72)+ALGBRC(75)+ALGBRC(78)+ALGBRC(90)+ALGBRC(104)+ALGBRC(95)
      RATES(1) = - (ALGBRC(106)+ALGBRC(1))/CONSTS(4)
      ALGBRC(107) = ALGBRC(43)+ALGBRC(74)+ALGBRC(77)+ALGBRC(49)+ALGBRC(80)+ALGBRC(85)+ 3.00000*ALGBRC(104)+ 3.00000*ALGBRC(95)
      RATES(2) = - ALGBRC(107)/( CONSTS(3)*CONSTS(13))
      ALGBRC(108) = (ALGBRC(71)+ALGBRC(72)+ALGBRC(76)+ALGBRC(88)+ALGBRC(73)+ALGBRC(45)+ALGBRC(50)+ALGBRC(79)+ALGBRC(84)+ALGBRC(82)) -  2.00000*ALGBRC(104)
      RATES(3) = - (ALGBRC(108)+ALGBRC(1))/( CONSTS(3)*CONSTS(13))
      ALGBRC(105) = 1.00000/(1.00000+ CONSTS(57)/CONSTS(12) ** 1.12000*(1.00000+ALGBRC(102)/CONSTS(55) ** 1.06000))
      ALGBRC(110) =  CONSTS(51)*ALGBRC(103)+ CONSTS(53)*ALGBRC(105)
      ALGBRC(100) = 1.00000/(1.00000+ CONSTS(56)/STATES(3) ** 1.12000*(1.00000+STATES(2)/CONSTS(54) ** 1.06000))
      ALGBRC(112) =  ALGBRC(101)*ALGBRC(99)+ CONSTS(52)*ALGBRC(100)
      RATES(28) =  ALGBRC(110)*(1.00000 - STATES(28)) -  ALGBRC(112)*STATES(28)
      ALGBRC(111) = 1.00000/(1.00000+CONSTS(64)/ALGBRC(30))
      ALGBRC(109) = 1.00000/(1.00000+CONSTS(63)/STATES(29))
      ALGBRC(115) = 1.00000/(1.00000+CONSTS(65)/STATES(5))
      ALGBRC(116) =  CONSTS(66)*1.00000*( CONSTS(60)*ALGBRC(109)*STATES(30) -  ALGBRC(115)*ALGBRC(111)*(1.00000 - STATES(30)))
      ALGBRC(118) =  - 0.400000*STATES(35)*CONSTS(76)
      RATES(5) = (( CONSTS(16)*(CONSTS(17) - STATES(5))+ALGBRC(118)) - ALGBRC(104)/( CONSTS(3)*CONSTS(13)))+ALGBRC(116)/( 4.00000*CONSTS(3)*CONSTS(13))
      ALGBRC(114) = 1.00000 - ALGBRC(111)
      ALGBRC(117) =  ALGBRC(115)*ALGBRC(111)+ CONSTS(62)*ALGBRC(114)
      ALGBRC(113) = 1.00000 - ALGBRC(109)
      ALGBRC(119) =  CONSTS(60)*ALGBRC(109)+ CONSTS(61)*ALGBRC(113)
      RATES(30) =  ALGBRC(117)*(1.00000 - STATES(30)) -  ALGBRC(119)*STATES(30)
      ALGBRC(53) =  ALGBRC(52)*ALGBRC(47)
      ALGBRC(120) =  280000.*ALGBRC(30)/1.00000 ** 2.00000+ CONSTS(69)*ALGBRC(53)
      ALGBRC(122) = (CONSTS(74) - STATES(33))+CONSTS(75)
      ALGBRC(123) =  CONSTS(75)*STATES(33)
      ALGBRC(124) = ( (ALGBRC(122) ** 2.00000+ 4.00000*ALGBRC(123)) ** (1.0 / 2) - ALGBRC(122))/2.00000
      ALGBRC(125) = 0.0800000/(1.00000+0.360000/ALGBRC(124))
      RATES(31) =  STATES(32)*ALGBRC(120) -  STATES(31)*ALGBRC(125)
      ALGBRC(126) =  0.000377000*ALGBRC(124)/1.00000 ** 2.00000
      RATES(32) =  ALGBRC(126)*(1.00000 - (STATES(31)+STATES(32))) -  (ALGBRC(120)+CONSTS(68))*STATES(32)
      ALGBRC(127) =  CONSTS(67)*(ALGBRC(124) - ALGBRC(30))*STATES(31)
      ALGBRC(128) =  CONSTS(70)*(STATES(29) - ALGBRC(124))
      RATES(33) = (ALGBRC(128) - ALGBRC(127))/( 2.00000*CONSTS(3)*CONSTS(72))
      ALGBRC(121) =  CONSTS(71)*(STATES(29) - ALGBRC(30))
      RATES(29) = ((- ALGBRC(116) - ALGBRC(128)) - ALGBRC(121))/( 2.00000*CONSTS(3)*CONSTS(73))
      ALGBRC(135) =  CONSTS(78)*STATES(34)*CONSTS(99) -  CONSTS(83)*STATES(35)
      ALGBRC(131) = ((1.00000 - STATES(34)) - STATES(35)) - STATES(36)
      ALGBRC(133) =  CONSTS(77)*ALGBRC(30)*ALGBRC(131) -  CONSTS(82)*STATES(34)
      RATES(34) = ALGBRC(133) - ALGBRC(135)
      ALGBRC(129) = CONSTS(86) - STATES(37)
      RATES(37) =  CONSTS(91)*(ALGBRC(129) - CONSTS(92))
      ALGBRC(137) =  CONSTS(79)*STATES(35) -  CONSTS(84)*STATES(36)*ALGBRC(30)
      ALGBRC(140) =  CONSTS(81)*RATES(37) ** 2.00000*STATES(35)
      RATES(35) = (ALGBRC(135) - ALGBRC(137)) - ALGBRC(140)
      ALGBRC(138) =  CONSTS(80)*STATES(36)
      ALGBRC(139) =  CONSTS(81)*RATES(37) ** 2.00000*STATES(36)
      RATES(36) = (ALGBRC(137) - ALGBRC(138)) - ALGBRC(139)
      ALGBRC(97) = (ALGBRC(48)+ALGBRC(55)+ALGBRC(89)) -  2.00000*ALGBRC(95)
      ALGBRC(141) =  CONSTS(76)*((ALGBRC(140)+ALGBRC(137)) - ALGBRC(133))
      RATES(4) = - (((ALGBRC(97) - ALGBRC(116)) - ALGBRC(127)) - ALGBRC(121))/( 2.00000*CONSTS(3)*CONSTS(13))+ALGBRC(141)
      RETURN
      END
      SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC)
      REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*)
      ALGBRC(3) = 1.00000/( 9.00000e+09*EXP(STATES(1)/5.00000)+ 8000.00*EXP(STATES(1)/100.000))
      ALGBRC(16) = 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(18) = 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(19) = 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(20) = 1.00000/( 0.000600000*EXP(- STATES(1)/6.70000)+ 1.20000*EXP(- STATES(1)/25.0000))
      ALGBRC(8) = 1.00000/( 20.0000*EXP(- STATES(1)/11.5000)+ 5.00000*EXP(- STATES(1)/300.000))
      ALGBRC(21) = 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(9) = 1.00000/( 200.000*EXP(- STATES(1)/13.0000)+ 20.0000*EXP(- STATES(1)/300.000))
      ALGBRC(22) = 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(10) = 1.00000/( 10.0000*EXP(STATES(1)/17.0000)+ 2.50000*EXP(STATES(1)/300.000))
      ALGBRC(23) = 1.00000/( 0.350000*EXP(- STATES(1)/17.0000)+ 2.00000*EXP(- STATES(1)/150.000))
      ALGBRC(11) = 1.00000/( 85.0000*EXP(- STATES(1)/10.5000)+ 370.000*EXP(- STATES(1)/62.0000))
      ALGBRC(24) = 1.00000/( 1450.00*EXP(STATES(1)/20.0000)+ 260.000*EXP(STATES(1)/100.000))
      ALGBRC(12) = 1.00000/( 11.0000*EXP(- STATES(1)/28.0000)+ 0.200000*EXP(- STATES(1)/400.000))
      ALGBRC(25) = 1.00000/( 4.40000*EXP(STATES(1)/16.0000)+ 0.200000*EXP(STATES(1)/500.000))
      ALGBRC(13) = ( 0.00380000*EXP(- (STATES(1)+13.5000)/11.3000))/(1.00000+ 0.0513350*EXP(- (STATES(1)+13.5000)/11.3000))
      ALGBRC(26) = ( 0.00380000*EXP((STATES(1)+13.5000)/11.3000))/(1.00000+ 0.0670830*EXP((STATES(1)+13.5000)/11.3000))
      ALGBRC(17) = 1.00000/( 0.270000*EXP(- STATES(1)/5.90000)+ 1.50000*EXP(- STATES(1)/65.0000))
      ALGBRC(29) = 1.00000/( 480.000*EXP(STATES(1)/7.00000)+ 2.20000*EXP(STATES(1)/65.0000))
      ALGBRC(15) = 1.00000/( 0.102700*EXP(- STATES(1)/8.00000)+ 0.250000*EXP(- STATES(1)/50.0000))
      ALGBRC(28) = 1.00000/( 26.0000*EXP(STATES(1)/17.0000)+ 0.0200000*EXP(STATES(1)/800.000))
      ALGBRC(31) = 1.00000/( 0.800000*EXP(- STATES(1)/400.000))
      ALGBRC(4) = ((1.00000 - STATES(10)) - STATES(14)) - STATES(15)
      ALGBRC(32) = 1.00000/( 0.00180000*EXP(- STATES(1)/7.40000)+ 2.00000*EXP(- STATES(1)/100.000))
      ALGBRC(35) = 1.00000/( 2.20000e+06*EXP(STATES(1)/7.40000)+ 11.0000*EXP(STATES(1)/100.000))
      ALGBRC(14) = (CONSTS(14) - STATES(4))+CONSTS(15)
      ALGBRC(27) =  CONSTS(15)*STATES(4)
      ALGBRC(30) = ( (ALGBRC(14) ** 2.00000+ 4.00000*ALGBRC(27)) ** (1.0 / 2) - ALGBRC(14))/2.00000
      ALGBRC(36) =  3.70000*ALGBRC(30)
      ALGBRC(2) = ((1.00000 - STATES(8)) - STATES(6)) - STATES(9)
      ALGBRC(38) = 1.00000/( 0.000102700*EXP(- STATES(1)/8.00000)+ 5.00000*EXP(- STATES(1)/400.000))
      ALGBRC(34) = 1.00000/( 1300.00*EXP(STATES(1)/20.0000)+ 0.0400000*EXP(STATES(1)/800.000))
      ALGBRC(39) = 0.0400000/(1.00000+( CONSTS(22)*ALGBRC(29)*ALGBRC(32))/( CONSTS(21)*ALGBRC(17)*ALGBRC(35)))
      ALGBRC(42) = 0.0400000 - ALGBRC(39)
      ALGBRC(41) = 0.0100000/(1.00000+( CONSTS(19)*ALGBRC(28)*ALGBRC(38))/( ALGBRC(31)*ALGBRC(15)*ALGBRC(34)))
      ALGBRC(44) = 0.0100000 - ALGBRC(41)
      ALGBRC(37) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(11), ( (( 2.00000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(ALGBRC(30) -  CONSTS(11)*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(52) =  0.0676000*ALGBRC(37)
      ALGBRC(54) = ALGBRC(30) -  0.300000*ALGBRC(52)
      ALGBRC(56) =  ALGBRC(54)*STATES(10)
      ALGBRC(58) = ALGBRC(56)+ ALGBRC(30)*(1.00000 - STATES(10))
      ALGBRC(61) =  CONSTS(27)*ALGBRC(58)
      ALGBRC(63) = ((1.00000 - STATES(16)) - STATES(11)) - STATES(12)
      ALGBRC(57) =  (( CONSTS(1)*CONSTS(2))/CONSTS(3))*log(CONSTS(12)/STATES(3))
      ALGBRC(67) = 1.00000/( 8000.00*EXP(((STATES(1) - ALGBRC(57)) - 97.0000)/8.50000)+ 7.00000*EXP(((STATES(1) - ALGBRC(57)) - 97.0000)/300.000))
      ALGBRC(60) = ( 0.750000*EXP( 0.0350000*((STATES(1) - ALGBRC(57)) - 10.0000)))/(1.00000+EXP( 0.0150000*((STATES(1) - ALGBRC(57)) - 140.000)))
      ALGBRC(62) = ( 3.00000*EXP( - 0.0480000*((STATES(1) - ALGBRC(57)) - 10.0000))*(1.00000+EXP( 0.0640000*((STATES(1) - ALGBRC(57)) - 38.0000))))/(1.00000+EXP( 0.0300000*((STATES(1) - ALGBRC(57)) - 70.0000)))
      ALGBRC(65) = ALGBRC(62)/(ALGBRC(60)+ALGBRC(62))
      ALGBRC(69) = ( ALGBRC(65) ** 4.00000*1.00000)/( 0.000140000*EXP(- ((STATES(1) - ALGBRC(57)) - 97.0000)/9.10000)+ 0.200000*EXP(- ((STATES(1) - ALGBRC(57)) - 97.0000)/500.000))
      ALGBRC(94) =  1.00000*EXP(( (CONSTS(50) - 1.00000)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(96) =  ALGBRC(94)*CONSTS(98)+ CONSTS(44)*CONSTS(101)
      ALGBRC(91) = 1.00000/(1.00000+ CONSTS(45)/STATES(2) ** 3.00000*(1.00000+ALGBRC(30)/CONSTS(47)))
      ALGBRC(92) = 1.00000/(1.00000+ (CONSTS(47)/ALGBRC(30))*(1.00000+STATES(2)/CONSTS(45) ** 3.00000))
      ALGBRC(93) =  1.00000*EXP(( CONSTS(50)*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(98) =  ALGBRC(93)*ALGBRC(91)+ CONSTS(43)*ALGBRC(92)
      ALGBRC(1) = TERNRY(VOI.GE.CONSTS(5).AND.VOI.LE.CONSTS(6).AND.(VOI - CONSTS(5)) -  INT((VOI - CONSTS(5))/CONSTS(7))*CONSTS(7).LE.CONSTS(8), CONSTS(9), 0.00000)
      ALGBRC(33) = TERNRY(STATES(1).EQ.0.00000, - CONSTS(10), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(2) -  CONSTS(10)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(80) =  CONSTS(38)*ALGBRC(33)
      ALGBRC(40) = TERNRY(STATES(1).EQ.0.00000, STATES(3), ( (( CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))*(STATES(3) -  CONSTS(12)*EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))/(1.00000 - EXP(( - CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))))
      ALGBRC(79) =  0.400000*CONSTS(38)*ALGBRC(40)
      ALGBRC(81) = ALGBRC(79)+ALGBRC(80)
      ALGBRC(89) =  CONSTS(42)*ALGBRC(37)
      ALGBRC(82) = TERNRY(STATES(1).EQ.- 3.00000,  CONSTS(96)*ALGBRC(40)*13.0077, ( CONSTS(96)*ALGBRC(40)*(STATES(1)+3.00000))/(1.00000 - EXP(- (STATES(1)+3.00000)/13.0000)))
      ALGBRC(83) = 1.00000/(1.00000+0.00120000/ALGBRC(30) ** 3.00000)
      ALGBRC(85) =  CONSTS(39)*ALGBRC(33)*ALGBRC(83)
      ALGBRC(84) =  CONSTS(39)*ALGBRC(40)*ALGBRC(83)
      ALGBRC(86) = ALGBRC(84)+ALGBRC(85)
      ALGBRC(87) = 0.800000/(1.00000+STATES(5)/0.100000 ** 2.00000)
      ALGBRC(88) =  CONSTS(97)*(STATES(1) - ALGBRC(57))*ALGBRC(87)
      ALGBRC(90) = ALGBRC(81)+ALGBRC(89)+ALGBRC(82)+ALGBRC(86)+ALGBRC(88)
      ALGBRC(43) =  CONSTS(18)*ALGBRC(33)*STATES(6)*STATES(7)
      ALGBRC(45) =  0.100000*CONSTS(18)*ALGBRC(40)*STATES(6)*STATES(7)
      ALGBRC(46) = ALGBRC(43)+ALGBRC(45)
      ALGBRC(47) = ( STATES(10)*(STATES(11)+STATES(12))*STATES(13))/(1.00000+1.40000/STATES(5) ** 3.00000)
      ALGBRC(49) =  1.85000e-05*CONSTS(20)*ALGBRC(33)*ALGBRC(47)
      ALGBRC(50) =  0.000365000*CONSTS(20)*ALGBRC(40)*ALGBRC(47)
      ALGBRC(48) =  CONSTS(20)*ALGBRC(37)*ALGBRC(47)
      ALGBRC(51) = ALGBRC(49)+ALGBRC(48)+ALGBRC(50)
      ALGBRC(55) =  CONSTS(30)*ALGBRC(37)*STATES(17)*STATES(18)
      ALGBRC(64) = ALGBRC(60)/(ALGBRC(60)+ALGBRC(62))
      ALGBRC(66) =  2.00000*ALGBRC(65) ** 2.00000*ALGBRC(64) ** 2.00000
      ALGBRC(68) =  (8.00000/3.00000)*ALGBRC(65) ** 3.00000*ALGBRC(64)
      ALGBRC(70) = ALGBRC(65) ** 4.00000
      ALGBRC(71) =  CONSTS(94)*(STATES(1) - ALGBRC(57))*(ALGBRC(70)+ALGBRC(68)+ALGBRC(66))*STATES(19)
      ALGBRC(72) =  CONSTS(95)*(STATES(1) - ALGBRC(57))*( 0.600000*STATES(20)+ 0.400000*STATES(21))*STATES(22)
      ALGBRC(73) =  CONSTS(33)*ALGBRC(40)*STATES(23) ** 2.00000*( 0.900000*STATES(24)+0.100000)
      ALGBRC(74) =  CONSTS(34)*ALGBRC(33)*STATES(23) ** 2.00000*( 0.900000*STATES(24)+0.100000)
      ALGBRC(75) = ALGBRC(74)+ALGBRC(73)
      ALGBRC(76) =  CONSTS(36)*ALGBRC(40)*STATES(25) ** 3.00000*STATES(26)
      ALGBRC(77) =  CONSTS(37)*ALGBRC(33)*STATES(25) ** 3.00000*STATES(26)
      ALGBRC(78) = ALGBRC(77)+ALGBRC(76)
      ALGBRC(102) =  CONSTS(10)*EXP(( - 0.820000*CONSTS(3)*STATES(1))/( CONSTS(1)*CONSTS(2)))
      ALGBRC(103) = 1.00000/(1.00000+ CONSTS(55)/ALGBRC(102) ** 1.06000*(1.00000+CONSTS(12)/CONSTS(57) ** 1.12000))
      ALGBRC(99) = 1.00000/(1.00000+ CONSTS(54)/STATES(2) ** 1.06000*(1.00000+STATES(3)/CONSTS(56) ** 1.12000))
      ALGBRC(101) = 0.370000/(1.00000+CONSTS(58)/STATES(5))
      ALGBRC(104) =  CONSTS(59)*CONSTS(4)*1.00000*( ALGBRC(101)*ALGBRC(99)*STATES(28) -  CONSTS(51)*ALGBRC(103)*(1.00000 - STATES(28)))
      ALGBRC(95) =  CONSTS(49)*CONSTS(4)*1.00000*( ALGBRC(93)*ALGBRC(91)*STATES(27) -  ALGBRC(94)*CONSTS(98)*(1.00000 - STATES(27)))
      ALGBRC(106) = ALGBRC(46)+ALGBRC(51)+ALGBRC(55)+ALGBRC(71)+ALGBRC(72)+ALGBRC(75)+ALGBRC(78)+ALGBRC(90)+ALGBRC(104)+ALGBRC(95)
      ALGBRC(107) = ALGBRC(43)+ALGBRC(74)+ALGBRC(77)+ALGBRC(49)+ALGBRC(80)+ALGBRC(85)+ 3.00000*ALGBRC(104)+ 3.00000*ALGBRC(95)
      ALGBRC(108) = (ALGBRC(71)+ALGBRC(72)+ALGBRC(76)+ALGBRC(88)+ALGBRC(73)+ALGBRC(45)+ALGBRC(50)+ALGBRC(79)+ALGBRC(84)+ALGBRC(82)) -  2.00000*ALGBRC(104)
      ALGBRC(105) = 1.00000/(1.00000+ CONSTS(57)/CONSTS(12) ** 1.12000*(1.00000+ALGBRC(102)/CONSTS(55) ** 1.06000))
      ALGBRC(110) =  CONSTS(51)*ALGBRC(103)+ CONSTS(53)*ALGBRC(105)
      ALGBRC(100) = 1.00000/(1.00000+ CONSTS(56)/STATES(3) ** 1.12000*(1.00000+STATES(2)/CONSTS(54) ** 1.06000))
      ALGBRC(112) =  ALGBRC(101)*ALGBRC(99)+ CONSTS(52)*ALGBRC(100)
      ALGBRC(111) = 1.00000/(1.00000+CONSTS(64)/ALGBRC(30))
      ALGBRC(109) = 1.00000/(1.00000+CONSTS(63)/STATES(29))
      ALGBRC(115) = 1.00000/(1.00000+CONSTS(65)/STATES(5))
      ALGBRC(116) =  CONSTS(66)*1.00000*( CONSTS(60)*ALGBRC(109)*STATES(30) -  ALGBRC(115)*ALGBRC(111)*(1.00000 - STATES(30)))
      ALGBRC(118) =  - 0.400000*STATES(35)*CONSTS(76)
      ALGBRC(114) = 1.00000 - ALGBRC(111)
      ALGBRC(117) =  ALGBRC(115)*ALGBRC(111)+ CONSTS(62)*ALGBRC(114)
      ALGBRC(113) = 1.00000 - ALGBRC(109)
      ALGBRC(119) =  CONSTS(60)*ALGBRC(109)+ CONSTS(61)*ALGBRC(113)
      ALGBRC(53) =  ALGBRC(52)*ALGBRC(47)
      ALGBRC(120) =  280000.*ALGBRC(30)/1.00000 ** 2.00000+ CONSTS(69)*ALGBRC(53)
      ALGBRC(122) = (CONSTS(74) - STATES(33))+CONSTS(75)
      ALGBRC(123) =  CONSTS(75)*STATES(33)
      ALGBRC(124) = ( (ALGBRC(122) ** 2.00000+ 4.00000*ALGBRC(123)) ** (1.0 / 2) - ALGBRC(122))/2.00000
      ALGBRC(125) = 0.0800000/(1.00000+0.360000/ALGBRC(124))
      ALGBRC(126) =  0.000377000*ALGBRC(124)/1.00000 ** 2.00000
      ALGBRC(127) =  CONSTS(67)*(ALGBRC(124) - ALGBRC(30))*STATES(31)
      ALGBRC(128) =  CONSTS(70)*(STATES(29) - ALGBRC(124))
      ALGBRC(121) =  CONSTS(71)*(STATES(29) - ALGBRC(30))
      ALGBRC(135) =  CONSTS(78)*STATES(34)*CONSTS(99) -  CONSTS(83)*STATES(35)
      ALGBRC(131) = ((1.00000 - STATES(34)) - STATES(35)) - STATES(36)
      ALGBRC(133) =  CONSTS(77)*ALGBRC(30)*ALGBRC(131) -  CONSTS(82)*STATES(34)
      ALGBRC(129) = CONSTS(86) - STATES(37)
      ALGBRC(137) =  CONSTS(79)*STATES(35) -  CONSTS(84)*STATES(36)*ALGBRC(30)
      ALGBRC(140) =  CONSTS(81)*RATES(37) ** 2.00000*STATES(35)
      ALGBRC(138) =  CONSTS(80)*STATES(36)
      ALGBRC(139) =  CONSTS(81)*RATES(37) ** 2.00000*STATES(36)
      ALGBRC(97) = (ALGBRC(48)+ALGBRC(55)+ALGBRC(89)) -  2.00000*ALGBRC(95)
      ALGBRC(141) =  CONSTS(76)*((ALGBRC(140)+ALGBRC(137)) - ALGBRC(133))
      ALGBRC(59) =  CONSTS(29)*ALGBRC(56)
      ALGBRC(130) =  CONSTS(76)*(STATES(35)+STATES(36))
      ALGBRC(132) =  CONSTS(90)*ALGBRC(130)
      ALGBRC(134) =  ALGBRC(132)*ALGBRC(129)
      ALGBRC(136) = - CONSTS(100)+ALGBRC(134)
      RETURN
      END
      REAL FUNCTION TERNRY(TEST, VALA, VALB)
      LOGICAL TEST
      REAL VALA, VALB
      IF (TEST) THEN
        TERNRY = VALA
      ELSE
        TERNRY = VALB
      ENDIF
      RETURN
      END