f4=figure(4)
box on
yyaxis left
hold on
plot(MX,pkdosat_ns_1d,'-','color',[.4 .6 .6],'linewidth',2)
plot(MX,DO3_1d,'-','color','b','linewidth',2)
plot(MX,poDO4_1d,'-.','color','b','linewidth',2)
ylim([90 450])
ylabel('oxygen')
yyaxis right
hold on
plot(date_nf_off,newfrond,'.','markersize',37,'Color','r')
ylabel('New Fronds')
ylim([-1.5 7])
ax=gca;
ax.YAxis(1).Color='k'
ax.YAxis(2).Color='k'
set(gca,'fontsize',16)
x = datenum('June-27-2018'):7:datenum('Oct-3-2018');
set(gca, 'XTick', x);
datetick('x','mm/dd','keepticks');