#Project 1: Toy stochastic climate models

from pylab import *
import numpy
import datetime
import time
import glob, os
import math
import matplotlib.colors as colors
from scipy import stats

#-------------------------------------------------------------------------------
#----------------------------MAIN SCRIPT STARTS HERE----------------------------
#-------------------------------------------------------------------------------

#Get the time series
file = open('SSH_SST_Pacific.txt', 'r')
lines = file.readlines()
file.close()

time    = np.zeros(len(lines) - 3)
ssh     = np.zeros(len(time))
temp    = np.zeros(len(time))

for time_i in range(len(time)):
    #Read in the data
    line = lines[time_i + 3].split()

    #Save the corresponding data in the relevant array
    time[time_i]    = datetime.datetime(int(line[2]), int(line[1]), int(line[0])).toordinal()
    ssh[time_i]     = float(line[3])
    temp[time_i]    = float(line[4])
	
#-------------------------------------------------------------------------------
#----------------------------------EXERCISE 1-----------------------------------
#-------------------------------------------------------------------------------