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