Documentation of plot_30s30n_EOF_regressions

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

Help text

load gr_pcs.mat lam lds pcs per lims cthp gr 
grpcs = pcs;
for i = 1:900;
  ind = 12*(i-1)+[1:12];
  tem(i,:) = mean(grpcs(ind,:));
end
pcs = tem;

Cross-Reference Information

Listing of script plot_30s30n_EOF_regressions

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat

cd ~/matlab/CSIRO/Heat/Old_routines

nfrm = -12:2:12; lags = nfrm;
tim = 101:1000;
lev = 1:3;
lims = [100 305 -35 35];

tit = 'LP9';

varn1 = 'taux';
varn2 = 'tauy';

[var1, lat, lon] = getheat(lims, 1:3, tim);
%var1 = getnc(varn1, lims, 1, tim);
%[lat, lon, depth] = getll(varn1, lims);

[var2, lat2, lon2] = getheat(lims, 4:7, tim);
%var2 = getnc(varn2, lims, 1, tim);
%[lat2, lon2, depth2] = getll(varn2, lims);

%lev = 4:7;
%[var1] = getheat(lims, lev, tim, 'u');
%[var2, lat2, lon2] = getheat(lims, lev, tim, 'v');

reg1 = regress_eof(detrend(var1), detrend(pcs), nfrm);
reg2 = regress_eof(detrend(var2), detrend(pcs), nfrm);

%  Plot the data

default_global;  XAX = lon; YAX = lat; FRAME = [105 300 -30 30];

num = 1
cint = 0.1; 
figure(num); figure_orient; clf;
for i = 8:10; ind = i-7;
  pos = [.2 .92-ind*.1415 .2920 .1072];
  subplot('position', pos);
%  subplot(6,2,2*i-2);
    gcont(-1e-8*reg1(i,:,:), cint);
    dc2(reg1(i,:,:));
    ylabel(['Lag = ' num2str(lags(i))]);
  if ind == 1;
%    title([tit ' 0-80m HC']);
  end
%    xlabel(['Contour Interval:  ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']);
%    xlabel(['Contour Interval:  ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
  end
end
%subplot(6,2,1)
%  title([tit ' TAUX']);
%  title([tit ' 0-80m HC']);

default_global;  XAX = lon2; YAX = lat2; FRAME = [105 300 -30 30];

cint = 0.2;
figure(num); figure_orient;
for i = 8:10; ind = i-7;
  pos = [.558 .92-ind*.1415 .2920 .1072];
  subplot('position', pos);
%  subplot(6,2,2*i);
    gcont(-1e-8*reg2(i,:,:), cint);
    dc2(reg2(i,:,:));
%    ylabel(['Lag = ' num2str(lags(i))]);
  if ind == 1;
%    title([tit ' 80-270m HC']);
  end
%    xlabel(['Contour Interval:  ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
%    xlabel(['Contour Interval:  ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']);
  end
end
%subplot(6,2,2)
%  title([tit ' Zonal Wind Stress']);
%  title([tit ' 80-270m HC']);

cd ~/Thesis/Chap4

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 sfc_T_taux_PC1.ps




%  Plot HC80, HC80-270, TAUX on one figure

figure(1); figure_landscape;

default_global;  XAX = lon; YAX = lat; FRAME = [110 299 -30 30];
default_global;  XAX = lon2; YAX = lat2; FRAME = [110 299 -30 30];

cint = 1; clev = [-12:cint:-cint cint:cint:12];
for i = 1:5;
  ind = i+1;
  subplot(5,3,3*i);
    gcont(1e2*reg2(ind,:,:), clev);
    dc2(reg1(i,:,:));
    ylabel(['Lag = ' num2str(lags(ind))]);
  if i == 5;
    xlabel(['Contour Interval:  ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
%    xlabel(['Contour Interval:  ' num2str(10*cint) 'x10^7 J m^-^2 std^-^1']);
  end
end
subplot(5,3,3)
  title([tit ' Zonal Wind Stress']);
  title([tit ' 80-270m HC']);
  title([tit ' 0-80m HC']);


cd ~/Thesis/Chap3




%  Plot 60s to 60n data

lev = 1:3;
[var1] = getheat(lims, lev, tim, 'wl');
[var2, lat2, lon2] = getheat(lims, lev, tim, 'wl');

reg1 = regress_eof(var1, pcs, nfrm);
reg2 = regress_eof(var2, pcs, nfrm);

get_global; FRAME = [105 300 -60 60];
lm = getnc('temp', lims, 1, 1);
[latt, lont] = getll('temp', lims);


lind = 0;
ind = find(lags == lind);

figure(1); figure_landscape;
subplot(2,2,1);
  XAX = lon2; YAX = lat2;
  gquiv(reg1(ind,:,:), reg2(ind,:,:), -0.75, 1);
  XAX = lont; YAX = latt;
  dc2(lm);
  title('LP9:  0-80m < U'', V'' >');

subplot(2,2,2);
  XAX = lon2; YAX = lat2;
  gquiv(reg1(ind,:,:), reg2(ind,:,:), -0.25, 1);
  XAX = lont; YAX = latt;
  dc2(lm);
  title('LP9:  80-270m < U'', V'' >');

cint = 0.05; clev = [-5:cint:-cint cint:cint:5];
subplot(2,2,3);
  xlabel('Contour Interval:  5x10^-^6 cm s^-1 std^1')
  XAX = lon2; YAX = lat2;
  gcont(reg1(ind,:,:), clev);
  XAX = lont; YAX = latt;
  dc2(lm);
  title('LP9:  0-80m W''');

cint = 0.025; clev = [-5:cint:-cint cint:cint:5];
subplot(2,2,4);
  xlabel('Contour Interval:  2.5x10^-^6 cm s^-1 std^1')
  XAX = lon2; YAX = lat2;
  gcont(reg2(ind,:,:), clev);
  XAX = lont; YAX = latt;
  dc2(lm);
  title('LP9:  80-270m W''');

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots
%  print -dps2 U_V_W_zerolag_regression.ps