function marked_times(maskfile,datafile,variable) %output times of marks in of EPIC mask netcdf's nc = netcdf(maskfile) marks = nc{variable}(:); close (nc) nc = netcdf(datafile) time = nc{'time'}(:); time2 = nc{'time2'}(:); close (nc) fprintf('Mark no.\tRec. no\tDate and Time\n') whichrecs = find (marks); for i = 1:length(whichrecs) time_num = ep_datenum([time(whichrecs(i)) time2(whichrecs(i))]); time_string = datestr(time_num,0); fprintf('%8.1f\t%7.1f\t%20s\n',i,whichrecs(i),time_string) end