OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
h3d_nodal_vector.F File Reference
#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "intstamp_c.inc"
#include "scr03_c.inc"
#include "scr16_c.inc"
#include "param_c.inc"
#include "sphcom.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine h3d_nodal_vector (elbuf_tab, nodal_vector, ifunc, iparg, geo, mass, pm, itab, node_id, info1, info2, is_written_node, h3d_part, ipartc, iparttg, ixc, ixtg, temp, iflow, rflow, ixs, ixq, nv46, monvol, diag_sms, ms, pdama2, x, volmon, stifr, stifn, a, d, v, cont, fcontg, fint, fext, keyword, fncont, fncontg, ftcont, ftcontg, fncont2, dr, dxancg, fanreac, fcluster, mcluster, vr, fopt, npby, vgaz, ipari, igrnod, weight, nodglob, fcont_max, fncontp2, ftcontp2, ar, ipartsp, ipartr, ipartp, ipartt, iparts, ipartq, kxsp, ixr, ixp, ixt, n_h3d_part_list, nodal_vector_fvm, is_written_node_fvm, airbags_total_fvm_in_h3d, smonvol, svolmon, ispmd, fvdata_p, airbags_node_id_shift, w, sw, x_c)

Function/Subroutine Documentation

◆ h3d_nodal_vector()

subroutine h3d_nodal_vector ( type (elbuf_struct_), dimension(ngroup) elbuf_tab,
nodal_vector,
integer ifunc,
integer, dimension(nparg,*) iparg,
geo,
mass,
pm,
integer, dimension(numnod) itab,
integer, dimension(*) node_id,
integer info1,
integer info2,
integer, dimension(numnod) is_written_node,
integer, dimension(*) h3d_part,
integer, dimension(*) ipartc,
integer, dimension(*) iparttg,
integer, dimension(nixc,numelc) ixc,
integer, dimension(nixtg,numeltg) ixtg,
temp,
integer, dimension(*) iflow,
rflow,
integer, dimension(nixs,numels) ixs,
integer, dimension(nixq,numelq) ixq,
integer nv46,
integer, dimension(smonvol) monvol,
diag_sms,
ms,
pdama2,
x,
volmon,
stifr,
stifn,
a,
d,
v,
cont,
fcontg,
fint,
fext,
character(len=ncharline100) keyword,
fncont,
fncontg,
ftcont,
ftcontg,
fncont2,
dr,
dxancg,
fanreac,
fcluster,
mcluster,
vr,
fopt,
integer, dimension(nnpby,*) npby,
vgaz,
integer, dimension(npari,ninter) ipari,
type (group_), dimension(ngrnod) igrnod,
integer, dimension(*) weight,
integer, dimension(*) nodglob,
fcont_max,
fncontp2,
ftcontp2,
dimension(3,numnod), intent(in) ar,
integer, dimension(numsph), intent(in) ipartsp,
integer, dimension(numelr), intent(in) ipartr,
integer, dimension(numelp), intent(in) ipartp,
integer, dimension(numelt), intent(in) ipartt,
integer, dimension(numels), intent(in) iparts,
integer, dimension(numelq), intent(in) ipartq,
integer, dimension(nisp,numsph), intent(in) kxsp,
integer, dimension(nixr,numelr), intent(in) ixr,
integer, dimension(nixp,numelp), intent(in) ixp,
integer, dimension(nixt,numelt), intent(in) ixt,
integer, intent(in) n_h3d_part_list,
dimension(3,airbags_total_fvm_in_h3d), intent(inout) nodal_vector_fvm,
integer, dimension(airbags_total_fvm_in_h3d), intent(inout) is_written_node_fvm,
integer, intent(in) airbags_total_fvm_in_h3d,
integer, intent(in) smonvol,
integer, intent(in) svolmon,
integer, intent(in) ispmd,
type(fvbag_data), dimension(nfvbag), intent(in) fvdata_p,
integer, intent(in) airbags_node_id_shift,
dimension(3,sw/3), intent(in) w,
integer, intent(in) sw,
dimension(3,numnod) x_c )

Definition at line 40 of file h3d_nodal_vector.F.

