Passar dos emails à prática

Os emails já nos haviam preparado para o facto de que verdades inconvenientes estariam a ser apagadas. Vejamos o que disse Phil Jones, na primeira pessoa:

From: Phil Jones, 2/2/2005 09:41 AM
The two MMs have been after the CRU station data for years. If they ever hear there is a Freedom of Information Act now in the UK, I think I'll delete the file rather than send to anyone.

Ora, saiu ontem no Sunday Times a admissão pela parte da Universidade de East Anglia, onde trabalha o tal Phil Jones, de que os dados de base que servem para o cálculo das temperaturas mundiais, foram efectivamente enviados para o lixo. Tal significa que ficaram apenas os dados tratados (leia-se manipulados), tornando impossível que outros cientistas possam confirmar a suposta subida de temperaturas nos últimos 150 anos.

Aparentemente, a desculpa é a de que os dados originais, em papel e banda magnética, foram deitados fora quando a CRU mudou para instalações novas! Dadas as manipulações já detectadas, e as ameaças explícitas nos emails, como podemos confiar neste gang do Climategate?

www.timesonline.co.uk/tol/news/environment/article6936328.ece

Comecem os cortes de estrada

Ainda não chegamos ao Inverno, e já começa o corte de estradas neste País. Certamente, para muitos será mais uma consequência do Aquecimento Global... Fica o registo do corte na Serra da Estrela, para estragar o turismo, como é habitual, os escuteiros por lá desprevenidos, e o IP4 também cortado! Fica também o registo do Instituto de Meteorologia:

"Queda de neve acima dos 1400 metros, descendo a cota para os 800/1000 metros durante a manhã e para os 600/800 metros para o final do dia."

www.correiomanha.pt/noticia.aspx?channelid=00000021-0000-0000-0000-000000000021&contentid=7872E84E-5FD0-4397-99BD-308B7447CFDC&h=6
http://tsf.sapo.pt/PaginaInicial/Vida/Interior.aspx?content_id=1434178
www.publico.clix.pt/Local/primeira-queda-de-neve-deste-inverno-corta-o-ip4_1411925

Aquecimento Global a acelerar nos media

O Aquecimento Global está a acelerar nos media, a poucos dias do início da Cimeira de Copenhaga. Os excerptos das notícias abaixo dão conta de uma subida exponencial das temperaturas previstas para 2100, que já duplicam as piores previsões do IPCC! Tal subida é apenas correlacionável com a estupidez dos cientistas e jornalistas que as publicam:

• (24/03/2006) O modelo prevê que as temperaturas no Árctico aumentem três a cinco graus Celsius, durante o Verão, até 2100.
  
• (18/11/2009) O documento de 150 páginas conclui que é necessário examinar com mais detalhes as maneiras como o aquecimento da superfície terrestre - que se prosseguir, ou acelerar as tendências recentes, pode aumentar de quatro a seis graus até 2100
  
• (25/11/2009) De acordo com um grupo de 24 cientistas internacionais, a temperatura global deverá aumentar entre dois e sete graus até 2100.
  
• (28/11/2009) O Árctico é a região do planeta onde a temperatura está a aumentar mais rapidamente, com uma taxa de aquecimento duas vezes superior à do resto do planeta. A previsão é de um aumento de até 9ºC durante o século XXI.

Arrogância ecologista malarenta

A arrogância dos ecologistas não tem limites. Na sequência do ClimateGate, a televisão Russia Today, provavelmente a que de forma mais isenta tem transmitido o tema, fez um pequeno debate entre Piers Corbyn, dono do Weather Action, e Aleksey Kokorin, Coordenador do programa climático, da WWF Russia. O debate tem interesse na sua totalidade, mas o que interessa mais neste contexto é o que decorre a partir do minuto 9:06.

A parte interessante da discussão começa quando Piers Corbyn, aos 7:50, refere que os recursos que se estão a consumir para Copenhaga seriam mais úteis combatendo a malária. O moderador surpreendido, ataca Corbyn, referindo que a malária seria uma consequência do Aquecimento Global. Corbyn refere que isso não faz sentido, e aos 8:00 refere que toda a gente na Rússia sabe que existiu malária no passado. Aos 9.10, Aleksey refere que nunca houve malária na Rússia antes, que tem como resposta uma sonora gargalhada de Corbyn, afirmando que tal é uma mentira!

