Sources Esope | Cast3M

denchp
C DENCHP    SOURCE    CHAT      05/01/12    22:41:53     5004
C calcule une taille de maille en fonction d'un chpoint de densite
      function denchp(ichp,x,y)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
-INC CCREEL
 
-INC PPARAM
-INC CCOPTIO
      real*8 x,y
      segment ichp
        real*8 xr(npc),yr(npc),den(npc)
      endsegment
      denf=0.
      denom=0.
      denm=xpetit
      do 20 i=1,den(/1)
        denm=max(denm,den(i))
  20  continue
      denm=max(1.d0,denm)
      xcorr=den(/1)*denm*1D17
      do 10 i=1,den(/1)
        dist=sqrt((xr(i)-x)**2+(yr(i)-y)**2)
*        if (abs(dist).lt.1d-30) dist=1d-30
        if (abs(dist).lt.xpetit*xcorr) dist=xpetit*xcorr
        denf=denf+den(i)/dist
        denom=denom+1/dist
  10  continue
      denchp=denf/denom
      if (iimpi.ne.0)
     # write (6,*) ' denchp retour ',denf,denom,x,y,den(/1)
      return
      end