MODULE cosp_math_constants
  USE cosp_kinds, only: wp
  IMPLICIT NONE
  REAL(wp), PARAMETER :: pi =  acos(-1.0_wp)

END MODULE cosp_math_constants