O que se segue é um Aleksey passado, a perguntar desalmadamente quando é que o Piers esteve na Rússia pela última vez, a afirmar que é ele que vive lá, a perguntar ao Piers se é um médico ou um meteorologista, ... A salganhada só terminou com o moderador a interrompê-los a ambos!

Quando vi isto, não sabia a resposta. Também não precisava de ter estado na Rússia para saber. Na Internet, em 5 minutos, encontrei todas as respostas que precisava. E a verdade é muito inconveniente para este Aleksey, com números que nos dizem ter existido mais de 30 milhões de casos de málária, na antiga União Soviética, durante a década de 20, do século passado! Será que estes ecologistas da treta não enxergam? Ou será que são simplesmente estúpidos?

www.ncbi.nlm.nih.gov/pmc/articles/PMC2537933/pdf/bullwho00505-0074.pdf
www.who.int/library/collections/historical/en/index4.html
www.nature.com/nature/journal/v139/n3505/abs/139022c0.html
www.uel.ac.uk/mosquito/issue12/russia.htm
http://query.nytimes.com/gst/abstract.html?res=9803EFDE113FE432A25753C2A9649C946095D6CF

Autópsia do Aquecimento Global

Enquanto se discute se o Aquecimento Global morreu ou não, para mim ele já está na morgue. Prontinho para a autópsia. O Aquecimento Global teve uma vida curta, tendo tido alguns irmãos mais velhos, mas nasceu verdadeiramente a 23 de Junho de 1988. O pai foi o James Hansen, que provocou o parto desligando os sistemas do ar condicionado do Congresso dos Estados Unidos. O miúdo cresceu rapidamente durante a década de 90, e antes dos 10 anos já era um superdotado, tendo terminado o primeiro ciclo com Kyoto. Como teenager, cresceu intelectualmente, tendo atingido a maioridade em 2006, com o filme "A Verdade Inconveniente". Em 2009 morreu de morte súbita, aparentemente de enfarte alarmista.

A morte ainda não foi reconhecida pelos seus fãs, que já anseiam pelo regresso do seu Elvis. Também ainda não apareceu na secção de necrologia da maioria dos media, pois estes continuam subjugados ao lápis azul da censura politicamente correcta. Quando este bebé nasceu, era assim que se vivia atrás do muro: as pessoas começavam a saber o que se passava do outro lado, mas não passava nada na imprensa oficial. Tal como quando o Titanic se afundou, a banda continua a tocar, na expectativa de que nada possa terminar a viagem com rumo para Copenhaga!

O que se seguirá só pode ter um caminho: Temperature Data Wants To Be Free! Já começaram a surgir algumas vontades nesse sentido, de que se destaca a vontade de tornar este processo open-source. O processo do open-source já provou ao mundo que é uma forma muito interessante de evolução, e todos poderiam contribuir para o esforço, mutando-se o código genético do Aquecimento Global, feito em Fortran, por formas de codificação mais modernas. Certamente, a geração do Aquecimento Global perecerá, incluindo o meio-irmão do Peer-Review, que será substituído na hierarquia por um processo mais aberto e distribuído!

Somos os maiores (outra vez)

O nosso primeiro só sabe dizer que somos os maiores. Quando a esmola é grande, o pobre desconfia sempre, e ele voltou à carga hoje: "Esta é a maior central fotovoltaica do mundo em ambiente urbano e fica como símbolo de um país que fez há alguns anos uma escolha e uma opção: a energia renovável"

José Sócrates efectuou esta afirmação na inauguração da central fotovoltaica do Mercado Abastecedor da Região de Lisboa. A dita, segundo várias fontes (eg. Wikipedia), tem uma potência de 6MW e dará para abastecer 3 mil lares / 12 mil pessoas. Custou a módica quantia de 31 milhões de euros, que os portugueses pagarão pelo deficit tarifário.

Mas uma pesquisa rápida pela Internet desmistifica rapidamente esta ideia de sermos os maiores. Aqui ao lado, em Espanha, em Saragoça, a General Motors colocou no topo da sua fábrica, 85.000 módulos solares fotovoltaicos, em 183.000 metros quadrados, com uma potência de 10MW, com capacidade de satisfazer o consumo de 4.600 lares. Tudo isto há mais de um ano!

