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 HP10_detrend_L1-7_CEOF_yr101-1000.mat
cd ~/matlab/CSIRO/Heat/Old_routines
nfrm = 6; lags = 180/nfrm * [0:nfrm-1];
tim = 101:1000;
lev = 1:3;
lims = [106 312 -31 31];
tit = 'HP10';
varn1 = 'temp';
varn2 = 'tauy';
[var1, 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_ceof(var1, pcs, nfrm);
reg2 = regress_ceof(var2, 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*reg1(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