58C-----------------------------------------------
59C D e s c r i p t i o n
60C-----------------------------------------------
61C This subroutine is storing relevant vector values requested
62C for H3D export by /H3D/NODA/... keywords
63C-----------------------------------------------
64C M o d u l e s
65C-----------------------------------------------
66 USE elbufdef_mod
67 USE groupdef_mod
68 USE outmax_mod
69 USE fvbag_mod , only : fvbag_data, nfvbag
70 USE names_and_titles_mod, ONLY: ncharline100
71 use element_mod , only : nixs,nixq,nixc,nixt,nixr,nixp,nixtg
72C-----------------------------------------------
73C I m p l i c i t T y p e s
74C-----------------------------------------------
75#include "implicit_f.inc"
76C-----------------------------------------------
77C C o m m o n B l o c k s
78C-----------------------------------------------
79#include "com01_c.inc"
80#include "com04_c.inc"
81#include "intstamp_c.inc"
82#include "scr03_c.inc"
83#include "scr16_c.inc"
84#include "param_c.inc"
85#include "sphcom.inc"
86C-----------------------------------------------
87C D u m m y A r g u m e n t s
88C-----------------------------------------------
89 INTEGER,INTENT(IN) :: SMONVOL, SVOLMON, ISPMD,AIRBAGS_NODE_ID_SHIFT, AIRBAGS_TOTAL_FVM_IN_H3D
90 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
91 INTEGER,INTENT(IN) :: SW
92 my_real
93 . nodal_vector(3*numnod),mass(*),geo(npropg,numgeo),
94 . pm(npropm,nummat),temp(*),rflow(*),volmon(svolmon), diag_sms(*),ms(numnod),
95 . pdama2(2,*),x(3,numnod),stifr(*),stifn(numnod),a(3,numnod),d(3,numnod),v(3,numnod), cont(3,*),
96 . fcontg(3,*), fint(3,numnod), fext(3,numnod),fncont(3,*),fncontg(3,*),
97 . ftcont(3,*),ftcontg(3,*),fncont2(3,*), dr(3,numnod),dxancg(3,*),
98 . fanreac(6,*),fcluster(3,*),mcluster(3,*),vr(3,numnod),fopt(6,*),vgaz(3,*),
99 . fcont_max(3,*),fncontp2(3,*),ftcontp2(3,*)
100 my_real,INTENT(IN) :: w(3,sw/3)
101 INTEGER IPARG(NPARG,*),IFUNC,NODE_ID(*),
102 . INFO1,INFO2,IS_WRITTEN_NODE(NUMNOD),H3D_PART(*),ITAB(NUMNOD),
103 . IXC(NIXC,NUMELC),IXTG(NIXTG,NUMELTG),IPARTC(*),IPARTTG(*),IFLOW(*),
104 . IXS(NIXS,NUMELS),IXQ(NIXQ,NUMELQ),NV46,MONVOL(SMONVOL),NPBY(NNPBY,*),
105 . IPARI(NPARI,NINTER),WEIGHT(*),NODGLOB(*)
106 CHARACTER(LEN=NCHARLINE100):: KEYWORD
107 TYPE (GROUP_) , DIMENSION(NGRNOD) :: IGRNOD
108 my_real , INTENT(IN) :: ar(3,numnod)
109 my_real , INTENT(IN) , DIMENSION(3,NUMNOD) :: x_c
110 INTEGER ,INTENT(IN) :: IPARTSP(NUMSPH),IPARTR(NUMELR),IPARTP(NUMELP),
111 . IPARTT(NUMELT),IPARTS(NUMELS),IPARTQ(NUMELQ)
112 INTEGER ,INTENT(IN) :: KXSP(NISP,NUMSPH),IXR(NIXR,NUMELR),IXP(NIXP,NUMELP),
113 . IXT(NIXT,NUMELT)
114 INTEGER ,INTENT(IN) :: N_H3D_PART_LIST
115 my_real, INTENT(INOUT) :: nodal_vector_fvm(3,airbags_total_fvm_in_h3d)
116 INTEGER, INTENT(INOUT) :: IS_WRITTEN_NODE_FVM(AIRBAGS_TOTAL_FVM_IN_H3D)
117 TYPE(FVBAG_DATA),INTENT(IN) :: FVDATA_P(NFVBAG)
118C-----------------------------------------------
119C L o c a l V a r i a b l e s
120C-----------------------------------------------
121 INTEGER I,J,K,N,IOK_PART(NUMNOD)
122 my_real value(3)
123C-----------------------------------------------
124
125 VALUE = zero
126
127 DO i=1,numnod
128 node_id(i) = itab(i)
129 iok_part(i) = 0
130 is_written_node(i) = 0
131 ENDDO
132
133 IF(n_h3d_part_list /= 0)THEN
134 DO i=1,numsph
135 IF ( h3d_part(ipartsp(i)) == 1) THEN
136 IF(kxsp(2,i) > 0 )iok_part(kxsp(2,i)) = 1
137 ENDIF
138 ENDDO
139
140 DO i=1,numelr
141 IF ( h3d_part(ipartr(i)) == 1) THEN
142 DO j=2,4
143 IF(ixr(j,i) > 0 )iok_part(ixr(j,i)) = 1
144 ENDDO
145 ENDIF
146 ENDDO
147
148 DO i=1,numelp
149 IF ( h3d_part(ipartp(i)) == 1) THEN
150 DO j=2,4
151 IF(ixp(j,i) > 0 )iok_part(ixp(j,i)) = 1
152 ENDDO
153 ENDIF
154 ENDDO
155
156 DO i=1,numelt
157 IF ( h3d_part(ipartt(i)) == 1) THEN
158 DO j=2,4
159 IF(ixt(j,i) > 0 )iok_part(ixt(j,i)) = 1
160 ENDDO
161 ENDIF
162 ENDDO
163
164 DO i=1,numelc
165 IF ( h3d_part(ipartc(i)) == 1) THEN
166 DO j=2,5
167 IF(ixc(j,i) > 0 )iok_part(ixc(j,i)) = 1
168 ENDDO
169 ENDIF
170 ENDDO
171
172 DO i=1,numeltg
173 IF ( h3d_part(iparttg(i)) == 1) THEN
174 DO j=2,4
175 IF(ixtg(j,i) > 0 )iok_part(ixtg(j,i)) = 1
176 ENDDO
177 ENDIF
178 ENDDO
179
180 DO i=1,numels
181 IF ( h3d_part(iparts(i)) == 1) THEN
182 DO j=2,9
183 IF(ixs(j,i) > 0 )iok_part(ixs(j,i)) = 1
184 ENDDO
185 ENDIF
186 ENDDO
187
188 DO i=1,numelq
189 IF ( h3d_part(ipartq(i)) == 1) THEN
190 DO j=2,5
191 IF(ixq(j,i) > 0 )iok_part(ixq(j,i)) = 1
192 ENDDO
193 ENDIF
194 ENDDO
195 ELSE
196 iok_part(1:numnod) = 1
197 ENDIF
198C
199C--------------------------------------------------
200 IF(keyword == 'VEL') THEN
201C--------------------------------------------------
202 DO i=1,numnod
203 value(1) = v(1,i)
204 value(2) = v(2,i)
205 value(3) = v(3,i)
206 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
207 ENDDO
208
209 !------------------------------------------
210 ! /MONVOL/FVMBAG*
211 ! set velocity from polyhedra centroids
212 !------------------------------------------
213 IF(airbags_total_fvm_in_h3d > 0)THEN
214 is_written_node_fvm(1:airbags_total_fvm_in_h3d)=0
215 IF(nfvbag>0 .AND. ispmd == 0 .AND. airbags_total_fvm_in_h3d > 0)THEN
216 CALL anim_nodal_vector_fvmbags('VEL ' , nodal_vector_fvm,monvol,volmon ,fvdata_p, nfvbag ,
217 . smonvol, svolmon ,airbags_total_fvm_in_h3d,is_written_node_fvm,
218 . airbags_node_id_shift )
219 ENDIF
220 ENDIF
221
222C--------------------------------------------------
223 ELSEIF(keyword == 'DIS') THEN
224C--------------------------------------------------
225 DO i=1,numnod
226 value(1) = d(1,i)
227 value(2) = d(2,i)
228 value(3) = d(3,i)
229 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
230 ENDDO
231C--------------------------------------------------
232 ELSEIF(keyword == 'ACC') THEN
233C--------------------------------------------------
234 DO i=1,numnod
235 value(1) = a(1,i)
236 value(2) = a(2,i)
237 value(3) = a(3,i)
238 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
239 ENDDO
240C--------------------------------------------------
241 ELSEIF(keyword == 'CONT'.AND.keyword /= 'CONT/TMAX')THEN
242C--------------------------------------------------
243 IF(nintstamp==0)THEN
244 DO i=1,numnod
245 value(1) = cont(1,i)
246 value(2) = cont(2,i)
247 value(3) = cont(3,i)
248 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
249 ENDDO
250 ELSE
251 DO i=1,numnod
252 k=nodglob(i)
253 value(1) = cont(1,i) + fcontg(1,k)
254 value(2) = cont(2,i) + fcontg(2,k)
255 value(3) = cont(3,i) + fcontg(3,k)
256 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
257 ENDDO
258 END IF
259C------------------------------------------------------------------------------
260 ELSEIF(keyword == 'CONT/TMAX')THEN
261C------------------------------------------------------------------------------
262 IF(nintstamp==0.OR.nspmd==1)THEN
263 DO i=1,numnod
264 value(1) = fcont_max(1,i)
265 value(2) = fcont_max(2,i)
266 value(3) = fcont_max(3,i)
267 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
268 ENDDO
269 ELSE
270 DO i=1,numnod
271 value(1) = zero
272 value(2) = zero
273 value(3) = zero
274 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
275 ENDDO
276 ENDIF
277C--------------------------------------------------
278 ELSEIF(keyword == 'FINT') THEN
279C--------------------------------------------------
280 DO i=1,numnod
281 value(1) = fint(1,i)
282 value(2) = fint(2,i)
283 value(3) = fint(3,i)
284 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
285 ENDDO
286C--------------------------------------------------
287 ELSEIF(keyword == 'FEXT') THEN
288C--------------------------------------------------
289 DO i=1,numnod
290 value(1) = fext(1,i)
291 value(2) = fext(2,i)
292 value(3) = fext(3,i)
293 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
294 ENDDO
295C--------------------------------------------------
296 ELSEIF(keyword == 'FOPT/FORCE') THEN
297C--------------------------------------------------
298 DO i=1,numnod
299 is_written_node(i) = 0
300 ENDDO
301 DO n=1,nrbody
302 i = npby(1,n)
303 IF (i>0) THEN
304 IF (weight(i)==1) THEN
305 value(1) = fopt(1,nsect+n)
306 value(2) = fopt(2,nsect+n)
307 value(3) = fopt(3,nsect+n)
308 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
309 ENDIF
310 ENDIF
311 ENDDO
312C--------------------------------------------------
313 ELSEIF(keyword == 'FOPT/MOMENT') THEN
314C--------------------------------------------------
315 DO i=1,numnod
316 is_written_node(i) = 0
317 ENDDO
318 DO n=1,nrbody
319 i = npby(1,n)
320 IF (i>0) THEN
321 IF (weight(i)==1) THEN
322 value(1) = fopt(4,nsect+n)
323 value(2) = fopt(5,nsect+n)
324 value(3) = fopt(6,nsect+n)
325 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
326 ENDIF
327 ENDIF
328 ENDDO
329C--------------------------------------------------
330 ELSEIF(keyword == 'VROT') THEN
331C--------------------------------------------------
332 IF(iroddl/=0) THEN
333 DO i=1,numnod
334 value(1) = vr(1,i)
335 value(2) = vr(2,i)
336 value(3) = vr(3,i)
337 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
338 ENDDO
339 ENDIF
340C--------------------------------------------------
341 ELSEIF(keyword == 'FVEL') THEN
342C--------------------------------------------------
343C--------------------------------------------------
344 ELSEIF(keyword == 'FRES') THEN
345C--------------------------------------------------
346C--------------------------------------------------
347 ELSEIF(keyword == 'PCONT/NORMAL') THEN
348C--------------------------------------------------
349 IF(nintstamp==0)THEN
350 DO i=1,numnod
351 value(1) = fncont(1,i)
352 value(2) = fncont(2,i)
353 value(3) = fncont(3,i)
354 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
355 ENDDO
356 ELSE
357 DO i=1,numnod
358 k=nodglob(i)
359 value(1) = fncont(1,i)+fncontg(1,k)
360 value(2) = fncont(2,i)+fncontg(2,k)
361 value(3) = fncont(3,i)+fncontg(3,k)
362 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
363 ENDDO
364 ENDIF
365C------------------------------------------------------------------------------
366 ELSEIF(keyword == 'MAXPCONT/NORMAL')THEN
367C------------------------------------------------------------------------------
368 IF(nintstamp==0.OR.nspmd==1)THEN
369 DO i=1,numnod
370 value(1) = fncont_max(3*(i-1)+1)
371 value(2) = fncont_max(3*(i-1)+2)
372 value(3) = fncont_max(3*(i-1)+3)
373 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
374 ENDDO
375 ELSE
376 DO i=1,numnod
377 value(1) = zero
378 value(2) = zero
379 value(3) = zero
380 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
381 ENDDO
382 ENDIF
383C--------------------------------------------------
384 ELSEIF(keyword == 'PCONT/TANGENT') THEN
385C--------------------------------------------------
386 IF(nintstamp==0)THEN
387 DO i=1,numnod
388 value(1) = ftcont(1,i)
389 value(2) = ftcont(2,i)
390 value(3) = ftcont(3,i)
391 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
392 ENDDO
393 ELSE
394 DO i=1,numnod
395 k=nodglob(i)
396 value(1) = ftcont(1,i)+ftcontg(1,k)
397 value(2) = ftcont(2,i)+ftcontg(2,k)
398 value(3) = ftcont(3,i)+ftcontg(3,k)
399 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
400 ENDDO
401 ENDIF
402C------------------------------------------------------------------------------
403 ELSEIF(keyword == 'MAXPCONT/TANGENT')THEN
404C------------------------------------------------------------------------------
405 IF(nintstamp==0.OR.nspmd==1)THEN
406 DO i=1,numnod
407 value(1) = ftcont_max(3*(i-1)+1)
408 value(2) = ftcont_max(3*(i-1)+2)
409 value(3) = ftcont_max(3*(i-1)+3)
410 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
411 ENDDO
412 ELSE
413 DO i=1,numnod
414 value(1) = zero
415 value(2) = zero
416 value(3) = zero
417 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
418 ENDDO
419 ENDIF
420C--------------------------------------------------
421 ELSEIF(keyword == 'CONT2') THEN
422C--------------------------------------------------
423 DO i=1,numnod
424 value(1) = fncont2(1,i)
425 value(2) = fncont2(2,i)
426 value(3) = fncont2(3,i)
427 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
428 ENDDO
429C--------------------------------------------------
430 ELSEIF(keyword == 'CONT2/MOMENT') THEN
431C--------------------------------------------------
432 DO i=1,numnod
433 value(1) = mcont2(1,i)
434 value(2) = mcont2(2,i)
435 value(3) = mcont2(3,i)
436 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
437 ENDDO
438C--------------------------------------------------
439 ELSEIF(keyword == 'PCONT2/NORMAL') THEN
440C--------------------------------------------------
441 DO i=1,numnod
442 value(1) = fncontp2(1,i)
443 value(2) = fncontp2(2,i)
444 value(3) = fncontp2(3,i)
445 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
446 ENDDO
447C--------------------------------------------------
448 ELSEIF(keyword == 'PCONT2/TANGENT') THEN
449C--------------------------------------------------
450 DO i=1,numnod
451 value(1) = ftcontp2(1,i)
452 value(2) = ftcontp2(2,i)
453 value(3) = ftcontp2(3,i)
454 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
455 ENDDO
456C--------------------------------------------------
457 ELSEIF(keyword == 'CONT2/TMAX') THEN
458C--------------------------------------------------
459 DO i=1,numnod
460 value(1) = fcont2_max(3*(i-1)+1)
461 value(2) = fcont2_max(3*(i-1)+2)
462 value(3) = fcont2_max(3*(i-1)+3)
463 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
464 ENDDO
465C--------------------------------------------------
466 ELSEIF(keyword == 'CONT2/TMIN') THEN
467C--------------------------------------------------
468 DO i=1,numnod
469 value(1) = fcont2_min(3*(i-1)+1)
470 value(2) = fcont2_min(3*(i-1)+2)
471 value(3) = fcont2_min(3*(i-1)+3)
472 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
473 ENDDO
474C--------------------------------------------------
475 ELSEIF(keyword == 'MAXPCONT2/NORMAL') THEN
476C--------------------------------------------------
477 DO i=1,numnod
478 value(1) = fncont2_max(3*(i-1)+1)
479 value(2) = fncont2_max(3*(i-1)+2)
480 value(3) = fncont2_max(3*(i-1)+3)
481 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
482 ENDDO
483C--------------------------------------------------
484 ELSEIF(keyword == 'MAXPCONT2/TANGENT') THEN
485C--------------------------------------------------
486 DO i=1,numnod
487 value(1) = ftcont2_max(3*(i-1)+1)
488 value(2) = ftcont2_max(3*(i-1)+2)
489 value(3) = ftcont2_max(3*(i-1)+3)
490 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
491 ENDDO
492C--------------------------------------------------
493 ELSEIF(keyword == 'MINPCONT2/NORMAL') THEN
494C--------------------------------------------------
495 DO i=1,numnod
496 value(1) = fncont2_min(3*(i-1)+1)
497 value(2) = fncont2_min(3*(i-1)+2)
498 value(3) = fncont2_min(3*(i-1)+3)
499 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
500 ENDDO
501C--------------------------------------------------
502 ELSEIF(keyword == 'MINPCONT2/TANGENT') THEN
503C--------------------------------------------------
504 DO i=1,numnod
505 value(1) = ftcont2_min(3*(i-1)+1)
506 value(2) = ftcont2_min(3*(i-1)+2)
507 value(3) = ftcont2_min(3*(i-1)+3)
508 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
509 ENDDO
510C--------------------------------------------------
511 ELSEIF(keyword == 'DROT')THEN
512C--------------------------------------------------
513 IF( (idrot == 1 .OR. isecut>0 .OR. iisrot>0 .OR. impose_dr>0) .AND. iroddl/=0) THEN
514 DO i=1,numnod
515 value(1) = dr(1,i)
516 value(2) = dr(2,i)
517 value(3) = dr(3,i)
518 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
519 ENDDO
520 ENDIF
521C--------------------------------------------------
522 ELSEIF (keyword == 'DXANC') THEN
523C--------------------------------------------------
524 DO i=1,numnod
525 value(1) = dxancg(1,i)
526 value(2) = dxancg(2,i)
527 value(3) = dxancg(3,i)
528 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
529 ENDDO
530C--------------------------------------------------
531 ELSEIF (keyword == 'GVEL') THEN
532C--------------------------------------------------
533 IF(ialelag > 0 ) THEN
534 DO i=1,numnod
535 value(1) = vgaz(1,i)
536 value(2) = vgaz(2,i)
537 value(3) = vgaz(3,i)
538 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
539 ENDDO
540 ENDIF
541C--------------------------------------------------
542 ELSEIF(keyword == 'FREAC') THEN
543C--------------------------------------------------
544 DO i=1,numnod
545 value(1)=fanreac(1,i)
546 value(2)=fanreac(2,i)
547 value(3)=fanreac(3,i)
548 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
549 ENDDO
550C--------------------------------------------------
551 ELSEIF(keyword == 'MREAC') THEN
552C--------------------------------------------------
553 DO i=1,numnod
554 value(1)=fanreac(4,i)
555 value(2)=fanreac(5,i)
556 value(3)=fanreac(6,i)
557 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
558 ENDDO
559C--------------------------------------------------
560 ELSEIF(keyword == 'CLUSTER/FORCE') THEN
561C--------------------------------------------------
562 DO i=1,numnod
563 value(1)=fcluster(1,i)
564 value(2)=fcluster(2,i)
565 value(3)=fcluster(3,i)
566 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
567 ENDDO
568C--------------------------------------------------
569 ELSEIF(keyword == 'CLUSTER/MOMENT') THEN
570C--------------------------------------------------
571 DO i=1,numnod
572 value(1)=mcluster(1,i)
573 value(2)=mcluster(2,i)
574 value(3)=mcluster(3,i)
575 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
576 ENDDO
577C--------------------------------------------------
578 ELSEIF(keyword == 'ZVEL') THEN
579 CALL h3d_velvecc22(elbuf_tab,iparg,1,ixs,ixq,itab,iok_part,is_written_node,nodal_vector)
580C--------------------------------------------------
581C--------------------------------------------------
582 ELSEIF(keyword == 'ZFVEL') THEN
583 CALL h3d_velvecz22(elbuf_tab,iparg,ipari,igrnod,x,ixs,ixq,itab,1,iok_part,is_written_node,nodal_vector)
584C--------------------------------------------------
585C--------------------------------------------------
586 ELSEIF(keyword == 'ZMOM') THEN
587 CALL h3d_velvecc22(elbuf_tab,iparg,2,ixs,ixq,itab,iok_part,is_written_node,nodal_vector)
588C--------------------------------------------------
589C--------------------------------------------------
590 ELSEIF(keyword == 'ZFP') THEN
591 CALL h3d_velvecz22(elbuf_tab,iparg,ipari,igrnod,x,ixs,ixq,itab,2,iok_part,is_written_node,nodal_vector)
592C--------------------------------------------------
593C--------------------------------------------------
594 ELSEIF(keyword == 'ZFINT') THEN
595 CALL h3d_velvecc22(elbuf_tab,iparg,3,ixs,ixq,itab,iok_part,is_written_node,nodal_vector)
596C--------------------------------------------------
597 ELSEIF(keyword == 'VEL/TMAX') THEN
598C--------------------------------------------------
599 DO i=1,numnod
600 j = 3*(i-1)+1
601 value(1) = tm_vel(j)
602 value(2) = tm_vel(j+1)
603 value(3) = tm_vel(j+2)
604 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
605 ENDDO
606C--------------------------------------------------
607 ELSEIF(keyword == 'DIS/TMAX') THEN
608C--------------------------------------------------
609 DO i=1,numnod
610 j = 3*(i-1)+1
611 value(1) = tm_dis(j)
612 value(2) = tm_dis(j+1)
613 value(3) = tm_dis(j+2)
614 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
615 ENDDO
616C--------------------------------------------------
617 ELSEIF(keyword == 'AROT') THEN
618C--------------------------------------------------
619 IF(iroddl/=0) THEN
620 DO i=1,numnod
621 value(1) = ar(1,i)
622 value(2) = ar(2,i)
623 value(3) = ar(3,i)
624 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
625 ENDDO
626 ENDIF
627C--------------------------------------------------
628 ELSEIF(keyword == 'VEL/GRID') THEN
629C--------------------------------------------------
630 IF(sw > 0) THEN
631 DO i=1,numnod
632 value(1) = w(1,i)
633 value(2) = w(2,i)
634 value(3) = w(3,i)
635 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
636 ENDDO
637 ENDIF
638C--------------------------------------------------
639 ELSEIF(keyword == 'SHELL_OFFSET') THEN
640C--------------------------------------------------
641 DO i=1,numnod
642 value(1:3) = x_c(1:3,i) - x(1:3,i)
643 CALL h3d_write_vector(iok_part,is_written_node,nodal_vector,i,0,0,VALUE)
644 ENDDO
645C--------------------------------------------------
646 ENDIF
647C--------------------------------------------------
648 RETURN
anim_nodal_vector_fvmbags
subroutine anim_nodal_vector_fvmbags(key, wa4, monvol, volmon, fvdata, nfvbag, smonvol, svolmon, airbags_total_fvm_in_h3d, is_written_node_fvm, airbags_node_id_shift)
Definition anim_nodal_vector_fvmbags.F:33
my_real
#define my_real
Definition cppsort.cpp:32
h3d_velvecc22
subroutine h3d_velvecc22(elbuf_tab, iparg, iflg, ixs, ixq, itab, iok_part, is_written_node, nodal_vector)
Definition h3d_velvecc22.F:39
h3d_velvecz22
subroutine h3d_velvecz22(elbuf_tab, iparg, ipari, igrnod, x, ixs, ixq, itab, iflg, iok_part, is_written_node, nodal_vector)
Definition h3d_velvecz22.F:40
h3d_write_vector
subroutine h3d_write_vector(iok_part, is_written, vector, i, offset, nft, value)
Definition h3d_write_vector.F:36
fvbag_mod
Definition fvbag_mod.F:44
fvbag_mod::airbags_total_fvm_in_h3d
integer airbags_total_fvm_in_h3d
Definition fvbag_mod.F:191
fvbag_mod::nfvbag
integer nfvbag
Definition fvbag_mod.F:127
groupdef_mod
Definition groupdef_mod.F:668
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::ncharline100
integer, parameter ncharline100
Definition names_and_titles_mod.F:1007
outmax_mod
Definition outmax_mod.F:52
fvbag_mod::fvbag_data
Definition fvbag_mod.F:52