1 The SAS System 14:12 Monday, September 26, 2022 NOTE: Unable to open SASUSER.REGSTRY. WORK.REGSTRY will be opened instead. NOTE: All registry changes will be lost at the end of the session. WARNING: Unable to copy SASUSER registry to WORK registry. Because of this, you will not see registry customizations during this session. NOTE: Copyright (c) 2016 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software 9.4 (TS1M7) Licensed to US ENVIRONMENTAL PROTECTION AGENCY, Site 70072220. NOTE: This session is executing on the Linux 3.10.0-1160.76.1.el7.x86_64 (LIN X64) platform. NOTE: Analytical products: SAS/STAT 15.2 NOTE: Additional host information: Linux LIN X64 3.10.0-1160.76.1.el7.x86_64 #1 SMP Tue Jul 26 14:15:37 UTC 2022 x86_64 Red Hat Enterprise Linux Server release 7.9 (Maipo) NOTE: SAS initialization used: real time 0.07 seconds cpu time 0.03 seconds 1 %let base=.; 2 3 libname sasd "&base./sasd"; NOTE: Libref SASD was successfully assigned as follows: Engine: V9 Physical Name: /work/ROMO/users/hou/400.1/EPA-R6-2016/2.postproc/xTESTING/1.master/yQA.zeroRGN/4/sasd 4 5 options nosymbolgen mprint mlogic; 6 7 ********************************************************************; 8 9 %MACRO program; 10 11 /****** 12 camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.origRGN.txt 13 ------------------------------------------------ 14 Tracer,site,ss_rayleigh,mon,day,year,fsrh,flrh,fssrh,frh,date,NO2,Soil,Amm_NO3,OMC,LAC,CM,Amm_SO4,Sea_Salt,L_oc,S_oc,L_so 14 ! 4,S_so4,L_no3,S_no3,E_oc,E_so4,E_no3,E_ec,E_soil,E_cm,E_no2,E_salt,T_bext,dv,nc_ext,elta_dv,E_oc_p,E_so4_p,E_no3_p,E_ec_p 14 ! ,E_soil_p,E_cm_p,E_no2_p,T_bext_p,dv_p 15 S001,BAND1,9,1,1,2002,2.66,2.1,2.91,2.2,01/01/02,3.61064,1.06337,2.0022,0.521507,0.180249,1.78806,0.810426,0,0.0135984775 15 ! 5245,0.50790852244755,0.0328395150738,0.7775864849262,0.200440242,1.801759758,1.50509457592308,4.88145842173203,13.649149 15 ! 286892,1.80249,1.06337,1.072836,1.1915112,0,25.1659094845471,9.22905187433472,14.00281,10.2862043215348,97.9658240067896, 15 ! 84.4270021606255,95.0189608455167,97.534170968475,99.8000938526513,99.6533428450409,99.9357867233515,64.8390190988065,68. 15 ! 0524330572417 16 S001,BAND1,9,1,2,2002,2.66,2.1,2.91,2.2,01/02/02,3.54629,1.20898,2.97722,0.504017,0.186733,2.17892,1.24482,0,0.0127016568 16 ! 1445,0.49131534318555,0.07747884162,1.16734115838,0.44319194642,2.53402805358,1.45316306748768,7.61226718236936,20.923820 16 ! 8402129,1.86733,1.20898,1.307352,1.1702757,0,35.54318879007,12.6816344968246,14.00281,12.6364348279769,97.6840539780506,8 16 ! 3.3087113489762,95.70045658671,97.0757649798813,99.6809168487447,99.4128087088635,99.9199803895603,71.1169852697451,78.81 16 ! 65065198656 17 S001,BAND1,9,1,3,2002,2.66,2.1,2.91,2.2,01/03/02,1.72125,0.394685,0.601157,0.231319,0.0889029,0.958739,0.141538,0,0.00267 17 ! 542398805,0.22864357601195,0.0010016502722,0.1405363497278,0.01806948693245,0.58308751306755,0.656522099160565,0.83251535 17 ! 3350862,3.91595488846978,0.889029,0.394685,0.5752434,0.5680125,0,7.83196224098121,-2.44372008898346,14.00281,4.4424573220 17 ! 8035,94.6295837592383,53.8650013553335,82.049760920751,93.5638781541259,98.7821329442299,98.2648810865991,99.616292798111 2 The SAS System 14:12 Monday, September 26, 2022 17 ! ,39.4220753460531,-35.5982366539213 18 F012,WIMO1,11,12,29,2002,3.2,2.4,2.63,2.8,12/29/02,6.46951e-06,0.000279461,0.0001433461,1.475693e-05,9.94636e-06,7.86728e 18 ! -16,0.002008757,0,1.0888349151245e-11,1.47569191116508e-05,2.0175523425245e-07,0.00200855524476575,1.0274052192605e-09,0. 18 ! 000143345072594781,4.13194399315522e-05,0.0141425531434495,0.0011009027329678,9.94636e-05,0.000279461,4.720368e-16,2.1349 18 ! 383e-06,0,0.0156658348546493,-64.5885815472788,16.62253,0.00942002053164881,0.00215516961928577,0.0733644249947294,0.0062 18 ! 6028548126112,0.00561297494963404,0.0118339960449035,1.09712264828479e-14,0.000366023954602803,0.0257524548719343,-357.72 18 ! 498112837 19 F012,WIMO1,11,12,30,2002,3.2,2.4,2.63,2.8,12/30/02,2.331593e-06,0.0001278892,6.31389e-05,6.94981e-06,4.77547e-06,7.2161e- 19 ! 16,0.00102921,0,2.414992951805e-12,6.94980758500705e-06,5.2963661205e-08,0.00102915703633879,1.993260346605e-10,6.3138700 19 ! 6739653e-05,1.94594759694767e-05,0.0072458756772022,0.000484907660926718,4.77547e-05,0.0001278892,4.32966e-16,7.6942569e- 19 ! 07,0,0.00792665613978883,-71.4010909741591,16.62253,0.00476748532208049,0.000780565645773459,0.0481776873972115,0.0021368 19 ! 1011984107,0.00189122281758525,0.0041277616218083,7.47790907461121e-15,8.37040480771992e-05,0.0121124291081846,-380.07926 19 ! 0406046 20 F012,WIMO1,11,12,31,2002,3.2,2.4,2.63,2.8,12/31/02,1.259447e-08,9.55557e-08,4.04931e-07,7.86377e-09,5.82946e-09,6.8661e-1 20 ! 6,1.245352e-06,0,3.091943930645e-18,7.86376999690806e-09,7.75450801952e-14,1.24535192245492e-06,8.19845573805e-15,4.04930 20 ! 991801544e-07,2.20185560102034e-08,8.76727842740196e-06,3.10987011738496e-06,5.82946e-08,9.55557e-08,4.11966e-16,4.156175 20 ! 1e-09,0,1.20571735763091e-05,-136.284358506168,16.62253,7.25351026886756e-06,1.43765890186667e-06,0.000167925611678289,1. 20 ! 79751279705781e-05,3.5073281671159e-06,3.7946334257281e-06,6.72382388656038e-15,6.81740283641875e-07,2.46894135312848e-05 20 ! ,-859.366335088252 21 ******/ 22 23 %let desc=camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv; 24 25 data sasd.v1a_bextdv_orig; 26 length n 8.; 27 infile "&base./&desc..origRGN.txt" DLM=',' DSD MISSOVER FIRSTOBS=2; 28 input Tracer :$9. site :$5. ssray mon day year fsrh flrh fssrh frh date :$8. 29 oNO2 oSoil oAmm_NO3 oOMC oLAC oCM oAmm_SO4 oSea_Salt oL_oc oS_oc oL_so4 oS_so4 oL_no3 oS_no3 30 oE_oc oE_so4 oE_no3 oE_ec oE_soil oE_cm oE_no2 oE_salt oT_bext odv onc_ext odelta_dv 31 oE_oc_p oE_so4_p oE_no3_p oE_ec_p oE_soil_p oE_cm_p oE_no2_p oT_bext_p odv_p; 32 n = _N_; 33 run; 34 35 data sasd.v1b_bextdv_zero; 36 length n 8.; 37 infile "&base./&desc..zeroRGN.txt" DLM=',' DSD MISSOVER FIRSTOBS=2; 38 input Tracer :$9. site :$5. ssray mon day year fsrh flrh fssrh frh date :$8. 39 zNO2 zSoil zAmm_NO3 zOMC zLAC zCM zAmm_SO4 zSea_Salt zL_oc zS_oc zL_so4 zS_so4 zL_no3 zS_no3 40 zE_oc zE_so4 zE_no3 zE_ec zE_soil zE_cm zE_no2 zE_salt zT_bext zdv znc_ext zdelta_dv 41 zE_oc_p zE_so4_p zE_no3_p zE_ec_p zE_soil_p zE_cm_p zE_no2_p zT_bext_p zdv_p; 42 n = _N_; 43 run; 44 45 ******************; 46 47 data sasd.v2a_bextdv_mrg1; 48 length n 8. Tracer $9. site $5. ssray mon day year fsrh flrh fssrh frh 8. date $8. 49 oNO2 zNO2 50 oSoil zSoil 51 oAmm_NO3 zAmm_NO3 52 oOMC zOMC 53 oLAC zLAC 54 oCM zCM 55 oAmm_SO4 zAmm_SO4 56 oSea_Salt zSea_Salt 57 oL_oc zL_oc 58 oS_oc zS_oc 59 oL_so4 zL_so4 3 The SAS System 14:12 Monday, September 26, 2022 60 oS_so4 zS_so4 61 oL_no3 zL_no3 62 oS_no3 zS_no3 63 oE_oc zE_oc 64 oE_so4 zE_so4 65 oE_no3 zE_no3 66 oE_ec zE_ec 67 oE_soil zE_soil 68 oE_cm zE_cm 69 oE_no2 zE_no2 70 oE_salt zE_salt 71 oT_bext zT_bext 72 odv zdv 73 onc_ext znc_ext 74 odelta_dv zdelta_dv 75 oE_oc_p zE_oc_p 76 oE_so4_p zE_so4_p 77 oE_no3_p zE_no3_p 78 oE_ec_p zE_ec_p 79 oE_soil_p zE_soil_p 80 oE_cm_p zE_cm_p 81 oE_no2_p zE_no2_p 82 oT_bext_p zT_bext_p 83 odv_p zdv_p 8.; 84 merge sasd.v1a_bextdv_orig(in=a) 85 sasd.v1b_bextdv_zero(in=b); 86 by n Tracer site mon day year date; 87 run; 88 89 data sasd.v2b_bextdv_mrg1(keep=n Tracer site ssray mon day year date 90 oNO2 zNO2 dfNO2 pctNO2 91 oE_so4 zE_so4 dfE_so4 pctE_so4 92 oE_no3 zE_no3 dfE_no3 pctE_no3 93 oE_no2 zE_no2 dfE_no2 pctE_no2 94 oT_bext zT_bext dfT_bext pctT_bext 95 odv zdv dfdv pctdv 96 odelta_dv zdelta_dv dfdelta_dv pctdelta_dv 97 oE_so4_p zE_so4_p dfE_so4_p pctE_so4_p 98 oE_no3_p zE_no3_p dfE_no3_p pctE_no3_p 99 oE_no2_p zE_no2_p dfE_no2_p pctE_no2_p 100 oT_bext_p zT_bext_p dfT_bext_p pctT_bext_p 101 odv_p zdv_p dfdv_p pctdv_p ); 102 length n 8. Tracer $9. site $5. ssray mon day year 8. date $8. 103 oNO2 zNO2 dfNO2 pctNO2 104 oE_so4 zE_so4 dfE_so4 pctE_so4 105 oE_no3 zE_no3 dfE_no3 pctE_no3 106 oE_no2 zE_no2 dfE_no2 pctE_no2 107 oT_bext zT_bext dfT_bext pctT_bext 108 odv zdv dfdv pctdv 109 odelta_dv zdelta_dv dfdelta_dv pctdelta_dv 110 oE_so4_p zE_so4_p dfE_so4_p pctE_so4_p 111 oE_no3_p zE_no3_p dfE_no3_p pctE_no3_p 112 oE_no2_p zE_no2_p dfE_no2_p pctE_no2_p 113 oT_bext_p zT_bext_p dfT_bext_p pctT_bext_p 114 odv_p zdv_p dfdv_p pctdv_p 8.; 115 set sasd.v2a_bextdv_mrg1; 116 dfNO2 = zNO2 - oNO2 ; 117 dfE_so4 = zE_so4 - oE_so4 ; 4 The SAS System 14:12 Monday, September 26, 2022 118 dfE_no3 = zE_no3 - oE_no3 ; 119 dfE_no2 = zE_no2 - oE_no2 ; 120 dfT_bext = zT_bext - oT_bext ; 121 dfdv = zdv - odv ; 122 dfdelta_dv = zdelta_dv - odelta_dv ; 123 dfE_so4_p = zE_so4_p - oE_so4_p ; 124 dfE_no3_p = zE_no3_p - oE_no3_p ; 125 dfE_no2_p = zE_no2_p - oE_no2_p ; 126 dfT_bext_p = zT_bext_p - oT_bext_p ; 127 dfdv_p = zdv_p - odv_p ; 128 pctNO2 = 100. * (dfNO2 / oNO2 ); 129 pctE_so4 = 100. * (dfE_so4 / oE_so4 ); 130 pctE_no3 = 100. * (dfE_no3 / oE_no3 ); 131 pctE_no2 = 100. * (dfE_no2 / oE_no2 ); 132 pctT_bext = 100. * (dfT_bext / oT_bext ); 133 pctdv = 100. * (dfdv / odv ); 134 pctdelta_dv = 100. * (dfdelta_dv / odelta_dv ); 135 pctE_so4_p = 100. * (dfE_so4_p / oE_so4_p ); 136 pctE_no3_p = 100. * (dfE_no3_p / oE_no3_p ); 137 pctE_no2_p = 100. * (dfE_no2_p / oE_no2_p ); 138 pctT_bext_p = 100. * (dfT_bext_p / oT_bext_p ); 139 pctdv_p = 100. * (dfdv_p / odv_p ); 140 format 141 pctNO2 142 pctE_so4 143 pctE_no3 144 pctE_no2 145 pctT_bext 146 pctdv 147 pctdelta_dv 148 pctE_so4_p 149 pctE_no3_p 150 pctE_no2_p 151 pctT_bext_p 152 pctdv_p 8.2; 153 run; 154 155 ******************; 156 157 *proc means data=sasd.v2b_bextdv_mrg1 min p5 p10 p25 p50 p75 p90 p95 max mean; 158 proc means data=sasd.v2b_bextdv_mrg1 min p1 p5 p10 p90 p95 p99 max; 159 var pctT_bext pctT_bext_p pctdv pctdv_p pctdelta_dv; 160 run; 161 /* 162 proc univariate data=sasd.v2b_bextdv_mrg1; 163 var pctT_bext; 164 histogram / vscale=count midpoints=-100 to 0 by 2; 165 var pctdv; 166 histogram / vscale=count midpoints=-7500 to 25000 by 500; 167 var pctdelta_dv; 168 histogram / vscale=count midpoints=-100 to 0 by 2; 169 run; 170 */ 171 172 proc univariate data=sasd.v2b_bextdv_mrg1; 173 var pctT_bext pctdv pctdelta_dv; 174 histogram / midpoints=-100 to 10 by 2; 175 run; 5 The SAS System 14:12 Monday, September 26, 2022 176 177 178 ******************; 179 180 proc export data=sasd.v2b_bextdv_mrg1 181 outfile="&base./v2b_bextdv_mrg1.csv" 182 DBMS=CSV REPLACE; 183 run; 184 185 ******************; 186 %MEND program; 187 %program MLOGIC(PROGRAM): Beginning execution. MLOGIC(PROGRAM): %LET (variable name is DESC) MPRINT(PROGRAM): data sasd.v1a_bextdv_orig; MPRINT(PROGRAM): length n 8.; MPRINT(PROGRAM): infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.origRGN.txt" DLM=',' DSD MISSOVER FIRSTOBS=2; MPRINT(PROGRAM): input Tracer :$9. site :$5. ssray mon day year fsrh flrh fssrh frh date :$8. oNO2 oSoil oAmm_NO3 oOMC oLAC oCM oAmm_SO4 oSea_Salt oL_oc oS_oc oL_so4 oS_so4 oL_no3 oS_no3 oE_oc oE_so4 oE_no3 oE_ec oE_soil oE_cm oE_no2 oE_salt oT_bext odv onc_ext odelta_dv oE_oc_p oE_so4_p oE_no3_p oE_ec_p oE_soil_p oE_cm_p oE_no2_p oT_bext_p odv_p; MPRINT(PROGRAM): n = _N_; MPRINT(PROGRAM): run; NOTE: The infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.origRGN.txt" is: Filename=/work/ROMO/users/hou/400.1/EPA-R6-2016/2.postproc/xTESTING/1.master/yQA.zeroRGN/4/camx_v630.EPA-R6-2016.TX-BART.CAMxR un4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.origRGN.txt, Owner Name=hou,Group Name=romo, Access Permission=-rw-r--r--, Last Modified=15Sep2022:11:05:06, File Size (bytes)=124201356 NOTE: 204035 records were read from the infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.origRGN.txt". The minimum record length was 363. The maximum record length was 688. NOTE: The data set SASD.V1A_BEXTDV_ORIG has 204035 observations and 47 variables. NOTE: DATA statement used (Total process time): real time 1.28 seconds cpu time 1.14 seconds MPRINT(PROGRAM): data sasd.v1b_bextdv_zero; MPRINT(PROGRAM): length n 8.; MPRINT(PROGRAM): infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.zeroRGN.txt" DLM=',' DSD MISSOVER FIRSTOBS=2; MPRINT(PROGRAM): input Tracer :$9. site :$5. ssray mon day year fsrh flrh fssrh frh date :$8. zNO2 zSoil zAmm_NO3 zOMC zLAC zCM zAmm_SO4 zSea_Salt zL_oc zS_oc zL_so4 zS_so4 zL_no3 zS_no3 zE_oc zE_so4 zE_no3 zE_ec zE_soil zE_cm zE_no2 zE_salt zT_bext zdv znc_ext zdelta_dv zE_oc_p zE_so4_p zE_no3_p zE_ec_p zE_soil_p zE_cm_p zE_no2_p zT_bext_p zdv_p; MPRINT(PROGRAM): n = _N_; MPRINT(PROGRAM): run; NOTE: The infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.zeroRGN.txt" is: Filename=/work/ROMO/users/hou/400.1/EPA-R6-2016/2.postproc/xTESTING/1.master/yQA.zeroRGN/4/camx_v630.EPA-R6-2016.TX-BART.CAMxR 6 The SAS System 14:12 Monday, September 26, 2022 un4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.zeroRGN.txt, Owner Name=hou,Group Name=romo, Access Permission=-rw-r--r--, Last Modified=22Sep2022:13:57:04, File Size (bytes)=116286949 NOTE: 204035 records were read from the infile "./camx_v630.EPA-R6-2016.TX-BART.CAMxRun4.36-12.PiG-PSAT.2002001-2002365.data.template.bext_dv.zeroRGN.txt". The minimum record length was 363. The maximum record length was 646. NOTE: The data set SASD.V1B_BEXTDV_ZERO has 204035 observations and 47 variables. NOTE: DATA statement used (Total process time): real time 1.18 seconds cpu time 1.08 seconds MPRINT(PROGRAM): ******************; MPRINT(PROGRAM): data sasd.v2a_bextdv_mrg1; MPRINT(PROGRAM): length n 8. Tracer $9. site $5. ssray mon day year fsrh flrh fssrh frh 8. date $8. oNO2 zNO2 oSoil zSoil oAmm_NO3 zAmm_NO3 oOMC zOMC oLAC zLAC oCM zCM oAmm_SO4 zAmm_SO4 oSea_Salt zSea_Salt oL_oc zL_oc oS_oc zS_oc oL_so4 zL_so4 oS_so4 zS_so4 oL_no3 zL_no3 oS_no3 zS_no3 oE_oc zE_oc oE_so4 zE_so4 oE_no3 zE_no3 oE_ec zE_ec oE_soil zE_soil oE_cm zE_cm oE_no2 zE_no2 oE_salt zE_salt oT_bext zT_bext odv zdv onc_ext znc_ext odelta_dv zdelta_dv oE_oc_p zE_oc_p oE_so4_p zE_so4_p oE_no3_p zE_no3_p oE_ec_p zE_ec_p oE_soil_p zE_soil_p oE_cm_p zE_cm_p oE_no2_p zE_no2_p oT_bext_p zT_bext_p odv_p zdv_p 8.; MPRINT(PROGRAM): merge sasd.v1a_bextdv_orig(in=a) sasd.v1b_bextdv_zero(in=b); MPRINT(PROGRAM): by n Tracer site mon day year date; MPRINT(PROGRAM): run; NOTE: There were 204035 observations read from the data set SASD.V1A_BEXTDV_ORIG. NOTE: There were 204035 observations read from the data set SASD.V1B_BEXTDV_ZERO. NOTE: The data set SASD.V2A_BEXTDV_MRG1 has 204035 observations and 82 variables. NOTE: DATA statement used (Total process time): real time 0.48 seconds cpu time 0.30 seconds MPRINT(PROGRAM): data sasd.v2b_bextdv_mrg1(keep=n Tracer site ssray mon day year date oNO2 zNO2 dfNO2 pctNO2 oE_so4 zE_so4 dfE_so4 pctE_so4 oE_no3 zE_no3 dfE_no3 pctE_no3 oE_no2 zE_no2 dfE_no2 pctE_no2 oT_bext zT_bext dfT_bext pctT_bext odv zdv dfdv pctdv odelta_dv zdelta_dv dfdelta_dv pctdelta_dv oE_so4_p zE_so4_p dfE_so4_p pctE_so4_p oE_no3_p zE_no3_p dfE_no3_p pctE_no3_p oE_no2_p zE_no2_p dfE_no2_p pctE_no2_p oT_bext_p zT_bext_p dfT_bext_p pctT_bext_p odv_p zdv_p dfdv_p pctdv_p ); MPRINT(PROGRAM): length n 8. Tracer $9. site $5. ssray mon day year 8. date $8. oNO2 zNO2 dfNO2 pctNO2 oE_so4 zE_so4 dfE_so4 pctE_so4 oE_no3 zE_no3 dfE_no3 pctE_no3 oE_no2 zE_no2 dfE_no2 pctE_no2 oT_bext zT_bext dfT_bext pctT_bext odv zdv dfdv pctdv odelta_dv zdelta_dv dfdelta_dv pctdelta_dv oE_so4_p zE_so4_p dfE_so4_p pctE_so4_p oE_no3_p zE_no3_p dfE_no3_p pctE_no3_p oE_no2_p zE_no2_p dfE_no2_p pctE_no2_p oT_bext_p zT_bext_p dfT_bext_p pctT_bext_p odv_p zdv_p dfdv_p pctdv_p 8.; MPRINT(PROGRAM): set sasd.v2a_bextdv_mrg1; MPRINT(PROGRAM): dfNO2 = zNO2 - oNO2 ; MPRINT(PROGRAM): dfE_so4 = zE_so4 - oE_so4 ; MPRINT(PROGRAM): dfE_no3 = zE_no3 - oE_no3 ; MPRINT(PROGRAM): dfE_no2 = zE_no2 - oE_no2 ; MPRINT(PROGRAM): dfT_bext = zT_bext - oT_bext ; MPRINT(PROGRAM): dfdv = zdv - odv ; MPRINT(PROGRAM): dfdelta_dv = zdelta_dv - odelta_dv ; MPRINT(PROGRAM): dfE_so4_p = zE_so4_p - oE_so4_p ; MPRINT(PROGRAM): dfE_no3_p = zE_no3_p - oE_no3_p ; MPRINT(PROGRAM): dfE_no2_p = zE_no2_p - oE_no2_p ; MPRINT(PROGRAM): dfT_bext_p = zT_bext_p - oT_bext_p ; MPRINT(PROGRAM): dfdv_p = zdv_p - odv_p ; MPRINT(PROGRAM): pctNO2 = 100. * (dfNO2 / oNO2 ); MPRINT(PROGRAM): pctE_so4 = 100. * (dfE_so4 / oE_so4 ); 7 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): pctE_no3 = 100. * (dfE_no3 / oE_no3 ); MPRINT(PROGRAM): pctE_no2 = 100. * (dfE_no2 / oE_no2 ); MPRINT(PROGRAM): pctT_bext = 100. * (dfT_bext / oT_bext ); MPRINT(PROGRAM): pctdv = 100. * (dfdv / odv ); MPRINT(PROGRAM): pctdelta_dv = 100. * (dfdelta_dv / odelta_dv ); MPRINT(PROGRAM): pctE_so4_p = 100. * (dfE_so4_p / oE_so4_p ); MPRINT(PROGRAM): pctE_no3_p = 100. * (dfE_no3_p / oE_no3_p ); MPRINT(PROGRAM): pctE_no2_p = 100. * (dfE_no2_p / oE_no2_p ); MPRINT(PROGRAM): pctT_bext_p = 100. * (dfT_bext_p / oT_bext_p ); MPRINT(PROGRAM): pctdv_p = 100. * (dfdv_p / odv_p ); MPRINT(PROGRAM): format pctNO2 pctE_so4 pctE_no3 pctE_no2 pctT_bext pctdv pctdelta_dv pctE_so4_p pctE_no3_p pctE_no2_p pctT_bext_p pctdv_p 8.2; MPRINT(PROGRAM): run; NOTE: There were 204035 observations read from the data set SASD.V2A_BEXTDV_MRG1. NOTE: The data set SASD.V2B_BEXTDV_MRG1 has 204035 observations and 56 variables. NOTE: DATA statement used (Total process time): real time 0.31 seconds cpu time 0.21 seconds MPRINT(PROGRAM): ******************; MPRINT(PROGRAM): *proc means data=sasd.v2b_bextdv_mrg1 min p5 p10 p25 p50 p75 p90 p95 max mean; MPRINT(PROGRAM): proc means data=sasd.v2b_bextdv_mrg1 min p1 p5 p10 p90 p95 p99 max; MPRINT(PROGRAM): var pctT_bext pctT_bext_p pctdv pctdv_p pctdelta_dv; MPRINT(PROGRAM): run; NOTE: There were 204035 observations read from the data set SASD.V2B_BEXTDV_MRG1. NOTE: The PROCEDURE MEANS printed page 1. NOTE: PROCEDURE MEANS used (Total process time): real time 0.64 seconds cpu time 0.88 seconds MPRINT(PROGRAM): proc univariate data=sasd.v2b_bextdv_mrg1; MPRINT(PROGRAM): var pctT_bext pctdv pctdelta_dv; MPRINT(PROGRAM): histogram / midpoints=-100 to 10 by 2; MPRINT(PROGRAM): run; ERROR: Expecting page 1, got page -1 instead. ERROR: Page validation error while reading SASUSER.PROFILE.CATALOG. NOTE: This SAS session is using a registry in WORK. All changes will be lost at the end of this session. NOTE: Unable to open SASUSER.PROFILE. WORK.PROFILE will be opened instead. NOTE: All profile changes will be lost at the end of the session. WARNING: The specified value of 37.9167 inches for HSIZE= is larger than 8.6776 inches which is the maximum for the device Xcolor. HSIZE is ignored. WARNING: The specified value of 9.3333 inches for VSIZE= is larger than 4.2257 inches which is the maximum for the device Xcolor. VSIZE is ignored. WARNING: The MIDPOINTS= list was extended to accommodate the data. WARNING: Because the data range of at least one histogram cell is much greater than that of the key cell, the number of bins required is very large. NOTE: The PROCEDURE UNIVARIATE printed pages 2-7. NOTE: PROCEDURE UNIVARIATE used (Total process time): real time 46.88 seconds cpu time 0.97 seconds MPRINT(PROGRAM): ******************; 8 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): proc export data=sasd.v2b_bextdv_mrg1 outfile="./v2b_bextdv_mrg1.csv" DBMS=CSV REPLACE; MPRINT(PROGRAM): ADLM; MPRINT(PROGRAM): run; 188 /********************************************************************** 189 * PRODUCT: SAS 190 * VERSION: 9.4 191 * CREATOR: External File Interface 192 * DATE: 26SEP22 193 * DESC: Generated SAS Datastep Code 194 * TEMPLATE SOURCE: (None Specified.) 195 ***********************************************************************/ 196 data _null_; MPRINT(PROGRAM): data _null_; 197 %let _EFIERR_ = 0; /* set the ERROR detection macro variable */ 198 %let _EFIREC_ = 0; /* clear export record count macro variable */ 199 file './v2b_bextdv_mrg1.csv' delimiter=',' DSD DROPOVER lrecl=32767; MPRINT(PROGRAM): file './v2b_bextdv_mrg1.csv' delimiter=',' DSD DROPOVER lrecl=32767; 200 if _n_ = 1 then /* write column names or labels */ 201 do; MPRINT(PROGRAM): if _n_ = 1 then do; 202 put 203 "n" 204 ',' 205 "Tracer" 206 ',' 207 "site" 208 ',' 209 "ssray" 210 ',' 211 "mon" 212 ',' 213 "day" 214 ',' 215 "year" 216 ',' 217 "date" 218 ',' 219 "oNO2" 220 ',' 221 "zNO2" 222 ',' 223 "dfNO2" 224 ',' 225 "pctNO2" 226 ',' 227 "oE_so4" 228 ',' 229 "zE_so4" 230 ',' 231 "dfE_so4" 232 ',' 233 "pctE_so4" 234 ',' 235 "oE_no3" 236 ',' 237 "zE_no3" 238 ',' 9 The SAS System 14:12 Monday, September 26, 2022 239 "dfE_no3" 240 ',' 241 "pctE_no3" 242 ',' 243 "oE_no2" 244 ',' 245 "zE_no2" 246 ',' 247 "dfE_no2" 248 ',' 249 "pctE_no2" 250 ',' 251 "oT_bext" 252 ',' 253 "zT_bext" 254 ',' 255 "dfT_bext" 256 ',' 257 "pctT_bext" 258 ',' 259 "odv" 260 ',' 261 "zdv" 262 ',' 263 "dfdv" 264 ',' 265 "pctdv" 266 ',' 267 "odelta_dv" 268 ',' 269 "zdelta_dv" 270 ',' 271 "dfdelta_dv" 272 ',' 273 "pctdelta_dv" 274 ',' 275 "oE_so4_p" 276 ',' 277 "zE_so4_p" 278 ',' 279 "dfE_so4_p" 280 ',' 281 "pctE_so4_p" 282 ',' 283 "oE_no3_p" 284 ',' 285 "zE_no3_p" 286 ',' 287 "dfE_no3_p" 288 ',' 289 "pctE_no3_p" 290 ',' 291 "oE_no2_p" 292 ',' 293 "zE_no2_p" 294 ',' 295 "dfE_no2_p" 296 ',' 10 The SAS System 14:12 Monday, September 26, 2022 297 "pctE_no2_p" 298 ',' 299 "oT_bext_p" 300 ',' 301 "zT_bext_p" 302 ',' 303 "dfT_bext_p" 304 ',' 305 "pctT_bext_p" 306 ',' 307 "odv_p" 308 ',' 309 "zdv_p" 310 ',' 311 "dfdv_p" 312 ',' 313 "pctdv_p" 314 ; MPRINT(PROGRAM): put "n" ',' "Tracer" ',' "site" ',' "ssray" ',' "mon" ',' "day" ',' "year" ',' "date" ',' "oNO2" ',' "zNO2" ',' "dfNO2" ',' "pctNO2" ',' "oE_so4" ',' "zE_so4" ',' "dfE_so4" ',' "pctE_so4" ',' "oE_no3" ',' "zE_no3" ',' "dfE_no3" ',' "pctE_no3" ',' "oE_no2" ',' "zE_no2" ',' "dfE_no2" ',' "pctE_no2" ',' "oT_bext" ',' "zT_bext" ',' "dfT_bext" ',' "pctT_bext" ',' "odv" ',' "zdv" ',' "dfdv" ',' "pctdv" ',' "odelta_dv" ',' "zdelta_dv" ',' "dfdelta_dv" ',' "pctdelta_dv" ',' "oE_so4_p" ',' "zE_so4_p" ',' "dfE_so4_p" ',' "pctE_so4_p" ',' "oE_no3_p" ',' "zE_no3_p" ',' "dfE_no3_p" ',' "pctE_no3_p" ',' "oE_no2_p" ',' "zE_no2_p" ',' "dfE_no2_p" ',' "pctE_no2_p" ',' "oT_bext_p" ',' "zT_bext_p" ',' "dfT_bext_p" ',' "pctT_bext_p" ',' "odv_p" ',' "zdv_p" ',' "dfdv_p" ',' "pctdv_p" ; 315 end; MPRINT(PROGRAM): end; 316 set SASD.V2B_BEXTDV_MRG1 end=EFIEOD; MPRINT(PROGRAM): set SASD.V2B_BEXTDV_MRG1 end=EFIEOD; 317 format n best12. ; MPRINT(PROGRAM): format n best12. ; 318 format Tracer $9. ; MPRINT(PROGRAM): format Tracer $9. ; 319 format site $5. ; MPRINT(PROGRAM): format site $5. ; 320 format ssray best12. ; MPRINT(PROGRAM): format ssray best12. ; 321 format mon best12. ; MPRINT(PROGRAM): format mon best12. ; 322 format day best12. ; MPRINT(PROGRAM): format day best12. ; 323 format year best12. ; MPRINT(PROGRAM): format year best12. ; 324 format date $8. ; MPRINT(PROGRAM): format date $8. ; 325 format oNO2 best12. ; MPRINT(PROGRAM): format oNO2 best12. ; 326 format zNO2 best12. ; MPRINT(PROGRAM): format zNO2 best12. ; 327 format dfNO2 best12. ; MPRINT(PROGRAM): format dfNO2 best12. ; 328 format pctNO2 8.2 ; MPRINT(PROGRAM): format pctNO2 8.2 ; 329 format oE_so4 best12. ; MPRINT(PROGRAM): format oE_so4 best12. ; 330 format zE_so4 best12. ; MPRINT(PROGRAM): format zE_so4 best12. ; 331 format dfE_so4 best12. ; 11 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): format dfE_so4 best12. ; 332 format pctE_so4 8.2 ; MPRINT(PROGRAM): format pctE_so4 8.2 ; 333 format oE_no3 best12. ; MPRINT(PROGRAM): format oE_no3 best12. ; 334 format zE_no3 best12. ; MPRINT(PROGRAM): format zE_no3 best12. ; 335 format dfE_no3 best12. ; MPRINT(PROGRAM): format dfE_no3 best12. ; 336 format pctE_no3 8.2 ; MPRINT(PROGRAM): format pctE_no3 8.2 ; 337 format oE_no2 best12. ; MPRINT(PROGRAM): format oE_no2 best12. ; 338 format zE_no2 best12. ; MPRINT(PROGRAM): format zE_no2 best12. ; 339 format dfE_no2 best12. ; MPRINT(PROGRAM): format dfE_no2 best12. ; 340 format pctE_no2 8.2 ; MPRINT(PROGRAM): format pctE_no2 8.2 ; 341 format oT_bext best12. ; MPRINT(PROGRAM): format oT_bext best12. ; 342 format zT_bext best12. ; MPRINT(PROGRAM): format zT_bext best12. ; 343 format dfT_bext best12. ; MPRINT(PROGRAM): format dfT_bext best12. ; 344 format pctT_bext 8.2 ; MPRINT(PROGRAM): format pctT_bext 8.2 ; 345 format odv best12. ; MPRINT(PROGRAM): format odv best12. ; 346 format zdv best12. ; MPRINT(PROGRAM): format zdv best12. ; 347 format dfdv best12. ; MPRINT(PROGRAM): format dfdv best12. ; 348 format pctdv 8.2 ; MPRINT(PROGRAM): format pctdv 8.2 ; 349 format odelta_dv best12. ; MPRINT(PROGRAM): format odelta_dv best12. ; 350 format zdelta_dv best12. ; MPRINT(PROGRAM): format zdelta_dv best12. ; 351 format dfdelta_dv best12. ; MPRINT(PROGRAM): format dfdelta_dv best12. ; 352 format pctdelta_dv 8.2 ; MPRINT(PROGRAM): format pctdelta_dv 8.2 ; 353 format oE_so4_p best12. ; MPRINT(PROGRAM): format oE_so4_p best12. ; 354 format zE_so4_p best12. ; MPRINT(PROGRAM): format zE_so4_p best12. ; 355 format dfE_so4_p best12. ; MPRINT(PROGRAM): format dfE_so4_p best12. ; 356 format pctE_so4_p 8.2 ; MPRINT(PROGRAM): format pctE_so4_p 8.2 ; 357 format oE_no3_p best12. ; MPRINT(PROGRAM): format oE_no3_p best12. ; 358 format zE_no3_p best12. ; MPRINT(PROGRAM): format zE_no3_p best12. ; 359 format dfE_no3_p best12. ; MPRINT(PROGRAM): format dfE_no3_p best12. ; 360 format pctE_no3_p 8.2 ; 12 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): format pctE_no3_p 8.2 ; 361 format oE_no2_p best12. ; MPRINT(PROGRAM): format oE_no2_p best12. ; 362 format zE_no2_p best12. ; MPRINT(PROGRAM): format zE_no2_p best12. ; 363 format dfE_no2_p best12. ; MPRINT(PROGRAM): format dfE_no2_p best12. ; 364 format pctE_no2_p 8.2 ; MPRINT(PROGRAM): format pctE_no2_p 8.2 ; 365 format oT_bext_p best12. ; MPRINT(PROGRAM): format oT_bext_p best12. ; 366 format zT_bext_p best12. ; MPRINT(PROGRAM): format zT_bext_p best12. ; 367 format dfT_bext_p best12. ; MPRINT(PROGRAM): format dfT_bext_p best12. ; 368 format pctT_bext_p 8.2 ; MPRINT(PROGRAM): format pctT_bext_p 8.2 ; 369 format odv_p best12. ; MPRINT(PROGRAM): format odv_p best12. ; 370 format zdv_p best12. ; MPRINT(PROGRAM): format zdv_p best12. ; 371 format dfdv_p best12. ; MPRINT(PROGRAM): format dfdv_p best12. ; 372 format pctdv_p 8.2 ; MPRINT(PROGRAM): format pctdv_p 8.2 ; 373 do; MPRINT(PROGRAM): do; 374 EFIOUT + 1; MPRINT(PROGRAM): EFIOUT + 1; 375 put n @; MPRINT(PROGRAM): put n @; 376 put Tracer $ @; MPRINT(PROGRAM): put Tracer $ @; 377 put site $ @; MPRINT(PROGRAM): put site $ @; 378 put ssray @; MPRINT(PROGRAM): put ssray @; 379 put mon @; MPRINT(PROGRAM): put mon @; 380 put day @; MPRINT(PROGRAM): put day @; 381 put year @; MPRINT(PROGRAM): put year @; 382 put date $ @; MPRINT(PROGRAM): put date $ @; 383 put oNO2 @; MPRINT(PROGRAM): put oNO2 @; 384 put zNO2 @; MPRINT(PROGRAM): put zNO2 @; 385 put dfNO2 @; MPRINT(PROGRAM): put dfNO2 @; 386 put pctNO2 @; MPRINT(PROGRAM): put pctNO2 @; 387 put oE_so4 @; MPRINT(PROGRAM): put oE_so4 @; 388 put zE_so4 @; MPRINT(PROGRAM): put zE_so4 @; 389 put dfE_so4 @; 13 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): put dfE_so4 @; 390 put pctE_so4 @; MPRINT(PROGRAM): put pctE_so4 @; 391 put oE_no3 @; MPRINT(PROGRAM): put oE_no3 @; 392 put zE_no3 @; MPRINT(PROGRAM): put zE_no3 @; 393 put dfE_no3 @; MPRINT(PROGRAM): put dfE_no3 @; 394 put pctE_no3 @; MPRINT(PROGRAM): put pctE_no3 @; 395 put oE_no2 @; MPRINT(PROGRAM): put oE_no2 @; 396 put zE_no2 @; MPRINT(PROGRAM): put zE_no2 @; 397 put dfE_no2 @; MPRINT(PROGRAM): put dfE_no2 @; 398 put pctE_no2 @; MPRINT(PROGRAM): put pctE_no2 @; 399 put oT_bext @; MPRINT(PROGRAM): put oT_bext @; 400 put zT_bext @; MPRINT(PROGRAM): put zT_bext @; 401 put dfT_bext @; MPRINT(PROGRAM): put dfT_bext @; 402 put pctT_bext @; MPRINT(PROGRAM): put pctT_bext @; 403 put odv @; MPRINT(PROGRAM): put odv @; 404 put zdv @; MPRINT(PROGRAM): put zdv @; 405 put dfdv @; MPRINT(PROGRAM): put dfdv @; 406 put pctdv @; MPRINT(PROGRAM): put pctdv @; 407 put odelta_dv @; MPRINT(PROGRAM): put odelta_dv @; 408 put zdelta_dv @; MPRINT(PROGRAM): put zdelta_dv @; 409 put dfdelta_dv @; MPRINT(PROGRAM): put dfdelta_dv @; 410 put pctdelta_dv @; MPRINT(PROGRAM): put pctdelta_dv @; 411 put oE_so4_p @; MPRINT(PROGRAM): put oE_so4_p @; 412 put zE_so4_p @; MPRINT(PROGRAM): put zE_so4_p @; 413 put dfE_so4_p @; MPRINT(PROGRAM): put dfE_so4_p @; 414 put pctE_so4_p @; MPRINT(PROGRAM): put pctE_so4_p @; 415 put oE_no3_p @; MPRINT(PROGRAM): put oE_no3_p @; 416 put zE_no3_p @; MPRINT(PROGRAM): put zE_no3_p @; 417 put dfE_no3_p @; MPRINT(PROGRAM): put dfE_no3_p @; 418 put pctE_no3_p @; 14 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): put pctE_no3_p @; 419 put oE_no2_p @; MPRINT(PROGRAM): put oE_no2_p @; 420 put zE_no2_p @; MPRINT(PROGRAM): put zE_no2_p @; 421 put dfE_no2_p @; MPRINT(PROGRAM): put dfE_no2_p @; 422 put pctE_no2_p @; MPRINT(PROGRAM): put pctE_no2_p @; 423 put oT_bext_p @; MPRINT(PROGRAM): put oT_bext_p @; 424 put zT_bext_p @; MPRINT(PROGRAM): put zT_bext_p @; 425 put dfT_bext_p @; MPRINT(PROGRAM): put dfT_bext_p @; 426 put pctT_bext_p @; MPRINT(PROGRAM): put pctT_bext_p @; 427 put odv_p @; MPRINT(PROGRAM): put odv_p @; 428 put zdv_p @; MPRINT(PROGRAM): put zdv_p @; 429 put dfdv_p @; MPRINT(PROGRAM): put dfdv_p @; 430 put pctdv_p ; MPRINT(PROGRAM): put pctdv_p ; 431 ; MPRINT(PROGRAM): ; 432 end; MPRINT(PROGRAM): end; 433 if _ERROR_ then call symputx('_EFIERR_',1); /* set ERROR detection macro variable */ MPRINT(PROGRAM): if _ERROR_ then call symputx('_EFIERR_',1); 434 if EFIEOD then call symputx('_EFIREC_',EFIOUT); MPRINT(PROGRAM): if EFIEOD then call symputx('_EFIREC_',EFIOUT); 435 run; MPRINT(PROGRAM): run; NOTE: The file './v2b_bextdv_mrg1.csv' is: Filename=/work/ROMO/users/hou/400.1/EPA-R6-2016/2.postproc/xTESTING/1.master/yQA.zeroRGN/4/v2b_bextdv_mrg1.csv, Owner Name=hou,Group Name=romo, Access Permission=-rw-r--r--, Last Modified=26Sep2022:14:13:02 NOTE: 204036 records were written to the file './v2b_bextdv_mrg1.csv'. The minimum record length was 317. The maximum record length was 547. NOTE: There were 204035 observations read from the data set SASD.V2B_BEXTDV_MRG1. NOTE: DATA statement used (Total process time): real time 2.58 seconds cpu time 2.48 seconds 204035 records created in ./v2b_bextdv_mrg1.csv from SASD.V2B_BEXTDV_MRG1. NOTE: "./v2b_bextdv_mrg1.csv" file was successfully created. NOTE: PROCEDURE EXPORT used (Total process time): real time 2.68 seconds cpu time 2.53 seconds 15 The SAS System 14:12 Monday, September 26, 2022 MPRINT(PROGRAM): ; MPRINT(PROGRAM): ******************; MLOGIC(PROGRAM): Ending execution. ERROR: Errors printed on page 7. NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414 NOTE: The SAS System used: real time 53.60 seconds cpu time 7.15 seconds