30
31
32
33 USE multi_fvm_mod
35
36
37
38#include "implicit_f.inc"
39
40
41
42
43
44#include "com04_c.inc"
45
46
47
48 TYPE(FVM_INIVEL_STRUCT), INTENT(IN) :: FVM_INIVEL(*)
49 TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
50
51 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
52 TYPE (GROUP_) , DIMENSION(NGRQUAD) :: IGRQUAD
53 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
54
55
56
57 INTEGER I, IGRS, NOSYS, J
58
59 DO i = 1, ninvel
60 IF (fvm_inivel(i)%FLAG) THEN
61
62 IF (fvm_inivel(i)%GRBRICID /= -1) THEN
63 igrs = fvm_inivel(i)%GRBRICID
64 DO j=1,igrbric(igrs)%NENTITY
65 nosys=igrbric(igrs)%ENTITY(j)
66 multi_fvm%VEL(1, nosys) = fvm_inivel(i)%VX
67 multi_fvm%VEL(2, nosys) = fvm_inivel(i)%VY
68 multi_fvm%VEL(3, nosys) = fvm_inivel(i)%VZ
69 ENDDO
70 ENDIF
71
72 IF (fvm_inivel(i)%GRQUADID /= -1) THEN
73 igrs = fvm_inivel(i)%GRQUADID
74 DO j=1,igrquad(igrs)%NENTITY
75 nosys=igrquad(igrs)%ENTITY(j)
76 multi_fvm%VEL(1, nosys) = fvm_inivel(i)%VX
77 multi_fvm%VEL(2, nosys) = fvm_inivel(i)%VY
78 multi_fvm%VEL(3, nosys) = fvm_inivel(i)%VZ
79 ENDDO
80 ENDIF
81
82 IF (fvm_inivel(i)%GRSH3NID /= -1) THEN
83 igrs = fvm_inivel(i)%GRSH3NID
84 DO j=1,igrsh3n(igrs)%NENTITY
85 nosys=igrsh3n(igrs)%ENTITY(j)
86 multi_fvm%VEL(1, nosys) = fvm_inivel(i)%VX
87 multi_fvm%VEL(2, nosys) = fvm_inivel(i)%VY
88 multi_fvm%VEL(3, nosys) = fvm_inivel(i)%VZ
89 ENDDO
90 ENDIF
91 ENDIF
92 ENDDO