% read the formatted file for model interpolated SST.

fid = fopen('interp.dat','r');

a = fscanf(fid,'%f',199*106);
sst = reshape(a,199,106);
x = 230:0.333:296;
y = -15:0.333:20;

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