%script to check time base of raw MIDAS burst data

nc = netcdf('930TBST.CDF');
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);
i450 = find(tdifs==450);
next450 = i450+1;
t450 = base + time(next450)./86400;
str450 = datestr(t450,0)