Wednesday, December 23, 2009

Phew! Got the limited coverage error sorted out

I had been worried about reproducing the error range shown in Brohan et al. section 6.1 because I was getting much larger error bars than in the paper.

But I wrote some tests and dug into the code and discovered a problem where my regridder/subsampler was ignoring some data because of a difference in the way longitude was reported between my gridding program and the NetCDF file from the reanalysis data.

Bottom line: I now get values much, much closer to the published information and all my worries about the temperature for 2009 appearing inside the error bars for the 1860s are gone.

Phew.

Here's the new chart showing the smoothed trend line with 95% confidence bars.


And here for comparison is the chart from the paper (mine is northern hemisphere and this is global but you can clearly see the similarity now).


Now, I can retire from amateur climatology... I've worked through Brohan et al. using the Met Office data, found a few problems, but mostly reproduced their results.

Labels:

If you enjoyed this blog post, you might enjoy my travel book for people interested in science and technology: The Geek Atlas. Signed copies of The Geek Atlas are available.

<$BlogCommentBody$>

<$BlogCommentDateTime$> <$BlogCommentDeleteIcon$>

Post a Comment

Links to this post:

<$BlogBacklinkControl$> <$BlogBacklinkTitle$> <$BlogBacklinkDeleteIcon$>
<$BlogBacklinkSnippet$>
Create a Link

<< Home