% read the interp.dat file output by sst_set in the RCM.

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

a = fscanf(fid,'%f');
a(a==0) = NaN;
sst = reshape(a(1:199*106),199,106);
x = -130:0.333:-64;
y = -15 :0.333: 20;


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