Global Index (short | long) | Local contents | Local Index (short | long)
height=load('geop500mb98111000.txt');
trying to load a more convenient way^M
| This script calls | |
|---|---|
clean
cd /home/disk/synoptic/class/atm442/matlab_data_nov98
load geop500mb98111012.txt ; %load geopotential height data set
latr = [20 : 1.25 : 60]; %latitude range of phi grid
lonr = [220 : 1.25 : 300]; %longitude range of phi grid
phi=geop500mb98111012*9.81 ; % phi is the geopotential
height=geop500mb98111012 ; % height is an easier name to use
gr = 0.8; % define gray for coastline and state boundaries
load coast; %data for coastlines
figure(1); fl(1); clf % plot a mercator projection of US
clow=4680; chigh=6000; cint=60; %contour line limits for 850 phi
conts = [clow:cint:chigh]; %define contour lines
V = [clow:cint*2:chigh]; % label every 2th line
latw = [20 60]; lonw = [220 300]; %limits of lat and long on map
mapp = 'mercator';
axesm('mapprojection',mapp,'maplatlimit',latw,'maplonlimit',lonw)
plotm(coast,'color',[gr gr gr],'linewidth',2) %plots coastline
displaym(worldlo('POline'))
displaym(usalo('stateborder'))
gridm; plabel; mlabel
[c, h] = pncontm(lonr, latr, height, conts, 0, 'k');
clabelm(c, h, V);
load urel500mb98111012.txt
uwnd = urel500mb98111012;
load vrel500mb98111012.txt
vwnd = vrel500mb98111012;
hold on;
quiverm2(thin(latr'*ones(size(lonr)), 2), ...
thin(ones(size(latr'))*lonr, 2), ...
thin(uwnd, 2), thin(vwnd, 2), 'k');
hold off;
vor = sph_curl(uwnd, vwnd, lonr, latr);
lev2 = 1e-4*[0.5:.5:3];
color_under(latr, lonr, vor, lev2, '+');
color_under(latr, lonr, vor, -lev2, '-');