% figure(1); clf; % giffig
 
% load wcoastline.dat
mapaxis=[-122.5 -121.75 36.6 37.1]; 
% mapaxis =[-122.08 -121.78 36.57 37]; 
[fillLat,fillLon] = get_landfill(mapaxis);

fill(fillLon,fillLat,[.8 .8 .8]); axis(mapaxis);
% plot(wcoastline(:,1),wcoastline(:,2),'r')
hold on; plot_iso([200 500 800],mapaxis,'b');
llgrid(20); set(gca,'tickdir','out','fontsize',8,'Xaxislocation','top');

%% transect
% load c:\muse\aviris\lowalt1
% y=y-.022; 	% flightline 2;
% [dist phas]=sw_dist(y,x,'km'); 
% p=plot(x,y,'r'); set(p,'linewidth',2);

% NHlat=[37+2.7/60 36+53.9/60]; NHlon=[-122-21.9/60 -122-10.4/60];
% p=plot(NHlon,NHlat,'b'); set(p,'linewidth',2);


