LD = load('Pul MODESTAMMEQ.mat');
MODESTAMMEQ = LD.MODESTAMMEQ;
dv = regexp(MODESTAMMEQ.Jan1989, '\d*\-\w*\-\d*','match');
MODESTAMMEQ.Var1 = datetime([dv{:}], 'InputFormat',"dd-MMM-yyyy").';
y = MODESTAMMEQ.Var2*10*0.35;
load('Pul MODESTA_PALINE_DELTA');
errorbar(x,MODESTAMMEQ.Var5(:,1),MODESTAMMEQ.Var4(:,1),'.r','MarkerSize',6)
plot(x,[y, MODESTAMMEQ.Var5(:,1)]);
title('Modesta AWS and Stakes')
title('Modesta AWS and Stakes')
load('Pul MODESTA_PALINE_DELTA');
errorbar(x(5432:6513,1),MODESTAMMEQ.Var5(5432:6513,1),MODESTAMMEQ.Var4(5432:6513,1),'.r','MarkerSize',6);
errorbar(x(7639:10569,1),MODESTAMMEQ.Var5(7639:10569,1),MODESTAMMEQ.Var4(7639:10569,1),'.b','MarkerSize',6);
errorbar(x(10570:11265,1),MODESTAMMEQ.Var5(10570:11265,1),MODESTAMMEQ.Var4(10570:11265,1),'.m','MarkerSize',6);
plot(MODESTAMMEQ.Var1(:,1),y);
title('Modesta AWS and Stakes')
load('Pul DATI_MAR_ANNUALI');
plot(MODESTAMMEQ.Var1(:,1),[y, MODESTAMMEQ.Var5(:,1)], 'g' , 'DisplayName','AWS')
plot(datetime(DATIMARannuali.Year,1,1), DATIMARannuali.SMB_priestleymm,'--co', 'MarkerEdgeColor','k','MarkerFaceColor','c', 'DisplayName','ECMWF');
plot(datetime(DATIECMWFannuali.Year,1,1), DATIECMWFannuali.SMB_priestleymm,'m-*', 'DisplayName','MAR');
legend('Location','best')