Documentation of zonal_depth_eq

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

Help text

 print -dps2 eq_tcline_neglag.ps

Cross-Reference Information

Listing of script zonal_depth_eq

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_detrend_L1-7_yr101-550.mat
lims = [105 310 -5 5];                                     
tim = 101:550;
nfrm = -10:10;
lev = 1:10;
temp = getnc('temp', lims, lev, tim);
[lat, lon, depth] = getll('temp', lims);
temp = squeeze(mean(temp, 3));
treg = regress_eof(temp, pcs, nfrm);    
default_global;
figure(1);
for i = 1:5;
  subplot(5, 2, 2*i-1);
    tem = squeeze(treg(i,:,:));
    pncont(lon, -1*depth, tem, [-1:.02:1], 0, 'k');
    dcv(lon, -1*depth, tem);
    ylabel(['Lag = ' num2str(nfrm(i))]);
    set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0);
    manual_ticks('x', 120:30:300, label_eastwest(120:30:300));
    axis([120 280 -510 0]);
  subplot(5, 2, 2*i);
    tem = squeeze(treg(i+5,:,:));
    pncont(lon, -1*depth, tem, [-1:.02:1], 0, 'k');
    dcv(lon, -1*depth, tem);
    ylabel(['Lag = ' num2str(nfrm(i+5))]);
    set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0);
    manual_ticks('x', 120:30:300, label_eastwest(120:30:300));
    axis([120 280 -550 0]);
end
subplot(5,2,9);
  xlabel('Contour Interval:  0.02 K std^-^1');
subplot(5,2,10);
  xlabel('Contour Interval:  0.02 K std^-^1');
subplot(5,2,1);
  title('Equatorial Temperature');
subplot(5,2,2);
  title('Equatorial Temperature');
figure(2);% figure_orient;
for i = 1:5;
  subplot(5, 2, 2*i-1);
    tem = squeeze(treg(i+10,:,:));
    pncont(lon, -1*depth, tem, [-1:.02:1], 0, 'k');
    dcv(lon, -1*depth, tem);
    ylabel(['Lag = ' num2str(nfrm(i+10))]);
    set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0);
    manual_ticks('x', 120:30:300, label_eastwest(120:30:300));
    axis([120 280 -500 0]);
  subplot(5, 2, 2*i);
    tem = squeeze(treg(i+15,:,:));
    pncont(lon, -1*depth, tem, [-1:.02:1], 0, 'k');
    dcv(lon, -1*depth, tem);
    ylabel(['Lag = ' num2str(nfrm(i+15))]);
    set(gca, 'YTick', -500:125:0, 'YTickLabel', 500:-125:0);
    manual_ticks('x', 120:30:300, label_eastwest(120:30:300));
    axis([120 280 -500 0]);
end
subplot(5,2,9);
  xlabel('Contour Interval:  0.02 K std^-^1');
subplot(5,2,10);
  xlabel('Contour Interval:  0.02 K std^-^1');
subplot(5,2,1);
  title('Equatorial Temperature');
subplot(5,2,2);
  title('Equatorial Temperature');
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_EOF
figure(1);
figure(2);
% print -dps2 eq_tcline_poslag.ps