!-------------------------- read table of P(k)
PROGRAM READTABLE
CHARACTER :: LINE*120
real      :: dk=0.


Open(1,file='pk_cambWMAP.dat')

read(1,'(a)') LINE

   read(1,*,iostat=iFlag)x0,Pk0

do
   read(1,*,iostat=iFlag)x,Pk
   If(iFlag.ne.0)Exit
   if(dk==0.)Then
         dk =x/x0
   else
         dd =x/x0
         write(*,'(5g12.4)') x,dd,dk,dd/dk
   endif
         x0 = x
enddo


stop
end PROGRAM READTABLE