# GISS Temp Anomaly Trend - Lead & Sea C - 1880 - 2008
## Monthly anomaly data
link1<- "ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_and_ocean.90S.90N.df_1901-2000mean.dat"
mo_GISS_data <- read.table(link1,
             skip = 0, sep = "", dec=".",
             row.names = NULL, header = FALSE,
             as.is = T,
             colClasses = rep("numeric",3),
             comment.char = "#", na.strings = c("*", "-",-99.9, -999.9),
             col.names = c("Yr", "Mo",  "Mo_Anomaly") )
mo_GISS_ts <- ts(mo_GISS_data$Mo_Anomaly, start = c(1880,1), frequency= 12)
par(las=1); par(ps = c(8)) ;  par(mfrow=c(1,1))
plot(mo_GISS_ts, type = "h", col = "grey", xlab = "",
ylab = expression(paste("GISS Temperature Anomaly - ",degree,"C")),
        xlim=c(1880, 2009),
     main = expression(paste("GISS Land & Sea Temperature Anomalies - ",degree, "C (1901 - 2000 Baseline)")) ,
     xaxs="i", yaxs="i",
     sub="Data Source: ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/monthly.land_and_ocean.90S.90N.df_1901-2000mean.dat")
# Annual Anomaly Data
link2 <- "ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/annual.land_and_ocean.90S.90N.df_1901-2000mean.dat"
 ann_GISS_data <- read.table(link2,
             skip = 0, sep = "", dec=".",
             row.names = NULL, header = FALSE,
             as.is = T,
             colClasses = rep("numeric",2),
             comment.char = "#", na.strings = c("*", "-",-99.9, -999.9),
             col.names = c("Yr", "Ann_Anomaly") )
ann_GISS_ts <- ts(ann_GISS_data$Ann_Anomaly, start = c(1880,1), frequency= 1)
points(ann_GISS_ts, type = "s", col = "red", lwd=1.5)
points(lowess(mo_GISS_ts, f=0.15),type = "l", lwd=1.5,  col = "blue")
lines(c(1882,1885),c(0.8, 0.8), col = "grey", lwd=2)
lines(c(1882, 1885), c(0.75, 0.75), col = "red", lwd = 2)
lines(c(1882, 1885), c(0.7, 0.7), col = "blue",lwd = 2)
text( 1886, 0.8, "Monthly Anomaly", pos =4)
text( 1886, 0.75, "Annual Anomaly",pos = 4)
text( 1886, 0.7, "Lowess Smooth Fit, f= 0.15", pos = 4)

#lines(c(1998, 1998), c(-0.3, -0.03), col = "brown")
#lines(c(2008, 2008), c(-0.3, -0.03), col = "brown")
#text(1990, -0.5, "Recent period considered by \n G Will, B Carter, D Whitehouse \n Ignores long term trend ")
#arrows(1990, -.42,2003, -0.3, length = 0.0, col = "brown")