Documentation of vertical_ocean_structure

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

Help text

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

Cross-Reference Information

Listing of script vertical_ocean_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'}(:);
  ctlim = [120 290 -5 45];
  [xk, yk] = keep_var(ctlim, lon, lat);
  tim = 301:400;

  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);

top = squeeze(mean2(mean2(shiftdim(temp, 2))));
top = detrend(top)';
clims = [-.6:.1:.6]
top(10,1:100:1000) = max(clims)*ones(1,10);
top(10,2:100:1000) = min(clims)*ones(1,10);

figure(1);  figure_landscape;
for i = 1:3
  subplot(3,1,i);
     tim = (i+2)*100 + [1:100];
     contourf(tim, -1*depth, top(:,tim),clims)
     grid on
     colorbar2
end

[npx, npy] = keep_var([165 210 22.5 37.5], lon, lat);
[wbx, wby] = keep_var([120 165 5 20], lon, lat);
[eqx, eqy] = keep_var([170 230 -5 5], lon, lat);

nptemp = detrend(squeeze(mean2(mean2(shiftdim(temp(:,:,npy,npx), 2)))))';
wbtemp = detrend(squeeze(mean2(mean2(shiftdim(temp(:,:,wby,wbx), 2)))))';
eqtemp = detrend(squeeze(mean2(mean2(shiftdim(temp(:,:,eqy,eqx), 2)))))';

clims = [-1.5:0.25:1.125];
nptemp(10,1:100:1000) = min(clims)*ones(1,10);
nptemp(10,2:100:1000) = max(clims)*ones(1,10);
wbtemp(10,1:100:1000) = min(clims)*ones(1,10);
wbtemp(10,2:100:1000) = max(clims)*ones(1,10);
eqtemp(10,1:100:1000) = min(clims)*ones(1,10);
eqtemp(10,2:100:1000) = max(clims)*ones(1,10);

figure(1);  figure_landscape;
%for i = 1:2:9;

i = 3;
tim = i*100 + [1:100];
subplot(3,1,1);
  contourf(tim, -0.01*depth, nptemp(:,tim), clims);
  colorbar2
subplot(3,1,2);
  contourf(tim, -0.01*depth, wbtemp(:,tim), clims);
  colorbar2
subplot(3,1,3);
  contourf(tim, -0.01*depth, eqtemp(:,tim), clims);
  colorbar2