MODULE constants IMPLICIT NONE ! ** Double precision (i.e accuracy) ** INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(12) END MODULE constants