diff --git a/src/gray.f b/src/gray.f
index 781a860..94acbe8 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -257,7 +257,7 @@ c              ierr=0
 
           if(i.gt.1.and.psinv.ge.0.and.psinv.lt.psdbnd) 
      .     iov(j,k)=1
-          if(iov(j,k).eq.1.and.psinv.ge.1.0d0) iov(j,k)=2
+          if(iov(j,k).eq.1.and.psinv.ge.psdbnd) iov(j,k)=2
 c         iov=0 initially, iov=1 first entrance in plasma, 
 c         iov=2 first exit from plasma, iov=3 after 2nd entrance into plasma