32
33
34
36
37
38
39#include "implicit_f.inc"
40
41
42
43#include "param_c.inc"
44#include "com01_c.inc"
45#include "com04_c.inc"
46#include "com_xfem1.inc"
47
48
49
50 INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),LCNE_CRKXFEM
51 TYPE (XFEM_LVSET_) , DIMENSION(NLEVMAX) :: CRKLVSET
52 TYPE (XFEM_SKY_) , DIMENSION(NLEVMAX) :: CRKSKY
53 TYPE (XFEM_AVX_) , DIMENSION(NLEVMAX) :: CRKAVX
54 TYPE (XFEM_EDGE_) , DIMENSION(NXLAYMAX):: CRKEDGE
55 TYPE (XFEM_PHANTOM_), DIMENSION(NXLAYMAX):: XFEM_PHANTOM
56
57
58
59 INTEGER KNOD2ELC(NUMNOD+1)
60 INTEGER I,J,N,LFT,LLT,,NEL,ITY,IXFEM,NG
61
62 knod2elc = 0
63 DO ng=1,ngroup
64 nel=iparg(2,ng)
65 nft=iparg(3,ng)
66 ity=iparg(5,ng)
67 ixfem=iparg(54,ng)
68 lft=1
70 IF (ixfem == 0 ) cycle
71 IF (ity==3) THEN
72 DO i=lft,llt
73 DO j=2,5
74 n = ixc(j,i+nft)
75 knod2elc(n) = knod2elc(n) + 1
76 END DO
77 END DO
78 ELSE IF(ity==7) THEN
79 DO i=lft,llt
80 DO j=2,4
81 n = ixtg(j,i+nft)
82 knod2elc(n) = knod2elc(n) + 1
83 END DO
84 END DO
85 END IF
86 ENDDO
87
88 ienrnod = 0
89 DO i=1,numnod
90 ienrnod =
max(ienrnod,knod2elc(i))
91 END DO
92
93 ienrnod =
max(ienrnod,4)
94 ienrnod = ienrnod * nlevmax
95
96
97
98
99
100
101 DO i=1,nlevmax
102 ALLOCATE(crklvset(i)%ELCUT(ecrkxfe))
103 ALLOCATE(crklvset(i)%HFI0(4,ecrkxfe))
104 ALLOCATE(crklvset(i)%ENRICH0(4,ecrkxfe))
105 ALLOCATE(crklvset(i)%EDGE0(4,ecrkxfe))
106
107 ALLOCATE(crklvset(i)%EDGE(4,ecrkxfec))
108 ALLOCATE(crklvset(i)%EDGETG(3,ecrkxfetg))
109 ALLOCATE(crklvset(i)%ICUTEDGE(numedges))
110 ALLOCATE(crklvset(i)%RATIOEDGE(numedges))
111
112 ALLOCATE(crkavx(i)%A(3,lcne_crkxfem))
113 ALLOCATE(crkavx(i)%AR(3,lcne_crkxfem))
114 ALLOCATE(crkavx(i)%V(3,lcne_crkxfem))
115 ALLOCATE(crkavx(i)%VR(3,lcne_crkxfem))
116 ALLOCATE(crkavx(i)%X(3,lcne_crkxfem))
117 ALLOCATE(crkavx(i)%U(3,lcne_crkxfem))
118 ALLOCATE(crkavx(i)%XX(3,4*ecrkxfe))
119
120 ALLOCATE(crksky(i)%FSKY(8,lcne_crkxfem))
121
122 crklvset(i)%ELCUT= 0
123 crklvset(i)%HFI0 = 0
124 crklvset(i)%ENRICH0 = 0
125 crklvset(i)%EDGE0 = 0
126
127 crklvset(i)%EDGE= 0
128 crklvset(i)%EDGETG = 0
129 crklvset(i)%ICUTEDGE = 0
130 crklvset(i)%RATIOEDGE = 0
131
132 crkavx(i)%A = 0
133 crkavx(i)%AR = 0
134 crkavx(i)%V = 0
135 crkavx(i)%VR = 0
136 crkavx(i)%X = 0
137 crkavx(i)%U = 0
138 crkavx(i)%XX = 0
139 crksky(i)%FSKY = 0
140 ENDDO
141
142 DO i=1,nxlaymax
143 ALLOCATE(crkedge(i)%LAYCUT(ecrkxfe))
144 ALLOCATE(crkedge(i)%EDGEICRK(numedges))
145 ALLOCATE(crkedge(i)%EDGEIFI(2,numedges))
146 ALLOCATE(crkedge(i)%EDGEENR(2,numedges))
147 ALLOCATE(crkedge(i)%EDGETIP(2,numedges))
148 crkedge(i)%LAYCUT = 0
149 crkedge(i)%EDGEICRK = 0
150 crkedge(i)%EDGEIFI = 0
151 crkedge(i)%EDGEENR = 0
152 crkedge(i)%EDGETIP = 0
153 ENDDO
154
155 DO i=1,nxlaymax
156 ALLOCATE(xfem_phantom(i)%ELCUT(ecrkxfe))
157 ALLOCATE(xfem_phantom(i)%IFI(lcnecrkxfem))
158 ALLOCATE(xfem_phantom(i)%ITRI(2,ecrkxfe))
159 ALLOCATE(xfem_phantom(i)%TAGXP(1,1,1))
160 xfem_phantom(i)%ELCUT = 0
161 xfem_phantom(i)%TAGXP = 0
162 xfem_phantom(i)%IFI = 0
163 xfem_phantom(i)%ITRI = 0
164 ENDDO
165
166 RETURN
1.15.0