colors = flipud(jet(20));

subtitle_comment = config.comment;

for ii=1:4
   eval(sprintf('ax%d = subplot(4,1,%d);', ii, ii))
   pos = get(gca, 'pos');
   set(gca, 'pos', pos + [0 0 0 .03]);
   if (ii < 4)
      set(gca, 'xtick', [])
   end
end

%%%% uabs from gps
axes(ax1)
nbins = length(data.bins);
plot([min(data.dday) max(data.dday)]',[0 0],'b');hold on
plot(data.dday,data.uship_g,'g','linewidth',2);hold on
for swathi = 1:length(colors(:,1))
   ii = (swathi - 1)*2 + [1:2];
   plot(data.dday, -data.umeas(ii,:), 'color',colors(swathi,:))
   hold on
end
xl(min(data.dday), max(data.dday))
yl(-8,8)
ylabel('shipu and -measu')



%%% measured velocity from package
axes(ax2);
plot([min(data.dday) max(data.dday)]',[0 0],'b');hold on
plot(nav_dday, nav_umeas, 'r')
hold on
plot(data.dday, data.uship_g, 'g', 'linewidth', 2);hold on
plot(data.dday, data.uship_sm, 'k', 'linewidth', 1);
ylabel('ship U and refl_g: m/s')
xl(min(data.dday), max(data.dday))
yl(-.5,.5)

axes(ax3)
plot([min(data.dday) max(data.dday)]',[0 0],'b');hold on
plot(data.dday,data.vship_g,'g','linewidth',2);hold on
for swathi = 1:length(colors(:,1))
   ii = (swathi - 1)*2 + [1:2];
   plot(data.dday, -data.vmeas(ii,:), 'color',colors(swathi,:))
   hold on
end
xl(min(data.dday), max(data.dday))
yl(-8,8)
ylabel('shipv and -measv')


%%% measured velocity from package
axes(ax4);
plot([min(data.dday) max(data.dday)]',[0 0],'b');hold on
plot(nav_dday, nav_vmeas, 'r')
hold on
plot(data.dday, data.vship_g, 'g', 'linewidth', 2);hold on
plot(data.dday, data.vship_sm, 'k', 'linewidth', 1);
ylabel('ship V and refl_g: m/s')
xl(min(data.dday), max(data.dday))
yl(-.5,.5)


date_it
subtitle_comment = sprintf('%s (dday %3.1f - %3.1f)',...
                           subtitle_comment, min(data.dday), max(data.dday));

subtitle(subtitle_comment)

