load 'MLOprec5y' M = MLOprec5y; mat1 = M (1:31,:); years1 = mat1 (:,1); T1 = mat1 (:,2); bb1 = robustfit (years1,T1); load 'MLOprec5y' M = MLOprec5y; mat2 = M (2:32,:); years2 = mat2 (:,1); T2 = mat2 (:,2); bb2 = robustfit (years2,T2); load 'MLOprec5y' M = MLOprec5y; mat3 = M (3:33,:); years3 = mat3 (:,1); T3 = mat3 (:,2); bb3 = robustfit (years3,T3); load 'MLOprec5y' M = MLOprec5y; mat4 = M (4:34,:); years4 = mat4 (:,1); T4 = mat4 (:,2); bb4 = robustfit (years4,T4); load 'MLOprec5y' M = MLOprec5y; mat5 = M (5:35,:); years5 = mat5 (:,1); T5 = mat5 (:,2); bb5 = robustfit (years5,T5); load 'MLOprec5y' M = MLOprec5y; mat6 = M (6:36,:); years6 = mat6 (:,1); T6 = mat6 (:,2); bb6 = robustfit (years6,T6); load 'MLOprec5y' M = MLOprec5y; mat7 = M (7:37,:); years7 = mat7 (:,1); T7 = mat7 (:,2); bb7 = robustfit (years7,T7); load 'MLOprec5y' M = MLOprec5y; mat8 = M (8:38,:); years8 = mat8 (:,1); T8 = mat8 (:,2); bb8 = robustfit (years8,T8); load 'MLOprec5y' M = MLOprec5y; mat9 = M (9:39,:); years9 = mat9 (:,1); T9 = mat9 (:,2); bb9 = robustfit (years9,T9); load 'MLOprec5y' M = MLOprec5y; mat10 = M (10:40,:); years10 = mat10 (:,1); T10 = mat10 (:,2); bb10 = robustfit (years10,T10); load 'MLOprec5y' M = MLOprec5y; mat11 = M (11:41,:); years11 = mat11 (:,1); T11 = mat11 (:,2); bb11 = robustfit (years11,T11); load 'MLOprec5y' M = MLOprec5y; mat12 = M (12:42,:); years12 = mat12 (:,1); T12 = mat12 (:,2); bb12 = robustfit (years12,T12); load 'MLOprec5y' M = MLOprec5y; mat13 = M (13:43,:); years13 = mat13 (:,1); T13 = mat13 (:,2); bb13 = robustfit (years13,T13); load 'MLOprec5y' M = MLOprec5y; mat14 = M (14:44,:); years14 = mat14 (:,1); T14 = mat14 (:,2); bb14 = robustfit (years14,T14); load 'MLOprec5y' M = MLOprec5y; mat15 = M (15:45,:); years15 = mat15 (:,1); T15 = mat15 (:,2); bb15 = robustfit (years15,T15); load 'MLOprec5y' M = MLOprec5y; mat16 = M (16:46,:); years16 = mat16 (:,1); T16 = mat16 (:,2); bb16 = robustfit (years16,T16); load 'MLOprec5y' M = MLOprec5y; mat17 = M (17:47,:); years17 = mat17 (:,1); T17 = mat17 (:,2); bb17 = robustfit (years17,T17); load 'MLOprec5y' M = MLOprec5y; mat18 = M (18:48,:); years18 = mat18 (:,1); T18 = mat18 (:,2); bb18 = robustfit (years18,T18); 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 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 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)','30 yrs trends');