OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i25ass_e2s.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "parit_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i25asse0_e2s (jlt, cs_loc, n1, n2, m1, m2, hs1, hs2, hm1, hm2, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, a, stifn, stif, nedge, nin, jtask, pene)
subroutine i25asse05_e2s (jlt, cs_loc, n1, n2, m1, m2, hs1, hs2, hm1, hm2, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, a, stifn, nedge, k1, k2, k3, k4, c1, c2, c3, c4, viscn, nin, jtask, pene)
subroutine i25asse2_e2s (jlt, cs_loc, n1, n2, m1, m2, hs1, hs2, hm1, hm2, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fskyi, isky, niskyfie, stif, nedge, nin, noint, pene, edge_id, tagip)
subroutine i25asse25_e2s (jlt, cs_loc, n1, n2, m1, m2, hs1, hs2, hm1, hm2, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, isky, niskyfie, nedge, k1, k2, k3, k4, c1, c2, c3, c4, nin, noint, pene, tagip)

Function/Subroutine Documentation

◆ i25asse05_e2s()

subroutine i25asse05_e2s ( integer jlt,
integer, dimension(*) cs_loc,
integer, dimension(*) n1,
integer, dimension(*) n2,
integer, dimension(*) m1,
integer, dimension(*) m2,
hs1,
hs2,
hm1,
hm2,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
a,
stifn,
integer nedge,
k1,
k2,
k3,
k4,
c1,
c2,
c3,
c4,
viscn,
integer nin,
integer jtask,
pene )

Definition at line 122 of file i25ass_e2s.F.

130C-----------------------------------------------
131C M o d u l e s
132C-----------------------------------------------
133 USE tri7box
134C-----------------------------------------------
135C I m p l i c i t T y p e s
136C-----------------------------------------------
137#include "implicit_f.inc"
138C-----------------------------------------------
139C D u m m y A r g u m e n t s
140C-----------------------------------------------
141 INTEGER JLT, NEDGE, NIN,
142 + CS_LOC(*),
143 + N1(*),N2(*),M1(*),M2(*),JTASK
144 my_real
145 . hs1(*),hs2(*),hm1(*),hm2(*),
146 . fx1(*),fy1(*),fz1(*),
147 . fx2(*),fy2(*),fz2(*),
148 . fx3(*),fy3(*),fz3(*),
149 . fx4(*),fy4(*),fz4(*),
150 . k1(*),k2(*),k3(*),k4(*),
151 . c1(*),c2(*),c3(*),c4(*),
152 . a(3,*), stifn(*), viscn(*), pene(*)
153C-----------------------------------------------
154C L o c a l V a r i a b l e s
155C-----------------------------------------------
156 INTEGER I, J1,NODFI,ISHIFT
157C-----------------------------------------------
158C
159 nodfi = nlskyfie(nin)
160 ishift = nodfi*(jtask-1)
161C
162 DO i=1,jlt
163 IF(pene(i)==zero) cycle
164C
165 IF(cs_loc(i)<=nedge) THEN
166 j1=n1(i)
167 a(1,j1)=a(1,j1)+fx1(i)
168 a(2,j1)=a(2,j1)+fy1(i)
169 a(3,j1)=a(3,j1)+fz1(i)
170 stifn(j1)=stifn(j1)+k1(i)
171 viscn(j1)=viscn(j1)+c1(i)
172C
173 j1=n2(i)
174 a(1,j1)=a(1,j1)+fx2(i)
175 a(2,j1)=a(2,j1)+fy2(i)
176 a(3,j1)=a(3,j1)+fz2(i)
177 stifn(j1)=stifn(j1)+k2(i)
178 viscn(j1)=viscn(j1)+c2(i)
179 ELSE
180 j1=n1(i)
181 afie(nin)%P(1,j1+ishift)=afie(nin)%P(1,j1+ishift)+fx1(i)
182 afie(nin)%P(2,j1+ishift)=afie(nin)%P(2,j1+ishift)+fy1(i)
183 afie(nin)%P(3,j1+ishift)=afie(nin)%P(3,j1+ishift)+fz1(i)
184 stnfie(nin)%P(j1+ishift)=stnfie(nin)%P(j1+ishift)+k1(i)
185 vscfie(nin)%P(j1+ishift)=vscfie(nin)%P(j1+ishift)+c1(i)
186C
187 j1=n2(i)
188 afie(nin)%P(1,j1+ishift)=afie(nin)%P(1,j1+ishift)+fx2(i)
189 afie(nin)%P(2,j1+ishift)=afie(nin)%P(2,j1+ishift)+fy2(i)
190 afie(nin)%P(3,j1+ishift)=afie(nin)%P(3,j1+ishift)+fz2(i)
191 stnfie(nin)%P(j1+ishift)=stnfie(nin)%P(j1+ishift)+k2(i)
192 vscfie(nin)%P(j1+ishift)=vscfie(nin)%P(j1+ishift)+c2(i)
193 END IF
194 END DO
195C
196 DO i=1,jlt
197 IF(pene(i)==zero) cycle
198C
199 j1=m1(i)
200 a(1,j1)=a(1,j1)+fx3(i)
201 a(2,j1)=a(2,j1)+fy3(i)
202 a(3,j1)=a(3,j1)+fz3(i)
203 stifn(j1)=stifn(j1)+k3(i)
204 viscn(j1)=viscn(j1)+c3(i)
205C
206 j1=m2(i)
207 a(1,j1)=a(1,j1)+fx4(i)
208 a(2,j1)=a(2,j1)+fy4(i)
209 a(3,j1)=a(3,j1)+fz4(i)
210 stifn(j1)=stifn(j1)+k4(i)
211 viscn(j1)=viscn(j1)+c4(i)
212 ENDDO
213C
214 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
tri7box::stnfie
type(real_pointer), dimension(:), allocatable stnfie
Definition tri7box.F:449
tri7box::afie
type(real_pointer2), dimension(:), allocatable afie
Definition tri7box.F:459
tri7box::vscfie
type(real_pointer), dimension(:), allocatable vscfie
Definition tri7box.F:449
tri7box::nlskyfie
integer, dimension(:), allocatable nlskyfie
Definition tri7box.F:512

