Global Index (short | long) | Local contents | Local Index (short | long)
[clevt, cintt] = get_contour_int(top, 12); [clevb, cintb] = get_contour_int(top, 12);
| This script calls | |
|---|---|
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[lat1, lon1, lev1, hyam, hybm, P0] = getll('wct.nc');
[lat2, lon2, lev2] = getll('wct_globsst_t31.nc');
[whgt1, wprc1, wprl1, wu1, wv1, wps1, wpsl1] = ...
getnc('wct.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[chgt1, cprc1, cprl1, cu1, cv1, cps1, cpsl1] = ...
getnc('cct.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[whgt2, wprc2, wprl2, wu2, wv2, wps2, wpsl2] = ...
getnc('wct_globsst_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[chgt2, cprc2, cprl2, cu2, cv2, cps2, cpsl2] = ...
getnc('cct_globsst_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
whgt1 = squeeze(mean(whgt1));
chgt1 = squeeze(mean(chgt1));
whgt2 = squeeze(mean(whgt2));
chgt2 = squeeze(mean(chgt2));
wprc1 = squeeze(mean(wprc1));
cprc1 = squeeze(mean(cprc1));
wprc2 = squeeze(mean(wprc2));
cprc2 = squeeze(mean(cprc2));
wprl1 = squeeze(mean(wprl1));
cprl1 = squeeze(mean(cprl1));
wprl2 = squeeze(mean(wprl2));
cprl2 = squeeze(mean(cprl2));
wu1 = squeeze(mean(wu1));
cu1 = squeeze(mean(cu1));
wu2 = squeeze(mean(wu2));
cu2 = squeeze(mean(cu2));
wv1 = squeeze(mean(wv1));
cv1 = squeeze(mean(cv1));
wv2 = squeeze(mean(wv2));
cv2 = squeeze(mean(cv2));
wps1 = squeeze(mean(wps1));
cps1 = squeeze(mean(cps1));
wps2 = squeeze(mean(wps2));
cps2 = squeeze(mean(cps2));
wpsl1 = squeeze(mean(wpsl1));
cpsl1 = squeeze(mean(cpsl1));
wpsl2 = squeeze(mean(wpsl2));
cpsl2 = squeeze(mean(cpsl2));
get_global
FRAME = [0 360 -90 90];
top = squeeze((wprc1 + wprl1) - (cprc1 + cprl1)) * (3600 * 24 * 1000) / 2;
bot = squeeze((wprc2 + wprl2) - (cprc2 + cprl2)) * (3600 * 24 * 1000) / 2;
units = ['mm (day)^-^1'];
head = ['Precipitation: WARM - COLD'];
FRAME = [0 360 -40 40];
cnum = 12;
top = squeeze(wps1 - cps1) / 200;
bot = squeeze(wps2 - cps2) / 200;
units = ['mb'];
head = ['Surface Pressure: WARM - COLD'];
FRAME = [0 360 -90 90];
top = squeeze(wpsl1 - cpsl1) / 200;
bot = squeeze(wpsl2 - cpsl2) / 200;
units = ['mb'];
head = ['SLP: WARM - COLD'];
FRAME = [0 360 -90 90];
num = 18;
top = squeeze(wu1(num,:,:) - cu1(num,:,:)) / 2;
bot = squeeze(wu2(num,:,:) - cu2(num,:,:)) / 2;
units = ['M/S'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb Zonal Wind: WARM - COLD'];
FRAME = [0 360 -90 90];
cnum = 12;
num = 17;
divw1 = sph_div(squeeze(wu1(num,:,:)), squeeze(wv1(num,:,:)), lon1, lat1);
divc1 = sph_div(squeeze(cu1(num,:,:)), squeeze(cv1(num,:,:)), lon1, lat1);
divw2 = sph_div(squeeze(wu2(num,:,:)), squeeze(wv2(num,:,:)), lon2, lat2);
divc2 = sph_div(squeeze(cu2(num,:,:)), squeeze(cv2(num,:,:)), lon2, lat2);
top = (divw1 - divc1) / 2;
bot = (divw2 - divc2) / 2;
head = [num2str(10 * round(lev1(num) / 10)) 'mb Zonal Wind: WARM - COLD'];
units = ['S^-^1'];
cnum = 10;
FRAME = [0 360 -40 40];
num = 8;
curlw1 = sph_curl(squeeze(wu1(num,:,:)), squeeze(wv1(num,:,:)), lon1, lat1);
curlc1 = sph_curl(squeeze(cu1(num,:,:)), squeeze(cv1(num,:,:)), lon1, lat1);
curlw2 = sph_curl(squeeze(wu2(num,:,:)), squeeze(wv2(num,:,:)), lon2, lat2);
curlc2 = sph_curl(squeeze(cu2(num,:,:)), squeeze(cv2(num,:,:)), lon2, lat2);
top = (curlw1 - curlc1) / 2;
bot = (curlw2 - curlc2) / 2;
head = [num2str(10 * round(lev1(num) / 10)) 'mb Relative Vort.: WARM - COLD'];
units = ['S^-^1'];
cnum = 12;
FRAME = [0 360 -90 90];
[clevt, cintt] = get_contour_int(top, cnum);
[clevb, cintb] = get_contour_int(top, cnum);
figure(1)
sp(1)
XAX = lon1;
YAX = lat1;
mcont(top, clevt, 'giso', [0 180]);
title(['CT\_run, ' head ', T42 version']);
xlabel(['Units: ' num2str(cintt) ' ' units]);
sp(2)
XAX = lon2;
YAX = lat2;
mcont(bot, clevb, 'giso', [0 180]);
title(['CT\_run, ' head ', t31x15 version']);
xlabel(['Units: ' num2str(cintb) ' ' units]);
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots
num = 8;
top = squeeze(whgt1(num,:,:) - chgt1(num,:,:)) / 2;
bot = squeeze(whgt2(num,:,:) - chgt2(num,:,:)) / 2;
units = ['M'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: WARM - COLD'];
FRAME = [0 360 -90 90];
num = 11;
top = squeeze(wu1(num,:,:) - cu1(num,:,:)) / 2;
bot = squeeze(wu2(num,:,:) - cu2(num,:,:)) / 2;
units = ['M/S'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: WARM - COLD'];
FRAME = [0 360 -90 90];
top1 = atlev(whgt1, [250 500], wps1, hyam, hybm, P0);
top2 = atlev(chgt1, [250 500], cps1, hyam, hybm, P0);
bot1 = atlev(whgt2, [250 500], wps2, hyam, hybm, P0);
bot2 = atlev(chgt2, [250 500], cps2, hyam, hybm, P0);
top = top1 - top2;
bot = bot1 - bot2;
top = squeeze(top(1,:,:))/2;
bot = squeeze(bot(1,:,:))/2;
units = ['M'];
head = ['250 mb Height: WARM - COLD'];
FRAME = [0 360 -90 90];
cintt = 10;
clevt = [-500:cintt:500];
cintb = cintt;
clevb = clevt;
figure(1)
sp(1)
XAX = lon1;
YAX = lat1;
mcont(top, clevt, 'stereo', [90 270]);
title([head ', T42 version']);
xlabel(['Units: ' num2str(cintt) ' ' units]);
sp(2)
XAX = lon2;
YAX = lat2;
mcont(bot, clevb, 'stereo', [90 270]);
title([head ', t31x15 version']);
xlabel(['Units: ' num2str(cintb) ' ' units]);
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots