%For comparison between HWRF POM output

clear all;close all
%load('eastcoastmask_new.mat')

%file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110826T070000-avhrr_3day1200to1500-coldcomp.nc';
%file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110827T070000-avhrr_3day1200to1500-coldcomp.nc';
%file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110828T070000-avhrr_3day1200to1400-coldcomp.nc';
%file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110829T070000-avhrr_3day1200to1400-coldcomp.nc';
%file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110830T070000-avhrr_3day1300to1400-coldcomp.nc';
file_nc     = '/Volumes/boardwalk/seroka/data/SST/rutgers/irene/sport_20110831T070000-avhrr_3day1200to1400-coldcomp.nc';
xgrid       = nc_varget(file_nc,'lon');
ygrid       = nc_varget(file_nc,'lat');
sst         = nc_varget(file_nc,'SST');

figure;
set(gcf,'Paperposition',[0 0 6 6]);
pcolor(xgrid,ygrid,sst); shading flat; caxis([14 30]);
axis([-77.5 -68 34 43])
colormap(jet); hc=colorbar; set(hc,'tickdir','out');
hold on
grid

%plot coastline
%plot(lon,lat,'k')

set(gca,'tickdir','out','fontsize',12);
ylabel('Latitude (^oN)');
xlabel('Longitude (^oW)');

title('AVHRR+NASA SPoRT SST 20110831 07Z','fontsize',12);

cbfreeze(hc)
freezeColors
MAB_no_contour_generator(1)
grid on;
set(gca,'layer','top')
    
project_mercator
wysiwyg

eval('print -dpng /Users/seroka/Research/data/HWRF/sport_20110831T070000.png');
