OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
st_qaprint_ebcs.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| st_qaprint_ebcs ../starter/source/output/qaprint/st_qaprint_ebcs.F
25!||--- called by ------------------------------------------------------
26!|| st_qaprint_driver ../starter/source/output/qaprint/st_qaprint_driver.F
27!||--- calls -----------------------------------------------------
28!||--- uses -----------------------------------------------------
29!|| restmod ../starter/share/modules1/restart_mod.F
30!||====================================================================
31 SUBROUTINE st_qaprint_ebcs(EBCS_TAB)
32C-----------------------------------------------
33C M o d u l e s
34C-----------------------------------------------
35 USE qa_out_mod
36 USE ale_ebcs_mod
37 USE restmod
38 USE ebcs_mod
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C D u m m y A r g u m e n t s
46C-----------------------------------------------
47 TYPE(t_ebcs_tab), TARGET, INTENT(IN) :: EBCS_TAB
48C--------------------------------------------------
49C L o c a l V a r i a b l e s
50C-----------------------------------------------
51 LOGICAL :: OK_QA
52 CHARACTER (LEN=255) :: VARNAME
53 INTEGER :: II, JJ
54 DOUBLE PRECISION :: TMPVAL
55 CLASS(t_ebcs), POINTER :: EBCS
56 CHARACTER(LEN=NCHARTITLE) :: TITLE
57C-----------------------------------------------
58C S o u r c e L i n e s
59C-----------------------------------------------
60 IF(nebcs == 0)RETURN
61
62 ok_qa = myqakey('/EBCS')
63
64 IF (ok_qa) THEN
65 DO ii = 1, nebcs
66 ebcs => ebcs_tab%tab(ii)%poly
67 title = ebcs%title
68 IF(len_trim(title) /= 0)THEN
69 CALL qaprint(ebcs%title,0,0.0_8)
70 ELSE
71 CALL qaprint('EBCS_FAKE_TITLE',0,0.0_8)
72 ENDIF
73 CALL qaprint('EBCS ID ',ebcs%ebcs_id,0.0_8)
74 CALL qaprint('EBCS TYPE ',ebcs%type,0.0_8)
75 CALL qaprint('EBCS SURF_ID ',ebcs%surf_id,0.0_8)
76 CALL qaprint('EBCS NB_NODE ',ebcs%nb_node,0.0_8)
77 CALL qaprint('EBCS NB_ELEM ',ebcs%nb_elem,0.0_8)
78 SELECT TYPE (twf => ebcs)
79
80 TYPE IS (t_ebcs_pres)
81 CALL qaprint('IPRES ',twf%ipres,0.0_8)
82 CALL qaprint('IRHO ',twf%irho,0.0_8)
83 CALL qaprint('IENER ',twf%iener,0.0_8)
84 CALL qaprint('IVX ',twf%ivx,0.0_8)
85 CALL qaprint('IVY ',twf%ivy,0.0_8)
86 CALL qaprint('IVZ ',twf%ivz,0.0_8)
87 tmpval = twf%c
88 CALL qaprint('C ', 0, tmpval)
89 tmpval = twf%pres
90 CALL qaprint('PRES ', 0, tmpval)
91 tmpval = twf%rho
92 CALL qaprint('RHO ', 0, tmpval)
93 tmpval = twf%ener
94 CALL qaprint('ENER ', 0, tmpval)
95 tmpval = twf%vx
96 CALL qaprint('VX ', 0, tmpval)
97 tmpval = twf%vy
98 CALL qaprint('VY ', 0, tmpval)
99 tmpval = twf%vz
100 CALL qaprint('VZ ', 0, tmpval)
101 tmpval = twf%lcar
102 CALL qaprint('LCAR ',0, tmpval)
103 tmpval = twf%r1
104 CALL qaprint('R1 ', 0, tmpval)
105 tmpval = twf%r2
106 CALL qaprint('R2 ', 0, tmpval)
107
108 TYPE IS (t_ebcs_valvin)
109 CALL qaprint('IPRES ',twf%ipres,0.0_8)
110 CALL qaprint('IRHO ',twf%irho,0.0_8)
111 CALL qaprint('IENER ',twf%iener,0.0_8)
112 tmpval = twf%c
113 CALL qaprint('C ', 0, tmpval)
114 tmpval = twf%pres
115 CALL qaprint('PRES ', 0, tmpval)
116 tmpval = twf%rho
117 CALL qaprint('RHO ', 0, tmpval)
118 tmpval = twf%ener
119 CALL qaprint('ENER ', 0, tmpval)
120 tmpval = twf%lcar
121 CALL qaprint('LCAR ', 0, tmpval)
122 tmpval = twf%r1
123 CALL qaprint('R1 ', 0, tmpval)
124 tmpval = twf%r2
125 CALL qaprint('R2 ', 0, tmpval)
126
127 TYPE IS (t_ebcs_valvout)
128 CALL qaprint('IPRES ',twf%ipres,0.0_8)
129 CALL qaprint('IRHO ',twf%irho,0.0_8)
130 CALL qaprint('IENER ',twf%iener,0.0_8)
131 tmpval = twf%c
132 CALL qaprint('C ', 0, tmpval)
133 tmpval = twf%pres
134 CALL qaprint('PRES ', 0, tmpval)
135 tmpval = twf%rho
136 CALL qaprint('RHO ', 0, tmpval)
137 tmpval = twf%ener
138 CALL qaprint('ENER ', 0, tmpval)
139 tmpval = twf%lcar
140 CALL qaprint('LCAR ', 0, tmpval)
141 tmpval = twf%r1
142 CALL qaprint('R1 ', 0, tmpval)
143 tmpval = twf%r2
144 CALL qaprint('R2 ', 0, tmpval)
145
146 TYPE IS (t_ebcs_gradp0)
147 CALL qaprint('IPRES ',twf%ipres,0.0_8)
148 CALL qaprint('IRHO ',twf%irho,0.0_8)
149 CALL qaprint('IENER ',twf%iener,0.0_8)
150 tmpval = twf%c
151 CALL qaprint('C ', 0, tmpval)
152 tmpval = twf%pres
153 CALL qaprint('PRES ', 0, tmpval)
154 tmpval = twf%rho
155 CALL qaprint('RHO ', 0, tmpval)
156 tmpval = twf%ener
157 CALL qaprint('ENER ', 0, tmpval)
158 tmpval = twf%lcar
159 CALL qaprint('LCAR ', 0, tmpval)
160 tmpval = twf%r1
161 CALL qaprint('R1 ', 0, tmpval)
162 tmpval = twf%r2
163 CALL qaprint('R2 ', 0, tmpval)
164
165 TYPE IS (t_ebcs_vel)
166 CALL qaprint('IRHO ',twf%irho,0.0_8)
167 CALL qaprint('IENER ',twf%iener,0.0_8)
168 CALL qaprint('IVX ',twf%ivx,0.0_8)
169 CALL qaprint('IVY ',twf%ivy,0.0_8)
170 CALL qaprint('IVZ ',twf%ivz,0.0_8)
171 tmpval = twf%c
172 CALL qaprint('C ', 0, tmpval)
173 tmpval = twf%rho
174 CALL qaprint('RHO ', 0, tmpval)
175 tmpval = twf%ener
176 CALL qaprint('ENER ', 0, tmpval)
177 tmpval = twf%vx
178 CALL qaprint('VX ', 0, tmpval)
179 tmpval = twf%vy
180 CALL qaprint('VY ', 0, tmpval)
181 tmpval = twf%vz
182 CALL qaprint('VZ ', 0, tmpval)
183 tmpval = twf%lcar
184 CALL qaprint('LCAR ', 0, tmpval)
185 tmpval = twf%r1
186 CALL qaprint('R1 ', 0, tmpval)
187 tmpval = twf%r2
188 CALL qaprint('R2 ', 0, tmpval)
189
190 TYPE IS (t_ebcs_normv)
191 CALL qaprint('IRHO ',twf%irho,0.0_8)
192 CALL qaprint('IENER ',twf%iener,0.0_8)
193 CALL qaprint('IVIMP ',twf%ivimp,0.0_8)
194 tmpval = twf%c
195 CALL qaprint('C ', 0, tmpval)
196 tmpval = twf%rho
197 CALL qaprint('RHO ', 0, tmpval)
198 tmpval = twf%ener
199 CALL qaprint('ENER ', 0, tmpval)
200 tmpval = twf%vimp
201 CALL qaprint('VIMP ',0, tmpval)
202 tmpval = twf%lcar
203 CALL qaprint('LCAR ', 0, tmpval)
204 tmpval = twf%r1
205 CALL qaprint('R1 ', 0, tmpval)
206 tmpval = twf%r2
207 CALL qaprint('R2 ', 0, tmpval)
208
209 TYPE IS (t_ebcs_inip)
210 tmpval = twf%c
211 CALL qaprint('C ', 0, tmpval)
212 tmpval = twf%rho
213 CALL qaprint('RHO ', 0, tmpval)
214 tmpval = twf%lcar
215 CALL qaprint('LCAR ',0, tmpval)
216
217 TYPE IS (t_ebcs_iniv)
218 tmpval = twf%c
219 CALL qaprint('C ', 0, tmpval)
220 tmpval = twf%rho
221 CALL qaprint('RHO ', 0, tmpval)
222 tmpval = twf%lcar
223 CALL qaprint('LCAR ', 0, tmpval)
224
225 TYPE IS (t_ebcs_monvol)
226 CALL qaprint('MONVOLID ',twf%monvol_id,0.0_8)
227 CALL qaprint('SENSORID ',twf%sensor_id,0.0_8)
228
229 TYPE IS (t_ebcs_inlet)
230 CALL qaprint('FORMULATION ',twf%fvm_inlet_data%formulation,0.0_8)
231 CALL qaprint('VECTOR_VELOCITY ',twf%fvm_inlet_data%vector_velocity,0.0_8)
232 DO jj = 1, 3
233 IF (twf%fvm_inlet_data%func_vel(jj) /= 0) THEN
234 WRITE(varname,'(A,I0)') 'IVEL', jj
235 CALL qaprint(trim(varname),twf%fvm_inlet_data%func_vel(jj),0.0_8)
236 ENDIF
237 ENDDO
238 DO jj = 1, 21
239 IF (twf%fvm_inlet_data%func_alpha(jj) /= 0) THEN
240 WRITE(varname,'(A,I0)') 'IALPHA', jj
241 CALL qaprint(trim(varname),twf%fvm_inlet_data%func_alpha(jj),0.0_8)
242 ENDIF
243 ENDDO
244 DO jj = 1, 21
245 IF (twf%fvm_inlet_data%func_rho(jj) /= 0) THEN
246 WRITE(varname,'(A,I0)') 'IRHO', jj
247 CALL qaprint(trim(varname),twf%fvm_inlet_data%func_rho(jj),0.0_8)
248 ENDIF
249 ENDDO
250 DO jj = 1, 21
251 IF (twf%fvm_inlet_data%func_pres(jj) /= 0) THEN
252 WRITE(varname,'(A,I0)') 'Ipres', jj
253 CALL qaprint(trim(varname),twf%fvm_inlet_data%func_pres(jj),0.0_8)
254 ENDIF
255 ENDDO
256 DO jj = 1, 3
257 tmpval = twf%fvm_inlet_data%val_vel(jj)
258 IF (twf%fvm_inlet_data%val_vel(jj) /= zero) THEN
259 WRITE(varname,'(A,I0)') 'FVEL', jj
260 CALL qaprint(trim(varname),0,tmpval)
261 ENDIF
262 ENDDO
263 DO jj = 1, 21
264 tmpval = twf%fvm_inlet_data%val_alpha(jj)
265 IF (twf%fvm_inlet_data%val_alpha(jj) /= zero) THEN
266 WRITE(varname,'(A,I0)') 'FALPHA', jj
267 CALL qaprint(trim(varname),0,tmpval)
268 ENDIF
269 ENDDO
270 DO jj = 1, 21
271 tmpval = twf%fvm_inlet_data%val_rho(jj)
272 IF (twf%fvm_inlet_data%val_rho(jj) /= zero) THEN
273 WRITE(varname,'(A,I0)') 'FRHO', jj
274 CALL qaprint(trim(varname),0,tmpval)
275 ENDIF
276 ENDDO
277 DO jj = 1, 21
278 tmpval = twf%fvm_inlet_data%val_pres(jj)
279 IF (twf%fvm_inlet_data%val_pres(jj) /= zero) THEN
280 WRITE(varname,'(A,I0)') 'Fpres', jj
281 CALL qaprint(trim(varname),0,tmpval)
282 ENDIF
283 ENDDO
284
285 TYPE IS (t_ebcs_fluxout)
286 !no parameter for this option
287
288 TYPE IS (t_ebcs_nrf)
289 tmpval = twf%tcar_p
290 CALL qaprint('TCAR_P ', 0, tmpval)
291 tmpval = twf%tcar_vf
292 CALL qaprint('TCAR_VF ', 0, tmpval)
293
294 TYPE IS (t_ebcs_propellant)
295 IF (twf%sensor_id /= 0) THEN
296 WRITE(varname,'(A)') 'sensor_id '
297 CALL qaprint(trim(varname),twf%sensor_id ,0.0_8)
298 ENDIF
299 IF (twf%submat_id /= 0) THEN
300 WRITE(varname,'(A)') 'submat_id '
301 CALL qaprint(trim(varname),twf%submat_id ,0.0_8)
302 ENDIF
303 IF (twf%ienthalpy /= 0) THEN
304 WRITE(varname,'(A)') 'ienthalpy '
305 CALL qaprint(trim(varname),twf%ienthalpy ,0.0_8)
306 ENDIF
307 tmpval = twf%a
308 CALL qaprint('a ', 0, TMPVAL)
309 TMPVAL = twf%n
310 CALL QAPRINT('n ', 0, TMPVAL)
311 TMPVAL = twf%t
312 CALL QAPRINT('tburn ', 0, TMPVAL)
313 TMPVAL = twf%rho0s
314 CALL QAPRINT('rho0s ', 0, TMPVAL)
315 !---function f
316 IF (twf%ffunc_id /= 0) THEN
317 WRITE(VARNAME,'(a)') 'ffunc_id '
318 CALL QAPRINT(TRIM(VARNAME),twf%ffunc_id ,0.0_8)
319 ENDIF
320 IF (twf%fscaleX /= 0) THEN
321 TMPVAL = twf%fscaleX
322 CALL QAPRINT('fscalex ', 0, TMPVAL)
323 ENDIF
324 IF (twf%fscaleY /= zero) THEN
325 TMPVAL = twf%fscaleY
326 CALL QAPRINT('fscaley ', 0, TMPVAL)
327 ENDIF
328 !---function g
329 IF (twf%gfunc_id /= 0) THEN
330 WRITE(VARNAME,'(a)') 'gfunc_id '
331 CALL QAPRINT(TRIM(VARNAME),twf%gfunc_id ,0.0_8)
332 ENDIF
333 IF (twf%gscaleX /= 0) THEN
334 TMPVAL = twf%gscaleX
335 CALL QAPRINT('gscalex ', 0, TMPVAL)
336 ENDIF
337 IF (twf%gscaleY /= zero) THEN
338 TMPVAL = twf%gscaleY
339 CALL QAPRINT('gscaley ', 0, TMPVAL)
340 ENDIF
341 !---function h
342 IF (twf%hfunc_id /= 0) THEN
343 WRITE(VARNAME,'(a)') 'hfunc_id '
344 CALL QAPRINT(TRIM(VARNAME),twf%hfunc_id ,0.0_8)
345 ENDIF
346 IF (twf%hscaleX /= 0) THEN
347 TMPVAL = twf%hscaleX
348 CALL QAPRINT('hscalex ', 0, TMPVAL)
349 ENDIF
350 IF (twf%hscaleY /= zero) THEN
351 TMPVAL = twf%hscaleY
352 CALL QAPRINT('hscaley ', 0, TMPVAL)
353 ENDIF
354
355 END SELECT
356
357 ENDDO
358 ENDIF
359C-----------------------------------------------
360 RETURN
361 END
362
integer nebcs
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 ...
Definition qa_out_mod.F:694
subroutine qaprint(name, idin, value)
@purpose print one entry to QA extract file example of call for real print CALL QAPRINT('MY_LABEL',...
Definition qa_out_mod.F:390
subroutine st_qaprint_ebcs(ebcs_tab)