34 use element_mod , only : nixs
35
36
37
38
39
40#include "implicit_f.inc"
41
42
43
44#include "scrcut_c.inc"
45
46
47
49 . xcut(7,*),xyz0(3,*),d(3,*),cbuf(*)
50 INTEGER ICUT(44,*),IXS(NIXS,*),ICBUF(*),NODCUT,NELCUT,LEN,NBF
51
52
53
54 INTEGER IC,ITYP,NUMTOT,NUMEL,
55 .MIC1,MIC2,,MIC4,MIC5,MIC6,MAC1,MAC2,MAC3
56 mic1=1
57 mic2=mic1+10*len
58 mic3=mic2+12*len
59 mic4=mic3+6*len
60 mic5=mic4+ncuts
61 mic6=mic5+ncuts
62 mac1=1
63 mac2=mac1+18*len
64 mac3=mac2+6*len
65 DO ic=1,ncuts
66 ityp=icut(1,ic)
67 CALL cutcon(ityp,xcut(1,ic),xcut(4,ic),xyz0,ixs,d,
68 . numtot,numel,
69 . icbuf(mic1),icbuf(mic2),icbuf(mic3),
70 . cbuf(mac1),cbuf(mac2),nodcut)
71 nodcut=nodcut+numtot
72 nelcut=nelcut+numel
73 icbuf(mic4)=nelcut+nbf
74 icbuf(mic5)=nodcut
75 mic1=mic1+5*numel
76 mic2=mic2+numtot*2
77 mic4=mic4+1
78 mic5=mic5+1
79 mac1=mac1+3*numtot
80 mac2=mac2+numtot
81 ENDDO
82 RETURN
subroutine cutcon(ityp, x0, vn, xyz0, ixs, d, numtot, numel, nc, nvois, na, xyz, al, nodcut)