Documentation of regress_hp_lp_ct

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

Help text

  cut down on the size of the data

Cross-Reference Information

Listing of script regress_hp_lp_ct

clear
cd /home/disk/hayes2/dvimont/ocean/centralpacific/stoc
filin = 'sal.cdf';
nc = netcdf(filin, 'nowrite');
  ct = nc{'h'}(:,1,50:53,52:102);
  h = nc{'h'}(:,1,1:102,1:102);
end

ct = squeeze(mean2(mean2(shiftdim(ct, 2))));
h = squeeze(h);
h = (h(:,1:2:101,1:2:101) + h(:,2:2:102,1:2:101) + ...
     h(:,1:2:101,2:2:102) + h(:,2:2:102,2:2:102))/4;

%save h_ct_small.mat ct h

h2 = detrend(h);
ct2 = detrend(ct);

ct2 = rave(ct2, 3) - rave(ct2, 40);
reg = regress_eof(h2, ct2, 0);

figure(1); fl(1); clf;
  pncont(1:2:101, 1:2:101, squeeze(reg), [-3:.5:5], 0, 'k')
  axis([0 102 0 102]);

ct2 = detrend(ct);
ct2 = rave(ct2, 40);
reg = regress_eof(h2, ct2, 10);

figure(2); fl(1); clf;
  pncont(1:2:101, 1:2:101, squeeze(reg), [-3:.15:5], 0, 'k')
  axis([0 102 0 102]);