diff --git a/src/gray.f b/src/gray.f
index cc31738..0144309 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -660,7 +660,7 @@ c
       use coreprofiles, only : psdbnd,read_profiles,tene_scal,set_prfspl
       use interp_eqprof, only : rmxm,zbmin,zbmax,rcen,psia,rmaxis,
      . zmaxis,rbbbs,zbbbs,read_eqdsk,change_cocos,eq_scal,set_eqspl,
-     . btrcen,psinr=>psiq,qpsi=>q
+     . btrcen,psinr=>psiq,qpsi=>q,equinum_fpol
       use reflections, only : rlim,zlim,nlim,alloc_lim
       use beamdata, only : nrayr,nrayth,nstep
       implicit none
diff --git a/src/interp_eqprof.f90 b/src/interp_eqprof.f90
index e681cae..01733c9 100644
--- a/src/interp_eqprof.f90
+++ b/src/interp_eqprof.f90
@@ -440,7 +440,7 @@ contains
     real(wp_), dimension(1) :: rrs,ffspl
     real(wp_), dimension(nsf) :: wrkfd
 !      
-    if(psinv.le.1.0_wp_.and.psinv.gt.0.0_wp_) then
+    if(psinv.le.1.0_wp_.and.psinv.ge.0.0_wp_) then
       rrs(1)=psinv
       call splev(tfp,nsf,cfp,3,rrs,ffspl,1,ier)
       fpolv=ffspl(1)