% read the binary file for daily 2001 SST.

fid = fopen('EPIC_oi_binary.dat','r','ieee-be');

nrec = 61;
nlon = 276;
nlat = 148;

sst = NaN+zeros(nlon,nlat,nrec);
len = fread(fid,1,'uint32');

for rec = 1:nrec;
  a = fread(fid,len/4,'float');
  sst(:,:,rec) = reshape(a,nlon,nlat);
end

x = -130.875:0.25:-62.125;
y = -15.875:0.25:20.875;

% PLOT
pcolor(-x,y,squeeze(sst(:,:,3))'); shading flat
colorbar
set(gca,'xdir','reverse','tickdir','out');
ylabel('latitude')
xlabel('W longitude')
title('SST (K)')
axis equal; axis tight
