nc = netcdf('sst_fill.nc');
lon = nc{'lon'}(:);
lat = nc{'lat'}(:);
sst = nc{'sst'}(:);
sst(sst==nc{'sst'}.missing_value(:)) = NaN;
am = nc{'available_mask'}(:); % FYI the same every day

M = moviein(size(sst,1));
for j = 1:size(sst,1) % really slow!
  [c,h] = contourf(360-lon,lat,squeeze(sst(j,:,:))',15:0.5:33);
  set(h(:),'linestyle','none');
  title(['SST (C), ' datestr(j+datenum('31-aug-2001'))]);
  colorbar
  set(gca,'xdir','reverse');
  xlabel('degrees W longitude');
  ylabel('degrees N latitude');
  M(:,j) = getframe
end

t = sst;
t(~am) = NaN;
hovsst = squeeze(nanmean(t(:,:,80:100),3));
pcolor(360-lon(1:200),1:61,hovsst(:,1:200));
shading flat;
ylabel('day after 31 Aug 2001');
xlabel('degrees West longitude');
set(gca,'xdir','reverse','tickdir','out')
title('6-11 degrees N mean SST Hovmoller')
% probabaly should 3-day bandpass filter this dataset
% worry also about aliasing of the diurnal cycle on the TRMM orbit.
%
% Justin's SST from 3-day running averages looks similar.  Is there a 6-12
% day cycle? Is it real or aliased? The amplitude is 0.5C, which is on the
% order of what to expect from a diurnal skin temperature variation. For
% now (9/21/04) we ignore it--let the model integrate it!
