Documentation of ct_lpct_regs

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

Help text

  Load Data

Cross-Reference Information

Listing of script ct_lpct_regs

clean
cd /home/disk/tao/dvimont/Papers/obs/matlab/Data
load cti2.mat

ct = detrend(ct);

ctlp = rave(rave(ct, 25), 37);
cthp = ct - ctlp;

%  Define limits
lims = [0 360 -75 90];
tim = 1:639;
tim1 = [11:15]; tim2 = tim1+5-12; tim3 = tim1+11-12;

%  Start with SLP
cd /home/disk/tao/data/nmc.reanalysis/monthly
%filin = 'skt.mon.mean.nc';
%varn = 'skt';
%filin = 'slp.mon.mean.nc';
%varn = 'slp';
filin = 'hgt.mon.mean.nc';
varn = 'hgt'; cint = 2;

[dat, lat1, lon1] = getnc2(filin, varn, lims, 850, tim);
dat = squeeze(dat);
dat = annave(dat);

[ntim, nlat, nlon] = size(dat); nyr = floor(ntim/12);

dat = detrend(dat);

datlp = rave(rave(dat, 25), 37);

[pat1lp, c1lp] = regress_eof(datlp, ctlp, 0);
[pat1hp, c1hp] = regress_eof(dat-datlp, cthp, 0);
[pat1, c1] = regress_eof(dat, ct, 0);

%  Plot regressions
figure_tall(2); clf;
global_axes(6, 7/3, 0, .4, 1.5);

subplot2(1,1); cla;
global_latlon(lat1, lon1, lims);
  map_axis('giso');
  [c, h] = map_contour_pn(pat1, cint);
  tem = c1; tem(abs(tem)<0.2) = NaN;
  hold on;
    h2 = map_surface(tem, -0.5);
  hold off;
  caxis([-1 1]);
  fill_landmap('under', 0.8*[1 1 1]);
  gridm on; tightmap
  title2('Unfiltered CTI; Lag 0');
  set(gca, 'fontsize', 10);
  
subplot2(1,2); cla;
global_latlon(lat1, lon1, lims);
  map_axis('giso');
  [c, h] = map_contour_pn(pat1hp, cint);
  tem = c1hp; tem(abs(tem)<0.2) = NaN;
  hold on;
    h2 = map_surface(tem, -0.5);
  hold off;
  caxis([-1 1]);
  fill_landmap('under', 0.8*[1 1 1]);
  gridm on; tightmap
  title2('Highpass-filtered CTI; Lag 0');
  set(gca, 'fontsize', 10);

subplot2(1,3); cla;
global_latlon(lat1, lon1, lims);
  map_axis('giso');
  [c, h] = map_contour_pn(pat1lp, cint/2);
  tem = c1lp; tem(abs(tem)<0.2) = NaN;
  hold on;
    h2 = map_surface(tem, -0.5);
  hold off;
  caxis([-1 1]);
  fill_landmap('under', 0.8*[1 1 1]);
  gridm on; tightmap
  title2('Lowpass-filtered CTI; Lag 0');
  set(gca, 'fontsize', 10);

global_axes(4.5, .2, 0, .4, 8.2+1.5);
cb = subplot2(1,1);
  colorbar(cb);
  set(cb, 'XTick', 1:2:11, 'XTickLabel', [-1:.4:1]);
  xlabel('Correlation', 'fontsize', 9);
  

cd /home/disk/tao/dvimont/pub/david
%print -dpsc2 CTI_SLP_regressions.ps
%print -dpsc2 CTI_SKT_regressions.ps
%print -dpsc2 CTI_Z850_regressions.ps