Documentation of regress_vertical_temp

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

Help text

[ntim, nlev, nlat, nlon] = size(temp1);

Cross-Reference Information

Listing of script regress_vertical_temp

clear
cd /home/disk/hayes2/dvimont/csiro/data
lonlim = [110 285];
latlim = [[-37.5:15:22.5]' [-22.5:15:37.5]'];
filin = ['temp_A_L1-10.nc'];
nc = netcdf(filin, 'nowrite');
  depth = nc{'depth'}(:);
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
for i = 1:5;
  ctlim = [lonlim latlim(i, :)];
  [xk, yk] = keep_var(ctlim, lon, lat);
  temp = nc{'temp'}(:,:,yk,xk);
  mv = nc{'temp'}.missing_value(:);
  temp(find(temp == mv)) = NaN * ones(size(find(temp == mv)));
  eval(['temp' num2str(i) ' = shiftdim(squeeze(mean2(shiftdim(temp, 2))), 1);']);
  eval(['lat' num2str(i) ' = lat(yk); lon' num2str(i) ' = lon(xk);']);
end;
nc = close(nc);
[ntim, nlev, nlon] = size(temp1);

cd /home/disk/tao/dvimont/matlab/CSIRO/Data
load butter_4.5_ctstar.mat
ctlow = detrend(ctlow);
ctlow = (ctlow - mean(ctlow)) / std(ctlow);

for i = 1:5;
  eval(['temp' num2str(i) ' = reshape(temp' num2str(i) ', ntim, nlev*nlon);']);
  eval(['temp' num2str(i) ' = detrend(temp' num2str(i) ');']);
end

figure(1);  figure_orient; clf;

lattit = [-30:15:30];

for i = 1:6;
  figure(i); figure_orient; clf;
  for j = 1:5;
    eval(['temp = temp' num2str(j) ';']);
    eval(['lon = lon' num2str(j) ';']);
    tit1 = num2str(lattit(i));
    for k = 1:2;

      spnum = (5-j)*2 + k;
      subplot(5,2,spnum);

      l = 2*(i-4) + (k-1);
      if l <= 0;
        tem = ctlow((1-l):ntim)' * temp(1:(ntim+l),:) / (ntim+l);
        ptit = ['lag_' num2str(abs(l))];
        tit = ['lag = ' num2str(abs(l))];
      else
        tem = ctlow(1:(ntim-l))' * temp((l+1):ntim,:) / (ntim-l);
        ptit = ['lead_' num2str(l)];
        tit = ['lead = ' num2str(l)];
      end
 
      tem = reshape(tem, nlev, nlon);

      subplot(5,2,spnum);

      pncont(lon, -0.01*depth, tem, [-.2:.01:.2], 'k');
      grid on
      title(['<' tit1 ', LPCT> :  ' tit]);
      set(gca, 'XTick', 150:30:300, 'YTick', -500:100:0, 'YTickLabel', 500:-100:0);
      axis([150 280 -545 -5])

    end
  end
end

for i = 6:-1:1; figure(i); end;