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