Documentation of wind_stress

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

Help text

  load LP10_L1-7_CEOF.mat; tit = 'LP10';
  load HP8_L1-7_CEOF.mat; tit = 'HP8';
  load RAW_L1-7_CEOF.mat; tit = 'RAW';

Cross-Reference Information

Listing of script wind_stress

clear
varn = 'temp';
lims = [110 300 -60 60];
lev = 1;
tim = 101:550;
nfrm = -6:5;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
  load BP10-60yr_L1-7_EOF.mat; tit = 'LP10'; %pcs = -1*pcs;

%[taux, tauy] = getnc('temp', 'tauy', lims, lev, tim);
taux = getnc(varn, lims, lev, tim);
%taux = getheat(1:7, tim, lims);
[lat, lon] = getll(varn, lims);

[b, a] = butter(6, 2/10);
[b2, a2] = butter(6, 2/60);
taux = filtfilt(b, a, taux) - filtfilt(b2, a2, taux);

[xreg, xcor] = regress_eof(taux, pcs, nfrm);
%[yreg, ycor] = regress_eof(tauy, pcs, nfrm);
%[xreg, xcor] = regress_ceof(taux, pcs, nfrm);
%[yreg, ycor] = regress_ceof(tauy, pcs, nfrm);

%  Plot the data

default_global;

figure(1);
for i = 1:6;
  tem = squeeze(xcor(i,:,:));
  subplot(3,2,i);
    gcont(tem, [-1:.1:1]);
    dc2(tem);
    title(['Lag:  ' num2str(nfrm(i))]);
end
figure(2);
for i = 7:12;
  tem = squeeze(xcor(i,:,:));
  subplot(3,2,i-6);
    gcont(tem, [-1:.1:1]);
    dc2(tem);
    title(['Lag:  ' num2str(nfrm(i))]);
end

%  Look at convergence

[nfrm, nlat, nlon] = size(xreg);
for i = 1:nfrm;
  [divtau(i,:,:), lat2, lon2] = sph_div1(squeeze(xreg(i,:,:)), ...
					 squeeze(yreg(i,:,:)), lat, lon, 0);
end

XAX = lon2; YAX = lat2;

figure(2); figure_orient;
for i = 1:6;
  tem = 1e8*squeeze(divtau(i,:,:));
  subplot(3,2,i);
    gcont(tem, [-10:.5:10]);
    dc2(tem);
end