32
33
34
37
38
39
40#include "implicit_f.inc"
41
42
43
44#include "com04_c.inc"
45
46
47
48 TYPE (DRAPE_), DIMENSION(NUMELC_DRAPE + NUMELTG_DRAPE), TARGET :: DRAPE
49 TYPE (DRAPEG_) :: DRAPEG
50 TYPE (DRAPE_PLY_) ,POINTER :: DRAPE_PLY
51
52
53
54 INTEGER I,II,MY_ID,NPT,,TEMP_INT,IE_DRAPE,KK,NSLICE,IP
55 CHARACTER (LEN=255) :: VARNAME
56 DOUBLE PRECISION TEMP_DOUBLE
57
59 DO i = 1,numelc+ numeltg
60 ie_drape = drapeg%INDX(i)
61 IF(ie_drape == 0) cycle
62 npt = drape(ie_drape)%NPLY
63 IF (npt > 0) THEN
64 IF (i <= numelc) THEN
65 WRITE(varname,'(A)') 'SHELL INTERNAL NUMBER'
66 temp_int = i
67 CALL qaprint(varname(1:len_trim(varname)),temp_int,0.0_8)
68 ELSE
69 WRITE(varname,'(A)') 'SH3N INTERNAL NUMBER'
70 temp_int = i
71 CALL qaprint(varname(1:len_trim(varname)),temp_int,0.0_8)
72 ENDIF
73 DO ii = 1, npt
74 ip = drape(ie_drape)%INDX_PLY(ii)
75 IF(ip == 0) cycle
76 drape_ply => drape(ie_drape)%DRAPE_PLY(ip)
77 nslice = drape_ply%NSLICE
78 WRITE(varname,'(A,I10)') 'IDDRAPE',ip
79 temp_int = drape_ply%IPID
80 CALL qaprint(varname(1:len_trim(varname)),temp_int,0.0_8)
81 DO kk =1,nslice
82 DO j = 1,2
83 WRITE(varname,'(A,I10,I10,I10)') 'IDRAPE',ii,kk,j
84 temp_int = drape_ply%IDRAPE(kk,j)
85 CALL qaprint(varname(1:len_trim(varname)),temp_int,0.0_8)
86 ENDDO
87 DO j = 1,2
88 WRITE(varname,'(A,I10,I10,I10)') 'RDRAPE',ii,kk,j
89 temp_double = drape_ply%RDRAPE(kk,j)
91 ENDDO
92 ENDDO
93 ENDDO
94 ENDIF
95 ENDDO
96 END IF
97
98 RETURN
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',...