◆ i25asse0_e2s()

subroutine i25asse0_e2s ( integer jlt,
integer, dimension(*) cs_loc,
integer, dimension(*) n1,
integer, dimension(*) n2,
integer, dimension(*) m1,
integer, dimension(*) m2,
hs1,
hs2,
hm1,
hm2,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
a,
stifn,
stif,
integer nedge,
integer nin,
integer jtask,
pene )

Definition at line 30 of file i25ass_e2s.F.

36C-----------------------------------------------
37C M o d u l e s
38C-----------------------------------------------
39 USE tri7box
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 INTEGER JLT, NEDGE, NIN,
48 + CS_LOC(*),
49 + N1(*),N2(*),M1(*),M2(*),JTASK
50 my_real
51 . hs1(*),hs2(*),hm1(*),hm2(*),
52 . fx1(*),fy1(*),fz1(*),
53 . fx2(*),fy2(*),fz2(*),
54 . fx3(*),fy3(*),fz3(*),
55 . fx4(*),fy4(*),fz4(*),
56 . a(3,*), stifn(*), stif(*), pene(*)
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER I, J1,NODFI,ISHIFT
61C-----------------------------------------------
62C
63 nodfi = nlskyfie(nin)
64 ishift = nodfi*(jtask-1)
65C
66 DO i=1,jlt
67 IF(pene(i)==zero) cycle
68C
69 IF(cs_loc(i)<=nedge) THEN
70 j1=n1(i)
71 a(1,j1)=a(1,j1)+fx1(i)
72 a(2,j1)=a(2,j1)+fy1(i)
73 a(3,j1)=a(3,j1)+fz1(i)
74 stifn(j1) = stifn(j1) + stif(i)*abs(hs1(i))
75C
76 j1=n2(i)
77 a(1,j1)=a(1,j1)+fx2(i)
78 a(2,j1)=a(2,j1)+fy2(i)
79 a(3,j1)=a(3,j1)+fz2(i)
80 stifn(j1) = stifn(j1) + stif(i)*abs(hs2(i))
81 ELSE
82 j1=n1(i)
83 afie(nin)%P(1,j1+ishift)=afie(nin)%P(1,j1+ishift)+fx1(i)
84 afie(nin)%P(2,j1+ishift)=afie(nin)%P(2,j1+ishift)+fy1(i)
85 afie(nin)%P(3,j1+ishift)=afie(nin)%P(3,j1+ishift)+fz1(i)
86 stnfie(nin)%P(j1+ishift) = stnfie(nin)%P(j1+ishift) + stif(i)*abs(hs1(i))
87C
88 j1=n2(i)
89 afie(nin)%P(1,j1+ishift)=afie(nin)%P(1,j1+ishift)+fx2(i)
90 afie(nin)%P(2,j1+ishift)=afie(nin)%P(2,j1+ishift)+fy2(i)
91 afie(nin)%P(3,j1+ishift)=afie(nin)%P(3,j1+ishift)+fz2(i)
92 stnfie(nin)%P(j1+ishift) = stnfie(nin)%P(j1+ishift) + stif(i)*abs(hs2(i))
93 END IF
94 END DO
95C
96 DO i=1,jlt
97 IF(pene(i)==zero) cycle
98C
99 j1=m1(i)
100 a(1,j1)=a(1,j1)+fx3(i)
101 a(2,j1)=a(2,j1)+fy3(i)
102 a(3,j1)=a(3,j1)+fz3(i)
103 stifn(j1) = stifn(j1) + stif(i)*abs(hm1(i))
104C
105 j1=m2(i)
106 a(1,j1)=a(1,j1)+fx4(i)
107 a(2,j1)=a(2,j1)+fy4(i)
108 a(3,j1)=a(3,j1)+fz4(i)
109 stifn(j1) = stifn(j1) + stif(i)*abs(hm2(i))
110 ENDDO
111C
112 RETURN

