Skip to content
Permalink
Branch: master
Find file Copy path
Fetching contributors…
Cannot retrieve contributors at this time
91 lines (77 sloc) 2 KB
function sdwspace
% SDWSPACE
%
% Simons, Dahlen and Wieczorek (2006) - Figure 5.1
% Plots spatial functions for spherical polar cap
% for various angular orders with eigenvalue grey shading.
%
% Last modified by fjsimons-at-alum.mit.edu, August 19th, 2004
TH=40; TH=20
nth=128;
L=18; SN=4;
M=4;
clf
[ah,ha]=krijetem(subnum(M+1,SN));
for m=0:M
[E1,V1{m+1},N1,th1,C1,ngl1,ngl2]=sdwcap(TH,L,m,nth);
% Note that we are requesting but not using C2, never mind the error
[E2,V2{m+1},N2,th2,C2]=sdwcapt(TH,L,m,nth);
E1m(m+1)=max(abs(E1(:)));
E2m(m+1)=max(abs(E2(:)));
for ondex=1:SN
axes(ah(m*SN+ondex))
% pm(m+1,ondex)=plot(th1,E1(:,ondex),'k-');
% hold on
% pg(m+1,ondex)=plot(th2,E2(:,ondex),'-','Color',grey(6));
pm(m+1,ondex)=plot(th1,E1(:,ondex),'-','Color',grey(6));
hold on
pg(m+1,ondex)=plot(th2,E1(1:size(E2,1),ondex),'-','Color','k');
drawnow
% set(ah(m*SN+ondex),'Color',grey(V1(ondex)*10))
% if V1(ondex)<0.4
% set(pm(m+1,ondex),'Color','w')
% end
end
end
set(ah,'xlim',[0 180],'xtick',[0 TH 180],'xgrid','on','ygrid','on',...
'ylim',[-1.1 1.1]*max([E1m E2m]))
nolabels(ah(1:end-SN),1)
nolabels(ha(M+2:end),2)
%nolabels(ha(1:M+1),2)
longticks(ah,1/2)
% Eigenvalue labels
nf=9;
for m=0:M
for ondex=1:SN
axes(ah(m*SN+ondex))
t{m+1,ondex}=sprintf('%s = %9.6f','\lambda',V1{m+1}(ondex));
if V1{m+1}(ondex)<0.975
lox='lr';
else
lox='ur';
end
[bh(m+1,ondex),th(m+1,ondex)]=...
boxtex(lox,ah(m*SN+ondex),t{m+1,ondex},nf);
end
end
set(th,'FontS',nf-1)
for index=SN*M+1:SN*(M+1)
axes(ah(index))
deggies(ah(index),1)
xl1(index)=xlabel(sprintf('colatitude %s','\theta'));
end
for index=1:M+1
serre(ah([1:SN]+(index-1)*SN),1/3)
axes(ha(index))
ylb(index)=ylabel(sprintf('m = %i',index-1));
end
for index=1:SN
axes(ah(index))
tlb(index)=title(sprintf('%s = %i','\alpha',index));
end
set([pm(:) ; pg(:)],'LineW',1)
set([ah],'Fonts',nf)
fig2print(gcf,'portrait')
% set(gcf,'Color','w','Inv','off')
shrink(ah,1,1/1.1)
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.