http://diariodigital.sapo.pt/news.asp?section_id=12&id_news=422735
www.oje.pt/noticias/economia/central-fotovoltaica-no-marl-e-a-maior-do-mundo-em-ambiente-urbano-diz-socrates

Bugs do Climategate

Depois dos emails, começou a ser escrutinado o código revelado a partir do Climate Research Unit. Este código, que há muito é solicitado, ao abrigo do Freedom of Information Act, foi conhecido no âmbito do ClimateGate. Steve McIntyre há muito que pede estes dados, para replicar as conclusões do estudo, mas agora sabemos que havia muitas Verdades Inconvenientes nele escondido...

Agora que o código é conhecido, as primeiras conclusões são de horror! Sendo que a teoria do Aquecimento Global está baseada nestes dados/algoritmos, codificados em programas informáticos com a qualidade atribuível a um caloiro de um qualquer curso de Informática, não há nada como começar a dar uma vista de olhos! E não me venham com a ideia de que são expressões manhosas de uns quantos emails. Estamos a falar dos programas dos supostos modelos, e das notas dos respectivos programadores, que preveêm um futuro Mundo muito negro. Espalhem pelos V/ amigos informáticos: eles vão adorar estes bugs...

• FOIA\documents\osborn-tree6\mann\oldprog\maps12.pro
  FOIA\documents\osborn-tree6\mann\oldprog\maps15.pro
  FOIA\documents\osborn-tree6\mann\oldprog\maps24.pro
  ; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions
; of growing season temperatures. Uses "corrected" MXD - but shouldn't usually
; plot past 1960 because these will be artificially adjusted to look closer to
; the real temperatures.
  
  
• FOIA\documents\harris-tree\recon_esper.pro
  ; Computes regressions on full, high and low pass Esper et al. (2002) series,
