      MODULE constants
        IMPLICIT NONE

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