diff --git a/src/density.irp.f b/src/density.irp.f
index 86f5669..9361f57 100644
--- a/src/density.irp.f
+++ b/src/density.irp.f
@@ -53,6 +53,8 @@ BEGIN_PROVIDER [ real, density_alpha_value_p ]
    exc(4) = exc(1)*exc(2)
    if (exc(4) /= 0) then
     exc(4) = exc(4)/abs(exc(4))
+   else
+    exc(4) = 1
    endif
 
    phase = dble(exc(4))
@@ -100,6 +102,8 @@ BEGIN_PROVIDER [ real, density_beta_value_p ]
    exc(4) = exc(1)*exc(2)
    if (exc(4) /= 0) then
     exc(4) = exc(4)/abs(exc(4))
+   else
+    exc(4) = 1
    endif
 
    phase = dble(exc(4))
diff --git a/src/test_1d.irp.f b/src/test_1d.irp.f
index 3873c92..2fe848a 100644
--- a/src/test_1d.irp.f
+++ b/src/test_1d.irp.f
@@ -13,6 +13,7 @@ subroutine run
  do i=-40,40
   point(3) = real(i)/20.
   TOUCH point
-  print *,  point(3), eplf_value_p, eplf_up_up, eplf_up_dn
+! print *,  point(3), eplf_value_p, eplf_up_up, eplf_up_dn
+  print *,  point(3), density_alpha_value_p
  enddo
 end