44
45
46
53
54
55
56#include "implicit_f.inc"
57
58
59
60#include "com04_c.inc"
61#include "scr17_c.inc"
62
63
64
65 INTEGER, INTENT(IN) :: ITAB(*),IPART(LIPART1,*),IXC(NIXC,*),
66 . IXTG(NIXTG,*),IXQ(NIXQ,*),IXP(NIXP,*),IXT(NIXT,*),IXR(NIXR,*),
67 . IXS(NIXS,*)
68
69 TYPE (SET_) , DIMENSION(NSETS) :: SET
70 TYPE (SUBMODEL_DATA),INTENT(IN)::LSUBMODEL(*)
71 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRNOD) :: IGRNOD
72 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRPART) :: IGRPART
73 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRBRIC) :: IGRBRIC
74 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRSHEL) :: IGRSH4N
75 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRSH3N) :: IGRSH3N
76 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRQUAD) :: IGRQUAD
77 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRBEAM) :: IGRBEAM
78 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRTRUS) :: IGRTRUSS
79 TYPE (GROUP_) , INTENT(IN), DIMENSION(NGRSPRI) :: IGRSPRING
80 TYPE (SURF_) , INTENT(IN), DIMENSION(NSURF) :: IGRSURF
81 TYPE (SURF_) , INTENT(IN), DIMENSION(NSLIN) :: IGRSLIN
82
83
84
85 INTEGER I,J,K,N,ID,IGS,CLAUSES_MAX,ISET_TYPE,ITMP,ICODE,IDS_MAX,IDS,
86 . OPT_,OPT_D,OPT_O,OPT_G,OPT_B,OPT_A,OPT_E,OPT_I,OPT_C,
87 . IGR,NOD,NB_NODE,NB_PART,IP,IE,NB_SOLID,NB_SH4N,NB_SH3N,
88 . NB_QUAD,NB_TRIA,NB_BEAM,NB_TRUSS,NB_SPRING,NB_SURF_SEG,
89 . NB_LINE_SEG,ISEG,NB_NODENS
90 CHARACTER(LEN = nchartitle) :: TITLE
91 CHARACTER(LEN = ncharfield) :: KEYSET,SET_TYPE
92 CHARACTER(LEN = ncharkey) :: KEY
93 CHARACTER (LEN=255) :: VARNAME
94 DOUBLE PRECISION TEMP_DOUBLE
95 LOGICAL IS_AVAILABLE
96
97 is_available = .false.
98
99
100
101
103
105
107
109
111 WRITE(varname,'(A)') trim(title)
112 CALL qaprint(title(1:len_trim(varname)),
id, 0.0_8)
113
115
116! issue 'SET_TYPE' ---> read one more character than the SET_TYPE
117
118 itmp = len(trim(set_type))
119 IF (itmp > 0) THEN
120 icode = iachar(set_type(itmp:itmp))
121 IF (icode == 0) set_type(itmp:itmp)=' '
122 ENDIF
123
124 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//trim(set_type)
125 CALL qaprint(varname(1:len_trim(varname)),
id, 0.0_8)
126
127 CALL hm_get_intv(
'iset_Type', iset_type,is_available,lsubmodel)
128 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'iset_Type'
129 CALL qaprint(varname(1:len_trim(varname)), iset_type, 0.0_8)
130
131 CALL hm_get_intv(
'clausesmax',clauses_max,is_available,lsubmodel)
132
133 DO j=1,clauses_max
135
136
137
138 itmp = len(trim(keyset))
139 icode = iachar(keyset(itmp:itmp))
140 IF (icode == 0) keyset(itmp:itmp)=' '
141
142
143
152
153
154
155
156
157 IF (opt_d == 1) THEN
158 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_D'
159 CALL qaprint(varname(1:len_trim(varname)),opt_d,0.0_8
160 ENDIF
161 IF (opt_o == 1) THEN
162 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_O'
163 CALL qaprint(varname(1:len_trim(varname)),opt_o,0.0_8)
164 ENDIF
165 IF (opt_g == 1) THEN
166 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_G'
167 CALL qaprint(varname(1:len_trim(varname)),opt_g,0.0_8)
168 ENDIF
169 IF (opt_b == 1) THEN
170 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_B'
171 CALL qaprint(varname(1:len_trim(varname
172 ENDIF
173 IF (opt_a == 1) THEN
174 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_A'
175 CALL qaprint(varname(1:len_trim(varname)),opt_a,0.0_8)
176 ENDIF
177 IF (opt_e == 1) THEN
178 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_E'
179 CALL qaprint(varname(1:len_trim(varname)),opt_e,0.0_8
180 ENDIF
181 IF (opt_i == 1) THEN
182 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_I'
183 CALL qaprint(varname(1:len_trim(varname)),opt_i,0.0_8)
184 ENDIF
185 IF (opt_c == 1) THEN
186 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'opt_C'
187 CALL qaprint(varname(1:len_trim(varname)),opt_c,0.0_8)
188 ENDIF
189
191 DO k=1,ids_max
193 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//trim(keyset)//
'_',k
194 CALL qaprint(varname(1:len_trim(varname)),ids,0.0_8)
195 ENDDO
196
197 ENDDO
198
199
200
201
202
203
204
205 IF( set(igs)%SET_ACTIV == 0 ) cycle
206
207 nb_node = set(igs)%NB_NODE
208
209 IF (nb_node > 0) THEN
210 igr = set(igs)%SET_GRNOD_ID
211 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRNOD'
212 CALL qaprint(varname(1:len_trim(varname)),igrnod(igr)%ID,0.0_8)
213 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRNOD_NB_NODE'
214 CALL qaprint(varname(1:len_trim(varname)),nb_node,0.0_8)
215 DO n = 1,nb_node
216 nod = igrnod(igr)%ENTITY(n)
217 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'NODE'//
'_',n
218 CALL qaprint(varname(1:len_trim(varname)),itab(nod),0.0_8)
219 ENDDO
220 ENDIF
221
222
223
224 IF( set(igs)%SET_ACTIV == 0 ) cycle
225
226 nb_nodens = set(igs)%NB_NODENS
227
228 IF (nb_nodens > 0) THEN
229 igr = set(igs)%SET_GRNOD_ID
230 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRNOD'
231 CALL qaprint(varname(1:len_trim(varname)),igrnod(igr)%ID,0.0_8)
232 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRNOD_NB_NODENS'
233 CALL qaprint(varname(1:len_trim(varname)),nb_nodens,0.0_8)
234 DO n = 1,nb_nodens
235 nod = igrnod(igr)%ENTITY(n)
236 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'NODENS'//
'_',n
237 CALL qaprint(varname(1:len_trim(varname)),itab(nod),0.0_8)
238 ENDDO
239 ENDIF
240
241
242
243 nb_part = set(igs)%NB_PART
244 IF (nb_part > 0) THEN
245 igr = set(igs)%SET_GRPART_ID
246 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRPART'
247 CALL qaprint(varname(1:len_trim(varname)),igrpart(igr)%ID,0.0_8)
248 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRPART_NB_PART'
249 CALL qaprint(varname(1:len_trim(varname)),nb_part,0.0_8)
250 DO n = 1,nb_part
251 ip = igrpart(igr)%ENTITY(n)
252 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'PART'//
'_',n
253 CALL qaprint(varname(1:len_trim(varname)),ipart(4,ip),0.0_8)
254 ENDDO
255 ENDIF
256
257
258
259
260 nb_solid = set(igs)%NB_SOLID
261 IF (nb_solid > 0) THEN
262 igr = set(igs)%SET_GRSOLID_ID
263 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRBRIC'
264 CALL qaprint(varname(1:len_trim(varname)),igrbric(igr)%ID,0.0_8)
265 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRBRIC_NB_SOLID'
266 CALL qaprint(varname(1:len_trim(varname)),nb_solid,0.0_8)
267 DO n = 1,nb_solid
268 ie = igrbric(igr)%ENTITY(n)
269 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'SOLID'//
'_',n
270 CALL qaprint(varname(1:len_trim(varname)),ixs(nixs,ie),0.0_8)
271 ENDDO
272 ENDIF
273
274
275 nb_sh4n = set(igs)%NB_SH4N
276 IF (nb_sh4n > 0) THEN
277 igr = set(igs)%SET_GRSH4N_ID
278 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRSH4N'
279 CALL qaprint(varname(1:len_trim(varname)),igrsh4n(igr)%ID,0.0_8)
280 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRSH4N_NB_SH4N'
281 CALL qaprint(varname(1:len_trim(varname)),nb_sh4n,0.0_8)
282 DO n = 1,nb_sh4n
283 ie = igrsh4n(igr)%ENTITY(n)
284 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'SHELL'//
'_',n
285 CALL qaprint(varname(1:len_trim(varname)),ixc(nixc,ie),0.0_8)
286 ENDDO
287 ENDIF
288
289
290 nb_sh3n = set(igs)%NB_SH3N
291 IF (nb_sh3n > 0) THEN
292 igr = set(igs)%SET_GRSH3N_ID
293 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRSH3N'
294 CALL qaprint(varname(1:len_trim(varname)),igrsh3n(igr)%ID,0.0_8)
295 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRSH3N_NB_SH3N'
296 CALL qaprint(varname(1:len_trim(varname)),nb_sh3n,0.0_8)
297 DO n = 1,nb_sh3n
298 ie = igrsh3n(igr)%ENTITY(n)
299 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'SH3N'//
'_',n
300 CALL qaprint(varname(1:len_trim(varname)),ixtg(nixtg,ie),0.0_8)
301 ENDDO
302 ENDIF
303
304 ! quad
305 nb_quad = set(igs)%NB_QUAD
306 IF (nb_quad > 0) THEN
307 igr = set(igs)%SET_GRQUAD_ID
308 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRQUAD'
309 CALL qaprint(varname(1:len_trim(varname)),igrquad(igr)%ID,0.0_8)
310 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRQUAD_NB_QUAD'
311 CALL qaprint(varname(1:len_trim(varname)),nb_quad,0.0_8)
312 DO n = 1,nb_quad
313 ie = igrquad(igr)%ENTITY(n)
314 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'QUAD'//
'_',n
315 CALL qaprint(varname(1:len_trim(varname)),ixq(nixq
316 ENDDO
317 ENDIF
318
319
320 nb_tria = set(igs)%NB_TRIA
321 IF (nb_tria > 0) THEN
322 igr = set(igs)%SET_GRTRIA_ID
323 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRTRIA'
324 CALL qaprint(varname(1:len_trim(varname)),igrsh3n(igr)%ID,0.0_8)
325 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRTRIA_NB_TRIA'
326 CALL qaprint(varname(1:len_trim(varname)),nb_tria,0.0_8
327 DO n = 1,nb_tria
328 ie = igrsh3n(igr)%ENTITY(n)
329 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'TRIA'//
'_',n
330 CALL qaprint(varname(1:len_trim(varname)),ixtg(nixtg,ie),0.0_8)
331 ENDDO
332 ENDIF
333
334
335 nb_beam = set(igs)%NB_BEAM
336 IF (nb_beam > 0) THEN
337 igr = set(igs)%SET_GRBEAM_ID
338 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRBEAM'
339 CALL qaprint(varname(1:len_trim(varname)),igrbeam(igr)%ID,0.0_8)
340 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRBEAM_NB_BEAM'
341 CALL qaprint(varname(1:len_trim(varname)),nb_beam,0.0_8)
342 DO n = 1,nb_beam
343 ie = igrbeam(igr)%ENTITY(n)
344 WRITE(varname,
'(A,I0,A,I0)')
'SET_',
id,
'_'//
'BEAM'//
'_',n
345 CALL qaprint(varname(1:len_trim(varname)),ixp(nixp,ie),0.0_8)
346 ENDDO
347 ENDIF
348
349
350 nb_truss = set(igs)%NB_TRUSS
351 IF (nb_truss > 0) THEN
352 igr = set(igs)%SET_GRTRUSS_ID
353 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRTRUSS'
354 CALL qaprint(varname(1:len_trim(varname)),igrtruss(igr)%ID,0.0_8)
355 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'GRTRUSS_NB_TRUSS'
356 CALL qaprint(varname(1:len_trim(varname)),nb_truss,0.0_8)
357 DO n = 1,nb_truss
358 ie = igrtruss(igr)%ENTITY(n)
359 WRITE(varname,'(A,I0,A,I0)')'SET_''//'truss'//'_',N
360 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXT(NIXT,IE),0.0_8)
361 ENDDO
362 ENDIF ! IF (NB_TRUSS > 0)
363!
364 ! spring
365 NB_SPRING = SET(IGS)%NB_SPRING
366 IF (NB_SPRING > 0) THEN
367 IGR = SET(IGS)%SET_GRSPRING_ID
368 WRITE(VARNAME,'(a,i0,a)
') 'set_',ID,'_
'//'grspring
'
369 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IGRSPRING(IGR)%ID,0.0_8)
370 WRITE(VARNAME,'(a,i0,a)
') 'set_',ID,'_
'//'grspring_nb_spring
'
371 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),NB_SPRING,0.0_8)
372 DO N = 1,NB_SPRING
373 IE = IGRSPRING(IGR)%ENTITY(N)
374 WRITE(VARNAME,'(a,i0,a,i0)
') 'set_',ID,'_
'//'spring
'//'_
',N
375 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXR(NIXR,IE),0.0_8)
376 ENDDO
377 ENDIF ! IF (NB_SPRING > 0)
378!---
379! --- New /SET grsurf --
380!---
381 ! surface segments
382 NB_SURF_SEG = SET(IGS)%HAS_SURF_SEG
383 IF (NB_SURF_SEG > 0) THEN
384 IGR = SET(IGS)%SET_NSURF_ID
385 WRITE(VARNAME,'(a,i0,a)
') 'set_',ID,'_
'//'surface
'
386 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IGRSURF(IGR)%ID,0.0_8)
387 WRITE(VARNAME,'(a,i0,a)')
'SET_',
id,
'_'//
'SURFACE_NB_SEG'
388 CALL qaprint(varname(1:len_trim(varname)),nb_surf_seg
389 DO n = 1,nb_surf_seg
390 IF(set(igs)%NB_ELLIPSE == 0 .AND. set(igs)%NB_PLANE== 0)THEN
391 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'SURFACE_SEG_',n,
'_NODE_1'
392 CALL qaprint(varname(1:len_trim(varname
393 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'SURFACE_SEG_',n,
'_NODE_2'
394 CALL qaprint(varname(1:len_trim(varname)),itab(igrsurf(igr)%NODES(n,2)),0.0_8)
395 IF(igrsurf(igr)%NODES(n,3) > 0)THEN
396 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'SURFACE_SEG_',n,
'_NODE_3'
397 CALL qaprint(varname(1:len_trim(varname)),itab(igrsurf
398 ENDIF
399 IF(igrsurf(igr)%NODES(n,4) > 0)THEN
400 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id'_''SURFACE_SEG_''_NODE_4'
401 CALL qaprint(varname(1:len_trim(varname)),itab(igrsurf(igr)%NODES(n,4)),0.0_8)
402 ENDIF
403 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'SURFACE_SEG_',n,
'_ELTYP'
404 CALL qaprint(varname(1:len_trim(varname)),igrsurf(igr)%ELTYP(n),0.0_8)
405 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'SURFACE_SEG_',n,
'_ELEM'
406 IF (igrsurf(igr)%ELTYP(n) == 3 ) THEN
407 CALL qaprint(varname(1:len_trim(varname)),ixc(nixc,igrsurf(igr)%ELEM(n)),0.0_8)
408 ELSEIF (igrsurf(igr)%ELTYP(n) == 7 ) THEN
409 CALL qaprint(varname(1:len_trim(varname)),ixtg(nixtg,igrsurf(igr)%ELEM(n)),0.0_8)
410 ELSEIF (igrsurf(igr)%ELTYP(n) == 1 ) THEN
411 CALL qaprint(varname(1:len_trim(varname)),ixs(nixs,igrsurf(igr)%ELEM(n)),0.0_8)
412 ENDIF
413
414 ELSE IF(set(igs)%NB_ELLIPSE == 1)THEN
415 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_Xc = '
416 temp_double = set(igs)%ELLIPSE_XC
417 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
418 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_Yc = '
419 temp_double = set(igs)%ELLIPSE_YC
420 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
421 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_Zc = '
422 temp_double = set(igs)%ELLIPSE_ZC
423 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
424 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_a = '
425 temp_double = set(igs)%ELLIPSE_A
426 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
427 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_b = '
428 temp_double = set(igs)%ELLIPSE_B
429 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
430 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_c = '
431 temp_double = set(igs)%ELLIPSE_C
432 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
433 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_N'
434 CALL qaprint(varname(1:len_trim(varname)),set(igs)%ELLIPSE_N,0.0_8)
435 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Ellipse_Skew_ID'
436 CALL qaprint(varname(1:len_trim(varname)),set(igs)%ELLIPSE_ID_MADYMO,0.0_8)
437 ELSE IF(set(igs)%NB_PLANE == 1)THEN
438 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'Plane_XM = '
439 temp_double = set(igs)%PLANE_XM
440 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
441 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'PLANE_YM = '
442 temp_double = set(igs)%PLANE_YM
443 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
444 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'PLANE_ZM = '
445 temp_double = set(igs)%PLANE_ZM
446 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
447 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'PLANE_XM1 = '
448 temp_double = set(igs)%PLANE_XM1
449 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
450 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'PLANE_YM1 = '
451 temp_double = set(igs)%PLANE_YM1
452 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
453 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'PLANE_ZM1 = '
454 temp_double = set(igs)%PLANE_ZM1
455 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
456 ENDIF
457
458 ENDDO
459 ENDIF
460
461
462
463
464 nb_line_seg = set(igs)%HAS_LINE_SEG
465 IF (nb_line_seg > 0) THEN
466 igr = set(igs)%SET_NSLIN_ID
467 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'LINE'
468 CALL qaprint(varname(1:len_trim(varname)),igrslin(igr)%ID,0.0_8)
469 WRITE(varname,
'(A,I0,A)')
'SET_',
id,
'_'//
'LINE_NB_SEG'
470 CALL qaprint(varname(1:len_trim(varname)),nb_line_seg,0.0_8)
471 DO n = 1,nb_line_seg
472 WRITE(varname,
'(A,I0,A,I0,A)')
'SET_',
id,
'_'//
'LINE_SEG_',n,'_node_1
'
473 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),ITAB(IGRSLIN(IGR)%NODES(N,1)),0.0_8)
474 WRITE(VARNAME,'(a,i0,a,i0,a)
') 'set_',ID,'_
'//'line_seg_
',N,'_node_2
'
475 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),ITAB(IGRSLIN(IGR)%NODES(N,2)),0.0_8)
476 WRITE(VARNAME,'(a,i0,a,i0,a)
') 'set_',ID,'_
'//'line_seg_
',N,'_eltyp
'
477 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IGRSLIN(IGR)%ELTYP(N),0.0_8)
478 WRITE(VARNAME,'(a,i0,a,i0,a)
') 'set_',ID,'_
'//'line_seg_
',N,'_elem
'
479 IF (IGRSLIN(IGR)%ELTYP(N) == 3 ) THEN ! SH4N
480 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXC(NIXC,IGRSLIN(IGR)%ELEM(N)),0.0_8)
481 ELSEIF (IGRSLIN(IGR)%ELTYP(N) == 7 ) THEN ! SH3N
482 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXTG(NIXTG,IGRSLIN(IGR)%ELEM(N)),0.0_8)
483 ELSEIF (IGRSLIN(IGR)%ELTYP(N) == 1 ) THEN ! SOLID
484 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXS(NIXS,IGRSLIN(IGR)%ELEM(N)),0.0_8)
485 ELSEIF (IGRSLIN(IGR)%ELTYP(N) == 2 ) THEN ! QUAD
486 CALL QAPRINT(VARNAME(1:LEN_TRIM(VARNAME)),IXQ(NIXQ,IGRSLIN(IGR)%ELEM(N)),0.0_8)
487 ENDIF ! IF (IGRSLIN(IGR)%ELTYP(N) == 3 )
488 ENDDO
489 ENDIF ! IF (NB_LINE_SEG > 0)
490!---
491 ENDDO ! DO KK = 1, NSETS
492 ENDIF ! IF (NSETS > 0)
493 ENDIF ! IF (MYQAKEY('/set'))
494
495 RETURN
subroutine hm_get_int_array_2indexes(name, ival, index1, index2, is_available, lsubmodel)
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
subroutine hm_get_string(name, sval, size, is_available)
subroutine hm_get_string_index(name, sval, index, size, is_available)
subroutine hm_option_start(entity_type)
integer, parameter nchartitle
integer, parameter ncharkey
integer, parameter ncharfield
integer, parameter ncharline
logical function myqakey(value)
@purpose Check if a given value is part of the values set by env variable Useful to make a condition ...
subroutine qaprint(name, idin, value)
@purpose print one entry to QA extract file example of call for real print CALL QAPRINT('MY_LABEL',...