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