ccc-gistemp

[pdjakow]

Polecam.

Chyba każdy wie, co to jest GISTEMP. A jeśli nie, to w skrócie - jest to oprogramowanie służące do wyciągania anomalii globalnych w oparciu o szereg danych wejściowych (np. GHCN, SST). Wiemy też, że wokół tego typu serii danych narosło ostatnio wiele kontrowersji - głównie chodzi tu o HadCRUT z CRU.

Do rzeczy. Oryginalny GISS obliczany jest przez NASA. Jednak kod źródłowy programu został opublikowany i nie wiele czasu upłynęło, a pojawił się "port" oryginalnego GISSa - CCC-GISTEMP. Przepisany z fortrana na pythona. Projekt ma na celu udostępnienie tego typu rzeczy szerszemu ogółowi, nadanie im przejrzystości (stąd CCC - Clear Climate Code).

Strona projektu dostępna jest tutaj. Najnowszą wersją jest 0.3.0, jednak w najbliższych dniach pojawi się wersja 0.4.0.
Źródła (python) pobrać można stąd.

Użycie jest bardzo proste - wystarczy uruchomić plik tools/run.py i czekać na wynik (oczywiście trzeba wpierw zainstalować pythona). Nie radzę jednak uruchamiać na początku miesiąca, gdyż wyskoczą bzdury dla miesiąca ostatniego - jest to spowodowane faktem, że baza GHCN najczęściej jest już uaktualniona, natomiast reszta jeszcze nie. Ale można to obejść:

• Ściągamy ze strony http://http://ftp.emc.ncep.noaa.gov/cmb/sst/oimonth_v2 plik oiv2mon.month.gz, gdzie month to ostatni miesiąc (teraz to 201001)
• Wsadzamy go do katalogu input (jeśli katalogu niema, tworzymy go)
• Ściągamy bazę GHCN - plik v2.mean.Z
• Pobieramy oryginalnego GISTEMP: http://http://data.giss.nasa.gov/gistemp/sources/GISTEMP_sources.tar.gz
• Wyciągamy z niego plik oisstv2_mod4.clim i wrzucamy do naszego katalogu input
• Przetwarzamy tylko krok 4 CCC-GISTEMP:  tools/run.py -s 4
• Uruchamiamy CCC od początku: tools/run.py

I tyle - możemy wyliczyć średnie odchylenia globalne GISTEMP o niemal 2 tygodnie wcześniej, niż będą dostępne oficjalnie.