OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
printgroup.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!|| printgroup ../starter/source/output/outfile/printgroup.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!||--- uses -----------------------------------------------------
30!|| format_mod ../starter/share/modules1/format_mod.F90
31!|| message_mod ../starter/share/message_module/message_mod.F
32!|| r2r_mod ../starter/share/modules1/r2r_mod.F
33!||====================================================================
34 SUBROUTINE printgroup(
35 1 ITAB ,ITABM1 ,IGRNOD ,NINICRACK,BUFSF ,
36 2 IGRSURF ,IGRSLIN ,INICRACK ,IXS ,IXQ ,
37 3 IXC ,IXT ,IXP ,IXR ,IXTG ,
38 4 IXS10 ,IXS20 ,IXS16 ,IPART ,
39 5 KXX ,IXIG3D ,KXIG3D ,
40 6 IGRQUAD ,IGRBRIC ,IGRSH4N ,IGRSH3N ,IGRTRUSS,
41 7 IGRBEAM ,IGRSPRING,IGRPART )
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE groupdef_mod
46 USE message_mod
47 USE r2r_mod
48 USE optiondef_mod
49 USE names_and_titles_mod , ONLY : nchartitle
50 USE format_mod , ONLY : fmw_10i, fmw_6i, fmw_2i_x_f, fmw_4i
51 use element_mod , only : nixs,nixq,nixc,nixt,nixp,nixr,nixtg
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C D u m m y A r g u m e n t s
58C-----------------------------------------------
59 INTEGER ITAB(*),ITABM1(*),NINICRACK,IXS(NIXS,*),
60 . IXC(NIXC,*),IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),
61 . IXTG(NIXTG,*),IXS10(6,*),IXS20(12,*),IXS16(8,*),
62 . IPART(LIPART1,*),IXQ(NIXQ,*),
63 . KXX(NIXX,*),KXIG3D(NIXIG3D,*),IXIG3D(*)
64 my_real BUFSF(*)
65C-----------------------------------------------
66 TYPE (GROUP_) , DIMENSION(NGRNOD) :: IGRNOD
67 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
68 TYPE (GROUP_) , DIMENSION(NGRQUAD) :: IGRQUAD
69 TYPE (GROUP_) , DIMENSION(NGRSHEL) :: IGRSH4N
70 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
71 TYPE (GROUP_) , DIMENSION(NGRTRUS) :: IGRTRUSS
72 TYPE (GROUP_) , DIMENSION(NGRBEAM) :: IGRBEAM
73 TYPE (GROUP_) , DIMENSION(NGRSPRI) :: IGRSPRING
74 TYPE (GROUP_) , DIMENSION(NGRPART) :: IGRPART
75 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
76 TYPE (SURF_) , DIMENSION(NSLIN) :: IGRSLIN
77 TYPE (INICRACK_), DIMENSION(NINICRACK) :: INICRACK
78C-----------------------------------------------
79C C o m m o n B l o c k s
80C-----------------------------------------------
81#include "scr17_c.inc"
82#include "scr03_c.inc"
83#include "scr23_c.inc"
84#include "com01_c.inc"
85#include "com04_c.inc"
86#include "units_c.inc"
87#include "param_c.inc"
88#include "r2r_c.inc"
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER I,J,K,L,II,KK,NEL,NNOD,NSEG,IAD0,DGR,ITYP,SRFTYP,NE,
93 . KR,IADR,NSEGIGE,ID
94 CHARACTER(LEN=NCHARTITLE) :: TITR
95 CHARACTER MESS*40,FMTT*15
96C-----------------------------------------------
97 IF (IPRI >= 0) then
98C-----------------------------------------------
99 IF (ngrbric > 0) THEN
100c MESS='BRIC. ELEMENT GROUP'
101c WRITE(ISTDO,'(2A)')' .. ',MESS
102 DO i = 1,ngrbric
103 id = igrbric(i)%ID
104 titr = igrbric(i)%TITLE
105 nel = igrbric(i)%NENTITY
106 IF (nsubdom>0) THEN
107C-----------> Multidomaines -> on prend le nb d'elts avant split----
108 IF ((nel==0).AND.(igrbric(i)%R2R_ALL/=0)) GOTO 150
109 ENDIF
110C----------------------------------------------------------------------
111 IF (nel == 0 .AND. igrbric(i)%SET_GROUP == 0) THEN
112 CALL ancmsg(msgid=534,
113 . msgtype=msgwarning,
114 . anmode=aninfo_blind_1,
115 . c1='BRICK',c2='BRICK',c3=trim(titr),
116 . i1=id)
117 ENDIF
118 IF (ipri > 0) THEN
119 WRITE(iout,'(//)')
120 WRITE(iout,'(A,I10,3A,I6,A)')' BRIC. GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
121 WRITE(iout,'(A)')' --------------------'
122 WRITE(iout,fmt=fmw_10i)(ixs(nixs,igrbric(i)%ENTITY(j)),j=1,nel)
123 END IF
124150 CONTINUE
125 ENDDO
126 ENDIF
127C-----------------------------------------------
128 IF (ngrquad > 0) THEN
129c MESS='QUAD ELEMENT GROUP'
130c WRITE(ISTDO,'(2A)')' .. ',MESS
131 DO i = 1,ngrquad
132 id = igrquad(i)%ID
133 titr = igrquad(i)%TITLE
134 nel =igrquad(i)%NENTITY
135 IF (nsubdom>0) THEN
136C-----------> Multidomaines -> on prend le nb d'elts avant split----
137 IF ((nel==0).AND.(igrquad(i)%R2R_ALL/=0)) GOTO 160
138 ENDIF
139C----------------------------------------------------------------------
140 IF (nel == 0 .AND. igrquad(i)%SET_GROUP == 0) THEN
141 CALL ancmsg(msgid=534,
142 . msgtype=msgwarning,
143 . anmode=aninfo_blind_1,
144 . c1='QUAD',c2='QUAD',c3=trim(titr),
145 . i1=id)
146 ENDIF
147 IF (ipri > 0) THEN
148 WRITE(iout,'(//)')
149 WRITE(iout,'(A,I10,3A,I6,A)')' QUAD GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
150 WRITE(iout,'(A)')' --------------------'
151 WRITE(iout,fmt=fmw_10i)(ixq(nixq,igrquad(i)%ENTITY(j)),j=1,nel)
152 END IF
153160 CONTINUE
154 ENDDO
155 ENDIF
156C-----------------------------------------------
157 IF (ngrshel > 0) THEN
158c MESS ='SHEL. ELEMENT GROUP'
159c WRITE(ISTDO,'(2A)')' .. ',MESS
160 DO i = 1,ngrshel
161 id = igrsh4n(i)%ID
162 titr = igrsh4n(i)%TITLE
163 nel =igrsh4n(i)%NENTITY
164 IF (nsubdom>0) THEN
165C-----------> Multidomaines -> on prend le nb d'elts avant split----
166 IF ((nel==0).AND.(igrsh4n(i)%R2R_ALL/=0)) GOTO 170
167 ENDIF
168C----------------------------------------------------------------------
169 IF (nel == 0 .AND. igrsh4n(i)%SET_GROUP == 0) THEN
170 CALL ancmsg(msgid=534,
171 . msgtype=msgwarning,
172 . anmode=aninfo_blind_1,
173 . c1='SHELL',c2='SHELL',c3=trim(titr),
174 . i1=id)
175 ENDIF
176 IF (ipri > 0) THEN
177 WRITE(iout,'(A)') ' '
178 WRITE(iout,'(A)') ' '
179 WRITE(iout,'(A,I10,3A,I6,A)')' SHEL. GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
180 WRITE(iout,'(A)')' --------------------'
181 WRITE(iout,fmt=fmw_10i) (ixc(nixc,igrsh4n(i)%ENTITY(j)),j=1,nel)
182 END IF
183170 CONTINUE
184 ENDDO
185 ENDIF
186C-----------------------------------------------
187 IF (ngrtrus > 0) THEN
188c MESS ='TRUS ELEMENT GROUP'
189c WRITE(ISTDO,'(2A)')' .. ',MESS
190 DO i = 1,ngrtrus
191 id = igrtruss(i)%ID
192 titr = igrtruss(i)%TITLE
193 nel =igrtruss(i)%NENTITY
194 IF (nsubdom>0) THEN
195C-----------> Multidomaines -> on prend le nb d'elts avant split----
196 IF ((nel==0).AND.(igrtruss(i)%R2R_ALL/=0)) GOTO 180
197 ENDIF
198C----------------------------------------------------------------------
199 IF (nel == 0.AND. igrtruss(i)%SET_GROUP == 0) THEN
200 CALL ancmsg(msgid=534,
201 . msgtype=msgwarning,
202 . anmode=aninfo_blind_1,
203 . c1='TRUS',c2='TRUS',c3=trim(titr),
204 . i1=id)
205 ENDIF
206 IF (ipri > 0) THEN
207 WRITE(iout,'(//)')
208 WRITE(iout,'(A,I10,3A,I6,A)')' TRUSS GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
209 WRITE(iout,'(A)')' --------------------'
210 WRITE(iout,fmt=fmw_10i)(ixt(nixt,igrtruss(i)%ENTITY(j)),j=1,nel)
211 END IF
212180 CONTINUE
213 ENDDO
214 ENDIF
215C-----------------------------------------------
216 IF (ngrbeam > 0) THEN
217c MESS ='BEAM ELEMENT GROUP'
218c WRITE(ISTDO,'(2A)')' .. ',MESS
219 DO i = 1,ngrbeam
220 id = igrbeam(i)%ID
221 nel =igrbeam(i)%NENTITY
222 titr = igrbeam(i)%TITLE
223 IF (nsubdom>0) THEN
224C-----------> Multidomaines -> on prend le nb d'elts avant split----
225 IF ((nel==0).AND.(igrbeam(i)%R2R_ALL/=0)) GOTO 190
226 ENDIF
227C----------------------------------------------------------------------
228 IF (nel == 0 .AND. igrbeam(i)%SET_GROUP == 0) THEN
229 CALL ancmsg(msgid=534,
230 . msgtype=msgwarning,
231 . anmode=aninfo_blind_1,
232 . c1='BEAM',c2='BEAM',c3=trim(titr),
233 . i1=id)
234 ENDIF
235 IF (ipri > 0) THEN
236 WRITE(iout,'(//)')
237 WRITE(iout,'(A,I10,3A,I6,A)')' BEAM GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
238 WRITE(iout,'(A)')' --------------------'
239 WRITE(iout,fmt=fmw_10i)(ixp(nixp,igrbeam(i)%ENTITY(j)),j=1,nel)
240 END IF
241190 CONTINUE
242 ENDDO
243 ENDIF
244C-----------------------------------------------
245 IF (ngrspri > 0) THEN
246c MESS ='SPRI ELEMENT GROUP'
247c WRITE(ISTDO,'(2A)')' .. ',MESS
248 DO i = 1,ngrspri
249 id = igrspring(i)%ID
250 titr = igrspring(i)%TITLE
251 nel =igrspring(i)%NENTITY
252 IF (nsubdom>0) THEN
253C-----------> Multidomaines -> on prend le nb d'elts avant split----
254 IF ((nel==0).AND.(igrspring(i)%R2R_ALL/=0)) GOTO 200
255 ENDIF
256C----------------------------------------------------------------------
257 IF (nel == 0 .AND. igrspring(i)%SET_GROUP == 0) THEN
258 CALL ancmsg(msgid=534,
259 . msgtype=msgwarning,
260 . anmode=aninfo_blind_1,
261 . c1='SPRI',c2='SPRI',c3=trim(titr),
262 . i1=id)
263 ENDIF
264 IF (ipri > 0) THEN
265 WRITE(iout,'(//)')
266 WRITE(iout,'(A,I10,3A,I6,A)')' SPRING GROUP:',id,',',trim(titr),',',nel,' ELEMENTS: '
267 WRITE(iout,'(A)')' --------------------'
268 WRITE(iout,fmt=fmw_10i)(ixr(nixr,igrspring(i)%ENTITY(j)),j=1,nel)
269 END IF
270200 CONTINUE
271 ENDDO
272 ENDIF
273C-----------------------------------------------
274 IF (ngrsh3n > 0) THEN
275c MESS ='SHELL3N ELEMENT GROUP'
276c WRITE(ISTDO,'(2A)')' .. ',MESS
277 DO i = 1,ngrsh3n
278 id = igrsh3n(i)%ID
279 titr = igrsh3n(i)%TITLE
280 nel =igrsh3n(i)%NENTITY
281 IF (nsubdom>0) THEN
282C-----------> Multidomaines -> on prend le nb d'elts avant split----
283 IF ((nel==0).AND.(igrsh3n(i)%R2R_ALL/=0)) GOTO 210
284 ENDIF
285C----------------------------------------------------------------------
286 IF (nel == 0 .AND. igrsh3n(i)%SET_GROUP == 0) THEN
287 IF(n2d==0)
288 . CALL ancmsg(msgid=534,
289 . msgtype=msgwarning,
290 . anmode=aninfo_blind_1,
291 . c1='SH3N',c2='SH3N',c3=trim(titr),
292 . i1=id)
293 IF(n2d/=0)
294 . CALL ancmsg(msgid=534,
295 . msgtype=msgwarning,
296 . anmode=aninfo_blind_1,
297 . c1='TRIA',c2='TRIA',c3=trim(titr),
298 . i1=id)
299 ENDIF
300 IF (ipri > 0) THEN
301 WRITE(iout,'(A)') ' '
302 WRITE(iout,'(A)') ' '
303 IF(n2d==0)WRITE(iout,'(A,I10,3A,I6,A)')' SH3N. GROUP:', id,',',trim(titr),',',nel,' ELEMENTS: '
304 IF(n2d/=0)WRITE(iout,'(A,I10,3A,I6,A)')' TRIA. GROUP:', id,',',trim(titr),',',nel,' ELEMENTS: '
305 WRITE(iout,'(A)')' --------------------'
306 WRITE(iout,fmt=fmw_10i)(ixtg(nixtg,igrsh3n(i)%ENTITY(j)),j=1,nel)
307 END IF
308210 CONTINUE
309 ENDDO
310 ENDIF
311C-----------------------------------------------
312 IF (ngrpart > 0) THEN
313c MESS ='PART GROUP'
314c WRITE(ISTDO,'(2A)')' .. ',MESS
315 DO i = 1,ngrpart
316 id = igrpart(i)%ID
317 titr = igrpart(i)%TITLE
318 nel =igrpart(i)%NENTITY
319 IF (nsubdom>0) THEN
320C-----------> Multidomaines -> on prend le nb de parts avant split-----
321 IF ((nel==0).AND.(igrpart(i)%R2R_ALL/=0)) GOTO 220
322 ENDIF
323C----------------------------------------------------------------------
324 IF (nel == 0 .AND. igrpart(i)%SET_GROUP == 0) THEN
325 CALL ancmsg(msgid=534,
326 . msgtype=msgwarning,
327 . anmode=aninfo_blind_1,
328 . c1='PART',c2='PART',c3=trim(titr),
329 . i1=id)
330 ENDIF
331 IF (ipri > 0) THEN
332 WRITE(iout,'(//)')
333 WRITE(iout,'(A,I10,3A,I6,A)')' PART GROUP:',id,',',trim(titr),',',nel,' PARTS: '
334 WRITE(iout,'(A)')' ---------------------'
335 WRITE(iout,fmt=fmw_10i)(ipart(4,igrpart(i)%ENTITY(j)),j=1,nel)
336 END IF
337220 CONTINUE
338 ENDDO
339 ENDIF
340C-----------------------------------------------
341 IF (nsurf > 0)THEN
342 DO i=1,nsurf
343 srftyp=igrsurf(i)%TYPE
344 IF (srftyp == 100) THEN
345 IF (igrsurf(i)%ID_MADYMO == 0) THEN
346 CALL ancmsg(msgid=692,
347 . msgtype=msgwarning,
348 . anmode=aninfo_blind_1,
349 . i1=igrsurf(i)%ID)
350 END IF
351 ELSEIF (srftyp == 101) THEN
352 CONTINUE
353 ELSEIF (srftyp == 200) THEN
354 CONTINUE
355 ELSE
356 nseg=igrsurf(i)%NSEG
357 nsegige=igrsurf(i)%NSEG_IGE
358C-----------> Multidomaines -> on prend le nb de segments avant split----
359 IF (nsubdom>0) nseg=isurf_r2r(2,i)
360C----------------------------------------------------------------------
361 IF (nseg == 0 .AND. nsegige == 0 .AND. igrsurf(i)%SET_GROUP == 0) THEN
362 CALL ancmsg(msgid=692,
363 . msgtype=msgwarning,
364 . anmode=aninfo_blind_1,
365 . i1=igrsurf(i)%ID)
366 ENDIF
367 ENDIF
368 ENDDO
369C---
370 IF (ipri > 0) THEN
371 DO i=1,nsurf
372C-----------> Multidomaines --------------------------------------------
373 IF (nsubdom>0) THEN
374 IF((igrsurf(i)%NSEG==0).AND.(isurf_r2r(2,i)/=0)) GOTO 300
375 ENDIF
376C----------------------------------------------------------------------
377 srftyp=igrsurf(i)%TYPE
378 IF (srftyp == 100) THEN
379 WRITE(iout,'(A)') ' '
380 WRITE(iout,'(A)') ' '
381 titr = igrsurf(i)%TITLE
382 WRITE(iout,'(A,I10,1X,A)')' EXTERNAL COUPLING TO SURFACE:',igrsurf(i)%ID,trim(titr)
383 WRITE(iout,'(A)') ' ----------------------------'
384 WRITE(iout,1000) igrsurf(i)%ID_MADYMO
385 WRITE(iout,'(10X,A)')'CORRESPONDING DATA WILL BE READ IN RADIOSS ENGINE.'
386 ELSEIF (srftyp == 101) THEN
387 WRITE(iout,'(A)') ' '
388 WRITE(iout,'(A)') ' '
389 titr = igrsurf(i)%TITLE
390 WRITE(iout,'(A,I10,1X,A)')' ELLIPSOIDAL SURFACE:',igrsurf(i)%ID,trim(titr)
391 WRITE(iout,'(A)') ' --------------------'
392 dgr=bufsf(igrsurf(i)%IAD_BUFR+36)
393 WRITE(iout,1010) bufsf(igrsurf(i)%IAD_BUFR+4),bufsf(igrsurf(i)%IAD_BUFR+5),
394 . bufsf(igrsurf(i)%IAD_BUFR+6),bufsf(igrsurf(i)%IAD_BUFR+1),
395 . bufsf(igrsurf(i)%IAD_BUFR+2),bufsf(igrsurf(i)%IAD_BUFR+3),
396 . dgr,igrsurf(i)%ID_MADYMO,
397 . (bufsf(igrsurf(i)%IAD_BUFR+7+j-1),j=1,3),
398 . (bufsf(igrsurf(i)%IAD_BUFR+7+j-1),j=4,6),
399 . (bufsf(igrsurf(i)%IAD_BUFR+7+j-1),j=7,9)
400 igrsurf(i)%ID_MADYMO=0
401 ELSEIF (srftyp == 200) THEN
402 WRITE(iout,'(A)') ' '
403 WRITE(iout,'(A)') ' '
404 titr = igrsurf(i)%TITLE
405 WRITE(iout,'(A,I10,A,A)')' INFINITE PLANE SURFACE:',igrsurf(i)%ID,',',trim(titr)
406 WRITE(iout,'(A)') ' --------------------'
407 WRITE(iout,'(A)') ' NORMAL PLANE VECTOR :'
408 WRITE(iout,1020)(bufsf(igrsurf(i)%IAD_BUFR+j),j=1,6)
409 ELSE
410 nseg=igrsurf(i)%NSEG
411!! IF (NSEG/=0) THEN
412 WRITE(iout,'(A)') ' '
413 WRITE(iout,'(A)') ' '
414 titr = igrsurf(i)%TITLE
415 WRITE(iout,'(A,I10,3A,I10,A)')' SURFACE:',igrsurf(i)%ID,',',trim(titr),',',nseg,' SEGMENTS: '
416 WRITE(iout,'(A)')' -------------------'
417 WRITE(iout,'(A)')' NODE 1 NODE 2 NODE 3 NODE 4 ELTYP ELEMENT'
418 DO j=1,nseg
419 ityp=igrsurf(i)%ELTYP(j)
420 IF (ityp == 1) THEN
421 ne=ixs(nixs,igrsurf(i)%ELEM(j))
422 ELSEIF (ityp == 2) THEN
423 ne=ixq(nixq,igrsurf(i)%ELEM(j))
424 ELSEIF (ityp == 3) THEN
425 ne=ixc(nixc,igrsurf(i)%ELEM(j))
426 ELSEIF (ityp == 7) THEN
427 ne=ixtg(nixtg,igrsurf(i)%ELEM(j))
428 ELSE
429 ne=0
430 ENDIF
431 IF(n2d == 0)THEN
432 WRITE(iout,'(6I10)')(itab(igrsurf(i)%NODES(j,l)),l=1,4),ityp,ne
433 ELSE
434 WRITE(iout,'(2I10,16X,2I10)')(itab(igrsurf(i)%NODES(j,l)),l=1,2),ityp,ne
435 ENDIF
436 ENDDO
437!! ENDIF ! IF (NSEG/=0)
438 IF (nsegige/=0) THEN
439 nsegige=igrsurf(i)%NSEG_IGE
440 WRITE(iout,'(a)') ' '
441 WRITE(IOUT,'(a)') ' '
442 TITR = IGRSURF(I)%TITLE
443 WRITE(IOUT,'(a,i10,3a,i10,a)')' surface:',IGRSURF(I)%ID,',',TRIM(TITR),',',NSEGIGE,' segments isogeometriques: '
444 WRITE(IOUT,'(a)')' -------------------'
445 WRITE(IOUT,'(a)')' cp 1 cp 2 cp 3 cp 4 eltyp element'
446 DO J=1,NSEGIGE
447 ITYP=IGRSURF(I)%ELTYP_IGE(J)
448 IF (ITYP == 101) THEN
449 NE=KXIG3D(5,IGRSURF(I)%ELEM_IGE(J))
450 ELSE
451 NE=0
452 ENDIF
453 IF(N2D == 0)THEN
454 WRITE(IOUT,FMT=FMW_6I)(IGRSURF(I)%NODES_IGE(J,L),L=1,4),ITYP,NE
455 ELSE
456 WRITE(IOUT,'(2i10,16x,2i10)')(IGRSURF(I)%NODES_IGE(J,L),L=1,2),ITYP,NE
457 ENDIF
458 ENDDO
459 ENDIF ! IF (NSEGIGE/=0)
460 ENDIF ! IF (SRFTYP == 100)
461300 CONTINUE
462 ENDDO ! DO I=1,NSURF
463 ENDIF ! IF (IPRI > 0)
464 ENDIF ! IF (NSURF > 0)
465C-----------------------------------------------
466 IF (NSLIN > 0) THEN
467 DO I=1,NSLIN
468 NSEG=IGRSLIN(I)%NSEG
469C-----------> Multidomaines -> on prend le nb de segments avant split----
470 IF (NSUBDOM>0) NSEG=IGRSLIN(I)%NSEG_R2R_ALL
471C-----------> Type19 interface -> warning deactivated for lines generated in translator
472 TITR = IGRSLIN(I)%TITLE
473 IF (TITR(1:19)=='type19 interface nb') CYCLE
474C----------------------------------------------------------------------
475.AND. IF (NSEG == 0 IGRSLIN(I)%SET_GROUP == 0) THEN
476 CALL ANCMSG(MSGID=691,
477 . MSGTYPE=MSGWARNING,
478 . ANMODE=ANINFO_BLIND_1,
479 . I1=IGRSLIN(I)%ID)
480 END IF
481 END DO
482C---
483 DO I=1,NSLIN
484 NSEG=IGRSLIN(I)%NSEG
485 IF (NSUBDOM>0) THEN
486C-----------> Multidomaines --------------------------------------------
487.AND. IF ((NSEG==0)(IGRSLIN(I)%NSEG_R2R_ALL/=0)) GOTO 310
488 ENDIF
489C----------------------------------------------------------------------
490 IF (IPRI > 0) THEN
491 WRITE(IOUT,'(a)') ' '
492 WRITE(IOUT,'(a)') ' '
493 TITR = IGRSLIN(I)%TITLE
494 WRITE(IOUT,'(a,i10,3a,i10,a)')' line: ',IGRSLIN(I)%ID,',',TRIM(TITR),',',NSEG,' segments: '
495 WRITE(IOUT,'(a)')' ----------------------'
496 WRITE(IOUT,'(a)')' node 1 node 2 eltyp element'
497 DO J=1,NSEG
498 ITYP=IGRSLIN(I)%ELTYP(J)
499 IF(ITYP == 1)THEN
500 NE=IXS(NIXS,IGRSLIN(I)%ELEM(J))
501 ELSEIF(ITYP == 2)THEN
502 NE=IXQ(NIXQ,IGRSLIN(I)%ELEM(J))
503 ELSEIF(ITYP == 3)THEN
504 NE=IXC(NIXC,IGRSLIN(I)%ELEM(J))
505 ELSEIF(ITYP == 7)THEN
506 NE=IXTG(NIXTG,IGRSLIN(I)%ELEM(J))
507 ELSEIF(ITYP == 4)THEN
508 NE=IXT(NIXT,IGRSLIN(I)%ELEM(J))
509 ELSEIF(ITYP == 5)THEN
510 NE=IXP(NIXP,IGRSLIN(I)%ELEM(J))
511 ELSEIF(ITYP == 6)THEN
512 NE=IXR(NIXR,IGRSLIN(I)%ELEM(J))
513 ELSEIF(ITYP == 8)THEN
514 NE=KXX(NIXX,IGRSLIN(I)%ELEM(J))
515 ELSE
516 NE=0
517 ENDIF
518 WRITE(IOUT,FMT=FMW_4I)(ITAB(IGRSLIN(I)%NODES(J,L)),L=1,2),ITYP,NE
519 ENDDO
520 END IF
521310 CONTINUE
522 ENDDO
523 ENDIF
524C-----------------------------------------------
525 IF (NINICRACK > 0) THEN
526 DO I=1,NINICRACK
527 NSEG=INICRACK(I)%NSEG
528 IF (NSEG==0) THEN
529 CALL ANCMSG(MSGID=693,
530 . MSGTYPE=MSGWARNING,
531 . ANMODE=ANINFO_BLIND_1,
532 . I1=INICRACK(I)%ID)
533 END IF
534 END DO
535C
536 DO I=1,NINICRACK
537 NSEG=INICRACK(I)%NSEG
538 IF (IPRI > 0) THEN
539 WRITE(IOUT,'(a)') ' '
540 WRITE(IOUT,'(a)') ' '
541 TITR = INICRACK(I)%TITLE
542 WRITE(IOUT,'(a,i10,3a)')' inicrack: ',INICRACK(I)%ID,',',TRIM(TITR)
543 WRITE(IOUT,'(a)')' ----------------------'
544 WRITE(IOUT,'(a)')' node 1 node 2 ratio'
545 DO J=1,NSEG
546 WRITE(IOUT,FMT=FMW_2I_X_F)(ITAB(INICRACK(I)%SEG(J)%NODES(L)),L=1,2),INICRACK(I)%SEG(J)%RATIO
547 ENDDO
548 END IF
549 ENDDO
550 ENDIF
551C-----------------------------------------------
552 IF (NGRNOD > 0)THEN
553 DO I=1,NGRNOD
554 ID = IGRNOD(I)%ID
555 NNOD=IGRNOD(I)%NENTITY
556 IF (NSUBDOM>0) NNOD=IGRNOD(I)%R2R_ALL
557C-----------> Multidomaines -> on prend le nb de noeuds avant split----
558C----------------------------------------------------------------------
559 IF (NNOD==0) THEN
560 CALL ANCMSG(MSGID=690,
561 . MSGTYPE=MSGWARNING,
562 . ANMODE=ANINFO_BLIND_1,
563 . I1=ID)
564 END IF
565 ENDDO
566C
567 DO I=1,NGRNOD
568 ID = IGRNOD(I)%ID
569 NNOD=IGRNOD(I)%NENTITY
570 IF (NSUBDOM>0) THEN
571C-----------> Multidomaines --------------------------------------------
572.AND. IF ((NNOD==0)(IGRNOD(I)%R2R_ALL/=0)) GOTO 350
573 ENDIF
574C----------------------------------------------------------------------
575 IF (IPRI > 0) THEN
576 WRITE(IOUT,'(a)') ' '
577 WRITE(IOUT,'(a)') ' '
578 TITR = IGRNOD(I)%TITLE
579 WRITE(IOUT,'(a,i10,3a,i10,a)')' node group:',ID,',',TRIM(TITR),',',NNOD,' nodes: '
580 WRITE(IOUT,'(a)')' -------------------'
581 WRITE(IOUT,FMT=FMW_10I)(ITAB(IGRNOD(I)%ENTITY(J)),J=1,NNOD)
582 END IF
583350 CONTINUE
584 ENDDO
585 ENDIF
586C-----------
587 ENDIF
588C-----------
5891000 FORMAT(/10X,'surface definition : '
590 . /10X,'-------------------- '
591 . /10X,'reference to madymo coupling surface : ',I8)
5921010 FORMAT(/10X,'ellipsoid definition : '
593 . /10X,'---------------------- '
594 . /10X,'ellipsoid center ',1P3G20.13
595 . /10X,'ellipsoid semi-axes length ',1P3G20.13
596 . /10X,'ellipsoid degree ',I10
597 . /10X,'skew number for positioning ',I10
598 . /10X,'computed ellipsoid rotation ',
599 . /10X,' with respect to skew : ',1P3G20.13
600 . /10X,' ',1P3G20.13
601 . /10X,' ',1P3G20.13)
6021020 FORMAT(/10X,'plane definition : '
603 . /10X,'---------------------- '
604 . /48X,'x-ord y-ord z-ord'
605 . /10X,'tail coordinates ',1P3G20.13
606 . /10X,'head coordinates ',1P3G20.13)
607C-----------
608
609 RETURN
610 END
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
groupdef_mod
Definition groupdef_mod.F:668
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
optiondef_mod
Definition optiondef_mod.F:90
r2r_mod
Definition r2r_mod.F:125
r2r_mod::isurf_r2r
integer, dimension(:,:), allocatable isurf_r2r
Definition r2r_mod.F:143
printgroup
subroutine printgroup(itab, itabm1, igrnod, ninicrack, bufsf, igrsurf, igrslin, inicrack, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, ixs10, ixs20, ixs16, ipart, kxx, ixig3d, kxig3d, igrquad, igrbric, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, igrpart)
Definition printgroup.F:42
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:895