#2008 Area and Extent source("http://www.climateaudit.org/scripts/seaice/functions.seaice.txt") download.file("http://www.climateaudit.org/data/seaice/Arctic.tab","temp.dat",mode="wb");load("temp.dat") dim(Arctic) #may 2009 #f15 stopped in Feb 2009 after being up since Jan 1, 2008 #f13 unavailable for first 3 months of 2008; substituted back Arctic=update.nsidc(working=Arctic,end0= Sys.Date() ) 11003 2009_02_14_f15.tab 14289 2009 02 14 f15 45 12.68360 14.25993 NA 0.3107784 14289 11004 2009_02_15_f15.tab 14290 2009 02 15 f15 46 12.16574 13.53851 NA 0.3107784 14290 11005 2009_02_16_f15.tab 14291 2009 02 16 f15 47 NA NA NA NA 14291 11006 2009_02_17_f15.tab 14292 2009 02 17 f15 48 NA NA NA NA 14292 k=11005 day=working$dd[k];month=working$mm[k];year=working$year[k];sat=working$sat[k] loc= "ftp://sidads.colorado.edu/pub/DATASETS/seaice/polar-stereo/nasateam/near-real-time/north" ; orig.file=paste("nt_",year,pastedate(month),pastedate(day),"_",sat,"_nrt_n.bin",sep="") url=file.path(loc,orig.file) #ftp://sidads.colorado.edu/pub/DATASETS/seaice/polar-stereo/nasateam/near-real-time/north/nt_20081215_f13_nrt_n.bin #ftp://sidads.colorado.edu/pub/DATASETS/seaice/polar-stereo/nasateam/near-real-time/north/nt_20090216_f15_nrt_n.bin" #DAILY ANOMALY temp= Arctic$year<2009 #anomaly.seaice=tapply(Arctic$extent[temp],Arctic$jday[temp],mean,na.rm=TRUE) anomaly.seaice=tapply(Arctic$area[temp]+Arctic$hole[temp],Arctic$jday[temp],mean,na.rm=TRUE) Arctic$anomaly=factor(Arctic$jday) levels(Arctic$anomaly)=anomaly.seaice Arctic$anomaly=as.numeric(as.character(Arctic$anomaly)) #Arctic$anomaly=Arctic$extent-Arctic$anomaly Arctic$anomaly=Arctic$area+Arctic$hole-Arctic$anomaly temp= Arctic$mm=="01" &Arctic$jday==1 & Arctic$year %%10 ==0 sum(temp) index=Arctic$julian[temp] index=c(index,max(Arctic$julian)+1) temp=!is.na(Arctic$anomaly) par(mar=c(3,4,2,1)) plot(Arctic$julian[temp],Arctic$anomaly[temp],type="l",axes=FALSE) axis(side=1,at=index,labels=paste(seq(1980,2010,10) ) ) axis(side=2,las= 1);box() temp=!is.na(Arctic$hole) par(mar=c(3,4,2,1)) plot(Arctic$julian[temp],Arctic$hole[temp],type="l",axes=FALSE) axis(side=1,at=index,labels=paste(seq(1980,2010,10) ) ) axis(side=2,las= 1);box() plot(Arctic$julian[temp],Arctic$hole[temp]+Arctic$area[temp],type="l",axes=FALSE) index=(1:nrow(Arctic)) #area prior to 1987 is about 0.555 greater than reported due to changed hole temp=(Arctic$year==2008) plot(Arctic$jday[temp],Arctic$extent[temp],type="l",xlab="",ylab="million km2",ylim=c(4,16)) lines(Arctic$jday[temp],Arctic$area[temp],lty=3,col=2) #1990,2007 and 2008 temp=(Arctic$year==2008) plot(Arctic$jday[temp],Arctic$extent[temp],type="l",xlab="",ylab="million km2",ylim=c(0,17),xlim=c(1,366)) lines(Arctic$jday[temp],Arctic$area[temp],lty=3,col=1) for(year in 2000:2007){ temp=(Arctic$year==year) lines(Arctic$jday[temp],Arctic$extent[temp],lty=1,col=2+(2007-year)) lines(Arctic$jday[temp],Arctic$area[temp],lty=3,col=2+(2007-year)) } N0=max((1:nrow(Arctic))[!is.na(Arctic$hole)]);N0 points(Arctic$jday[N0],Arctic$area[N0],pch=19,cex=.6) points(Arctic$jday[N0],Arctic$extent[N0],pch=19,cex=.6) temp=(daily$year==2008) lines(daily$dd[temp],daily$ice[temp],lty=2) temp=(Arctic$year==2008) plot(Arctic$jday[temp],Arctic$extent[temp],type="l",xlab="",ylab="million km2",ylim=c(0,17),xlim=c(120,260)) lines(Arctic$jday[temp],Arctic$area[temp],lty=3,col=1) for(year in 1980:2007){ temp=(Arctic$year==year) lines(Arctic$jday[temp],Arctic$extent[temp],lty=1,col=2+(2007-year)) lines(Arctic$jday[temp],Arctic$area[temp],lty=3,col=2+(2007-year)) } ##COMPAR Days temp=(Arctic$jday==203) Arctic[temp,] tapply(Arctic$area,Arctic$year,min,na.rm=TRUE) N=length(1979:2008) year0=1979;temp=(Arctic$year==year0)&!is.na(Arctic$extent) plot(Arctic$jday[temp],Arctic$extent[temp],type="l",ylim=c(3,20),col=tim.colors(34)[1]) for (year0 in 1980:2008) { temp=(Arctic$year==year0)&!is.na(Arctic$extent) lines(Arctic$jday[temp],Arctic$extent[temp],col=tim.colors(34)[year0-1978]) } temp=(Arctic$year==2008)&!is.na(Arctic$extent) lines(Arctic$jday[temp],Arctic$extent[temp],col=1,lwd=3) #AREA N=length(1979:2008) year0=1979;temp=(Arctic$year==year0)&!is.na(Arctic$area) plot(Arctic$jday[temp],Arctic$area[temp],type="l",ylim=c(3,20),col=tim.colors(34)[1]) for (year0 in 1980:2008) { temp=(Arctic$year==year0)&!is.na(Arctic$area) lines(Arctic$jday[temp],Arctic$area[temp],col=tim.colors(34)[year0-1978]) } temp=(Arctic$year==2008)&!is.na(Arctic$area) lines(Arctic$jday[temp],Arctic$area[temp],col=1,lwd=3) ##HOLE temp=!is.na(Arctic$hole) plot(Arctic$julian[temp],Arctic$hole[temp])