;CP_get_times, Directory_file, indir
;
;INPUTS:
;	Directory_file		the file containing the names of all the 
;				fits files for a given date
;	indir			the directory where the .fits files are stored
;	
;
;

pro get_timestamp, Directory_file=Directory_file, $
	       Readdir=Readdir, Writedir=Writedir,$
	       date=date

;Time will be coded as 2005*1000+Day_Number+hr/24+min/60+second/3600
year_offset=double(2005000)
time_gap=0.0001
num_of_telescope_shifts=0

;set directory path delimiter
    if (!version.os_family eq 'unix') then begin
       delim='/'
    endif else begin
       delim='\'
    endelse

;Defaults for keywords if not entered at command line
if (not keyword_set(Date)) then Date='10'  
if (not keyword_set(Rawdir)) then $
Rawdir='/nfs/slac/g/ki/ki08/lsst/CPanalysis/2005-05-'+Date+'/ShackHartman/raw/'
if (not keyword_set(Writedir)) then $
Writedir='/nfs/slac/g/ki/ki08/lsst/CPanalysis/2005-05-'+Date+'/ShackHartman/'
if (not keyword_set(Directory_file)) then Directory_file=Writedir+'Directory_files.txt'

rawfiles=file_search(Rawdir+'*5132*_raw.fits',count=num_of_files)
                                                                                
;The files are in Writedir are listed in the file Directory_files.txt
time_array=dblarr(num_of_files)
uttime_array=strarr(num_of_files)

get_lun,time_file
openw,time_file,Writedir+'/5-'+Date+'-Time_stamp.txt'

;Create a timestap for each of the files and write it to a time file
hours=fltarr(num_of_files)
minutes=fltarr(num_of_files)
seconds=fltarr(num_of_files)


for i=0,num_of_files-1 do begin

	file_name=rawfiles(i)
	fits_read,file_name,image,header,/header_only
	uttime=SXPAR(header,'DATE-OBS')
        time_arr=strsplit(uttime,':',/extract)
        date_arr=strsplit(uttime,'-',/extract)
        time2_arr=strsplit(time_arr(0),'T',/extract)
 	date2_arr=strsplit(date_arr(2),'T',/extract)
        printf,time_file,FORMAT='(%"%d\t%d\t%d\t%d\t%d\t%d\t%s")',$
	       date_arr(0),date_arr(1),date2_arr(0),$
	       time2_arr(1),time_arr(1),time_arr(2),$
	       file_name

	time_with_year=double(DATE_CONV(uttime,'R'))
	time_wo_year=time_with_year-year_offset
	time_array(i)=time_wo_year
	uttime_array(i)=uttime

;If the lapse in time between two successive images is greater than 
;time_gap, the telescope pointing has probably changed

endfor

close,time_file
free_lun,time_file

stop

end