32
33
34
37
38
39
40#include "implicit_f.inc"
41
42
43
44#include "com04_c.inc"
45#include "param_c.inc"
46#include "scr17_c.inc"
47#include "remesh_c.inc"
48
49
50
51 INTEGER, INTENT(IN) :: IPART(LIPART1,*), IPADMESH(KIPADMESH,*)
52 INTEGER, INTENT(IN) :: SH4TREE(KSH4TREE,*), SH3TREE(KSH3TREE,*),
53 . SH4TRIM(*),SH3TRIM(*) ,IXC(NIXC,*), IXTG(NIXTG,*)
54 my_real,
INTENT(IN) :: padmesh(kpadmesh,*)
55
56
57
58
59 INTEGER I, MY_ID, N , IP
60 CHARACTER (LEN=255) :: VARNAME
61 DOUBLE PRECISION TEMP_DOUBLE
62
63
64
66
67 IF (nadmesh > 0) THEN
68
69 WRITE(varname,'(A)') 'NADMESH'
70 CALL qaprint(varname(1:len_trim(varname)),nadmesh,0.0_8)
71
72 IF(levelmax /= 0) THEN
73 WRITE(varname,'(A)') 'LEVELMAX'
74 CALL qaprint(varname(1:len_trim(varname)),levelmax,0.0_8)
75 ENDIF
76
77 IF(iadmrule /= 0) THEN
78 WRITE(varname,'(A)') 'IADMRULE'
79 CALL qaprint(varname(1:len_trim(varname)),iadmrule,0.0_8)
80 ENDIF
81
82 IF(iadmstat /= 0) THEN
83 WRITE(varname,'(A)') 'IADMSTAT'
84 CALL qaprint(varname(1:len_trim(varname)),iadmstat,0.0_8)
85 ENDIF
86
87 IF(istatcnd /= 0) THEN
88 WRITE(varname,'(A)') 'ISTATCND'
89 CALL qaprint(varname(1:len_trim(varname)),istatcnd,0.0_8)
90 ENDIF
91
92 IF(iadmerrt /= 0) THEN
93 WRITE(varname,'(A)') 'IADMERRT'
94 CALL qaprint(varname(1:len_trim(varname)),iadmerrt,0.0_8)
95 ENDIF
96
97 IF(dtadmesh /= zero) THEN
98 WRITE(varname,'(A)') 'DTADMESH'
99 temp_double = dtadmesh
100 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
101 ENDIF
102
103 IF(tadmesh /= zero) THEN
104 WRITE(varname,'(A)') 'TADMESH'
105 temp_double = tadmesh
106 CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
107 ENDIF
108
109
110 DO ip=1,npart
111
112 my_id = ipart(4,ip)
113 IF(ipart(10,ip) /= 0) THEN
114
115 CALL qaprint(
'PART_ID_ADMESH_', my_id,0.0_8)
116
117 CALL qaprint(
'LEVELMAX',ipart(10,ip),0.0_8)
118
119 DO i=1,kipadmesh
120 WRITE(varname,'(A,I0)') 'IPADMESH_',i
121 IF(ipadmesh(i,ip)/= 0)
CALL qaprint(varname(1:len_trim(varname)),ipadmesh(i,ip),0.0_8)
122 ENDDO
123
124 DO i=1,kpadmesh
125 WRITE(varname,'(A,I0)') 'PADMESH_',i
126 temp_double = padmesh(i,ip)
127 IF(temp_double/= zero)
CALL qaprint(varname(1:len_trim(varname)),0,temp_double)
128 ENDDO
129 ENDIF
130
131 ENDDO
132
133 DO n=1,numelc
134 CALL qaprint(
'ELEM_ADMESH_4N_ID_', ixc(nixc,n),0.0_8)
135 DO i=1,ksh4tree
136 WRITE(varname,'(A,I0)') 'SH4TREE_',i
137 IF(sh4tree(i,n)/= 0)
CALL qaprint(varname(1:len_trim(varname)),sh4tree(i,n),0.0_8)
138 ENDDO
139 ENDDO
140
141 DO n=1,numeltg
142 CALL qaprint(
'ELEM_ADMESH_3N_ID_', ixtg(nixtg,n),0.0_8)
143 DO i=1,ksh3tree
144 WRITE(varname,'(A,I0)') 'SH3TREE_',i
145 IF(sh3tree(i,n)/= 0)
CALL qaprint(varname(1:len_trim(varname)),sh3tree(i,n),0.0_8)
146 ENDDO
147 ENDDO
148
149 IF(iadmstat/=0)THEN
150
151 DO n=1,lsh4trim
152 IF(sh4trim(n)/= 0) THEN
153 CALL qaprint(
'ELEM_ADMESH_STAT_4N_ID_', ixc(nixc,n),0.0_8)
154 WRITE(varname,'(A,I0)') 'SH4TRIM_',n
155 CALL qaprint(varname(1:len_trim(varname)),sh4trim(n),0.0_8)
156 ENDIF
157 ENDDO
158
159 DO n=1,lsh3trim
160 IF(sh3trim(n)/= 0) THEN
161 CALL qaprint(
'ELEM_ADMESH_STAT_3N_ID_', ixtg(nixtg,n),0.0_8)
162 WRITE(varname,'(A,I0)') 'SH3TRIM_',n
163 CALL qaprint(varname(1:len_trim(varname)),sh3trim(n),0.0_8)
164 ENDIF
165 ENDDO
166
167 ENDIF
168
169 ENDIF
170
171 ENDIF
172
173 RETURN
integer, parameter nchartitle
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',...