Imaginary unit was not written in double precision

src/prec_const_mod.F90

@@ -23,31 +23,31 @@ MODULE prec_const
 
 
   ! Some useful constants, to avoid recomputing then too often
-  REAL(dp), PARAMETER :: PI=3.141592653589793238462643383279502884197_dp
-  REAL(dp), PARAMETER :: PIO2=1.57079632679489661923132169163975144209858_dp
-  REAL(dp), PARAMETER :: TWOPI=6.283185307179586476925286766559005768394_dp
-  REAL(dp), PARAMETER :: SQRT2=1.41421356237309504880168872420969807856967_dp
-  REAL(dp), PARAMETER :: INVSQRT2=0.7071067811865475244008443621048490392848359377_dp
-  REAL(dp), PARAMETER :: SQRT3=1.73205080756887729352744634150587236694281_dp
-  REAL(dp), PARAMETER :: onetwelfth=0.08333333333333333333333333333333333333333333333_dp
-  REAL(dp), PARAMETER :: onetwentyfourth=0.04166666666666666666666666666666666666666666667_dp
-  REAL(dp), PARAMETER :: onethird=0.33333333333333333333333333333333333333333333333_dp
-  REAL(dp), PARAMETER :: onesixth=0.1666666666666666666666666666666666666666666667_dp
-  REAL(dp), PARAMETER :: fivesixths=0.8333333333333333333333333333333333333333333333_dp
-  REAL(dp), PARAMETER :: sevensixths=1.1666666666666666666666666666666666666666666667_dp
-  REAL(dp), PARAMETER :: elevensixths=1.833333333333333333333333333333333333333333333_dp
-  REAL(dp), PARAMETER :: nineeighths=1.125_dp
-  REAL(dp), PARAMETER :: onesixteenth=0.0625_dp
-  REAL(dp), PARAMETER :: ninesixteenths=0.5625_dp
-  REAL(dp), PARAMETER :: thirteentwelfths = 1.083333333333333333333333333333333333333333333_dp
-  COMPLEX(dp),  PARAMETER :: imagu = (0,1)
+  REAL(dp),    PARAMETER :: PI=3.141592653589793238462643383279502884197_dp
+  REAL(dp),    PARAMETER :: PIO2=1.57079632679489661923132169163975144209858_dp
+  REAL(dp),    PARAMETER :: TWOPI=6.283185307179586476925286766559005768394_dp
+  REAL(dp),    PARAMETER :: SQRT2=1.41421356237309504880168872420969807856967_dp
+  REAL(dp),    PARAMETER :: INVSQRT2=0.7071067811865475244008443621048490392848359377_dp
+  REAL(dp),    PARAMETER :: SQRT3=1.73205080756887729352744634150587236694281_dp
+  REAL(dp),    PARAMETER :: onetwelfth=0.08333333333333333333333333333333333333333333333_dp
+  REAL(dp),    PARAMETER :: onetwentyfourth=0.04166666666666666666666666666666666666666666667_dp
+  REAL(dp),    PARAMETER :: onethird=0.33333333333333333333333333333333333333333333333_dp
+  REAL(dp),    PARAMETER :: onesixth=0.1666666666666666666666666666666666666666666667_dp
+  REAL(dp),    PARAMETER :: fivesixths=0.8333333333333333333333333333333333333333333333_dp
+  REAL(dp),    PARAMETER :: sevensixths=1.1666666666666666666666666666666666666666666667_dp
+  REAL(dp),    PARAMETER :: elevensixths=1.833333333333333333333333333333333333333333333_dp
+  REAL(dp),    PARAMETER :: nineeighths=1.125_dp
+  REAL(dp),    PARAMETER :: onesixteenth=0.0625_dp
+  REAL(dp),    PARAMETER :: ninesixteenths=0.5625_dp
+  REAL(dp),    PARAMETER :: thirteentwelfths = 1.083333333333333333333333333333333333333333333_dp
+  COMPLEX(dp), PARAMETER :: imagu = (0._dp,1._dp)
   CONTAINS
 
     SUBROUTINE INIT_PREC_CONST
 
       IMPLICIT NONE
       integer :: ierr,me
-      
+
       REAL(sp) :: a = 1_sp
       REAL(dp) :: b = 1_dp