Global Index (short | long) | Local contents | Local Index (short | long)
Start off with the typical variables.
| This script calls | |
|---|---|
clean
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[cpc, cpl, ch, cps] = getnc('cgr.nc', 'PRECC', 'PRECL', 'Z3', 'PS');
[wpc, wpl, wh, wps] = getnc('wgr.nc', 'PRECC', 'PRECL', 'Z3', 'PS');
[ipc, ipl, ih, ips] = getnc('india.nc', 'PRECC', 'PRECL', 'Z3', 'PS');
[clpc, clpl, clh, clps] = getnc('clim.nc', 'PRECC', 'PRECL', 'Z3', 'PS');
[lat, lon, lev, hyam, hybm, P0] = getll('cct.nc');
whos
cp = squeeze(mean(cpc+cpl))*(24*3600*1000);
wp = squeeze(mean(wpc+wpl))*(24*3600*1000);
ip = squeeze(mean(ipc+ipl))*(24*3600*1000);
clp = squeeze(mean(clpc+clpl))*(24*3600*1000);
ch = squeeze(mean(ch));
wh = squeeze(mean(wh));
ih = squeeze(mean(ih));
clh = squeeze(mean(clh));
cps = squeeze(mean(cps));
wps = squeeze(mean(wps));
ips = squeeze(mean(ips));
clps = squeeze(mean(clps));
levs = [1000 925 850 500 250];
chgt = atlev(ch, levs, cps, hyam, hybm, P0);
whgt = atlev(wh, levs, wps, hyam, hybm, P0);
ihgt = atlev(ih, levs, ips, hyam, hybm, P0);
clhgt = atlev(clh, levs, clps, hyam, hybm, P0);
% Try some plots: NOTE: The variable 'weight' is calculated
% somewhere else, and must be used here to scale the india run.
% In this case, just scroll down and you'll find it.
lev = 500;
ind = find(levs == lev);
top = squeeze(whgt(ind, :, :) - clhgt(ind, :, :))/1.0;
%bot = squeeze(whgt(ind, :, :) + chgt(ind, :, :) - 2*clhgt(ind, :, :));
%bot = squeeze(ihgt(ind, :, :) - clhgt(ind, :, :))/(1.4*weight);
%top = top - bot;
bot = squeeze(chgt(ind, :, :) - clhgt(ind, :, :))/1.0;
cind = 10;
origin = [90 270];
%origin = [-90 270];
%get_global
default_global
figure(1)
sp(1)
mcont(top, [-500:cind:500], origin)
title(['CT: WARM - CLIM ' num2str(lev) 'mb HGT']);
xlabel(['Contour Interval: ' num2str(cind) 'm']);
sp(2)
mcont(bot, [-500:cind:500], origin)
title(['CT: COLD - CLIM ' num2str(lev) 'mb HGT']);
xlabel(['Contour Interval: ' num2str(cind) 'm']);
% NOTE: Here's one way to get the variable 'weight'.
top = wp - clp;
%bot = wp + cp - 2*clp;
bot = ip - clp;
cint = 2;
FRAME = [30 120 -5 15];
figure(1)
sp(1)
mcont(top, [-20:cint:20], 'giso', [0 180]);
title(['CT: WARM - CLIM Precipitation'])
xlabel(['Contour Interval: ' num2str(cint) 'mm (day)^-^1'])
sp(2)
mcont(bot, [-20:cint:20], 'giso', [0 180]);
title(['CT: COLD - CLIM Precipitation'])
xlabel(['Contour Interval: ' num2str(cint) 'mm (day)^-^1'])
[ix, iy] = keep_var([60 105 0 10], lon, lat);
weight = mean(mean(bot(ix, iy))) / mean(mean(top(ix, iy)));
% Here's another:
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[wzt, wct] = getnc('wct.nc', 'ZMDT', 'CMFDT');
[izt, ict] = getnc('india.nc', 'ZMDT', 'CMFDT');
[clzt, clct] = getnc('clim.nc', 'ZMDT', 'CMFDT');
tem1 = (wzt + wct - clzt - clct);
tem1 = squeeze(mean(mean(mean(mean(tem1(:,8:18,iy,ix))))));
tem2 = (izt + ict - clzt - clct);
tem2 = squeeze(mean(mean(mean(mean(tem2(:,8:18,iy,ix))))));
weight = tem2 / tem1
% Check out SST's:
wts = getnc('wgr.nc', 'TS');
cts = getnc('cgr.nc', 'TS');
clts = getnc('clim.nc', 'TS');
top = squeeze(mean(wts-clts))/2;
bot = squeeze(mean(cts-clts))/2;
cint = 0.2;
default_global;
figure(1)
sp(1)
mcont(top, [-5:cint:5], 'giso', [0 180]);
title(['CT: WARM - CLIM Precipitation'])
xlabel(['Contour Interval: ' num2str(cint) 'mm (day)^-^1'])
sp(2)
mcont(bot, [-5:cint:5], 'giso', [0 180]);
title(['CT: COLD - CLIM Precipitation'])
xlabel(['Contour Interval: ' num2str(cint) 'mm (day)^-^1'])
filin = 'hgt.mon.mean.nc';
varn = 'hgt'
cd /home/disk/tao/dvimont/matlab/CCM/CT
[ppat, npat, pat, clim] = pn_reg_ct(filin, varn, 500);
cd /home/disk/tao/data/nmc.reanalysis/monthly
[lat, lon] = getll(filin);
if (strcmp(varn, 'prate'));
ppat = ppat*(24*3600);
npat = npat*(24*3600);
end
get_global
default_global
%FRAME = [0 360 -40 40];
cint = 10
figure(1)
sp(1)
% mcont(ppat, [-5:cint:5], 'giso', [0 180]);
mcont(ppat, [-500:cint:500]);
title(['NMC: WARM CT* Weighted Average 500mb HGT'])
xlabel(['Contour Interval ' num2str(cint) 'm'])
sp(2)
% mcont(npat, [-5:cint:5], 'giso', [0 180]);
mcont(npat, [-500:cint:500]);
title(['NMC: COLD CT* Weighted Average 500mb HGT'])
xlabel(['Contour Interval ' num2str(cint) 'm'])
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots/WA