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"))

bm1<-read.table(file="bmfile1.txt",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","spd","dir"))

bm2<-read.table(file="bmfile2.txt",header=F,na.strings="-999",col.name=c("yr","mo","dy","hr","spd","dir"))

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



x11(width=6,height=3)
par(mar=c(3,3,2,3))
par(mgp=c(3,0.5,0))
par(cex=.9)
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 Consistency",font=2,side=2,line=2)
mtext("Hour",font=2,side=1,line=1.5,)
axis(2,las=2)
axis(1,c(1,4,8,12,16,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,ylim=c(0,360))
axis(4,c(0,90,180,270,360),las=2)
mtext("Wind Direction (degrees)",side=4,font=2,line=2)
legend(3,365,legend=c("Avg Wspd","Avg Wdir","Consistency"),lty=c(1,NA,2),pch=c(4,4,20),col=c(1,4,3),cex=0.9)

#hly tmp and spd from aes_pg:

h<-hlyp$tmp/10
hr<-hlyp$hr
s<-hlyp$spd*1000/3600

i<-1
for(i in 1:length(s))
if (s[i]==0){
s[i]<-NA}


x11(height=4,width=4)
par(mar=c(3,3,2,3))
par(mgp=c(3,0.5,0))
par(cex=.9)
par(tcl=-0.3)

mat<-(matrix(1:2,2,1))
layout(mat)
layout.show(2)

boxplot(split(h,hr),ylim=c(0,30),axes=F)
axis(2,las=2)
axis(1,c(0,2,4,6,8,10,12,14,16,18,20,22,24))
box()
mtext("Temperature (Degrees)",font=2,side=2,line=2)
mtext("Hour",font=2,side=1,line=1.5,)

boxplot(split(s,hr),axes=F)
axis(2,las=2)
axis(1,c(0,2,4,6,8,10,12,14,16,18,20,22,24))
box()
mtext("Wind Speed (m/s)",font=2,side=2,line=2)
mtext("Hour",font=2,side=1,line=1.5,)


#Length of day graphs:

hist(all$n,breaks=seq(from=0,to=12,by=1),axes=F)
axis(2,las=2)
axis(1,at=c(0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5),labels=c(1,2,3,4,5,6,7,8,9,10,11,12))

,pch=c(3),labels=c(1,2,3,4,5,6,7,8,9,10,11,12))

x11(width=3,height=3)
hist(rte$n,ylim=c(0,40),xlim=c(0,12))
x11(width=3,height=3)
hist(rtw$n)
x11(width=3,height=3)
hist(roh$n)
x11(width=3,height=3)
hist(all$n)