Vous êtes ici : Accueil / Alexandre Ganachaud / xbt / cmpxbttsg / xbttsg_pltbpoints.m

xbttsg_pltbpoints.m

Par Webmaster Legos Dernière modification 03/07/2012 18:33

Objective-C source code icon xbttsg_pltbpoints.m — Objective-C source code, 1 kB (1319 bytes)

Contenu du fichier

%plot bad points geographically along lines

corder=['b+';'b.';'bx';'bs';'g+';'g.';'gx';'c+';'c.';...
        'cx';'m+';'m.';'mx'];

for idev=1:length(gdevthxb)
  figure(idev);clf
  m_proj('Gall-Peters','long',[0 360],'lat',[-50 80]);
  hold on
  m_coast('line','color','k')
  m_grid('box','fancy','xticklabels',[],'yticklabels',[]);
  for iship =1: length(thxbt)
    clear xbtmean
    thg=thxbt(iship).thg;
    xbt=thxbt(iship).xbt;
    nxbt=length(xbt);
    tlat=cat(1,thg.thlat);
    tlon=cat(1,thg.thlon);
    for ix=1:nxbt
      xbtmean(ix)=mean(xbt(ix).xsst);
    end
    gixdev=find(abs(cat(1,thg.thsst)-xbtmean')>gdevthxb(idev));
    l2shift=find(tlon<0);
    tlon(l2shift)=tlon(l2shift)+360;
    [X,Y]=m_ll2xy(tlon,tlat,'clip','on');
    pl(iship)=plot(X,Y,corder(iship,:),'markersize',2);
    plot(X(gixdev),Y(gixdev),'o','marker','o','markersize',4,...
         'MarkerFaceColor','r',...
            'MarkerEdgeColor','r');
  end %on iship
  title(sprintf('Deviation XBT-TSG > %d deg.',gdevthxb(idev)));
  legend(pl,cat(1,thxbt.shipid))
  drawnow;
  land;setlargefig;
%  eval(sprintf('print -depsc fig_map_tsgxbtdev%3.1f.eps',...
%               gdevthxb(idev)))
  eval(sprintf('print -djpeg90 fig_map_tsgxbtdev%3.1f.jpg',...
               gdevthxb(idev)))
end % idev

Actions sur le document

logo cnes logo IRD Logo université de Toulouse Logo université Paul Sabatier Logo CNRS
Logo bibliothèque OBS Logo Observatoire Midi Pyrénées