library("Cairo") #glsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land.90S.90N.df_1901-2000mean.dat" #nhsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land.00N.90N.df_1901-2000mean.dat" #shsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land.90S.00N.df_1901-2000mean.dat" #glsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.ocean.90S.90N.df_1901-2000mean.dat" #nhsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.ocean.00N.90N.df_1901-2000mean.dat" #shsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.ocean.90S.00N.df_1901-2000mean.dat" #glsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_ocean.90S.90N.df_1901-2000mean.dat" #nhsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_ocean.00N.90N.df_1901-2000mean.dat" #shsrc="ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_ocean.90S.00N.df_1901-2000mean.dat" # land and ocean nhsrc="../../data/ncdc/data/monthly.land_ocean.00N.90N.df_1901-2000mean.dat" shsrc="../../data/ncdc/data/monthly.land_ocean.90S.00N.df_1901-2000mean.dat" glsrc="../../data/ncdc/data/monthly.land_ocean.90S.90N.df_1901-2000mean.dat" # land only #nhsrc="../../data/ncdc/data/monthly.land.00N.90N.df_1901-2000mean.dat" #shsrc="../../data/ncdc/data/monthly.land.90S.00N.df_1901-2000mean.dat" #glsrc="../../data/ncdc/data/monthly.land.90S.90N.df_1901-2000mean.dat" # ocean only # nhsrc="../../data/ncdc/data/monthly.ocean.00N.90N.df_1901-2000mean.dat" # shsrc="../../data/ncdc/data/monthly.ocean.90S.00N.df_1901-2000mean.dat" # glsrc="../../data/ncdc/data/monthly.ocean.90S.90N.df_1901-2000mean.dat" fgl<-read.table(glsrc,fill=TRUE,na.strings="-999.0000") fnh<-read.table(nhsrc,fill=TRUE,na.strings="-999.0000") fsh<-read.table(shsrc,fill=TRUE,na.strings="-999.0000") gl <- aggregate(fgl[,3],list(fgl[,1]), FUN=mean, na.rm=T) nh <- aggregate(fnh[,3],list(fnh[,1]), FUN=mean, na.rm=T) sh <- aggregate(fsh[,3],list(fsh[,1]), FUN=mean, na.rm=T) # fourier just from 1880-2010 spgl <- spectrum(gl[,2], log="no") spnh <- spectrum(nh[,2], log="no") spsh <- spectrum(sh[,2], log="no") 1/spgl$freq[spgl$spec>0.05] 1/spnh$freq[spnh$spec>0.05] 1/spsh$freq[spsh$spec>0.05] # Land and Ocean # 1/spgl$freq[spgl$spec>0.05] # [1] 135.00000 67.50000 45.00000 33.75000 19.28571 # 1/spnh$freq[spnh$spec>0.05] # [1] 135.000000 67.500000 45.000000 33.750000 19.285714 9.000000 5.192308 # 1/spsh$freq[spsh$spec>0.05] # [1] 135.0 67.5 45.0 # Land Only # 1/spgl$freq[spgl$spec>0.05] # [1] 135.000000 67.500000 45.000000 33.750000 19.285714 5.192308 # [7] 4.821429 3.970588 3.552632 2.045455 2.014925 # 1/spnh$freq[spnh$spec>0.05] # [1] 135.000000 67.500000 45.000000 33.750000 19.285714 9.642857 # [7] 9.000000 8.437500 7.500000 5.400000 5.192308 4.821429 # [13] 4.655172 4.218750 3.970588 3.552632 2.500000 2.454545 # [19] 2.045455 2.014925 # 1/spsh$freq[spsh$spec>0.05] # [1] 135.000000 67.500000 33.750000 22.500000 19.285714 15.000000 # [7] 13.500000 6.428571 4.218750 3.648649 3.552632 3.214286 # Ocean Only # 1/spgl$freq[spgl$spec>0.05] # [1] 135.0 67.5 45.0 # 1/spnh$freq[spnh$spec>0.05] # [1] 67.50 45.00 33.75 9.00 # 1/spsh$freq[spsh$spec>0.05] # [1] 135 45 27