Télécharger denchp.eso

Retour à la liste

Numérotation des lignes :

denchp
  1. C DENCHP SOURCE CHAT 05/01/12 22:41:53 5004
  2. C calcule une taille de maille en fonction d'un chpoint de densite
  3. function denchp(ichp,x,y)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. -INC CCREEL
  7.  
  8. -INC PPARAM
  9. -INC CCOPTIO
  10. real*8 x,y
  11. segment ichp
  12. real*8 xr(npc),yr(npc),den(npc)
  13. endsegment
  14. denf=0.
  15. denom=0.
  16. denm=xpetit
  17. do 20 i=1,den(/1)
  18. denm=max(denm,den(i))
  19. 20 continue
  20. denm=max(1.d0,denm)
  21. xcorr=den(/1)*denm*1D17
  22. do 10 i=1,den(/1)
  23. dist=sqrt((xr(i)-x)**2+(yr(i)-y)**2)
  24. * if (abs(dist).lt.1d-30) dist=1d-30
  25. if (abs(dist).lt.xpetit*xcorr) dist=xpetit*xcorr
  26. denf=denf+den(i)/dist
  27. denom=denom+1/dist
  28. 10 continue
  29. denchp=denf/denom
  30. if (iimpi.ne.0)
  31. # write (6,*) ' denchp retour ',denf,denom,x,y,den(/1)
  32. return
  33. end
  34.  
  35.  
  36.  
  37.  
  38.  

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