Documentation of hofmoeller_depth_temp

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

Help text

  ctlim = [-.1 360 -6 6];

Cross-Reference Information

Listing of script hofmoeller_depth_temp

clear
cd /home/disk/hayes2/dvimont/csiro/data
filin = ['temp_A_L1-10.nc'];
nc = netcdf(filin, 'nowrite');
  depth = nc{'depth'}(:);
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  ctlim = [120 290 -45 45];
  [xk, yk] = keep_var(ctlim, lon, lat);
  tim = 466:515;

  temp = nc{'temp'}(tim,:,yk,xk);
  mv = nc{'temp'}.missing_value(:);

nc = close(nc);
[ntim, nlev, nlat, nlon] = size(temp);
temp(find(temp == mv)) = NaN * ones(size(find(temp == mv)));
lat = lat(yk); lon = lon(xk);

temp = reshape(temp, ntim, nlev*nlat*nlon);
temp = detrend(temp);
[bl,al]=butter(6,(2/4.5));
[bh,ah]=butter(6,(2/12));
%temp = filtfilt(bl, al, temp) - filtfilt(bh, ah, temp);
temp = filtfilt(bl, al, temp);
temp = reshape(temp, ntim, nlev, nlat, nlon);

%  Look at longitude/time sections

lims = flipud([-8 8]);
tim = 1:50;
for j = 1:4
figure(j); clf; figure_orient(1);
colormap('pink');
for i = 1:5;
  [xk, yk] = keep_var([ctlim(1:2) lims(i,:)], lon, lat);
  tem = squeeze(mean2(shiftdim(temp(tim,((2*j)-1),yk,xk),1)));

  subplot(1,5,i);
  contourf(lon(xk), (tim+465), tem', 10);
%  hold on;
%  contour(lon(xk), (tim+465), tem', [0 0], '-k');
%  hold off;
%  colorbar('horiz');
  title(['<'num2str(lims(i,:)) '>']);%, ' num2str(depth((2*j)-1)/100) 'm']);
%  xlabel('longitude'); ylabel('time');

end
  subplot(1,5,1); ylabel('Year');
  subplot(1,5,3); 
  xlabel(['Time-Longitude Hofmoeller Diagram, Depth = '...
           num2str(depth((2*j)-1)/100) 'm']);
end
for i=5:-1:1; figure(i); colormap('pink'); end;

cd /home/disk/tao/dvimont/matlab/CSIRO/Plots3


%  Look at latitude/time sections

lims = [120 150; 150 180; 180 210; 210 240; 240 270];
for j = 1:5
figure(j); set(gcf, 'Position', [2 222 950 750]);
for i = 1:5;

  [xk, yk] = keep_var([lims(j,:) ctlim(3:4)], lon, lat);
  tem = squeeze(mean2(shiftdim(temp2(:,((2*i)-1),yk,xk),3)));

  subplot(5,1,i);
  contourf(tim, lat(yk), tem', 10);
%  colorbar('horiz');
  ylabel(['<'num2str(lims(j,:)) '>, ' num2str(depth((2*i)-1)/100) 'm']);
%  xlabel('time'); 
%  ylabel('latitude');

end
end

%  Time-depth plots

limy = [25 35];
limx = [120 150; 150 180; 180 210; 210 240; 240 270];
limy = flipud([-35 -25; -25 -15; -15 -5; -5 5; 5 15; 15 25; 25 35]);
%limx = flipud([150 170; 150 170; 150 170; 135 170; 120 140; 120 140; 120 140]);
limx = flipud([265 290; 265 290; 255 280; 255 280; 255 280; 225 250; 215 240]);
for j = 1:7;
  [xk, yk] = keep_var([limx(j,:) limy(j,:)], lon, lat);
  tem = squeeze(mean2(mean2(shiftdim(temp2(:,:,yk,xk), 2))));

  figure(i); figure_orient;
  subplot(7,1,j);
  contourf(465+tim, -1*depth/100, tem', 10);
  
  set(gca, 'YTickLabel', -1*str2num(get(gca, 'YTickLabel')));

  ylabel(['<' num2str(limy(j,1)) ' to ' num2str(limy(j,2)) '>']);
end
subplot(7,1,1);
title(['Time - Depth Hofmoeller Diagram for Eastern Boundary']);

% lat: ' num2str(limy(1)) ...
%       ' to ' num2str(limy(2))]);
subplot(7,1,7);
xlabel('Years'); 

cd /home/disk/tao/dvimont/matlab/CSIRO/Plots3