####################################################
# 2014-03-27 updated to work with data from
#           state subfolders
#
# 2014-02-05 Read data from tables 
# with monthly mean data and create a
# line, or barplot
###################################################
# clear memory
rm(list=ls())

###################################################
# User's program control parameters
###################################################
# state id
state<-"NY"
# station
stationname<-"USW00094789"
# variable
varname<-"tavg"
# first and last year for averaging
yr1<-1981
yr2<-2010
type<-"bar"


###################################################
# Main program 
###################################################
# 2014-03-27 updated to work with state subfolder
# (also corrected the file name to match 
#  file names from climatology.R)
infile<-paste("data/",state,"/",stationname,"_",varname,"_mon_mean_climc_",yr1,"-",yr2,".csv",sep='')


print ("plot monthly mean climatology")
print(paste("Station :",stationname))
print(paste("variable :",varname))
buffer<-read.csv(infile)
month<-buffer[,1]
xmean<-buffer[,2]
nsum<-buffer[,3]



if (type=="line") {
  
  plot(xmean,typ='l',ylab='tavg [C]')
}

if (type=="bar") {
  
  barplot(xmean,names.arg=buffer[,1],ylab='tavg [C]')
}