Documentation of get_windstress

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

Help text

  Currently, this is the model wind stress:

    taustr = cos(3.1415927*(j-51-0.5)/100.)

Cross-Reference Information

Listing of script get_windstress

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1-7_EOF_yr101-1000.mat
lims = [110 300 -90 90];
lev = 1;
tim = 101:1000;
taux = getnc('taux', lims, 1, tim);
[lat, lon] = getll('taux', lims);
mtaux = squeeze(mean2(taux));
zon_taux = squeeze(mean2(shiftdim(mtaux, 1)));
figure(1); figure_landscape;
subplot(2,1,1);
plot(lat, zon_taux);
rtaux = regress_eof(taux, pcs, 0);
rtaux = squeeze(rtaux);
zon_rtaux = squeeze(mean2(shiftdim(rtaux, 1)));
subplot(2,1,2);
plot(lat, -zon_rtaux);
grid on

taustr = zeros(102);
for i = 11:92;
  for j = 11:92;
%    taustr(j,i) = cos(pi*(j-51-0.5)/100);
%    taustr(j,i) = exp(-1*((j-51-0.5)/8.413)^2);
    taustr(j,i) = exp(-1*((j-51-0.5)/8.413)^2)*...
                  exp(-1*((i-51-0.5)/20)^2);
  end
end

latgrd = [-50.5:1:50.5]*180/101

subplot(2,1,1)
plot(latgrd, mean(taustr, 2));
axis([-90 90 -0.5 1.5]);
set(gca, 'XTick', -90:10:90)
grid on

subplot(2,1,2);
contourf(taustr);