Télécharger vcrit.eso

Retour à la liste

Numérotation des lignes :

  1. C VCRIT SOURCE JC220346 16/11/29 21:15:40 9221
  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.  
  26.  
  27.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales