module set_precision ! working precision = double precision integer, parameter :: wp = kind(0.0d0) end module set_precision