Documentation of regress_sv

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script regress_sv

clear
lev = 1:9;
tim = 101:550;
lims = [100 300 8 70];
nfrm = -20:2:20;

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_detrend_L1-7_yr101-550.mat
pcs = -1*pcs;

%  Load ocean velocities

[v, latu, lonu, depthu, middepth] = getheat(lims, lev, tim, 'v');
u = getheat(lims, lev, tim, 'u');

ureg = regress_eof(u/100, pcs, nfrm);
vreg = regress_eof(v/100, pcs, nfrm);

[taux, tauy] = getnc('taux', 'tauy', lims, lev, tim);
[latt, lont] = getll('taux', lims);

xreg = regress_eof(taux/10, pcs, nfrm);
yreg = regress_eof(tauy/10, pcs, nfrm);

strf1 = NaN * ones(length(nfrm), length(latu), length(lonu)+1);
strf2 = NaN * ones(length(nfrm), length(latt), length(lont)+1);
H = middepth(max(lev)+1) - middepth(min(lev));
for i = 1:length(nfrm);
  strf1(i,:,:) = vel_to_strfnc(squeeze(vreg(i,:,:)), latu, lonu);
  [strf2(i,:,:), lat_out2, lon_out2] = tau_to_strf(...
      squeeze(xreg(i,:,:)), squeeze(yreg(i,:,:)), latt, lont, ...
      squeeze(strf1(i,:,:)), H);
end

figure(1)
default_global; XAX = lon_out2; YAX = lat_out2; FRAME = lims;

cint = 500; clev = [-5000:cint:5000];
for i = 1:5;
  subplot(5,2,2*i-1);
    gcont(strf1(i+6,:,:), clev);
    ylabel(['Lag = ' num2str(nfrm(i+6))]);
    dc2(strf1(i+6,:,:));
    axis([110 280 8 60])
  subplot(5,2,2*i);
    gcont(strf2(i+6,:,:), clev);
    ylabel(['Lag = ' num2str(nfrm(i+6))]);
    dc2(strf2(i+6,:,:));
    axis([110 280 8 60])
end

figure(2)
default_global; XAX = lon_out2; YAX = lat_out2; FRAME = lims;

cint = 500; clev = [-5000:cint:5000];
for i = 1:5;
  subplot(5,2,2*i-1);
    gcont(strf1(i+11,:,:), clev);
    ylabel(['Lag = ' num2str(nfrm(i+6))]);
    dc2(strf1(i+6,:,:));
    axis([110 280 8 60])
  subplot(5,2,2*i);
    gcont(strf2(i+11,:,:), clev);
    ylabel(['Lag = ' num2str(nfrm(i+6))]);
    dc2(strf2(i+6,:,:));
    axis([110 280 8 60])
end