indx=data~=-9999;
agg=data(indx);
x=datenum(2018,1,1): datenum(2018,12,31);
f=31+1:31+28;
plot(x(f(indx(f))),(data(f((indx(f))))))
datetick('x','dd-mmm-yyyy')