pro newest_center_fit,Date=Date,stats_directory=stats_directory,noab=noab

If (not keyword_set(noab)) then noab=0
If (not keyword_set(Date)) then Date='10'
If (not keyword_set(stats_directory)) then stats_directory='stats'

If noab then begin
stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$
          '/ShackHartman/stats/'
postscript_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$
          '/ShackHartman/distorted_postscripts_no_ab/'
endif

stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$
          '/ShackHartman/stats/'
offset_files=file_search(stats_directory+'*5132*newest_offsets.txt',count=num_stats_files)
readcol,stats_directory+'Delta_offsets.txt',xavgcoord,deltax,yavgcoord,deltay,$
	format='(X,f,f,f,f)'
bad_delta=where(deltax eq -1)

ycenters=fltarr(num_stats_files)
xcenters=fltarr(num_stats_files)

get_lun,newcenters
openw,newcenters,stats_directory+'Distorted_Centers.txt'

for file_number=0,num_stats_files-1 do begin

   readcol,offset_files(file_number),xcoord,ycoord,format='(X,f,X,f,X)'
   x_center_arr=xcoord-deltax
   y_center_arr=ycoord-deltay
   bad_coord=where(xcoord eq -10) 
   x_center_arr(bad_coord)=0
   y_center_arr(bad_coord)=0
   x_center_arr(bad_delta)=0
   y_center_arr(bad_delta)=0
   xcenters(file_number)=avg(x_center_arr(where(x_center_arr ne 0)))
   ycenters(file_number)=avg(y_center_arr(where(y_center_arr ne 0)))
   printf,newcenters,format='(%"%f\t%f")',xcenters(file_number),ycenters(file_number)
   
endfor

free_lun,newcenters
close,newcenters

stop

end