###############
##GHCN LAND 
################

   #ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/anom/anom-grid2-1880-current.dat.gz ##STALE
              #updated to June 2009 2434 MB 6/8/2009 7:00:00
    ##ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/grid/grid_mean_temp_1880_current.dat.gz  ##STALE 
               #updated to June 2009 2434 MB 6/8/2009 7:00:00
   ## ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/anom/README_GRID_TEMP.txt

  #readme: ftp://ftp.ncdc.noaa.gov/pub/data/er-ghcn-sst/er-ghcn-sst.readme.txt anomaly reference not specified
  #readme:  http://www.ncdc.noaa.gov/oa/climate/research/sst/ersstv3.php  1971-2000 reference

	# download.file("ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/grid/grid_1880_2008.dat.gz", "anom.gz",mode="wb");
	# download.file("ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/anom/anom-grid2-1880-current.dat.gz", "anom.gz",mode="wb");
	# download.file("ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/grid/grid_mean_temp_1880_current.dat.gz", "anom.gz",mode="wb");
 

get.ghcn=function() {
	year=date();n=nchar(year)
	year=as.numeric(substr(year,n-3,n))
	download.file(paste("ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/grid/grid_1880_",year,".dat.gz",sep=""),"temp.dat",mode="wb");
	handle <- gzfile("temp.dat");
	fred=readLines(handle,n=-1)
	fred=gsub("-32768","    NA",fred)
	writeLines(fred,"temp.dat")
	x=scan("temp.dat")
	N=length(x); N # 4028482 item
		#N/2594 #1527 
	x=array(x,dim=c(2594,N/2594) )
	ghcn=ts( t( x[3:2594,]),start=1880,freq=12)/100 #eliminate month year #start at 1880 1
	  #tsp(noaa)
	  #1880.000 2007.167   12.000
	return(ghcn)
     }