attach("COmonthlyMet.rda") library( fields) png("fig3.png",width=480, height=360) quilt.plot( CO.loc,CO.tmax.MAM[103,] ) US( add=T) title( "Recorded MAM max temperatures (1997)") dev.off() png("fig4.png",width=480, height=360) matplot( CO.elev, cbind( CO.tmax.MAM[103,], CO.tmin.MAM[103,]), pch="o", type="p", col=c("red", "blue"), xlab="Elevation (m)", ylab="Temperature (C)") title("Recorded MAM max (red) and min (blue) temperatures 1997") dev.off() # time series plot of Boulder, CO July and Decembertemperatures png("fig5.png",width=720, height=240) # check that # CO.names[15] is BOULDER matplot( 1895:1997, CO.tmax[,c(7,12),15], type="l", lwd=2, lty=1, col=c("red", "green"), xlab="Years", ylab="Temperature (C)") title("Recorded July and December max temperatures") dev.off()