OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
w_bufbric_22.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!|| w_bufbric_22 ../starter/source/interfaces/inter3d1/w_bufbric_22.F
25!||--- called by ------------------------------------------------------
26!|| ddsplit ../starter/source/restart/ddsplit/ddsplit.F
27!||--- calls -----------------------------------------------------
28!|| arret ../starter/source/system/arret.F
29!||--- uses -----------------------------------------------------
30!||====================================================================
31 SUBROUTINE w_bufbric_22()
32C-----------------------------------------------
33C M o d u l e s
34C-----------------------------------------------
35 USE i22bufbric_mod
36 USE i22tri_mod
37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C C o m m o n B l o c k s
43C-----------------------------------------------
44#include "inter22.inc"
45#include "units_c.inc"
46C-----------------------------------------------
47C L o c a l V a r i a b l e s
48C-----------------------------------------------
49 my_real :: bid(1)
50 INTEGER :: NIN, I, J, K, IERR1, ISIZ
51 INTEGER :: Isectype(14)
52 CHARACTER*14 :: string
53C-----------------------------------------------
54C S o u r c e L i n e s
55C-----------------------------------------------
56
57 nb = 0
58
59 IF(.NOT.ALLOCATED(list_b_old))THEN
60 ALLOCATE(list_b_old(int22,i22len),stat=ierr1)
61 IF (ierr1/=0) GOTO 1000
62 list_b_old(:,:) = 0
63 ENDIF
64
65 IF(.NOT.ALLOCATED(nbold))THEN
66 ALLOCATE(nbold(int22) ,stat=ierr1)
67 IF (ierr1/=0) GOTO 1000
68 nbold(:) = 0
69 ENDIF
70
71 IF(.NOT.ALLOCATED(brick_list))THEN
72 ALLOCATE(brick_list(int22,i22len),stat=ierr1)
73 IF (ierr1/=0) GOTO 1000
74
75 brick_list(:,:)%ID = -1
76 brick_list(:,:)%NG = -2
77 brick_list(:,:)%IDLOC = -3
78 brick_list(:,:)%ICODE = -4
79 brick_list(:,:)%OLD_ICODE = -5
80 brick_list(:,:)%IDBLE = -6
81 brick_list(:,:)%NBITS = -7
82 brick_list(:,:)%NPQTS = -8
83 brick_list(:,:)%NBCUT = -9
84 brick_list(:,:)%MainID = -10
85 brick_list(:,:)%WasCut = -11
86 brick_list(:,:)%NewInBuffer = 0
87 brick_list(:,:)%OldMainStrongNode = -13
88 brick_list(:,:)%MLW = -14
89 brick_list(:,:)%ITASK = -15
90
91 DO j=1,8 ; brick_list(:,:)%SecId_Cell(j) = -16 ;ENDDO
92 DO j=1,8 ; brick_list(:,:)%NODE(j)%OLD_WhichCell = -17 ;ENDDO
93 DO j=1,8 ; brick_list(:,:)%NODE(j)%WhichCell = -18 ;ENDDO
94 DO j=1,9 ; brick_list(:,:)%POLY(j)%NumNOD = -19 ;ENDDO
95 DO j=1,9 ; brick_list(:,:)%POLY(j)%NumPOINT = -20 ;ENDDO
96
97 DO i=1,6
98 DO j=1,9 ; brick_list(:,:)%POLY(j)%FACE(i)%NumPOINT = -21 ;ENDDO
99 enddo;
100
101 brick_list(:,:)%Seg_add_LFT = -22
102 brick_list(:,:)%Seg_add_LLT = -23
103
104 DO i=1,5
105 DO j=1,6 ; brick_list(:,:)%Adjacent_Brick(j,i) = -24 ;ENDDO
106 ENDDO
107
108 DO i=1,6
109 DO j=1,9
110 DO k=1,5 ; brick_list(:,:)%POLY(j)%FACE(i)%Adjacent_Cell(k) = -25 ;enddo;
111 ENDDO
112 ENDDO
113
114 DO i=1,9
115 DO j=1,6 ; brick_list(:,:)%POLY(i)%FACE(j)%NAdjCell = -26 ;enddo;
116 ENDDO
117
118 DO i=1,8
119 DO j=1,9 ; brick_list(:,:)%POLY(j)%ListNodID(i) = -27 ;ENDDO
120 ENDDO
121
122 DO j=1,9 ; brick_list(:,:)%POLY(j)%IsMain = -28 ;ENDDO
123
124 DO j=1,8 ; brick_list(:,:)%NODE(j)%NodWasMain = -29 ;ENDDO
125
126 DO i=1,4
127 DO j=1,9 ; brick_list(:,:)%POLY(j)%WhereIsMain(i) = -30 ;ENDDO
128 ENDDO
129
130 DO i=1,8 ; brick_list(:,:)%NODE(i)%WhereWasMain = -31 ;ENDDO
131
132 DO i=1,3
133 DO j=1,5 ; brick_list(:,:)%MergeTarget(i,j) = -32 ;ENDDO
134 ENDDO
135
136 brick_list(:,:)%NTarget = -33
137
138 DO i=1,6 ; brick_list(:,:)%ClosedSurf(i) = -34 ;ENDDO
139
140 DO i=1,9 ; brick_list(:,:)%POLY(i)%ID_FREE_NODE = -35 ;ENDDO
141
142 DO i=1,6
143 DO j=1,9
144 DO k=1,5 ; brick_list(:,:)%POLY(j)%FACE(i)%Adjacent_FLUX(k) = -102.d00 ;ENDDO
145 ENDDO
146 ENDDO
147
148 DO i=1,6
149 DO j=1,9
150 DO k=1,5 ; brick_list(:,:)%POLY(j)%FACE(i)%Adjacent_upwFLUX(k)= -103.d00 ;ENDDO
151 ENDDO
152 ENDDO
153
154 DO i=1,9 ; brick_list(:,:)%POLY(i)%Adjacent_FLU1 = -104.d00 ;ENDDO
155
156 DO i=1,6
157 DO j=1,3 ; brick_list(:,:)%N(i,j) = -105.d00 ;ENDDO
158 ENDDO
159
160 DO j=1,9 ; brick_list(:,:)%POLY(j)%DDVOL = -106.d00 ;ENDDO
161
162 DO j=1,9 ; brick_list(:,:)%POLY(j)%PHI = -107.d00 ;ENDDO
163 DO j=1,9 ; brick_list(:,:)%POLY(j)%DPHI = -108.d00 ;ENDDO
164
165 DO i=1,9
166 DO j=1,4 ; brick_list(:,:)%POLY(i)%VFRACm(j) = -109.d00 ;ENDDO
167 ENDDO
168
169 brick_list(:,:)%Vnew_SCell = -110.d00
170 brick_list(:,:)%Vold_SCell = -111.d00
171 DO j=1,9 ; brick_list(:,:)%POLY(j)%Vnew = -112.d00 ;ENDDO
172 DO j=1,9 ; brick_list(:,:)%POLY(j)%Vold = -113.d00 ;ENDDO
173 DO j=1,9 ; brick_list(:,:)%POLY(j)%OLD_Vnew = -114.d00 ;ENDDO
174
175 DO i=1,3
176 brick_list(:,:)%SCellCenter(i) = -114.5d00
177 ENDDO
178
179 DO i=1,3
180 DO j=1,9 ; brick_list(:,:)%POLY(j)%CellCenter(i) = -115.d00 ;ENDDO
181 ENDDO
182
183 DO k=1,9
184 DO i=1,3
185 DO j=1,6 ; brick_list(:,:)%POLY(k)%FACE(j)%Center(i) = -115.5d00;ENDDO
186 ENDDO
187 ENDDO
188
189 DO j=1,6 ;
190 DO i=1,9 ; brick_list(:,:)%POLY(i)%FACE(j)%Surf = -116.d00 ;enddo;
191 ENDDO
192
193 DO k=1,9
194 DO j=1,3
195 DO i=1,6 ; brick_list(:,:)%POLY(k)%FACE(i)%Vel(j) = -117.d00 ;enddo;
196 ENDDO
197 ENDDO
198
199 DO j=1,6 ; brick_list(:,:)%Face_Brick(j) = -118.d00 ;ENDDO
200
201 brick_list(:,:)%UncutVol = -119.d00
202
203 DO k=1,9
204 DO i=1,6
205 DO j=1,3 ; brick_list(:,:)%POLY(k)%FACE(i)%F_FACE(j) = -121.d00 ;ENDDO
206 ENDDO
207 ENDDO
208
209 DO j=1,3 ; brick_list(:,:)%FCELL(j) = -122.d00 ;ENDDO
210
211 DO j=1,3 ; brick_list(:,:)%FEXT_CELL(j) = -123.d00 ;ENDDO
212
213 DO j=0,6 ; brick_list(:,:)%SIG(j) = -124.d00 ;ENDDO
214
215 DO i=0,2
216 DO j=1,9 ; brick_list(:,:)%POLY(j)%DVOL(i) = -125.d00 ;ENDDO
217 ENDDO
218 brick_list(:,:)%DVOL = -125.d00
219
220
221 DO j=1,2
222 DO i=1,6 ; brick_list(:,:)%Poly9woNodes(i,j) = -126.d00 ;ENDDO
223 ENDDO
224
225 brick_list(:,:)%ADJ_ELEMS%Num = -30
226 DO j=1,48
227 brick_list(:,:)%ADJ_ELEMS%Num_inv(j) = -31
228 brick_list(:,:)%ADJ_ELEMS%IV(j) = -32
229 brick_list(:,:)%ADJ_ELEMS%IB(j) = -33
230 brick_list(:,:)%ADJ_ELEMS%ICELL(j) = -34
231 brick_list(:,:)%ADJ_ELEMS%SecndFACE(j) = -35
232 ENDDO !next J
233
234 DO i=1,14
235 brick_list(:,:)%PCUT(i)%NP = -130.d00
236 brick_list(:,:)%PCUT(i)%N(1) = -131.d00
237 brick_list(:,:)%PCUT(i)%N(2) = -132.d00
238 brick_list(:,:)%PCUT(i)%N(3) = -133.d00
239 brick_list(:,:)%PCUT(i)%B(1) = -134.d00
240 brick_list(:,:)%PCUT(i)%B(2) = -135.d00
241 brick_list(:,:)%PCUT(i)%B(3) = -136.d00
242 brick_list(:,:)%PCUT(i)%SCUT(1) = -137.d00
243 brick_list(:,:)%PCUT(i)%SEFF(1) = -138.d00
244 DO j=1,8
245 brick_list(:,:)%PCUT(i)%P(1,j) = -139.d00
246 brick_list(:,:)%PCUT(i)%P(2,j) = -140.d00
247 brick_list(:,:)%PCUT(i)%P(3,j) = -141.d00
248 ENDDO !next J
249 ENDDO !next I
250
251 brick_list(:,:)%SecndList%VOL_Unmerged = -142.d00
252 brick_list(:,:)%SecndList%Num = -40
253 brick_list(:,:)%SecndList%NumSecndNodes = -41
254 DO j=1,24
255 brick_list(:,:)%SecndList%FM(j) = -42
256 brick_list(:,:)%SecndList%FV(j) = -43
257 brick_list(:,:)%SecndList%IV(j) = -44
258 brick_list(:,:)%SecndList%IBV(j) = -45
259 brick_list(:,:)%SecndList%ICELLv(j) = -46
260 brick_list(:,:)%SecndList%VOL(j) = -1.0d00
261 brick_list(:,:)%SecndList%NumNOD_Cell(j) = -47
262 do k=1,8
263 brick_list(:,:)%SecndList%ListNodID(j,k) = -48
264 enddo
265 ENDDO !next J
266
267 DO j=1,12
268 brick_list(:,:)%Edge(j)%NODE(1) = -49
269 brick_list(:,:)%Edge(j)%NODE(2) = -50
270 brick_list(:,:)%Edge(j)%NBCUT = 0
271 brick_list(:,:)%Edge(j)%CUTSHELL(1) = -59
272 brick_list(:,:)%Edge(j)%CUTSHELL(2) = -59
273 brick_list(:,:)%Edge(j)%CUTCOOR(1) = -60.0d00
274 brick_list(:,:)%Edge(j)%CUTCOOR(2) = -60.0d00
275 brick_list(:,:)%Edge(j)%CUTPOINT(1,1) = -61.0d00
276 brick_list(:,:)%Edge(j)%CUTPOINT(2,1) = -61.0d00
277 brick_list(:,:)%Edge(j)%CUTPOINT(3,1) = -61.0d00
278 brick_list(:,:)%Edge(j)%CUTPOINT(1,2) = -61.0d00
279 brick_list(:,:)%Edge(j)%CUTPOINT(2,2) = -61.0d00
280 brick_list(:,:)%Edge(j)%CUTPOINT(3,2) = -61.0d00
281 brick_list(:,:)%Edge(j)%CUTVEL(1,1) = -62.0d00
282 brick_list(:,:)%Edge(j)%CUTVEL(2,1) = -62.0d00
283 brick_list(:,:)%Edge(j)%CUTVEL(3,1) = -62.0d00
284 brick_list(:,:)%Edge(j)%CUTVEL(1,2) = -62.0d00
285 brick_list(:,:)%Edge(j)%CUTVEL(2,2) = -62.0d00
286 brick_list(:,:)%Edge(j)%CUTVEL(3,2) = -62.0d00
287 brick_list(:,:)%Edge(j)%VECTOR(1) = -63.0d00
288 brick_list(:,:)%Edge(j)%VECTOR(2) = -63.0d00
289 brick_list(:,:)%Edge(j)%VECTOR(3) = -63.0d00
290 brick_list(:,:)%Edge(j)%LEN = -64.0d00
291 ENDDO
292
293 brick_list(:,:)%SECTYPE(1) = 'init-1 '
294 brick_list(:,:)%SECTYPE(2) = 'init-2 '
295 brick_list(:,:)%SECTYPE(3) = 'init-3 '
296 brick_list(:,:)%SECTYPE(4) = 'init-4 '
297 brick_list(:,:)%SECTYPE(5) = 'init-5 '
298 brick_list(:,:)%SECTYPE(6) = 'init-6 '
299 brick_list(:,:)%SECTYPE(7) = 'init-7 '
300 brick_list(:,:)%SECTYPE(8) = 'init-8 '
301 ENDIF !(.NOT.ALLOCATED(BRICK_LIST))
302
303 brick_list(:,:)%bakMAT%rho = -65d00
304 brick_list(:,:)%bakMAT%rhoE = -66d00
305 brick_list(:,:)%bakMAT%rhoU(1) = -67d00
306 brick_list(:,:)%bakMAT%rhoU(2) = -67d00
307 brick_list(:,:)%bakMAT%rhoU(3) = -67d00
308 brick_list(:,:)%bakMAT%sig(1) = -68d00
309 brick_list(:,:)%bakMAT%sig(2) = -68d00
310 brick_list(:,:)%bakMAT%sig(3) = -68d00
311 brick_list(:,:)%bakMAT%sig(4) = -68d00
312 brick_list(:,:)%bakMAT%sig(5) = -68d00
313 brick_list(:,:)%bakMAT%sig(6) = -68d00
314 brick_list(:,:)%bakMAT%ssp = -69d00
315
316 ALLOCATE(old_secndlist(int22,i22len) )
317 DO nin=1,int22
318 !=====%OLD_SecndList(NIN,1:I22LEN)
319 old_secndlist(nin,1:i22len)%VOL_Unmerged = zero
320 old_secndlist(nin,1:i22len)%Num = 0
321 old_secndlist(nin,1:i22len)%NumSecndNodes = 0
322 DO j=1,24
323 old_secndlist(nin,1:i22len)%FM(j) = 0
324 old_secndlist(nin,1:i22len)%FV(j) = 0
325 old_secndlist(nin,1:i22len)%IV(j) = 0
326 old_secndlist(nin,1:i22len)%IBV(j) = 0
327 old_secndlist(nin,1:i22len)%ICELLv(j) = 0
328 old_secndlist(nin,1:i22len)%VOL(j) = zero
329 old_secndlist(nin,1:i22len)%SURF_v(j) = zero
330 old_secndlist(nin,1:i22len)%NumNOD_Cell(j) = 0
331 DO k=1,8
332 old_secndlist(nin,1:i22len)%ListNodID(j,k) = 0
333 ENDDO
334 ENDDO !next J
335 enddo!next NIN
336
337
338
339
340 DO nin=1,int22
341 !=====general data
342 CALL write_i_c( nbold(nin) , 1)
343 CALL write_i_c( list_b_old(nin,1:i22len) , i22len)
344 !=====cartesian cell data - INTEGER
345 CALL write_i_c( brick_list(nin,1:i22len)%ID , i22len)
346 CALL write_i_c( brick_list(nin,1:i22len)%NG , i22len)
347 CALL write_i_c( brick_list(nin,1:i22len)%IDLOC , i22len)
348 CALL write_i_c( brick_list(nin,1:i22len)%ICODE , i22len)
349 CALL write_i_c( brick_list(nin,1:i22len)%OLD_ICODE , i22len)
350 CALL write_i_c( brick_list(nin,1:i22len)%IDBLE , i22len)
351 CALL write_i_c( brick_list(nin,1:i22len)%NBITS , i22len)
352 CALL write_i_c( brick_list(nin,1:i22len)%NPQTS , i22len)
353 CALL write_i_c( brick_list(nin,1:i22len)%NBCUT , i22len)
354 CALL write_i_c( brick_list(nin,1:i22len)%MainID , i22len)
355 CALL write_i_c( brick_list(nin,1:i22len)%WasCut , i22len)
356 CALL write_i_c( brick_list(nin,1:i22len)%NewInBuffer , i22len)
357 CALL write_i_c( brick_list(nin,1:i22len)%OldMainStrongNode , i22len)
358 CALL write_i_c( brick_list(nin,1:i22len)%MLW , i22len)
359 CALL write_i_c( brick_list(nin,1:i22len)%ITASK , i22len)
360 DO j=1,8 ; CALL write_i_c( brick_list(nin,1:i22len)%SecID_Cell(j) , i22len) ; ENDDO
361 CALL write_i_c( brick_list(nin,1:i22len)%Seg_add_LFT , i22len)
362 CALL write_i_c( brick_list(nin,1:i22len)%Seg_add_LLT , i22len)
363 DO i=1,5
364 DO j=1,6 ; CALL write_i_c( brick_list(nin,1:i22len)%Adjacent_Brick(j,i) , i22len) ; ENDDO
365 ENDDO
366 DO i=1,3
367 DO j=1,5 ; CALL write_i_c ( brick_list(nin,1:i22len)%MergeTarget(i,j) , i22len) ; ENDDO
368 ENDDO
369 CALL write_i_c( brick_list(nin,1:i22len)%NTarget , i22len)
370 DO i=1,6 ; CALL write_i_c ( brick_list(nin,1:i22len)%ClosedSurf(i) , i22len) ; ENDDO
371 CALL write_i_c( brick_list(nin,1:i22len)%IsMErgeTarget , i22len)
372 !=====cartesian cell data - REAL
373 DO i=1,6
374 DO j=1,3 ; CALL write_db ( brick_list(nin,1:i22len)%N(i,j) , i22len) ; ENDDO
375 ENDDO
376 CALL write_db ( brick_list(nin,1:i22len)%Vnew_SCell , i22len)
377 CALL write_db ( brick_list(nin,1:i22len)%Vold_SCell , i22len)
378 DO i=1,3 ; CALL write_db ( brick_list(nin,1:i22len)%SCellCenter(i) , i22len) ; ENDDO
379 DO j=1,6 ; CALL write_db ( brick_list(nin,1:i22len)%Face_Brick(j) , i22len) ; ENDDO
380 CALL write_db ( brick_list(nin,1:i22len)%UncutVol , i22len)
381 DO j=1,3 ; CALL write_db ( brick_list(nin,1:i22len)%FCELL(j) , i22len) ; ENDDO
382 DO j=1,3 ; CALL write_db ( brick_list(nin,1:i22len)%FEXT_CELL(j) , i22len) ; ENDDO
383 DO j=0,6 ; CALL write_db ( brick_list(nin,1:i22len)%SIG(j) , i22len) ; ENDDO
384 CALL write_db ( brick_list(nin,1:i22len)%RHOC , i22len)
385 CALL write_db ( brick_list(nin,1:i22len)%MACH , i22len)
386 CALL write_db ( brick_list(nin,1:i22len)%DVOL , i22len)
387 DO j=1,2
388 DO i=1,6 ; CALL write_db ( brick_list(nin,1:i22len)%Poly9woNodes(i,j) , i22len) ; ENDDO
389 ENDDO
390 !=====%POLY(1:9) - INTEGER
391 DO j=1,9 ; CALL write_i_c( brick_list(nin,1:i22len)%POLY(j)%NumNOD , i22len) ; ENDDO
392 DO j=1,9 ; CALL write_i_c( brick_list(nin,1:i22len)%POLY(j)%NumPOINT , i22len) ; ENDDO
393 DO i=1,8
394 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%ListNodID(i) , i22len) ; ENDDO
395 ENDDO
396 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%IsMain , i22len) ; ENDDO
397 DO i=1,4
398 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%WhereIsMain(i) , i22len) ; ENDDO
399 ENDDO
400 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%ID_FREE_NODE , i22len) ; ENDDO
401 !=====%POLY(1:9) - REAL
402 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%Adjacent_FLU1 , i22len) ; ENDDO
403 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%DDVOL , i22len) ; ENDDO
404 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%PHI , i22len) ; ENDDO
405 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%dPHI , i22len) ; ENDDO
406 DO j=1,9
407 DO i=1,4 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%VFRACm(i) , i22len) ; ENDDO
408 ENDDO
409 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%Vnew , i22len) ; ENDDO
410 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%Vold , i22len) ; ENDDO
411 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%OLD_Vnew , i22len) ; ENDDO
412 DO i=1,3
413 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%CellCenter(i) , i22len) ; ENDDO
414 ENDDO
415 DO i=0,2
416 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%DVOL(i) , i22len) ; ENDDO
417 ENDDO
418 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%cutFace , i22len) ; ENDDO
419 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%DM , i22len) ; ENDDO
420 !=====%POLY(1:9)%FACE0
421 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%FACE0%NumPOINT , i22len) ; ENDDO
422 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%FACE0%NAdjCell , i22len) ; ENDDO
423 DO i=1,3
424 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE0%Center(i) , i22len) ; ENDDO
425 ENDDO
426 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE0%Surf , i22len) ; ENDDO
427 DO i=1,3
428 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE0%F_FACE(i) , i22len) ; ENDDO
429 ENDDO
430 DO i=1,9
431 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE0%U_N(i) , i22len) ; ENDDO
432 ENDDO
433 !=====%POLY(1:9)%FACE(1:6)
434 DO k=1,6
435 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%NumPOINT , i22len) ; ENDDO
436 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%NAdjCell , i22len) ; ENDDO
437 DO i=1,5
438 DO j=1,9 ; CALL write_i_c ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Adjacent_Cell(i) , i22len) ; ENDDO
439 ENDDO
440 DO i=1,5
441 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Adjacent_FLUX(i) , i22len) ; ENDDO
442 ENDDO
443 DO i=1,5
444 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Adjacent_upwFLUX(i), i22len) ; ENDDO
445 ENDDO
446 DO i=1,3
447 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Center(i) , i22len) ; ENDDO
448 ENDDO
449 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Surf , i22len) ; ENDDO
450 DO i=1,3
451 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%Vel(i) , i22len) ; ENDDO
452 ENDDO
453 DO i=1,3
454 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%F_FACE(i) , i22len) ; ENDDO
455 ENDDO
456 DO j=1,9 ; CALL write_db ( brick_list(nin,1:i22len)%POLY(j)%FACE(k)%U_N , i22len) ; ENDDO
457 enddo!next K
458 !=====%NODE
459 DO j=1,8 ; CALL write_i_c( brick_list(nin,1:i22len)%NODE(j)%OLD_WhichCell , i22len) ; ENDDO
460 DO j=1,8 ; CALL write_i_c( brick_list(nin,1:i22len)%NODE(j)%WhichCell , i22len) ; ENDDO
461 DO j=1,8 ; CALL write_i_c( brick_list(nin,1:i22len)%NODE(j)%NodWasMain , i22len) ; ENDDO
462 DO i=1,8 ; CALL write_i_c( brick_list(nin,1:i22len)%NODE(i)%WhereWasMain , i22len) ; ENDDO
463 !=====%ADJ_ELEMS
464 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%Num , i22len)
465 DO j=1,48
466 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%Num_inv(j) , i22len)
467 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%IV(j) , i22len)
468 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%IB(j) , i22len)
469 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%ICELL(j) , i22len)
470 CALL write_i_c ( brick_list(nin,1:i22len)%ADJ_ELEMS%SecndFACE(j) , i22len)
471 ENDDO !next J
472 !=====%PCUT
473 DO i=1,16
474 CALL write_i_c ( brick_list(nin,1:i22len)%PCUT(i)%NP , i22len)
475 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%N(1) , i22len)
476 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%N(2) , i22len)
477 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%N(3) , i22len)
478 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%B(1) , i22len)
479 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%B(2) , i22len)
480 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%B(3) , i22len)
481 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%SCUT(1) , i22len)
482 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%SEFF(1) , i22len)
483 DO j=1,8
484 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%P(1,j) , i22len)
485 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%P(2,j) , i22len)
486 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%P(3,j) , i22len)
487 ENDDO !next J
488 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%Vel(1) , i22len)
489 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%Vel(2) , i22len)
490 CALL write_db ( brick_list(nin,1:i22len)%PCUT(i)%Vel(3) , i22len)
491 ENDDO !next J
492 !=====%SecndList
493 CALL write_db ( brick_list(nin,1:i22len)%SecndList%VOL_Unmerged , i22len)
494 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%Num , i22len)
495 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%NumSecndNodes , i22len)
496 DO j=1,24
497 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%FM(j) , i22len)
498 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%FV(j) , i22len)
499 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%IV(j) , i22len)
500 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%IBV(j) , i22len)
501 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%ICELLv(j) , i22len)
502 CALL write_db ( brick_list(nin,1:i22len)%SecndList%VOL(j) , i22len)
503 CALL write_db ( brick_list(nin,1:i22len)%SecndList%SURF_v(j) , i22len)
504 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%NumNOD_Cell(j) , i22len)
505 DO k=1,8
506 CALL write_i_c ( brick_list(nin,1:i22len)%SecndList%ListNodID(j,k) , i22len)
507 ENDDO
508 ENDDO !next J
509 !=====%Edge
510 DO j=1,12
511 CALL write_i_c ( brick_list(nin,1:i22len)%Edge(j)%NODE(1) , i22len)
512 CALL write_i_c ( brick_list(nin,1:i22len)%Edge(j)%NODE(2) , i22len)
513 CALL write_i_c ( brick_list(nin,1:i22len)%Edge(j)%NBCUT , i22len)
514 CALL write_i_c ( brick_list(nin,1:i22len)%Edge(j)%CUTSHELL(1) , i22len)
515 CALL write_i_c ( brick_list(nin,1:i22len)%Edge(j)%CUTSHELL(2) , i22len)
516 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTCOOR(1) , i22len)
517 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTCOOR(2) , i22len)
518 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(1,1) , i22len)
519 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(2,1) , i22len)
520 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(3,1) , i22len)
521 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(1,2) , i22len)
522 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(2,2) , i22len)
523 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTPOINT(3,2) , i22len)
524 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(1,1) , i22len)
525 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(2,1) , i22len)
526 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(3,1) , i22len)
527 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(1,2) , i22len)
528 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(2,2) , i22len)
529 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%CUTVEL(3,2) , i22len)
530 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%VECTOR(1) , i22len)
531 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%VECTOR(2) , i22len)
532 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%VECTOR(3) , i22len)
533 CALL write_db ( brick_list(nin,1:i22len)%Edge(j)%LEN , i22len)
534 ENDDO
535 !=====%SECTYPE
536 DO j=1,8
537 DO i=1,i22len
538 string = brick_list(nin,i)%SECTYPE(j)
539 do k=1,14
540 CALL write_i_c( ichar(string(k:k)), 1)
541 enddo
542 ENDDO
543 ENDDO
544 !=====%bakMAT
545 CALL write_db ( brick_list(nin,1:i22len)%bakMAT%rho , i22len)
546 CALL write_db ( brick_list(nin,1:i22len)%bakMAT%rhoE , i22len)
547 DO i=1,3 ; CALL write_db ( brick_list(nin,1:i22len)%bakMAT%rhoU(i) , i22len) ; ENDDO
548 DO i=1,6 ; CALL write_db ( brick_list(nin,1:i22len)%bakMAT%sig(i) , i22len) ; ENDDO
549 CALL write_db ( brick_list(nin,1:i22len)%bakMAT%ssp , i22len)
550 isiz = max(i22law37, i22law51)
551 bid = zero
552 IF(isiz/=0)THEN
553 DO k=1,i22len
554 DO i=1,isiz
555 !CALL WRITE_DB ( BRICK_LIST(NIN,K)%bakMAT%UVAR(I) , I22LEN)
556 CALL write_db ( bid , 1)
557 ENDDO
558 ENDDO
559 ENDIF
560
561
562 !=====%OLD_SecndList(NIN,1:I22LEN)
563 CALL write_db (old_secndlist(nin,1:i22len)%VOL_Unmerged , i22len)
564 CALL write_i_c (old_secndlist(nin,1:i22len)%Num , i22len)
565 CALL write_i_c (old_secndlist(nin,1:i22len)%NumSecndNodes , i22len)
566 DO j=1,24
567 CALL write_i_c ( old_secndlist(nin,1:i22len)%FM(j) , i22len)
568 CALL write_i_c ( old_secndlist(nin,1:i22len)%FV(j) , i22len)
569 CALL write_i_c ( old_secndlist(nin,1:i22len)%IV(j) , i22len)
570 CALL write_i_c ( old_secndlist(nin,1:i22len)%IBV(j) , i22len)
571 CALL write_i_c ( old_secndlist(nin,1:i22len)%ICELLv(j) , i22len)
572 CALL write_db ( old_secndlist(nin,1:i22len)%VOL(j) , i22len)
573 CALL write_db ( old_secndlist(nin,1:i22len)%SURF_v(j) , i22len)
574 CALL write_i_c ( old_secndlist(nin,1:i22len)%NumNOD_Cell(j) , i22len)
575 DO k=1,8
576 CALL write_i_c ( old_secndlist(nin,1:i22len)%ListNodID(j,k) , i22len)
577 ENDDO
578 ENDDO !next J
579 DEALLOCATE(old_secndlist)
580
581
582 ENDDO !next NIN
583
584C-----------------------------------------------
585 1000 CONTINUE
586
587 IF (ierr1/=0) THEN
588 WRITE(iout,*)' ** ERROR IN MEMORY ALLOCATION (INTER22)'
589 WRITE(istdo,*)' ** ERROR IN MEMORY ALLOCATION (INTER22)'
590 CALL arret(2)
591 END IF
592C-----------------------------------------------
593 END SUBROUTINE
594
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
i22bufbric_mod
Definition cut-cell-search_mod.F:95
i22bufbric_mod::brick_list
type(brick_entity), dimension(:,:), allocatable, target brick_list
Definition cut-cell-search_mod.F:105
i22tri_mod
Definition cut-cell-search_mod.F:212
i22tri_mod::nb
integer nb
Definition cut-cell-search_mod.F:487
w_bufbric_22
subroutine w_bufbric_22()
Definition w_bufbric_22.F:32
arret
subroutine arret(nn)
Definition arret.F:87
write_db
subroutine write_db(a, n)
Definition write_db.F:140
write_i_c
void write_i_c(int *w, int *len)
Definition write_routtines.c:645