###################################################################################
# Date: 2014-04-30
# Time series analysis
# listening to Climate history
###################################################################################
rm(list=ls())
library(audio)

# interpolated 10-year resolution
# 65 North insolation during Jun-August

buffer<-read.table("data/timeseries_insolation_JJA_800000.asc")
time<-buffer[,1]
x<-buffer[,2]
n<-length(x)

plot(time,x,xlab="time [years]",
     ylab="Jun-Jul-Aug Insolation [W/m^2]",
     typ='l')

wait<-readline("press enter to listen to the time series")  
rate<-22050


play(x,rate=rate)

#