◆ i25asse25_e2s()

subroutine i25asse25_e2s ( integer jlt,
integer, dimension(*) cs_loc,
integer, dimension(*) n1,
integer, dimension(*) n2,
integer, dimension(*) m1,
integer, dimension(*) m2,
hs1,
hs2,
hm1,
hm2,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
integer, dimension(*) isky,
integer niskyfie,
integer nedge,
k1,
k2,
k3,
k4,
c1,
c2,
c3,
c4,
integer nin,
integer noint,
pene,
integer, dimension(4*mvsiz), intent(in) tagip )

Definition at line 393 of file i25ass_e2s.F.

401C-----------------------------------------------
402C M o d u l e s
403C-----------------------------------------------
404 USE tri7box
405 USE message_mod
406C-----------------------------------------------
407C I m p l i c i t T y p e s
408C-----------------------------------------------
409#include "implicit_f.inc"
410#include "comlock.inc"
411C-----------------------------------------------
412C G l o b a l P a r a m e t e r s
413C-----------------------------------------------
414#include "mvsiz_p.inc"
415C-----------------------------------------------
416C C o m m o n B l o c k s
417C-----------------------------------------------
418#include "parit_c.inc"
419C-----------------------------------------------
420C D u m m y A r g u m e n t s
421C-----------------------------------------------
422 INTEGER JLT, NEDGE,NISKYFIE,NIN,NOINT,
423 + CS_LOC(*),ISKY(*),
424 + N1(*),N2(*),M1(*),M2(*)
425 my_real
426 . hs1(*),hs2(*),hm1(*),hm2(*),
427 . fx1(*),fy1(*),fz1(*),
428 . fx2(*),fy2(*),fz2(*),
429 . fx3(*),fy3(*),fz3(*),
430 . fx4(*),fy4(*),fz4(*),
431 . k1(*),k2(*),k3(*),k4(*),
432 . c1(*),c2(*),c3(*),c4(*),
433 . fskyi(lskyi,nfskyi), pene(*)
434 INTEGER , INTENT(IN) :: TAGIP(4*MVSIZ)
435C-----------------------------------------------
436C L o c a l V a r i a b l e s
437C-----------------------------------------------
438 INTEGER I, J1, NISKYL1, NISKYL,IGP,IGM, NISKYFIEL
439C
440 niskyl1 = 0
441 DO i = 1, jlt
442 IF(pene(i)==zero.AND.tagip(i)==0) cycle
443 IF (hm1(i)/=zero.OR.tagip(i)==1) niskyl1 = niskyl1 + 1
444 IF (hm2(i)/=zero.OR.tagip(i)==1) niskyl1 = niskyl1 + 1
445 ENDDO
446
447 igp = 0
448 igm = 0
449 DO i=1,jlt
450 IF(pene(i) == zero.AND.tagip(i)==0) cycle
451 IF(cs_loc(i)<=nedge) THEN
452 igp = igp+2
453 ELSE
454 igm = igm+1
455 ENDIF
456 ENDDO
457
458#include "lockon.inc"
459 niskyl = nisky
460 nisky = nisky + niskyl1 + igp
461 niskyfiel = niskyfie
462 niskyfie = niskyfie + igm
463#include "lockoff.inc"
464C
465 IF (niskyl+niskyl1+igp > lskyi) THEN
466 CALL ancmsg(msgid=26,anmode=aninfo)
467 CALL arret(2)
468 ENDIF
469 IF (niskyfiel+igm > nlskyfie(nin)) THEN
470 CALL ancmsg(msgid=26,anmode=aninfo)
471 CALL arret(2)
472 ENDIF
473C
474 DO i=1,jlt
475 IF(pene(i)==zero.AND.tagip(i)==0) cycle
476C
477 IF(cs_loc(i)<=nedge) THEN
478 niskyl = niskyl + 1
479 fskyi(niskyl,1)=fx1(i)
480 fskyi(niskyl,2)=fy1(i)
481 fskyi(niskyl,3)=fz1(i)
482 fskyi(niskyl,4)=k1(i)
483 fskyi(niskyl,5)=c1(i)
484 isky(niskyl) = n1(i)
485C
486 niskyl = niskyl + 1
487 fskyi(niskyl,1)=fx2(i)
488 fskyi(niskyl,2)=fy2(i)
489 fskyi(niskyl,3)=fz2(i)
490 fskyi(niskyl,4)=k2(i)
491 fskyi(niskyl,5)=c2(i)
492 isky(niskyl) = n2(i)
493 ELSE
494 niskyfiel = niskyfiel + 1
495 fskyfie(nin)%P(1,niskyfiel)=fx1(i)
496 fskyfie(nin)%P(2,niskyfiel)=fy1(i)
497 fskyfie(nin)%P(3,niskyfiel)=fz1(i)
498 fskyfie(nin)%P(4,niskyfiel)=k1(i)
499 fskyfie(nin)%P(5,niskyfiel)=c1(i)
500 fskyfie(nin)%P(6,niskyfiel)=fx2(i)
501 fskyfie(nin)%P(7,niskyfiel)=fy2(i)
502 fskyfie(nin)%P(8,niskyfiel)=fz2(i)
503 fskyfie(nin)%P(9,niskyfiel)=k2(i)
504 fskyfie(nin)%P(10,niskyfiel)=c2(i)
505 iskyfie(nin)%P(niskyfiel) = cs_loc(i)-nedge
506 END IF
507 END DO
508C
509 DO i=1,jlt
510 IF(pene(i)==zero.AND.tagip(i)==0) cycle
511C
512 IF (hm1(i)/=zero.OR.tagip(i)==1) THEN
513 niskyl = niskyl + 1
514 fskyi(niskyl,1)=fx3(i)
515 fskyi(niskyl,2)=fy3(i)
516 fskyi(niskyl,3)=fz3(i)
517 fskyi(niskyl,4)=k3(i)
518 fskyi(niskyl,5)=c3(i)
519 isky(niskyl) = m1(i)
520 ENDIF
521 ENDDO
522 DO i=1,jlt
523 IF(pene(i)==zero.AND.tagip(i)==0) cycle
524C
525 IF (hm2(i)/=zero.OR.tagip(i)==1) THEN
526 niskyl = niskyl + 1
527 fskyi(niskyl,1)=fx4(i)
528 fskyi(niskyl,2)=fy4(i)
529 fskyi(niskyl,3)=fz4(i)
530 fskyi(niskyl,4)=k4(i)
531 fskyi(niskyl,5)=c4(i)
532 isky(niskyl) = m2(i)
533 ENDIF
534 ENDDO
535C
536 RETURN
message_mod
Definition message_mod.F:1249
tri7box::fskyfie
type(real_pointer2), dimension(:), allocatable fskyfie
Definition tri7box.F:459
tri7box::iskyfie
type(int_pointer), dimension(:), allocatable iskyfie
Definition tri7box.F:480
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
arret
subroutine arret(nn)
Definition arret.F:87

