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