I went through the CRU code below, and while I try to stay out of politics, this appears to be a programming issue with Politics embedded in it.<\/div>\n

I would hazzard a guess that this started as a cover-up from issues arising from bad programming choices, and unvetted data being used for a while…. \u00c2\u00a0And then when someone noticed “it was too late to fix” without embarassing people… \u00c2\u00a0Probably the wrong people.<\/div>\n

But never the less, if you are curious, take a look…<\/div>\n

The original information is from\u00c2\u00a0Climategate: hide the decline \u00e2\u20ac\u201c codified<\/a><\/strong><\/div>\n
\n
\n
WUWT blogging ally\u00c2\u00a0 Ecotretas<\/strong><\/a> writes in to say that he has made a compendium of programming code segments that show comments by the programmer that suggest places where data may be corrected, modified, adjusted, or busted. Some the\u00c2\u00a0 HARRY_READ_ME comments are quite revealing. For those that don\u00e2\u20ac\u2122t understand computer programming, don\u00e2\u20ac\u2122t fret, the comments by the programmer tell the story quite well even if the code itself makes no sense to you.<\/p>\n
To say that the CRU code might be \u00e2\u20ac\u0153buggy\u00e2\u20ac\u009d would be\u00e2\u20ac\u00a6well I\u00e2\u20ac\u2122ll just let CRU\u00e2\u20ac\u2122s programmer tell you in his own words.<\/p><\/blockquote>\n
\n
FOIAdocumentsosborn-tree6mannoldprogmaps12.pro<\/span>FOIAdocumentsosborn-tree6mannoldprogmaps15.pro<\/span>FOIAdocumentsosborn-tree6mannoldprogmaps24.pro<\/span>; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions \n; of growing season temperatures.\u00c2\u00a0Uses \"corrected\" MXD - but shouldn't usually \n; plot past 1960 because these will be artificially adjusted to look closer to \n; the real temperatures. <\/span><\/strong><\/code><\/li>\n
FOIAdocumentsharris-treerecon_esper.pro<\/span>; Computes regressions on full, high and low pass Esper et al. (2002) series, \n; anomalies against full NH temperatures and other series. \n; CALIBRATES IT AGAINST THE LAND-ONLY TEMPERATURES NORTH OF 20 N \n; \n; Specify period over which to compute the regressions (stop in 1960 to avoid \n; the decline<\/strong><\/code><\/li>\n
FOIAdocumentsharris-treecalibrate_nhrecon.pro<\/span>; \n; Specify period over which to compute the regressions (stop in 1960 to avoid<\/strong> \n;\u00c2\u00a0the decline that affects tree-ring density records<\/strong>) \n;<\/code><\/li>\n
FOIAdocumentsharris-treerecon1.pro<\/span> \nFOIAdocumentsharris-treerecon2.pro<\/span>FOIAdocumentsharris-treerecon_jones.pro<\/span>; \n; Specify period over which to compute the regressions (stop in 1940 to avoid<\/strong> \n;\u00c2\u00a0the decline<\/strong> \n;<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>17. Inserted debug statements into anomdtb.f90, discovered that \na\u00c2\u00a0sum-of-squared variable is becoming very, very negative<\/strong>! Key \noutput from the debug statements: \n(..) \nforrtl: error (75): floating point exception \nIOT trap (core dumped) \n..so the data value is unbfeasibly large,\u00c2\u00a0but why does the \nsum-of-squares parameter OpTotSq go negative?!!<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>22. Right, time to stop pussyfooting around the niceties of Tim's labyrinthine software \nsuites - let's have a go at producing CRU TS 3.0! since failing to do that will be the \ndefinitive failure of the entire project..<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>getting seriously fed up with the state of the Australian data. so many new stations have been \nintroduced, so many false references.. so many changes that aren't documented.<\/strong> Every time a \ncloud forms I'm presented with a bewildering selection of similar-sounding sites, some with \nreferences, some with WMO codes, and some with both. And if I look up the station metadata with \none of the local references, chances are the WMO code will be wrong (another station will have \nit) and the lat\/lon will be wrong too.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>I am very sorry to report that the rest of the databases seem to be in nearly as poor a state as \nAustralia was.\u00c2\u00a0There are hundreds if not thousands of pairs of dummy stations<\/strong>, one with no WMO \nand one with, usually overlapping and with the same station name and very similar coordinates. I \nknow it could be old and new stations, but why such large overlaps if that's the case? Aarrggghhh! \nThere truly is no end in sight.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>28. With huge reluctance, I have dived into 'anomdtb' -\u00c2\u00a0and already I have \nthat familiar Twilight Zone sensation.<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>Wrote 'makedtr.for' to tackle the thorny problem of the tmin and tmax databases not \nbeing kept in step. Sounds familiar, if worrying.\u00c2\u00a0am I the first person to attempt \nto get the CRU databases in working order?!!<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>Well, dtr2cld is not the world's most complicated program. Wheras cloudreg is, and\u00c2\u00a0I \nimmediately found a mistake!<\/strong> Scanning forward to 1951 was done with a loop that,\u00c2\u00a0for \ncompletely unfathomable reasons, didn't include months! So we read 50 grids instead \nof 600!!!<\/strong> That may have had something to do with it. I also noticed, as I was correcting \nTHAT, that I reopened the DTR and CLD data files when I should have been opening the \nbloody station files!!<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>Back to the gridding. I am seriously worried that our flagship gridded data product is produced by \nDelaunay triangulation - apparently linear as well. As far as I can see,\u00c2\u00a0this renders the station \ncounts totally meaningless<\/strong>. It also means that we cannot say exactly how the gridded data is arrived \nat from a statistical perspective - since we're using an off-the-shelf product that isn't documented \nsufficiently to say that. Why this wasn't coded up in Fortran I don't know - time pressures perhaps? \nWas too much effort expended on homogenisation, that there wasn't enough time to write a gridding \nprocedure?\u00c2\u00a0Of course, it's too late for me to fix it too.<\/strong> Meh.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>Here, the expected 1990-2003 period is MISSING - so the correlations aren't so hot! Yet \nthe WMO codes and station names \/locations are identical (or close). What the hell is \nsupposed to happen here?\u00c2\u00a0Oh yeah - there is no 'supposed', I can make it up. So I have :-)<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>Well,\u00c2\u00a0it's been a real day of revelations, never mind the week<\/strong>. This morning I \ndiscovered that proper angular weighted interpolation was coded into the IDL \nroutine, but that its use was discouraged because it was slow! Aaarrrgghh. \nThere is even an option to tri-grid at 0.1 degree resolution and then 'rebin' \nto 720x360 - also deprecated! And now, just before midnight (so it counts!), \nhaving gone back to the tmin\/tmax work,\u00c2\u00a0I've found that most if not all of the \nAustralian bulletin stations have been unceremoniously dumped into the files \nwithout the briefest check for existing stations.<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>As we can see, even I'm cocking it up!<\/strong> Though recoverably. DTR, TMN and TMX need to be written as (i7.7).\/code><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a>OH FUCK THIS.<\/strong> It's Sunday evening, I've worked all weekend, and just when I thought it was done I'm \nhitting yet\u00c2\u00a0another problem that's based on the hopeless state of our databases.<\/strong> There is no uniform \ndata integrity, it's just a catalogue of issues that continues to grow as they're found.<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannmxdgrid2ascii.pro<\/span>printf,1,\u00e2\u20ac\u2122Osborn et al. (2004) gridded reconstruction of warm-season\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122(April-September) temperature anomalies (from the 1961-1990 mean).\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122Reconstruction is based on tree-ring density records.\u00e2\u20ac\u2122 \nprintf,1 \nprintf,1,\u00e2\u20ac\u2122NOTE:\u00c2\u00a0recent decline in tree-ring density has been ARTIFICIALLY<\/strong>\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122REMOVED<\/strong> to facilitate calibration. THEREFORE, post-1960 values\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122will be much closer to observed temperatures then they should be,\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122which will incorrectly imply the reconstruction is more skilful\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122than it actually is. See Osborn et al. (2004).\u00e2\u20ac\u2122<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modesdata4sweden.pro<\/span> \nFOIAdocumentsosborn-tree6summer_modesdata4sweden.proprintf,1,'IMPORTANT NOTE:' \nprintf,1,'The data after 1960 should not be used. The tree-ring density' \nprintf,1,'records tend to show a decline after 1960 relative to the summer' \nprintf,1,'temperature in many high-latitude locations. In this data set' \nprintf,1,'this \"decline\" has been artificially removed in an ad-hoc way, and' \nprintf,1,'this means that\u00c2\u00a0data after 1960 no longer represent tree-ring<\/strong> \nprintf,1,'density variations, but have been modified to look more like the<\/strong> \nprintf,1,'observed temperatures.<\/strong>'<\/code><\/li>\n
FOIAdocumentsosborn-tree6combined_wavelet_col.pro<\/span>; \n; Remove missing data from start & end (end in 1960 due to decline<\/strong>) \n; \nkl=where((yrmxd ge 1402) and (yrmxd le 1960),n) \nsst=prednh(kl)<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannmxd_pcr_localtemp.pro<\/span>; Tries to reconstruct Apr-Sep temperatures, on a box-by-box basis, from the \n; EOFs of the MXD data set. This is PCR, although PCs are used as predictors \n; but not as predictands. This PCR-infilling must be done for a number of \n; periods, with different EOFs for each period (due to different spatial \n; coverage). *BUT* don\u00e2\u20ac\u2122t do special PCR for the modern period (post-1976), \n; since they won\u00e2\u20ac\u2122t be used\u00c2\u00a0due to the decline\/correction problem<\/strong>. \n; Certain boxes that appear to reconstruct well are \u00e2\u20ac\u0153manually\u00e2\u20ac\u009d removed because \n; they are isolated and away from any trees.<\/code><\/li>\n
FOIAdocumentsosborn-tree6briffa_sep98_d.pro<\/span>;mknormal,yyy,timey,refperiod=[1881,1940] \n; \n;\u00c2\u00a0Apply a VERY ARTIFICAL correction for decline!!<\/strong> \n; \nyrloc=[1400,findgen(19)*5.+1904] \nvaladj=[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,$ \n2.6,2.6,2.6]*0.75 ; fudge factor \n(...) \n; \n;\u00c2\u00a0APPLY ARTIFICIAL CORRECTION<\/strong> \n; \nyearlyadj=interpol(valadj,yrloc,x) \ndensall=densall+yearlyadj<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modespl_decline.pro<\/span>; \n; Plots density \u00e2\u20ac\u02dcdecline\u00e2\u20ac\u2122 as a time series of the difference between \n; temperature and density averaged over the region north of 50N, \n; and an associated pattern in the difference field. \n; The difference data set is computed using only boxes and years with \n; both temperature and density in them \u00e2\u20ac\u201c i.e., the grid changes in time. \n; The pattern is computed by correlating and regressing the *filtered* \n; time series against the unfiltered (or filtered) difference data set. \n; \n;*** MUST ALTER FUNCT_DECLINE.PRO TO MATCH THE COORDINATES OF THE<\/strong> \n; START OF THE DECLINE *** ALTER THIS EVERY TIME YOU CHANGE ANYTHING ***<\/strong><\/code><\/li>\n
FOIAdocumentsosborn-tree6mannoldprogmaps12.pro<\/span>; \n; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions \n; of growing season temperatures.\u00c2\u00a0Uses \u00e2\u20ac\u0153corrected\u00e2\u20ac\u009d MXD<\/strong> \u00e2\u20ac\u201c but shouldn\u00e2\u20ac\u2122t usually \n; plot past 1960\u00c2\u00a0because these will be artificially adjusted to look closer to<\/strong> \n;\u00c2\u00a0the real temperatures.<\/strong> \n;<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannoldprogcalibrate_correctmxd.pro<\/span>; We have previously (calibrate_mxd.pro) calibrated the high-pass filtered \n; MXD over 1911-1990, applied the calibration to unfiltered MXD data (which \n; gives a zero mean over 1881-1960) after extending the calibration to boxes \n; without temperature data (pl_calibmxd1.pro). We have identified and \n;\u00c2\u00a0artificially removed (i.e. corrected) the decline<\/strong> in this calibrated \n; data set. We now recalibrate this corrected calibrated dataset against \n; the unfiltered 1911-1990 temperature data, and apply the same calibration \n; to the corrected and uncorrected calibrated MXD data.<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modescalibrate_correctmxd.pro<\/span>; No need to verify the correct and uncorrected versions, since these \n; should be identical prior to 1920 or 1930 or\u00c2\u00a0whenever the decline \n; was corrected onwards from.<\/strong><\/code><\/li>\n
FOIAdocumentsosborn-tree5densplus188119602netcdf.pro<\/span>;\u00c2\u00a0we know the file starts at yr 440, but we want nothing till 1400,<\/strong> so we \n; can skill lines (1400-440)\/10 + 1 header line \n; we now want all lines (10 yr per line) from 1400 to 1980, which is \n; (1980-1400)\/10 + 1 lines \n(...) \n;\u00c2\u00a0we know the file starts at yr 1070, but we want nothing till 1400,<\/strong> so we \n; can skill lines (1400-1070)\/10 + 1 header line \n; we now want all lines (10 yr per line) from 1400 to 1991, which is \n; (1990-1400)\/10 + 1 lines (since 1991 is on line beginning 1990)<\/code><\/li>\n<\/ul>\n
\n
\n
FOIAdocumentsosborn-tree6mannoldprogmaps12.pro<\/span>
\nFOIAdocumentsosborn-tree6mannoldprogmaps15.pro<\/span>
\nFOIAdocumentsosborn-tree6mannoldprogmaps24.pro<\/span>
\n; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions \n; of growing season temperatures.Uses \"corrected\" MXD - but shouldn't usually \n; plot past 1960 because these will be artificially adjusted to look closer to \n; the real temperatures.<\/strong><\/code><\/li>\n
FOIAdocumentsharris-treerecon_esper.pro<\/span> \n; Computes regressions on full, high and low pass Esper et al. (2002) series, \n; anomalies against full NH temperatures and other series. \n; CALIBRATES IT AGAINST THE LAND-ONLY TEMPERATURES NORTH OF 20 N \n; \n; Specify period over which to compute the regressions (stop in 1960 to avoid \n; the decline<\/strong><\/code><\/li>\n
FOIAdocumentsharris-treecalibrate_nhrecon.pro<\/span> \n; \n; Specify period over which to compute the regressions (stop in 1960 to avoid<\/strong> \n;the decline that affects tree-ring density records<\/strong>) \n;<\/code><\/li>\n
FOIAdocumentsharris-treerecon1.pro<\/span> \nFOIAdocumentsharris-treerecon2.pro<\/span>FOIAdocumentsharris-treerecon_jones.pro<\/span> \n; \n; Specify period over which to compute the regressions (stop in 1940 to avoid<\/strong> \n;the decline<\/strong> \n;<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \n17. Inserted debug statements into anomdtb.f90, discovered that \nasum-of-squared variable is becoming very, very negative<\/strong>! Key \noutput from the debug statements: \n(..) \nforrtl: error (75): floating point exception \nIOT trap (core dumped) \n..so the data value is unbfeasibly large,but why does the \nsum-of-squares parameter OpTotSq go negative?!!<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \n22. Right, time to stop pussyfooting around the niceties of Tim's labyrinthine software \nsuites - let's have a go at producing CRU TS 3.0! since failing to do that will be the \ndefinitive failure of the entire project..<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \ngetting seriously fed up with the state of the Australian data. so many new stations have been \nintroduced, so many false references.. so many changes that aren't documented.<\/strong>Every time a \ncloud forms I'm presented with a bewildering selection of similar-sounding sites, some with \nreferences, some with WMO codes, and some with both. And if I look up the station metadata with \none of the local references, chances are the WMO code will be wrong (another station will have \nit) and the lat\/lon will be wrong too.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nI am very sorry to report that the rest of the databases seem to be in nearly as poor a state as \nAustralia was.There are hundreds if not thousands of pairs of dummy stations<\/strong>, one with no WMO \nand one with, usually overlapping and with the same station name and very similar coordinates. I \nknow it could be old and new stations, but why such large overlaps if that's the case? Aarrggghhh! \nThere truly is no end in sight.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \n28. With huge reluctance, I have dived into 'anomdtb' -and already I have \nthat familiar Twilight Zone sensation.<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nWrote 'makedtr.for' to tackle the thorny problem of the tmin and tmax databases not \nbeing kept in step. Sounds familiar, if worrying.am I the first person to attempt \nto get the CRU databases in working order?!!<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nWell, dtr2cld is not the world's most complicated program. Wheras cloudreg is, andI \nimmediately found a mistake!<\/strong>Scanning forward to 1951 was done with a loop that,for \ncompletely unfathomable reasons, didn't include months! So we read 50 grids instead \nof 600!!!<\/strong>That may have had something to do with it. I also noticed, as I was correcting \nTHAT, that I reopened the DTR and CLD data files when I should have been opening the \nbloody station files!!<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nBack to the gridding. I am seriously worried that our flagship gridded data product is produced by \nDelaunay triangulation - apparently linear as well. As far as I can see,this renders the station \ncounts totally meaningless<\/strong>. It also means that we cannot say exactly how the gridded data is arrived \nat from a statistical perspective - since we're using an off-the-shelf product that isn't documented \nsufficiently to say that. Why this wasn't coded up in Fortran I don't know - time pressures perhaps? \nWas too much effort expended on homogenisation, that there wasn't enough time to write a gridding \nprocedure?Of course, it's too late for me to fix it too.<\/strong>Meh.<\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nHere, the expected 1990-2003 period is MISSING - so the correlations aren't so hot! Yet \nthe WMO codes and station names \/locations are identical (or close). What the hell is \nsupposed to happen here?Oh yeah - there is no 'supposed', I can make it up. So I have :-)<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nWell,it's been a real day of revelations, never mind the week<\/strong>. This morning I \ndiscovered that proper angular weighted interpolation was coded into the IDL \nroutine, but that its use was discouraged because it was slow! Aaarrrgghh. \nThere is even an option to tri-grid at 0.1 degree resolution and then 'rebin' \nto 720x360 - also deprecated! And now, just before midnight (so it counts!), \nhaving gone back to the tmin\/tmax work,I've found that most if not all of the \nAustralian bulletin stations have been unceremoniously dumped into the files \nwithout the briefest check for existing stations.<\/strong><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nAs we can see, even I'm cocking it up!<\/strong>Though recoverably. DTR, TMN and TMX need to be written as (i7.7).\/code><\/code><\/li>\n
FOIAdocumentsHARRY_READ_ME.txt<\/a> \nOH FUCK THIS.<\/strong>It's Sunday evening, I've worked all weekend, and just when I thought it was done I'm \nhitting yetanother problem that's based on the hopeless state of our databases.<\/strong>There is no uniform \ndata integrity, it's just a catalogue of issues that continues to grow as they're found.<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannmxdgrid2ascii.pro<\/span> \nprintf,1,\u00e2\u20ac\u2122Osborn et al. (2004) gridded reconstruction of warm-season\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122(April-September) temperature anomalies (from the 1961-1990 mean).\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122Reconstruction is based on tree-ring density records.\u00e2\u20ac\u2122 \nprintf,1 \nprintf,1,\u00e2\u20ac\u2122NOTE:recent decline in tree-ring density has been ARTIFICIALLY<\/strong>\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122REMOVED<\/strong>to facilitate calibration. THEREFORE, post-1960 values\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122will be much closer to observed temperatures then they should be,\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122which will incorrectly imply the reconstruction is more skilful\u00e2\u20ac\u2122 \nprintf,1,\u00e2\u20ac\u2122than it actually is. See Osborn et al. (2004).\u00e2\u20ac\u2122<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modesdata4sweden.pro<\/span> \nFOIAdocumentsosborn-tree6summer_modesdata4sweden.pro \nprintf,1,'IMPORTANT NOTE:' \nprintf,1,'The data after 1960 should not be used. The tree-ring density' \nprintf,1,'records tend to show a decline after 1960 relative to the summer' \nprintf,1,'temperature in many high-latitude locations. In this data set' \nprintf,1,'this \"decline\" has been artificially removed in an ad-hoc way, and' \nprintf,1,'this means thatdata after 1960 no longer represent tree-ring<\/strong> \nprintf,1,'density variations, but have been modified to look more like the<\/strong> \nprintf,1,'observed temperatures.<\/strong>'<\/code><\/li>\n
FOIAdocumentsosborn-tree6combined_wavelet_col.pro<\/span> \n; \n; Remove missing data from start & end (end in 1960 due to decline<\/strong>) \n; \nkl=where((yrmxd ge 1402) and (yrmxd le 1960),n) \nsst=prednh(kl)<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannmxd_pcr_localtemp.pro<\/span> \n; Tries to reconstruct Apr-Sep temperatures, on a box-by-box basis, from the \n; EOFs of the MXD data set. This is PCR, although PCs are used as predictors \n; but not as predictands. This PCR-infilling must be done for a number of \n; periods, with different EOFs for each period (due to different spatial \n; coverage). *BUT* don\u00e2\u20ac\u2122t do special PCR for the modern period (post-1976), \n; since they won\u00e2\u20ac\u2122t be useddue to the decline\/correction problem<\/strong>. \n; Certain boxes that appear to reconstruct well are \u00e2\u20ac\u0153manually\u00e2\u20ac\u009d removed because \n; they are isolated and away from any trees.<\/code><\/li>\n
FOIAdocumentsosborn-tree6briffa_sep98_d.pro<\/span>;mknormal,yyy,timey,refperiod=[1881,1940] \n; \n;Apply a VERY ARTIFICAL correction for decline!!<\/strong> \n; \nyrloc=[1400,findgen(19)*5.+1904] \nvaladj=[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,$ \n2.6,2.6,2.6]*0.75 ; fudge factor \n(...) \n; \n;APPLY ARTIFICIAL CORRECTION<\/strong> \n; \nyearlyadj=interpol(valadj,yrloc,x) \ndensall=densall+yearlyadj<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modespl_decline.pro<\/span> \n; \n; Plots density \u00e2\u20ac\u02dcdecline\u00e2\u20ac\u2122 as a time series of the difference between \n; temperature and density averaged over the region north of 50N, \n; and an associated pattern in the difference field. \n; The difference data set is computed using only boxes and years with \n; both temperature and density in them \u00e2\u20ac\u201c i.e., the grid changes in time. \n; The pattern is computed by correlating and regressing the *filtered* \n; time series against the unfiltered (or filtered) difference data set. \n; \n;*** MUST ALTER FUNCT_DECLINE.PRO TO MATCH THE COORDINATES OF THE<\/strong> \n;START OF THE DECLINE *** ALTER THIS EVERY TIME YOU CHANGE ANYTHING ***<\/strong><\/code><\/li>\n
FOIAdocumentsosborn-tree6mannoldprogmaps12.pro<\/span> \n; \n; Plots 24 yearly maps of calibrated (PCR-infilled or not) MXD reconstructions \n; of growing season temperatures.Uses \u00e2\u20ac\u0153corrected\u00e2\u20ac\u009d MXD<\/strong>\u00e2\u20ac\u201c but shouldn\u00e2\u20ac\u2122t usually \n; plot past 1960because these will be artificially adjusted to look closer to<\/strong> \n;the real temperatures.<\/strong> \n;<\/code><\/li>\n
FOIAdocumentsosborn-tree6mannoldprogcalibrate_correctmxd.pro<\/span> \n; We have previously (calibrate_mxd.pro) calibrated the high-pass filtered \n; MXD over 1911-1990, applied the calibration to unfiltered MXD data (which \n; gives a zero mean over 1881-1960) after extending the calibration to boxes \n; without temperature data (pl_calibmxd1.pro). We have identified and \n;artificially removed (i.e. corrected) the decline<\/strong>in this calibrated \n; data set. We now recalibrate this corrected calibrated dataset against \n; the unfiltered 1911-1990 temperature data, and apply the same calibration \n; to the corrected and uncorrected calibrated MXD data.<\/code><\/li>\n
FOIAdocumentsosborn-tree6summer_modescalibrate_correctmxd.pro<\/span> \n; No need to verify the correct and uncorrected versions, since these \n; should be identical prior to 1920 or 1930 orwhenever the decline \n; was corrected onwards from.<\/strong><\/code><\/li>\n
FOIAdocumentsosborn-tree5densplus188119602netcdf.pro<\/span> \n;we know the file starts at yr 440, but we want nothing till 1400,<\/strong>so we \n; can skill lines (1400-440)\/10 + 1 header line \n; we now want all lines (10 yr per line) from 1400 to 1980, which is \n; (1980-1400)\/10 + 1 lines \n(...) \n;we know the file starts at yr 1070, but we want nothing till 1400,<\/strong>so we \n; can skill lines (1400-1070)\/10 + 1 header line \n; we now want all lines (10 yr per line) from 1400 to 1991, which is \n; (1990-1400)\/10 + 1 lines (since 1991 is on line beginning 19<\/code><\/li>\n<\/ul>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"I went through the CRU code below, and while I try to stay out of politics, this appears to be a programming issue with Politics embedded in it. I would hazzard a guess that this started as a cover-up from issues arising from bad programming choices, and unvetted data being used for a while…. \u00c2\u00a0And […]<\/a><\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0},"categories":[54,47,19,48],"tags":[91],"_links":{"self":[{"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/posts\/2039"}],"collection":[{"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/comments?post=2039"}],"version-history":[{"count":0,"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/posts\/2039\/revisions"}],"wp:attachment":[{"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/media?parent=2039"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/categories?post=2039"},{"taxonomy":"post_tag","embeddable":true,"href":"http:\/\/www.schollnick.net\/wordpress\/wp-json\/wp\/v2\/tags?post=2039"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}