Reanalasis Graphs: original files are on stratus in /home/lindsayv/grads/reanal_comps

all<-read.table(file="peakall.out",header=F,na.strings="-999",col.name=c("dat","jd","yr","mo","dy","n","tmp","tmp3","ppt","spd","dir"))

noc<-read.table(file="noClass.both",header=F,na.strings="-999",col.name=c("dat","jd","yr","mo","dy","n","tmp","tmp3","ppt","spd","dir"))

rtw<-read.table(file="rtw.both",header=F,na.strings="-999",col.name=c("dat","jd","yr","mo","dy","n","tmp","tmp3","ppt","spd","dir"))

rte<-read.table(file="rte.both",header=F,na.strings="-999",col.name=c("dat","jd","yr","mo","dy","n","tmp","tmp3","ppt","spd","dir"))

roh<-read.table(file="roh.both",header=F,na.strings="-999",col.name=c("dat","jd","yr","mo","dy","n","tmp","tmp3","ppt","spd","dir"))

allhly<-read.table(file="peakall.hly",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","n","tmp","tmp3","ppt","spd","dir","sol"))

nochly<-read.table(file="noClass.hly",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","n","tmp","tmp3","ppt","spd","dir","sol"))

rtwhly<-read.table(file="rtw.hly",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","n","tmp","tmp3","ppt","spd","dir","sol"))

rtehly<-read.table(file="rte.hly",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","n","tmp","tmp3","prs","spd","dir","sol"))

rohhly<-read.table(file="roh.hly",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","n","tmp","tmp3","ppt","spd","dir","sol"))

### Run wind_vec_mean.txt on hourly data and then:



x11(width=5,height=6)
par(mar=c(3,3,2,3))
par(mgp=c(3,0.5,0))
par(cex=1.3)
par(tcl=-0.3)
plot(sMean,ylim=c(0,3),xlim=c(0,24),ann=F,pch=4,axes=F)
points(sMean,type="l")
mgp=c(3,0.5,0)
mtext("Wind Speed(m/s) and Consistancy",side=2,line=1.2,adj=0.2,cex=1.6)
mtext("Hours",side=1,line=1.2,cex=1.6)
axis(2)
axis(1,c(1,5,10,15,20,24))
box()
points(cons,col=3,pch=20)
points(cons,type="l",lty=2,col=3)
par(new=T)
par(mgp=c(3,0.5,0))
plot(d,axes=F,col=4,ann=F,pch=4)
axis(4,col=4)
mtext("Wind Direction (degrees)",side=4,line=1.2,adj=0.35,cex=1.6)
legend(1.5,325,legend=c("Avg Wspd","Avg Wdir","Consistancy"),lty=c(1,NA,2),pch=c(4,4,20),col=c(1,4,3))