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 35 entries in the algebraic variable array. C There are a total of 16 entries in each of the rate and state variable arrays. C There are a total of 134 entries in the constant variable array. C C C VOI is time in component environment (second). C STATES(1) is ADP_m in component ADP_m (millimolar). C ALGBRC(33) is V_ANT in component calcium_dynamics (flux). C ALGBRC(31) is V_ATPase in component oxidative_phosphorylation (flux). C ALGBRC(12) is V_SL in component V_SL (flux). C STATES(2) is NADH in component NADH (millimolar). C ALGBRC(27) is V_O2 in component oxidative_phosphorylation (flux). C ALGBRC(10) is V_IDH in component V_IDH (flux). C ALGBRC(11) is V_KGDH in component V_KGDH (flux). C ALGBRC(14) is V_MDH in component V_MDH (flux). C STATES(3) is ISOC in component ISOC (millimolar). C ALGBRC(9) is V_ACO in component V_ACO (flux). C STATES(4) is alpha_KG in component alpha_KG (millimolar). C ALGBRC(16) is V_AAT in component V_AAT (flux). C STATES(5) is SCoA in component SCoA (millimolar). C STATES(6) is Suc in component Suc (millimolar). C ALGBRC(13) is V_SDH in component V_SDH (flux). C STATES(7) is FUM in component FUM (millimolar). C ALGBRC(17) is V_FH in component V_FH (flux). C STATES(8) is MAL in component MAL (millimolar). C STATES(9) is OAA in component OAA (millimolar). C ALGBRC(8) is V_CS in component V_CS (flux). C STATES(10) is ASP in component ASP (millimolar). C ALGBRC(19) is V_C_ASP in component V_C_ASP (flux). C STATES(11) is Ca_m in component Ca_m (micromolar). C CONSTS(1) is f in component Ca_m (dimensionless). C ALGBRC(34) is V_uni in component calcium_dynamics (flux). C ALGBRC(35) is V_NaCa in component calcium_dynamics (flux). C ALGBRC(2) is Ca_i in component Ca_i (micromolar). C CONSTS(2) is stim_start in component Ca_i (second). C CONSTS(3) is stim_end in component Ca_i (second). C ALGBRC(1) is stim_period in component Ca_i (second). C CONSTS(4) is stim_duration in component Ca_i (second). C CONSTS(5) is pulse_value in component Ca_i (micromolar). C CONSTS(6) is Na_i in component Na_i (millimolar). C CONSTS(7) is ATP_i in component ATP_i (millimolar). C ALGBRC(3) is ATP_m in component ATP_m (millimolar). C CONSTS(8) is Cm in component ATP_m (millimolar). C ALGBRC(4) is ADP_i in component ADP_i (millimolar). C CONSTS(9) is pulse_value in component ADP_i (millimolar). C CONSTS(10) is GLU in component GLU (millimolar). C CONSTS(11) is Mg in component Mg (millimolar). C CONSTS(12) is H in component H (millimolar). C CONSTS(13) is Pi in component Pi (millimolar). C CONSTS(14) is CoA in component CoA (millimolar). C CONSTS(15) is AcCoA in component AcCoA (millimolar). C CONSTS(16) is FAD in component FAD (millimolar). C CONSTS(17) is FADH2 in component FADH2 (millimolar). C ALGBRC(5) is NAD in component NAD (millimolar). C CONSTS(18) is C_PN in component NAD (millimolar). C CONSTS(19) is NADPH in component NADPH (millimolar). C STATES(12) is O2_m in component O2_m (millimolar). C CONSTS(20) is shunt in component O2_m (dimensionless). C ALGBRC(21) is VTr_ROS in component VTr_ROS (flux). C STATES(13) is O2_i in component O2_i (millimolar). C ALGBRC(15) is V_SOD in component V_SOD (flux). C STATES(14) is H2O2 in component H2O2 (millimolar). C ALGBRC(20) is V_CAT in component V_CAT (flux). C ALGBRC(22) is V_GPX in component V_GPX (flux). C STATES(15) is GSH in component GSH (millimolar). C ALGBRC(24) is V_GR in component V_GR (flux). C ALGBRC(6) is GSSG in component GSSG (millimolar). C CONSTS(21) is GT in component GSSG (millimolar). C ALGBRC(7) is CIT in component CIT (millimolar). C CONSTS(22) is C_Kint in component CIT (millimolar). C STATES(16) is delta_psi_m in component mitochondrial_membrane (volt). C CONSTS(23) is R in component mitochondrial_membrane (volt_coulomb_per_mole_kelvin). C CONSTS(24) is T in component mitochondrial_membrane (kelvin). C CONSTS(25) is F in component mitochondrial_membrane (coulomb_per_mole). C CONSTS(26) is C_mito in component mitochondrial_membrane (millimolar_per_volt). C ALGBRC(26) is V_He in component oxidative_phosphorylation (flux). C ALGBRC(28) is V_He_F in component oxidative_phosphorylation (flux). C ALGBRC(32) is V_Hu in component oxidative_phosphorylation (flux). C ALGBRC(29) is V_Hleak in component oxidative_phosphorylation (flux). C CONSTS(27) is Km_AcCoA in component V_CS (millimolar). C CONSTS(28) is Km_OAA in component V_CS (millimolar). C CONSTS(29) is Kcat_CS in component V_CS (first_order_rate_constant). C CONSTS(30) is ET_CS in component V_CS (millimolar). C CONSTS(31) is Kf_ACO in component V_ACO (first_order_rate_constant). C CONSTS(32) is KE_ACO in component V_ACO (dimensionless). C CONSTS(33) is Kh_1 in component V_IDH (millimolar). C CONSTS(34) is Kh_2 in component V_IDH (millimolar). C CONSTS(35) is Km_ISOC in component V_IDH (millimolar). C CONSTS(36) is Ka_ADP in component V_IDH (millimolar). C CONSTS(37) is Ka_Ca in component V_IDH (micromolar). C CONSTS(38) is Km_NAD in component V_IDH (millimolar). C CONSTS(39) is Ki_NADH in component V_IDH (millimolar). C CONSTS(40) is Kcat_IDH in component V_IDH (first_order_rate_constant). C CONSTS(41) is ET_IDH in component V_IDH (millimolar). C CONSTS(42) is ni in component V_IDH (dimensionless). C CONSTS(43) is Km_alpha_KG in component V_KGDH (millimolar). C CONSTS(44) is Kcat_KGDH in component V_KGDH (first_order_rate_constant). C CONSTS(45) is ET_KGDH in component V_KGDH (millimolar). C CONSTS(46) is Kd_Mg in component V_KGDH (millimolar). C CONSTS(47) is Kd_Ca in component V_KGDH (micromolar). C CONSTS(48) is n_alpha_KG in component V_KGDH (dimensionless). C CONSTS(49) is Km_NAD in component V_KGDH (millimolar). C CONSTS(50) is kf_SL in component V_SL (second_order_rate_constant). C CONSTS(51) is Ke_SL in component V_SL (millimolar). C CONSTS(52) is Kisdh_OAA in component V_SDH (millimolar). C CONSTS(53) is Kcat_SDH in component V_SDH (first_order_rate_constant). C CONSTS(54) is ET_SDH in component V_SDH (millimolar). C CONSTS(55) is Km_Suc in component V_SDH (millimolar). C CONSTS(56) is Ki_FUM in component V_SDH (millimolar). C CONSTS(57) is Km_MAL in component V_MDH (millimolar). C CONSTS(58) is Kcat_MDH in component V_MDH (first_order_rate_constant). C CONSTS(59) is ET_MDH in component V_MDH (millimolar). C CONSTS(60) is Ki_OAA in component V_MDH (millimolar). C CONSTS(133) is fh_a in component V_MDH (dimensionless). C CONSTS(134) is fh_i in component V_MDH (dimensionless). C CONSTS(61) is Km_NAD in component V_MDH (millimolar). C CONSTS(62) is kh1 in component V_MDH (millimolar). C CONSTS(63) is kh2 in component V_MDH (millimolar). C CONSTS(64) is kh3 in component V_MDH (millimolar). C CONSTS(65) is kh4 in component V_MDH (millimolar). C CONSTS(66) is k_offset in component V_MDH (dimensionless). C CONSTS(67) is Ke_FH in component V_FH (dimensionless). C CONSTS(68) is kf_FH in component V_FH (first_order_rate_constant). C CONSTS(69) is Ke_AAT in component V_AAT (dimensionless). C CONSTS(70) is kf_AAT in component V_AAT (second_order_rate_constant). C CONSTS(71) is k_C_ASP in component V_C_ASP (first_order_rate_constant). C CONSTS(72) is k1_SOD in component V_SOD (second_order_rate_constant). C CONSTS(73) is k3_SOD in component V_SOD (second_order_rate_constant). C CONSTS(74) is k5_SOD in component V_SOD (first_order_rate_constant). C CONSTS(75) is ET_SOD in component V_SOD (millimolar). C CONSTS(76) is Ki_H2O2 in component V_SOD (millimolar). C CONSTS(77) is k1_CAT in component V_CAT (second_order_rate_constant). C CONSTS(78) is ET_CAT in component V_CAT (millimolar). C CONSTS(79) is fr in component V_CAT (per_millimolar). C CONSTS(80) is phi1_GPX in component V_GPX (millimolar_second). C CONSTS(81) is phi2_GPX in component V_GPX (millimolar_second). C CONSTS(82) is ET_GPX in component V_GPX (millimolar). C CONSTS(83) is k1_GR in component V_GR (first_order_rate_constant). C CONSTS(84) is KM_GSSG in component V_GR (millimolar). C CONSTS(85) is KM_NADPH in component V_GR (millimolar). C CONSTS(86) is ET_GR in component V_GR (millimolar). C ALGBRC(18) is V_IMAC in component V_IMAC (flux). C CONSTS(87) is Gmax in component V_IMAC (millimolar_per_second_volt). C CONSTS(88) is GL in component V_IMAC (millimolar_per_second_volt). C CONSTS(89) is a in component V_IMAC (dimensionless). C CONSTS(90) is b in component V_IMAC (dimensionless). C CONSTS(91) is kappa in component V_IMAC (per_volt). C CONSTS(92) is delta_psi_bm in component V_IMAC (volt). C CONSTS(93) is Kcc in component V_IMAC (millimolar). C CONSTS(94) is j in component VTr_ROS (dimensionless). C CONSTS(95) is rho_res in component oxidative_phosphorylation (millimolar). C CONSTS(96) is rho_res_F in component oxidative_phosphorylation (millimolar). C CONSTS(97) is ra in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(98) is rc1 in component oxidative_phosphorylation (first_order_rate_constant). C ALGBRC(23) is Ares in component oxidative_phosphorylation (volt). C CONSTS(132) is Ares_F in component oxidative_phosphorylation (volt). C CONSTS(99) is r1 in component oxidative_phosphorylation (dimensionless). C CONSTS(100) is r2 in component oxidative_phosphorylation (dimensionless). C CONSTS(101) is r3 in component oxidative_phosphorylation (dimensionless). C CONSTS(102) is rb in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(103) is rc2 in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(104) is Kres in component oxidative_phosphorylation (dimensionless). C CONSTS(105) is Kres_F in component oxidative_phosphorylation (dimensionless). C CONSTS(106) is gH in component oxidative_phosphorylation (millimolar_per_second_volt). C CONSTS(107) is delta_psi_B in component oxidative_phosphorylation (volt). C CONSTS(108) is g in component oxidative_phosphorylation (dimensionless). C ALGBRC(25) is delta_mu_H in component oxidative_phosphorylation (volt). C CONSTS(109) is delta_pH in component oxidative_phosphorylation (dimensionless). C CONSTS(110) is rho_F1 in component oxidative_phosphorylation (millimolar). C CONSTS(111) is pa in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(112) is pc1 in component oxidative_phosphorylation (first_order_rate_constant). C ALGBRC(30) is AF1 in component oxidative_phosphorylation (volt). C CONSTS(113) is p1 in component oxidative_phosphorylation (dimensionless). C CONSTS(114) is p2 in component oxidative_phosphorylation (dimensionless). C CONSTS(115) is p3 in component oxidative_phosphorylation (dimensionless). C CONSTS(116) is pb in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(117) is pc2 in component oxidative_phosphorylation (first_order_rate_constant). C CONSTS(118) is KF1 in component oxidative_phosphorylation (millimolar). C CONSTS(119) is h in component calcium_dynamics (dimensionless). C CONSTS(120) is delta_psi_0 in component calcium_dynamics (volt). C CONSTS(121) is Vmax_ANT in component calcium_dynamics (flux). C CONSTS(122) is L in component calcium_dynamics (dimensionless). C CONSTS(123) is na in component calcium_dynamics (dimensionless). C CONSTS(124) is Vmax_uni in component calcium_dynamics (flux). C CONSTS(125) is K_act in component calcium_dynamics (micromolar). C CONSTS(126) is K_trans in component calcium_dynamics (micromolar). C CONSTS(127) is n in component calcium_dynamics (dimensionless). C CONSTS(128) is Vmax_NaCa in component calcium_dynamics (flux). C CONSTS(129) is KNa in component calcium_dynamics (millimolar). C CONSTS(130) is KCa in component calcium_dynamics (micromolar). C CONSTS(131) is b in component calcium_dynamics (dimensionless). C RATES(1) is d/dt ADP_m in component ADP_m (millimolar). C RATES(2) is d/dt NADH in component NADH (millimolar). C RATES(3) is d/dt ISOC in component ISOC (millimolar). C RATES(4) is d/dt alpha_KG in component alpha_KG (millimolar). C RATES(5) is d/dt SCoA in component SCoA (millimolar). C RATES(6) is d/dt Suc in component Suc (millimolar). C RATES(7) is d/dt FUM in component FUM (millimolar). C RATES(8) is d/dt MAL in component MAL (millimolar). C RATES(9) is d/dt OAA in component OAA (millimolar). C RATES(10) is d/dt ASP in component ASP (millimolar). C RATES(11) is d/dt Ca_m in component Ca_m (micromolar). C RATES(12) is d/dt O2_m in component O2_m (millimolar). C RATES(13) is d/dt O2_i in component O2_i (millimolar). C RATES(14) is d/dt H2O2 in component H2O2 (millimolar). C RATES(15) is d/dt GSH in component GSH (millimolar). C RATES(16) is d/dt delta_psi_m in component mitochondrial_membrane (volt). C SUBROUTINE initConsts(CONSTS, RATES, STATES) REAL CONSTS(*), RATES(*), STATES(*) STATES(1) = 0.1 STATES(2) = 0.01 STATES(3) = 0.01 STATES(4) = 0.01 STATES(5) = 0.01 STATES(6) = 0.01 STATES(7) = 0.01 STATES(8) = 0.01 STATES(9) = 0.01 STATES(10) = 0.01 STATES(11) = 0.01 CONSTS(1) = 0.0003 CONSTS(2) = 0 CONSTS(3) = 10000 CONSTS(4) = 0.4 CONSTS(5) = 1 CONSTS(6) = 10.0 CONSTS(7) = 6.5 CONSTS(8) = 15.0 CONSTS(9) = 0.15 CONSTS(10) = 20 CONSTS(11) = 0.4 CONSTS(12) = 2.5E-5 CONSTS(13) = 20.0 CONSTS(14) = 0.02 CONSTS(15) = 0.0002 CONSTS(16) = 0.01 CONSTS(17) = 1.24 CONSTS(18) = 10.0 CONSTS(19) = 1.0 STATES(12) = 0.01 CONSTS(20) = 0.05 STATES(13) = 0.01 STATES(14) = 0.01 STATES(15) = 0.01 CONSTS(21) = 1.0 CONSTS(22) = 1.0 STATES(16) = 0.01 CONSTS(23) = 8.315 CONSTS(24) = 310.16 CONSTS(25) = 96480 CONSTS(26) = 1.812 CONSTS(27) = 1.26E-2 CONSTS(28) = 6.4E-4 CONSTS(29) = 3.2 CONSTS(30) = 0.4 CONSTS(31) = 12.5 CONSTS(32) = 2.22 CONSTS(33) = 8.1E-5 CONSTS(34) = 5.98E-5 CONSTS(35) = 1.52 CONSTS(36) = 6.2E-2 CONSTS(37) = 1.41 CONSTS(38) = 0.923 CONSTS(39) = 0.19 CONSTS(40) = 1.94 CONSTS(41) = 0.109 CONSTS(42) = 1 CONSTS(43) = 1.94 CONSTS(44) = 0.15 CONSTS(45) = 0.5 CONSTS(46) = 0.0308 CONSTS(47) = 1.27 CONSTS(48) = 1.2 CONSTS(49) = 38.7 CONSTS(50) = 0.127 CONSTS(51) = 3.115 CONSTS(52) = 0.15 CONSTS(53) = 1.0 CONSTS(54) = 0.5 CONSTS(55) = 3.0E-2 CONSTS(56) = 1.3 CONSTS(57) = 1.493 CONSTS(58) = 2.775E1 CONSTS(59) = 0.154 CONSTS(60) = 3.1E-3 CONSTS(61) = 0.2244 CONSTS(62) = 1.13E-5 CONSTS(63) = 26.7 CONSTS(64) = 6.68E-9 CONSTS(65) = 5.62E-6 CONSTS(66) = 3.99E-2 CONSTS(67) = 1.0 CONSTS(68) = 0.83 CONSTS(69) = 6.6 CONSTS(70) = 0.644 CONSTS(71) = 0.01 CONSTS(72) = 2.4E6 CONSTS(73) = 4.8E4 CONSTS(74) = 5.0E-1 CONSTS(75) = 0.4E-3 CONSTS(76) = 0.5 CONSTS(77) = 1.7E4 CONSTS(78) = 0.001 CONSTS(79) = 50.0 CONSTS(80) = 0.15 CONSTS(81) = 0.5 CONSTS(82) = 0.00141 CONSTS(83) = 0.0308 CONSTS(84) = 1.94 CONSTS(85) = 38.7 CONSTS(86) = 1.27E-3 CONSTS(87) = 7.82 CONSTS(88) = 0.0782 CONSTS(89) = 1.0E-3 CONSTS(90) = 1.0E4 CONSTS(91) = 70.0 CONSTS(92) = 0.004 CONSTS(93) = 0.01 CONSTS(94) = 0.12 CONSTS(95) = 0.0006 CONSTS(96) = 0.0045 CONSTS(97) = 6.394E-10 CONSTS(98) = 2.656E-19 CONSTS(99) = 2.077E-18 CONSTS(100) = 1.728E-9 CONSTS(101) = 1.059E-26 CONSTS(102) = 1.762E-13 CONSTS(103) = 8.632E-27 CONSTS(104) = 1.35E18 CONSTS(105) = 5.765E13 CONSTS(106) = 0.01 CONSTS(107) = 0.05 CONSTS(108) = 0.85 CONSTS(109) = -0.6 CONSTS(110) = 0.06 CONSTS(111) = 1.656E-5 CONSTS(112) = 9.651E-14 CONSTS(113) = 1.346E-8 CONSTS(114) = 7.739E-7 CONSTS(115) = 6.65E-15 CONSTS(116) = 3.373E-7 CONSTS(117) = 4.585E-14 CONSTS(118) = 1.71E6 CONSTS(119) = 0.5 CONSTS(120) = 0.091 CONSTS(121) = 0.05 CONSTS(122) = 110.0 CONSTS(123) = 2.8 CONSTS(124) = 0.000625 CONSTS(125) = 0.38 CONSTS(126) = 19.0 CONSTS(127) = 3.0 CONSTS(128) = 0.005 CONSTS(129) = 9.4 CONSTS(130) = 3.75E-1 CONSTS(131) = 0.5 CONSTS(132) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(105)*CONSTS(17)/CONSTS(16) ** 0.500000) CONSTS(133) = 1.00000/(1.00000+CONSTS(12)/CONSTS(62)+CONSTS(12) ** 2.00000/( CONSTS(62)*CONSTS(63)))+CONSTS(66) CONSTS(134) = 1.00000/(1.00000+CONSTS(64)/CONSTS(12)+( CONSTS(64)*CONSTS(65))/CONSTS(12) ** 2.00000) ** 2.00000 RETURN END SUBROUTINE computeRates(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(5) = CONSTS(18) - STATES(2) ALGBRC(10) = ( CONSTS(40)*CONSTS(41))/(1.00000+CONSTS(12)/CONSTS(33)+CONSTS(34)/CONSTS(12)+CONSTS(35)/STATES(3) ** CONSTS(42)/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))+ (CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39))+( CONSTS(35)/STATES(3) ** CONSTS(42)*(CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39)))/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))) ALGBRC(7) = CONSTS(22) - (STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)) ALGBRC(9) = CONSTS(31)*(ALGBRC(7) - STATES(3)/CONSTS(32)) RATES(3) = ALGBRC(9) - ALGBRC(10) ALGBRC(3) = CONSTS(8) - STATES(1) ALGBRC(12) = CONSTS(50)*( STATES(5)*STATES(1) - ( STATES(6)*ALGBRC(3)*CONSTS(14))/CONSTS(51)) ALGBRC(11) = ( CONSTS(44)*CONSTS(45))/(1.00000+CONSTS(43)/STATES(4) ** CONSTS(48)/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))+(CONSTS(49)/ALGBRC(5))/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))) RATES(5) = ALGBRC(11) - ALGBRC(12) ALGBRC(13) = ( CONSTS(53)*CONSTS(54))/(1.00000+ (CONSTS(55)/STATES(6))*(1.00000+STATES(9)/CONSTS(52))*(1.00000+STATES(7)/CONSTS(56))) RATES(6) = ALGBRC(12) - ALGBRC(13) ALGBRC(16) = CONSTS(70)*( STATES(9)*CONSTS(10) - ( STATES(4)*STATES(10))/CONSTS(69)) RATES(4) = ALGBRC(16)+ALGBRC(10)+- ALGBRC(11) ALGBRC(17) = CONSTS(68)*(STATES(7) - STATES(8)/CONSTS(67)) RATES(7) = ALGBRC(13) - ALGBRC(17) ALGBRC(14) = ( CONSTS(58)*CONSTS(59)*CONSTS(133)*CONSTS(134))/(1.00000+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))+CONSTS(61)/ALGBRC(5)+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))*(CONSTS(61)/ALGBRC(5))) RATES(8) = ALGBRC(17) - ALGBRC(14) ALGBRC(8) = ( CONSTS(29)*CONSTS(30))/(1.00000+CONSTS(27)/CONSTS(15)+CONSTS(28)/STATES(9)+ (CONSTS(27)/CONSTS(15))*(CONSTS(28)/STATES(9))) RATES(9) = ALGBRC(14) - (ALGBRC(8)+ALGBRC(16)) ALGBRC(19) = CONSTS(71)*STATES(10) RATES(10) = ALGBRC(16) - ALGBRC(19) ALGBRC(18) = (CONSTS(89)+CONSTS(90)/(1.00000+CONSTS(93)/STATES(13)))*(CONSTS(88)+CONSTS(87)/(1.00000+EXP( CONSTS(91)*(CONSTS(92) - STATES(16)))))*STATES(16) ALGBRC(21) = CONSTS(94)*(ALGBRC(18)/STATES(16))*(STATES(16) - (( CONSTS(23)*CONSTS(24))/CONSTS(25))*arbitrary_log(STATES(12)/STATES(13), 10)) ALGBRC(15) = ( 2.00000*CONSTS(72)*CONSTS(74)*(CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))*CONSTS(75)*STATES(13))/( CONSTS(74)*( 2.00000*CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))+ STATES(13)*CONSTS(72)*CONSTS(73)*(1.00000+STATES(14)/CONSTS(76))) RATES(13) = ALGBRC(21) - ALGBRC(15) ALGBRC(20) = 2.00000*CONSTS(77)*CONSTS(78)*STATES(14)*EXP( - CONSTS(79)*STATES(14)) ALGBRC(22) = ( CONSTS(82)*STATES(14)*STATES(15))/( CONSTS(80)*STATES(15)+ CONSTS(81)*STATES(14)) RATES(14) = ALGBRC(15) - (ALGBRC(20)+ALGBRC(22)) ALGBRC(6) = 2.00000*(CONSTS(21) - STATES(15)) ALGBRC(24) = ( CONSTS(83)*CONSTS(86))/(1.00000+CONSTS(84)/ALGBRC(6)+CONSTS(85)/CONSTS(19)+ (CONSTS(84)/ALGBRC(6))*(CONSTS(85)/CONSTS(19))) RATES(15) = ALGBRC(24) - ALGBRC(22) ALGBRC(23) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(104)*STATES(2)/ALGBRC(5) ** 0.500000) ALGBRC(25) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*CONSTS(109)+STATES(16) ALGBRC(27) = 0.500000*CONSTS(95)*((( (CONSTS(97)+ CONSTS(98)*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - CONSTS(97)*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))+ CONSTS(103)*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) RATES(2) = - ALGBRC(27)+ALGBRC(10)+ALGBRC(11)+ALGBRC(14) RATES(12) = CONSTS(20)*ALGBRC(27) - ALGBRC(21) ALGBRC(4) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, CONSTS(9), 0.0500000) ALGBRC(33) = CONSTS(121)*((1.00000 - ( 0.0500000*CONSTS(7)*0.450000*0.800000*STATES(1))/( 0.450000*ALGBRC(4)*0.0500000*ALGBRC(3)))/( (1.00000+ (( 0.0500000*CONSTS(7))/( 0.450000*ALGBRC(4)))*EXP(( - CONSTS(119)*CONSTS(25)*CONSTS(120))/( CONSTS(23)*CONSTS(24))))*(1.00000+( 0.450000*0.800000*STATES(1))/( 0.0500000*ALGBRC(3))))) ALGBRC(30) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(118)*(ALGBRC(3)/( STATES(1)*CONSTS(13)))) ALGBRC(31) = - CONSTS(110)*(( ( 100.000*CONSTS(111)+ CONSTS(112)*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - ( CONSTS(111)*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))+ CONSTS(117)*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) RATES(1) = ALGBRC(33) - (ALGBRC(31)+ALGBRC(12)) ALGBRC(1) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, 0.500000, 4.00000) ALGBRC(2) = TERNRY(VOI.GE.CONSTS(2).AND.VOI.LE.CONSTS(3).AND.(VOI - CONSTS(2)) - INT((VOI - CONSTS(2))/ALGBRC(1))*ALGBRC(1).LE.CONSTS(4), CONSTS(5), 0.100000) ALGBRC(34) = CONSTS(124)*(( (ALGBRC(2)/CONSTS(126))*1.00000+ALGBRC(2)/CONSTS(126) ** 3.00000*(( 2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))/(1.00000+ALGBRC(2)/CONSTS(126) ** 4.00000+ (CONSTS(122)/1.00000+ALGBRC(2)/CONSTS(125) ** CONSTS(123))*(1.00000 - EXP(( -2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))))) ALGBRC(35) = CONSTS(128)*(( EXP(( CONSTS(131)*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))*EXP(log(ALGBRC(2)/STATES(11))))/( 1.00000+CONSTS(129)/CONSTS(6) ** CONSTS(127)*(1.00000+CONSTS(130)/STATES(11)))) RATES(11) = CONSTS(1)*1.00000*(ALGBRC(34) - ALGBRC(35)) ALGBRC(26) = 6.00000*CONSTS(95)*(( CONSTS(97)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))) - (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(28) = 6.00000*CONSTS(96)*(( CONSTS(97)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))) - (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(32) = -3.00000*CONSTS(110)*(( 100.000*CONSTS(111)*(1.00000+EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24)))) - (CONSTS(111)+CONSTS(116))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(29) = CONSTS(106)*ALGBRC(25) RATES(16) = (ALGBRC(26)+ALGBRC(28)+- (ALGBRC(32)+ALGBRC(33)+ALGBRC(29)+ALGBRC(35)+ 2.00000*ALGBRC(34)))/CONSTS(26) RETURN END SUBROUTINE computeVariables(VOI, CONSTS, RATES, STATES, ALGBRC) REAL VOI, CONSTS(*), RATES(*), STATES(*), ALGBRC(*) ALGBRC(5) = CONSTS(18) - STATES(2) ALGBRC(10) = ( CONSTS(40)*CONSTS(41))/(1.00000+CONSTS(12)/CONSTS(33)+CONSTS(34)/CONSTS(12)+CONSTS(35)/STATES(3) ** CONSTS(42)/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))+ (CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39))+( CONSTS(35)/STATES(3) ** CONSTS(42)*(CONSTS(38)/ALGBRC(5))*(1.00000+STATES(2)/CONSTS(39)))/( (1.00000+STATES(1)/CONSTS(36))*(1.00000+STATES(11)/CONSTS(37)))) ALGBRC(7) = CONSTS(22) - (STATES(3)+STATES(4)+STATES(5)+STATES(6)+STATES(7)+STATES(8)+STATES(9)) ALGBRC(9) = CONSTS(31)*(ALGBRC(7) - STATES(3)/CONSTS(32)) ALGBRC(3) = CONSTS(8) - STATES(1) ALGBRC(12) = CONSTS(50)*( STATES(5)*STATES(1) - ( STATES(6)*ALGBRC(3)*CONSTS(14))/CONSTS(51)) ALGBRC(11) = ( CONSTS(44)*CONSTS(45))/(1.00000+CONSTS(43)/STATES(4) ** CONSTS(48)/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))+(CONSTS(49)/ALGBRC(5))/( (1.00000+CONSTS(11)/CONSTS(46))*(1.00000+STATES(11)/CONSTS(47)))) ALGBRC(13) = ( CONSTS(53)*CONSTS(54))/(1.00000+ (CONSTS(55)/STATES(6))*(1.00000+STATES(9)/CONSTS(52))*(1.00000+STATES(7)/CONSTS(56))) ALGBRC(16) = CONSTS(70)*( STATES(9)*CONSTS(10) - ( STATES(4)*STATES(10))/CONSTS(69)) ALGBRC(17) = CONSTS(68)*(STATES(7) - STATES(8)/CONSTS(67)) ALGBRC(14) = ( CONSTS(58)*CONSTS(59)*CONSTS(133)*CONSTS(134))/(1.00000+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))+CONSTS(61)/ALGBRC(5)+ (CONSTS(57)/STATES(8))*(1.00000+STATES(9)/CONSTS(60))*(CONSTS(61)/ALGBRC(5))) ALGBRC(8) = ( CONSTS(29)*CONSTS(30))/(1.00000+CONSTS(27)/CONSTS(15)+CONSTS(28)/STATES(9)+ (CONSTS(27)/CONSTS(15))*(CONSTS(28)/STATES(9))) ALGBRC(19) = CONSTS(71)*STATES(10) ALGBRC(18) = (CONSTS(89)+CONSTS(90)/(1.00000+CONSTS(93)/STATES(13)))*(CONSTS(88)+CONSTS(87)/(1.00000+EXP( CONSTS(91)*(CONSTS(92) - STATES(16)))))*STATES(16) ALGBRC(21) = CONSTS(94)*(ALGBRC(18)/STATES(16))*(STATES(16) - (( CONSTS(23)*CONSTS(24))/CONSTS(25))*arbitrary_log(STATES(12)/STATES(13), 10)) ALGBRC(15) = ( 2.00000*CONSTS(72)*CONSTS(74)*(CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))*CONSTS(75)*STATES(13))/( CONSTS(74)*( 2.00000*CONSTS(72)+ CONSTS(73)*(1.00000+STATES(14)/CONSTS(76)))+ STATES(13)*CONSTS(72)*CONSTS(73)*(1.00000+STATES(14)/CONSTS(76))) ALGBRC(20) = 2.00000*CONSTS(77)*CONSTS(78)*STATES(14)*EXP( - CONSTS(79)*STATES(14)) ALGBRC(22) = ( CONSTS(82)*STATES(14)*STATES(15))/( CONSTS(80)*STATES(15)+ CONSTS(81)*STATES(14)) ALGBRC(6) = 2.00000*(CONSTS(21) - STATES(15)) ALGBRC(24) = ( CONSTS(83)*CONSTS(86))/(1.00000+CONSTS(84)/ALGBRC(6)+CONSTS(85)/CONSTS(19)+ (CONSTS(84)/ALGBRC(6))*(CONSTS(85)/CONSTS(19))) ALGBRC(23) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(104)*STATES(2)/ALGBRC(5) ** 0.500000) ALGBRC(25) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*CONSTS(109)+STATES(16) ALGBRC(27) = 0.500000*CONSTS(95)*((( (CONSTS(97)+ CONSTS(98)*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - CONSTS(97)*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))+ CONSTS(103)*EXP(( ALGBRC(23)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(4) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, CONSTS(9), 0.0500000) ALGBRC(33) = CONSTS(121)*((1.00000 - ( 0.0500000*CONSTS(7)*0.450000*0.800000*STATES(1))/( 0.450000*ALGBRC(4)*0.0500000*ALGBRC(3)))/( (1.00000+ (( 0.0500000*CONSTS(7))/( 0.450000*ALGBRC(4)))*EXP(( - CONSTS(119)*CONSTS(25)*CONSTS(120))/( CONSTS(23)*CONSTS(24))))*(1.00000+( 0.450000*0.800000*STATES(1))/( 0.0500000*ALGBRC(3))))) ALGBRC(30) = (( CONSTS(23)*CONSTS(24))/CONSTS(25))*log( CONSTS(118)*(ALGBRC(3)/( STATES(1)*CONSTS(13)))) ALGBRC(31) = - CONSTS(110)*(( ( 100.000*CONSTS(111)+ CONSTS(112)*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24))))*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24))) - ( CONSTS(111)*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))+ CONSTS(117)*EXP(( ALGBRC(30)*CONSTS(25))/( CONSTS(23)*CONSTS(24)))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24)))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(1) = TERNRY(VOI.GE.100.000.AND.VOI.LT.300.000, 0.500000, 4.00000) ALGBRC(2) = TERNRY(VOI.GE.CONSTS(2).AND.VOI.LE.CONSTS(3).AND.(VOI - CONSTS(2)) - INT((VOI - CONSTS(2))/ALGBRC(1))*ALGBRC(1).LE.CONSTS(4), CONSTS(5), 0.100000) ALGBRC(34) = CONSTS(124)*(( (ALGBRC(2)/CONSTS(126))*1.00000+ALGBRC(2)/CONSTS(126) ** 3.00000*(( 2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24))))/(1.00000+ALGBRC(2)/CONSTS(126) ** 4.00000+ (CONSTS(122)/1.00000+ALGBRC(2)/CONSTS(125) ** CONSTS(123))*(1.00000 - EXP(( -2.00000*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))))) ALGBRC(35) = CONSTS(128)*(( EXP(( CONSTS(131)*CONSTS(25)*(STATES(16) - CONSTS(120)))/( CONSTS(23)*CONSTS(24)))*EXP(log(ALGBRC(2)/STATES(11))))/( 1.00000+CONSTS(129)/CONSTS(6) ** CONSTS(127)*(1.00000+CONSTS(130)/STATES(11)))) ALGBRC(26) = 6.00000*CONSTS(95)*(( CONSTS(97)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))) - (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*ALGBRC(23))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(28) = 6.00000*CONSTS(96)*(( CONSTS(97)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))) - (CONSTS(97)+CONSTS(102))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(99)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( 6.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(100)+ CONSTS(101)*EXP(( CONSTS(25)*CONSTS(132))/( CONSTS(23)*CONSTS(24))))*EXP(( CONSTS(108)*6.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(32) = -3.00000*CONSTS(110)*(( 100.000*CONSTS(111)*(1.00000+EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24)))) - (CONSTS(111)+CONSTS(116))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))/( (1.00000+ CONSTS(113)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*CONSTS(107))/( CONSTS(23)*CONSTS(24)))+ (CONSTS(114)+ CONSTS(115)*EXP(( CONSTS(25)*ALGBRC(30))/( CONSTS(23)*CONSTS(24))))*EXP(( 3.00000*CONSTS(25)*ALGBRC(25))/( CONSTS(23)*CONSTS(24))))) ALGBRC(29) = CONSTS(106)*ALGBRC(25) RETURN END REAL FUNCTION TERNRY(TEST, VALA, VALB) LOGICAL TEST REAL VALA, VALB IF (TEST) THEN TERNRY = VALA ELSE TERNRY = VALB ENDIF RETURN END