load 'MLOprec5y'
M = MLOprec5y;

mat1 = M (1:11,:);

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

bb1 = robustfit (years1,T1);



load 'MLOprec5y'
M = MLOprec5y;

mat2 = M (2:12,:);

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

bb2 = robustfit (years2,T2);


load 'MLOprec5y'
M = MLOprec5y;

mat3 = M (3:13,:);

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

bb3 = robustfit (years3,T3);


load 'MLOprec5y'
M = MLOprec5y;

mat4 = M (4:14,:);

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

bb4 = robustfit (years4,T4);



load 'MLOprec5y'
M = MLOprec5y;

mat5 = M (5:15,:);

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

bb5 = robustfit (years5,T5);


load 'MLOprec5y'
M = MLOprec5y;

mat6 = M (6:16,:);

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

bb6 = robustfit (years6,T6);


load 'MLOprec5y'
M = MLOprec5y;

mat7 = M (7:17,:);

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

bb7 = robustfit (years7,T7);



load 'MLOprec5y'
M = MLOprec5y;

mat8 = M (8:18,:);

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

bb8 = robustfit (years8,T8);


load 'MLOprec5y'
M = MLOprec5y;

mat9 = M (9:19,:);

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

bb9 = robustfit (years9,T9);


load 'MLOprec5y'
M = MLOprec5y;

mat10 = M (10:20,:);

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

bb10 = robustfit (years10,T10);



load 'MLOprec5y'
M = MLOprec5y;

mat11 = M (11:21,:);

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

bb11 = robustfit (years11,T11);



load 'MLOprec5y'
M = MLOprec5y;

mat12 = M (12:22,:);

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

bb12 = robustfit (years12,T12);



load 'MLOprec5y'
M = MLOprec5y;

mat13 = M (13:23,:);

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

bb13 = robustfit (years13,T13);


load 'MLOprec5y'
M = MLOprec5y;

mat14 = M (14:24,:);

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

bb14 = robustfit (years14,T14);


load 'MLOprec5y'
M = MLOprec5y;

mat15 = M (15:25,:);

years15 = mat15 (:,1);
T15 = mat15 (:,2);

bb15 = robustfit (years15,T15);


load 'MLOprec5y'
M = MLOprec5y;

mat16 = M (16:26,:);

years16 = mat16 (:,1);
T16 = mat16 (:,2);

bb16 = robustfit (years16,T16);


load 'MLOprec5y'
M = MLOprec5y;

mat17 = M (17:27,:);

years17 = mat17 (:,1);
T17 = mat17 (:,2);

bb17 = robustfit (years17,T17);



load 'MLOprec5y'
M = MLOprec5y;

mat18 = M (18:28,:);

years18 = mat18 (:,1);
T18 = mat18 (:,2);

bb18 = robustfit (years18,T18);



load 'MLOprec5y'
M = MLOprec5y;

mat19 = M (19:29,:);

years19 = mat19 (:,1);
T19 = mat19 (:,2);

bb19 = robustfit (years19,T19);


load 'MLOprec5y'
M = MLOprec5y;

mat20 = M (20:30,:);

years20 = mat20 (:,1);
T20 = mat20 (:,2);

bb20 = robustfit (years20,T20);



load 'MLOprec5y'
M = MLOprec5y;

mat21 = M (21:31,:);

years21 = mat21 (:,1);
T21 = mat21 (:,2);

bb21 = robustfit (years21,T21);



load 'MLOprec5y'
M = MLOprec5y;

mat22 = M (22:32,:);

years22 = mat22 (:,1);
T22 = mat22 (:,2);

bb22 = robustfit (years22,T22);


load 'MLOprec5y'
M = MLOprec5y;

mat23 = M (23:33,:);

years23 = mat23 (:,1);
T23 = mat23 (:,2);

bb23 = robustfit (years23,T23);


load 'MLOprec5y'
M = MLOprec5y;

mat24 = M (24:34,:);

years24 = mat24 (:,1);
T24 = mat24 (:,2);

bb24 = robustfit (years24,T24);

load 'MLOprec5y'
M = MLOprec5y;

mat25 = M (25:35,:);

years25 = mat25 (:,1);
T25 = mat25 (:,2);

bb25 = robustfit (years25,T25);

load 'MLOprec5y'
M = MLOprec5y;

mat26 = M (26:36,:);

years26 = mat26 (:,1);
T26 = mat26 (:,2);

bb26 = robustfit (years26,T26);

load 'MLOprec5y'
M = MLOprec5y;

mat27 = M (27:37,:);

years27 = mat27 (:,1);
T27 = mat27 (:,2);

bb27 = robustfit (years27,T27);


load 'MLOprec5y'
M = MLOprec5y;

mat28 = M (28:38,:);

years28 = mat28 (:,1);
T28 = mat28 (:,2);

bb28 = robustfit (years28,T28);


load 'MLOprec5y'
M = MLOprec5y;

