C VCRIT     SOURCE    JC220346  16/11/29    21:15:40     9221           
      SUBROUTINE VCRIT(I)
*  calcul la taille locale
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
-INC TDEMAIT
      denf=0
      denom=0
      xi=xyz(1,i)
      yi=xyz(2,i)
      zi=xyz(3,i)
      do 219 ipts=1,nptcom
        dist=sqrt((xyz(1,ipts)-xi)**2+
     >            (xyz(2,ipts)-yi)**2+
     >            (xyz(3,ipts)-zi)**2)
        if (abs(dist).lt.1d-30) dist=1d-30
        poids=xyz(4,ipts)**2/dist
        denf=denf+xyz(4,ipts)*poids
        denom=denom+poids
 219  continue
      xmoy =denf/denom
      XYZ(4,I)=xmoy
      end