#NOAA 

get_noaa=function(firma="loti",hemi="gl") {
	aliasf=c(land="land",ocean="ocean",loti="land_ocean")
	aliash=c(gl="90S.90N",nh="00N.90N",sh="90S.00N")
	loc=paste("ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.",aliasf[firma],".",aliash[hemi],".df_1901-2000mean.dat",sep="")
	noaa.glb=read.table(loc)
	noaa.glb=ts(noaa.glb[,3],start=c(1880,1),freq=12)
	noaa.glb[noaa.glb== -999]=NA
	noaa.glb=window(noaa.glb,end= max(time(noaa.glb)[!is.na(noaa.glb)]) )
	return(noaa.glb)
}

 

noaa=get_noaa(firma="loti") 
	
	#tropical ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/all.land_and_ocean.20S.20N.dat
	# ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_and_ocean.20S.20N.df_1961-1990mean.dat NOT UP TO DATA
	#ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_and_ocean.20S.20N.df_1901-2000mean.dat  2006 11 

	#ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/grid/grid_1880-2008_RVose.dat.gz  for gridded
	#see spaghetti.tropical