41
42
43
44#include "implicit_f.inc"
45
46
47
48#include "scr03_c.inc"
49#include "vect01_c.inc"
50
51
52
53 INTEGER IXS(NIXS,*), NGL(*), MXT(*),NGEO(*),
54 . IX1(*), IX2(*), IX3(*), IX4(*)
55 double precision
56 . x1(*),x2(*),x3(*),x4(*),
57 . y1(*),y2(*),y3(*),y4(*),
58 . z1(*),z2(*),z3(*),z4(*)
59
61 . x(3,*),xrefs(8,3,*)
62
63
64
65 INTEGER I,N1,N2,N3,N4
66
67
68
71
72
73
74 DO i=lft,llt
75 mxt(i) =ixs(1,i)
76 ngeo(i)=ixs(nixs-1,i)
77 ngl(i) =ixs(nixs,i)
78 ix1(i) =ixs(2,i)
79 ix2(i) =ixs(4,i)
80 ix3(i) =ixs(7,i)
81 ix4(i) =ixs(6,i)
82 ENDDO
83
84 IF (nxref == 0) THEN
85 DO i=lft,llt
87
88 ix2(i)=ixs(6,i)
89 ix4(i)=ixs(4,i)
90 ixs(4,i)=ix2(i)
91 ixs(6,i)=ix4(i)
92 ixs(5,i)=ix2(i)
93 ixs(9,i)=ix4(i)
94 ENDIF
95 x1(i)=x(1,ix1(i))
96 y1(i)=x(2,ix1(i))
97 z1(i)=x(3,ix1(i))
98 x2(i)=x(1,ix2(i))
99 y2(i)=x(2,ix2(i))
100 z2(i)=x(3,ix2(i))
101 x3(i)=x(1,ix3(i))
102 y3(i)=x(2,ix3(i))
103 z3(i)=x(3,ix3(i))
104 x4(i)=x(1,ix4(i))
105 y4(i)=x(2,ix4(i))
106 z4(i)=x(3,ix4(i))
107 ENDDO
108 ELSE
109 DO i=lft,llt
111
112 ix2(i)=ixs(6,i)
113 ix4(i)=ixs(4,i)
114 ixs(4,i)=ix2(i)
115 ixs(6,i)=ix4(i)
116 ixs(5,i)=ix2(i)
117 ixs(9,i)=ix4(i)
118 x1(i) = xrefs(1,1,i)
119 y1(i) = xrefs(1,2,i)
120 z1(i) = xrefs(1,3,i)
121 x2(i) = xrefs(5,1,i)
122 y2(i) = xrefs(5,2,i)
123 z2(i) = xrefs(5,3,i)
124 x3(i) = xrefs(6,1,i)
125 y3(i) = xrefs(6,2,i)
126 z3(i) = xrefs(6,3,i)
127 x4(i) = xrefs(3,1,i)
128 y4(i) = xrefs(3,2,i)
129 z4(i) = xrefs(3,3,i)
130 ELSE
131 x1(i) = xrefs(1,1,i)
132 y1(i) = xrefs(1,2,i)
133 z1(i) = xrefs(1,3,i)
134 x2(i) = xrefs(3,1,i)
135 y2(i) = xrefs(3,2,i)
136 z2(i) = xrefs(3,3,i)
137 x3(i) = xrefs(6,1,i)
138 y3(i) = xrefs(6,2,i)
139 z3(i) = xrefs(6,3,i)
140 x4(i) = xrefs(5,1,i)
141 y4(i) = xrefs(5,2,i)
142 z4(i) = xrefs(5,3,i)
143 ENDIF
144 xrefs(1,1,i) = x1(i)
145 xrefs(1,2,i) = y1(i)
146 xrefs(1,3,i) = z1(i)
147 xrefs(2,1,i) = x2(i)
148 xrefs(2,2,i) = y2(i)
149 xrefs(2,3,i) = z2(i)
150 xrefs(3,1,i) = x3(i)
151 xrefs(3,2,i) = y3(i)
152 xrefs(3,3,i) = z3(i)
153 xrefs(4,1,i) = x4(i)
154 xrefs(4,2,i) = y4(i)
155 xrefs(4,3,i) = z4(i)
156 ENDDO
157 ENDIF
158
159 RETURN
function checkvolume_4n(x, ixs)
1.15.0