clear


% Read Scientific_Data_Sets

fname='A2006275183500.L2_LAC_SST';



longitude=hdfread(fname,'longitude');
latitude=hdfread(fname,'latitude');
sst=hdfread(fname,'sst');
qual_sst=hdfread(fname,'qual_sst');
cntl_pt_cols=hdfread(fname,'cntl_pt_cols');

lonsat=double(longitude);
latsat=double(latitude);
SST=.005*double(sst(:,cntl_pt_cols));  % down sample to control points as well as add slope of .005. 
flag=double(qual_sst(:,cntl_pt_cols));  % down sample to control points



[mm,nn]=size(SST);

ind=find(flag > 2);  % keep only quality flags of 0,1,2 (3 and 4) denote cloud and land contaminated. 
SST(ind)=NaN;

% flag==4 only filters land
% flag==3 only filters cloud




% determine latitude and longitude limits. 

minlat1=min(min(latsat));
maxlat1=max(max(latsat));
minlon1=min(min(lonsat));
maxlon1=max(max(lonsat));
latlim1=[minlat1 maxlat1];
lonlim1=[minlon1 maxlon1];


% lay on a mercator projection map with limits defined above 
axesm('mercator','MapLatLimit',latlim1,'MapLonLimit',lonlim1,'MeridianLabel','on', ...
    'ParallelLabel','on');

% place on  map axis. 
surfm(latsat,lonsat,SST); colorbar

hold on

% loading and plotting coastlines

load coast;
latc=lat;
lonc=long;

symbols = makesymbolspec('Line',{'Default','Color','black'})
geoshow(latc,lonc,'SymbolSpec',symbols)


% loading state boundaries.

load lons 
load lats

symbols = makesymbolspec('Line',{'Default','Color','black'})
geoshow(lats,lons,'SymbolSpec',symbols)