; anomalies against full NH temperatures and other series.
; CALIBRATES IT AGAINST THE LAND-ONLY TEMPERATURES NORTH OF 20 N
;
; Specify period over which to compute the regressions (stop in 1960 to avoid
; the decline
  
  
• FOIA\documents\harris-tree\calibrate_nhrecon.pro
  ;
; Specify period over which to compute the regressions (stop in 1960 to avoid
; the decline that affects tree-ring density records)
;
  
  
• FOIA\documents\harris-tree\recon1.pro
  FOIA\documents\harris-tree\recon2.proFOIA\documents\harris-tree\recon_jones.pro
  ;
; Specify period over which to compute the regressions (stop in 1940 to avoid
; the decline
;
  
  
• FOIA\documents\HARRY_READ_ME.txt
  17. Inserted debug statements into anomdtb.f90, discovered that
a sum-of-squared variable is becoming very, very negative! Key
output from the debug statements:
(..)
forrtl: error (75): floating point exception
IOT trap (core dumped)
..so the data value is unbfeasibly large, but why does the
sum-of-squares parameter OpTotSq go negative?!!
  
  
• FOIA\documents\HARRY_READ_ME.txt
  22. Right, time to stop pussyfooting around the niceties of Tim's labyrinthine software
suites - let's have a go at producing CRU TS 3.0! since failing to do that will be the
definitive failure of the entire project..
  
  
• FOIA\documents\HARRY_READ_ME.txt
  getting seriously fed up with the state of the Australian data. so many new stations have been
introduced, so many false references.. so many changes that aren't documented. Every time a
cloud forms I'm presented with a bewildering selection of similar-sounding sites, some with
references, some with WMO codes, and some with both. And if I look up the station metadata with
one of the local references, chances are the WMO code will be wrong (another station will have
it) and the lat/lon will be wrong too.
  
  
• FOIA\documents\HARRY_READ_ME.txt
  I am very sorry to report that the rest of the databases seem to be in nearly as poor a state as
Australia was. There are hundreds if not thousands of pairs of dummy stations, one with no WMO
and one with, usually overlapping and with the same station name and very similar coordinates. I
know it could be old and new stations, but why such large overlaps if that's the case? Aarrggghhh!
There truly is no end in sight.
  
  
• FOIA\documents\HARRY_READ_ME.txt
  28. With huge reluctance, I have dived into 'anomdtb' - and already I have
that familiar Twilight Zone sensation.
  
  
• FOIA\documents\HARRY_READ_ME.txt
  Wrote 'makedtr.for' to tackle the thorny problem of the tmin and tmax databases not
being kept in step. Sounds familiar, if worrying. am I the first person to attempt
to get the CRU databases in working order?!!
  
  
• FOIA\documents\HARRY_READ_ME.txt
  Well, dtr2cld is not the world's most complicated program. Wheras cloudreg is, and I
immediately found a mistake! Scanning forward to 1951 was done with a loop that, for
completely unfathomable reasons, didn't include months! So we read 50 grids instead
of 600!!! That may have had something to do with it. I also noticed, as I was correcting
THAT, that I reopened the DTR and CLD data files when I should have been opening the
bloody station files!!
  
  
• FOIA\documents\HARRY_READ_ME.txt
  Back to the gridding. I am seriously worried that our flagship gridded data product is produced by
Delaunay triangulation - apparently linear as well. As far as I can see, this renders the station
counts totally meaningless. It also means that we cannot say exactly how the gridded data is arrived
at from a statistical perspective - since we're using an off-the-shelf product that isn't documented
sufficiently to say that. Why this wasn't coded up in Fortran I don't know - time pressures perhaps?
Was too much effort expended on homogenisation, that there wasn't enough time to write a gridding
procedure? Of course, it's too late for me to fix it too. Meh.
  
  
• FOIA\documents\HARRY_READ_ME.txt
  Here, the expected 1990-2003 period is MISSING - so the correlations aren't so hot! Yet
the WMO codes and station names /locations are identical (or close). What the hell is
supposed to happen here? Oh yeah - there is no 'supposed', I can make it up. So I have :-)
  
  
• FOIA\documents\HARRY_READ_ME.txt
  Well, it's been a real day of revelations, never mind the week. This morning I
discovered that proper angular weighted interpolation was coded into the IDL
routine, but that its use was discouraged because it was slow! Aaarrrgghh.
There is even an option to tri-grid at 0.1 degree resolution and then 'rebin'
to 720x360 - also deprecated! And now, just before midnight (so it counts!),
having gone back to the tmin/tmax work, I've found that most if not all of the
Australian bulletin stations have been unceremoniously dumped into the files
without the briefest check for existing stations.
  
  
• FOIA\documents\HARRY_READ_ME.txt
  As we can see, even I'm cocking it up! Though recoverably. DTR, TMN and TMX need to be written as (i7.7)./code>


• FOIA\documents\HARRY_READ_ME.txt
OH FUCK THIS. It's Sunday evening, I've worked all weekend, and just when I thought it was done I'm
hitting yet another problem that's based on the hopeless state of our databases. There is no uniform
data integrity, it's just a catalogue of issues that continues to grow as they're found.


• FOIA\documents\osborn-tree6\mann\mxdgrid2ascii.pro
printf,1,’Osborn et al. (2004) gridded reconstruction of warm-season’
printf,1,’(April-September) temperature anomalies (from the 1961-1990 mean).’
printf,1,’Reconstruction is based on tree-ring density records.’
printf,1
printf,1,’NOTE: recent decline in tree-ring density has been ARTIFICIALLY’
printf,1,’REMOVED to facilitate calibration. THEREFORE, post-1960 values’
printf,1,’will be much closer to observed temperatures then they should be,’
printf,1,’which will incorrectly imply the reconstruction is more skilful’
printf,1,’than it actually is. See Osborn et al. (2004).’


• FOIA\documents\osborn-tree6\summer_modes\data4sweden.pro
FOIA\documents\osborn-tree6\summer_modes\data4sweden.pro
printf,1,'IMPORTANT NOTE:'
printf,1,'The data after 1960 should not be used. The tree-ring density'
printf,1,'records tend to show a decline after 1960 relative to the summer'
printf,1,'temperature in many high-latitude locations. In this data set'
printf,1,'this "decline" has been artificially removed in an ad-hoc way, and'
printf,1,'this means that data after 1960 no longer represent tree-ring
printf,1,'density variations, but have been modified to look more like the
printf,1,'observed temperatures.'


• FOIA\documents\osborn-tree6\combined_wavelet_col.pro
;
; Remove missing data from start & end (end in 1960 due to decline)
;
kl=where((yrmxd ge 1402) and (yrmxd le 1960),n)
sst=prednh(kl)


• FOIA\documents\osborn-tree6\mann\mxd_pcr_localtemp.pro
; Tries to reconstruct Apr-Sep temperatures, on a box-by-box basis, from the
; EOFs of the MXD data set. This is PCR, although PCs are used as predictors
; but not as predictands. This PCR-infilling must be done for a number of
; periods, with different EOFs for each period (due to different spatial
; coverage). *BUT* don’t do special PCR for the modern period (post-1976),
; since they won’t be used due to the decline/correction problem.
; Certain boxes that appear to reconstruct well are “manually” removed because
; they are isolated and away from any trees.


• FOIA\documents\osborn-tree6\briffa_sep98_d.pro
;mknormal,yyy,timey,refperiod=[1881,1940]
;
; Apply a VERY ARTIFICAL correction for decline!!
;
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 ; fudge factor


• FOIA\documents\harris-tree\briffa_sep98_e.pro
;****** APPLIES A VERY ARTIFICIAL CORRECTION FOR DECLINE*********
;
(...)
mknormal,densadj,x,refperiod=[1881,1960],refmean=refmean,refsd=refsd
mknormal,densall,x,refperiod=[1881,1960],refmean=refmean,refsd=refsd
;
; APPLY ARTIFICIAL CORRECTION
;
yearlyadj=interpol(valadj,yrloc,x)
densall=densall+yearlyadj


• FOIA\documents\osborn-tree6\summer_modes\pl_decline.pro
;
; Plots density ‘decline’ as a time series of the difference between
; temperature and density averaged over the region north of 50N,
; and an associated pattern in the difference field.
; The difference data set is computed using only boxes and years with
; both temperature and density in them – i.e., the grid changes in time.
; The pattern is computed by correlating and regressing the *filtered*
; time series against the unfiltered (or filtered) difference data set.
;
;*** MUST ALTER FUNCT_DECLINE.PRO TO MATCH THE COORDINATES OF THE
; START OF THE DECLINE *** ALTER THIS EVERY TIME YOU CHANGE ANYTHING ***


• FOIA\documents\osborn-tree6\mann\oldprog\maps12.pro
;
; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions
; of growing season temperatures. Uses “corrected” MXD – but shouldn’t usually
; plot past 1960 because these will be artificially adjusted to look closer to
; the real temperatures.
;


• FOIA\documents\osborn-tree6\mann\oldprog\calibrate_correctmxd.pro
; We have previously (calibrate_mxd.pro) calibrated the high-pass filtered
; MXD over 1911-1990, applied the calibration to unfiltered MXD data (which
; gives a zero mean over 1881-1960) after extending the calibration to boxes
; without temperature data (pl_calibmxd1.pro). We have identified and
; artificially removed (i.e. corrected) the decline in this calibrated
; data set. We now recalibrate this corrected calibrated dataset against
; the unfiltered 1911-1990 temperature data, and apply the same calibration
; to the corrected and uncorrected calibrated MXD data.


• FOIA\documents\osborn-tree6\summer_modes\calibrate_correctmxd.pro
; No need to verify the correct and uncorrected versions, since these
; should be identical prior to 1920 or 1930 or whenever the decline
; was corrected onwards from.


• FOIA\documents\osborn-tree5\densplus188119602netcdf.pro
; we know the file starts at yr 440, but we want nothing till 1400, so we
; can skill lines (1400-440)/10 + 1 header line
; we now want all lines (10 yr per line) from 1400 to 1980, which is
; (1980-1400)/10 + 1 lines
(...)
; we know the file starts at yr 1070, but we want nothing till 1400, so we
; can skill lines (1400-1070)/10 + 1 header line
; we now want all lines (10 yr per line) from 1400 to 1991, which is
; (1990-1400)/10 + 1 lines (since 1991 is on line beginning 1990)


Nota: Editado para corrigir falha de localização de código, originalmente referenciando FOIA\documents\osborn-tree6\briffa_sep98_d.pro mas na verdade em FOIA\documents\harris-tree\briffa_sep98_e.pro