load 'uh5ydpm'
M = uh5ydpm;

mat1 = M (1:11,:);

years1 = mat1 (:,1);
dp1 = mat1 (:,2);

bb1 = robustfit (years1,dp1);



load 'uh5ydpm'
M = uh5ydpm;

mat2 = M (2:12,:);

years2 = mat2 (:,1);
dp2 = mat2 (:,2);

bb2 = robustfit (years2,dp2);


load 'uh5ydpm'
M = uh5ydpm;

mat3 = M (3:13,:);

years3 = mat3 (:,1);
dp3 = mat3 (:,2);

bb3 = robustfit (years3,dp3);


load 'uh5ydpm'
M = uh5ydpm;

mat4 = M (4:14,:);

years4 = mat4 (:,1);
dp4 = mat4 (:,2);

bb4 = robustfit (years4,dp4);



load 'uh5ydpm'
M = uh5ydpm;

mat5 = M (5:15,:);

years5 = mat5 (:,1);
dp5 = mat5 (:,2);

bb5 = robustfit (years5,dp5);


load 'uh5ydpm'
M = uh5ydpm;

mat6 = M (6:16,:);

years6 = mat6 (:,1);
dp6 = mat6 (:,2);

bb6 = robustfit (years6,dp6);


load 'uh5ydpm'
M = uh5ydpm;

mat7 = M (7:17,:);

years7 = mat7 (:,1);
dp7 = mat7 (:,2);

bb7 = robustfit (years7,dp7);



load 'uh5ydpm'
M = uh5ydpm;

mat8 = M (8:18,:);

years8 = mat8 (:,1);
dp8 = mat8 (:,2);

bb8 = robustfit (years8,dp8);


load 'uh5ydpm'
M = uh5ydpm;

mat9 = M (9:19,:);

years9 = mat9 (:,1);
dp9 = mat9 (:,2);

bb9 = robustfit (years9,dp9);


load 'uh5ydpm'
M = uh5ydpm;

mat10 = M (10:20,:);

years10 = mat10 (:,1);
dp10 = mat10 (:,2);

bb10 = robustfit (years10,dp10);



load 'uh5ydpm'
M = uh5ydpm;

mat11 = M (11:21,:);

years11 = mat11 (:,1);
dp11 = mat11 (:,2);

bb11 = robustfit (years11,dp11);



load 'uh5ydpm'
M = uh5ydpm;

mat12 = M (12:22,:);

years12 = mat12 (:,1);
dp12 = mat12 (:,2);

bb12 = robustfit (years12,dp12);



load 'uh5ydpm'
M = uh5ydpm;

mat13 = M (13:23,:);

years13 = mat13 (:,1);
dp13 = mat13 (:,2);

bb13 = robustfit (years13,dp13);


load 'uh5ydpm'
M = uh5ydpm;

mat14 = M (14:24,:);

years14 = mat14 (:,1);
dp14 = mat14 (:,2);

bb14 = robustfit (years14,dp14);

% 
% load 'uh5ydpm'
% M = uh5ydpm;
% 
% mat15 = M (15:25,:);
% 
% years15 = mat15 (:,1);
% dp15 = mat15 (:,2);
% 
% bb15 = robustfit (years15,dp15);



trendchangesplot = figure;
hold on
trenchangesplot = plot (years1,years1*bb1(2)+bb1(1),'blue','Linewidth',2);
hold on
grid on
trenchangesplot = plot (years2,years2*bb2(2)+bb2(1),'blue','Linewidth',2);
hold on
grid on
trenchangesplot = plot (years3,years3*bb3(2)+bb3(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years4,years4*bb4(2)+bb4(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years5,years5*bb5(2)+bb5(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years6,years6*bb6(2)+bb6(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years7,years7*bb7(2)+bb7(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years8,years8*bb8(2)+bb8(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years9,years9*bb9(2)+bb9(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years10,years10*bb10(2)+bb10(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years11,years11*bb11(2)+bb11(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years12,years12*bb12(2)+bb12(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years13,years13*bb13(2)+bb13(1),'blue','Linewidth',2);
hold on
%trenchangesplot = plot (years14,years14*bb14(2)+bb14(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('dew point (Celsius)','Fontsize',12);


load 'uh5ydpm'
MM = uh5ydpm;

load 'bb'
bb = bb;

years = MM (:,1);
dp = MM (:,2);

trendchengesmean = figure;
hold on
trendchangesmean = plot (years, dp, 'cyan','Linewidth',2);
hold on
trenschangesmean = plot (years, years*bb(2) + bb(1),'magenta','Linewidth',2);
hold on
grid on
trenchangesmean = plot (years1,years1*bb1(2)+bb1(1),'blue','Linewidth',2);
hold on
grid on
trenchangesmean = plot (years2,years2*bb2(2)+bb2(1),'blue','Linewidth',2);
hold on
grid on
trenchangesmean = plot (years3,years3*bb3(2)+bb3(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years4,years4*bb4(2)+bb4(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years5,years5*bb5(2)+bb5(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years6,years6*bb6(2)+bb6(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years7,years7*bb7(2)+bb7(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years8,years8*bb8(2)+bb8(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years9,years9*bb9(2)+bb9(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years10,years10*bb10(2)+bb10(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years11,years11*bb11(2)+bb11(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years12,years12*bb12(2)+bb12(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years13,years13*bb13(2)+bb13(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years14,years14*bb14(2)+bb14(1),'blue','Linewidth',2);
hold on
%trenchangesmean = plot (years15,years15*bb15(2)+bb15(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('dew point (Celsius)','Fontsize',12);
hold on
axis ([1982 2012 -22 -12]);
hold on
legend ('5 yr running average','-0.09 + 171.19 ((C) / yr)','10 yrs trends');