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