###LOAD GHCN V2 STATION INVENTORY
	#this can be collated once and saved
# [1] "40371951001 KEMANO,BC                       53.57 -127.93   70  793R   -9MVxxCO10x-9TUNDRA          B"
# [2] "40371951002 KILDALA,BC                      53.83 -128.48   30  554R   -9MVxxCO 1x-9COOL CONIFER    A"
# [3] "40371951003 KITIMAT MISSION,BC              53.98 -128.65    6  211R   -9MVxxCO 1x-9COOL CONIFER    A"
 
	url<-"d:/climate/data/station/ghcn/stations.tab"
	test<-try(load(url) ) 
	if(class(test)=="try-error") {
		url<-"ftp://ftp.ncdc.noaa.gov/pub/data/ghcn/v2/v2.temperature.inv"
		test<-readLines(url);length(test) #7280
		stations<-data.frame(as.numeric(paste(substr(test,1,8),"0",substr(test,9,11),sep="")),as.character(sub(" +$","",substr(test,12,42) )))
		names(stations)<-c("id","site")
		stations$wmo<-as.numeric(substr(test,4,8))
		stations$version<-stations$id%%100
		stations$lat<-as.numeric(substr(test,43,49))
		stations$long<-as.numeric(substr(test,50,58))
		stations$pop<-as.numeric(substr(test,69,73))
		stations$country<- as.numeric(substr(test,1,3))
		}
	
	#id0<-29612;stations[stations$wmo==id0,]