Documentation of hp_lp_SST_hoffmoeller

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

Function Synopsis

maxt = hp_lp(t1);

Help text

if i == 1;

Cross-Reference Information

Listing of function hp_lp_SST_hoffmoeller

function maxt = hp_lp(t1);
get_cdtem;
cdir /home/disk/hayes2/dvimont/csiro/data/Individual_levels
figure(1); fo; clf;
for i = 1:2;
  if i == 1;
    filin = 'temp_M_L1_1000_years_new.nc';
  elseif i == 2;
    filin = 'temp_M_L5_1000_years_new.nc';
  end
lims = [100 300 -5 5];
tim = t1 + [1:600];
nc = netcdf(filin, 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  [xk, yk] = keep_var(lims, lon, lat);
  temp = nc{'temp'}(tim, 1, yk, xk);
  mv = nc{'temp'}.missing_value(:);
nc = close(nc);
temp(temp == mv) = NaN;
temp = squeeze(temp);
temp = squeeze(mean2(shiftdim(temp, 1)));
temp = temp';
[temp, clim] = annave(temp);
temp = detrend(temp);
temp = rave(temp, 5);
maxt(i) = max([max(max(temp)) -1*min(min(temp))]);
tind = [1:600];
  subplot(1,2,i);
  pcolor(lon(xk), tind, temp(tind,:));
  shading interp
  axis([100 280 min(tind)-1 max(tind)+1]);
  set(gca, 'YTick', 0:60:12000, 'YTickLabel', (t1/12):5:1000, ...
	   'XTick', 0:30:360)
  caxis([-maxt(i) maxt(i)]);
%elseif i == 2;
%  subplot(1,2,2);
%  pcolor(lon(xk), tind, temp(tind,:));
%  shading interp
%  axis([100 280 min(tind)-1 max(tind)+1]);
%  set(gca, 'YTick', 0:60:12000, 'YTickLabel', (t1/12):5:1000, ... 
%	   'XTick', 0:30:360)
%end

end

back

end