Documentation of plot_lag_regress2

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script plot_lag_regress2

clear

biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
%  load LP10_L1-7_CEOF.mat; tit = 'LP10';
  load BP10-60yr_L1-7_EOF.mat; tit = 'BP10-60';
  load LP10_detrend_L1-7_yr101-550.mat; tit = 'LP10';
elseif biff2 == 2;
  load HP8_L1-7_CEOF.mat; tit = 'HP8';
elseif biff2 == 3;
  load RAW_L1-7_CEOF.mat; tit = 'RAW';
end

cd ~/matlab/CSIRO/Heat/Old_routines

lags = [-10:10]; nfrm = length(lags);
tim = 101:550;
lev = 1:7;
lims = [106 312 -44 63];

%  Start with SFC temp and TAUX

[temp, taux] = getnc('temp', 'taux', lims, 1, tim);
[latt, lont, deptht] = getll('temp', lims);
[latx, lonx] = getll('taux', lims);

szt = size(temp);
szx = size(taux);
temp = reshape(temp, szt(1), prod(szt(2:3)));
taux = reshape(taux, szx(1), prod(szx(2:3)));

[b, a] = butter(6, 2/10);
temp = filtfilt(b, a, detrend(temp));
taux = filtfilt(b, a, detrend(taux));
temp = reshape(temp, szt);
taux = reshape(taux, szx);

pcs = -1*pcs;
[treg, tcoef] = regress_eof(temp, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);

%  Plot Temperature

lag = 1;

dofx = get_dof(pcs(:,1));
score = tscore(dofx, 2.5);
tval = tcoef * sqrt(dofx - 1) ./ sqrt(1 - tcoef.^2);

default_global;  XAX = lont; YAX = latt; FRAME = [110 299 -40 60];

cint = 0.025; clev = [-5:cint:-cint 0 cint:cint:5];
figure(1); figure_orient;
for i = 1:4;
  subplot(4,2,2*i-1);
    ind = i + lag;
    gcont(treg(ind,:,:), clev);
    dc2(treg(ind,:,:));
    tem = zeros(size(tval(ind,:,:)));
    tem(find(tval(ind,:,:) > score)) = NaN;
    shade_ccoef(tem);
    ylabel(['Lag = ' num2str(lags(ind))]);
    if i == 4;
      xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);
%      xlabel(['Contour Interval:  ' num2str(cint) ' cm s^-^1 std^-^1']);
    end
end
subplot(4,2,1); 
  title([tit ':  TEMP']);
%  title([tit ':  0:80m V']);
%  title([tit ':  80-270m V']);

%  Plot TauX

dofx = get_dof(pcs(:,1));
score = tscore(dofx, 2.5);
tval = xcoef * sqrt(dofx - 1) ./ sqrt(1 - xcoef.^2);

default_global;  XAX = lonx; YAX = latx; FRAME = [110 299 -40 60];

cint = 0.01; clev = [-15:cint:-cint 0 cint:cint:15];
figure(1); figure_orient;
for i = 1:4;
  subplot(4,2,2*i);
    ind = i + lag;
    gcont(xreg(ind,:,:), clev);
    dc2(xreg(ind,:,:));
    tem = zeros(size(tval(ind,:,:)));
    tem(find(tval(ind,:,:) > score)) = NaN;
    shade_ccoef(tem);
    if i == 4;
      xlabel(['Contour Interval:  ' num2str(cint) ' dynes cm^-^2 std^-^1']);
%      xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^1']);
%      xlabel(['Contour Interval:  ' num2str(cint) ' J m^-^2 std^-^1']);
    end
end
subplot(4,2,2); 
  title([tit ':  TAUX']);
%  title([tit ':  HFLX']);
%  title([tit ':  80-270m HC']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Mid_Lats/Plot_lag6_lag3
%cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Mid_Lats/Plot_lag8_lag5

%  Look at Heat Flux and V - component of ocean current (top 80m)

lev = 1:3;
[temp, latt, lont, depth, middepth] = getheat(lev, tim, lims, 'v');
taux = getnc('heat', lims, 1, tim);
[latx, lonx] = getll('heat', lims);

szt = size(temp);
szx = size(taux);
temp = reshape(temp, szt(1), prod(szt(2:3)));
taux = reshape(taux, szx(1), prod(szx(2:3)));

[b, a] = butter(6, 2/10);
temp = filtfilt(b, a, detrend(temp));
taux = filtfilt(b, a, detrend(taux));
temp = reshape(temp, szt);
taux = reshape(taux, szx);

[treg, tcoef] = regress_eof(temp, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);


%  Look at Lower Level Heat Content Heat and V - component of ocean current (80-270m)

lev = 4:7;
[temp, latt, lont, depth, middepth] = getheat(lev, tim, lims, 'v');
[taux, latx, lonx] = getheat(lev, tim, lims);

szt = size(temp);
szx = size(taux);
temp = reshape(temp, szt(1), prod(szt(2:3)));
taux = reshape(taux, szx(1), prod(szx(2:3)));

[b, a] = butter(6, 2/10);
temp = filtfilt(b, a, detrend(temp));
taux = filtfilt(b, a, detrend(taux));
temp = reshape(temp, szt);
taux = reshape(taux, szx);

[treg, tcoef] = regress_eof(temp, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);