Blog

Blog

20 de mai. de 2012

Climategate 2: Código Usado no Modelo Causa Aumento Artificial da Temperatura



Eu ontem escrevi sobre o que já é considerado por muitos a maior fraude científica do milênio : Climategate: entitulado na cola do famoso Watergate, o escândalo em que emails e documentos hackeados de um dos mais reconhecidos centro de aquecimento global do planeta e que mostram a face suja e podre da pseudo-ciência do aquecimento global. Vou em breve publicar a tradução de mais alguns trechos quentes.

O que me levou a escrever este post, no entanto, foi o código de software que faz parte do modelo que prevê o aumento catastrófico da temperatura e que é a base dos relatórios do IPCC e de toda a histeria em torno do aquecimento global.

Primeiramente, o código é um lixo, como nós programadores costumamos dizer. Uma vergonha para qualquer programador e o qual não deveria NUNCA ter parado em mãos de terceiros. Felizmente para a humanidade foi exatamente o que aconteceu.

Ao longo do código o(s) programadore(s) lamentam a falta de dados precisos, o precário estado do banco de dados do centro, e como foi preciso ajustar diversos algorítmos para que mostrasse os resultados "corretos".

Muitas pessoas dizem que os trechos dos emails foram tirados fora de contexto anos após terem sido escritos. Comentários de código, por sua vez, tem o objetivo de documenta-lo de forma que possa ser mantido anos após ter sido escrito.

Segue abaixo alguns trechos do código:


;
; Isso funciona com REMTS SENDO UM ARRAY bi-dimensional (Nseries, ntime) de múltiplas TimeSeries
, Cuja influência será removida. Infelizmente, o IDL5.4 p_correlate
; Falha com SÉRIES > 1 mantida constante, por isso tenho de remover a sua
; INFLUÊNCIA DE AMBOS INDTS E Depts Usando Regressão Linear Múltipla e entao usar a funcao correlata nos resíduos.

;
pro maps12, yrstart, doinfill = doinfill
;
; Plota 24 mapas anuais de reconstruções MXD calibradas (PCR-infilled ou não)
; de temperaturas crescentes. Usa MXD "corrigido" - mas geralmente não deve
; plotar antes de 1960, porque estes últimos serão artificialmente ajustados para parecer mais perto de temperaturas reais



; Calcula regressões em passa-alta e passa-baixa séries Esper et al. (2002)
; ,anomalias contra temperaturas NH completa e outras séries.
; CALIBRA CONTRA TEMPERATURAS DE TERRA AO NORTE DA LATITUDE 20N
; Especificar o período acima do qual se de deseja calcular as regressões
;(parar em 1960 para evitar O declínio)

http://www.anenglishmanscastle.com/HARRY_READ_ME.txt

Estou ficando seriamente cansado com o estado dos dados da australia. tantas novas estações foram introduzidas, tantas referências falsas .. tantas mudanças que não são documentadas. Cada vez que uma nuvem se forma eu sou apresentado com uma seleção assombrosa de sites semelhantes de sondagem, alguns com referências, alguns com os códigos da OMM (Organizacao metereológica mundial), e alguns com ambos.

28. Com grande relutância, eu mergulhei na 'anomdtb "- e já estou com aquela sensação familiar de Twilight Zone.

eu sou a primeira pessoa a tentar por os dados CRU (Climate Research Unit) em ordem?!

Bem, dtr2cld não é o programa mais complicados do mundo. Considerando cloudreg, que é, e eu imediatamente encontrei um erro! A leitura para a frente a 1951 foi feito com um loop que, por razões completamente insondáveis, não incluem meses! Assim, lemos 50 grids ao invés
de 600!

...Por isso que isto não foi codificado em Fortran não sei - por pressões de tempo talvez? Foi muito esforço dispendido na homogeneização, que não houve tempo suficiente para escrever um procedimento de gridding? Naturalmente, é tarde demais para eu corrigi-lo também. Meh.

Aqui, o esperado período 1990-2003 esta FALTANDO- assim as correlações não são tão quentes! No entanto, os códigos da OMM e os nomes das estações / posições são idênticos (ou quase). Que diabos devia acontecer aqui? Ah, sim - nao devia acontecer, entao eu posso inventar. E então foi o que eu fiz :-)

Como podemos ver, até eu estou fazendo merda! Embora possa recuperar. DTR, TMN e TMX precisa ser escrito como (i7.7). / Code>

AHH FODA-SE. É noite de domingo, eu trabalhei todos os fins de semana, e só quando eu pensei que havia acabado eu recebo ainda outro problema que é baseado no estado desesperador de nossos bancos de dados. Não há integridade de dados uniforme, é apenas um catálogo de problemas que continua a crescer a medida que eles são encontrados.


;mknormal,yyy,timey,refperiod=[1881,1940]
;
; Aplicar uma correção MUITO ARTIFICIAL para declínio!
yrloc=[1400,findgen(19)*5.+1904]
valadj=[0.,0.,0.,0.,0.,-0.1,-0.25,-0.3,0.,-0.1,0.3,0.8,1.2,1.7,2.5,2.6,2.6,$
2.6,2.6,2.6]*0.75 ; fator lorota

FOIA\documents\osborn-tree6\summer_modes\pl_decline.pro
; *** Deve alterar FUNCT_DECLINE.PRO PARA IGUALAR COM AS COORDENADAS DO INÍCIO DO DECLÍNIO. ALTERAR ISTO CADA VEZ QUE VOCE MUDAR QUALQUER COISA ***

O código original pode ser encontrado aqui.

Neste post, um blogger tenta rodar o código, que parece necessitar várias intervenções de operação.

O BRASIL IRÁ TORNAR LEI O TRATADO DE COPENHAGUE MESMO QUE OS OUTROS PAÍSES NAO O FAÇAM. CONTATE O SEU SENADOR/DEPUTADO E MOSTRE A ELE ESTAS NOTÍCIAS.
Abaixo é um vídeo de um programador sobre o código:


Update:
Este blog faz uma análise sobre o código Fortran divulgado:
http://www.di2.nu/blog.htm