Documentation of plot_30s30n_regressions

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

Help text

  Plot the data

Cross-Reference Information

Listing of script plot_30s30n_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 = -6:6; lags = nfrm;
tim = 101:1000;
lev = 1:3;
lims = [106 312 -35 35];
tit = '0-80m LP9';
varn1 = 'taux';
varn2 = 'tauy';
[var3, lat, lon] = getheat(lims, 4:7, tim);
var1 = getnc(varn1, lims, 1, tim);
[lat, lon, depth] = getll(varn1, lims);
[var2, lat2, lon2] = getheat(lims, 1:3, tim);
var2 = getnc(varn2, lims, 1, tim);
[lat2, lon2, depth2] = getll(varn2, lims);
reg1 = regress_eof(var1, pcs, nfrm);
reg2 = regress_eof(var2, pcs, nfrm);
reg3 = regress_eof(var3, pcs, nfrm);

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

cint = 0.25; clev = [-10:cint:-cint cint:cint:10];
figure(1); figure_orient;
for i = 1:6;
  subplot(6,2,2*i-1);
    gcont(1e-8*reg3(i,:,:), clev);
    dc2(reg1(i,:,:));
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);
  end
end
subplot(6,2,1)
  title([tit ' 12.5m Temp']);
  title([tit ' 80-270m HC']);

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

cint = 1; clev = [-10:cint:-cint cint:cint:10];
figure(1); figure_orient;
for i = 1:6;
  subplot(6,2,2*i);
    gcont(1e2*reg2(i,:,:), clev);
    dc2(reg2(i,:,:));
    ylabel(['Phase = ' num2str(lags(i))]);
  if i == 6;
    xlabel(['Contour Interval:  ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']);
  end
end
subplot(6,2,2)
  title([tit ' Zonal Wind Stress']);
  title([tit ' 0-80m HC']);

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







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

figure(2); figure_landscape(1);

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

cint = 0.5; clev = [-12:cint:-cint 0 cint:cint:12];
ind1 = 0;

  for i = 1:6;
    ind = i+ind1;
    subplot(6,3,3*i);
      gcont(-1e2*reg1(ind,:,:), clev);
      XAX = lon2; YAX = lat2;
      dc2(reg2(i,:,:));
      XAX = lon; YAX = lat;
      ylabel(['Lag = ' num2str(lags(ind))]);
    if i == 6;
      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(6,3,3)
      title([tit ' Zonal Wind Stress']);
      title([tit ' 80-270m HC']);
      title([tit ' 0-80m HC']);

cd ~/Thesis/Chap3


for i = 1:18;
  subplot(6,3,i);
  set(gca, 'Fontsize', 8);
end