filename='mw_ir.fusion.2015.239.rt';
ia=findstr('mw_ir', filename);
if length(ia) < 1,
      xdim=1440
      ydim=720
else
      xdim=4096
      ydim=2048
end;
xcell=[1:xdim];
ycell=[1:ydim];
xmax=max(xcell)
ymax=max(ycell)
dx=360./xmax
dy=180./ymax
xx  = dx*xcell-dx/2.    
yy  = dy*ycell-(90+dy/2.) 

centigrade_scale_sst   = .15;
centigrade_offset_sst  = -3;
centigrade_scale_error = .005;
centigrade_offset_error= 0;


fid=fopen(filename,'r');
sst = fread(fid,[xdim ydim], 'uchar');
error = fread(fid,[xdim ydim], 'uchar');
mask = fread(fid,[xdim ydim], 'uchar');
fclose(fid);
%scale SST
good = find(sst <= 250);
sst(good) = (sst(good) * centigrade_scale_sst) + centigrade_offset_sst;

%scale error
good = find(error <= 250);
%set flags in SST and error fields using the mask data
land=find(bitget(mask,1)==1);sst(land)=255;error(land)=255;
ice =find(bitget(mask,2)==1);sst(ice)=252;error(ice)=252;
bad =find(bitget(mask,5)==1);sst(bad)=254;error(bad)=254;
yyn=fliplr(yy);
figure(1);
subplot(1,1,1);

imagesc(xx,fliplr(yy),sst',[-2,34])