mat29 = M (29:39,:);

years29 = mat29 (:,1);
T29 = mat29 (:,2);

bb29 = robustfit (years29,T29);


load 'MLOprec5y'
M = MLOprec5y;

mat30 = M (30:40,:);

years30 = mat30 (:,1);
T30 = mat30 (:,2);

bb30 = robustfit (years30,T30);


load 'MLOprec5y'
M = MLOprec5y;

mat31 = M (31:41,:);

years31 = mat31 (:,1);
T31 = mat31 (:,2);

bb31 = robustfit (years31,T31);


load 'MLOprec5y'
M = MLOprec5y;

mat32 = M (32:42,:);

years32 = mat32 (:,1);
T32 = mat32 (:,2);

bb32 = robustfit (years32,T32);


load 'MLOprec5y'
M = MLOprec5y;

mat33 = M (33:43,:);

years33 = mat33 (:,1);
T33 = mat33 (:,2);

bb33 = robustfit (years33,T33);

load 'MLOprec5y'
M = MLOprec5y;

mat34 = M (34:44,:);

years34 = mat34 (:,1);
T34 = mat34 (:,2);

bb34 = robustfit (years34,T34);

load 'MLOprec5y'
M = MLOprec5y;

mat35 = M (35:45,:);

years35 = mat35 (:,1);
T35 = mat35 (:,2);

bb35 = robustfit (years35,T35);

load 'MLOprec5y'
M = MLOprec5y;

mat36 = M (36:46,:);

years36 = mat36 (:,1);
T36 = mat36 (:,2);

bb36 = robustfit (years36,T36);


load 'MLOprec5y'
M = MLOprec5y;

mat37 = M (37:47,:);

years37 = mat37 (:,1);
T37 = mat37 (:,2);

bb37 = robustfit (years37,T37);


load 'MLOprec5y'
M = MLOprec5y;

mat38 = M (38:48,:);

years38 = mat38 (:,1);
T38 = mat38 (:,2);

bb38 = robustfit (years38,T38);



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
trenchangesplot = plot (years15,years15*bb15(2)+bb15(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years16,years16*bb16(2)+bb16(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years17,years17*bb17(2)+bb17(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years18,years18*bb18(2)+bb18(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years19,years19*bb19(2)+bb19(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years20,years20*bb20(2)+bb20(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years21,years21*bb21(2)+bb21(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years22,years22*bb22(2)+bb22(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years23,years23*bb23(2)+bb23(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years24,years24*bb24(2)+bb24(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years25,years25*bb25(2)+bb25(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years26,years26*bb26(2)+bb26(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years27,years27*bb27(2)+bb27(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years28,years28*bb28(2)+bb28(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years29,years29*bb29(2)+bb29(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years30,years30*bb30(2)+bb30(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years31,years31*bb31(2)+bb31(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years32,years32*bb32(2)+bb32(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years33,years33*bb33(2)+bb33(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years34,years34*bb34(2)+bb34(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years35,years35*bb35(2)+bb35(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years36,years36*bb36(2)+bb36(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years37,years37*bb37(2)+bb37(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years38,years38*bb38(2)+bb38(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('precipitation (mm)','Fontsize',12);


load 'MLOprec5y'
MM = MLOprec5y;

load 'bbMLOprec'
bb = bbMLOprec;

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

trendchengesmean = figure;
hold on
trendchangesmean = plot (years, T, '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
trenchangesmean = plot (years16,years16*bb16(2)+bb16(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years17,years17*bb17(2)+bb17(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years18,years18*bb18(2)+bb18(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years19,years19*bb19(2)+bb19(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years20,years20*bb20(2)+bb20(1),'blue','Linewidth',2);
hold on
trenchangesmean = plot (years21,years21*bb21(2)+bb21(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years22,years22*bb22(2)+bb22(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years23,years23*bb23(2)+bb23(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years24,years24*bb24(2)+bb24(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years25,years25*bb25(2)+bb25(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years26,years26*bb26(2)+bb26(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years27,years27*bb27(2)+bb27(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years28,years28*bb28(2)+bb28(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years29,years29*bb29(2)+bb29(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years30,years30*bb30(2)+bb30(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years31,years31*bb31(2)+bb31(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years32,years32*bb32(2)+bb32(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years33,years33*bb33(2)+bb33(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years34,years34*bb34(2)+bb34(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years35,years35*bb35(2)+bb35(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years36,years36*bb36(2)+bb36(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years37,years37*bb37(2)+bb37(1),'blue','Linewidth',2);
hold on
trenchangesplot = plot (years38,years38*bb38(2)+bb38(1),'blue','Linewidth',2);
hold on
xlabel ('years','Fontsize',12);
hold on
ylabel ('precipitation (mm)','Fontsize',12);
hold on
axis ([1957 2012 0 1000]);
hold on
legend ('5 yr running average','3.45 yrs + 7.29 e3 (mm / yr)','10 yrs trends');