pro read_dig,file,x,y

Sz=size(file)

if Sz(n_elements(Sz)-2) ne 7 or n_params() ne 3 then message,'Incorrect call'
if file eq '' then message,'Incorrect call'

	openr,lun,file,/get_lun
stat=fstat(lun)
x=(y=fltarr(stat.size/8))
readu,lun,x,y
	free_lun,lun

	end