% Script: load CRUTEM4, CO2, emi
% Data obtainable independently from the three indicated files
% Origin changed to 1971 values for CO2 and TSI

Y = (1851:2019)';

% 1.  CLIMATE DATASET
% Global CRUTEM4 for 1851-2019 from  Hadley Climate center
% https://crudata.uea.ac.uk/cru/data/temperature/CRUTEM4-gl.dat
% Annual data taken from rightmost column (column 14) less 1st row (1850)
CRUTEM4 = [
-0.004
-0.137
-0.326
-0.125
-0.522
-0.587
-0.580
-0.555
-0.193
-0.649
-0.516
-0.756
-0.275
-0.825
-0.453
-0.440
-0.603
-0.424
-0.393
-0.549
-0.659
-0.355
-0.384
-0.483
-0.832
-0.497
-0.252
-0.023
-0.553
-0.363
-0.440
-0.336
-0.617
-0.724
-0.693
-0.532
-0.578
-0.594
-0.276
-0.446
-0.536
-0.633
-0.708
-0.484
-0.574
-0.395
-0.310
-0.396
-0.329
-0.194
-0.157
-0.398
-0.432
-0.533
-0.405
-0.184
-0.606
-0.457
-0.423
-0.333
-0.415
-0.434
-0.306
-0.063
-0.060
-0.346
-0.652
-0.417
-0.275
-0.262
-0.099
-0.248
-0.263
-0.314
-0.212
-0.064
-0.223
-0.146
-0.469
-0.108
-0.072
-0.076
-0.317
-0.025
-0.163
-0.117
-0.032
0.152
-0.002
-0.041
-0.029
-0.046
0.010
0.107
-0.119
-0.078
-0.005
0.000
-0.108
-0.301
-0.107
-0.063
0.129
-0.145
-0.192
-0.426
-0.078
0.041
0.019
-0.068
0.039
0.027
0.024
-0.283
-0.233
-0.132
-0.090
-0.207
-0.114
-0.054
-0.202
-0.225
0.138
-0.278
-0.065
-0.354
0.082
-0.047
0.039
0.125
0.300
-0.020
0.317
-0.047
-0.042
0.106
0.264
0.361
0.234
0.490
0.380
0.094
0.177
0.331
0.570
0.228
0.490
0.841
0.556
0.474
0.668
0.762
0.758
0.660
0.863
0.805
0.895
0.682
0.731
0.914
0.707
0.761
0.826
0.866
1.151
1.292
1.120
0.941
1.158
];

% 2. CO2 DATASET
% Scripps Atmospheric CO2 Data
% http://scrippsco2.ucsd.edu/assets/data/atmospheric/merged_ice_core_mlo_spo/spline_merged_ice_core_yearly.csv
% Data taken from the second column for the years 1851-2017
CO2 = [
287.190
287.263
287.272
287.222
287.128
287.002
286.853
286.694
286.534
286.386
286.259
286.163
286.108
286.105
286.161
286.278
286.454
286.686
286.954
287.239
287.529
287.815
288.087
288.331
288.534
288.698
288.836
288.958
289.075
289.199
289.339
289.508
289.715
289.972
290.289
290.664
291.089
291.553
292.046
292.558
293.077
293.591
294.085
294.548
294.964
295.312
295.575
295.767
295.909
296.026
296.139
296.276
296.458
296.697
296.999
297.351
297.738
298.144
298.552
298.952
299.353
299.769
300.207
300.660
301.126
301.595
302.058
302.506
302.930
303.325
303.699
304.060
304.417
304.777
305.135
305.491
305.841
306.185
306.522
306.853
307.185
307.526
307.882
308.257
308.651
309.052
309.448
309.821
310.152
310.422
310.622
310.762
310.854
310.916
310.965
311.022
311.099
311.208
311.357
311.543
311.764
312.025
312.338
312.709
313.143
313.632
314.171
314.757
315.382
316.040
316.722
317.426
318.150
318.900
319.683
320.507
321.376
322.290
323.250
324.256
325.305
326.397
327.532
328.709
329.929
331.195
332.513
333.880
335.291
336.735
338.204
339.691
341.194
342.711
344.236
345.765
347.293
348.811
350.307
351.772
353.212
354.638
356.071
357.535
359.048
360.616
362.242
363.924
365.657
367.431
369.246
371.102
372.997
374.924
376.874
378.845
380.834
382.846
384.889
386.974
389.109
391.302
393.555
395.861
398.212
400.595
402.997
405.404
408.793
];

emi = [
54
57
59
69
71
76
77
78
83
91
95
97
104
112
119
122
130
135
142
147
156
173
184
174
188
191
194
196
210
236
243
256
272
275
277
281
295
327
327
356
372
374
370
383
406
419
440
465
507
534
552
566
617
624
663
707
784
750
785
819
836
879
943
850
838
901
955
936
806
932
803
845
970
963
975
983
1062
1065
1145
1053
940
847
893
973
1027
1130
1209
1142
1192
1299
1334
1342
1391
1383
1160
1238
1392
1469
1419
1630
1767
1795
1841
1865
2042
2177
2270
2330
2454
2569
2580
2686
2833
2995
3130
3288
3393
3566
3780
4053
4208
4376
4614
4623
4596
4864
5016
5074
5357
5301
5138
5094
5075
5258
5417
5583
5725
5936
6066
6074
6142
6078
6070
6174
6305
6448
6556
6576
6561
6733
6893
6994
7376
7743
8042
8336
8503
8776
8697
9128
9503
9673
9773
9855
10060
10160
10260
10360
10460
];
