It looks like you're using an Ad Blocker.
Please white-list or disable AboveTopSecret.com in your ad-blocking tool.
Thank you.
Some features of ATS will be disabled while you continue to use an ad-blocker.
People are talking about the emails being smoking guns but I find the remarks in the code and the code more of a smoking gun. The code is so hacked around to give predetermined results that it shows the bias of the coder. In other words make the code ignore inconvenient data to show what I want it to show.
function mkp2correlation,indts,depts,remts,t,filter=filter,refperiod=refperiod,$
datathresh=datathresh
;
; THIS WORKS WITH REMTS BEING A 2D ARRAY (nseries,ntime) OF MULTIPLE TIMESERIES
; WHOSE INFLUENCE IS TO BE REMOVED. UNFORTUNATELY THE IDL5.4 p_correlate
; FAILS WITH >1 SERIES TO HOLD CONSTANT, SO I HAVE TO REMOVE THEIR INFLUENCE
; FROM BOTH INDTS AND DEPTS USING MULTIPLE LINEAR REGRESSION AND THEN USE THE
; USUAL correlate FUNCTION ON THE RESIDUALS.
;
pro maps12,yrstart,doinfill=doinfill
;
; 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.
;
;
; Plots (1 at a time) 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.
From 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
;
; Computes regressions on full, high and low pass MEAN timeseries of MXD
; anomalies against full NH temperatures.
; THIS IS FOR THE AGE-BANDED (ALL BANDS) STUFF OF HARRY’S
;
; Specify period over which to compute the regressions (stop in 1940 to avoid
; the decline
;
recon_mann.pro:
; Computes regressions on full, high and low pass MEAN timeseries of MXD
; anomalies against full NH temperatures.
; THIS IS FOR THE Mann et al. reconstruction
; CALIBRATES IT AGAINST THE LAND-ONLY TEMPERATURES NORTH OF 20 N
; IN FACT, I NOW HAVE AN ANNUAL LAND-ONLY NORTH OF 20N VERSION OF MANN,
; SO I CAN CALIBRATE THIS TOO – WHICH MEANS I’m ONLY ALTERING THE SEASON
briff_sep98_e.pro:
;
; PLOTS ‘ALL’ REGION MXD timeseries from age banded and from hugershoff
; standardised datasets.
; Reads Harry’s regional timeseries and outputs the 1600-1992 portion
; with missing values set appropriately. Uses mxd, and just the
; “all band” timeseries
;****** APPLIES A VERY ARTIFICIAL CORRECTION FOR DECLINE*********
****** APPLIES A VERY ARTIFICIAL CORRECTION FOR DECLINE*********
these will be artificially adjusted
(stop in 1940 to avoid
; the decline
"I've just completed Mike's Nature [the science journal] trick of adding in the real temps to each series for the last 20 years (ie, from 1981 onwards) and from 1961 for Keith's to hide the decline."
from README_GRIDDING.TXT..
“Use dist to specify the correlation decay distance for the climate
variable being interpolated – necessary information to determine where to add dummy or synthetic data.”
The fact is that we can’t account for the lack of warming at the moment and it is a travesty that we can’t. The CERES data published in the August BAMS 09 supplement on 2008 shows there should be even more warming: but the data are surely wrong. Our observing system is inadequate.
……Phil and I have recently submitted a paper using about a dozen NH records that fit this category, and many of which are available nearly 2K back–I think that trying to adopt a timeframe of 2K, rather than the usual 1K, addresses a good earlier point that Peck made w/ regard to the memo, that it would be nice to try to “contain” the putative “MWP”, even if we don’t yet have a hemispheric mean reconstruction available that far back….
Well I have my own article on where the heck is global warming ? We are asking that here in Boulder where we have broken records the past two days for the coldest days on record. ... The fact is that we can’t account for the lack of warming at the moment and it is a travesty that we can’t. The CERES data published in the August BAMS 09 supplement on 2008 shows there should be even more warming: but the data are surely wrong. Our observing system is inadequate.
We probably need to say more about this. Land warming since 1980 has been twice the ocean warming — and skeptics might claim that this proves that urban warming is real and important.
Here are some speculations on correcting SSTs to partly explain the 1940s warming blip. If you look at the attached plot you will see that theland also shows the 1940s blip (as I’m sure you know).
So, if we could reduce the ocean blip by, say, 0.15 degC, then this would be significant for the global mean – but we’d still have to explain the land blip. I’ve chosen 0.15 here deliberately. This still leaves an ocean blip, and i think one needs to have some form of ocean blip to explain the land blip (via either some common forcing, or ocean forcing land, or vice versa, or all of these). When you look at other blips, the land blips are 1.5 to 2 times (roughly) the ocean blips—higher sensitivity plus thermal inertia effects. My 0.15 adjustment leaves things consistent with this, so you can see where I am coming from.
Removing ENSO does not affect this.
It would be good to remove at least part of the 1940s blip, but we are still left with “why the blip”.
Let me go further. If you look at NH vs SH and the aerosol effect (qualitatively or with MAGICC) then with a reduced ocean blip we get continuous warming in the SH, and a cooling in the NH—just as one would expect with mainly NH aerosols.
The other interesting thing is (as Foukal et al. note – from MAGICC) that the 1910-40 warming cannot be solar. The Sun can get at most 10% of this with Wang et al solar, less with Foukal solar. So this may well be NADW, as Sarah and I noted in 1987 (and also Schlesinger later). A reduced SST blip in the 1940s makes the 1910-40 warming larger than the SH (which it currently is not)—but not really enough.
So ... why was the SH so cold around 1910? Another SST problem? (SH/NH data also attached.)
This stuff is in a report I am writing for EPRI, so I’d appreciate any comments you (and Ben) might have.
Tom.
The Korttajarvi record was oriented in the reconstruction in the way that McIntyre said. I took a look at the original reference – the temperature proxy we looked at is x-ray density, which the author interprets to be inversely related to temperature. We had higher values as warmer in the reconstruction, so it looks to me like we got it wrong, unless we decided to reinterpret the record which I don’t remember.
As we all know, this isn't about truth at all, its about plausibly deniable accusations.
Yeah, it wasn’t so much 1998 and all that that I was concerned about, used to dealing with that, but the possibility that we might be going through a longer – 10 year – period of relatively stable temperatures beyond what you might expect from La Nina etc. Speculation, but if I see this as a possibility then others might also. Anyway, I’ll maybe cut the last few points off the filtered curve before I give the talk again as that’s trending down as a result of the end effects and the recent cold-ish years.
“This whole project is SUCH A MESS. No wonder I needed therapy!!”
— HARRY_READ_ME.txt
; calculate 1961-1990 synthetic normal from adjusted tmn
print,'Calculating synthetic frs normal'
for iy=nor1,nor2 do begin
tmpfl=strip(string(tmp_prefix,iy))
dtrfl=strip(string(dtr_prefix,iy))
rdbin,tmpgrd,tmpfl,gridsize=2.5,/quiet
rdbin,dtrgrd,dtrfl,gridsize=2.5,/quiet
tmn(nland)=(tmpgrd(nland)-(0.5*dtrgrd(nland)))/10.0
frssyn(nland)=frssyn(nland)+frscal(tmn(nland))
endfor
frssyn(nland)=frssyn(nland)/(nor2-nor1+1)
for im=0,11 do begin
temp=frssyn(*,*,im)
nfin=where(temp gt 0)
temp(nfin)=(temp(nfin)/100.0)*days(im)
frssyn(*,*,im)=temp
endfor
frssyn(nsea)=-999.9
; Calculate synthetic frs from tmin, convert to anomalies
; relative to synthetic mean frs, and apply to normal frs
print,'Calculating synthetic anomalies'
; HUGREG=Hugershoff regions, ABDREG=age-banded regions, HUGGRID=Hugershoff grid
; The calibrated (uncorrected) versions of all these data sets are used.
; However, the same adjustment is then applied to the corrected version of
; the grid Hugershoff data, so that both uncorrected and corrected versions
; are available with the appropriate low frequency variability. There is some
; ambiguity during the modern period here, however, because the corrected
; version has already been artificially adjusted to reproduce the largest
; scales of observed temperature over recent decades - so a new adjustment
; would be unwelcome. Therefore, the adjustment term is scaled back towards
; zero when being applied to the corrected data set, so that it is linearly
; interpolated from its 1950 value to zero at 1970 and kept at zero thereafter.
; Combines the directly calibrated MXD data set with the PCR-based
; reconstruction of gridded temperatures. There are various PCR models to
; use, according to period and spatial coverage of MXD data. We always
; use the later model (based on most MXD data), but we have to decide whether
; a grid box that was successfully reconstructed using an earlier subset of
; the MXD should be used throughout (or at all) if later subsets failed to
; successfully reconstruct it. **For now, I'm using them throughout.**
;
; Restore MXD gridded dataset
;
print,'Reading in MXD data'
restore,filename='calibmxd5_abdlow.idlsave'
; g,mxdyear,mxdnyr,fdcalibu,fdcalibc,mxdfd2,timey,fdseas
;
; Use the "corrected" calibrated version
;
fdcalibpcr=fdcalibc
timeyr=mxdyear
nyr=mxdnyr
;
; Now process each PCR version in turn
;
for iper = 0 , 6 do begin
;
case iper of
0: perst='14001976'
1: perst='14531976'
2: perst='15831976'
3: perst='16601976'
4: perst='16971976'
5: perst='17431976'
6: perst='18221976'
endcase
;
; Restore the next PCR-based reconstruction
;
fn='calibmxd5_abdlow_pcr'+perst+'.idlsave'
print,fn
restore,filename=fn
; Gets: mxdyear,mxdnyr,nx,ny,xlon,ylat,fdstatus,recontemp
;
; Check that the grids match
;
abserr=total(abs(g.x-xlon))
if abserr ne 0. then message,'Grids do not match!'
abserr=total(abs(g.y-ylat))
if abserr ne 0. then message,'Grids do not match!'
;
; Identify period of overlap
;
ist=where(timeyr eq mxdyear(0)) & ist=ist(0)
;
; Put in the reconstruction, replacing any data already there
;
for iyr = 0 , mxdnyr-1 do begin
fd=reform(recontemp(*,*,iyr))
kl=where(finite(fd),nkeep)
if nkeep gt 0 then begin
oldfd=fdcalibpcr(*,*,ist+iyr)
oldfd(kl)=fd(kl)
fdcalibpcr(*,*,ist+iyr)=oldfd
endif
endfor
;
endfor
;
; Finally, replace the original calibrated MXD values, since these are
; preferable to the PCR-based reconstructions
;
; Computes EOFs of infilled calibrated MXD gridded dataset.
; Can use corrected or uncorrected MXD data (i.e., corrected for the decline).
; Do not usually rotate, since this loses the common volcanic and global
; warming signal, and results in regional-mean series instead.
; Generally use the correlation matrix EOFs.
;
; 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.
; Plots a HovMueller diagram (longitude-time) of meridionally averaged
; growing season reconstructions. Uses "corrected" MXD - but shouldn't usually
; plot past 1960 because these will be artificially adjusted to look closer to
; the real temperatures.
; 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.
Originally posted by IgnoranceIsntBlisss
Perhaps you should head over to Climate Audit to put McIntyre and the others (many actual programmers and climate scientists) in their places, damn silly foolhearted feebleminded denier conspiracy kook nutters they all must be.
And then some:
"I've just completed Mike's Nature [the science journal] trick of adding in the real temps to each series for the last 20 years (ie, from 1981 onwards) and from 1961 for Keith's to hide the decline."
“This whole project is SUCH A MESS. No wonder I needed therapy!!”
— HARRY_READ_ME.txt
What do you say, mel?
it would be nice to try to “contain” the putative “MWP”, even if we don’t yet have a hemispheric mean reconstruction available that far back
Originally posted by pumpkinorange
reply to post by IgnoranceIsntBlisss
Editorial: Hiding Evidence of Global Cooling:
www.washingtontimes.com...