Documentation of ceof_regress_velocity

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

Help text

  Get PCs

Cross-Reference Information

Listing of script ceof_regress_velocity

clear
lev = 1;
lims = [106 310 -63 63];
tim = 101:550;
nfrm = 6;

for biff2 = 3:3;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
  load LP10_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
  ptit = 'LP10';
elseif biff2 == 2;
  load HP8_L1-7_CEOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
  ptit = 'HP8';
elseif biff2 == 3;
  load RAW_L1-7_CEOF.mat; tit = 'Unfiltered Data';
  ptit = 'RAW';
end

%  Get data

[u, v] = getnc('taux', 'tauy', lims, lev, tim);
[u, clim] = remove_mean(u);
[v, clim] = remove_mean(v);

[lat, lon] = getll('taux', lims);

%  Get regressions

ureg = regress_ceof(u, pcs, nfrm);
vreg = regress_ceof(v, pcs, nfrm);

%  Plot the data

default_global;  FRAME = [110 299 -60 60];
lag = 0;

figure(biff2); figure_orient(1);
for i = 1:nfrm
  temu = squeeze(ureg(i,:,:));
  temv = squeeze(vreg(i,:,:));
  
    subplot(3,2,i);
    [h, hx] = gquiv(temu, temv, -0.075, 2, 'dynes cm^-^2');
    dc
    title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
end
for i = 1:4;
  subplot(3,2,i);
  xlabel([]);
end
subplot(3,2,3);
  ylabel([tit ':  Wind Stress regressed on HC (L1-7) CPC1']);%;  '...
%          'Depth = ' num2str(depth(lev)) 'm']);
end

if 0;
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Plot_Ceof
eval(['print -dps2 ' ptit '_WSTRESS_L1-7.ps'])

%_oce_vel_L' num2str(lev) '.ps']);

end;