Esope source files

   1 : C VCRIT     SOURCE    PV        07/10/05    21:17:04     5890           
   2 :       SUBROUTINE VCRIT(I)
   3 : *  calcul la taille locale
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 : -INC TDEMAIT
   7 :       denf=0
   8 :       denom=0
   9 :       xi=xyz(1,i)
  10 :       yi=xyz(2,i)
  11 :       zi=xyz(3,i)
  12 :       do 219 ipts=1,nptcom
  13 :         dist=sqrt((xyz(1,ipts)-xi)**2+
  14 :      >            (xyz(2,ipts)-yi)**2+
  15 :      >            (xyz(3,ipts)-zi)**2)
  16 :         if (abs(dist).lt.1d-30) dist=1d-30
  17 :         poids=xyz(4,ipts)**2/dist
  18 :         denf=denf+xyz(4,ipts)*poids
  19 :         denom=denom+poids
  20 :  219  continue
  21 :       xmoy =denf/denom
  22 :       XYZ(4,I)=xmoy
  23 :       end
  24 :  
  25 :