C PROCON    SOURCE    PV090527  23/07/13    21:15:04     11708          
       FUNCTION PROCON(XX,YY,NSTRS)

       IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)

       DIMENSION XX(*),YY(*)
       data ISAV /0/
       save ISAV

       data sqxgr /0/
       save sqxgr
-INC CCREEL

       PROCO = 0.D0

       IF(ISAV .EQ. 0)THEN
        sqxgr=sqrt(xgrand*xzprec)
        isav=1
       ENDIF

       IF (NSTRS .GT. 3) GOTO 1
       DO I=1,NSTRS
         XI    = min(sqxgr,max(XX(i),-sqxgr))
         YI    = min(sqxgr,max(YY(i),-sqxgr))
         PROCO = PROCO + XI*YI
       ENDDO
       PROCON = PROCO
       RETURN

 1     DO I=4,NSTRS
         XI    = min(sqxgr,max(XX(i),-sqxgr))
         YI    = min(sqxgr,max(YY(i),-sqxgr))
         PROCO = PROCO + XI*YI
       ENDDO
       PROCO = PROCO * 2.D0

       DO I=1,3
         XI    = min(sqxgr,max(XX(i),-sqxgr))
         YI    = min(sqxgr,max(YY(i),-sqxgr))
         PROCO = PROCO + XI*YI
       ENDDO
       PROCON = PROCO
       END