#1871 - 2016 June - Sep
rain=scan("https://iridl.ldeo.columbia.edu/SOURCES/.IITM/.All_India/.v1871-2016/.Rainfall/.PCPN/T/(Jun-Sep)/seasonalAverage/4/mul/data.ch")


### NINO12, NINO3, NINO4, WPI, TROPGRAD and IOD - Jun-Sep


nino3 = scan("http://iridl.ldeo.columbia.edu/SOURCES/.Indices/.nino/.EXTENDED/.NINO3/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/data.ch")

nino34 = scan("https://iridl.ldeo.columbia.edu/SOURCES/.Indices/.nino/.EXTENDED/.NINO34/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/data.ch")

nino4 = scan("http://iridl.ldeo.columbia.edu/SOURCES/.Indices/.nino/.EXTENDED/.NINO4/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/data.ch")

nino12 = scan("http://iridl.ldeo.columbia.edu/SOURCES/.Indices/.nino/.EXTENDED/.NINO12/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/data.ch")

amodt = scan("http://iridl.ldeo.columbia.edu/SOURCES/.KAPLAN/.EXTENDED/.v3/.ssta/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/X/280/360/RANGE/Y/24/60/RANGE/%5BX/Y%5Daverage/detrend-bfl/data.ch")

iod = scan("https://iridl.ldeo.columbia.edu/SOURCES/.KAPLAN/.EXTENDED/.v3/.ssta/T/(Jun%201871)/(Dec%202016)/RANGE/T/4/boxAverage/T/12/STEP/Y/-10/10/RANGE/X/50/70/RANGE/%5BX/Y%5Daverage/SOURCES/.KAPLAN/.EXTENDED/.v2/.ssta/T/(Jun%201871)/(Nov%202016)/RANGE/T/4/boxAverage/T/12/STEP/Y/-10/0/RANGE/X/90/110/RANGE/%5BX/Y%5Daverage/sub/data.ch")

#Definition: area averaged SST: 60E-170E, 15S-15N 
#Hoerling et al., 2010
wpi = scan("https://iridl.ldeo.columbia.edu/SOURCES/.KAPLAN/.EXTENDED/.v3/.ssta_c9120/T/(Jun%201871)/(Sep%202016)/RANGE/T/4/boxAverage/T/12/STEP/Y/-15/15/RANGE/X/60/170/RANGE/%5BX/Y%5Daverage/data.ch")

## Trop gradient nino12 - wpi
tropgrad = nino12 - wp

## Predictor matrix
Xpred = cbind(nino12, nino3, nino4, wpi, tropgrad, iod)

##########################################################################