if 0

gps = gemreader0;
load sethgutman/seth
load ../aerigoes/aerigoes_values
load ../goes/goes_values
%load ../sondes_uwvai/uncor/uwvai
load ../sondes_uwvai/cor/andros/uwvai_cor

end

h=plot(dn(gps.x),gps.y/10,'b',...
       dn(gps.nss1970),gps.pwv/10,'ob',...
       dn(ag.nss1970),ag.pwv,'or',...
       dn(goes.nss1970),goes.pwv,'sc',...
       dn(uwvai_cor.nss1970),uwvai_cor.pwv_ucor,'or',...
       dn(uwvai_cor.nss1970),uwvai_cor.pwv_cor,'or',...
       dn(seth.nss1970),seth.ipwv,'^k');
set(h(2),'MarkerFaceColor','b','MarkerEdgeColor','k')
set(h(3),'MarkerFaceColor',[0 .8 0],'MarkerEdgeColor','k')
set(h(4),'MarkerFaceColor',[1 .8 0],'MarkerEdgeColor','k')
set(h(5),'MarkerFaceColor','r','MarkerEdgeColor','k','MarkerSize',8)
set(h(6),'MarkerFaceColor',[1 .7 .7],'MarkerEdgeColor','k','MarkerSize',8)
set(h(7),'MarkerFaceColor',[.8 .8 .8],'MarkerEdgeColor','k','MarkerSize',7)
set(gca,'FontSize',18)
datetick;
ylabel('pwv (cm)')

for i = 1:60
  close
  clf

  ind0 = find(gps.x >= nosecsince1970(1998,8,i,0) & ...
              gps.x < nosecsince1970(1998,8,i+1,0));  
  ind1 = find(gps.nss1970 >= nosecsince1970(1998,8,i,0) & ...
              gps.nss1970 < nosecsince1970(1998,8,i+1,0));  
  ind2 = find(uwvai_cor.nss1970 >= nosecsince1970(1998,8,i,0) & ...
              uwvai_cor.nss1970 < nosecsince1970(1998,8,i+1,0));
  ind3 = find(ag.nss1970 >= nosecsince1970(1998,8,i,0) & ...
              ag.nss1970 < nosecsince1970(1998,8,i+1,0));
  ind4 = find(goes.nss1970 >= nosecsince1970(1998,8,i,0) & ...
              goes.nss1970 < nosecsince1970(1998,8,i+1,0));
  ind5 = find(seth.nss1970 >= nosecsince1970(1998,8,i,0) & ...
              seth.nss1970 < nosecsince1970(1998,8,i+1,0));

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

    bd = nosecsince1970(1998,8,i,0);    
    h=plot((gps.x(ind0)-bd)/3600,gps.y(ind0)/10,'b');
    hold on
    h=plot((gps.nss1970(ind1)-bd)/3600,gps.pwv(ind1)/10,'ob');
    set(h,'MarkerFaceColor','b','MarkerEdgeColor','k')
    h=plot((seth.nss1970(ind5)-bd)/3600,seth.ipwv(ind5),'^k');
    set(h,'MarkerFaceColor',[.7 .7 .7],'MarkerEdgeColor','k','MarkerSize',8)
    h=plot((ag.nss1970(ind3)-bd)/3600,ag.pwv(ind3),'og');
    set(h,'MarkerFaceColor',[0 .8 0],'MarkerEdgeColor','k')
    h=plot((goes.nss1970(ind4)-bd)/3600,goes.pwv(ind4),'sc');
    set(h,'MarkerFaceColor',[1 .8 0],'MarkerEdgeColor','k','MarkerSize',8)
    h=plot((uwvai_cor.nss1970(ind2)-bd)/3600,uwvai_cor.pwv_ucor(ind2),'or');
    set(h,'MarkerFaceColor','r','MarkerEdgeColor','k','MarkerSize',10)
    h=plot((uwvai_cor.nss1970(ind2)-bd)/3600,uwvai_cor.pwv_cor(ind2),'or');
    set(h,'MarkerFaceColor',[1 .7 .7],'MarkerEdgeColor','k','MarkerSize',10)
    set(gca,'FontSize',18)
    ylabel('pwv (cm)')
    xlabel('decimal hour')
    set(gca,'XLim',[0 24])
    set(gca,'YLim',[0 8])
    title(['Andros Island pwv: ' datestr(dn(nosecsince1970(1998,8,i,0)))])

    h=text(1.5,2.1,'GPS, processing b');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,2.1,'^k','MarkerFaceColor',[.7 .7 .7],'MarkerSize',8);
    h=text(1.5,2.6,'GPS, processing a');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,2.6,'ok','MarkerFaceColor','b');
    h=plot([1.0 1.0],[2.4 2.8],'b');
    h=text(1.5,1.6,'UW Vaisala sonde');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,1.6,'ok','MarkerFaceColor','r','MarkerSize',10);
    h=text(1.5,1.1,'UW Vaisala sonde with packaging correction');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,1.1,'ok','MarkerFaceColor',[1 .7 .7],'MarkerSize',10);
    h=text(1.5,0.6,'AERI/GOES');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,0.6,'ok','MarkerFaceColor',[0 .8 0]);
    h=text(1.5,0.1,'GOES');set(h,'FontWeight','b','FontSize',14,'Color','k');
    h=plot(1.0,0.1,'sk','MarkerFaceColor',[1 .8 0],'MarkerSize',8);

    drawnow
    
    fname = ['../plots/PS/' datestr(dn(bd)) '.ps'];
    eval(['print -dpsc ' fname])
    fname = ['../plots/JPG/' datestr(dn(bd)) '.jpg'];
    printimg(fname,11,8)
  end
end