Subversion Repositories lagranto.ecmwf

Rev

Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
3 michaesp 1
% -------------------------------------------------------------------------
2
% Plot horizontal starting positions
3
% -------------------------------------------------------------------------
4
 
5
% Load starting points
6
[ start.lon start.lat start.p ] = textread('TEST','%f %f %f',-1); 
7
 
8
% Open a new figure and set the geographical projection and region
9
figure(1);
10
clf;
11
load coast
12
h=axesm('MapProjection','stereo','origin',[ 90 70 ]);
13
gridm;
14
h=plotm(lat,long,'Color','k','LineWidth',1.5)
15
%axis([-1.5 1.5 -2 0]);
16
 
17
 
18
% Plot starting points
19
for i=1:length(start.lon)
20
  linem(start.lat(i),start.lon(i),'marker','o', ...
21
        'markersize',4,'color','w','MarkerEdgeColor',[.4 .4 .4], ...
22
        'MarkerFaceColor','b');
23
end
24
 
25
 
26
% Save figure
27
figname = [ 'startf.eps' ];
28
set(gcf, 'PaperPosition', [2 1 15 10]);
29
print('-depsc2','-r0',figname);
30