testplot2

0001 dolandscape
0002 [ub,x,y]=testhisread('test.ubot.dat');       
0003 
0004 % create a suitable color map:
0005 %
0006 c=jet(32);
0007 cmap=c(20:29,:);
0008 
0009 %plot contours:
0010 
0011 extcontour(x/1000,y/1000,ub.*100,...
0012               [0 10 20 30 40 50],'fill','label');...
0013 colormap(cmap);...
0014 
0015 set(gca,'xlim',[690 770]);...
0016 set(gca,'ylim',[150 195]);...
0017 set(gca,'aspect',[1.77777 1]);
0018 
0019 hold on
0020 
0021 set(gca,'tickdir','out')
0022 
0023 titl='Bottom Orbital Velocity (cm/sec)'
0024 title(titl);
0025 xlabel('km');ylabel('km');
0026 
0027 %put a wind arrow on plot:
0028 
0029 dir=270.0;
0030 wind=8*exp(sqrt(-1)*dir*(pi/180));
0031 
0032 % draw a label box in which to put the wind arrow
0033 %
0034 %  ll=[755 180]; %lower left of label box
0035 %  ur=[765 190]; %upper right
0036 %  xbox=[ll(1) ur(1) ur(1) ll(1) ll(1)];
0037 %  ybox=[ll(2) ll(2) ur(2) ur(2) ll(2)];
0038 %  patch(xbox,ybox,'black');...
0039 
0040 %  draw and label the wind arrow
0041 %
0042   cg4=[.8 .8 .8];
0043   fatarrow(760,187,wind,cg4);...
0044   text(760,181,'Wind: 20 m/s','hor','center','fontsize',20,...
0045       'fontweight','bold','color',cg4);
0046