OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
printgroup.F File Reference
#include "implicit_f.inc"
#include "scr17_c.inc"
#include "scr03_c.inc"
#include "scr23_c.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "units_c.inc"
#include "param_c.inc"
#include "r2r_c.inc"

Go to the source code of this file.

Functions/Subroutines

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)

Function/Subroutine Documentation

◆ printgroup()

subroutine printgroup ( integer, dimension(*) itab,
integer, dimension(*) itabm1,
type (group_), dimension(ngrnod) igrnod,
integer ninicrack,
bufsf,
type (surf_), dimension(nsurf) igrsurf,
type (surf_), dimension(nslin) igrslin,
type (inicrack_), dimension(ninicrack) inicrack,
integer, dimension(nixs,*) ixs,
integer, dimension(nixq,*) ixq,
integer, dimension(nixc,*) ixc,
integer, dimension(nixt,*) ixt,
integer, dimension(nixp,*) ixp,
integer, dimension(nixr,*) ixr,
integer, dimension(nixtg,*) ixtg,
integer, dimension(6,*) ixs10,
integer, dimension(12,*) ixs20,
integer, dimension(8,*) ixs16,
integer, dimension(lipart1,*) ipart,
integer, dimension(nixx,*) kxx,
integer, dimension(*) ixig3d,
integer, dimension(nixig3d,*) kxig3d,
type (group_), dimension(ngrquad) igrquad,
type (group_), dimension(ngrbric) igrbric,
type (group_), dimension(ngrshel) igrsh4n,
type (group_), dimension(ngrsh3n) igrsh3n,
type (group_), dimension(ngrtrus) igrtruss,
type (group_), dimension(ngrbeam) igrbeam,
type (group_), dimension(ngrspri) igrspring,
type (group_), dimension(ngrpart) igrpart )

Definition at line 34 of file printgroup.F.

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