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 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
Source
Derived from workspace Matsuoka, Sarai, Kuratomi, Ono, Noma, 2003 at changeset dd1892d2a56a.
Collaboration
To begin collaborating on this work, please use your git client and issue this command:
License

This work is licensed under a Creative Commons Attribution 3.0 Unported License.