clf

  ind0 = find(gps.x >= nosecsince1970(1998,9,13,0) & ...
              gps.x < nosecsince1970(1998,9,15,0));
  ind1 = find(gps.nss1970 >= nosecsince1970(1998,9,13,0) & ...
              gps.nss1970 < nosecsince1970(1998,9,15,0));
  ind2 = find(uwvai_cor.nss1970 >= nosecsince1970(1998,9,13,0) & ...
              uwvai_cor.nss1970 < nosecsince1970(1998,9,15,0));
  ind3 = find(ag.nss1970 >= nosecsince1970(1998,9,13,0) & ...
              ag.nss1970 < nosecsince1970(1998,9,15,0));
  ind4 = find(goes.nss1970 >= nosecsince1970(1998,9,13,0) & ...
              goes.nss1970 < nosecsince1970(1998,9,15,0));
  ind5 = find(seth.nss1970 >= nosecsince1970(1998,9,13,0) & ...
              seth.nss1970 < nosecsince1970(1998,9,15,0));

  if length(ind1) >= 1 | length(ind2) >= 1 | length(ind3) >= 1 | ...
     length(ind4) >= 1 | length(ind5) >= 1

    h=plot(dn(gps.x(ind0)),gps.y(ind0)/10,'b');
    hold on
    h=plot(dn(gps.nss1970(ind1)),gps.pwv(ind1)/10,'ob');
    set(h,'MarkerFaceColor','b','MarkerEdgeColor','k')
    h=plot(dn(seth.nss1970(ind5)),seth.ipwv(ind5),'^k');
    set(h,'MarkerFaceColor',[.7 .7 .7],'MarkerEdgeColor','k','MarkerSize',8)
    h=plot(dn(ag.nss1970(ind3)),ag.pwv(ind3),'og');
    set(h,'MarkerFaceColor',[0 .8 0],'MarkerEdgeColor','k')
    h=plot(dn(goes.nss1970(ind4)),goes.pwv(ind4),'sc');
    set(h,'MarkerFaceColor',[1 .8 0],'MarkerEdgeColor','k','MarkerSize',8)
    h=plot(dn(uwvai_cor.nss1970(ind2)),uwvai_cor.pwv_ucor(ind2),'oc');
    set(h,'MarkerFaceColor','r','MarkerEdgeColor','k','MarkerSize',10)
    h=plot(dn(uwvai_cor.nss1970(ind2)),uwvai_cor.pwv_cor(ind2),'or');
    set(h,'MarkerFaceColor','m','MarkerEdgeColor','k','MarkerSize',10)
    set(gca,'FontSize',18)
    hold off
    datetick
    ylabel('pwv (cm)')
    set(gca,'YLim',[-0.2 6])
    title(['Andros Island pwv'])
    
    v=dn(nosecsince1970(1998,9,13,3,0,0));
    v2=dn(nosecsince1970(1998,9,13,4,0,0));
    
    hold on
    h=text(v2,2.1,'GPS, processing b');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,2.1,'^k','MarkerFaceColor',[.7 .7 .7],'MarkerSize',8);
    h=text(v2,2.6,'GPS, processing a');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,2.6,'ok','MarkerFaceColor','b');
    h=plot([v v],[2.4 2.8],'b');
    h=text(v2,1.6,'UW Vaisala sonde');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,1.6,'ok','MarkerFaceColor','r','MarkerSize',10);
    h=text(v2,1.1,'UW Vaisala sonde with packaging correction');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,1.1,'ok','MarkerFaceColor','m','MarkerSize',10);
    h=text(v2,0.6,'AERI/GOES');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,0.6,'ok','MarkerFaceColor',[0 .8 0]);
    h=text(v2,0.1,'GOES');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(v,0.1,'sk','MarkerFaceColor',[1 .8 0],'MarkerSize',8);

    drawnow
  end

ax=axis;
axis([ax(1) ax(2)-0.15*(ax(2)-ax(1)) ax(3) ax(4)])
grid
set(gca,'Position',[.2 .2 .6 .6])
orient landscape
print -depsc sept13-14.ps