Global Index (short | long) | Local contents | Local Index (short | long)
o = owd; v = vwd; u = uwd;
| This script calls | |
|---|---|
clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
var1 = 'OMEGA'; var2 = 'CMFDT'; var3 = 'U'; var4 = 'PS'; var5='ZMDT';
[ow1, vw1, uw1, pw1, tw1] = ...
getnc('wgr_5093.nc', var1, var2, var3, var4, var5);
[oc1, vc1, uc1, pc1, tc1] = ...
getnc('cgr_5093.nc', var1, var2, var3, var4, var5);
[ow2, vw2, uw2, pw2, tw2] = ...
getnc('wgr_SOM_hard.nc', var1, var2, var3, var4, var5);
[oc2, vc2, uc2, pc2, tc2] = ...
getnc('cgr_SOM_hard.nc', var1, var2, var3, var4, var5);
[lat, lon, lev, hyam, hybm, P0] = getll('wgr_5093.nc');
[ntim, nlev, nlat, nlon] = size(ow1);
vw1 = vw1+tw1; clear tw1;
vc1 = vc1+tc1; clear tc1;
vw2 = vw2+tw2; clear tw2;
vc2 = vc2+tc2; clear tc2;
[xk, yk] = keep_var([0 360 -22.5 -7.5], lon, lat);
tim = [1:4 11 12];
owd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(ow1(tim,:,yk,xk))),1))),1);
ocd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(oc1(tim,:,yk,xk))),1))),1);
ows = shiftdim(squeeze(mean(shiftdim(squeeze(mean(ow2(tim,:,yk,xk))),1))),1);
ocs = shiftdim(squeeze(mean(shiftdim(squeeze(mean(oc2(tim,:,yk,xk))),1))),1);
vwd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(vw1(tim,:,yk,xk))),1))),1);
vcd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(vc1(tim,:,yk,xk))),1))),1);
vws = shiftdim(squeeze(mean(shiftdim(squeeze(mean(vw2(tim,:,yk,xk))),1))),1);
vcs = shiftdim(squeeze(mean(shiftdim(squeeze(mean(vc2(tim,:,yk,xk))),1))),1);
uwd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(uw1(tim,:,yk,xk))),1))),1);
ucd = shiftdim(squeeze(mean(shiftdim(squeeze(mean(uc1(tim,:,yk,xk))),1))),1);
uws = shiftdim(squeeze(mean(shiftdim(squeeze(mean(uw2(tim,:,yk,xk))),1))),1);
ucs = shiftdim(squeeze(mean(shiftdim(squeeze(mean(uc2(tim,:,yk,xk))),1))),1);
pwd = squeeze(mean(shiftdim(squeeze(mean(pw1)),1)));
pcd = squeeze(mean(shiftdim(squeeze(mean(pc1)),1)));
pws = squeeze(mean(shiftdim(squeeze(mean(pw2)),1)));
pcs = squeeze(mean(shiftdim(squeeze(mean(pc2)),1)));
lev2 = [25 50:50:1000];
owd = atlev(owd, lev2, pwd, hyam, hybm, P0);
ocd = atlev(ocd, lev2, pwd, hyam, hybm, P0);
ows = atlev(ows, lev2, pwd, hyam, hybm, P0);
ocs = atlev(ocs, lev2, pwd, hyam, hybm, P0);
vwd = atlev(vwd, lev2, pwd, hyam, hybm, P0);
vcd = atlev(vcd, lev2, pwd, hyam, hybm, P0);
vws = atlev(vws, lev2, pwd, hyam, hybm, P0);
vcs = atlev(vcs, lev2, pwd, hyam, hybm, P0);
uwd = atlev(uwd, lev2, pwd, hyam, hybm, P0);
ucd = atlev(ucd, lev2, pwd, hyam, hybm, P0);
uws = atlev(uws, lev2, pwd, hyam, hybm, P0);
ucs = atlev(ucs, lev2, pwd, hyam, hybm, P0);
for i = 1:4
if i == 1;
figure(1);% figure_landscape;
o = owd - ocd; v = vwd - vcd; u = uwd - ucd;
elseif i == 2;
o = ows - ocs; v = vws - vcs; u = uws - ucs;
% o = ocd; v = vcd; u = ucd;
elseif i == 3;
figure(2);% figure_landscape;
o = ows - owd; v = vws - vwd; u = uws - uwd;
% o = ows; v = vws; u = uws;
elseif i == 4;
o = ocs - ocd; v = vcs - vcd; u = ucs - ucd;
% o = ocs; v = vcs; u = ucs;
end
sp(rem(i+1,2)+1); cla
quiver(lon, -0.3*lev2', u, -500*o, 1.5, '-k');
hold on
pncont(lon, -0.3*lev2', v*1e6, [-20:2:20], 0, 'k')
hold off
set(gca, 'XTick', [0:30:360], 'YTick', -300:30:0, 'YTickLabel', 1000:-100:0)
axis([0 360 -300 -25])
if i == 1;
titlab = ['(WGR\_DOM - CGR\_DOM)'];
elseif i == 2;
titlab = ['(WGR\_SOM - CGR\_SOM)'];
elseif i == 3;
titlab = ['(WGR\_SOM - WGR\_DOM)'];
elseif i == 4;
titlab = ['(CGR\_SOM - CGR\_DOM)'];
end
title(['Zonal Cross Section: Y = -22.5:-7.5, ' titlab]);
xlabel('T-Tendency due to Convection Contoured, Contour Interval: 2e^-^6 K/S')
end
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots/DOM_SOM