Monday, September 26, 2011

Lindzen-Choi 2011: paper, data, scripts

Steve McIntyre was initially slightly dissatisfied with the access to the data and code used by Richard Lindzen and Yong-Sang Choi in their 2011 paper; see TRF for texts on improvements relatively to LC09 and the publication in the Asia-Pacific journal.

However, Steve just announced via Climate Audit that he was given all the files by Richard and his colleague and you may download the data files as well as the scripts in IDL. Here is my concise, compact, and clear zipped copy of the folder:
Lindzen Choi 2011: paper in PDF, data as DAT, scripts PRO in IDL

Unzipped folder (at Steve's server)
If you download the 767 kB ZIP file from the first link above and unpack it, you will see the full Lindzen Choi 2011 paper, scripts, and data. This is pretty much the format in which all papers in (not only) climate science should be accessible, I think.

You got an "A" in transparency from me, Dick. And Yong-Sang, too. ;-)

The paper is 236-Lindzen-Choi-2011-from-lindzen-mit-website.pdf.

The readme file, both in TXT and PDF form, explains the content of the four subdirectories, and here is my simple description:
  • the "scripts" folder contains all the four IDL scripts plus a JPG image with some graphs (sample feedback estimates for LW): "main" runs the remaining three, "obs" loads observations, "amip" loads models, "lc11" computes regression etc.
  • the "obs" folder contains all the four observational files: tropical sea temperature, 36-day averages ("oisst" in the name), monthly averages with "amipsst" in the name, "erbe" 36-day averaged mean flux, and "ceres" monthly averaged mean flux
  • each of the two "amip" folders, "olr" and "osr" (long-wave vs short-wave radiation), contains 11 files; each of them corresponds to one of the 11 climate models and predicts fluxes to be compared with "amipsst" from "obs"
All the data files have a very simple format with lots of lines such as
1 198501 26.767
2 198502 26.999
The obscure IDL language is easily readable (and you may download its free version) but I would still prefer if someone had a Mathematica version of this. A short time ago, I wanted to spend hours by translating the scripts to Mathematica (which would be straightforward) but it seems as a waste of time.

If I will be doing something with the files, it will be my calculations based on my algorithms. If you can, you're invited to make your checks and your own version of Lindzen Choi 2011. Recall that the paper's main object to calculate is the response of the energy fluxes to the changes of the sea surface temperature in the tropics, perhaps with some delay (–3...+3 months) added.

No comments:

Post a Comment