;-------------------------------------------------------------------- ; Interpolates average time between exposures for MOSES-I dataset. ; Compares average cadence to actual time between exposures. If ; actual time is within +-(200ms) of average cadence, returns zero. ; If actual time deviates more than +-(200ms) from average cadence ; returns deviation in seconds. Uses files exptime.pro and index.sav ; Created on 02 June 2009. PRO cadence_int, comp = comp ; Call exptime.pro to get exposure times and lengths exptime, exptime = a, explength = b ; Determine average cadence t = (a(1,26)-total(b))/26 ; Find actual cadence cadence = dblarr(26) for i=0,25 do begin cadence(i) = a(0,i+1) - a(1,i) end ; Compare average cadence to actual cadence comp = fltarr(26) for i=0,25 do begin IF (cadence(i) gt t-.2 and cadence(i) lt t+.2) THEN BEGIN comp(i) = 0 ENDIF ELSE BEGIN comp(i) = cadence(i)-t ENDELSE end ; Output results print, 'Average Cadence: ',t,' seconds' print, ' ' print, 'If zero, then actual cadence is +-200ms from average cadence' print, 'If nonzero, then gives deviation from average cadence in seconds' print, comp end