http://beforeitsnews.com/story/0000000000000790
This is the actual source code (from the CRU file briffa_Sep98_d.pro), which is organized into numbered lines of code written in Interactive Data Language (IDL). The computer uses an interpreter to read each line, parses out the various commands and does exactly what they tell it to do.
1 ;
2 ; Apply a VERY ARTIFICAL correction for decline!!
3 ;
4 yrloc=[1400,findgen(19)*5.+1904]
5 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
6 if n_elements(yrloc) ne n_elements(valadj) then message,'Oooops!'
7
8
yearlyadj=interpol(valadj,yrloc,timey)
Going through, line by line is not difficult:
The first 3 lines are comments. The person writing the comment is explaining that he will apply something VERY ARTIFICIAL to the original temperature readings. They are not off to a good start -- if the evidence is so compelling, why apply anything artificial? Why not just show the actual readings? Because the actual readings don't behave the way you want them to, that's why.
Line 4 creates a 20 element array (a series of ordered numbers) for the variable yrloc, which represents the years in the series of data points, starting with the year 1400 and then 1904, plus another 18 years, each incremented by 5 years from the previous year in the series. This series of ordered numbers is created by the words and numbers on the right hand side of the "=" .
On the right side of the "=" , the term findgen makes the numbers floating point, i.e. have a decimal. In computer programming like this, you can't mix and match different number types (fixed = number with no decimal, e.g. 5, floating = number WITH a decimal, e.g. 5.0). So everything here is going to be floating point and in this case 5.0, which increments the series of numbers in the array.
Here is how the array gets populated with numbers. The first number is 1400, the next is 1904, those are the numbers plugged into the formula. The remaining numbers are generated as follows:
1904 + 5 = 1909
1909 + 5 = 1914
1914 + 5 = 1919...and so on, until you have done that 19 times - that's what the (19) means
This results in the following sequence of numbers for yrloc:
1400, 1904, 1909, 1914, 1919, 1924, 1929, 1935, 1939, 1944, 1949, 1954, 1959, 1964, 1969, 1975, 1979, 1984, 1989, 1994
Line 5, as it so eloquently calls itself is the "fudge factor" (the computer model author's word, not mine. These numbers are added into the model to obtain the correct result. Let's see what these "fudge factor" numbers look like on a graph. No temperature numbers, just the "fudge factor". These numbers are plugged into the variable array "valadj" as follows:
0, 0, 0, 0, 0, -0.1, -0.25, etc. I have plotted them here:
This is the actual source code (from the CRU file briffa_Sep98_d.pro), which is organized into numbered lines of code written in Interactive Data Language (IDL). The computer uses an interpreter to read each line, parses out the various commands and does exactly what they tell it to do.
1 ;
2 ; Apply a VERY ARTIFICAL correction for decline!!
3 ;
4 yrloc=[1400,findgen(19)*5.+1904]
5 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
6 if n_elements(yrloc) ne n_elements(valadj) then message,'Oooops!'
7
8
yearlyadj=interpol(valadj,yrloc,timey)
Going through, line by line is not difficult:
The first 3 lines are comments. The person writing the comment is explaining that he will apply something VERY ARTIFICIAL to the original temperature readings. They are not off to a good start -- if the evidence is so compelling, why apply anything artificial? Why not just show the actual readings? Because the actual readings don't behave the way you want them to, that's why.
Line 4 creates a 20 element array (a series of ordered numbers) for the variable yrloc, which represents the years in the series of data points, starting with the year 1400 and then 1904, plus another 18 years, each incremented by 5 years from the previous year in the series. This series of ordered numbers is created by the words and numbers on the right hand side of the "=" .
On the right side of the "=" , the term findgen makes the numbers floating point, i.e. have a decimal. In computer programming like this, you can't mix and match different number types (fixed = number with no decimal, e.g. 5, floating = number WITH a decimal, e.g. 5.0). So everything here is going to be floating point and in this case 5.0, which increments the series of numbers in the array.
Here is how the array gets populated with numbers. The first number is 1400, the next is 1904, those are the numbers plugged into the formula. The remaining numbers are generated as follows:
1904 + 5 = 1909
1909 + 5 = 1914
1914 + 5 = 1919...and so on, until you have done that 19 times - that's what the (19) means
This results in the following sequence of numbers for yrloc:
1400, 1904, 1909, 1914, 1919, 1924, 1929, 1935, 1939, 1944, 1949, 1954, 1959, 1964, 1969, 1975, 1979, 1984, 1989, 1994
Line 5, as it so eloquently calls itself is the "fudge factor" (the computer model author's word, not mine. These numbers are added into the model to obtain the correct result. Let's see what these "fudge factor" numbers look like on a graph. No temperature numbers, just the "fudge factor". These numbers are plugged into the variable array "valadj" as follows:
0, 0, 0, 0, 0, -0.1, -0.25, etc. I have plotted them here:
