Nothing
get_gsod_stn <-
function(){
USAF=NULL
END=NULL
BEGIN=NULL
GSOD_history=""
load(paste(path.package("EcoHydRology"), "data/GSOD_history.rda", sep = "/"))
# data(GSOD_history,package="EcoHydRology")
tmpdir=readline("Please enter a temp directory to store the datafiles in, remember this directory when parsing with build_gsod_forcing_data()? \n")
dir.create(tmpdir)
stnlist=GSOD_history
ANSWER="n"
while(ANSWER=="n"){
stnreq=readline("What is the station name you want to look for? \n")
test=grep(stnreq,stnlist$STATION.NAME)
print(stnlist[test,])
ANSWER <- readline("Good enough? ")
}
ANSWER <- readline("Given the start and end years, would you like to update the stnlist file from NDCD? (may take a few minutes) ")
## a better version would check the answer less cursorily, and
## perhaps re-prompt
if (substr(ANSWER, 1, 1) == "y")
stnlist=read.csv(url("ftp://ftp.ncdc.noaa.gov/pub/data/gsod/ish-history.csv"))
#
# Add one more request to make sure the number of years is ok as some stations have really long histories.
#
#
working_tmp <- stnlist$STATION.NAME
has_stn<- working_tmp %in% grep(stnreq, working_tmp, value=TRUE)
dl_wbans=subset(stnlist, has_stn , USAF:END)
start_year=min(trunc(subset(stnlist, has_stn, BEGIN)/10000))
end_year=max(trunc(subset(stnlist, has_stn , END)/10000))
start_year_rec<- readline(paste("Date range available is",start_year,"to",end_year,", what start year would you like?"))
end_year_rec<- readline(paste("Date range available is",start_year,"to",end_year,", what end year would you like?"))
for (wbans in row.names(dl_wbans)){
start_year=max(start_year_rec,substr(dl_wbans[wbans,"BEGIN"],1,4))
end_year=min(end_year_rec,substr(dl_wbans[wbans,"END"],1,4))
print(c(start_year,end_year))
if(start_year > end_year) next()
for (year in start_year:end_year){
dlurl=paste("ftp://ftp.ncdc.noaa.gov/pub/data/gsod/",year,"/",dl_wbans[wbans,1],"-",dl_wbans[wbans,2],"-",year,".op.gz",sep="")
dlfile=paste(tmpdir,"/",dl_wbans[wbans,1],"-",dl_wbans[wbans,2],"-",year,".op.gz",sep="")
print(dlurl)
try(download.file(dlurl,dlfile),silent=T)
}
}
files=dir(tmpdir,"*op.gz")
print(files)
ANSWER <- readline("Good enough? ")
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.