Global Index (short | long) | Local contents | Local Index (short | long)
mquiv(top1, top2, 1, 'stereo', [90 270]);
| This script calls | |
|---|---|
clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[lat, lon, lev, hyam, hybm, P0] = getll('c_cycle.nc');
[whgt1, wprc1, wprl1, wu1, wv1, wps1, wpsl1] = ...
getnc('w_cycle.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[chgt1, cprc1, cprl1, cu1, cv1, cps1, cpsl1] = ...
getnc('c_cycle.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[whgt2, wprc2, wprl2, wu2, wv2, wps2, wpsl2] = ...
getnc('wgr_5093.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
[chgt2, cprc2, cprl2, cu2, cv2, cps2, cpsl2] = ...
getnc('cgr_5093.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
'PSL');
tim = [1 2 3 4 11 12];
whgt1 = squeeze(mean(whgt1));
chgt1 = squeeze(mean(chgt1));
whgt2 = squeeze(mean(whgt2(tim,:,:,:)));
chgt2 = squeeze(mean(chgt2(tim,:,:,:)));
wprc1 = squeeze(mean(wprc1));
cprc1 = squeeze(mean(cprc1));
wprc2 = squeeze(mean(wprc2(tim,:,:)));
cprc2 = squeeze(mean(cprc2(tim,:,:)));
wprl1 = squeeze(mean(wprl1));
cprl1 = squeeze(mean(cprl1));
wprl2 = squeeze(mean(wprl2(tim,:,:)));
cprl2 = squeeze(mean(cprl2(tim,:,:)));
wu1 = squeeze(mean(wu1));
cu1 = squeeze(mean(cu1));
wu2 = squeeze(mean(wu2(tim,:,:,:)));
cu2 = squeeze(mean(cu2(tim,:,:,:)));
wv1 = squeeze(mean(wv1));
cv1 = squeeze(mean(cv1));
wv2 = squeeze(mean(wv2(tim,:,:,:)));
cv2 = squeeze(mean(cv2(tim,:,:,:)));
wps1 = squeeze(mean(wps1));
cps1 = squeeze(mean(cps1));
wps2 = squeeze(mean(wps2(tim,:,:)));
cps2 = squeeze(mean(cps2(tim,:,:)));
wpsl1 = squeeze(mean(wpsl1));
cpsl1 = squeeze(mean(cpsl1));
wpsl2 = squeeze(mean(wpsl2(tim,:,:)));
cpsl2 = squeeze(mean(cpsl2(tim,:,:)));
lev2 = [100 150 250 500 850 925 970 1000];
whgt1 = atlev(whgt1, lev2, wps1, hyam, hybm, P0);
chgt1 = atlev(chgt1, lev2, cps1, hyam, hybm, P0);
whgt2 = atlev(whgt2, lev2, wps2, hyam, hybm, P0);
chgt2 = atlev(chgt2, lev2, cps2, hyam, hybm, P0);
get_global
default_global
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 -60 60];
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];
cnum = 10;
num = 18;
top = squeeze(wu1(num,:,:) - cu1(num,:,:)) / 2;
bot = squeeze(wu2(num,:,:) - cu2(num,:,:)) / 2;
units = ['M/S'];
head = [num2str(10 * round(lev(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(bot, cnum);
figure(1)
sp(1)
mcont(top, clevt, 'giso', [0 180]);
title([head ', GR\_TPAC run']);
xlabel(['Units: ' num2str(cintt) ' ' units]);
sp(2)
mcont(bot, clevb, 'giso', [0 180]);
title([head ', CT\_TPAC run']);
xlabel(['Units: ' num2str(cintb) ' ' units]);
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots
lev = 250;
cint = 20; clev = [-100:cint:100];
num= find(lev2==lev);
top = squeeze(whgt1(num,:,:) - chgt1(num,:,:)) / 2;
bot = squeeze(whgt2(num,:,:) - chgt2(num,:,:)) / 2;
units = ['M'];
head = [num2str(lev) 'mb Height: WARM - COLD'];
FRAME = [0 360 -90 90];
num = 11;
top = squeeze(whgt1(num,:,:) - whgt2(num,:,:)) / 2;
bot = squeeze(chgt1(num,:,:) - chgt2(num,:,:)) / 2;
units = ['M'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: GLOB - TPAC'];
FRAME = [0 360 -90 90];
num = 8;
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];
cintt = 2;
cintb = 2;
num = 10;
top = squeeze(wu1(num,:,:) - wu2(num,:,:)) / 2;
bot = squeeze(cu1(num,:,:) - cu2(num,:,:)) / 2;
units = ['M/S'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb Height: GLOB - TPAC'];
FRAME = [0 360 -90 90];
cintt = 20; clevt = [-140:cintt:140];
cintb = 10; clevb = [-150:cintb:150];
figure(1)
sp(1)
clma
mcont(top, clevt, 'stereo', [90 270]);
title([head ', CT_CYCLE']);
xlabel(['Units: ' num2str(cintt) ' ' units]);
sp(2)
clma
mcont(bot, clevb, 'stereo', [90 270]);
title([head ', GR_DOM']);
xlabel(['Units: ' num2str(cintb) ' ' units]);
cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[wuu1, wvv1, wuv1] = getnc('wgr_TropPac.nc', ...
'UU', 'VV', 'UV');
[cuu1, cvv1, cuv1] = getnc('cgr_TropPac.nc', ...
'UU', 'VV', 'UV');
[wuu2, wvv2, wuv2] = getnc('wct_TropPac_t31.nc', ...
'UU', 'VV', 'UV');
[cuu2, cvv2, cuv2] = getnc('cct_TropPac_t31.nc', ...
'UU', 'VV', 'UV');
wuu1 = squeeze(mean(wuu1));
wvv1 = squeeze(mean(wvv1));
wuv1 = squeeze(mean(wuv1));
cuu1 = squeeze(mean(cuu1));
cvv1 = squeeze(mean(cvv1));
cuv1 = squeeze(mean(cuv1));
wuu2 = squeeze(mean(wuu2));
wvv2 = squeeze(mean(wvv2));
wuv2 = squeeze(mean(wuv2));
cuu2 = squeeze(mean(cuu2));
cvv2 = squeeze(mean(cvv2));
cuv2 = squeeze(mean(cuv2));
weke1 = (wuu1 + wvv1 - wu1.^2 - wv1.^2);
ceke1 = (cuu1 + cvv1 - cu1.^2 - cv1.^2);
weke2 = (wuu2 + wvv2 - wu2.^2 - wv2.^2);
ceke2 = (cuu2 + cvv2 - cu2.^2 - cv2.^2);
num = 11;
top = squeeze(weke1(num,:,:) - ceke1(num,:,:)) / 2;
bot = squeeze(weke2(num,:,:) - ceke2(num,:,:)) / 2;
units = ['M^2/S^2'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb EKE: WARM - COLD'];
FRAME = [0 360 -90 90];
cnum = 12;
[clevt, cintt] = get_contour_int(top, cnum);
[clevb, cintb] = get_contour_int(bot, cnum);
figure(1)
sp(1)
XAX = lon1;
YAX = lat1;
mcont(top, clevt, 'stereo', [90 270]);
title([head ', GLOBAL SST']);
xlabel(['Units: ' num2str(cintt) ' ' units]);
sp(2)
XAX = lon2;
YAX = lat2;
mcont(bot, clevb, 'stereo', [90 270]);
title([head ', TPAC SST']);
xlabel(['Units: ' num2str(cintb) ' ' units]);
wep1_1 = wvv1 - wv1.^2 - wuu1 + wu1.^2;
wep1_2 = wu1.*wv1-wuv1;
cep1_1 = cvv1 - cv1.^2 - cuu1 + cu1.^2;
cep1_2 = cu1.*cv1-cuv1;
wep2_1 = wvv2 - wv2.^2 - wuu2 + wu2.^2;
wep2_2 = wu2.*wv2-wuv2;
cep2_1 = cvv2 - cv2.^2 - cuu2 + cu2.^2;
cep2_2 = cu2.*cv2-cuv2;
clf
num = 11;
top1 = squeeze(wep1_1(num,:,:)-cep1_1(num,:,:))/2;
top2 = squeeze(wep1_2(num,:,:)-cep1_2(num,:,:))/2;
bot1 = squeeze(wep2_1(num,:,:)-cep2_1(num,:,:))/2;
bot2 = squeeze(wep2_2(num,:,:)-cep2_2(num,:,:))/2;
units = ['M^2/S^2'];
head = [num2str(10 * round(lev1(num) / 10)) 'mb E-VECT: WARM - COLD'];
FRAME = [120 270 0 90];
num = 8;
top1 = squeeze(wep1_1(num,:,:));
top2 = squeeze(wep1_2(num,:,:));
bot1 = squeeze(cep1_1(num,:,:));
bot2 = squeeze(cep1_2(num,:,:));
units = ['M^2/S^2'];
head1 = [num2str(10 * round(lev(num) / 10)) 'mb E-VECT: WARM'];
head2 = [num2str(10 * round(lev(num) / 10)) 'mb E-VECT: COLD'];
FRAME = [0 150 -20 80];
figure(1)
sp(1)
gquiv(top1, top2, -0.25);
title([head1 ', TPAC SST']);
xlabel(['Units: ' units]);
dc
sp(2)
gquiv(bot1, bot2, -0.25);
% mquiv(bot1, bot2, 1, 'stereo', [90 270]);
title([head2 ', TPAC SST']);
xlabel(['Units: ' units]);
dc