Posts Tagged ‘R code’

BOE: Rates of Change in Antarctic Ice Core Temperatures

2012 June 21 2 comments

Another line of argument for a catastrophic impact of climate change is in the rate of change of global temperatures and/or ocean acidification. In this post, I’m going to poke a bit at the rates of change in global temperature via ice core data. This step – from globally to particular Antarctic locations – is questionable and I’ll look closer to see if this is justified in a latter post. With that said, what do the ice cores hold ….?

Read more…


Colorado River Basin Storage: Lake Mead and Lake Powell

2012 June 6 Comments off

Some of my first comments on WUWT were in regards to Lake Powell water levels. Even then I understood that to correctly document water levels on the Colorado River Basin, you had to look, as a minimum, at both Mead and Powell. There was a new river operations approach begun in 2005, formalized in 2007, which has had the net effect of slowly increasing the water levels in Lake Powell, but at the expense of Lake Mead. Indeed, there was little change in the river system as a whole from 2005-2010, with combined water storage for the two reservoirs bouncing around 50%. This winter was the first year in the last eight or so with a sustained improvement in storage – although whether that can be maintained is still in question with forecasts of widespread low water flows.

Read more…

wxcalc: building a library of weather related functions

2012 June 4 1 comment

As we saw in the previous thread, the activity of eolian sands can be characterized by the ratio of precipitation to potential evapotranspiration. In order to calculate the latter, I need access to some weather related functions such as the derivation of vapor pressure from dew point temperature. Searching for the equations I needed, I stumbled over the “Weather Calculator” site (by Tim Brice and Todd Hall) hosted on a El Paso, Texas NOAA NWS site. While the calculator isn’t of much use for me, the equations behind it were available in a set of pdf files.

Read more…

Steig Siple dD: 11th – 20th Century with Lowess

2012 May 13 Comments off

Steve McIntyre has a post about Steig’s Siple dD series in which he posts a chart of the 20th century to undermine Schmidt’s claim that there are exceptionally high values in late 20th Century. And, indeed, the 20th century appears to be relatively trendless when eyeballed on that chart. McIntyre included his code.

Read more…

NCDC DS-9640: CONUS Temperature Anomalies and a couple of ax-grinders

2012 May 12 11 comments

For some reason, two posts at The Blackboard, It’s “Fancy,” Sort of … (Shollenberger) and “To get what he wanted”: Upturned end points. (Lucia), seem to be having difficulty understanding the mechanics of another post at “Open Mind”, In the Classroom (Tamino). But there is nothing unusual or difficult about the methods Tamino used to create the charts which have generated so much smoke and apparent frustration at The Blackboard – resulting in an outbreak of mcintyretude: scorn, derision, insults, and the questioning of motives. Since this is mostly a quick walk through of some code to clear the smoke, I will leave the charts generated to post at the end.

Read more…

SST Regions: Tropical West Pacific, 23N-23S, 120E-180E

2011 March 9 Comments off
Next up is the Tropical West Pacific as defined by a bounding box 23N-23S, 120E-180E. There are some issues cropping HadSST2. I noticed it in the NW Pacific as well, but thought it was related to the fact that what should have been the upper row was all ‘NA.’ Gonna have to crack the raster ‘crop’ code open. However, I doubt that a row more or less on the boundaries will materially affect the results.

Read more…

SST Regions: NorthEast Pacific 23N-65N, 120E-180E

2011 March 8 Comments off
Next up is the NorthEast Pacific as defined by a bounding box 23N-65N, 100W-180W. Over the weekend, there were some modifications to my first presentation. Inverted the x-axis on the Fourier chart. Formatted the freq text on Fourier chart. Created a preset color scale for the maps. Colors and cuts from HadSST2 map. I also saved off the ERSST annual raster so that I don’t have to recalculate it for each run.

Read more…