Skip to content
Permalink
Branch: master
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
77 lines (69 sloc) 1.88 KB
function sdwfried
% SDWFRIED
%
% Makes fried-egg plots for +- angular orders
% Figure 5.4 from Simons et al., SIAM Review (2006)
%
% Run twice. Second time round the color axis is OK
%
% Last modified by fjsimons-at-alum.mit.edu, 04/13/2009
TH=40;
L=18;
[lrnk,mrnk,lval]=sdwelm(TH,L);
% Nearly double the amount of requested tapers
lrnk=gamini(lrnk,(mrnk~=0)+1);
lval=gamini(lval,(mrnk~=0)+1);
mrnk=gamini(mrnk,(mrnk~=0)+1);
% This trick only works for the single cap
mrnk((diff(mrnk)==0))=-mrnk((diff(mrnk)==0));
[ah,ha]=krijetem(subnum(8,4));
for index=1:length(ah)
fnpl=sprintf('%s/SDWFRIED-%i-%i-%i.mat',...
fullfile(getenv('IFILES'),'LOCALIZE'),TH,L,index);
if exist(fnpl,'file')==2
load(fnpl)
disp(sprintf('%s loaded by SDWFRIED',fnpl))
axes(ah(index))
cbb=gray(10);
cbb=kelicol;
ispl=data/max(abs(data(:)));
ispl(abs(ispl)<1/100)=NaN;
imagefnan([-1 1],[1 -1],ispl,cbb,[-1 1])
colormap(cbb); hold on
ax2(1)=circ(1); ax2(2)=circ(sin(TH*pi/180));
set(ax2(1:2),'LineW',1)
set(ax2(2),'LineS','--')
axis([-1.0100 1.0100 -1.0100 1.0100])
axis off
if index==0
[axlim,handl,XYZ]=plotcont([0 90],[360 -90],4);
delete(handl)
hold on
p1=plot(XYZ(:,1),XYZ(:,2),'k-','LineWidth',1);
end
tl(index)=title(sprintf('%s_%i = %6.3f ; m = %i',...
'\lambda',...
lrnk(index),lval(index),mrnk(index)));
drawnow
else
[E,V,N,th,C]=sdwcap(TH,L,mrnk(index),128,[],2);
E=E{lrnk(index)};
Emax=max(E(:)); thresh=100;
Ed=E;
Ed(abs(Ed)<Emax/thresh)=0;
axes(ah(index))
[data,ch,ph]=plotplm(Ed,[],[],5,[],TH);
tl(index)=...
title(sprintf('%s = %6.3f ; m = %i','\lambda',...
lval(index),mrnk(index)));
delete(ch)
drawnow
save(fnpl,'data')
end
end
seemax(ah,3)
fig2print(gcf,'tall')
set(ah,'CameraV',6.25)
% Last-minute cosmetics
movev(tl,-0.25)
figdisp
You can’t perform that action at this time.
You signed in with another tab or window. Reload to refresh your session. You signed out in another tab or window. Reload to refresh your session.