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