Documentation of regress_CPC_windstress

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

Help text

  Load CPCs

Cross-Reference Information

Listing of script regress_CPC_windstress

clear
cd /home/disk/hayes2/dvimont/csiro/data
tim = 101:550;
ctlim = [109 302 -62 62];
nc = netcdf('zonal_stress_A_1000_years.nc', 'nowrite');
  lat1 = nc{'latitude'}(:);
  lon1 = nc{'longitude'}(:);
  [xk, yk] = keep_var(ctlim, lon1, lat1);
  taux = nc{'smfzon'}(tim,yk,xk);
  mv = nc{'smfzon'}.missing_value(:);
nc = close(nc);
taux(find(taux == mv)) = NaN * ones(size(find(taux == mv)));
nc = netcdf('merid_stress_A_1000_years.nc', 'nowrite');
  tauy = nc{'smfmer'}(tim,yk,xk);
  mv = nc{'smfmer'}.missing_value(:);
nc = close(nc);
tauy(find(tauy == mv)) = NaN * ones(size(find(tauy == mv)));
lat1 = lat1(yk); lon1 = lon1(xk);
[taux, clim] = remove_mean(taux);
[tauy, clim] = remove_mean(tauy);
[ntim, nlat, nlon] = size(taux);
taux = reshape(taux, ntim, nlat*nlon);
tauy = reshape(tauy, ntim, nlat*nlon);

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
%load LP10_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 10 Years )';
%load LP8_L1-7_CEOF.mat; tit = 'Lowpass Filtered Data ( > 8 Years )';
load RAW_L1-7_CEOF.mat; tit = 'Unfiltered Data';

%  Load CPCs

%for biff2 = 1:3;
biff2 = 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 regressions

lag = 0; lg = lag*pi/180; lg2 = 1;
num = 1; lind = 1;
nfrm = 6;

%  Store all regressions under one variable

j = sqrt(-1);
timeseries = sqrt(2)*pcs(:,num)./std(pcs(:,num));
clear temtim ureg vreg
for i = 1:nfrm
  wgt = conj(exp(j * ((i-1) * pi/(lg2*nfrm) + lg) ));
  temtim(:,i) = squeeze(real(wgt .* timeseries));
  ureg(i, :) = temtim(:,i)' * taux ./ ntim;
  vreg(i, :) = temtim(:,i)' * tauy ./ ntim;
end
ureg = reshape(ureg, nfrm, nlat, nlon);
vreg = reshape(vreg, nfrm, nlat, nlon);

%  Plot the data

default_global; XAX = lon1; YAX = lat1; FRAME = [110 300 -60 60];

figure(1); figure_orient;
for i = 1:nfrm
  temu = squeeze(ureg(i,:,:))/10;
  temv = squeeze(vreg(i,:,:))/10;
  
  subplot(3,2,i);
    [h, hx] = gquiv(temu, temv, -0.01, 2, 'dynes cm^-^2');
    dc2(temu)
    title(['Phase = ' num2str((i-1)*180/nfrm + lag)]);
%    xlabel(['Max = ' num2str(round(c(2)*100)/100) ' cm s^-^1'])
end
subplot(3,2,3);
  ylabel([tit ':  Wind Stress regressed on HC CPC1;  '...
          'Depth = 0:270m']);

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

end