OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
state_admesh.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!|| state_admesh ../starter/source/model/remesh/state_admesh.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
32!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
33!|| nintlst2 ../starter/source/system/nintrr.F
34!||--- uses -----------------------------------------------------
35!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE state_admesh(IPART ,IPARTC,IPARTTG,IXC ,IXTG ,
40 . SH4TREE,SH3TREE,SH4TRIM,SH3TRIM,LSUBMODEL)
41C----------------------------------------------------------
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE message_mod
46 USE hm_option_read_mod
47 USE submodel_mod
48 USE names_and_titles_mod , ONLY : nchartitle, ncharkey
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C C o m m o n B l o c k s
55C-----------------------------------------------
56#include "param_c.inc"
57#include "com04_c.inc"
58#include "scr17_c.inc"
59#include "remesh_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER IPART(LIPART1,*), IPARTC(*), IPARTTG(*),
64 . ixc(nixc,*), ixtg(nixtg,*),
65 . sh4tree(ksh4tree,*), sh3tree(ksh3tree,*),
66 . sh4trim(*), sh3trim(*)
67 TYPE(submodel_data) LSUBMODEL(*)
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER ID,ID1,ID2,ID3,ID4,II,I1,I2,I3,I4,NLIST,N,LEVEL,NN,
72 . itrim,i
73 INTEGER IERROR, NINTLST2,ERRORADJ,NSHELL,NSH3N
74 INTEGER, DIMENSION(:),ALLOCATABLE :: LIST,INDEXL
75 INTEGER IX1(MAX(NUMELC,NUMELTG)),
76 . ix2(max(numelc,numeltg)),
77 . index(2*max(numelc,numeltg))
78 CHARACTER MESS*40
79 CHARACTER(LEN=NCHARTITLE) :: TITR
80 CHARACTER(LEN=NCHARKEY) :: KEY,KEY2
81 LOGICAL IS_AVAILABLE
82C-----------------------------------------------
83 DATA mess /'ADAPTIVE MESHING STATE DEFINITION '/
84C-----------------------------------------------
85 nlist =0
86
87C-------------------------------------------------------------
88C READING /ADMESH/STATE/SHELL : 4-NODE SHELLS => Counting
89C------------------------------------------------------------
90C
91 CALL hm_option_start('/ADMESH/STATE')
92C
93 DO n =1,nadmeshstat
94 titr = ''
95C
96 CALL hm_option_read_key(lsubmodel,
97 . option_titr = titr,
98 . keyword2 = key,
99 . keyword3 = key2)
100
101 IF(key2(1:len_trim(key2))=='SHELL')THEN
102
103C
104 is_available = .false.
105C
106C--------* EXTRACT DATAS (INTEGER VALUES) *------
107C
108 CALL hm_get_intv('NSHELL',nshell,is_available,lsubmodel)
109 nlist = nlist + nshell
110C
111 ENDIF
112
113 ENDDO
114C------
115 ALLOCATE(list(5*nlist),stat=ierror)
116 IF(ierror/=0) CALL ancmsg(msgid=268,anmode=aninfo,
117 . msgtype=msgerror,
118 . c1='LIST')
119 ALLOCATE(indexl(10*nlist),stat=ierror)
120 IF(ierror/=0) CALL ancmsg(msgid=268,anmode=aninfo,
121 . msgtype=msgerror,
122 . c1='INDEXL')
123
124 nlist=0
125
126C-------------------------------------------------------------------------------
127C READING /ADMESH/STAT/SHELL : 4-NODE SHELLS => Storing and Tri( local Ids)
128C-------------------------------------------------------------------------------
129C
130 CALL hm_option_start('/ADMESH/STATE')
131
132 DO n =1,nadmeshstat
133 titr = ''
134
135 CALL hm_option_read_key(lsubmodel,
136 . option_titr = titr,
137 . keyword2 = key,
138 . keyword3 = key2)
139
140 IF(key2(1:len_trim(key2))=='SHELL')THEN
141
142C
143 is_available = .false.
144C
145C--------* EXTRACT DATAS (INTEGER VALUES) *------
146C
147 CALL hm_get_intv('NSHELL',nshell,is_available,lsubmodel)
148
149 DO i=1,nshell
150
151 CALL hm_get_int_array_index('shell_ID',id,i,is_available,lsubmodel)
152 CALL hm_get_int_array_index('shell_ID1',id1,i,is_available,lsubmodel)
153 CALL hm_get_int_array_index('shell_ID2',id2,i,is_available,lsubmodel)
154 CALL hm_get_int_array_index('shell_ID3',id3,i,is_available,lsubmodel)
155 CALL hm_get_int_array_index('shell_ID4',id4,i,is_available,lsubmodel)
156
157 nlist=nlist+1
158 list(nlist)=id
159 nlist=nlist+1
160 list(nlist)=id1
161 nlist=nlist+1
162 list(nlist)=id2
163 nlist=nlist+1
164 list(nlist)=id3
165 nlist=nlist+1
166 list(nlist)=id4
167C
168 ENDDO
169 ENDIF
170 ENDDO
171
172 nn=nintlst2(list,nlist,indexl,ixc,nixc,numelc,
173 . mess,ix1,ix2,index,0)
174
175 nlist=0
176
177C---------------------------------------------------------------------
178C READING /ADMESH/STAT/SHELL : 4-NODE SHELLS => Storing in SH4TREE
179C---------------------------------------------------------------------
180
181C
182 CALL hm_option_start('/ADMESH/STATE')
183
184 DO n =1,nadmeshstat
185 titr = ''
186
187 CALL hm_option_read_key(lsubmodel,
188 . option_titr = titr,
189 . keyword2 = key,
190 . keyword3 = key2)
191
192 IF(key2(1:len_trim(key2))=='SHELL')THEN
193
194C
195 is_available = .false.
196C
197C--------* EXTRACT DATAS (INTEGER VALUES) *------
198C
199 CALL hm_get_intv('NSHELL',nshell,is_available,lsubmodel)
200
201 DO i=1,nshell
202
203
204 CALL hm_get_int_array_index('shell_ID',id,i,is_available,lsubmodel)
205 CALL hm_get_int_array_index('shell_ID1',id1,i,is_available,lsubmodel)
206 CALL hm_get_int_array_index('shell_ID2',id2,i,is_available,lsubmodel)
207 CALL hm_get_int_array_index('shell_ID3',id3,i,is_available,lsubmodel)
208 CALL hm_get_int_array_index('shell_ID4',id4,i,is_available,lsubmodel)
209 CALL hm_get_int_array_index('Actlev',level,i,is_available,lsubmodel)
210 CALL hm_get_int_array_index('IMapping',itrim,i,is_available,lsubmodel)
211
212 nlist=nlist+1
213 ii=list(nlist)
214 nlist=nlist+1
215 i1=list(nlist)
216 nlist=nlist+1
217 i2=list(nlist)
218 nlist=nlist+1
219 i3=list(nlist)
220 nlist=nlist+1
221 i4=list(nlist)
222 IF(i1+i2+i3+i4 /=0 .AND.
223 . (i2-i1/=1 .OR. i3-i1 /= 2 .OR. i4-i1 /= 3))THEN
224 CALL ancmsg(msgid=654,
225 . msgtype=msgerror,
226 . anmode=aninfo,
227 . i1=id1,
228 . i2=id2,
229 . i3=id3,
230 . i4=id4,
231 . i5=id)
232 END IF
233C
234 IF(i1+i2+i3+i4 /=0) THEN
235 erroradj =0
236 IF(ixc(2,ii) /= ixc(2,i1).OR.ixc(3,ii) /= ixc(3,i2)
237 . .OR.ixc(4,ii) /= ixc(4,i3).OR.ixc(5,ii) /= ixc(5,i4)) THEN
238 erroradj =1
239 ELSEIF(ixc(4,i1) /= ixc(5,i2).OR.ixc(5,i2) /= ixc(2,i3)
240 . .OR.ixc(2,i3) /= ixc(3,i4).OR.ixc(4,i1) /= ixc(3,i4)) THEN
241 erroradj =1
242 ELSEIF(ixc(3,i1) /= ixc(2,i2).OR.ixc(4,i2) /= ixc(3,i3)
243 . .OR.ixc(5,i3) /= ixc(4,i4).OR.ixc(5,i1) /= ixc(2,i4)) THEN
244 erroradj =1
245 ENDIF
246 ENDIF
247 IF(erroradj ==1.AND.abs(level)<levelmax) THEN
248 CALL ancmsg(msgid=1023,
249 . msgtype=msgerror,
250 . anmode=aninfo,
251 . i1=id1,
252 . i2=id2,
253 . i3=id3,
254 . i4=id4,
255 . i5=id)
256 END IF
257C
258 IF(level<-levelmax-1.OR.level>levelmax)THEN
259 CALL ancmsg(msgid=656,
260 . msgtype=msgerror,
261 . anmode=aninfo,
262 . i1=id)
263 END IF
264 sh4tree(2,ii)=i1
265 sh4tree(3,ii)=level
266 IF(i1/=0)THEN
267 sh4tree(1,i1)=ii
268 sh4tree(1,i2)=ii
269 sh4tree(1,i3)=ii
270 sh4tree(1,i4)=ii
271 END IF
272 sh4trim(ii)=itrim
273C
274 ENDDO
275 ENDIF
276 ENDDO
277
278 DEALLOCATE(list)
279 DEALLOCATE(indexl)
280
281C-----------------------------------------------
282C 3-NODE SHELLS
283C-----------------------------------------------
284 nlist=0
285
286C-----------------------------------------------------------
287C READING /ADMESH/STAT/SH3N : 3-NODE SHELLS => Counting
288C----------------------------------------------------------
289
290 CALL hm_option_start('/ADMESH/STATE')
291
292 DO n =1,nadmeshstat
293 titr = ''
294
295 CALL hm_option_read_key(lsubmodel,
296 . option_titr = titr,
297 . keyword2 = key,
298 . keyword3 = key2)
299
300 IF(key2(1:len_trim(key2))=='SH3N')THEN
301
302C
303 is_available = .false.
304C
305C--------* EXTRACT DATAS (INTEGER VALUES) *------
306C
307 CALL hm_get_intv('NSH3N',nsh3n,is_available,lsubmodel)
308 nlist = nlist + nsh3n
309C
310
311 ENDIF
312
313 ENDDO
314
315C------
316 ALLOCATE(list(5*nlist),stat=ierror)
317 IF(ierror/=0) CALL ancmsg(msgid=268,anmode=aninfo,
318 . msgtype=msgerror,
319 . c1='LIST')
320 ALLOCATE(indexl(10*nlist),stat=ierror)
321 IF(ierror/=0) CALL ancmsg(msgid=268,anmode=aninfo,
322 . msgtype=msgerror,
323 . c1='INDEXL')
324
325C---------
326
327 nlist=0
328
329C------------------------------------------------------------------------------
330C READING /ADMESH/STAT/SH3N : 3-NODE SHELLS => Storing and Tri( local Ids)
331C------------------------------------------------------------------------------
332
333 CALL hm_option_start('/ADMESH/STATE')
334
335 DO n =1,nadmeshstat
336 titr = ''
337
338 CALL hm_option_read_key(lsubmodel,
339 . option_titr = titr,
340 . keyword2 = key,
341 . keyword3 = key2)
342
343 IF(key2(1:len_trim(key2))=='SH3N')THEN
344C
345 is_available = .false.
346C
347C--------* EXTRACT DATAS (INTEGER VALUES) *------
348C
349 CALL hm_get_intv('NSH3N',nsh3n,is_available,lsubmodel)
350
351 DO i=1,nsh3n
352
353 CALL hm_get_int_array_index('sh3n_ID',id,i,is_available,lsubmodel)
354 CALL hm_get_int_array_index('sh3n_ID1',id1,i,is_available,lsubmodel)
355 CALL hm_get_int_array_index('sh3n_ID2',id2,i,is_available,lsubmodel)
356 CALL hm_get_int_array_index('sh3n_ID3',id3,i,is_available,lsubmodel)
357 CALL hm_get_int_array_index('sh3n_id4',ID4,I,IS_AVAILABLE,LSUBMODEL)
358
359 NLIST=NLIST+1
360 LIST(NLIST)=ID
361 NLIST=NLIST+1
362 LIST(NLIST)=ID1
363 NLIST=NLIST+1
364 LIST(NLIST)=ID2
365 NLIST=NLIST+1
366 LIST(NLIST)=ID3
367 NLIST=NLIST+1
368 LIST(NLIST)=ID4
369 ENDDO
370C
371
372 ENDIF
373
374 ENDDO
375
376 NN=NINTLST2(LIST,NLIST,INDEXL,IXTG,NIXTG,NUMELTG,
377 . MESS,IX1,IX2,INDEX,0)
378
379 NLIST=0
380C---------------------------------------------------------------------
381C READING /ADMESH/STAT/SH3N : 3-NODE SHELLS => Storing in SH3TREE
382C---------------------------------------------------------------------
383
384 CALL HM_OPTION_START('/admesh/state')
385
386 DO N =1,NADMESHSTAT
387 TITR = ''
388C
389C--------* EXTRACT DATAS OF /ADMESH/... LINE *------
390C
391 CALL HM_OPTION_READ_KEY(LSUBMODEL,
392 . OPTION_TITR = TITR,
393 . KEYWORD2 = KEY,
394 . KEYWORD3 = KEY2)
395
396 IF(KEY2(1:LEN_TRIM(KEY2))=='sh3n')THEN
397
398C
399 IS_AVAILABLE = .FALSE.
400C
401C--------* EXTRACT DATAS (INTEGER VALUES) *------
402C
403 CALL HM_GET_INTV('nsh3n',NSH3N,IS_AVAILABLE,LSUBMODEL)
404
405 DO I=1,NSH3N
406
407 CALL HM_GET_INT_ARRAY_INDEX('sh3n_id',ID,I,IS_AVAILABLE,LSUBMODEL)
408 CALL HM_GET_INT_ARRAY_INDEX('sh3n_id1',ID1,I,IS_AVAILABLE,LSUBMODEL)
409 CALL HM_GET_INT_ARRAY_INDEX('sh3n_id2',ID2,I,IS_AVAILABLE,LSUBMODEL)
410 CALL HM_GET_INT_ARRAY_INDEX('sh3n_id3',ID3,I,IS_AVAILABLE,LSUBMODEL)
411 CALL HM_GET_INT_ARRAY_INDEX('sh3n_id4',ID4,I,IS_AVAILABLE,LSUBMODEL)
412 CALL HM_GET_INT_ARRAY_INDEX('actlev',LEVEL,I,IS_AVAILABLE,LSUBMODEL)
413 CALL HM_GET_INT_ARRAY_INDEX('imapping',ITRIM,I,IS_AVAILABLE,LSUBMODEL)
414
415
416 NLIST=NLIST+1
417 II=LIST(NLIST)
418 NLIST=NLIST+1
419 I1=LIST(NLIST)
420 NLIST=NLIST+1
421 I2=LIST(NLIST)
422 NLIST=NLIST+1
423 I3=LIST(NLIST)
424 NLIST=NLIST+1
425 I4=LIST(NLIST)
426C
427.AND. IF(I1+I2+I3+I4 /=0
428.OR..OR. . (I2-I1/=1 I3-I1 /= 2 I4-I1 /= 3))THEN
429 CALL ANCMSG(MSGID=655,
430 . MSGTYPE=MSGERROR,
431 . ANMODE=ANINFO,
432 . I1=ID1,
433 . I2=ID2,
434 . I3=ID3,
435 . I4=ID4,
436 . I5=ID)
437 END IF
438C
439 IF(I1+I2+I3+I4 /=0) THEN
440 ERRORADJ =0
441.OR. IF(IXTG(2,II) /= IXTG(2,I1)IXTG(3,II) /= IXTG(3,I2)
442.OR. . IXTG(4,II) /= IXTG(4,I3)) THEN
443 ERRORADJ =1
444.OR. ELSEIF(IXTG(3,I1) /= IXTG(2,I2)IXTG(4,I2) /= IXTG(3,I3)
445.OR. . IXTG(2,I3) /= IXTG(3,I4)) THEN
446 ERRORADJ =1
447.OR. ELSEIF(IXTG(4,I1) /= IXTG(2,I3)IXTG(4,I1) /= IXTG(3,I4)
448.OR..OR..OR. . IXTG(4,I2) /= IXTG(2,I4)IXTG(4,I4) /= IXTG(2,I2)
449 . IXTG(4,I4) /= IXTG(3,I1)) THEN
450 ERRORADJ =1
451 ENDIF
452 ENDIF
453.AND. IF(ERRORADJ ==1ABS(LEVEL)<LEVELMAX) THEN
454 CALL ANCMSG(MSGID=1023,
455 . MSGTYPE=MSGERROR,
456 . ANMODE=ANINFO,
457 . I1=ID1,
458 . I2=ID2,
459 . I3=ID3,
460 . I4=ID4,
461 . I5=ID)
462 END IF
463C
464.OR. IF(LEVEL<-LEVELMAX-1LEVEL>LEVELMAX)THEN
465 CALL ANCMSG(MSGID=657,
466 . MSGTYPE=MSGERROR,
467 . ANMODE=ANINFO,
468 . I1=ID)
469 END IF
470 SH3TREE(2,II)=I1
471 SH3TREE(3,II)=LEVEL
472 IF(I1/=0)THEN
473 SH3TREE(1,I1)=II
474 SH3TREE(1,I2)=II
475 SH3TREE(1,I3)=II
476 SH3TREE(1,I4)=II
477 END IF
478 SH3TRIM(II)=ITRIM
479 ENDDO
480
481 ENDIF
482
483 ENDDO
484
485 DEALLOCATE(LIST)
486 DEALLOCATE(INDEXL)
487
488 RETURN
489 END
hm_get_int_array_index
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
Definition hm_get_int_array_index.F:164
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:523
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:251
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
hm_option_read_mod
Definition hm_option_read_mod.F:408
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1000
submodel_mod
Definition submodel_mod.F:712
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
state_admesh
subroutine state_admesh(ipart, ipartc, iparttg, ixc, ixtg, sh4tree, sh3tree, sh4trim, sh3trim, lsubmodel)
Definition state_admesh.F:41
submodel_mod::submodel_data
Definition submodel_mod.F:723