This module defines single/double precision kind parameters and sets the accuracy to which the floating point operations are to be done throughout the entire package. The user may change the precision of floating point operations by simply editing the WP parameter:
integer, parameter :: WP = DP ! working precision definition
The default value is set to correspond to the obsolescent Fortran 77's double precision data type.