Esope source files

   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 :