Documentation of plot_TPac_vs_NMC

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

Help text

clhgt = atlev(hgt3, lev1, ps2, hyam, hybm, P0);

Cross-Reference Information

Listing of script plot_TPac_vs_NMC

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
filtag = '_ense.nc'
[hgt1, prc1, prl1, u1, v1, ps1, psl1] = ...
     getnc(['wct' filtag], 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
	   'PSL');
[hgt2, prc2, prl2, u2, v2, ps2, psl2] = ...
     getnc(['cct' filtag], 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
	   'PSL');
[hgt3, prc3, prl3, u3, v3, ps3, psl3] = ...
     getnc(['clim' filtag], 'Z3', 'PRECC', 'PRECL', 'U', 'V', 'PS', ...
	   'PSL');
[lat, lon, lev, hyam, hybm, P0] = getll('wct_TropPac_t31.nc');
[ntim, nlev, nlat, nlon] = size(hgt3);
cd /home/disk/tao/data/nmc.reanalysis/monthly
[lat1, lon1, lev1] = getnc('hgt.mon.mean.nc', 'lat', 'lon', 'level');
cd /home/disk/tao/dvimont/matlab/CCM/CT
lev1 = [1000 850 500 250 100];
klev = [5:nlev];
whgt = atlev(hgt1(:,klev,:,:), lev1, ps1, hyam(klev), hybm(klev), P0);
chgt = atlev(hgt2(:,klev,:,:), lev1, ps2, hyam(klev), hybm(klev), P0);

[nmc_hgt500, clim] = reg_ct('hgt.mon.mean.nc', 'hgt', 500);
[nmc_hgt250, clim] = reg_ct('hgt.mon.mean.nc', 'hgt', 250);

get_global
FRAME = [0 360 -90 90];

lind = find(lev1 == 250);
top = nmc_hgt250;
bot = squeeze(mean(whgt(:,lind,:,:)-chgt(:,lind,:,:)))/2.8;
cint = 10;
hem = 1;

figure(1)
sp(1)
     XAX = lon1;
     YAX = lat1;
     mcont(top, [-90:cint:90], 'stereo', [hem*90 270]);
     title(['NH NMC:  ' num2str(lev1(lind)) 'mb Height regressed on CT*']);
     xlabel(['Contour Interval = ' num2str(cint) 'm'])
sp(2)
     XAX = lon;
     YAX = lat;
     mcont(bot, [-90:cint:90], 'stereo', [hem*90 270]);
     title(['NH CCM\_TPAC:  ' num2str(lev1(lind)) 'mb Height, WARM-COLD']);
     xlabel(['Contour Interval = ' num2str(cint) 'm'])

cd /home/disk/tao/dvimont/pub/david
%cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots

%cd /home/disk/tao/dvimont/matlab/CCM/CT
%[prec_nmc, clim] = reg_ct('prate.mon.mean.nc', 'prate', 0);
cd /home/disk/tao/dvimont/matlab/CCM/CT/Vars
load NMC_Prec_regr_on_CTstar.mat
cd /home/disk/tao/data/nmc.reanalysis/monthly
[lat1, lon1] = getnc('prate.mon.mean.nc', 'lat', 'lon');

top = prec_nmc * 3600 * 24;
bot = squeeze(mean(prl1 + prc1 - prl2 - prc2)) * (24*3600*1000)/2.8;
cint = 1;
FRAME = [0 360 -60 60];

figure(1)
sp(1)
     XAX = lon1;
     YAX = lat1;
     mcont(top, [-8:cint:8], 'giso', [0 180]);
     title(['NMC:  PRECIP regressed on CT*']);
     xlabel(['Contour Interval = ' num2str(cint) 'mm (day)^-^1'])
sp(2)
     XAX = lon;
     YAX = lat;
     mcont(bot, [-80:cint:80], 'giso', [0 180]);
     title(['CCM\_TPAC:  PRECIP, WARM-COLD']);
     xlabel(['Contour Interval = ' num2str(cint) 'mm (day)^-^1'])

cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots

cd /home/disk/tao/dvimont/matlab/CCM/CT
%filin = 'rain.monthly.7998.nc';
%[prec_merged, prec_clim_merged] = reg_ct(filin, 'data', 0);
filin = 'msu_prec.monthly.7997.orig.nc';
[prec_msu, prec_clim_msu] = reg_ct(filin, 'data', 0);
%cd /home/disk/tao/data/merged_precip
cd /home/disk/tao/data/MSU/monthly
[lat1, lon1] = getnc(filin, 'lat', 'lon');

cd /home/disk/tao/dvimont/matlab/CCM/CT

top = prec_msu/10;
bot = squeeze(mean(prl1 + prc1 - prl2 - prc2)) * (24*3600*1000)/2.8;
cint = 1;
FRAME = [0 360 -60 60];

figure(1)
sp(1)
     XAX = lon1;
     YAX = lat1;
     mcont(top, [-8:cint:8], 'giso', [0 180]);
     title(['MERGED PRECIP regressed on CT* (1979 to 1997)']);
     xlabel(['Contour Interval = ' num2str(cint) 'mm (day)^-^1'])
sp(2)
     XAX = lon;
     YAX = lat;
     mcont(bot, [-80:cint:80], 'giso', [0 180]);
     title(['CCM\_TPAC:  PRECIP, WARM-COLD']);
     xlabel(['Contour Interval = ' num2str(cint) 'mm (day)^-^1'])

cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots

cd /home/disk/tao/data/nmc.reanalysis/monthly
[lat1, lon1] = getnc('slp.mon.mean.nc', 'lat', 'lon');

cd /home/disk/tao/dvimont/matlab/CCM/CT
[slp_nmc, clim] = reg_ct('slp.mon.mean.nc', 'slp', 0);

top = slp_nmc;
bot = squeeze(mean(psl1 - psl2))/280;
cint = 0.5;
FRAME = [0 360 -90 90];

figure(1)
sp(1)
     XAX = lon1;
     YAX = lat1;
     mcont(top, [-8:cint:8], 'giso', [0 180]);
     title(['NMC: SLP regressed on CT*']);
     xlabel(['Contour Interval = ' num2str(cint) 'mb'])
sp(2)
     XAX = lon;
     YAX = lat;
     mcont(bot, [-80:cint:80], 'giso', [0 180]);
     title(['CCM\_TPAC:  SLP, WARM-COLD']);
     xlabel(['Contour Interval = ' num2str(cint) 'mb'])

cd /home/disk/tao/dvimont/matlab/CCM/CT/CT_Plots