Documentation of plot_tpgr

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

Help text

  top = squeeze(mean(tem1(:,ind,:,:)-tem3(:,ind,:,:)));

Cross-Reference Information

Listing of script plot_tpgr

clear
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
var = 'Z3'; var2 = 'PS';
var = 'PRECL'; var2 = 'PRECC';
var = 'U'; var2 = 'V';
var = 'TS';
var = 'PSL';
[tem1, ps1] = getnc('wgr_SOM_hard.nc', var, var2);
[tem2, ps2] = getnc('cgr_SOM_hard.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('clim_t31x15_ense.nc');
for i = 1:2;
  eval(['ps' num2str(i) ' = ps' num2str(i) '([1 2 3 4 11 12], :, :);']);
  eval(['tem' num2str(i) ' = tem' num2str(i) '([1 2 3 4 11 12], :, :, :);']);
end
if ndims(tem1) == 4;
  klev = [5:size(tem1,2)];
  newlevs = [100 150 250 500 850 1000];
  tem1 = atlev(tem1(:,klev,:,:), newlevs, ps1, hyam(klev), hybm(klev), P0);
  tem2 = atlev(tem2(:,klev,:,:), newlevs, ps2, hyam(klev), hybm(klev), P0);
  tem3 = atlev(tem3(:,klev,:,:), newlevs, ps3, hyam(klev), hybm(klev), P0);
end
get_global
default_global
cd /home/disk/tao/dvimont/matlab/CCM/GR
if ndims(tem1) == 4;
  mblev = 500;
  cint = 5;
  filin = 'hgt.mon.mean.nc'; nmcvar = 'hgt';
  ind = find(newlevs == mblev);
  [top, clim] = reg_gr(filin, nmcvar, mblev);
  [lat1, lon1] = getll(filin);
  bot = squeeze(mean(tem1(:,ind,:,:)-tem2(:,ind,:,:)))/2;
  ttl = [num2str(mblev) 'mb HGT:  NMC Regressed on GR'];
  btl = [num2str(mblev) 'mb HGT:  CCM GR Run, (WARM - COLD) / 2'];
else
  top = squeeze(mean(tem1 - tem2));
  bot = squeeze(mean(tem2 - tem3));
  ttl = ['Precipitation'];
  xl = ['mm / day'];
end

figure(2)
sp(1)
     XAX = lon1;
     YAX = lat1;
     mcont(top, [-100:cint:100]);
     title('NMC:  NDJFMAM 500mb Height regressed on GR')
%     title(['GR\_TPAC, ' ttl]);
     xlabel(['Contour Interval:  ' num2str(cint) 'm']);
sp(2)
     XAX = lon;
     YAX = lat;
     mcont(bot, [-100:cint:100]);
     title('CCM:  GR\_SOM\_Hard, NDJFMA 500mb Height, (WARM - COLD)/2')
%     title(['GR\_TPAC, ' btl]);
     xlabel(['Contour Interval:  ' num2str(cint) 'm']);

cd /home/disk/tao/dvimont/pub/david
cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots

ttl = ['Precipitation'];
xl = ['mm / day'];
ttl = ['SLP'];
xl = ['mb'];
cint = 0.5;
figure(1)
     FRAME = [0 360 -90 90];
     XAX = lon;
     YAX = lat;
sp(1)
     mcont(top, [-5:cint:5], 'giso', [0 180])
     title([ttl ', (WARM - CLIM)']);
     xlabel(['Units:  ' num2str(cint) ' ' xl]);
sp(2)
     mcont(bot, [-5:cint:5], 'giso', [0 180]);
     title([ttl ', (COLD - CLIM)']);
     xlabel(['Units:  ' num2str(cint) ' ' xl]);

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('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