module prec
!M=========
  integer, parameter :: wp = kind(1.0d0)
end module prec

