Documentation of plot_wpgr

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

Cross-Reference Information

Listing of script plot_wpgr

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
var = 'Z3'; var2 = 'PS';
var = 'PRECL'; var2 = 'PRECC';
var = 'U'; var2 = 'V';
var = 'TS';
[tem1, ps1] = getnc('wwp_TropPac.nc', var, var2);
[tem2, ps2] = getnc('cwp_TropPac.nc', var, var2);
[tem3, ps3] = getnc('clim_t31x15_ense.nc', var, var2);
if var == 'PRECL';
  tem1 = (tem1 + ps1) * (3600 * 24 * 1000);
  tem2 = (tem2 + ps2) * (3600 * 24 * 1000);
  tem3 = (tem3 + ps3) * (3600 * 24 * 1000);
end
[lat, lon, lev, hyam, hybm, P0] = getll('wgr_TropPac.nc');
if ndims(tem1) == 4;
  newlevs = [100 150 250 500 850 925 1000];
  tem1 = atlev(tem1, newlevs, ps1, hyam, hybm, P0);
  tem2 = atlev(tem2, newlevs, ps2, hyam, hybm, P0);
  tem3 = atlev(tem3, newlevs, ps3, hyam, hybm, P0);
end
get_global
default_global
cd /home/disk/tao/dvimont/matlab/CCM/GR
if ndims(tem1) == 4;
  mblev = 250;
  cint = 10;
  ind = find(newlevs == 500);
  ind2 = find(newlevs == 250);
  top = squeeze(mean(tem1(:,ind2,:,:)-tem2(:,ind2,:,:)))/2;
  bot = squeeze(mean(tem1(:,ind,:,:)-tem2(:,ind,:,:)))/2;
  ttl = ['250mb HGT:  (WARM - COLD)'];
  btl = ['500mb HGT:  (WARM - COLD)'];
else
  top = squeeze(mean(tem1 - tem2)) / 2;
  bot = squeeze(mean(tem2 - tem3));
  ttl = ['Precipitation'];
  xl = ['mm / day'];
end
figure(1)
sp(1)
     mcont(top, [-500:cint:500]);
     title(['WP\_run, ' ttl]);
     xlabel(['Contour Interval:  ' num2str(cint) 'm']);
sp(2)
     mcont(bot, [-500:cint:500]);
     title(['WP\_run, ' btl]);
     xlabel(['Contour Interval:  ' num2str(cint) 'm']);
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots
ttl = ['Precipitation'];
xl = ['mm / day'];
cint = 1;
figure(1)
     FRAME = [0 360 -50 50];
     XAX = lon;
     YAX = lat;
sp(1)
     mcont(top, [-5:cint:5], 'giso', [0 180])
     title(['WP\_run:  ' ttl ', (WARM - COLD) / 2']);
     xlabel(['Units:  ' num2str(cint) ' ' xl]);
sp(2)
     mcont(bot, [-5:.1:5], 'giso', [0 180]);
     title(['WP\_run:  T_s, (WARM - COLD) / 2']);
     xlabel(['Units:  ' num2str(0.1) ' K']);
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots
ind = 18;
topu = squeeze(mean(tem1(:, ind, :, :) - tem2(:, ind, :, :))) / 2;
topv = squeeze(mean(ps1(:, ind, :, :) - ps2(:, ind, :, :))) / 2;
top2 = squeeze(mean(tem1 - tem2)) / 2;
bot = squeeze(mean(tem1 - tem2)) / 2;
figure(1)
XAX = lon;
YAX = lat;
FRAME = [40 290 -45 60];
sp(1)
     mquiv(topu, topv, 1.5, 'giso', [0 165]);
     title('WP\_run:  Surface Winds, (WARM - COLD) / 2');
sp(2)
     mcont(bot, [-5:.5:5], 'giso', [0 165]);
     title('Precipitation, (WARM - COLD) / 2');
     xlabel('Contour Interval:  0.5 mm / day');
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots