Documentation of plot_PREC_for_SUN

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

Help text

 print -dps2 CT_OMEGA_PRECIP.ps

Cross-Reference Information

Listing of script plot_PREC_for_SUN

clear
cd /cdrom
[p1, p2] = getnc('wct.nc', 'PRECC', 'PRECL');
pw = (p1+p2)*24*3600*1000;
[p1, p2] = getnc('cct.nc', 'PRECC', 'PRECL');
pc = (p1+p2)*24*3600*1000;
[lat, lon, lev] = getll('wct.nc');
get_global
default_global
[xk, yk] = keep_var([50 290 -10 10], lon, lat);
nc = netcdf('wct.nc', 'nowrite');
  omw = nc{'OMEGA'}(:,:,yk,xk);
nc = close(nc);
nc = netcdf('cct.nc', 'nowrite');
  omc = nc{'OMEGA'}(:,:,yk,xk);
nc = close(nc);
figure(1)
top = squeeze(mean(omw-omc))/2;
top = squeeze(mean(shiftdim(top, 1)));
top = shiftdim(top, 1);
sp(1)
     pncont(lon(xk), -1*lev, top, [-.1:.005:.1], 'k')
     set(gca, 'XTick', [60:30:300])
     set(gca, 'XTickLabel', [' 60E';' 90E';'120E';'150E';' 180';'150W';'120W';' 90W'])
     set(gca, 'YTickLabel', [1000:-100:0])
     title('CT RUN:  OMEGA (10S:10N AVERAGED) -- (WARM - COLD)/2');
     xlabel('CONTOUR INTERVAL:  0.005 mb s^-^1')
     ylabel('PRESSURE (mb)')
     axis([50 290 -1000 0])
sp(2)
     bot = squeeze(mean(pw - pc)) / 2;
     FRAME = [50 290 -40 40];
     gcont(bot, [-15:15]);
     dc
     xlabel('CONTOUR INTERVAL:  1 mm (day)^-^1');
     title('CT RUN:  PRECIPITATION -- (WARM - COLD)/2');
cd /home/disk/tao/dvimont/public_html/research/images

cd /home/disk/hayes/dvimont/ccm3.6/data/CT
tsw = getnc('ct_warm_sst.nc', 'SST');
tsc = getnc('ct_cold_sst.nc', 'SST');

clf 
orient landscape
     bot = squeeze(mean(tsw - tsc)) / 2;
     FRAME = [0 360 -90 90];
     gcont(bot, [-15:0.2:15]);
     dc
     xlabel('CONTOUR INTERVAL:  0.2 K');
     title('CT RUN:  SSTA -- (WARM - COLD)/2');

cd /home/disk/tao/dvimont/public_html/research/images
% print -dps2 CT_SSTA.ps


cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/t31_ct/tpdata

[p1, p2, omwg] = getnc('wct_globsst_t31.nc', 'PRECC', 'PRECL', 'OMEGA');
pwg = (p1 + p2) * 24 * 3600 * 1000;
[p1, p2, omcg] = getnc('cct_globsst_t31.nc', 'PRECC', 'PRECL', 'OMEGA');
pcg = (p1 + p2) * 24 * 3600 * 1000;
[p1, p2, omwt] = getnc('wct_TropPac_t31.nc', 'PRECC', 'PRECL', 'OMEGA');
pwt = (p1 + p2) * 24 * 3600 * 1000;
[p1, p2, omct] = getnc('cct_TropPac_t31.nc', 'PRECC', 'PRECL', 'OMEGA');
pct = (p1 + p2) * 24 * 3600 * 1000;
[lat, lon, lev] = getll('cct_TropPac_t31.nc');
get_global
default_global

[xk, yk] = keep_var([50 290 -10 10], lon, lat);
top = squeeze(mean(omwt-omct))/2;
top = squeeze(mean(shiftdim(top(:,yk,xk), 1)));
top = shiftdim(top, 1);

figure(1)
sp(1)
     pncont(lon(xk), -1*lev, top, [-.1:.005:.1], 'k')
     set(gca, 'XTick', [60:30:300])
     set(gca, 'XTickLabel', [' 60E';' 90E';'120E';'150E';' 180';'150W';'120W';' 90W'])
     set(gca, 'YTickLabel', [1000:-100:0])
     title('TROPICAL SST CT RUN:  OMEGA (10S:10N AVERAGED) -- (WARM - COLD)/2');
     xlabel('CONTOUR INTERVAL:  0.005 mb s^-^1')
     ylabel('PRESSURE (mb)')
     axis([50 290 -1000 0])

sp(2)
     bot = squeeze(mean(pwt - pct)) / 2;
     FRAME = [50 290 -40 40];
     gcont(bot, [-15:15]);
     dc
     xlabel('CONTOUR INTERVAL:  1 mm (day)^-^1');
     title('TROPICAL CT RUN:  PRECIPITATION -- (WARM - COLD)/2');

cd /home/disk/tao/dvimont/public_html/research/images
% print -dps glob_ct_omega_prec.ps
% print -dps trop_ct_omega_prec.ps