--- DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/08/31 14:56:52 1.21 +++ DarthVader/TrackerLevel2/src/F77/functionspfa.f 2007/09/04 09:47:49 1.22 @@ -622,8 +622,8 @@ enddo if(DEBUG.EQ.1) $ print*,'pfaeta2 *** warning *** angle out of range: ',angle - if(angle.lt.angL(1))iang=1 - if(angle.gt.angR(nangbin))iang=nangbin + if(angle.le.angL(1))iang=1 + if(angle.ge.angR(nangbin))iang=nangbin 98 continue !jump here if ok @@ -740,8 +740,8 @@ enddo if(DEBUG.EQ.1) $ print*,'pfaeta3 *** warning *** angle out of range: ',angle - if(angle.lt.angL(1))iang=1 - if(angle.gt.angR(nangbin))iang=nangbin + if(angle.le.angL(1))iang=1 + if(angle.ge.angR(nangbin))iang=nangbin 98 continue !jump here if ok @@ -856,8 +856,8 @@ enddo if(DEBUG.EQ.1) $ print*,'pfaeta4 *** warning *** angle out of range: ',angle - if(angle.lt.angL(1))iang=1 - if(angle.gt.angR(nangbin))iang=nangbin + if(angle.le.angL(1))iang=1 + if(angle.ge.angR(nangbin))iang=nangbin 98 continue !jump here if ok @@ -1961,8 +1961,8 @@ enddo if(DEBUG.eq.1) $ print*,'pfacorr *** warning *** angle out of range: ',angle - if(angle.lt.angL(1))iang=1 - if(angle.gt.angR(nangbin))iang=nangbin + if(angle.le.angL(1))iang=1 + if(angle.ge.angR(nangbin))iang=nangbin 98 continue !jump here if ok pfacorr = fcorr(iview,lad,iang)