pro read_diff_movie3,indir,file_1,file_2

file_name_1=indir+'/'+file_1
;file_name_2=indir+'/'+file_2

fits_read,file_name_1,im,h,/header_only
row_number=sxpar(h,'NAXIS1')
column_number=sxpar(h,'NAXIS2')
total_frames=sxpar(h,'NAXIS3')
temp=sxpar(h,'MOLYTABL')
bias=sxpar(h,'VSUB')
threshold=1000
rebinsize=4.
numpixrd=row_number*column_number

file_parts=strsplit(file_1,'_',/extract)

if N_elements(file_parts) eq 6 then begin

	file_parts_1=strsplit(file_1,'_',/extract)
;	file_parts_2=strsplit(file_2,'_',/extract)

	file_end_1=strsplit(file_parts_1(5),'.',/extract)
;	file_end_2=strsplit(file_parts_2(5),'.',/extract)

	file_num_1=file_end_1(0)
;	file_num_2=file_end_2(0)
	
	tempstr=file_parts_1(1)
	
endif else begin

	file_parts_1=strsplit(file_1,'_',/extract)
 ;       file_parts_2=strsplit(file_2,'_',/extract)

        file_num_1=file_parts_1(4)
  ;      file_num_2=file_parts_2(4)

	tempstr=file_parts(1)

endelse 

Dir_parts=strsplit(indir,'/',/extract)

Dir=Dir_parts(4)+'_'+Dir_parts(6)

;Create MPEG Object to generate movie
poofy_movie=mpeg_open([row_number/rebinsize,column_number/rebinsize],quality=10)

for read_num=0,total_frames-2 do begin

print,'processing frame #',strtrim(string(read_num+1),2)
;Take in section of one file as the vector x and the section of another file
;as the vector x2
fits_read,file_name_1,x1,first=long(read_num*numpixrd),last=long((read_num+1)*numpixrd)-1
fits_read,file_name_1,x2,first=long((read_num+1)*numpixrd),last=long((read_num+2)*numpixrd)-1

;Since x and x2 are one-dimensional vectors, we need to turn them into
;2-dimensional arrays where the size of each dimension is size of dimension in device
im1=reform(float(x1),row_number,column_number)
im2=reform(float(x2),row_number,column_number)
diff=float(im2-im1)
diff=bytscl(congrid(diff,1024,1024),min=-threshold,max=threshold)

for i=0,30 do begin & mpeg_put,poofy_movie,IMAGE=diff,FRAME=read_num+i &
endfor

endfor

mpeg_save,poofy_movie,FILENAME=Dir+'_'+tempstr+'_'+file_num_1+'.mpg'
mpeg_close,poofy_movie

end