Global Index (short | long) | Local contents | Local Index (short | long)
| This script calls | |
|---|---|
clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[uuw, vvw, uvw, uw, vw] = ...
getnc('wct_globsst_t31.nc', 'UU', 'VV', 'UV', 'U', 'V');
[uuc, vvc, uvc, uc, vc] = ...
getnc('cct_globsst_t31.nc', 'UU', 'VV', 'UV', 'U', 'V');
[lat, lon, lev] = getll('wct_globsst_t31.nc');
ew1 = vvw - vw.^2 - uuw + uw.^2;
ec1 = vvc - vc.^2 - uuc + uc.^2;
ew2 = uvw - uw.*vw;
ec2 = uvc - uc.*vc;
get_global
num = 8;
top1 = squeeze(mean(ew1(:,num, :, :)));
top2 = squeeze(mean(ew2(:,num, :, :)));
bot1 = squeeze(mean(ec1(:,num, :, :)));
bot2 = squeeze(mean(ec2(:,num, :, :)));
FRAME = [0 360 -90 90];
figure(1);
subplot(1,1,1);
top = top1 - bot1;
bot = bot1 - bot2;
FRAME = [0 360 -90 90];
gquiv(top, bot, 2);
figure(1)
XAX = lon;
YAX = lat;
sp(1)
gquiv(top1, top2, 2);
dc
sp(2)
gquiv(bot1, bot2, 2);
dc
divw = sph_div(top1, top2, lon, lat);
divc = sph_div(bot1, bot2, lon ,lat);
figure(1)
sp(1)
gcont(divw, [-2:.5:2]*1e-4);
dc
sp(2)
gcont(divc, [-2:.5:2]*1e-4);
dc
figure(2)
gcont((divw-divc), [-2:.5:2]*1e-4)
dc