Documentation of thermocline_structure

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

Help text

  ctlim = [-.1 360 -35 -20];
  ctlim = [-.1 360 22.5 37.5];
  ctlim = [165 210 22.5 37.5];

Cross-Reference Information

Listing of script thermocline_structure

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'}(:);
  [xk, yk] = keep_var(ctlim, lon, lat);

  temp = nc{'temp'}(:,:,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)));

top = squeeze(shiftdim(mean(shiftdim(temp, 2)), 2));

cint = 1;
clims = [6:cint:29];

colormap('default');

top2 = squeeze(mean(top));
top2 = [top2(:,3:length(lon)) top2(:,1:3)];
lon2 = [lon(3:length(lon)); lon(1:3)+360];

figure(1); figure_landscape; clf;
sp(1)
%     a = contourf(lon2, -.01*depth, top2, [5:30]);
     pncont(lon2, -.01*depth, top2, [-2:30], 1, 'k');
%     clabel(a, 'manual');
     grid on;
     set(gca, 'XTick', 0:30:360, 'YTick', -550:50:0, 'YTickLabel', 550:-50:0)
     axis([10 375 -545 -5])
     title('Mean Vertical Temperature Profile, Averaged from 5s to 5n');
     xlabel(['Contour Interval:  ' num2str(cint) ' degree C'])
     ylabel(['Depth in meters']);
%     colorbar2

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

top = reshape(top, ntim, nlev*nlon);
top = detrend(top);

bot = ctlow' * top / ntim;
bot = reshape(bot, nlev, nlon);
bot = [bot(:,3:length(lon)) bot(:,1:3)];

sp(2)
%     a = contourf(lon2, -.01*depth, bot, [-.2:.05:.35]);
%     clabel(a, 'manual');
     pncont(lon2, -.01*depth, bot, [-.2:.05:.35], 1, 'k');
     grid on;
     axis([10 375 -545 -5])
     set(gca, 'XTick', 0:30:360, 'YTick', -550:50:0, 'YTickLabel', 550:-50:0)
     title('Temperature Profile Regressed on CTLP, Averaged from 22.5n to 37.5n');
     xlabel(['Contour Interval:  0.05 degree C'])
     ylabel(['Depth in meters']);
%     colorbar2

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

l = 7;
bot = ctlow(1:(ntim-l+1))' * top(l:ntim,:) / (ntim-l);
bot = reshape(bot, nlev, nlon);
bot = [bot(:,3:length(lon)) bot(:,1:3)];