OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
multi_face_data_elem.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!|| multi_face_elem_data ../engine/source/multifluid/multi_face_data_elem.F
25!||--- called by ------------------------------------------------------
26!|| alemain ../engine/source/ale/alemain.F
27!|| multi_timeevolution ../engine/source/multifluid/multi_timeevolution.F
28!||--- calls -----------------------------------------------------
29!|| centroid2 ../engine/source/multifluid/centroid.F
30!|| centroid2t ../engine/source/multifluid/centroid.F
31!|| centroid3 ../engine/source/multifluid/centroid.F
32!|| centroid3t ../engine/source/multifluid/centroid.F
33!|| qnorm2 ../engine/source/multifluid/qnorm2.F
34!|| qnorm2t ../engine/source/multifluid/qnorm2t.F
35!|| snorm3 ../engine/source/multifluid/snorm3.F
36!|| snorm3t ../engine/source/multifluid/snorm3t.F
37!||--- uses -----------------------------------------------------
38!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
39!||====================================================================
40 SUBROUTINE multi_face_elem_data(ITASK, IPARG, IXS, IXQ, IXTG, XGRID, WGRID, MULTI_FVM)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE multi_fvm_mod
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49C-----------------------------------------------
50C C o m m o n B l o c k s
51C-----------------------------------------------
52C NIXS, NIXQ, NIXTG, NPARG
53#include "param_c.inc"
54C NGROUP
55#include "com01_c.inc"
56C NTHREAD
57#include "task_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER, INTENT(IN) :: ITASK, IPARG(NPARG, *), IXS(NIXS, *), IXQ(NIXQ, *), IXTG(NIXTG, *)
62 my_real, INTENT(IN) :: xgrid(3, *), wgrid(3, *)
63 TYPE(multi_fvm_struct), INTENT(INOUT) :: MULTI_FVM
64C-----------------------------------------------
65C L o c a l V a r i a b l e s
66C-----------------------------------------------
67 INTEGER :: NG, MTN, NEL, ITY, ISOLNOD, NFT, JALE
68 LOGICAL :: COMPUTE_CENTROID
69
70 compute_centroid = (multi_fvm%MUSCL > 0) .OR. multi_fvm%NS_DIFF
71
72 DO ng = itask + 1, ngroup, nthread
73 mtn = iparg(1, ng)
74 IF (mtn == 151) THEN
75 nel = iparg(2, ng)
76 nft = iparg(3, ng)
77 ity = iparg(5, ng)
78 jale = iparg(7, ng)
79 isolnod = iparg(28, ng)
80 IF (ity == 1) THEN
81C Solid
82 IF (isolnod == 4) THEN
83C TETRA
84 CALL snorm3t(nel, nft, jale, ixs, xgrid, wgrid,
85 . multi_fvm%FACE_DATA%NORMAL(1:3, 1:6, 1 + nft : nel + nft),
86 . multi_fvm%FACE_DATA%WFAC(1:3, 1:6, 1 + nft : nel + nft),
87 . multi_fvm%FACE_DATA%SURF(1:6, 1 + nft : nel + nft))
88 IF (compute_centroid) THEN
89 CALL centroid3t(nel, multi_fvm%NELEM, nft, ixs, xgrid,
90 . multi_fvm%ELEM_DATA%CENTROID(:, 1 + nft : nel + nft),
91 . multi_fvm%FACE_DATA%CENTROID(:, :, 1 + nft : nel + nft))
92 ENDIF
93 ELSE
94C OTHER SOLIDS
95 CALL snorm3(nel, nft, jale, ixs, xgrid, wgrid,
96 . multi_fvm%FACE_DATA%NORMAL(1:3, 1:6, 1 + nft : nel + nft),
97 . multi_fvm%FACE_DATA%WFAC(1:3, 1:6, 1 + nft : nel + nft),
98 . multi_fvm%FACE_DATA%SURF(1:6, 1 + nft : nel + nft))
99C
100 IF (compute_centroid) THEN
101 CALL centroid3(nel, multi_fvm%NELEM, nft, ixs, xgrid,
102 . multi_fvm%ELEM_DATA%CENTROID(:, 1 + nft : nel + nft),
103 . multi_fvm%FACE_DATA%CENTROID(:, :, 1 + nft : nel + nft))
104 ENDIF
105 ENDIF
106 ELSE IF (ity == 2) THEN
107C QUAD
108 CALL qnorm2(nel, nft, jale, multi_fvm%SYM, ixq, xgrid, wgrid,
109 . multi_fvm%FACE_DATA%NORMAL(1:3, 1:4, 1 + nft : nel + nft),
110 . multi_fvm%FACE_DATA%WFAC(1:3, 1:4, 1 + nft : nel + nft),
111 . multi_fvm%FACE_DATA%SURF(1:4, 1 + nft : nel + nft))
112 IF (compute_centroid) THEN
113 CALL centroid2(nel, multi_fvm%NELEM, nft, ixq, xgrid,
114 . multi_fvm%ELEM_DATA%CENTROID(:, 1 + nft : nel + nft),
115 . multi_fvm%FACE_DATA%CENTROID(:, :, 1 + nft : nel + nft))
116 ENDIF
117 ELSE IF (ity == 7) THEN
118C TRIANGLE
119 CALL qnorm2t(nel, nft, jale, multi_fvm%SYM, ixtg, xgrid, wgrid,
120 . multi_fvm%FACE_DATA%NORMAL(1:3, 1:3, 1 + nft : nel + nft),
121 . multi_fvm%FACE_DATA%WFAC(1:3, 1:3, 1 + nft : nel + nft),
122 . multi_fvm%FACE_DATA%SURF(1:3, 1 + nft : nel + nft))
123 IF (compute_centroid) THEN
124 CALL centroid2t(nel, multi_fvm%NELEM, nft, ixtg, xgrid,
125 . multi_fvm%ELEM_DATA%CENTROID(:, 1 + nft : nel + nft),
126 . multi_fvm%FACE_DATA%CENTROID(:, :, 1 + nft : nel + nft))
127 ENDIF
128 ENDIF
129 ENDIF
130 ENDDO
131 END SUBROUTINE multi_face_elem_data
centroid2t
subroutine centroid2t(nel, length, nft, ixtg, xgrid, elem_centroid, face_centroid)
Definition centroid.F:247
centroid3t
subroutine centroid3t(nel, length, nft, ixs, xgrid, elem_centroid, face_centroid)
Definition centroid.F:104
centroid3
subroutine centroid3(nel, length, nft, ixs, xgrid, elem_centroid, face_centroid)
Definition centroid.F:31
centroid2
subroutine centroid2(nel, length, nft, ixq, xgrid, elem_centroid, face_centroid)
Definition centroid.F:189
my_real
#define my_real
Definition cppsort.cpp:32
multi_face_elem_data
subroutine multi_face_elem_data(itask, iparg, ixs, ixq, ixtg, xgrid, wgrid, multi_fvm)
Definition multi_face_data_elem.F:41
qnorm2
subroutine qnorm2(nel, nft, jale, sym, ixq, xgrid, wgrid, norm, wfac, surf)
Definition qnorm2.F:30
qnorm2t
subroutine qnorm2t(nel, nft, jale, sym, ixtg, xgrid, wgrid, norm, wfac, surf)
Definition qnorm2t.F:32
snorm3
subroutine snorm3(nel, nft, jale, ixs, xgrid, wgrid, norm, wfac, surf)
Definition snorm3.F:31
snorm3t
subroutine snorm3t(nel, nft, jale, ixs, xgrid, wgrid, norm, wfac, surf)
Definition snorm3t.F:33