diff --git a/src/gray.f b/src/gray.f
index 3e3296b..7b4b44e 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -6295,7 +6295,7 @@ c     of gaussian profile
      
       rhpp=frhopol(rhotpav)
       rhpj=frhopol(rhotjava)
-      dpdvp=pabs*2.0d0/(sqrt(pi)*drhotjava*fdvdrhot(rhpp))
+      dpdvp=pabs*2.0d0/(sqrt(pi)*drhotpav*fdvdrhot(rhpp))
       ajphip=currt*2.0d0/(sqrt(pi)*drhotjava*fdadrhot(rhpj))
       call ratioj(rhpj,ratjamx,ratjbmx,ratjplmx)