diff --git a/fvn_misc/fvn_misc.f90 b/fvn_misc/fvn_misc.f90
index 01fc571..e8ee974 100644
--- a/fvn_misc/fvn_misc.f90
+++ b/fvn_misc/fvn_misc.f90
@@ -111,11 +111,11 @@ subroutine fvn_z_muller (f,eps,eps1,kn,nguess,n,x,itmax,infer,ier)
       integer :: infer(kn+n)
 
 
-      data                zero/(0.0,0.0)/,p1/(0.1,0.0)/, &
-                          one/(1.0,0.0)/,four/(4.0,0.0)/, &
-                          p5/(0.5,0.0)/, &
-                          rzero/0.0/,rten/10.0/,rhun/100.0/, &
-                          ax/0.1/,ickmax/3/,rp01/0.01/
+      data                zero/(0.0d0,0.0d0)/,p1/(0.1d0,0.0d0)/, &
+                          one/(1.0d0,0.0d0)/,four/(4.0d0,0.0d0)/, &
+                          p5/(0.5d0,0.0d0)/, &
+                          rzero/0.0d0/,rten/10.0d0/,rhun/100.0d0/, &
+                          ax/0.1d0/,ickmax/3/,rp01/0.01d0/
 
             ier = 0
             if (n .lt. 1) then ! What the hell are doing here then ...
@@ -184,7 +184,7 @@ icloop:         do
                         bi = t1-t2+xx+t3
                         den = bi*bi-four*(xx*t1-t3*(t2-xx))
 !                                  use denominator of maximum amplitude 
-                        t1 = cdsqrt(den)
+                        t1 = sqrt(den)
                         qz = rhun*max(cdabs(bi),cdabs(t1))
                         t2 = bi + t1
                         tpq = cdabs(t2)+qz
@@ -313,4 +313,4 @@ icloop:         do
 
 
 
-end module fvn_misc
\ No newline at end of file
+end module fvn_misc