36
37
38
40 USE elbufdef_mod
41 USE multi_fvm_mod
42
43
44
45#include "implicit_f.inc"
46
47
48
49#include "com01_c.inc"
50#include "param_c.inc"
51
52
53
54 INTEGER, INTENT(IN) :: NG
55 TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
56 INTEGER, INTENT(IN) :: IPARG(NPARG, *)
57 INTEGER, INTENT(IN) :: ITASK
58 TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
59 my_real,
INTENT(INOUT) :: partsav(*), gresav(*)
60 INTEGER, INTENT(IN) :: IPARTS(*), IGRTH(*), GRTH(*)
61 my_real,
INTENT(IN) :: volnew(*),timestep
62 LOGICAL, INTENT(IN) :: PRED
63
64
65
66 TYPE(G_BUFEL_), POINTER :: GBUF
67 INTEGER :: II, I
68 my_real :: rho, etot, vel2, vol, sumflux(5)
70 INTEGER :: NEL, NFT
71
72
73
74 gbuf => elbuf_tab(ng)%GBUF
75 nel = iparg(2, ng)
76 nft = iparg(3, ng)
77
78
79 DO ii = 1, nel
80 i = ii + nft
81 vii(1) = multi_fvm%RHO(i)
82 vii(2) = multi_fvm%VEL(1, i)
83 vii(3) = multi_fvm%VEL(2, i)
84 vii(4) = multi_fvm%VEL(3, i)
85 vii(5) = multi_fvm%EINT(i)
86
87 vii(1:5) = vii(1:5) / volnew(ii)
88
89 multi_fvm%RHO(i) = vii(1)
90
91 multi_fvm%VEL(1, i) = vii(2) / vii(1)
92 multi_fvm%VEL(2, i) = vii(3) / vii(1)
93 multi_fvm%VEL(3, i) = vii(4) / vii(1)
94
95 vel2 = multi_fvm%VEL(1, i)**2 + multi_fvm%VEL(2, i)**2 + multi_fvm%VEL(3, i)**2
96
97 multi_fvm%EINT(i) = vii(5) - half * vel2 * vii(1)
98 ENDDO
99
100 IF (.NOT. pred) THEN
102 . nel, nft, iparts(1 + nft), gresav, igrth, grth)
103 IF(ALLOCATED(multi_fvm%ACC))THEN
104 DO ii = 1, nel
105 i = ii + nft
106
107
108 multi_fvm%ACC(1,i) = (multi_fvm%VEL(1,i)-multi_fvm%ACC(1,i))/timestep
109 multi_fvm%ACC(2,i) = (multi_fvm%VEL(2,i)-multi_fvm%ACC(2,i))/timestep
110 multi_fvm%ACC(3,i) = (multi_fvm%VEL(3,i)-multi_fvm%ACC(3,i))/timestep
111 ENDDO
112 ENDIF
113 ENDIF
114
115
116
subroutine multi_bilan(partsav, vol, multi_fvm, nel, nft, iparts, gresav, igrth, grth)
1.15.0