Documentation of pac_vs_global

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

     mcont(bot, clevb, 'giso', [0 180]);

Cross-Reference Information

Listing of script pac_vs_global

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/t31_ct/data
[hgt1, prc1, prl1, u1, v1, ps1, psl1] = ...
     getnc('wct_globsst_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
	   'PSL');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/t31_ct/tpdata
[hgt2, prc2, prl2, u2, v2, ps2, psl2] = ...
     getnc('wct_TropPac_t31.nc', 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
	   'PSL');
[lat, lon, lev] = getll('wct_TropPac_t31.nc');
hgt1 = squeeze(mean(hgt1));
hgt2 = squeeze(mean(hgt2));
prc1 = squeeze(mean(prc1));
prc2 = squeeze(mean(prc2));
prl1 = squeeze(mean(prl1));
prl2 = squeeze(mean(prl2));
u1 = squeeze(mean(u1));
u2 = squeeze(mean(u2));
v1 = squeeze(mean(v1));
v2 = squeeze(mean(v2));
ps1 = squeeze(mean(ps1));
ps2 = squeeze(mean(ps2));
psl1 = squeeze(mean(psl1));
psl2 = squeeze(mean(psl2));
get_global
default_global
FRAME = [0 360 -90 90];
num = 18;
cnum = 10
top = squeeze((prc1 + prl1) - (prc2 + prl2)) * (3600 * 24 * 1000) / 2;
bot1 = squeeze(u1(num,:,:) - u2(num,:,:));
bot2 = squeeze(v1(num,:,:) - v2(num,:,:));
units1 = ['mm (day)^-^1'];
head1 = ['Precipitation:  GLOB-TPAC'];
units2 = ['M/S'];
head2 = [num2str(10 * round(lev(num) / 10)) 'mb Winds:  GLOB-TPAC'];
FRAME = [0 360 -40 40];
[clevt, cintt] = get_contour_int(top, cnum);
figure(1)
sp(1)
     gcont(top, clevt);
     title([head1 ]);
     xlabel(['Units:  ' num2str(cintt) ' ' units1]);
     dc
sp(2)
     gquiv(bot1, bot2, 2)
     title([head2 ]);
     xlabel(['Units:  ' num2str(cintb) ' ' units2]);
     dc

num = 11;
cnum = 12;
top = squeeze(hgt1(num,:,:) - hgt2(num,:,:));
bot = squeeze(u1(num,:,:) - u2(num,:,:));
head1 = [num2str(10 * round(lev(num) / 10)) 'mb HGT:  GLOB-TPAC'];
head2 = [num2str(10 * round(lev(num) / 10)) 'mb Zonal Wind:  GLOB-TPAC'];
units1 = ['M'];
units2 = ['M/S'];
FRAME = [0 360 -90 90];

top = (prc1 + prl1) * (3600 * 24 * 1000) / 2;
bot = (prc2 + prl2) * (3600 * 24 * 1000) / 2;
cnum = 8.5;
head1 = ['Precipitation, GLOB'];
head2 = ['Precipitation, TPAC'];
units1 = ['mm/day'];
units2 = units1;
FRAME = [0 360 -90 90];

[clevt, cintt] = get_contour_int(top, cnum);
[clevb, cintb] = get_contour_int(bot, cnum);
figure(1)
sp(1)
     gcont(top, clevt);
     title([head1]);
     xlabel(['Units:  ' num2str(cintt) ' ' units1])
     dc
sp(2)
     gcont(bot, clevb);
     title([head2]);
     xlabel(['Units:  ' num2str(cintb) ' ' units2])
     dc