#% Global Average Temperature Anomaly with Sea Ice Temperature Inferred from Air Temperatures
#% Global Average Temperature Anomaly with Sea Ice Temperature Inferred from Water Temperatures
#Estimated Jan 1951-Dec 1980 monthly absolute temperature 

get_berk=function(dset="water") {
loc="http://berkeleyearth.lbl.gov/auto/Global/Land_and_Ocean_complete.txt"
    #75] "%                  Monthly          Annual          Five-year        Ten-year        Twenty-year" 
  # [76] "% Year, Month,  Anomaly, Unc.,   Anomaly, Unc.,   Anomaly, Unc.,   Anomaly, Unc.,   Anomaly, Unc."
 
download.file(loc,dest)
fred=readLines(dest)
m=grep("Year",fred) #76 2086 

if(dset=="air") {
 tem=read.table(dest,skip=m[1]+1,nrow=m[2]-m[1]-7,colClasses=rep("numeric",12))
 berk=ts(tem[,3],start=1850,freq=12)
 return(berk)
 }

if(dset=="water") {
 tem=read.table(dest,skip=m[2]+1,colClasses=rep("numeric",12))
 berk=ts(tem[,3],start=1850,freq=12)
 return(berk)
 }
}

# berk=get_berk(dset="water")
# tsp(berk)#
#[1] 1850.000 2016.917   12.000