##COLLATION STEIG

#http://faculty.washington.edu/steig/nature09data/ReadMe.html

	
	grid.info=1:5509;grid.info=data.frame(grid.info);names(grid.info)="id"
	grid.info$lat=scan("http://faculty.washington.edu/steig/nature09data/data/Tir_lats.txt")
	grid.info$long=scan("http://faculty.washington.edu/steig/nature09data/data/Tir_lons.txt")
	temp=grid.info$long>180
	grid.info$long[temp]= -( 360- grid.info$long[temp])


  	       
  #ant_recon
	download.file("http://faculty.washington.edu/steig/nature09data/ant_recon.txt","temp.dat",mode="wb")
	grid=scan("temp.dat",n= -1) #  3305400
	  #600*5509 #[1] 3305400
	recon=ts(t(array(grid,dim=c(5509,600))),start=c(1957,1),freq=12)	
	dimnames(recon)[[2]]=1:5509
	save(recon,file="d:/climate/data/steig/recon.tab")
	load("d:/climate/data/steig/recon.tab")

  #ant_recon_corr 
 	download.file("http://faculty.washington.edu/steig/nature09data/data/awsreconcorrected.txt","temp.dat",mode="wb")
	grid=scan("temp.dat",n= -1) #  3305400
	recon=ts(t(array(grid,dim=c(length(grid)/600,600))),start=c(1957,1),freq=12)	
	dimnames(recon)[[2]]=1:63
	save(recon,file="d:/climate/data/steig/recon.corr.tab") #600 5509 
	recon_corr=recon


  #ant_recon_detrend
	download.file("http://faculty.washington.edu/steig/nature09data/data/ant_recon_detrend.txt","temp.dat",mode="wb")
	   #doesn't exist
	grid=scan("temp.dat",n= -1) #  3305400
	  #600*5509 #[1] 3305400
	recon_detrend=ts(t(array(grid,dim=c(5509,600))),start=c(1957,1),freq=12)	
	dimnames(recon_detrend)[[2]]=1:5509
	save(recon_detrend,file="d:/climate/data/steig/recon_detrend.tab")

  #http://faculty.washington.edu/steig/nature09data/data/ant_recon_pca.txt
	download.file("http://faculty.washington.edu/steig/nature09data/data/ant_recon_detrend.txt","temp.dat",mode="wb")
	   #doesn't exist
	grid=scan("temp.dat",n= -1) #  3305400
	  #600*5509 #[1] 3305400
	recon_pca=ts(t(array(grid,dim=c(5509,600))),start=c(1957,1),freq=12)	
	dimnames(recon_pca)[[2]]=1:5509
	save(recon_pca,file="d:/climate/data/steig/recon_pca.tab")


  #http://faculty.washington.edu/steig/nature09data/ant_recon_aws.txt
	download.file("http://faculty.washington.edu/steig/nature09data/data/ant_recon_aws.txt","temp.dat",mode="wb")
	grid=scan("temp.dat",n= -1) #  37800
	  #600*63 #[1] 3305400
	recon_aws=ts(t(array(grid,dim=c(63,600))),start=c(1957,1),freq=12)	
	dimnames(recon_aws)[[2]]=Info$aws_grid$id
	save(recon_aws,file="d:/climate/data/steig/recon_aws.tab")

	lat=scan("http://faculty.washington.edu/steig/nature09data/data/aws_lats.txt")
	long=scan("http://faculty.washington.edu/steig/nature09data/data/aws_lons.txt")



	
	plot(Info$aws_grid$long,Info$aws_grid$lat)

	dimnames(recon_aws)[[2]]=paste(Info$aws_grid$id)
	
	x=recon_aws[,"89324"]
	y=Data$aws[,"89324"]
	y=Anom[,"89324"]

	ts.plot(x-y,type="h")
	






#AVHRR
#	AVHRR Polar Pathfinder Twice-Daily 1.25 km EASE-Grid Composites  
#	     April 1992 - January 1996 for the Southern Hemisphere.
#		Clear sky surface broadband albedo and skin temperature

#	AVHRR Polar Pathfinder Twice-Daily 5 km EASE-Grid Composites
#	 	July 1981 - August 1998
#		Clear sky surface broadband albedo and skin temperature

#	AVHRR Polar Pathfinder Twice-Daily 25 km EASE-Grid Composites
		# Clear sky surface broadband albedo and skin temperature
		# Average albedo and temperature
#		July 1981 - December 2000