[1] "2025-05-13 13:43:28 MDT" options [1,] "file" "151002c" [2,] "flight" "RF05" [3,] "project" "ARISTO2015" [4,] "png" "TRUE" [5,] "valveskip" "60" [6,] "wtwtco" "200" [7,] "spspco" "500" [8,] "o2gco" "500" [9,] "co2gco" "50" [10,] "hr" "TRUE" [11,] "newskip" "0" [12,] "recalcpd" "TRUE" [13,] "mrgdir" "/scr/raf/stephens/HUMERGE" [14,] "datdir" "/scr/raf/stephens/AO2DAT" [15,] "netdir" "/net/www/docs/homes/stephens/AO2" [16,] "acdatdir" "/scr/raf/Prod_Data" [17,] "logfile" "TRUE" [18,] "uwyoflag" "none" [19,] "tlagff0" "0" [20,] "tlagff1" "0" [21,] "tlagopt" "FALSE" [22,] "useoptlag" "FALSE" [23,] "optto" "oms" [24,] "acavail" "TRUE" [25,] "huavail" "FALSE" [26,] "vxavail" "FALSE" [27,] "usetlagcalc" "TRUE" [28,] "inlet" "himl" [29,] "pcabcorr" "FALSE" [30,] "pcabforce" "ave,ave" [31,] "pcabfiltfirst" "TRUE" [32,] "pao2corr" "FALSE" [33,] "pao2force" "ave,ave" [34,] "pao2filtfirst" "TRUE" [35,] "tpboxcorr" "FALSE" [36,] "tpboxforce" "ave,ave" [37,] "tpboxfiltfirst" "TRUE" [38,] "inpcorr" "FALSE" [39,] "inpopt" "FALSE" [40,] "inpvals" "46,17,6" [41,] "inpsteps" "2,2,2" [42,] "actype" "nc2asc" [43,] "outputdesc" "250508" [44,] "readmergenew" "FALSE" [45,] "c13cor" "TRUE" [46,] "outdir" "." [47,] "co2caltype" "full" [48,] "o2caltype" "full" [49,] "lih2oco" "-10" [50,] "nasasmth" "FALSE" [51,] "o2slpcorr" "FALSE" [52,] "co2slpcorr" "FALSE" [53,] "co2usewt" "TRUE" [54,] "skipeffects" "FALSE" [55,] "hrlinredo" "FALSE" [56,] "jogsdfilt" "FALSE" [57,] "jogsdfact" "2" [58,] "use120sec" "FALSE" [59,] "lowerdec" "TRUE" [60,] "agg10sec" "FALSE" [61,] "plotxsect" "FALSE" [62,] "reavecals" "FALSE" [63,] "locate" "FALSE" [64,] "icartt" "0" [65,] "icarttrev" "0" [66,] "mmsrev" "TRUE" [67,] "usespline" "4pt" [68,] "updnave" "151002c" [1] "Reading in flight and plotting intervals" Read 23 records [1] "Reading in /scr/raf/stephens/AO2DAT/151002c.mr" Read 7434 records [1] "Converting statout" [1] "Median(diff(mrdata$time) = 2.59 seconds" [1] "Median(diff(mrdata$time,lag=240)/240) = 2.62 seconds" [1] "New JOGEST = 2.62 seconds" [1] "New VALVE SKIP = 23 jogs" [1] "Reading in /scr/raf/stephens/AO2DAT/151002c.hr" Read 29739 records [1] "HR FILE MATCHES" [1] "Converting hrdata statout" [1] "Tue May 13 01:43:29 PM MDT 2025 done" [1] "Reading in cylinder log" Read 2575 records [1] "Finding cylinder IDs" Read 1 record Read 536 records [1] "position 1 cylinder OK30771 description: HS CO2 = 391.34 O2 = -452.7" [1] "position 2 cylinder OK30764 description: LS CO2 = 404.02 O2 = -717.6" [1] "position 3 cylinder OK30772 description: LT CO2 = 392.44 O2 = -584.5" [1] "position 4 cylinder OK30768 description: WT CO2 = 394.39 O2 = -580.8" [1] "position 5 cylinder OK30761 description: LP CO2 = 369.41 O2 = -827.3" [1] 391.3386 404.0181 392.4445 394.3873 369.4096 [1] -452.7198 -717.5633 -584.4989 -580.7611 -827.3404 [1] "Nominal linear calibrations:" [1] "co2scl = 1.068 co2off = -11.73 o2scla = -0.109 o2offa = -184.2" [1] "(dV/V)/(dXO2/XO2) = -1.923" [1] "(dV/V)/d(O2/N2) = -1.52" [1] "Reading in very high resolution timeseries" Read 148680 records [1] 148680 [1] "VHR FILE MATCHES" [1] "Recalculating pdo2 and pdref corrections" [1] "Tue May 13 01:43:29 PM MDT 2025 done approxing mr variables" [1] "Tue May 13 01:43:30 PM MDT 2025 done undoing flight corrections" [1] "o2scaleavenew = -1.52" [1] "o2scaleavenew_drelVdrelP = -1.923" [1] "median(o2scalenew) = -1.524" [1] "Empirical dV/V vs dP/P relationship =-1.228, r^2 =0.024, damping =0.639" [1] "Tue May 13 01:43:30 PM MDT 2025 done redoing corrections" [1] "JOGAVES MATCHES" [1] "MASKING 20 VHR ROWS TO MAKE JOGS EVEN FOR SHAPE CALCS" [1] "Plotting Jog Shapes" [1] "Tue May 13 01:43:47 PM MDT 2025 done recreating mr" [1] "Using SPLINE interpolation between jogs" [1] "RECALC TIME MATCHES" [1] "Tue May 13 01:43:48 PM MDT 2025 done recreating mr o2diffs" [1] "Screening 0 MR O2 values when lamp off" [1] "Screening 0 HR O2 values when lamp off" [1] "MAX O2 DIFFERENCE (ORIG-UNCORR) -43.527 PER MEG" [1] "MAX O2 DIFFERENCE (ORIG-PDCORR) -42.214 PER MEG" [1] "MAX O2 DIFFERENCE (ORIG-PDPRCORR) -43.189 PER MEG" [1] "MAX O2 DIFFERENCE (PDCORR-PDPRCORR) 13.866 PER MEG" [1] "Removing O2 and CO2 glitches" [1] 0 [1] 0 [1] 33 [1] 33 [1] 33 [1] "O2 glitches (based on SP-SP, WT-WT, or single jog spikes or steps) found at (row, time, o2diff, o2jog, posneg):" [,1] [,2] [,3] [,4] [,5] [1,] 177 275.5992 -598.6060 942893.7 -1 [2,] 329 275.6039 -564.8704 948094.0 -1 [3,] 395 275.6059 -904.1940 949986.5 -1 [4,] 463 275.6081 -593.0976 952533.7 -1 [5,] 529 275.6103 -865.5532 954205.9 -1 [6,] 593 275.6124 -633.8131 956428.8 -1 [7,] 723 275.6165 -612.0285 959769.3 -1 [8,] 787 275.6186 -752.0126 960386.1 -1 [9,] 851 275.6208 -609.6233 960781.9 -1 [10,] 915 275.6229 -921.1738 960569.1 -1 [11,] 977 275.6248 -691.4379 961003.9 -1 [12,] 983 275.6250 -98.6932 961621.8 -1 [13,] 985 275.6250 375.8897 962096.4 -1 [14,] 987 275.6251 568.8950 962323.3 -1 [15,] 1473 275.6398 -775.3556 966029.2 -1 [16,] 1693 275.6468 -779.2446 968913.7 -1 [17,] 1951 275.6548 -770.0807 972078.0 -1 [18,] 2235 275.6623 -782.3700 974880.9 -1 [19,] 2523 275.6700 -765.4259 977687.2 -1 [20,] 2711 275.6749 -744.1765 979641.2 -1 [21,] 2713 275.6750 -532.0866 979856.5 -1 [22,] 3349 275.6932 -797.2317 995156.5 -1 [23,] 3857 275.7092 -789.3711 1001038.8 -1 [24,] 4359 275.7250 -518.5791 1000239.6 -1 [25,] 4425 275.7271 -731.7735 1000210.6 -1 [26,] 4955 275.7417 -762.2785 1003282.7 -1 [27,] 5847 275.7650 -754.8653 1007660.7 -1 [28,] 6379 275.7807 -877.7358 1016408.5 -1 [29,] 6447 275.7827 -618.5035 1017248.2 -1 [30,] 6661 275.7891 -543.9452 1018574.5 -1 [31,] 6737 275.7915 -750.0805 1019027.2 -1 [32,] 7019 275.8005 -546.2413 1020965.4 -1 [33,] 7277 275.8089 -651.9775 1021552.9 -1 [1] "CO2 glitches found at:" [,1] [,2] [,3] [1,] 1704 275.6472 289.122 [2,] 1819 275.6508 960879.100 [3,] 6388 275.7810 9737.195 [1] "MAX O2 DIFFERENCE (ORIG-UNCORR) -43.527 PER MEG" [1] "MAX O2 DIFFERENCE (ORIG-PDCORR) -42.214 PER MEG" [1] "MAX O2 DIFFERENCE (ORIG-PDPRCORR) -43.189 PER MEG" [1] "MAX O2 DIFFERENCE (PDCORR-PDPRCORR) 13.866 PER MEG" [1] 275.6350 275.6367 [1] "Screened 54 O2 values from screeno2 period 275.635011574074 to 275.636747685185" [1] 275.6433 275.6450 [1] "Screened 55 O2 values from screeno2 period 275.643298611111 to 275.645034722222" [1] 275.6513 275.6530 [1] "Screened 55 O2 values from screeno2 period 275.651296296296 to 275.653032407407" [1] 275.6588 275.6605 [1] "Screened 65 O2 values from screeno2 period 275.658784722222 to 275.660520833333" [1] 275.6657 275.6674 [1] "Screened 65 O2 values from screeno2 period 275.665671296296 to 275.667407407407" [1] 275.6878 275.6895 [1] "Screened 66 O2 values from screeno2 period 275.6878125 to 275.689548611111" [1] 275.7051 275.7068 [1] "Screened 55 O2 values from screeno2 period 275.705104166667 to 275.706840277778" [1] 275.7177 275.7195 [1] "Screened 55 O2 values from screeno2 period 275.717719907407 to 275.719456018519" [1] 275.7301 275.7318 [1] "Screened 55 O2 values from screeno2 period 275.730081018518 to 275.73181712963" [1] 275.7605 275.7622 [1] "Screened 67 O2 values from screeno2 period 275.760462962963 to 275.762199074074" [1] 391.3386 404.0181 392.4445 394.3873 369.4096 [1] -452.7198 -717.5633 -584.4989 -580.7611 -827.3404 [1] "Nominal linear calibrations:" [1] "co2scl = 1.068 co2off = -11.73 o2scla = -0.109 o2offa = -184" [1] "(dV/V)/(dXO2/XO2) = -1.921" [1] "(dV/V)/d(O2/N2) = -1.519" [1] "Finding cal and sample periods" [1] "Median Cal: 2.98 minutes" [1] "Average Cal: 3.02 minutes" [1] "Shortest Cal: 1.58 minutes" [1] "Longest Cal: 5 minutes" [1] "Removed 6 calibration period(s) with less than 30 seconds of data" [1] "at: 20" [1] "at: 22" [1] "at: 30" [1] "at: 39" [1] "at: 42" [1] "at: 69" [1] "Reading in aircraft data" [1] "Plotting CO2 cabin P sensitivity" [1] "HS+LS CO2 vs. PACAB: m = -0.000428" [1] "HS CO2 vs. PACAB: m = -0.00161, r^2 = 0.03" [1] "LS CO2 vs. PACAB: m = 0.000758, r^2 = 0.02" [1] "HS+LS LIH2O vs. PACAB: m = -0.000156" [1] "HS LIH2O vs. PACAB: m = -0.000137, r^2 = 0.89" [1] "LS LIH2O vs. PACAB: m = -0.000176, r^2 = 0.91" [1] "HS+LS CO2 vs. LIH2O: m = 13.3" [1] "HS CO2 vs. LIH2O: m = 20.6, r^2 = 0.12" [1] "LS CO2 vs. LIH2O: m = 5.88, r^2 = 0.02" [1] "APCO2 vs. PACAB: m = -0.0267, r^2 = 0.17" [1] "Plotting O2 cabin P, H2O, PaO2, and TPBox sensitivity" [1] "HS+LS O2 vs. PACAB: m = -0.00128" [1] "HS O2 vs. PACAB: m = 0.0164, r^2 = 0" [1] "LS O2 vs. PACAB: m = -0.019, r^2 = 0.01" [1] "HS+LS O2 vs. LIH2O: m = -71.3" [1] "HS O2 vs. LIH2O: m = -262, r^2 = 0.01" [1] "LS O2 vs. LIH2O: m = 120, r^2 = 0.01" [1] "HS+LS O2 vs. PAO2: m = -3.6" [1] "HS O2 vs. PAO2: m = -7.33, r^2 = 0.05" [1] "LS O2 vs. PAO2: m = 0.123, r^2 = 0" [1] "HS+LS O2 vs. TPBOX: m = 0.762" [1] "HS O2 vs. TPBOX: m = 1.87, r^2 = 0.03" [1] "LS O2 vs. TPBOX: m = -0.343, r^2 = 0" [1] "APO vs. PACAB: m = 0.228, r^2 = 0.17" [1] "APO vs. PAO2: m = 8.26, r^2 = 0.02" [1] "APO vs. TPBOX: m = -3.69, r^2 = 0.02" [1] "Plotting cal meas P and T sensitivity" [1] "HS CO2 vs. LiCor Cell P: m = -3.5, r^2 = 0.12" [1] "LS CO2 vs. LiCor Cell P: m = -9.09, r^2 = 0.77" [1] "HS CO2 vs. CO2 Vol P: m = 0.584, r^2 = 0.08" [1] "LS CO2 vs. CO2 Vol P: m = 1.7, r^2 = 0.56" [1] "HS CO2 vs. SLOPE: m = 0.718, r^2 = 0.25" [1] "LS CO2 vs. SLOPE: m = -0.988, r^2 = 0.06" [1] "HS O2 vs. O2 Vol P: m = -21.2, r^2 = 0.11" [1] "LS O2 vs. O2 Vol P: m = -0.618, r^2 = 0" [1] "HS O2 vs. O2 Cell T: m = 8.15, r^2 = 0.08" [1] "LS O2 vs. O2 Cell T: m = -0.184, r^2 = 0" [1] "HS O2 vs. CO2 Vol P: m = 33.1, r^2 = 0.04" [1] "LS O2 vs. CO2 Vol P: m = 28.3, r^2 = 0.07" [1] "HS O2 vs. O2 Vol P / O2 Cell T: m = -17.9, r^2 = 0.06" [1] "LS O2 vs. O2 Vol P / O2 Cell T: m = -4.49, r^2 = 0.01" [1] "HS O2 vs. CO2 Vol P / O2 Cell T: m = 33.1, r^2 = 0.04" [1] "LS O2 vs. CO2 Vol P / O2 Cell T: m = 28.3, r^2 = 0.07" [1] "HS O2 vs. LIH2O: m = -860, r^2 = 0.03" [1] "LS O2 vs. LIH2O: m = 244, r^2 = 0.01" [1] "HS O2 vs. SLOPE: m = 0.705, r^2 = 0.07" [1] "LS O2 vs. SLOPE: m = -0.305, r^2 = 0.03" [1] "HS CO2 vs. Time: m = 18.577, r^2 = 0.93" [1] "LS CO2 vs. Time: m = 17.846, r^2 = 0.94" [1] "HS CO2 vs. Cylinder T: m = -4.166, r^2 = 0.53" [1] "LS CO2 vs. Cylinder T: m = -3.904, r^2 = 0.331" [1] "HS CO2 vs. Regulator T: m = -0.3, r^2 = 0.2" [1] "LS CO2 vs. Regulator T: m = -0.573, r^2 = 0.815" [1] "HS CO2 vs. Ref. Volumes T: m = -0.206, r^2 = 0.08" [1] "LS CO2 vs. Ref. Volumes T: m = -0.626, r^2 = 0.771" [1] "HS CO2 vs. Manifold T: m = -0.088, r^2 = 0.07" [1] "LS CO2 vs. Manifold T: m = -0.269, r^2 = 0.642" [1] "HS CO2 vs. Pump Box T: m = -0.128, r^2 = 0.19" [1] "LS CO2 vs. Pump Box T: m = -0.262, r^2 = 0.793" [1] "HS O2 vs. Cylinder T: m = 100.453, r^2 = 0.05" [1] "LS O2 vs. Cylinder T: m = 69.735, r^2 = 0.047" [1] "HS O2 vs. Regulator T: m = 20.24, r^2 = 0.14" [1] "LS O2 vs. Regulator T: m = 0.457, r^2 = 0" [1] "HS O2 vs. Ref. Volumes T: m = 28.168, r^2 = 0.23" [1] "LS O2 vs. Ref. Volumes T: m = 2.135, r^2 = 0.004" [1] "HS O2 vs. Manifold T: m = 9.184, r^2 = 0.11" [1] "LS O2 vs. Manifold T: m = 0.049, r^2 = 0" [1] "HS O2 vs. Pump Box T: m = 4.51, r^2 = 0.04" [1] "LS O2 vs. Pump Box T: m = -1.015, r^2 = 0.005" [1] "HS CO2 vs. WT P: m = -0.005, r^2 = 0.18" [1] "LS CO2 vs. WT P: m = -0.004, r^2 = 0.178" [1] "HS CO2 vs. SP P: m = -0.00512, r^2 = 0.19" [1] "HS CO2 vs. PACAB comb with SP P: m = -0.001, r^2 = 0.18" [1] "HS CO2 vs. SP P comb with PACAB: m = -0.004, r^2 = 0.18" [1] "LS CO2 vs. SP P: m = -0.00173, r^2 = 0.17" [1] "LS CO2 vs. PACAB comb with SP P: m = -0.21, r^2 = 0.71" [1] "LS CO2 vs. SP P comb with PACAB: m = 0.199, r^2 = 0.71" [1] "HS CO2 vs. WT FL: m = -0.078, r^2 = 0.12" [1] "LS CO2 vs. WT FL: m = -0.194, r^2 = 0.788" [1] "HS CO2 vs. SP FL: m = 0.005, r^2 = 0" [1] "LS CO2 vs. SP FL: m = -0.068, r^2 = 0.092" [1] "HS O2 vs. SLPDIFF: m = 0.872, r^2 = 0.2" [1] "LS O2 vs. SLPDIFF: m = -1.888, r^2 = 0.418" [1] "HS O2 vs. WT P: m = 0.079, r^2 = 0.01" [1] "LS O2 vs. WT P: m = -0.002, r^2 = 0" [1] "HS O2 vs. SP P: m = 0.09, r^2 = 0.01" [1] "LS O2 vs. SP P: m = -0.004, r^2 = 0" [1] "HS O2 vs. WT FL: m = 7.934, r^2 = 0.2" [1] "LS O2 vs. WT FL: m = 0.506, r^2 = 0.002" [1] "HS O2 vs. SP FL: m = 5.942, r^2 = 0.06" [1] "LS O2 vs. SP FL: m = -0.219, r^2 = 0" [1] "Plotting calibration shapes" [1] "Empirical LP O2 lag time(s): 65.1, 83.52, 64.57, 42.47, 66.68, 79.81, 69.06, 74.16, 71.08, 77.28, 66.92, 78.33, 69.88, -4, 66.71, 13.6, 73.01, 32.58, 66.19, 74.74" [1] "Empirical LP O2 lag flow(s): 96, 96.2, 95.9, 96.1, 96.1, 96.1, 96.2, 96.1, 96.2, 96.4, 97.8, 98.9, 100.4, 100.5, 100.5, 100.4, 100, 100.5, 100.2, 100.1" [1] "Empirical LP CO2 lag time(s): 59.65, -22.84, 51.53, 42.47, -15.15, -16.78, 62.15, 0.29, -4.6, 27, 71.41, -20.25, 72.48, 4.3, 50.47, 49.28, 73.01, 41.56, 23.16, -1.64" [1] "Empirical LP CO2 lag flow(s): 96, 96.2, 95.9, 96.1, 96.1, 96.1, 96.2, 96.1, 96.2, 96.4, 97.8, 98.9, 100.4, 100.5, 100.5, 100.4, 100, 100.5, 100.2, 100.1" [1] "Plotting calibration sds and slopes (and running sds)" [1] "Median all cal CO2 SD = 0.88" [1] "Median all HS CO2 SD = 0.89" [1] "Median all LS CO2 SD = 0.85" [1] "Median all LT CO2 SD = 0.86" [1] "Median all WT CO2 SD = 0.78" [1] "Median all HS or LS CO2 SD = 0.87" [1] "Median all cal O2 SD = 7.91" [1] "Median all HS O2 SD = 6.07" [1] "Median all LS O2 SD = 7.17" [1] "Median all LT O2 SD = 9.14" [1] "Median all WT O2 SD = 9.37" [1] "Median all HS or LS O2 SD = 6.76" [1] "Median all cal CO2 Slope = -0.03" [1] "Median all HS CO2 Slope = 0.03" [1] "Median all LS CO2 Slope = -0.02" [1] "Median all LT CO2 Slope = -0.09" [1] "Median all WT CO2 Slope = -0.11" [1] "Median all HS or LS CO2 Slope = 0.01" [1] "Median all cal O2 Slope = 1.23" [1] "Median all HS O2 Slope = 1.04" [1] "Median all LS O2 Slope = 0.31" [1] "Median all LT O2 Slope = 0.61" [1] "Median all WT O2 Slope = -2.78" [1] "Median all HS or LS O2 Slope = 0.69" [1] "Plotting calibration fits" [1] "Plotting calibration parameters" [1] "SD of O2 Calibration = 41.56 per meg" [1] "Initial O2 Calibration Offset = 5.13 per meg" [1] "SD of O2 Calibration (post-takeoff) = 48.73 per meg" [1] "Initial O2 Calibration Offset (from post-takeoff median) = 41.96 per meg" [1] "Last Preflight O2 Calibration Offset (from post-takeoff median) = 39.07 per meg" [1] "Plotting all variables" [1] "Plotting calibration data" [1] "using HIPPO4-5 inlet volumes" [1] "before any lag: 51328.17 70326.4" [1] "Time lag adjustments = 0 + 0 * P/Flow seconds" [1] "Median CO2 delay + tlagff = 36.4 at 99.2 sccm" [1] "Median O2 delay + tlagff = 37.1 at 99.2 sccm" [1] "For comparison to empirical LP lag times:" [1] "Median CO2 delay from 3-way + tlagff = 33.9 at 99.2 sccm" [1] "Median O2 delay from 3-way + tlagff = 34.6 at 99.2 sccm" [1] "after calc and ff lag: 51290.31 70290.18" [1] "2025-05-13 13:44:13 MDT" [1] "Calculating 1-sec interpolated values" [1] "Calculating 10-sec medians" [1] "Plotting aircraft data" [1] "Outputting data files" [1] "Plotting intervals" [1] "Excluding 563 interpolated O2 values from O2 vs. CO2 plot for tdiff>5" [1] "making html" [1] "making index html"