lonc =   0.125:0.25:359.875; % center of the patches
latc = -89.875:0.25:89.875;
lonrange = [0 360];
latrange = [-40 40];
lon = lonc(lonc>=lonrange(1) & lonc<=lonrange(2));
lat = latc(latc>=latrange(1) & latc<=latrange(2));

start = datenum('1-may-2001')-datenum('0-jan-2001');
endt  = datenum('31-oct-2001')-datenum('0-jan-2001');

for jday = start:endt
  fname = ['/home/sol/deszoeke/data/oi_sst/tmi/2001/tmi.fusion.2001.',sprintf('%03d',jday),'.v01'];
  fid = fopen(fname);
  temp = reshape(fread(fid,1440*720,'char'),1440,720);
  t = temp(:,latc>=latrange(1) & latc<=latrange(2))
  t(t>250) = NaN;
end