OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i7dst3.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| i7dst3 ../engine/source/interfaces/int07/i7dst3.F
25!||--- called by ------------------------------------------------------
26!|| i7main_lmult ../engine/source/interfaces/int07/i7main_lmult.F
27!|| i7mainf ../engine/source/interfaces/int07/i7mainf.F
28!||--- calls -----------------------------------------------------
29!|| i7cmaj ../engine/source/interfaces/int07/i7curv.f
30!||====================================================================
31 SUBROUTINE i7dst3(
32 1 JLT ,CAND_N ,CAND_E ,CN_LOC ,CE_LOC ,
33 2 X1 ,X2 ,X3 ,X4 ,Y1 ,
34 3 Y2 ,Y3 ,Y4 ,Z1 ,Z2 ,
35 4 Z3 ,Z4 ,XI ,YI ,ZI ,
36 5 NX1 ,NX2 ,NX3 ,NX4 ,NY1 ,
37 6 NY2 ,NY3 ,NY4 ,NZ1 ,NZ2 ,
38 7 NZ3 ,NZ4 ,LB1 ,LB2 ,LB3 ,
39 8 LB4 ,LC1 ,LC2 ,LC3 ,LC4 ,
40 9 P1 ,P2 ,P3 ,P4 ,IX1 ,
41 A IX2 ,IX3 ,IX4 ,NSVG ,STIF ,
42 B JLT_NEW,GAPV ,INACTI ,CAND_P ,
43 C INDEX ,VXI ,VYI ,
44 D VZI ,MSI ,KINI ,ICURV ,IRECT ,
45 E NNX1 ,NNX2 ,NNX3 ,NNX4 ,NNY1 ,
46 F NNY2 ,NNY3 ,NNY4 ,NNZ1 ,NNZ2 ,
47 G NNZ3 ,NNZ4 ,NOD_NORMAL,IADM,RCURVI ,
48 H ANGLMI ,INTTH ,TEMPI ,PHI ,AREASI ,
49 I IELECI ,NSMS ,CMAJ ,DRAD2 , INTFRIC,IPARTFRICSI,
50 K IPARTFRICMI,ITIED ,JLT_TIED,CAND_F,IORTHFRIC,
51 L IREP_FRICMI,DIR_FRICMI,DGAPLOAD)
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C G l o b a l P a r a m e t e r s
58C-----------------------------------------------
59#include "mvsiz_p.inc"
60#include "sms_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 INTEGER JLT, JLT_NEW, INACTI, INTFRIC, ITIED, JLT_TIED,IORTHFRIC,
65 . CAND_N(*), CN_LOC(MVSIZ), CAND_E(*), CE_LOC(MVSIZ),
66 . NSVG(MVSIZ), KINI(*)
67 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
68 . INDEX(MVSIZ),INTTH,IELECI(MVSIZ),
69 . NSMS(*),IPARTFRICSI(MVSIZ),IPARTFRICMI(MVSIZ),
70 . IREP_FRICMI(*)
71 my_real
72 . NX1(MVSIZ), NX2(MVSIZ), NX3(MVSIZ), NX4(MVSIZ),
73 . NY1(MVSIZ), NY2(MVSIZ), NY3(MVSIZ), NY4(MVSIZ),
74 . NZ1(MVSIZ), NZ2(MVSIZ), NZ3(MVSIZ), NZ4(MVSIZ),
75 . LB1(MVSIZ), LB2(MVSIZ), LB3(MVSIZ), LB4(MVSIZ),
76 . LC1(MVSIZ), LC2(MVSIZ), LC3(MVSIZ), LC4(MVSIZ),
77 . KB1(MVSIZ), KB2(MVSIZ), KB3(MVSIZ), KB4(MVSIZ),
78 . KC1(MVSIZ), KC2(MVSIZ), KC3(MVSIZ), KC4(MVSIZ),
79 . X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
80 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
81 . Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ), Z4(MVSIZ),
82 . XI(MVSIZ), YI(MVSIZ), ZI(MVSIZ), STIF(MVSIZ),
83 . P1(MVSIZ), P2(MVSIZ), P3(MVSIZ), P4(MVSIZ),
84 . gapv(mvsiz), cand_p(*),
85 . vxi(mvsiz), vyi(mvsiz), vzi(mvsiz), msi(mvsiz),
86 . tempi(mvsiz),phi(mvsiz),areasi(mvsiz),dir_fricmi(mvsiz,2)
87 INTEGER IRECT(4,*),ICURV, IADM
88 my_real
89 . NNX1(MVSIZ), NNX2(MVSIZ), NNX3(MVSIZ), NNX4(MVSIZ),
90 . NNY1(MVSIZ), NNY2(MVSIZ), NNY3(MVSIZ), NNY4(MVSIZ),
91 . nnz1(mvsiz), nnz2(mvsiz), nnz3(mvsiz), nnz4(mvsiz),
92 . nod_normal(3,*), cand_f(8,*)
93 my_real
94 . rcurvi(mvsiz), anglmi(mvsiz),cmaj(*)
95 my_real , INTENT(IN) :: drad2, dgapload
96C-----------------------------------------------
97C L o c a l V a r i a b l e s
98C-----------------------------------------------
99 INTEGER I, IG, J, K, KEEP(MVSIZ), IPERM,IORTHFRIC_LOC
100 my_real
101 . X0(MVSIZ), Y0(MVSIZ), Z0(MVSIZ),
102 . AL1(MVSIZ), AL2(MVSIZ), AL3(MVSIZ), AL4(MVSIZ),
103 . X01(MVSIZ), X02(MVSIZ), X03(MVSIZ), X04(MVSIZ),
104 . Y01(MVSIZ), Y02(MVSIZ), Y03(MVSIZ), Y04(MVSIZ),
105 . Z01(MVSIZ), Z02(MVSIZ), Z03(MVSIZ), Z04(MVSIZ),
106 . XI1(MVSIZ), XI2(MVSIZ), XI3(MVSIZ), XI4(MVSIZ),
107 . yi1(mvsiz), yi2(mvsiz), yi3(mvsiz), yi4(mvsiz),
108 . zi1(mvsiz), zi2(mvsiz), zi3(mvsiz), zi4(mvsiz),
109 . pene2(mvsiz),
110 . hlb1(mvsiz), hlc1(mvsiz), hlb2(mvsiz),hlc2(mvsiz),
111 . hlb3(mvsiz), hlc3(mvsiz), hlb4(mvsiz),hlc4(mvsiz)
112 my_real
113 . s2,d1,d2,d3,d4,
114 . x12,x23,x34,x41,xi0,sx1,sx2,sx3,sx4,sx0,
115 . y12,y23,y34,y41,yi0,sy1,sy2,sy3,sy4,sy0,
116 . z12,z23,z34,z41,zi0,sz1,sz2,sz3,sz4,sz0,
117C
118 . gap2, ds2,t1,t2,t3,xxx,yyy,zzz,
119 . al1num,al2num,al3num,al4num,al1den,al2den,al3den,al4den,
120 . x23d,y23d,z23d,x34d,y34d,z34d,x41d,y41d,z41d,
121 . x12d,y12d,z12d,gap2d,xi0d,yi0d,zi0d,s2d, la, hla, aaa,
122 . xi0v(mvsiz), yi0v(mvsiz), zi0v(mvsiz)
123 my_real
124 . rperm
125 LOGICAL ITERATE
126C--------------------------------------------------------
127C CAS DES PAQUETS MIXTES
128C--------------------------------------------------------
129 DO I=1,jlt
130 IF(ix3(i)/=ix4(i))THEN
131 x0(i) = fourth*(x1(i)+x2(i)+x3(i)+x4(i))
132 y0(i) = fourth*(y1(i)+y2(i)+y3(i)+y4(i))
133 z0(i) = fourth*(z1(i)+z2(i)+z3(i)+z4(i))
134 ELSE
135 x0(i) = x3(i)
136 y0(i) = y3(i)
137 z0(i) = z3(i)
138 ENDIF
139 ENDDO
140C
141C--------------------------------------------------------
142C UNIQUEMENT POUR PAQUET DE TRIANGLE
143C--------------------------------------------------------
144C--------------------------------------------------------
145C CAS DES PAQUETS MIXTES OU QUADRANGLE
146C--------------------------------------------------------
147C
148 DO i=1,jlt
149 cmaj(i) = zero
150C
151 x01(i) = x1(i) - x0(i)
152 y01(i) = y1(i) - y0(i)
153 z01(i) = z1(i) - z0(i)
154C
155 x02(i) = x2(i) - x0(i)
156 y02(i) = y2(i) - y0(i)
157 z02(i) = z2(i) - z0(i)
158C
159 x03(i) = x3(i) - x0(i)
160 y03(i) = y3(i) - y0(i)
161 z03(i) = z3(i) - z0(i)
162C
163 x04(i) = x4(i) - x0(i)
164 y04(i) = y4(i) - y0(i)
165 z04(i) = z4(i) - z0(i)
166C
167 xi0v(i) = x0(i) - xi(i)
168 yi0v(i) = y0(i) - yi(i)
169 zi0v(i) = z0(i) - zi(i)
170C
171 xi1(i) = x1(i) - xi(i)
172 yi1(i) = y1(i) - yi(i)
173 zi1(i) = z1(i) - zi(i)
174C
175 xi2(i) = x2(i) - xi(i)
176 yi2(i) = y2(i) - yi(i)
177 zi2(i) = z2(i) - zi(i)
178C
179 xi3(i) = x3(i) - xi(i)
180 yi3(i) = y3(i) - yi(i)
181 zi3(i) = z3(i) - zi(i)
182C
183 xi4(i) = x4(i) - xi(i)
184 yi4(i) = y4(i) - yi(i)
185 zi4(i) = z4(i) - zi(i)
186C
187 sx1 = yi0v(i)*zi1(i) - zi0v(i)*yi1(i)
188 sy1 = zi0v(i)*xi1(i) - xi0v(i)*zi1(i)
189 sz1 = xi0v(i)*yi1(i) - yi0v(i)*xi1(i)
190C
191 sx2 = yi0v(i)*zi2(i) - zi0v(i)*yi2(i)
192 sy2 = zi0v(i)*xi2(i) - xi0v(i)*zi2(i)
193 sz2 = xi0v(i)*yi2(i) - yi0v(i)*xi2(i)
194C
195 sx0 = y01(i)*z02(i) - z01(i)*y02(i)
196 sy0 = z01(i)*x02(i) - x01(i)*z02(i)
197 sz0 = x01(i)*y02(i) - y01(i)*x02(i)
198 s2 = 1./max(em30,sx0*sx0 + sy0*sy0 + sz0*sz0)
199C
200 lb1(i) = -(sx0*sx2 + sy0*sy2 + sz0*sz2) * s2
201 lc1(i) = (sx0*sx1 + sy0*sy1 + sz0*sz1) * s2
202C
203 sx3 = yi0v(i)*zi3(i) - zi0v(i)*yi3(i)
204 sy3 = zi0v(i)*xi3(i) - xi0v(i)*zi3(i)
205 sz3 = xi0v(i)*yi3(i) - yi0v(i)*xi3(i)
206C
207 sx0 = y02(i)*z03(i) - z02(i)*y03(i)
208 sy0 = z02(i)*x03(i) - x02(i)*z03(i)
209 sz0 = x02(i)*y03(i) - y02(i)*x03(i)
210 s2 = 1./max(em30,sx0*sx0 + sy0*sy0 + sz0*sz0)
211C
212 lb2(i) = -(sx0*sx3 + sy0*sy3 + sz0*sz3) * s2
213 lc2(i) = (sx0*sx2 + sy0*sy2 + sz0*sz2) * s2
214C
215 sx4 = yi0v(i)*zi4(i) - zi0v(i)*yi4(i)
216 sy4 = zi0v(i)*xi4(i) - xi0v(i)*zi4(i)
217 sz4 = xi0v(i)*yi4(i) - yi0v(i)*xi4(i)
218C
219 sx0 = y03(i)*z04(i) - z03(i)*y04(i)
220 sy0 = z03(i)*x04(i) - x03(i)*z04(i)
221 sz0 = x03(i)*y04(i) - y03(i)*x04(i)
222 s2 = one/max(em30,sx0*sx0 + sy0*sy0 + sz0*sz0)
223C
224 lb3(i) = -(sx0*sx4 + sy0*sy4 + sz0*sz4) * s2
225 lc3(i) = (sx0*sx3 + sy0*sy3 + sz0*sz3) * s2
226C
227 sx0 = y04(i)*z01(i) - z04(i)*y01(i)
228 sy0 = z04(i)*x01(i) - x04(i)*z01(i)
229 sz0 = x04(i)*y01(i) - y04(i)*x01(i)
230 s2 = one/max(em30,sx0*sx0 + sy0*sy0 + sz0*sz0)
231C
232 lb4(i) = -(sx0*sx1 + sy0*sy1 + sz0*sz1) * s2
233 lc4(i) = (sx0*sx4 + sy0*sy4 + sz0*sz4) * s2
234C
235 ENDDO
236#include "nofusion.inc"
237 DO i=1,jlt
238 aaa = one/max(em30,x01(i)*x01(i)+y01(i)*y01(i)+z01(i)*z01(i))
239 hlc1(i)= lc1(i)*abs(lc1(i))*aaa
240 hlb4(i)= lb4(i)*abs(lb4(i))*aaa
241 al1(i) = -(xi0v(i)*x01(i)+yi0v(i)*y01(i)+zi0v(i)*z01(i))*aaa
242 al1(i) = max(zero,min(one,al1(i)))
243 aaa = one/max(em30,x02(i)*x02(i)+y02(i)*y02(i)+z02(i)*z02(i))
244 hlc2(i)= lc2(i)*abs(lc2(i))*aaa
245 hlb1(i)= lb1(i)*abs(lb1(i))*aaa
246 al2(i) = -(xi0v(i)*x02(i)+yi0v(i)*y02(i)+zi0v(i)*z02(i))*aaa
247 al2(i) = max(zero,min(one,al2(i)))
248 aaa = one/max(em30,x03(i)*x03(i)+y03(i)*y03(i)+z03(i)*z03(i))
249 hlc3(i)= lc3(i)*abs(lc3(i))*aaa
250 hlb2(i)= lb2(i)*abs(lb2(i))*aaa
251 al3(i) = -(xi0v(i)*x03(i)+yi0v(i)*y03(i)+zi0v(i)*z03(i))*aaa
252 al3(i) = max(zero,min(one,al3(i)))
253 aaa = one/max(em30,x04(i)*x04(i)+y04(i)*y04(i)+z04(i)*z04(i))
254 hlc4(i)= lc4(i)*abs(lc4(i))*aaa
255 hlb3(i)= lb3(i)*abs(lb3(i))*aaa
256 al4(i) = -(xi0v(i)*x04(i)+yi0v(i)*y04(i)+zi0v(i)*z04(i))*aaa
257 al4(i) = max(zero,min(one,al4(i)))
258C
259 ENDDO
260C
261#include "nofusion.inc"
262 DO i=1,jlt
263 x12 = x2(i) - x1(i)
264 y12 = y2(i) - y1(i)
265 z12 = z2(i) - z1(i)
266 la = one - lb1(i) - lc1(i)
267C HLA, HLB1, HLC1 necessaires pour triangle angle obtu
268 aaa = one / max(em20,x12*x12+y12*y12+z12*z12)
269 hla= la*abs(la) * aaa
270 IF(la<zero.AND.
271 + hla<=hlb1(i).AND.hla<=hlc1(i))THEN
272 lb1(i) = (xi2(i)*x12+yi2(i)*y12+zi2(i)*z12) * aaa
273 lb1(i) = max(zero,min(one,lb1(i)))
274 lc1(i) = one - lb1(i)
275 ELSEIF(lb1(i)<zero.AND.
276 + hlb1(i)<=hlc1(i).AND.hlb1(i)<=hla)THEN
277 lb1(i) = zero
278 lc1(i) = al2(i)
279 ELSEIF(lc1(i)<zero.AND.
280 + hlc1(i)<=hla.AND.hlc1(i)<=hlb1(i))THEN
281 lc1(i) = zero
282 lb1(i) = al1(i)
283 ENDIF
284 ENDDO
285C
286#include "nofusion.inc"
287 DO i=1,jlt
288 x23 = x3(i) - x2(i)
289 y23 = y3(i) - y2(i)
290 z23 = z3(i) - z2(i)
291 la = one - lb2(i) - lc2(i)
292C HLA, HLB1, HLC1 necessaires pour triangle angle obtu
293 aaa = one / max(em20,x23*x23+y23*y23+z23*z23)
294 hla= la*abs(la) * aaa
295 IF(la<zero.AND.
296 + hla<=hlb2(i).AND.hla<=hlc2(i))THEN
297 lb2(i) = (xi3(i)*x23+yi3(i)*y23+zi3(i)*z23)*aaa
298 lb2(i) = max(zero,min(one,lb2(i)))
299 lc2(i) = one - lb2(i)
300 ELSEIF(lb2(i)<zero.AND.
301 + hlb2(i)<=hlc2(i).AND.hlb2(i)<=hla)THEN
302 lb2(i) = zero
303 lc2(i) = al3(i)
304 ELSEIF(lc2(i)<zero.AND.
305 + hlc2(i)<=hla.AND.hlc2(i)<=hlb2(i))THEN
306 lc2(i) = zero
307 lb2(i) = al2(i)
308 ENDIF
309 ENDDO
310C
311#include "nofusion.inc"
312 DO i=1,jlt
313 x34 = x4(i) - x3(i)
314 y34 = y4(i) - y3(i)
315 z34 = z4(i) - z3(i)
316 la = one - lb3(i) - lc3(i)
317C HLA, HLB1, HLC1 necessaires pour triangle angle obtu
318 aaa = one / max(em20,x34*x34+y34*y34+z34*z34)
319 hla= la*abs(la) * aaa
320 IF(la<zero.AND.
321 + hla<=hlb3(i).AND.hla<=hlc3(i))THEN
322 lb3(i) = (xi4(i)*x34+yi4(i)*y34+zi4(i)*z34)*aaa
323 lb3(i) = max(zero,min(one,lb3(i)))
324 lc3(i) = one - lb3(i)
325 ELSEIF(lb3(i)<zero.AND.
326 + hlb3(i)<=hlc3(i).AND.hlb3(i)<=hla)THEN
327 lb3(i) = zero
328 lc3(i) = al4(i)
329 ELSEIF(lc3(i)<zero.AND.
330 + hlc3(i)<=hla.AND.hlc3(i)<=hlb3(i))THEN
331 lc3(i) = zero
332 lb3(i) = al3(i)
333 ENDIF
334 ENDDO
335C
336#include "nofusion.inc"
337 DO i=1,jlt
338 x41 = x1(i) - x4(i)
339 y41 = y1(i) - y4(i)
340 z41 = z1(i) - z4(i)
341 la = one - lb4(i) - lc4(i)
342C HLA, HLB1, HLC1 necessaires pour triangle angle obtu
343 aaa = one / max(em20,x41*x41+y41*y41+z41*z41)
344 hla= la*abs(la) * aaa
345 IF(la<zero.AND.
346 + hla<=hlb4(i).AND.hla<=hlc4(i))THEN
347 lb4(i) = (xi1(i)*x41+yi1(i)*y41+zi1(i)*z41)*aaa
348 lb4(i) = max(zero,min(one,lb4(i)))
349 lc4(i) = one - lb4(i)
350 ELSEIF(lb4(i)<zero.AND.
351 + hlb4(i)<=hlc4(i).AND.hlb4(i)<=hla)THEN
352 lb4(i) = zero
353 lc4(i) = al1(i)
354 ELSEIF(lc4(i)<zero.AND.
355 + hlc4(i)<=hla.AND.hlc4(i)<=hlb4(i))THEN
356 lc4(i) = zero
357 lb4(i) = al4(i)
358 ENDIF
359 ENDDO
360c---------------------------------------------------------
361c courbure cubique
362c---------------------------------------------------------
363 IF(icurv == 3)THEN
364 CALL i7cmaj(jlt ,cmaj ,irect ,nod_normal,cand_e,
365 2 x1 ,x2 ,x3 ,x4 ,
366 3 y1 ,y2 ,y3 ,y4 ,
367 4 z1 ,z2 ,z3 ,z4 ,
368 5 nnx1 ,nnx2 ,nnx3 ,nnx4 ,
369 6 nny1 ,nny2 ,nny3 ,nny4 ,
370 7 nnz1 ,nnz2 ,nnz3 ,nnz4 )
371 ELSEIF(icurv /= 0)THEN
372#include "nofusion.inc"
373 DO i=1,jlt
374 xxx=max(x1(i),x2(i),x3(i),x4(i))
375 . -min(x1(i),x2(i),x3(i),x4(i))
376 yyy=max(y1(i),y2(i),y3(i),y4(i))
377 . -min(y1(i),y2(i),y3(i),y4(i))
378 zzz=max(z1(i),z2(i),z3(i),z4(i))
379 . -min(z1(i),z2(i),z3(i),z4(i))
380 cmaj(i) = half * max(xxx,yyy,zzz)
381 ENDDO
382 ENDIF
383c
384#include "nofusion.inc"
385 DO i=1,jlt
386 gap2=(gapv(i)+cmaj(i)+dgapload)*(gapv(i)+cmaj(i)+dgapload)
387 gap2 =max(gap2,drad2)
388C
389 nx1(i) = xi(i)-(x0(i) + lb1(i)*x01(i) + lc1(i)*x02(i))
390 ny1(i) = yi(i)-(y0(i) + lb1(i)*y01(i) + lc1(i)*y02(i))
391 nz1(i) = zi(i)-(z0(i) + lb1(i)*z01(i) + lc1(i)*z02(i))
392 p1(i) = nx1(i)*nx1(i) + ny1(i)*ny1(i) +nz1(i)*nz1(i)
393 d1 = max(zero, gap2 - p1(i))
394C
395 nx2(i) = xi(i)-(x0(i) + lb2(i)*x02(i) + lc2(i)*x03(i))
396 ny2(i) = yi(i)-(y0(i) + lb2(i)*y02(i) + lc2(i)*y03(i))
397 nz2(i) = zi(i)-(z0(i) + lb2(i)*z02(i) + lc2(i)*z03(i))
398 p2(i) = nx2(i)*nx2(i) + ny2(i)*ny2(i) +nz2(i)*nz2(i)
399 d2 = max(zero, gap2 - p2(i))
400C
401 nx3(i) = xi(i)-(x0(i) + lb3(i)*x03(i) + lc3(i)*x04(i))
402 ny3(i) = yi(i)-(y0(i) + lb3(i)*y03(i) + lc3(i)*y04(i))
403 nz3(i) = zi(i)-(z0(i) + lb3(i)*z03(i) + lc3(i)*z04(i))
404 p3(i) = nx3(i)*nx3(i) + ny3(i)*ny3(i) +nz3(i)*nz3(i)
405 d3 = max(zero, gap2 - p3(i))
406C
407 nx4(i) = xi(i)-(x0(i) + lb4(i)*x04(i) + lc4(i)*x01(i))
408 ny4(i) = yi(i)-(y0(i) + lb4(i)*y04(i) + lc4(i)*y01(i))
409 nz4(i) = zi(i)-(z0(i) + lb4(i)*z04(i) + lc4(i)*z01(i))
410 p4(i) = nx4(i)*nx4(i) + ny4(i)*ny4(i) +nz4(i)*nz4(i)
411 d4 = max(zero, gap2 - p4(i))
412c IF(ITAB(NSVG(I))==40142.AND.ITAB(IX1(I))==40016.AND.ITAB(IX2(I))==40017.AND.ITAB(IX3(I))==40080
413c . .AND.ITAB(IX4(I))==40079) THEN
414c write(*,*) ' P1(I) P2(I) P3(I) P4(I) ',P1(I) , P2(I) , P3(I) , P4(I)
415c PAUSE
416c ENDIF
417C !!!!!!!!!!!!!!!!!!!!!!!
418C PENE2 = GAP^2 - DIST^2 UTILISE POUR TESTER SI NON NUL
419C!!!!!!!!!!!!!!!!!!!!!!!!
420 pene2(i) = max(d1,d2,d3,d4)
421C
422 ENDDO
423C---------------------
424C PENE INITIALE
425C---------------------
426 keep(1:jlt)=0
427 IF(itied == 0)THEN
428 IF(inacti==5.OR.inacti==6)THEN
429 DO i=1,jlt
430 IF(pene2(i)==zero.OR.stif(i)==zero)THEN
431 cand_p(index(i))=zero
432 ENDIF
433 ENDDO
434 ENDIF
435#include "vectorize.inc"
436 DO i=1,jlt
437 IF( pene2(i)/=zero .AND. stif(i)/=zero)THEN
438 keep(i) =1
439 END IF
440 ENDDO
441 ELSE
442 IF(inacti==5.OR.inacti==6)THEN
443 DO i=1,jlt
444 IF((pene2(i)==zero.AND.cand_f(1,index(i))==zero).OR.
445 . stif(i)==zero)THEN
446 cand_p(index(i))=zero
447 ENDIF
448 ENDDO
449 ENDIF
450 DO i=1,jlt
451 IF( (pene2(i)/=zero.OR.cand_f(1,index(i))/=zero) .AND.
452 . stif(i)/=zero)THEN
453 IF( cand_f(1,index(i))/=zero )THEN
454 jlt_tied =jlt_tied+1
455 keep(i) =2
456 ELSE
457 keep(i) =1
458 END IF
459 END IF
460 ENDDO
461 END IF
462C
463
464c----copy IORTHFRIC in local flag
465
466 iorthfric_loc = iorthfric
467 IF(idtmins < 2)THEN
468C
469 IF(intfric == 0) THEN
470 IF(intth == 0 ) THEN
471 IF(icurv==0.AND.iadm==0 )THEN
472 DO i=1,jlt
473 IF(keep(i)/=0)THEN
474 jlt_new = jlt_new + 1
475 cn_loc(jlt_new) = cand_n(i)
476 ce_loc(jlt_new) = cand_e(i)
477 ix1(jlt_new) = ix1(i)
478 ix2(jlt_new) = ix2(i)
479 ix3(jlt_new) = ix3(i)
480 ix4(jlt_new) = ix4(i)
481 nsvg(jlt_new) = nsvg(i)
482 nx1(jlt_new) = nx1(i)
483 nx2(jlt_new) = nx2(i)
484 nx3(jlt_new) = nx3(i)
485 nx4(jlt_new) = nx4(i)
486 ny1(jlt_new) = ny1(i)
487 ny2(jlt_new) = ny2(i)
488 ny3(jlt_new) = ny3(i)
489 ny4(jlt_new) = ny4(i)
490 nz1(jlt_new) = nz1(i)
491 nz2(jlt_new) = nz2(i)
492 nz3(jlt_new) = nz3(i)
493 nz4(jlt_new) = nz4(i)
494 p1(jlt_new) = p1(i)
495 p2(jlt_new) = p2(i)
496 p3(jlt_new) = p3(i)
497 p4(jlt_new) = p4(i)
498 lb1(jlt_new) = lb1(i)
499 lb2(jlt_new) = lb2(i)
500 lb3(jlt_new) = lb3(i)
501 lb4(jlt_new) = lb4(i)
502 lc1(jlt_new) = lc1(i)
503 lc2(jlt_new) = lc2(i)
504 lc3(jlt_new) = lc3(i)
505 lc4(jlt_new) = lc4(i)
506 stif(jlt_new) = stif(i)
507 gapv(jlt_new) = gapv(i)
508 index(jlt_new)= index(i)
509C
510 kini(jlt_new) = kini(i)
511 vxi(jlt_new) = vxi(i)
512 vyi(jlt_new) = vyi(i)
513 vzi(jlt_new) = vzi(i)
514 msi(jlt_new) = msi(i)
515C
516 xi(jlt_new) = xi(i)
517 yi(jlt_new) = yi(i)
518 zi(jlt_new) = zi(i)
519 x1(jlt_new) = x1(i)
520 y1(jlt_new) = y1(i)
521 z1(jlt_new) = z1(i)
522 x2(jlt_new) = x2(i)
523 y2(jlt_new) = y2(i)
524 z2(jlt_new) = z2(i)
525 x3(jlt_new) = x3(i)
526 y3(jlt_new) = y3(i)
527 z3(jlt_new) = z3(i)
528 x4(jlt_new) = x4(i)
529 y4(jlt_new) = y4(i)
530 z4(jlt_new) = z4(i)
531C
532C
533 keep(jlt_new) = keep(i)
534C
535 ENDIF
536 ENDDO
537 ELSE !
538 DO i=1,jlt
539 IF(keep(i)/=0)THEN
540 jlt_new = jlt_new + 1
541 cn_loc(jlt_new) = cand_n(i)
542 ce_loc(jlt_new) = cand_e(i)
543 ix1(jlt_new) = ix1(i)
544 ix2(jlt_new) = ix2(i)
545 ix3(jlt_new) = ix3(i)
546 ix4(jlt_new) = ix4(i)
547 nsvg(jlt_new) = nsvg(i)
548 nx1(jlt_new) = nx1(i)
549 nx2(jlt_new) = nx2(i)
550 nx3(jlt_new) = nx3(i)
551 nx4(jlt_new) = nx4(i)
552 ny1(jlt_new) = ny1(i)
553 ny2(jlt_new) = ny2(i)
554 ny3(jlt_new) = ny3(i)
555 ny4(jlt_new) = ny4(i)
556 nz1(jlt_new) = nz1(i)
557 nz2(jlt_new) = nz2(i)
558 nz3(jlt_new) = nz3(i)
559 nz4(jlt_new) = nz4(i)
560 p1(jlt_new) = p1(i)
561 p2(jlt_new) = p2(i)
562 p3(jlt_new) = p3(i)
563 p4(jlt_new) = p4(i)
564 lb1(jlt_new) = lb1(i)
565 lb2(jlt_new) = lb2(i)
566 lb3(jlt_new) = lb3(i)
567 lb4(jlt_new) = lb4(i)
568 lc1(jlt_new) = lc1(i)
569 lc2(jlt_new) = lc2(i)
570 lc3(jlt_new) = lc3(i)
571 lc4(jlt_new) = lc4(i)
572 stif(jlt_new) = stif(i)
573 gapv(jlt_new) = gapv(i)
574 index(jlt_new)= index(i)
575 kini(jlt_new) = kini(i)
576 vxi(jlt_new) = vxi(i)
577 vyi(jlt_new) = vyi(i)
578 vzi(jlt_new) = vzi(i)
579 msi(jlt_new) = msi(i)
580C
581 xi(jlt_new) = xi(i)
582 yi(jlt_new) = yi(i)
583 zi(jlt_new) = zi(i)
584 x1(jlt_new) = x1(i)
585 y1(jlt_new) = y1(i)
586 z1(jlt_new) = z1(i)
587 x2(jlt_new) = x2(i)
588 y2(jlt_new) = y2(i)
589 z2(jlt_new) = z2(i)
590 x3(jlt_new) = x3(i)
591 y3(jlt_new) = y3(i)
592 z3(jlt_new) = z3(i)
593 x4(jlt_new) = x4(i)
594 y4(jlt_new) = y4(i)
595 z4(jlt_new) = z4(i)
596 rcurvi(jlt_new) = rcurvi(i)
597 anglmi(jlt_new) = anglmi(i)
598 keep(jlt_new) = keep(i)
599C
600 ENDIF
601 ENDDO
602 ENDIF
603 ELSE
604 IF(iadm == 0 )THEN
605 DO i=1,jlt
606 IF(keep(i)/=0)THEN
607 jlt_new = jlt_new + 1
608 cn_loc(jlt_new) = cand_n(i)
609 ce_loc(jlt_new) = cand_e(i)
610 ix1(jlt_new) = ix1(i)
611 ix2(jlt_new) = ix2(i)
612 ix3(jlt_new) = ix3(i)
613 ix4(jlt_new) = ix4(i)
614 nsvg(jlt_new) = nsvg(i)
615 nx1(jlt_new) = nx1(i)
616 nx2(jlt_new) = nx2(i)
617 nx3(jlt_new) = nx3(i)
618 nx4(jlt_new) = nx4(i)
619 ny1(jlt_new) = ny1(i)
620 ny2(jlt_new) = ny2(i)
621 ny3(jlt_new) = ny3(i)
622 ny4(jlt_new) = ny4(i)
623 nz1(jlt_new) = nz1(i)
624 nz2(jlt_new) = nz2(i)
625 nz3(jlt_new) = nz3(i)
626 nz4(jlt_new) = nz4(i)
627 p1(jlt_new) = p1(i)
628 p2(jlt_new) = p2(i)
629 p3(jlt_new) = p3(i)
630 p4(jlt_new) = p4(i)
631 lb1(jlt_new) = lb1(i)
632 lb2(jlt_new) = lb2(i)
633 lb3(jlt_new) = lb3(i)
634 lb4(jlt_new) = lb4(i)
635 lc1(jlt_new) = lc1(i)
636 lc2(jlt_new) = lc2(i)
637 lc3(jlt_new) = lc3(i)
638 lc4(jlt_new) = lc4(i)
639 stif(jlt_new) = stif(i)
640 gapv(jlt_new) = gapv(i)
641 index(jlt_new)= index(i)
642C
643 kini(jlt_new) = kini(i)
644 vxi(jlt_new) = vxi(i)
645 vyi(jlt_new) = vyi(i)
646 vzi(jlt_new) = vzi(i)
647 msi(jlt_new) = msi(i)
648C
649 xi(jlt_new) = xi(i)
650 yi(jlt_new) = yi(i)
651 zi(jlt_new) = zi(i)
652 x1(jlt_new) = x1(i)
653 y1(jlt_new) = y1(i)
654 z1(jlt_new) = z1(i)
655 x2(jlt_new) = x2(i)
656 y2(jlt_new) = y2(i)
657 z2(jlt_new) = z2(i)
658 x3(jlt_new) = x3(i)
659 y3(jlt_new) = y3(i)
660 z3(jlt_new) = z3(i)
661 x4(jlt_new) = x4(i)
662 y4(jlt_new) = y4(i)
663 z4(jlt_new) = z4(i)
664C
665 tempi(jlt_new) = tempi(i)
666 phi(jlt_new) = zero
667 areasi(jlt_new) = areasi(i)
668 ieleci(jlt_new) =ieleci(i)
669C
670C
671 keep(jlt_new) = keep(i)
672C
673 ENDIF
674 ENDDO
675 ELSE !
676 DO i=1,jlt
677 IF(keep(i)/=0)THEN
678 jlt_new = jlt_new + 1
679 cn_loc(jlt_new) = cand_n(i)
680 ce_loc(jlt_new) = cand_e(i)
681 ix1(jlt_new) = ix1(i)
682 ix2(jlt_new) = ix2(i)
683 ix3(jlt_new) = ix3(i)
684 ix4(jlt_new) = ix4(i)
685 nsvg(jlt_new) = nsvg(i)
686 nx1(jlt_new) = nx1(i)
687 nx2(jlt_new) = nx2(i)
688 nx3(jlt_new) = nx3(i)
689 nx4(jlt_new) = nx4(i)
690 ny1(jlt_new) = ny1(i)
691 ny2(jlt_new) = ny2(i)
692 ny3(jlt_new) = ny3(i)
693 ny4(jlt_new) = ny4(i)
694 nz1(jlt_new) = nz1(i)
695 nz2(jlt_new) = nz2(i)
696 nz3(jlt_new) = nz3(i)
697 nz4(jlt_new) = nz4(i)
698 p1(jlt_new) = p1(i)
699 p2(jlt_new) = p2(i)
700 p3(jlt_new) = p3(i)
701 p4(jlt_new) = p4(i)
702 lb1(jlt_new) = lb1(i)
703 lb2(jlt_new) = lb2(i)
704 lb3(jlt_new) = lb3(i)
705 lb4(jlt_new) = lb4(i)
706 lc1(jlt_new) = lc1(i)
707 lc2(jlt_new) = lc2(i)
708 lc3(jlt_new) = lc3(i)
709 lc4(jlt_new) = lc4(i)
710 stif(jlt_new) = stif(i)
711 gapv(jlt_new) = gapv(i)
712 index(jlt_new)= index(i)
713 kini(jlt_new) = kini(i)
714 vxi(jlt_new) = vxi(i)
715 vyi(jlt_new) = vyi(i)
716 vzi(jlt_new) = vzi(i)
717 msi(jlt_new) = msi(i)
718 xi(jlt_new) = xi(i)
719 yi(jlt_new) = yi(i)
720 zi(jlt_new) = zi(i)
721 x1(jlt_new) = x1(i)
722 y1(jlt_new) = y1(i)
723 z1(jlt_new) = z1(i)
724 x2(jlt_new) = x2(i)
725 y2(jlt_new) = y2(i)
726 z2(jlt_new) = z2(i)
727 x3(jlt_new) = x3(i)
728 y3(jlt_new) = y3(i)
729 z3(jlt_new) = z3(i)
730 x4(jlt_new) = x4(i)
731 y4(jlt_new) = y4(i)
732 z4(jlt_new) = z4(i)
733 rcurvi(jlt_new) = rcurvi(i)
734 anglmi(jlt_new) = anglmi(i)
735 tempi(jlt_new) = tempi(i)
736 phi(jlt_new) = zero
737 areasi(jlt_new) = areasi(i)
738 ieleci(jlt_new) =ieleci(i)
739C
740C
741 keep(jlt_new) = keep(i)
742C
743 ENDIF
744 ENDDO
745 ENDIF
746 ENDIF
747 ELSE ! INTFRIC
748
749 IF(intth == 0 ) THEN
750 IF(icurv==0.AND.iadm==0 )THEN
751 DO i=1,jlt
752 IF(keep(i)/=0)THEN
753 jlt_new = jlt_new + 1
754 cn_loc(jlt_new) = cand_n(i)
755 ce_loc(jlt_new) = cand_e(i)
756 ix1(jlt_new) = ix1(i)
757 ix2(jlt_new) = ix2(i)
758 ix3(jlt_new) = ix3(i)
759 ix4(jlt_new) = ix4(i)
760 nsvg(jlt_new) = nsvg(i)
761 nx1(jlt_new) = nx1(i)
762 nx2(jlt_new) = nx2(i)
763 nx3(jlt_new) = nx3(i)
764 nx4(jlt_new) = nx4(i)
765 ny1(jlt_new) = ny1(i)
766 ny2(jlt_new) = ny2(i)
767 ny3(jlt_new) = ny3(i)
768 ny4(jlt_new) = ny4(i)
769 nz1(jlt_new) = nz1(i)
770 nz2(jlt_new) = nz2(i)
771 nz3(jlt_new) = nz3(i)
772 nz4(jlt_new) = nz4(i)
773 p1(jlt_new) = p1(i)
774 p2(jlt_new) = p2(i)
775 p3(jlt_new) = p3(i)
776 p4(jlt_new) = p4(i)
777 lb1(jlt_new) = lb1(i)
778 lb2(jlt_new) = lb2(i)
779 lb3(jlt_new) = lb3(i)
780 lb4(jlt_new) = lb4(i)
781 lc1(jlt_new) = lc1(i)
782 lc2(jlt_new) = lc2(i)
783 lc3(jlt_new) = lc3(i)
784 lc4(jlt_new) = lc4(i)
785 stif(jlt_new) = stif(i)
786 gapv(jlt_new) = gapv(i)
787 index(jlt_new)= index(i)
788C
789 kini(jlt_new) = kini(i)
790 vxi(jlt_new) = vxi(i)
791 vyi(jlt_new) = vyi(i)
792 vzi(jlt_new) = vzi(i)
793 msi(jlt_new) = msi(i)
794C
795 xi(jlt_new) = xi(i)
796 yi(jlt_new) = yi(i)
797 zi(jlt_new) = zi(i)
798 x1(jlt_new) = x1(i)
799 y1(jlt_new) = y1(i)
800 z1(jlt_new) = z1(i)
801 x2(jlt_new) = x2(i)
802 y2(jlt_new) = y2(i)
803 z2(jlt_new) = z2(i)
804 x3(jlt_new) = x3(i)
805 y3(jlt_new) = y3(i)
806 z3(jlt_new) = z3(i)
807 x4(jlt_new) = x4(i)
808 y4(jlt_new) = y4(i)
809 z4(jlt_new) = z4(i)
810C
811 ipartfricsi(jlt_new)=ipartfricsi(i)
812 ipartfricmi(jlt_new)=ipartfricmi(i)
813 IF(iorthfric_loc >0) THEN
814 irep_fricmi(jlt_new)= irep_fricmi(i)
815 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
816 ENDIF
817C
818 keep(jlt_new) = keep(i)
819C
820 ENDIF
821 ENDDO
822 ELSE !
823 DO i=1,jlt
824 IF(keep(i)/=0)THEN
825 jlt_new = jlt_new + 1
826 cn_loc(jlt_new) = cand_n(i)
827 ce_loc(jlt_new) = cand_e(i)
828 ix1(jlt_new) = ix1(i)
829 ix2(jlt_new) = ix2(i)
830 ix3(jlt_new) = ix3(i)
831 ix4(jlt_new) = ix4(i)
832 nsvg(jlt_new) = nsvg(i)
833 nx1(jlt_new) = nx1(i)
834 nx2(jlt_new) = nx2(i)
835 nx3(jlt_new) = nx3(i)
836 nx4(jlt_new) = nx4(i)
837 ny1(jlt_new) = ny1(i)
838 ny2(jlt_new) = ny2(i)
839 ny3(jlt_new) = ny3(i)
840 ny4(jlt_new) = ny4(i)
841 nz1(jlt_new) = nz1(i)
842 nz2(jlt_new) = nz2(i)
843 nz3(jlt_new) = nz3(i)
844 nz4(jlt_new) = nz4(i)
845 p1(jlt_new) = p1(i)
846 p2(jlt_new) = p2(i)
847 p3(jlt_new) = p3(i)
848 p4(jlt_new) = p4(i)
849 lb1(jlt_new) = lb1(i)
850 lb2(jlt_new) = lb2(i)
851 lb3(jlt_new) = lb3(i)
852 lb4(jlt_new) = lb4(i)
853 lc1(jlt_new) = lc1(i)
854 lc2(jlt_new) = lc2(i)
855 lc3(jlt_new) = lc3(i)
856 lc4(jlt_new) = lc4(i)
857 stif(jlt_new) = stif(i)
858 gapv(jlt_new) = gapv(i)
859 index(jlt_new)= index(i)
860 kini(jlt_new) = kini(i)
861 vxi(jlt_new) = vxi(i)
862 vyi(jlt_new) = vyi(i)
863 vzi(jlt_new) = vzi(i)
864 msi(jlt_new) = msi(i)
865C
866 xi(jlt_new) = xi(i)
867 yi(jlt_new) = yi(i)
868 zi(jlt_new) = zi(i)
869 x1(jlt_new) = x1(i)
870 y1(jlt_new) = y1(i)
871 z1(jlt_new) = z1(i)
872 x2(jlt_new) = x2(i)
873 y2(jlt_new) = y2(i)
874 z2(jlt_new) = z2(i)
875 x3(jlt_new) = x3(i)
876 y3(jlt_new) = y3(i)
877 z3(jlt_new) = z3(i)
878 x4(jlt_new) = x4(i)
879 y4(jlt_new) = y4(i)
880 z4(jlt_new) = z4(i)
881 rcurvi(jlt_new) = rcurvi(i)
882 anglmi(jlt_new) = anglmi(i)
883 ipartfricsi(jlt_new)=ipartfricsi(i)
884 ipartfricmi(jlt_new)=ipartfricmi(i)
885 IF(iorthfric_loc >0) THEN
886 irep_fricmi(jlt_new)= irep_fricmi(i)
887 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
888 ENDIF
889C
890 keep(jlt_new) = keep(i)
891C
892 ENDIF
893 ENDDO
894 ENDIF
895 ELSE
896 IF(iadm == 0 )THEN
897 DO i=1,jlt
898 IF(keep(i)/=0)THEN
899 jlt_new = jlt_new + 1
900 cn_loc(jlt_new) = cand_n(i)
901 ce_loc(jlt_new) = cand_e(i)
902 ix1(jlt_new) = ix1(i)
903 ix2(jlt_new) = ix2(i)
904 ix3(jlt_new) = ix3(i)
905 ix4(jlt_new) = ix4(i)
906 nsvg(jlt_new) = nsvg(i)
907 nx1(jlt_new) = nx1(i)
908 nx2(jlt_new) = nx2(i)
909 nx3(jlt_new) = nx3(i)
910 nx4(jlt_new) = nx4(i)
911 ny1(jlt_new) = ny1(i)
912 ny2(jlt_new) = ny2(i)
913 ny3(jlt_new) = ny3(i)
914 ny4(jlt_new) = ny4(i)
915 nz1(jlt_new) = nz1(i)
916 nz2(jlt_new) = nz2(i)
917 nz3(jlt_new) = nz3(i)
918 nz4(jlt_new) = nz4(i)
919 p1(jlt_new) = p1(i)
920 p2(jlt_new) = p2(i)
921 p3(jlt_new) = p3(i)
922 p4(jlt_new) = p4(i)
923 lb1(jlt_new) = lb1(i)
924 lb2(jlt_new) = lb2(i)
925 lb3(jlt_new) = lb3(i)
926 lb4(jlt_new) = lb4(i)
927 lc1(jlt_new) = lc1(i)
928 lc2(jlt_new) = lc2(i)
929 lc3(jlt_new) = lc3(i)
930 lc4(jlt_new) = lc4(i)
931 stif(jlt_new) = stif(i)
932 gapv(jlt_new) = gapv(i)
933 index(jlt_new)= index(i)
934C
935 kini(jlt_new) = kini(i)
936 vxi(jlt_new) = vxi(i)
937 vyi(jlt_new) = vyi(i)
938 vzi(jlt_new) = vzi(i)
939 msi(jlt_new) = msi(i)
940C
941 xi(jlt_new) = xi(i)
942 yi(jlt_new) = yi(i)
943 zi(jlt_new) = zi(i)
944 x1(jlt_new) = x1(i)
945 y1(jlt_new) = y1(i)
946 z1(jlt_new) = z1(i)
947 x2(jlt_new) = x2(i)
948 y2(jlt_new) = y2(i)
949 z2(jlt_new) = z2(i)
950 x3(jlt_new) = x3(i)
951 y3(jlt_new) = y3(i)
952 z3(jlt_new) = z3(i)
953 x4(jlt_new) = x4(i)
954 y4(jlt_new) = y4(i)
955 z4(jlt_new) = z4(i)
956C
957 tempi(jlt_new) = tempi(i)
958 phi(jlt_new) = zero
959 areasi(jlt_new) = areasi(i)
960 ieleci(jlt_new) =ieleci(i)
961C
962 ipartfricsi(jlt_new)=ipartfricsi(i)
963 ipartfricmi(jlt_new)=ipartfricmi(i)
964 IF(iorthfric_loc >0) THEN
965 irep_fricmi(jlt_new)= irep_fricmi(i)
966 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
967 ENDIF
968C
969 keep(jlt_new) = keep(i)
970C
971 ENDIF
972 ENDDO
973 ELSE !
974 DO i=1,jlt
975 IF(keep(i)/=0)THEN
976 jlt_new = jlt_new + 1
977 cn_loc(jlt_new) = cand_n(i)
978 ce_loc(jlt_new) = cand_e(i)
979 ix1(jlt_new) = ix1(i)
980 ix2(jlt_new) = ix2(i)
981 ix3(jlt_new) = ix3(i)
982 ix4(jlt_new) = ix4(i)
983 nsvg(jlt_new) = nsvg(i)
984 nx1(jlt_new) = nx1(i)
985 nx2(jlt_new) = nx2(i)
986 nx3(jlt_new) = nx3(i)
987 nx4(jlt_new) = nx4(i)
988 ny1(jlt_new) = ny1(i)
989 ny2(jlt_new) = ny2(i)
990 ny3(jlt_new) = ny3(i)
991 ny4(jlt_new) = ny4(i)
992 nz1(jlt_new) = nz1(i)
993 nz2(jlt_new) = nz2(i)
994 nz3(jlt_new) = nz3(i)
995 nz4(jlt_new) = nz4(i)
996 p1(jlt_new) = p1(i)
997 p2(jlt_new) = p2(i)
998 p3(jlt_new) = p3(i)
999 p4(jlt_new) = p4(i)
1000 lb1(jlt_new) = lb1(i)
1001 lb2(jlt_new) = lb2(i)
1002 lb3(jlt_new) = lb3(i)
1003 lb4(jlt_new) = lb4(i)
1004 lc1(jlt_new) = lc1(i)
1005 lc2(jlt_new) = lc2(i)
1006 lc3(jlt_new) = lc3(i)
1007 lc4(jlt_new) = lc4(i)
1008 stif(jlt_new) = stif(i)
1009 gapv(jlt_new) = gapv(i)
1010 index(jlt_new)= index(i)
1011 kini(jlt_new) = kini(i)
1012 vxi(jlt_new) = vxi(i)
1013 vyi(jlt_new) = vyi(i)
1014 vzi(jlt_new) = vzi(i)
1015 msi(jlt_new) = msi(i)
1016 xi(jlt_new) = xi(i)
1017 yi(jlt_new) = yi(i)
1018 zi(jlt_new) = zi(i)
1019 x1(jlt_new) = x1(i)
1020 y1(jlt_new) = y1(i)
1021 z1(jlt_new) = z1(i)
1022 x2(jlt_new) = x2(i)
1023 y2(jlt_new) = y2(i)
1024 z2(jlt_new) = z2(i)
1025 x3(jlt_new) = x3(i)
1026 y3(jlt_new) = y3(i)
1027 z3(jlt_new) = z3(i)
1028 x4(jlt_new) = x4(i)
1029 y4(jlt_new) = y4(i)
1030 z4(jlt_new) = z4(i)
1031 rcurvi(jlt_new) = rcurvi(i)
1032 anglmi(jlt_new) = anglmi(i)
1033 tempi(jlt_new) = tempi(i)
1034 phi(jlt_new) = zero
1035 areasi(jlt_new) = areasi(i)
1036 ieleci(jlt_new) =ieleci(i)
1037 ipartfricsi(jlt_new)=ipartfricsi(i)
1038 ipartfricmi(jlt_new)=ipartfricmi(i)
1039 IF(iorthfric_loc >0) THEN
1040 irep_fricmi(jlt_new)= irep_fricmi(i)
1041 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
1042 ENDIF
1043C
1044 keep(jlt_new) = keep(i)
1045C
1046 ENDIF
1047 ENDDO
1048 ENDIF
1049 ENDIF
1050 ENDIF
1051
1052
1053C
1054C IDTMINS=2
1055C ---------
1056 ELSE
1057C
1058 IF(intfric == 0) THEN
1059 IF(intth == 0 ) THEN
1060 IF(icurv==0.AND.iadm==0 )THEN
1061 DO i=1,jlt
1062 IF(keep(i)/=0)THEN
1063 jlt_new = jlt_new + 1
1064 cn_loc(jlt_new) = cand_n(i)
1065 ce_loc(jlt_new) = cand_e(i)
1066 ix1(jlt_new) = ix1(i)
1067 ix2(jlt_new) = ix2(i)
1068 ix3(jlt_new) = ix3(i)
1069 ix4(jlt_new) = ix4(i)
1070 nsvg(jlt_new) = nsvg(i)
1071 nx1(jlt_new) = nx1(i)
1072 nx2(jlt_new) = nx2(i)
1073 nx3(jlt_new) = nx3(i)
1074 nx4(jlt_new) = nx4(i)
1075 ny1(jlt_new) = ny1(i)
1076 ny2(jlt_new) = ny2(i)
1077 ny3(jlt_new) = ny3(i)
1078 ny4(jlt_new) = ny4(i)
1079 nz1(jlt_new) = nz1(i)
1080 nz2(jlt_new) = nz2(i)
1081 nz3(jlt_new) = nz3(i)
1082 nz4(jlt_new) = nz4(i)
1083 p1(jlt_new) = p1(i)
1084 p2(jlt_new) = p2(i)
1085 p3(jlt_new) = p3(i)
1086 p4(jlt_new) = p4(i)
1087 lb1(jlt_new) = lb1(i)
1088 lb2(jlt_new) = lb2(i)
1089 lb3(jlt_new) = lb3(i)
1090 lb4(jlt_new) = lb4(i)
1091 lc1(jlt_new) = lc1(i)
1092 lc2(jlt_new) = lc2(i)
1093 lc3(jlt_new) = lc3(i)
1094 lc4(jlt_new) = lc4(i)
1095 stif(jlt_new) = stif(i)
1096 gapv(jlt_new) = gapv(i)
1097 index(jlt_new)= index(i)
1098C
1099 kini(jlt_new) = kini(i)
1100 vxi(jlt_new) = vxi(i)
1101 vyi(jlt_new) = vyi(i)
1102 vzi(jlt_new) = vzi(i)
1103 msi(jlt_new) = msi(i)
1104C
1105 xi(jlt_new) = xi(i)
1106 yi(jlt_new) = yi(i)
1107 zi(jlt_new) = zi(i)
1108 x1(jlt_new) = x1(i)
1109 y1(jlt_new) = y1(i)
1110 z1(jlt_new) = z1(i)
1111 x2(jlt_new) = x2(i)
1112 y2(jlt_new) = y2(i)
1113 z2(jlt_new) = z2(i)
1114 x3(jlt_new) = x3(i)
1115 y3(jlt_new) = y3(i)
1116 z3(jlt_new) = z3(i)
1117 x4(jlt_new) = x4(i)
1118 y4(jlt_new) = y4(i)
1119 z4(jlt_new) = z4(i)
1120C
1121 nsms(jlt_new) = nsms(i)
1122C
1123 keep(jlt_new) = keep(i)
1124C
1125 ENDIF
1126 ENDDO
1127 ELSE !
1128 DO i=1,jlt
1129 IF(keep(i)/=0)THEN
1130 jlt_new = jlt_new + 1
1131 cn_loc(jlt_new) = cand_n(i)
1132 ce_loc(jlt_new) = cand_e(i)
1133 ix1(jlt_new) = ix1(i)
1134 ix2(jlt_new) = ix2(i)
1135 ix3(jlt_new) = ix3(i)
1136 ix4(jlt_new) = ix4(i)
1137 nsvg(jlt_new) = nsvg(i)
1138 nx1(jlt_new) = nx1(i)
1139 nx2(jlt_new) = nx2(i)
1140 nx3(jlt_new) = nx3(i)
1141 nx4(jlt_new) = nx4(i)
1142 ny1(jlt_new) = ny1(i)
1143 ny2(jlt_new) = ny2(i)
1144 ny3(jlt_new) = ny3(i)
1145 ny4(jlt_new) = ny4(i)
1146 nz1(jlt_new) = nz1(i)
1147 nz2(jlt_new) = nz2(i)
1148 nz3(jlt_new) = nz3(i)
1149 nz4(jlt_new) = nz4(i)
1150 p1(jlt_new) = p1(i)
1151 p2(jlt_new) = p2(i)
1152 p3(jlt_new) = p3(i)
1153 p4(jlt_new) = p4(i)
1154 lb1(jlt_new) = lb1(i)
1155 lb2(jlt_new) = lb2(i)
1156 lb3(jlt_new) = lb3(i)
1157 lb4(jlt_new) = lb4(i)
1158 lc1(jlt_new) = lc1(i)
1159 lc2(jlt_new) = lc2(i)
1160 lc3(jlt_new) = lc3(i)
1161 lc4(jlt_new) = lc4(i)
1162 stif(jlt_new) = stif(i)
1163 gapv(jlt_new) = gapv(i)
1164 index(jlt_new)= index(i)
1165 kini(jlt_new) = kini(i)
1166 vxi(jlt_new) = vxi(i)
1167 vyi(jlt_new) = vyi(i)
1168 vzi(jlt_new) = vzi(i)
1169 msi(jlt_new) = msi(i)
1170C
1171 xi(jlt_new) = xi(i)
1172 yi(jlt_new) = yi(i)
1173 zi(jlt_new) = zi(i)
1174 x1(jlt_new) = x1(i)
1175 y1(jlt_new) = y1(i)
1176 z1(jlt_new) = z1(i)
1177 x2(jlt_new) = x2(i)
1178 y2(jlt_new) = y2(i)
1179 z2(jlt_new) = z2(i)
1180 x3(jlt_new) = x3(i)
1181 y3(jlt_new) = y3(i)
1182 z3(jlt_new) = z3(i)
1183 x4(jlt_new) = x4(i)
1184 y4(jlt_new) = y4(i)
1185 z4(jlt_new) = z4(i)
1186 rcurvi(jlt_new) = rcurvi(i)
1187 anglmi(jlt_new) = anglmi(i)
1188C
1189 nsms(jlt_new) = nsms(i)
1190C
1191 keep(jlt_new) = keep(i)
1192C
1193 ENDIF
1194 ENDDO
1195 ENDIF
1196 ELSE
1197 IF(iadm == 0 )THEN
1198 DO i=1,jlt
1199 IF(keep(i)/=0)THEN
1200 jlt_new = jlt_new + 1
1201 cn_loc(jlt_new) = cand_n(i)
1202 ce_loc(jlt_new) = cand_e(i)
1203 ix1(jlt_new) = ix1(i)
1204 ix2(jlt_new) = ix2(i)
1205 ix3(jlt_new) = ix3(i)
1206 ix4(jlt_new) = ix4(i)
1207 nsvg(jlt_new) = nsvg(i)
1208 nx1(jlt_new) = nx1(i)
1209 nx2(jlt_new) = nx2(i)
1210 nx3(jlt_new) = nx3(i)
1211 nx4(jlt_new) = nx4(i)
1212 ny1(jlt_new) = ny1(i)
1213 ny2(jlt_new) = ny2(i)
1214 ny3(jlt_new) = ny3(i)
1215 ny4(jlt_new) = ny4(i)
1216 nz1(jlt_new) = nz1(i)
1217 nz2(jlt_new) = nz2(i)
1218 nz3(jlt_new) = nz3(i)
1219 nz4(jlt_new) = nz4(i)
1220 p1(jlt_new) = p1(i)
1221 p2(jlt_new) = p2(i)
1222 p3(jlt_new) = p3(i)
1223 p4(jlt_new) = p4(i)
1224 lb1(jlt_new) = lb1(i)
1225 lb2(jlt_new) = lb2(i)
1226 lb3(jlt_new) = lb3(i)
1227 lb4(jlt_new) = lb4(i)
1228 lc1(jlt_new) = lc1(i)
1229 lc2(jlt_new) = lc2(i)
1230 lc3(jlt_new) = lc3(i)
1231 lc4(jlt_new) = lc4(i)
1232 stif(jlt_new) = stif(i)
1233 gapv(jlt_new) = gapv(i)
1234 index(jlt_new)= index(i)
1235C
1236 kini(jlt_new) = kini(i)
1237 vxi(jlt_new) = vxi(i)
1238 vyi(jlt_new) = vyi(i)
1239 vzi(jlt_new) = vzi(i)
1240 msi(jlt_new) = msi(i)
1241C
1242 xi(jlt_new) = xi(i)
1243 yi(jlt_new) = yi(i)
1244 zi(jlt_new) = zi(i)
1245 x1(jlt_new) = x1(i)
1246 y1(jlt_new) = y1(i)
1247 z1(jlt_new) = z1(i)
1248 x2(jlt_new) = x2(i)
1249 y2(jlt_new) = y2(i)
1250 z2(jlt_new) = z2(i)
1251 x3(jlt_new) = x3(i)
1252 y3(jlt_new) = y3(i)
1253 z3(jlt_new) = z3(i)
1254 x4(jlt_new) = x4(i)
1255 y4(jlt_new) = y4(i)
1256 z4(jlt_new) = z4(i)
1257C
1258 tempi(jlt_new) = tempi(i)
1259 phi(jlt_new) = zero
1260 areasi(jlt_new) = areasi(i)
1261 ieleci(jlt_new) =ieleci(i)
1262C
1263 nsms(jlt_new) = nsms(i)
1264C
1265 keep(jlt_new) = keep(i)
1266C
1267 ENDIF
1268 ENDDO
1269 ELSE !
1270 DO i=1,jlt
1271 IF(keep(i)/=0)THEN
1272 jlt_new = jlt_new + 1
1273 cn_loc(jlt_new) = cand_n(i)
1274 ce_loc(jlt_new) = cand_e(i)
1275 ix1(jlt_new) = ix1(i)
1276 ix2(jlt_new) = ix2(i)
1277 ix3(jlt_new) = ix3(i)
1278 ix4(jlt_new) = ix4(i)
1279 nsvg(jlt_new) = nsvg(i)
1280 nx1(jlt_new) = nx1(i)
1281 nx2(jlt_new) = nx2(i)
1282 nx3(jlt_new) = nx3(i)
1283 nx4(jlt_new) = nx4(i)
1284 ny1(jlt_new) = ny1(i)
1285 ny2(jlt_new) = ny2(i)
1286 ny3(jlt_new) = ny3(i)
1287 ny4(jlt_new) = ny4(i)
1288 nz1(jlt_new) = nz1(i)
1289 nz2(jlt_new) = nz2(i)
1290 nz3(jlt_new) = nz3(i)
1291 nz4(jlt_new) = nz4(i)
1292 p1(jlt_new) = p1(i)
1293 p2(jlt_new) = p2(i)
1294 p3(jlt_new) = p3(i)
1295 p4(jlt_new) = p4(i)
1296 lb1(jlt_new) = lb1(i)
1297 lb2(jlt_new) = lb2(i)
1298 lb3(jlt_new) = lb3(i)
1299 lb4(jlt_new) = lb4(i)
1300 lc1(jlt_new) = lc1(i)
1301 lc2(jlt_new) = lc2(i)
1302 lc3(jlt_new) = lc3(i)
1303 lc4(jlt_new) = lc4(i)
1304 stif(jlt_new) = stif(i)
1305 gapv(jlt_new) = gapv(i)
1306 index(jlt_new)= index(i)
1307 kini(jlt_new) = kini(i)
1308 vxi(jlt_new) = vxi(i)
1309 vyi(jlt_new) = vyi(i)
1310 vzi(jlt_new) = vzi(i)
1311 msi(jlt_new) = msi(i)
1312 xi(jlt_new) = xi(i)
1313 yi(jlt_new) = yi(i)
1314 zi(jlt_new) = zi(i)
1315 x1(jlt_new) = x1(i)
1316 y1(jlt_new) = y1(i)
1317 z1(jlt_new) = z1(i)
1318 x2(jlt_new) = x2(i)
1319 y2(jlt_new) = y2(i)
1320 z2(jlt_new) = z2(i)
1321 x3(jlt_new) = x3(i)
1322 y3(jlt_new) = y3(i)
1323 z3(jlt_new) = z3(i)
1324 x4(jlt_new) = x4(i)
1325 y4(jlt_new) = y4(i)
1326 z4(jlt_new) = z4(i)
1327 rcurvi(jlt_new) = rcurvi(i)
1328 anglmi(jlt_new) = anglmi(i)
1329 tempi(jlt_new) = tempi(i)
1330 phi(jlt_new) = zero
1331 areasi(jlt_new) = areasi(i)
1332 ieleci(jlt_new) =ieleci(i)
1333C
1334 nsms(jlt_new) = nsms(i)
1335C
1336 keep(jlt_new) = keep(i)
1337C
1338 ENDIF
1339 ENDDO
1340 ENDIF
1341 ENDIF
1342 ELSE ! INTFRIC
1343 IF(intth == 0 ) THEN
1344 IF(icurv==0.AND.iadm==0 )THEN
1345 DO i=1,jlt
1346 IF(keep(i)/=0)THEN
1347 jlt_new = jlt_new + 1
1348 cn_loc(jlt_new) = cand_n(i)
1349 ce_loc(jlt_new) = cand_e(i)
1350 ix1(jlt_new) = ix1(i)
1351 ix2(jlt_new) = ix2(i)
1352 ix3(jlt_new) = ix3(i)
1353 ix4(jlt_new) = ix4(i)
1354 nsvg(jlt_new) = nsvg(i)
1355 nx1(jlt_new) = nx1(i)
1356 nx2(jlt_new) = nx2(i)
1357 nx3(jlt_new) = nx3(i)
1358 nx4(jlt_new) = nx4(i)
1359 ny1(jlt_new) = ny1(i)
1360 ny2(jlt_new) = ny2(i)
1361 ny3(jlt_new) = ny3(i)
1362 ny4(jlt_new) = ny4(i)
1363 nz1(jlt_new) = nz1(i)
1364 nz2(jlt_new) = nz2(i)
1365 nz3(jlt_new) = nz3(i)
1366 nz4(jlt_new) = nz4(i)
1367 p1(jlt_new) = p1(i)
1368 p2(jlt_new) = p2(i)
1369 p3(jlt_new) = p3(i)
1370 p4(jlt_new) = p4(i)
1371 lb1(jlt_new) = lb1(i)
1372 lb2(jlt_new) = lb2(i)
1373 lb3(jlt_new) = lb3(i)
1374 lb4(jlt_new) = lb4(i)
1375 lc1(jlt_new) = lc1(i)
1376 lc2(jlt_new) = lc2(i)
1377 lc3(jlt_new) = lc3(i)
1378 lc4(jlt_new) = lc4(i)
1379 stif(jlt_new) = stif(i)
1380 gapv(jlt_new) = gapv(i)
1381 index(jlt_new)= index(i)
1382C
1383 kini(jlt_new) = kini(i)
1384 vxi(jlt_new) = vxi(i)
1385 vyi(jlt_new) = vyi(i)
1386 vzi(jlt_new) = vzi(i)
1387 msi(jlt_new) = msi(i)
1388C
1389 xi(jlt_new) = xi(i)
1390 yi(jlt_new) = yi(i)
1391 zi(jlt_new) = zi(i)
1392 x1(jlt_new) = x1(i)
1393 y1(jlt_new) = y1(i)
1394 z1(jlt_new) = z1(i)
1395 x2(jlt_new) = x2(i)
1396 y2(jlt_new) = y2(i)
1397 z2(jlt_new) = z2(i)
1398 x3(jlt_new) = x3(i)
1399 y3(jlt_new) = y3(i)
1400 z3(jlt_new) = z3(i)
1401 x4(jlt_new) = x4(i)
1402 y4(jlt_new) = y4(i)
1403 z4(jlt_new) = z4(i)
1404C
1405 nsms(jlt_new) = nsms(i)
1406C
1407 ipartfricsi(jlt_new)=ipartfricsi(i)
1408 ipartfricmi(jlt_new)=ipartfricmi(i)
1409 IF(iorthfric_loc >0) THEN
1410 irep_fricmi(jlt_new)= irep_fricmi(i)
1411 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
1412 ENDIF
1413C
1414 keep(jlt_new) = keep(i)
1415C
1416 ENDIF
1417 ENDDO
1418 ELSE !
1419 DO i=1,jlt
1420 IF(keep(i)/=0)THEN
1421 jlt_new = jlt_new + 1
1422 cn_loc(jlt_new) = cand_n(i)
1423 ce_loc(jlt_new) = cand_e(i)
1424 ix1(jlt_new) = ix1(i)
1425 ix2(jlt_new) = ix2(i)
1426 ix3(jlt_new) = ix3(i)
1427 ix4(jlt_new) = ix4(i)
1428 nsvg(jlt_new) = nsvg(i)
1429 nx1(jlt_new) = nx1(i)
1430 nx2(jlt_new) = nx2(i)
1431 nx3(jlt_new) = nx3(i)
1432 nx4(jlt_new) = nx4(i)
1433 ny1(jlt_new) = ny1(i)
1434 ny2(jlt_new) = ny2(i)
1435 ny3(jlt_new) = ny3(i)
1436 ny4(jlt_new) = ny4(i)
1437 nz1(jlt_new) = nz1(i)
1438 nz2(jlt_new) = nz2(i)
1439 nz3(jlt_new) = nz3(i)
1440 nz4(jlt_new) = nz4(i)
1441 p1(jlt_new) = p1(i)
1442 p2(jlt_new) = p2(i)
1443 p3(jlt_new) = p3(i)
1444 p4(jlt_new) = p4(i)
1445 lb1(jlt_new) = lb1(i)
1446 lb2(jlt_new) = lb2(i)
1447 lb3(jlt_new) = lb3(i)
1448 lb4(jlt_new) = lb4(i)
1449 lc1(jlt_new) = lc1(i)
1450 lc2(jlt_new) = lc2(i)
1451 lc3(jlt_new) = lc3(i)
1452 lc4(jlt_new) = lc4(i)
1453 stif(jlt_new) = stif(i)
1454 gapv(jlt_new) = gapv(i)
1455 index(jlt_new)= index(i)
1456 kini(jlt_new) = kini(i)
1457 vxi(jlt_new) = vxi(i)
1458 vyi(jlt_new) = vyi(i)
1459 vzi(jlt_new) = vzi(i)
1460 msi(jlt_new) = msi(i)
1461C
1462 xi(jlt_new) = xi(i)
1463 yi(jlt_new) = yi(i)
1464 zi(jlt_new) = zi(i)
1465 x1(jlt_new) = x1(i)
1466 y1(jlt_new) = y1(i)
1467 z1(jlt_new) = z1(i)
1468 x2(jlt_new) = x2(i)
1469 y2(jlt_new) = y2(i)
1470 z2(jlt_new) = z2(i)
1471 x3(jlt_new) = x3(i)
1472 y3(jlt_new) = y3(i)
1473 z3(jlt_new) = z3(i)
1474 x4(jlt_new) = x4(i)
1475 y4(jlt_new) = y4(i)
1476 z4(jlt_new) = z4(i)
1477 rcurvi(jlt_new) = rcurvi(i)
1478 anglmi(jlt_new) = anglmi(i)
1479C
1480 nsms(jlt_new) = nsms(i)
1481C
1482 ipartfricsi(jlt_new)=ipartfricsi(i)
1483 ipartfricmi(jlt_new)=ipartfricmi(i)
1484 IF(iorthfric_loc >0) THEN
1485 irep_fricmi(jlt_new)= irep_fricmi(i)
1486 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
1487 ENDIF
1488C
1489 keep(jlt_new) = keep(i)
1490C
1491 ENDIF
1492 ENDDO
1493 ENDIF
1494 ELSE
1495 IF(iadm == 0 )THEN
1496 DO i=1,jlt
1497 IF(keep(i)/=0)THEN
1498 jlt_new = jlt_new + 1
1499 cn_loc(jlt_new) = cand_n(i)
1500 ce_loc(jlt_new) = cand_e(i)
1501 ix1(jlt_new) = ix1(i)
1502 ix2(jlt_new) = ix2(i)
1503 ix3(jlt_new) = ix3(i)
1504 ix4(jlt_new) = ix4(i)
1505 nsvg(jlt_new) = nsvg(i)
1506 nx1(jlt_new) = nx1(i)
1507 nx2(jlt_new) = nx2(i)
1508 nx3(jlt_new) = nx3(i)
1509 nx4(jlt_new) = nx4(i)
1510 ny1(jlt_new) = ny1(i)
1511 ny2(jlt_new) = ny2(i)
1512 ny3(jlt_new) = ny3(i)
1513 ny4(jlt_new) = ny4(i)
1514 nz1(jlt_new) = nz1(i)
1515 nz2(jlt_new) = nz2(i)
1516 nz3(jlt_new) = nz3(i)
1517 nz4(jlt_new) = nz4(i)
1518 p1(jlt_new) = p1(i)
1519 p2(jlt_new) = p2(i)
1520 p3(jlt_new) = p3(i)
1521 p4(jlt_new) = p4(i)
1522 lb1(jlt_new) = lb1(i)
1523 lb2(jlt_new) = lb2(i)
1524 lb3(jlt_new) = lb3(i)
1525 lb4(jlt_new) = lb4(i)
1526 lc1(jlt_new) = lc1(i)
1527 lc2(jlt_new) = lc2(i)
1528 lc3(jlt_new) = lc3(i)
1529 lc4(jlt_new) = lc4(i)
1530 stif(jlt_new) = stif(i)
1531 gapv(jlt_new) = gapv(i)
1532 index(jlt_new)= index(i)
1533C
1534 kini(jlt_new) = kini(i)
1535 vxi(jlt_new) = vxi(i)
1536 vyi(jlt_new) = vyi(i)
1537 vzi(jlt_new) = vzi(i)
1538 msi(jlt_new) = msi(i)
1539C
1540 xi(jlt_new) = xi(i)
1541 yi(jlt_new) = yi(i)
1542 zi(jlt_new) = zi(i)
1543 x1(jlt_new) = x1(i)
1544 y1(jlt_new) = y1(i)
1545 z1(jlt_new) = z1(i)
1546 x2(jlt_new) = x2(i)
1547 y2(jlt_new) = y2(i)
1548 z2(jlt_new) = z2(i)
1549 x3(jlt_new) = x3(i)
1550 y3(jlt_new) = y3(i)
1551 z3(jlt_new) = z3(i)
1552 x4(jlt_new) = x4(i)
1553 y4(jlt_new) = y4(i)
1554 z4(jlt_new) = z4(i)
1555C
1556 tempi(jlt_new) = tempi(i)
1557 phi(jlt_new) = zero
1558 areasi(jlt_new) = areasi(i)
1559 ieleci(jlt_new) =ieleci(i)
1560 nsms(jlt_new) = nsms(i)
1561 ipartfricsi(jlt_new)=ipartfricsi(i)
1562 ipartfricmi(jlt_new)=ipartfricmi(i)
1563 IF(iorthfric_loc >0) THEN
1564 irep_fricmi(jlt_new)= irep_fricmi(i)
1565 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
1566 ENDIF
1567C
1568 keep(jlt_new) = keep(i)
1569C
1570 ENDIF
1571 ENDDO
1572 ELSE !
1573 DO i=1,jlt
1574 IF(keep(i)/=0)THEN
1575 jlt_new = jlt_new + 1
1576 cn_loc(jlt_new) = cand_n(i)
1577 ce_loc(jlt_new) = cand_e(i)
1578 ix1(jlt_new) = ix1(i)
1579 ix2(jlt_new) = ix2(i)
1580 ix3(jlt_new) = ix3(i)
1581 ix4(jlt_new) = ix4(i)
1582 nsvg(jlt_new) = nsvg(i)
1583 nx1(jlt_new) = nx1(i)
1584 nx2(jlt_new) = nx2(i)
1585 nx3(jlt_new) = nx3(i)
1586 nx4(jlt_new) = nx4(i)
1587 ny1(jlt_new) = ny1(i)
1588 ny2(jlt_new) = ny2(i)
1589 ny3(jlt_new) = ny3(i)
1590 ny4(jlt_new) = ny4(i)
1591 nz1(jlt_new) = nz1(i)
1592 nz2(jlt_new) = nz2(i)
1593 nz3(jlt_new) = nz3(i)
1594 nz4(jlt_new) = nz4(i)
1595 p1(jlt_new) = p1(i)
1596 p2(jlt_new) = p2(i)
1597 p3(jlt_new) = p3(i)
1598 p4(jlt_new) = p4(i)
1599 lb1(jlt_new) = lb1(i)
1600 lb2(jlt_new) = lb2(i)
1601 lb3(jlt_new) = lb3(i)
1602 lb4(jlt_new) = lb4(i)
1603 lc1(jlt_new) = lc1(i)
1604 lc2(jlt_new) = lc2(i)
1605 lc3(jlt_new) = lc3(i)
1606 lc4(jlt_new) = lc4(i)
1607 stif(jlt_new) = stif(i)
1608 gapv(jlt_new) = gapv(i)
1609 index(jlt_new)= index(i)
1610 kini(jlt_new) = kini(i)
1611 vxi(jlt_new) = vxi(i)
1612 vyi(jlt_new) = vyi(i)
1613 vzi(jlt_new) = vzi(i)
1614 msi(jlt_new) = msi(i)
1615 xi(jlt_new) = xi(i)
1616 yi(jlt_new) = yi(i)
1617 zi(jlt_new) = zi(i)
1618 x1(jlt_new) = x1(i)
1619 y1(jlt_new) = y1(i)
1620 z1(jlt_new) = z1(i)
1621 x2(jlt_new) = x2(i)
1622 y2(jlt_new) = y2(i)
1623 z2(jlt_new) = z2(i)
1624 x3(jlt_new) = x3(i)
1625 y3(jlt_new) = y3(i)
1626 z3(jlt_new) = z3(i)
1627 x4(jlt_new) = x4(i)
1628 y4(jlt_new) = y4(i)
1629 z4(jlt_new) = z4(i)
1630 rcurvi(jlt_new) = rcurvi(i)
1631 anglmi(jlt_new) = anglmi(i)
1632 tempi(jlt_new) = tempi(i)
1633 phi(jlt_new) = zero
1634 areasi(jlt_new) = areasi(i)
1635 ieleci(jlt_new) =ieleci(i)
1636C
1637 nsms(jlt_new) = nsms(i)
1638C
1639 ipartfricsi(jlt_new)=ipartfricsi(i)
1640 ipartfricmi(jlt_new)=ipartfricmi(i)
1641 IF(iorthfric_loc >0) THEN
1642 irep_fricmi(jlt_new)= irep_fricmi(i)
1643 dir_fricmi(jlt_new,1:2) = dir_fricmi(i,1:2)
1644 ENDIF
1645C
1646 keep(jlt_new) = keep(i)
1647C
1648 ENDIF
1649 ENDDO
1650 ENDIF
1651 ENDIF
1652 ENDIF
1653 END IF
1654C
1655 IF(itied/=0.AND.jlt_tied/=0)THEN
1656C---------------------------------------------
1657C Tri a bulle
1658C---------------------------------------------
1659 iterate=.true.
1660 DO WHILE (iterate .EQV. .true.)
1661 iterate=.false.
1662 DO j=1,jlt_new-1
1663 IF(keep(j) > keep(j+1) )THEN
1664
1665 iperm =keep(j)
1666 keep(j) =keep(j+1)
1667 keep(j+1) =iperm
1668
1669 iperm =cn_loc(j)
1670 cn_loc(j) =cn_loc(j+1)
1671 cn_loc(j+1)=iperm
1672 iperm =ce_loc(j)
1673 ce_loc(j) =ce_loc(j+1)
1674 ce_loc(j+1)=iperm
1675 iperm =ix1(j)
1676 ix1(j) =ix1(j+1)
1677 ix1(j+1) =iperm
1678 iperm =ix2(j)
1679 ix2(j) =ix2(j+1)
1680 ix2(j+1) =iperm
1681 iperm =ix3(j)
1682 ix3(j) =ix3(j+1)
1683 ix3(j+1) =iperm
1684 iperm =ix4(j)
1685 ix4(j) =ix4(j+1)
1686 ix4(j+1) =iperm
1687 iperm =nsvg(j)
1688 nsvg(j) =nsvg(j+1)
1689 nsvg(j+1) =iperm
1690 rperm =nx1(j)
1691 nx1(j) =nx1(j+1)
1692 nx1(j+1) =rperm
1693 rperm =nx2(j)
1694 nx2(j) =nx2(j+1)
1695 nx2(j+1) =rperm
1696 rperm =nx3(j)
1697 nx3(j) =nx3(j+1)
1698 nx3(j+1) =rperm
1699 rperm =nx4(j)
1700 nx4(j) =nx4(j+1)
1701 nx4(j+1) =rperm
1702 rperm =ny1(j)
1703 ny1(j) =ny1(j+1)
1704 ny1(j+1) =rperm
1705 rperm =ny2(j)
1706 ny2(j) =ny2(j+1)
1707 ny2(j+1) =rperm
1708 rperm =ny3(j)
1709 ny3(j) =ny3(j+1)
1710 ny3(j+1) =rperm
1711 rperm =ny4(j)
1712 ny4(j) =ny4(j+1)
1713 ny4(j+1) =rperm
1714 rperm =nz1(j)
1715 nz1(j) =nz1(j+1)
1716 nz1(j+1) =rperm
1717 rperm =nz2(j)
1718 nz2(j) =nz2(j+1)
1719 nz2(j+1) =rperm
1720 rperm =nz3(j)
1721 nz3(j) =nz3(j+1)
1722 nz3(j+1) =rperm
1723 rperm =nz4(j)
1724 nz4(j) =nz4(j+1)
1725 nz4(j+1) =rperm
1726 rperm =p1(j)
1727 p1(j) =p1(j+1)
1728 p1(j+1) =rperm
1729 rperm =p2(j)
1730 p2(j) =p2(j+1)
1731 p2(j+1) =rperm
1732 rperm =p3(j)
1733 p3(j) =p3(j+1)
1734 p3(j+1) =rperm
1735 rperm =p4(j)
1736 p4(j) =p4(j+1)
1737 p4(j+1) =rperm
1738 rperm =lb1(j)
1739 lb1(j) =lb1(j+1)
1740 lb1(j+1) =rperm
1741 rperm =lb2(j)
1742 lb2(j) =lb2(j+1)
1743 lb2(j+1) =rperm
1744 rperm =lb3(j)
1745 lb3(j) =lb3(j+1)
1746 lb3(j+1) =rperm
1747 rperm =lb4(j)
1748 lb4(j) =lb4(j+1)
1749 lb4(j+1) =rperm
1750 rperm =lc1(j)
1751 lc1(j) =lc1(j+1)
1752 lc1(j+1) =rperm
1753 rperm =lc2(j)
1754 lc2(j) =lc2(j+1)
1755 lc2(j+1) =rperm
1756 rperm =lc3(j)
1757 lc3(j) =lc3(j+1)
1758 lc3(j+1) =rperm
1759 rperm =lc4(j)
1760 lc4(j) =lc4(j+1)
1761 lc4(j+1) =rperm
1762
1763 rperm =stif(j)
1764 stif(j) =stif(j+1)
1765 stif(j+1) =rperm
1766 rperm =gapv(j)
1767 gapv(j) =gapv(j+1)
1768 gapv(j+1) =rperm
1769
1770 iperm =index(j)
1771 index(j) =index(j+1)
1772 index(j+1) =iperm
1773 iperm =kini(j)
1774 kini(j) =kini(j+1)
1775 kini(j+1) =iperm
1776
1777 rperm =vxi(j)
1778 vxi(j) =vxi(j+1)
1779 vxi(j+1) =rperm
1780 rperm =vyi(j)
1781 vyi(j) =vyi(j+1)
1782 vyi(j+1) =rperm
1783 rperm =vzi(j)
1784 vzi(j) =vzi(j+1)
1785 vzi(j+1) =rperm
1786 rperm =msi(j)
1787 msi(j) =msi(j+1)
1788 msi(j+1) =rperm
1789 rperm =xi(j)
1790 xi(j) =xi(j+1)
1791 xi(j+1) =rperm
1792 rperm =yi(j)
1793 yi(j) =yi(j+1)
1794 yi(j+1) =rperm
1795 rperm =zi(j)
1796 zi(j) =zi(j+1)
1797 zi(j+1) =rperm
1798 rperm =x1(j)
1799 x1(j) =x1(j+1)
1800 x1(j+1) =rperm
1801 rperm =y1(j)
1802 y1(j) =y1(j+1)
1803 y1(j+1) =rperm
1804 rperm =z1(j)
1805 z1(j) =z1(j+1)
1806 z1(j+1) =rperm
1807 rperm =x2(j)
1808 x2(j) =x2(j+1)
1809 x2(j+1) =rperm
1810 rperm =y2(j)
1811 y2(j) =y2(j+1)
1812 y2(j+1) =rperm
1813 rperm =z2(j)
1814 z2(j) =z2(j+1)
1815 z2(j+1) =rperm
1816 rperm =x3(j)
1817 x3(j) =x3(j+1)
1818 x3(j+1) =rperm
1819 rperm =y3(j)
1820 y3(j) =y3(j+1)
1821 y3(j+1) =rperm
1822 rperm =z3(j)
1823 z3(j) =z3(j+1)
1824 z3(j+1) =rperm
1825 rperm =x4(j)
1826 x4(j) =x4(j+1)
1827 x4(j+1) =rperm
1828 rperm =y4(j)
1829 y4(j) =y4(j+1)
1830 y4(j+1) =rperm
1831 rperm =z4(j)
1832 z4(j) =z4(j+1)
1833 z4(j+1) =rperm
1834C
1835C
1836 IF(icurv/=0.OR.iadm/=0)THEN
1837 rperm =rcurvi(j)
1838 rcurvi(j) =rcurvi(j+1)
1839 rcurvi(j+1) =rperm
1840 rperm =anglmi(j)
1841 anglmi(j) =anglmi(j+1)
1842 anglmi(j+1) =rperm
1843 END IF
1844
1845 IF(intth/=0)THEN
1846 rperm =tempi(j)
1847 tempi(j) =tempi(j+1)
1848 tempi(j+1) =rperm
1849C PHI(J) = ZERO
1850 rperm =areasi(j)
1851 areasi(j) =areasi(j+1)
1852 areasi(j+1) =rperm
1853 iperm =ieleci(j)
1854 ieleci(j) =ieleci(j+1)
1855 ieleci(j+1) =iperm
1856 END IF
1857
1858C
1859 IF(intfric/=0)THEN
1860 iperm =ipartfricsi(j)
1861 ipartfricsi(j) =ipartfricsi(j+1)
1862 ipartfricsi(j+1) =iperm
1863 iperm =ipartfricmi(j)
1864 ipartfricmi(j) =ipartfricmi(j+1)
1865 ipartfricmi(j+1) =iperm
1866 IF(iorthfric_loc >0) THEN
1867 iperm =irep_fricmi(j)
1868 irep_fricmi(j) = irep_fricmi(j+1)
1869 irep_fricmi(j+1)= iperm
1870 rperm = dir_fricmi(j,1)
1871 dir_fricmi(j,1) = dir_fricmi(j+1,1)
1872 dir_fricmi(j+1,1)= rperm
1873 rperm = dir_fricmi(j,2)
1874 dir_fricmi(j,2) = dir_fricmi(j+1,2)
1875 dir_fricmi(j+1,2)= rperm
1876 ENDIF
1877 END IF
1878
1879 IF(idtmins==2)THEN
1880 iperm =nsms(j)
1881 nsms(j) =nsms(j+1)
1882 nsms(j+1) =iperm
1883 END IF
1884
1885 iterate = .true.
1886 ENDIF
1887 ENDDO
1888 ENDDO
1889 END IF ! IF(ITIED/=0.AND.JLT_TIED/=0)THEN
1890
1891C
1892 RETURN
1893 END
i7dst3
subroutine i7dst3(jlt, cand_n, cand_e, cn_loc, ce_loc, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, ix1, ix2, ix3, ix4, nsvg, stif, jlt_new, gapv, inacti, cand_p, index, vxi, vyi, vzi, msi, kini, icurv, irect, nnx1, nnx2, nnx3, nnx4, nny1, nny2, nny3, nny4, nnz1, nnz2, nnz3, nnz4, nod_normal, iadm, rcurvi, anglmi, intth, tempi, phi, areasi, ieleci, nsms, cmaj, drad2, intfric, ipartfricsi, ipartfricmi, itied, jlt_tied, cand_f, iorthfric, irep_fricmi, dir_fricmi, dgapload)
Definition i7dst3.F:52
i7curv
subroutine i7curv(jlt, pene, n1, n2, n3, gapv, x, nod_normal, ix1, ix2, ix3, ix4, h1, h2, h3, h4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi)
Definition i7curv.F:443
i7cmaj
subroutine i7cmaj(jlt, cmaj, irect, nod_normal, cand_e, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, nnx1, nnx2, nnx3, nnx4, nny1, nny2, nny3, nny4, nnz1, nnz2, nnz3, nnz4)
Definition i7curv.F:272
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21