###MAKE AND SAVE MSU DATA AVERAGES FOR WEBSTER BASIN
   ##Webster et al 2005
 #We define the ocean basins that support tropical cyclone development as follows: 
  	#North Atlantic (90- to 20-W, 5- to 25-N), 
	#western North Pacific (120- to 180-E, 5- to 20-N), 
	#eastern North Pacific (90- to 120-W, 5- to 20-N), 
	#South Indian (50- to 115-E, 5--20-S), 
	#North Indian (55- to 90-E, 5--20-N), and 
	#Southwest Pacific (155- to 180-E, 5- to 20-S).

  #Figure 1 LEgend
	#NATL: 90- to 20-E, 5- to 25-N, June-October), 
	#the Western Pacific Ocean (WPAC: 120- to 180-E, 5-to 20-N, May-December), 
	#the East Pacific Ocean (EPAC: 90- to 120-W, 5- to 20-N, June-October),
	#the Southwest Pacific Ocean (SPAC: 155- to 180-E, 5- to 20-S,December-April),
	# the North Indian Ocean (NIO: 55- to 90-E, 5- to 20-N, April-May and September-November), 
	#and the South Indian Ocean (SIO: 50- to 115-E, 5- to 20-S, November-April).

	basin<-rep(list(NA),6);names(basin)<-c("NATL","WNP","ENP","SIO","NIO","SWP")
	basin$NATL<- c(t(outer( seq(7.5,17.5,5),seq(-87.5, -22.5,5),jones)))
	basin$WNP<- c(t(outer( seq(7.5,17.5,5),seq(122.5, 177.5,5),jones)))
	basin$ENP<- c(t(outer( seq(7.5,17.5,5),seq(-117.5, -92.5,5),jones)))
	basin$SIO<- c(t(outer( seq(-17.5,-7.5,5),seq(52.5, 112.5,5),jones)))
	basin$NIO<- c(t(outer( seq(7.5,17.5,5),seq(57.5, 87.5,5),jones)))
	basin$SWP<- c(t(outer( seq(-17.5,-7.5,5),seq(157.5, 177.5,5),jones)))

	season<-list(c(6:10),c(5:12),c(6:10),c(1:4,12),c(4:5,9:11),c(1:4,11:12))
	NHindex<-c(TRUE,TRUE,TRUE,FALSE,TRUE,FALSE)


  ##LOAD SATELLITE TEMPERATURE
	load("c:/climate/data/satellite/msu.tab"); dim(msu)# 332 2592

  ##MAKE AND SAVE BASIN AVERAGES 

	basin.average<-array(NA,dim=c(332,6))
	dimnames(basin.average)[[2]]<-names(basin)
	for ( i in 1:6) {basin.average[,i]<-apply(msu[,basin[[i]]],1,mean) }
	basin.average<-ts(basin.average[2:313,],start=c(1979,1),freq=12)

	year<-seq(tsp(basin.average)[1],tsp(basin.average)[2],1/tsp(basin.average)[3])
	nf<-layout(array(1:6,dim=c(6,1)),heights=c(1.1,1,1,1,1,1.3))
	par(mar=c(0,4,1,1))
	plot(year,basin.average[,1],type="l",axes=FALSE,ylab=names(basin)[1])
	axis(side=1,labels=FALSE);axis(side=2,las=1);box();abline(h=0,lty=2)
	par(mar=c(0,4,0,1))
	for (i in 2:5){	plot(year,basin.average[,i],type="l",axes=FALSE,ylab=names(basin)[i])
	axis(side=1,labels=FALSE);axis(side=2,las=1);box();abline(h=0,lty=2)}
	par(mar=c(3,4,0,1))
	plot(year,basin.average[,6],type="l",axes=FALSE,ylab=names(basin)[6])
	axis(side=1);axis(side=2,las=1);box();abline(h=0,lty=2)

	write.table(cbind(time(basin.average),basin.average),file="c:/climate/data/hurricane/basin.average.msu.txt",sep="\t")