Documentation of plot_30s30n_EOF_regressions

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

Help text

load LP9_detrend_L1-7_EOF_yr101-1000.mat

Cross-Reference Information

Listing of script plot_30s30n_EOF_regressions

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

cd ~/matlab/CSIRO/Heat/Old_routines

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

tit = 'HP10';

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

varn1 = 'temp';
varn2 = 'taux';

%[var1, lat, lon] = getheat(lims, 1:3, tim);
var1 = getnc(varn1, lims, 4, 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(var1, pcs, nfrm);
reg2 = regress_eof(var2, pcs, nfrm);

%  Plot the data

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

clf
cint = 0.05; 
figure(1); figure_orient;
for i = 1:5;
  pos = [.2 .92-i*.1415 .2920 .1072];
  subplot('position', pos);
    gcont(1e0*reg1(i,:,:), cint);
    dc2(reg1(i,:,:));
    ylabel(['Lag = ' num2str(lags(i))]);
  if i == 1;
    title(['< 0-80m HC, HPPC1 >']);
  elseif i == 5;
%    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(['< 0-80m HC, HPPC1 >']);

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

cint = 0.05;
figure(1); figure_orient;
for i = 1:5;
  pos = [.558 .92-i*.1415 .2920 .1072];
  subplot('position', pos);
    gcont(1e0*reg1(i,:,:), cint);
    dc2(reg1(i,:,:));
  if i == 1;
    title(['< 80-270m HC, HPPC1 >']);
  elseif i == 5;
%    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/Chap3

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots
%print -dps2 HC0-80_80-270_lag-2to2_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