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