31 use element_mod , only : nixr
32
33
34
35#include "implicit_f.inc"
36
37
38
39#include "vect01_c.inc"
40#include "param_c.inc"
41
42
43
44 INTEGER IX(NIXR,*),IGEO(NPROPGI,*)
46 . off(*), geo(npropg,*), x(3,*), al(*),ipos(5,*)
47
48
49
50 INTEGER , J, I1, I2, USENS
51
53 . ex, ey, ez
54
55 DO i=lft,llt
56 j=i+nft
57 usens=igeo(3,ix(1,j))
58 IF (usens <= 0) THEN
59
60 off(i)=one
61 ELSE
62 off(i)=-ten
63 ENDIF
64 ENDDO
65
66 DO j=1,5
67 DO i=lft,llt
68 ipos(j,i)= zero
69 ENDDO
70 ENDDO
71
72 DO i=lft,llt
73 j=i+nft
74 i1=ix(2,j)
75 i2=ix(3,j)
76 ex=x(1,i2)-x(1,i1)
77 ey=x(2,i2)-x(2,i1)
78 ez=x(3,i2)-x(3,i1)
79 al(i)=sqrt(ex*ex+ey*ey+ez*ez)
80 ENDDO
81
82 RETURN