Télécharger denchp.eso

Retour à la liste

Numérotation des lignes :

  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. -INC CCOPTIO
  8. real*8 x,y
  9. segment ichp
  10. real*8 xr(npc),yr(npc),den(npc)
  11. endsegment
  12. denf=0.
  13. denom=0.
  14. denm=xpetit
  15. do 20 i=1,den(/1)
  16. denm=max(denm,den(i))
  17. 20 continue
  18. denm=max(1.d0,denm)
  19. xcorr=den(/1)*denm*1D17
  20. do 10 i=1,den(/1)
  21. dist=sqrt((xr(i)-x)**2+(yr(i)-y)**2)
  22. * if (abs(dist).lt.1d-30) dist=1d-30
  23. if (abs(dist).lt.xpetit*xcorr) dist=xpetit*xcorr
  24. denf=denf+den(i)/dist
  25. denom=denom+1/dist
  26. 10 continue
  27. denchp=denf/denom
  28. if (iimpi.ne.0)
  29. # write (6,*) ' denchp retour ',denf,denom,x,y,den(/1)
  30. return
  31. end
  32.  
  33.  
  34.  
  35.  
  36.  

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