◆ i25asse2_e2s()

subroutine i25asse2_e2s ( integer jlt,
integer, dimension(*) cs_loc,
integer, dimension(*) n1,
integer, dimension(*) n2,
integer, dimension(*) m1,
integer, dimension(*) m2,
hs1,
hs2,
hm1,
hm2,
fx1,
fy1,
fz1,
fx2,
fy2,
fz2,
fx3,
fy3,
fz3,
fx4,
fy4,
fz4,
fskyi,
integer, dimension(*) isky,
integer niskyfie,
stif,
integer nedge,
integer nin,
integer noint,
pene,
integer, dimension(2,4*mvsiz) edge_id,
integer, dimension(4*mvsiz), intent(in) tagip )

Definition at line 228 of file i25ass_e2s.F.

235C-----------------------------------------------
236C M o d u l e s
237C-----------------------------------------------
238 USE tri7box
239 USE message_mod
240C-----------------------------------------------
241C I m p l i c i t T y p e s
242C-----------------------------------------------
243#include "implicit_f.inc"
244#include "comlock.inc"
245C-----------------------------------------------
246C G l o b a l P a r a m e t e r s
247C-----------------------------------------------
248#include "mvsiz_p.inc"
249C-----------------------------------------------
250C C o m m o n B l o c k s
251C-----------------------------------------------
252#include "parit_c.inc"
253C-----------------------------------------------
254C D u m m y A r g u m e n t s
255C-----------------------------------------------
256 INTEGER :: EDGE_ID(2,4*MVSIZ)
257 INTEGER JLT, NEDGE,NISKYFIE,NIN,NOINT,
258 + CS_LOC(*),ISKY(*),
259 + N1(*),N2(*),M1(*),M2(*)
260 INTEGER , INTENT(IN) :: TAGIP(4*MVSIZ)
261 my_real
262 . hs1(*),hs2(*),hm1(*),hm2(*),
263 . fx1(*),fy1(*),fz1(*),
264 . fx2(*),fy2(*),fz2(*),
265 . fx3(*),fy3(*),fz3(*),
266 . fx4(*),fy4(*),fz4(*),
267 . fskyi(lskyi,nfskyi), stif(*), pene(*)
268C-----------------------------------------------
269C L o c a l V a r i a b l e s
270C-----------------------------------------------
271 INTEGER I, J1, NISKYL1, NISKYL,IGP,IGM, NISKYFIEL
272C
273 niskyl1 = 0
274 DO i = 1, jlt
275 IF(pene(i)==zero.AND.tagip(i)==0) cycle
276 IF (hm1(i)/=zero.OR.tagip(i)==1) niskyl1 = niskyl1 + 1
277 IF (hm2(i)/=zero.OR.tagip(i)==1) niskyl1 = niskyl1 + 1
278 ENDDO
279
280 igp = 0
281 igm = 0
282 DO i=1,jlt
283 IF(pene(i)==zero.AND.tagip(i)==0) cycle
284 IF(cs_loc(i)<=nedge) THEN
285 igp = igp+2
286 ELSE
287 igm = igm+1
288 ENDIF
289 ENDDO
290
291#include "lockon.inc"
292 niskyl = nisky
293 nisky = nisky + niskyl1 + igp
294 niskyfiel = niskyfie
295 niskyfie = niskyfie + igm
296#include "lockoff.inc"
297
298 IF (niskyl+niskyl1+igp > lskyi) THEN
299 CALL ancmsg(msgid=26,anmode=aninfo)
300 CALL arret(2)
301 ENDIF
302 IF (niskyfiel+igm > nlskyfie(nin)) THEN
303 CALL ancmsg(msgid=26,anmode=aninfo)
304 CALL arret(2)
305 ENDIF
306C
307 DO i=1,jlt
308 IF(pene(i)==zero.AND.tagip(i)==0) cycle
309C
310 IF(cs_loc(i)<=nedge) THEN
311 niskyl = niskyl + 1
312 fskyi(niskyl,1)=fx1(i)
313 fskyi(niskyl,2)=fy1(i)
314 fskyi(niskyl,3)=fz1(i)
315 fskyi(niskyl,4)=stif(i)*abs(hs1(i))
316 isky(niskyl) = n1(i)
317C
318 niskyl = niskyl + 1
319 fskyi(niskyl,1)=fx2(i)
320 fskyi(niskyl,2)=fy2(i)
321 fskyi(niskyl,3)=fz2(i)
322 fskyi(niskyl,4)=stif(i)*abs(hs2(i))
323 isky(niskyl) = n2(i)
324
325#ifdef WITH_ASSERT
326 WRITE(6,"(2I20,X,A,6Z20)") edge_id(1,i),edge_id(2,i),"A",fx1(i),fy1(i),fz1(i),fx2(i),fy2(i),fz2(i)
327#endif
328
329 ELSE
330 niskyfiel = niskyfiel + 1
331 fskyfie(nin)%P(1,niskyfiel)=fx1(i)
332 fskyfie(nin)%P(2,niskyfiel)=fy1(i)
333 fskyfie(nin)%P(3,niskyfiel)=fz1(i)
334 fskyfie(nin)%P(4,niskyfiel)=stif(i)*abs(hs1(i))
335 fskyfie(nin)%P(5,niskyfiel)=fx2(i)
336 fskyfie(nin)%P(6,niskyfiel)=fy2(i)
337 fskyfie(nin)%P(7,niskyfiel)=fz2(i)
338 fskyfie(nin)%P(8,niskyfiel)=stif(i)*abs(hs2(i))
339 iskyfie(nin)%P(niskyfiel) = cs_loc(i)-nedge
340C ASSERT(CS_LOC(I) - NEDGE > 0)
341
342#ifdef WITH_ASSERT
343 WRITE(6,"(2I20,X,A,6Z20)") edge_id(1,i),edge_id(2,i),"B",fx1(i),fy1(i),fz1(i),fx2(i),fy2(i),fz2(i)
344#endif
345 END IF
346 END DO
347C
348 DO i=1,jlt
349 IF(pene(i)==zero.AND.tagip(i)==0) cycle
350C
351 IF (hm1(i)/=zero.OR.tagip(i)==1) THEN
352 niskyl = niskyl + 1
353 fskyi(niskyl,1)=fx3(i)
354 fskyi(niskyl,2)=fy3(i)
355 fskyi(niskyl,3)=fz3(i)
356 fskyi(niskyl,4)=stif(i)*abs(hm1(i))
357 isky(niskyl) = m1(i)
358#ifdef WITH_ASSERT
359 WRITE(6,"(2I20,X,A,3Z20)") edge_id(1,i),edge_id(2,i),"C",fx3(i),fy3(i),fz3(i)
360#endif
361 ENDIF
362 ENDDO
363 DO i=1,jlt
364 IF(pene(i)==zero.AND.tagip(i)==0) cycle
365C
366 IF (hm2(i)/=zero.OR.tagip(i)==1) THEN
367 niskyl = niskyl + 1
368 fskyi(niskyl,1)=fx4(i)
369 fskyi(niskyl,2)=fy4(i)
370 fskyi(niskyl,3)=fz4(i)
371 fskyi(niskyl,4)=stif(i)*abs(hm2(i))
372 isky(niskyl) = m2(i)
373#ifdef WITH_ASSERT
374 WRITE(6,"(2I20,X,A,3Z20)") edge_id(1,i),edge_id(2,i),"D",fx4(i),fy4(i),fz4(i)
375#endif
376 ENDIF
377 ENDDO
378C
379 RETURN