subroutine cplogp(wvec,wt,ne)

  implicit none

! arguments:
  real, dimension(:) :: wvec, wt
  integer ne

  wt(1:ne) = log10(abs(wvec(1:ne)) + 1.e-15)
  
end subroutine cplogp