openr,1,'plbs.br.kRA'&s=''&readf,1,s&print,s&d=fltarr(6,13)
readf,1,d&close,1
fac=fltarr(13)&bl=fac&dbl=bl
w=transpose(d(0,*)*10.)
for i=0,12 do begin
    fac(i)=.2651*d(0,i)/(d(3,i)*d(4,i)*d(5,i))
    bl(i)=fac(i)*(d(1,i)-d(2,i))
    dbl(i)=fac(i)*sqrt(2)*.001
endfor
end