There's been a lot of analysis of the Hadley/CRU emails disclosed in the data leak. There's been less - but very interesting indeed - analysis of the source code disclosed in the leak.
That got me wondering: what happens when you run the code? I mean, the whole IPCC Thermageddon story is based on this code, right? Running the code might be interesting.
So I ran it. "Interesting"? Boy, Howdy.
The first thing that became obvious is that this is not at all a manual process. There is a lot of operator intervention, with the program asking you for all sorts of input. So much so, in fact, that it looks like a programmer at CRU took a shot at trying to automate the process:
I'm not interested in grants, and it looks like the editor of GRL has already been properly dealt with, so I selected the first option and forged ahead.
Unfortunately, the code looks to be very fragile, with poor error handling. I'm in the middle of stack dumps, but thought that these results so far were interesting enough to post. Here's the (admittedly limited) error report that I have so far.
UPDATE 8 December 2009 15:42: Climate Model Clippy is now available as a coffee mug.
UPDATE 2 December 2009 16:23: Snark aside, I've distilled the (scientific) reasons to be an AGW skeptic here.
UPDATE 24 November 2009 18:05: Welcome visitors from A Libertarian's Perspective! Lots of Climate Change posts here. And more Climate Change Snark here.
UPDATE 24 November 2009 18:16: Bienvenue visiteurs de Changement Climatique!
UPDATE 24 November 2009 18:19: Welcome visitors from The Common Room, 180grader.dk, Whirlpool.net.au, OFN Community, and Ann Coulter!
UPDATE 24 November 2009 18:31: Welcome visitors from Bishop's Hill!
8 comments:
*APPLAUSE*
Bravo, encore, more !
Have twittered.
Lol! That's awesome.
I'm here all week. Try the veal.
;-)
Borepatch, thanks for making me spray coffee all over my flat panel.
That Clippy panel has got to be one of the funniest things I've seen in a long time!
Brilliant!!!
Excellent use of the Office Irritant.
Linked and distributing.
How about this piece of code, have you ran it yet?
function sdtest(n1,n2,iy,im,data,imiss,nyrs)
integer data(1701:2000,12)
sum=0.0
sum2=0.0
n=0
do jy=n1,n2
if(jy.ne.iy.and.data(jy,im).ne.imiss)then
sum=sum+0.1*real(data(jy,im))
sum2=sum2+(0.1*real(data(jy,im)))**2
n=n+1
endif
enddo
if(n.gt.nyrs)then
sum=sum/n
sum2=sum2/n
diff=sum2-sum**2
if(diff.lt.0)diff=0.0
sdtest=sqrt(diff)
else
sdtest=imiss
endif
return
end
Great clippy, Bore Patch. Update the info on my sideblog at americandigest.org.
The Clippy is one of my favorites. I laughed for an hour when I first saw it, and I still smile widely now. Thanks.
Post a Comment