Documentation of plot_hgt

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

Cross-Reference Information

Listing of script plot_hgt

cd /cdrom
filin = 'clim_dec82-jun98.nc';
nc = netcdf(filin, 'nowrite');
  hgt = nc{'Z3'}(:,11,:,:);
  ps  = nc{'PS'}(:);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
nc = close(nc);
filin = 'warmenso.nc';
nc = netcdf(filin, 'nowrite');
  whgt = nc{'Z3'}(:,11,:,:);
  wps  = nc{'PS'}(:);
nc = close(nc);
filin = 'coldenso.nc';
nc = netcdf(filin, 'nowrite');
  chgt = nc{'Z3'}(:,11,:,:);
  cps  = nc{'PS'}(:);
nc = close(nc);
hgt = squeeze(mean(hgt([1 2 12],:,:,:)));
whgt = squeeze(mean(whgt(2:4,:,:,:)));
chgt = squeeze(mean(chgt(2:4,:,:,:)));
cd ~/matlab/CCM/oi_plots
orient tall
global XAX YAX FRAME
FRAME = [0 360 -90 90];
XAX = lon;
YAX = lat;
tem = whgt - hgt;
subplot('Position', [.1 .53 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('OI DJF Warm Event 500mb heights (WARM - CLIM)')
     xlabel('Contour Interval:  20m')
tem = chgt - hgt;
subplot('Position', [.1 .0525 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('OI DJF Cold Event 500mb heights (COLD - CLIM)')
     xlabel('Contour Interval:  20m')
tem = whgt - chgt;
subplot('Position', [.1 .53 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('OI DJF 500mb heights (WARM - COLD)')
     xlabel('Contour Interval:  20m')
tem = chgt + whgt - 2*hgt;
subplot('Position', [.1 .0525 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('OI DJF 500mb heights (WARM + COLD - 2*CLIM)')
     xlabel('Contour Interval:  20m')
clear
cd /home/disk/tao/data/nmc.reanalysis/monthly
filin = 'hgt.mon.mean.nc';
nc = netcdf(filin, 'nowrite');
  hgt = nc{'hgt'}(289:486,11,:,:);
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  time = nc{'time'}(289:486);
nc = close(nc);
hgt = squeeze(hgt);
[ntim, nlat, nlon] = size(hgt);
tem = (time - time(1))/(24);
mo = [];
for i = 1:(ceil(ntim/12));
  mo = [mo 1:12];
end
mo = mo(1:ntim);
yr = floor(tem / 365) + 1982;
wyr = [1983 1988 1992 1995 1998];
cyr = [1984 1985 1989 1996 1997];
wind = intersect(find(ismember(yr, wyr)), find(ismember(mo, 1)));
cind = intersect(find(ismember(yr, cyr)), find(ismember(mo, 1)));
whgt = reshape(mean2(hgt([wind wind-1 wind+1],:,:)),nlat,nlon);
chgt = reshape(mean2(hgt([cind cind-1 cind+1],:,:)),nlat,nlon);
[anom, hgt] = annave(hgt);
hgt = reshape(hgt, 12, nlat, nlon);
hgt = squeeze(mean(hgt([1 2 12],:,:)));
global XAX YAX FRAME
cd ~/matlab/CCM/oi_plots
FRAME = [0 360 -90 90];
XAX = lon;
YAX = lat;
tem = whgt - hgt;
subplot('Position', [.1 .53 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('NMC DJF Warm Event 500mb heights (WARM - CLIM)')
     xlabel('Contour Interval:  20m')
tem = chgt - hgt;
subplot('Position', [.1 .0525 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('NMC DJF Cold Event 500mb heights (COLD - CLIM)')
     xlabel('Contour Interval:  20m')
tem = whgt - chgt;
subplot('Position', [.1 .53 .8 .42])
     mcont(tem, [-210:40:210], 1);
     title('NMC DJF 500mb heights (WARM - COLD)')
     xlabel('Contour Interval:  40m')
tem = chgt + whgt - 2*hgt;
subplot('Position', [.1 .0525 .8 .42])
     mcont(tem, [-110:20:110], 1);
     title('NMC DJF 500mb heights (WARM + COLD - 2*CLIM)')
     xlabel('Contour Interval:  20m')