Documentation of wstp_upwelling

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

Help text

  First look at mean:

Cross-Reference Information

Listing of script wstp_upwelling

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_detrend_L1-7_EOF_yr101-1000.mat
tim = 101:1000;
lims2 = [110 165 5 10]; tit = 'WSTP'; %  WSTP index
[taux, tauy] = getnc('taux', 'tauy', lims2, 1, tim);
[lat, lon] = getll('taux', lims2);
w = getnc('wl', lims2, 1:9, tim); w = 1e-6*w;
[latw, lonw, depthw] = getll('wl', lims2);

wreg = squeeze(mean(w)); 
wreg = [zeros(1, size(wreg, 2), size(wreg,3)); wreg];
depthw = [0; depthw];
txreg = squeeze(mean(taux));
tyreg = squeeze(mean(tauy));

%  Get Ekman upwelling

[nlat, nlon] = size(txreg);
f = (2*7.292e-5)*sin(lat*pi/180)*ones(1, nlon);
txreg = 0.1*txreg./f;
tyreg = 0.1*tyreg./f;
[wek, lat2, lon2] = sph_curl1(txreg, tyreg, lat, lon, 0);
wek = wek./1000;

%  Get thickness of w-layers

[tem1, tem2, depth] = getll('temp', lims2);
dz = diff([0; depth]);

%  compare w and wek

[nlev,  nlon] = size(wreg);

for i = 1:8;
lev = 1:(i+1);
  wtem = dz(lev)'*reshape(wreg(lev,:), length(lev), nlon);
  wave(i,:) = reshape(wtem, 1, nlon)./sum(dz(lev));
end

get_global; XAX = lon2; YAX = lat2; FRAME = lims2;

sptalk(3,2,1);
  gcont(0.5*1e6*wek, 0.2);
for i = 2:6;
sptalk(3,2,i);
  gcont(1e6*wave(i,:,:), 0.2);
end

squeeze(mean2(shiftdim(wave, 1)))
0.5*mean2(mean2(wek))




%  That doesn't show much, so try regressions:

clear

loadpcs;
tim = 101:1000;
pcs = hpcs;

lims2 = [135 180 3 16]; tit = 'WSTP'; %  WSTP index
lims2 = [131 166 5 10]; tit = 'WSTP'; %  WSTP index
lims2 = [110 150 5 13]; tit = 'WSTP'; %  WSTP index

[taux, tauy] = getnc('taux', 'tauy', lims2, 1, tim);
[lat, lon] = getll('taux', lims2);

w = getnc('wl', lims2, 1:9, tim); w = 1e-6*w;
[latw, lonw, depthw] = getll('wl', lims2);

%  Get regressions:

lags = -10:10;
wreg = regress_eof(w, pcs, lags);
txreg = regress_eof(taux, pcs, lags);
tyreg = regress_eof(tauy, pcs, lags);

[nlag, nlat, nlon] = size(txreg);
f = (2*7.292e-5)*sin(lat*pi/180)*ones(1, nlon);

clear wek
for i = 1:nlag
  temx = squeeze(0.1*txreg(i,:,:))./f;
  temy = squeeze(0.1*tyreg(i,:,:))./f;
  [wek(i,:,:), lat2, lon2] = sph_curl1(temx, temy, lat, lon, 0);
end
wek = wek./1000;

%  Average wreg

[tem1, tem2, depth] = getll('temp', lims2);
dz = diff([0; depth]);

[nlag, nlev, nlat, nlon] = size(wreg);
wreg = reshape(shiftdim(wreg, 1), nlev, nlat*nlon*nlag);
wreg = [zeros(1, nlat*nlon*nlag); wreg];

clear wtem
for i = 1:(nlev-1);
  lev = 1:(i+1);
  wtem(i,:) = dz(lev)'*wreg(lev,:)./sum(dz(lev));
end

wtem = shiftdim(reshape(wtem, size(wtem, 1),  nlat, nlon, nlag), 3);

figure(2); fl(1);
sptalk(1,1,1);
  tem1 = 0.5*squeeze(mean2(mean2(shiftdim(wek, 1))));
  tem2 = squeeze(mean2(mean2(shiftdim( ...
          wtem(:,:,:,:), 2))));
  plot(lags, tem1, 'x-k', lags, tem2, '^-k')



1.7052e-07