clear all;clf

% example on how to read and plot viirs 4km data
vfile='V20170602017090.L3m_MO_SNPP_CHL_chl_ocx_4km.nc';
a='''';
info=ncinfo(vfile);
for k=1:length(info.Variables)
    str=[info.Variables(k).Name '=ncread(' a vfile a ',' ...
         a info.Variables(k).Name a ');'];
    eval(str);
end

[x,y]=meshgrid(lon,lat);
st1=ncreadatt(vfile,'/','time_coverage_start');
st2=ncreadatt(vfile,'/','time_coverage_end');
t1=datenum(st1([1:4 6:7 9:10 12:13 15:16 18:19]),'yyyymmddHHMMSS');
t2=datenum(st2([1:4 6:7 9:10 12:13 15:16 18:19]),'yyyymmddHHMMSS');
t=(t1+t2)/2;