Permalink
Switch branches/tags
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
59 lines (51 sloc) 1.44 KB
function varargout=plotplates(c11,cmn,res)
% [handl,XY]=PLOTPLATES(c11,cmn,res)
%
% Plots plate boundaries contained within 'c11' and 'cmn'
%
% res : 0 regular (default)
% : 1 regular (for compatibility with PLOTCONT)
% : 2 Global Mollweide projection
%
% Last modified by fjsimons-at-alum.mit.edu, 24.11.2004
defval('c11',[0 90])
defval('cmn',[360 -90])
defval('res',0)
switch res
case {0,1}
pathname=getenv('IFILES');
fid=fopen(fullfile(pathname,'PLATES','plates.mtl'),'r','b');
plates=fread(fid,[1692 2 ],'uint16');
plates(plates==0)=NaN;
plates=plates/100-90;
lon=plates(:,1);
lat=plates(:,2);
lon(~(lon>=c11(1) & lon<=cmn(1)))=NaN;
lat(~(lat<=c11(2) & lat>=cmn(2)))=NaN;
case 2
load(fullfile(getenv('IFILES'),'COASTS','platm'))
lon=plxm;
lat=plym;
end
hold on
handl=plot(lon,lat,'k');
XY=[lon(:) lat(:)];
varnames={'handl' 'XY'};
for index=1:nargout
varargout{index}=eval(varnames{index});
end
%file=fullfile(getenv('IFILES'),'PLATES','nuvel1_plates.datp');
%[err,dat]=unix([ 'awk ''{if($1~/>/) {print "NaN NaN"} else {print}}'' ',...
% file]);
%dat=str2num(dat);
%[m,n]=size(dat);
% Look at [(1:m)' dat] find out where NaN's need be inserted.
%positions=[24 28 44 106 410];
%lat=insert(dat(:,1),NaN,positions);
%lon=insert(dat(:,2),NaN,positions);
%plot(lon,lat)
%plates=([lon' lat']+90)*100;
%pathname=getenv('IFILES');
%fid=fopen([pathname '/PLATES/plates.mtl'],'wb');
%fwrite(fid,plates,'uint16')
%fclose(fid)