subroutine curcor(n,curv,v,sl,pei,dcut)
*     *****************************************************
*     Correct the curvature dependence of the burning
*     velocity
*     *****************************************************
      integer n,j
      double precision sl,pei,cj,dcut
      double precision curv(0:8192),v(0:8192)
      do 10 j=1,n
         cj = pei*curv(j)
         v(j) = v(j) - sl*(dexp(-cj)-1.D0+cj)
 10   continue
      v(0) = v(n)
      
      call kfilter(n,dcut,v)
   
      end