double complex fnbt
      double precision btlap, btdx, btdz, btkl
      fnbt(btdx,btdz,btkl,btlap)
     &   = (0.0d0, 1.0d0)/( 2.0d0* btkl)*btlap / btdx**2 * btdz