%script to check time base of raw MIDAS burst data

nc = netcdf
time = nc{'time'}(:);
yr = nc{'time'}.base_date(1);
mo = nc{'time'}.base_date(2);
da = nc{'time'}.base_date(3);
hr = nc{'time'}.base_date(4);
mi = nc{'time'}.base_date(5);
sc = nc{'time'}.base_date(6);
base = datenum(yr,mo,da,hr,mi,sc);
close (nc)

tdifs = diff(time);
[theValues, theCounts]=abundance(tdifs);
plot(theValues, theCounts,'o')
zoomsafe   

fprintf('Gap    \t        \tRecords before and after first gap this size\n')
fprintf('size(s)\tCount   \trecord no. \t           \ttimes\n')

for i = 1:length(theValues)
   thisdif = find (tdifs == theValues(i));
   rec_before(i) = thisdif(1);
	time_before(i) = base + time(rec_before(i))./86400;
   str_before = datestr(time_before(i),0);
   rec_after(i) = thisdif(1) +1;
	time_after(i) = base + time(rec_after(i))./86400;
   str_after = datestr(time_after(i),0);
   fprintf('%7.1f\t%8.1f\t%8.1f\t%8.1f\t%20s\t%20s\n',...
      theValues(i),theCounts(i),rec_before(i),rec_after(i),...
      str_before,str_after)
end