Skip to content

Instantly share code, notes, and snippets.

@mikelove

mikelove/elnino.R

Created August 20, 2015 20:15
Show Gist options
  • Save mikelove/0490ef13647e73581382 to your computer and use it in GitHub Desktop.
Save mikelove/0490ef13647e73581382 to your computer and use it in GitHub Desktop.
Download ZIP
elnino
Raw
elnino.R
library(RColorBrewer)
palette(brewer.pal(4,"Set1"))
x <- read.table("snow.txt", header=TRUE)
x$EN <- factor(x$EN, c("n","m","s","vs"))
levels(x$EN) <- c("no/weak", "medium", "strong", "very strong")
plot(x$Tot, col=x$EN, pch=20, cex=3,
xaxt="n", xlab="",
ylab="total snowfall", main="Boston snowfall")
abline(h=1:5 * 20, col=rgb(0,0,0,.5))
axis(1, 1:nrow(x), as.character(x$Season), las=2)
legend("topleft", col=1:4, levels(x$EN), pch=20,
title="El Nino", bg="white", inset=.05)
Raw
elnino.txt
Season Jul Aug Sep Oct Nov Dec Jan Feb Mar Apr May Jun Tot EN
1950-51 0.0 0.0 0.0 0.0 T 2.7 13.9 9.2 3.9 T 0.0 0.0 29.7 n
1951-52 0.0 T 0.0 0.0 T 8.4 10.7 19.1 1.4 0.0 0.0 T 39.6 m
1952-53 0.0 0.0 0.0 T T 2.4 13.6 11.4 0.2 2.2 T 0.0 29.8 n
1953-54 0.0 0.0 0.0 0.0 T T 19.2 1.9 0.4 2.1 0.0 0.0 23.6 n
1954-55 T T 0.0 0.0 T 10.3 0.9 6.5 7.0 0.4 0.0 0.0 25.1 n
1955-56 T 0.0 0.0 0.0 2.5 1.8 7.7 14.5 31.2 3.2 0.0 0.0 60.9 n
1956-57 0.0 0.0 0.0 0.0 0.6 15.4 20.6 2.8 11.5 1.1 0.0 0.0 52.0 n
1957-58 0.0 0.0 0.0 T 0.1 T 6.6 23.9 12.0 2.1 T 0.0 44.7 s
1958-59 0.0 0.0 0.0 0.0 0.1 4.6 4.1 10.7 14.6 T 0.0 0.0 34.1 n
1959-60 0.0 0.0 0.0 0.0 0.6 5.5 10.2 2.3 22.3 T 0.0 0.0 40.9 n
1960-61 0.0 0.0 0.0 T T 16.9 18.7 14.9 9.0 2.0 0.0 0.0 61.5 n
1961-62 0.0 0.0 0.0 T 0.9 11.4 2.5 28.7 1.1 0.1 0.0 0.0 44.7 n
1962-63 0.0 0.0 0.0 T 0.9 5.3 6.5 4.6 13.6 T 0.0 0.0 30.9 n
1963-64 0.0 0.0 0.0 T 0.0 17.7 14.4 23.2 7.7 T 0.0 0.0 63.0 m
1964-65 0.0 0.0 0.0 T T 12.2 22.2 4.7 9.7 1.6 0.0 0.0 50.4 n
1965-66 0.0 0.0 0.0 0.0 T 2.3 26.4 12.1 3.3 T T 0.0 44.1 s
1966-67 0.0 0.0 0.0 0.0 0.0 9.9 0.5 23.5 22.9 3.3 T 0.0 60.1 n
1967-68 0.0 0.0 0.0 0.0 2.2 14.7 17.7 3.4 6.8 0.0 0.0 0.0 44.8 n
1968-69 0.0 0.0 0.0 0.0 0.4 5.1 0.9 41.3 6.1 T 0.0 0.0 53.8 m
1969-70 0.0 0.0 0.0 T T 12.6 7.4 10.5 18.2 0.1 0.0 0.0 48.8 n
1970-71 0.0 0.0 0.0 T T 27.9 12.0 8.1 7.4 1.9 0.0 0.0 57.3 n
1971-72 0.0 0.0 0.0 0.0 2.8 7.9 7.8 16.5 12.1 0.4 0.0 0.0 47.5 n
1972-73 0.0 0.0 0.0 T 0.6 3.3 3.6 2.5 0.3 T 0.0 0.0 10.3 s
1973-74 0.0 0.0 0.0 0.0 0.0 T 16.0 17.8 0.1 3.0 0.0 0.0 36.9 n
1974-75 0.0 0.0 0.0 0.0 2.0 3.6 2.2 17.0 1.8 1.0 0.0 0.0 27.6 n
1975-76 0.0 0.0 0.0 T 0.1 19.3 15.0 1.4 10.8 T 0.0 0.0 46.6 n
1976-77 0.0 0.0 0.0 0.0 1.0 17.2 23.2 5.9 10.7 T 0.5 0.0 58.5 n
1977-78 0.0 0.0 0.0 0.0 0.7 5.2 35.9 27.2 16.1 T 0.0 0.0 85.1 n
1978-79 0.0 0.0 0.0 0.0 4.2 5.8 10.5 6.6 T 0.4 0.0 0.0 27.5 n
1979-80 0.0 0.0 0.0 0.2 T 2.0 0.4 6.5 3.6 T 0.0 0.0 12.7 n
1980-81 0.0 0.0 0.0 0.0 2.4 5.6 11.9 1.9 0.5 0.0 0.0 0.0 22.3 n
1981-82 0.0 0.0 0.0 0.0 T 17.6 18.0 7.6 5.3 13.3 0.0 0.0 61.8 n
1982-83 0.0 0.0 0.0 0.0 T 5.5 4.7 22.3 0.2 T 0.0 0.0 32.7 vs
1983-84 0.0 0.0 0.0 0.0 T 2.6 21.1 0.3 19.0 T 0.0 0.0 43.0 n
1984-85 0.0 0.0 0.0 0.0 T 3.7 7.0 10.2 3.7 2.0 0.0 0.0 26.6 n
1985-86 0.0 0.0 0.0 0.0 3.0 1.3 0.8 10.4 2.6 T 0.0 0.0 18.1 n
1986-87 0.0 0.0 0.0 0.0 3.5 3.4 24.3 3.7 3.5 4.1 0.0 0.0 42.5 m
1987-88 0.0 0.0 0.0 0.0 9.0 7.5 17.0 14.1 5.0 T 0.0 0.0 52.6 s
1988-89 0.0 0.0 0.0 T 0.0 3.7 1.5 6.7 3.2 0.4 0.0 0.0 15.5 n
1989-90 0.0 0.0 0.0 0.0 4.5 6.2 7.0 16.9 4.1 0.5 0.0 0.0 39.2 n
1990-91 0.0 0.0 0.0 0.0 T 1.2 11.7 2.8 3.4 0.0 0.0 0.0 19.1 n
1991-92 0.0 0.0 0.0 0.0 T 5.8 0.4 4.0 10.8 1.0 0.0 0.0 22.0 m
1992-93 0.0 0.0 0.0 0.0 0.6 9.7 12.9 19.6 38.9 2.2 0.0 0.0 83.9 n
1993-94 0.0 0.0 0.0 0.0 T 11.6 33.7 36.2 14.8 0.0 0.0 0.0 96.3 n
1994-95 0.0 T 0.0 0.0 0.1 1.5 4.4 8.5 0.4 T 0.0 0.0 14.9 m
1995-96 0.0 0.0 0.0 0.0 4.1 24.1 39.8 15.5 16.8 7.3 0.0 0.0 107.6 n
1996-97 0.0 0.0 0.0 0.0 1.8 5.0 9.7 4.8 8.2 22.4 0.0 0.0 51.9 n
1997-98 0.0 0.0 0.0 0.0 3.9 7.8 10.0 0.3 3.6 0.0 0.0 0.0 25.6 vs
1998-99 0.0 0.0 0.0 0.0 0.0 0.8 16.8 7.4 11.4 0.0 0.0 0.0 36.4 n
1999-00 0.0 0.0 0.0 0.0 0.0 0.0 13.7 9.2 2.0 T 0.0 0.0 24.9 n
2000-01 0.0 0.0 0.0 T T 4.5 12.4 9.8 19.2 T 0.0 0.0 45.9 n
2001-02 0.0 0.0 0.0 0.0 0.0 5.0 7.9 0.5 1.4 0.3 0.0 T 15.1 n
2002-03 0.0 0.0 0.0 0.0 3.6 9.6 4.2 39.9 8.1 2.3 0.0 0.0 67.7 m
2003-04 0.0 0.0 0.0 0.0 0.0 17.2 4.9 2.4 10.6 0.0 0.0 0.0 35.1 n
2004-05 0.0 0.0 0.0 0.0 3.9 7.2 43.3 17.7 14.5 0.0 0.0 0.0 86.6 n
2005-06 0.0 0.0 0.0 1.1 0.0 10.7 8.1 20.0 0.0 0.0 0.0 0.0 39.9 n
2006-07 0.0 0.0 0.0 0.0 0.0 0.8 1.0 4.6 10.2 0.5 0.0 0.0 17.1 n
2007-08 0.0 0.0 0.0 0.0 0.0 26.9 8.3 15.0 1.0 0.0 0.0 0.0 51.2 n
2008-09 0.0 0.0 0.0 0.0 0.0 25.3 23.7 6.2 10.7 0.0 0.0 0.0 65.9 n
2009-10 0.0 0.0 0.0 0.1 0.0 15.2 13.2 7.0 0.2 0.0 0.0 0.0 35.7 m
2010-11 0.0 0.0 0.0 0.0 0.0 22.0 38.3 18.5 1.3 0.9 0.0 0.0 81.0 n
@mikelove
Copy link
Author

data provenance:
http://www.erh.noaa.gov/box/climate/bossnw.shtml
http://liveweatherblogs.com/index.php/community/groups/viewdiscussion/1457-boston-boston-ma-winter-outlook-2015-2016-strong-to-super-el-nino-impact-on-snow-temps?groupid=64

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment