options linesize = 75 pagesize = 60 nodate nonumber ;

data depo ;
infile 'depoRepCO.asp' firstobs = 6 ;
input SiteID  $QUOTE6. @9 Per $QUOTE8. Year  Crit1 Crit2 Crit3 Crit4 Ca Mg  
K  Na  NH4  NO3 InorgN Cl SO4 HLab  HField  Svol  Ppt  Pct  ValidF ValidL Days  @201 Date1 $QUOTE12. @217 Date2 $QUOTE12.  ;
drop Per Year Crit1-Crit4 ;
sdate = input(Date1, mmddyy10.) ;   * convert from char var to date var ;
edate = input(Date2, mmddyy10.) ;   
daysop = edate - sdate ;
format sdate edate date8. ;
run ;

* proc print data  = depo (obs=10) ;
* run ;

data sites ;
infile 'stateCO.asp' firstobs = 19 missover ;
input @13 SiteID $ @20 sitename $18. @40 strtdate mmddyy10. @53 stopdate mmddyy10. @68 elev ;
if strtdate ne . ;   * subsetting if: exclude observations meeting condition ;
format strtdate stopdate date8. ;
run ;

* proc print ;
* run ;

data combined ;
merge depo sites ;
run ;

proc print ;
run ;

proc plot data = combined ; * text plot ;
plot SO4 * elev = '.' /hpos = 40 vpos = 20 ; * syntax: y-var * x-var = 'plotting symbol' ;
title1 'Sulfuric Acid vs. Elevation' ;
title2 'Colorado NADP sites, Annual Data for 2000' ;
footnote 'Sulfuric Acid in kg/ha' ;
label SO4 = 'Sulfuric Acid' elev = 'Elevation in Meters' ;
run ;

proc gplot data = combined ;   * SAS/Graph plot ;
plot SO4 * elev = '.' ; * syntax: y-var * x-var = 'plotting symbol' ;
title1 'Sulfuric Acid vs. Elevation' ;
title2 'Colorado NADP sites, Annual Data for 2000' ;
footnote 'Sulfuric Acid in kg/ha' ;
label SO4 = 'Sulfuric Acid' elev = 'Elevation in Meters' ;
run ;

proc plot data = combined  ; * text plot ;
plot SO4 * elev = '.' days * elev = '*' /hpos = 40 vpos = 20 overlay ; 
* syntax: y-var * x-var = 'plotting symbol' ;
title1 'Sulfuric Acid vs. Elevation' ;
title2 'Colorado NADP sites, Annual Data for 2000' ;
footnote 'Sulfuric Acid in kg/ha' ;
label SO4 = 'Sulfuric Acid' elev = 'Elevation in Meters' ;
run ;

proc gplot data = combined ;   * SAS/Graph plot ;
plot SO4 * elev = '.' days * elev = '*' /overlay ; 
* syntax: y-var * x-var = 'plotting symbol' ;
title1 'Sulfuric Acid vs. Elevation' ;
title2 'Colorado NADP sites, Annual Data for 2000' ;
footnote 'Sulfuric Acid in kg/ha' ;
label SO4 = 'Sulfuric Acid' elev = 'Elevation in Meters' ;
run ;