OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10cndv.F File Reference
#include "implicit_f.inc"
#include "com04_c.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine s10cndv (icnds10, vnd, v)
subroutine s10mvnd (icnds10, vnd, v)
subroutine s10prtv (icnds10, v, itab)
subroutine s10cndi2a (icnds10, itagnd, a)
subroutine s10cndi2a1 (icnds10, itagnd, a)
subroutine s10getvdm (icnds10, v, vnd, vmd)

Function/Subroutine Documentation

◆ s10cndi2a()

subroutine s10cndi2a ( integer, dimension(3,*) icnds10,
integer, dimension(*) itagnd,
a )

Definition at line 159 of file s10cndv.F.

160C-----------------------------------------------
161C I m p l i c i t T y p e s
162C-----------------------------------------------
163#include "implicit_f.inc"
164C-----------------------------------------------
165C C o m m o n B l o c k s
166C-----------------------------------------------
167#include "com04_c.inc"
168C-----------------------------------------------
169C D u m m y A r g u m e n t s
170C-----------------------------------------------
171 INTEGER ICNDS10(3,*),ITAGND(*)
172 my_real
173 . a(3,*)
174C-----------------------------------------------
175C L o c a l V a r i a b l e s
176C-----------------------------------------------
177 INTEGER I,ND,N1,N2
178 my_real
179 . ax,ay,az
180C======================================================================|
181C----for Nd as main of int2 only; Ad-> acce absolu
182 DO i=1,ns10e
183 nd = icnds10(1,i)
184 IF (iabs(itagnd(iabs(nd)))>ns10e) THEN
185 n1 = icnds10(2,i)
186 n2 = icnds10(3,i)
187 ax = half*(a(1,n1) + a(1,n2))
188 ay = half*(a(2,n1) + a(2,n2))
189 az = half*(a(3,n1) + a(3,n2))
190 IF (nd<0) THEN
191 nd = -nd
192 a(1,nd) = ax
193 a(2,nd) = ay
194 a(3,nd) = az
195 ELSE
196 a(1,nd) = a(1,nd) + ax
197 a(2,nd) = a(2,nd) + ay
198 a(3,nd) = a(3,nd) + az
199 END IF
200 END IF !(IABS(ITAGND(IABS(ND)))>NS10E) THEN
201 END DO
202C----6---------------------------------------------------------------7---------8
203 RETURN
my_real
#define my_real
Definition cppsort.cpp:32

◆ s10cndi2a1()

subroutine s10cndi2a1 ( integer, dimension(3,*) icnds10,
integer, dimension(*) itagnd,
a )

Definition at line 210 of file s10cndv.F.

211C-----------------------------------------------
212C I m p l i c i t T y p e s
213C-----------------------------------------------
214#include "implicit_f.inc"
215C-----------------------------------------------
216C C o m m o n B l o c k s
217C-----------------------------------------------
218#include "com04_c.inc"
219C-----------------------------------------------
220C D u m m y A r g u m e n t s
221C-----------------------------------------------
222 INTEGER ICNDS10(3,*),ITAGND(*)
223C REAL
224 my_real
225 . a(3,*)
226C-----------------------------------------------
227C L o c a l V a r i a b l e s
228C-----------------------------------------------
229 INTEGER I,ND,N1,N2
230C REAL
231 my_real
232 . ax,ay,az
233C======================================================================|
234C----for Nd as main of int2 only; Ad-> acce relative
235 DO i=1,ns10e
236 nd = icnds10(1,i)
237 IF (nd>0.AND.iabs(itagnd(iabs(nd)))>ns10e) THEN
238 n1 = icnds10(2,i)
239 n2 = icnds10(3,i)
240 ax = half*(a(1,n1) + a(1,n2))
241 ay = half*(a(2,n1) + a(2,n2))
242 az = half*(a(3,n1) + a(3,n2))
243 a(1,nd) = a(1,nd) - ax
244 a(2,nd) = a(2,nd) - ay
245 a(3,nd) = a(3,nd) - az
246 END IF !(IABS(ITAGND(IABS(ND)))>NS10E) THEN
247 END DO
248C----6---------------------------------------------------------------7---------8
249 RETURN

◆ s10cndv()

subroutine s10cndv ( integer, dimension(3,*) icnds10,
vnd,
v )

Definition at line 28 of file s10cndv.F.

29C-----------------------------------------------
30C I m p l i c i t T y p e s
31C-----------------------------------------------
32#include "implicit_f.inc"
33C-----------------------------------------------
34C C o m m o n B l o c k s
35C-----------------------------------------------
36#include "com04_c.inc"
37C-----------------------------------------------
38C D u m m y A r g u m e n t s
39C-----------------------------------------------
40 INTEGER ICNDS10(3,*)
41 my_real
42 . v(3,*),vnd(3,*)
43C-----------------------------------------------
44C L o c a l V a r i a b l e s
45C-----------------------------------------------
46 INTEGER I,ND,N1,N2
47 my_real
48 . vx,vy,vz
49C======================================================================|
50 DO i=1,ns10e
51 nd = icnds10(1,i)
52 n1 = icnds10(2,i)
53 n2 = icnds10(3,i)
54 vx = half*(v(1,n1) + v(1,n2))
55 vy = half*(v(2,n1) + v(2,n2))
56 vz = half*(v(3,n1) + v(3,n2))
57 IF (nd<0) THEN
58 nd = -nd
59 v(1,nd) = vx
60 v(2,nd) = vy
61 v(3,nd) = vz
62 ELSE
63 v(1,nd) = v(1,nd)-vnd(1,i) + vx
64 v(2,nd) = v(2,nd)-vnd(2,i) + vy
65 v(3,nd) = v(3,nd)-vnd(3,i) + vz
66 END IF
67 vnd(1,i) = vx
68 vnd(2,i) = vy
69 vnd(3,i) = vz
70 END DO
71C----6---------------------------------------------------------------7---------8
72 RETURN

◆ s10getvdm()

subroutine s10getvdm ( integer, dimension(3,*) icnds10,
v,
vnd,
vmd )

Definition at line 256 of file s10cndv.F.

257C-----------------------------------------------
258C I m p l i c i t T y p e s
259C-----------------------------------------------
260#include "implicit_f.inc"
261C-----------------------------------------------
262C C o m m o n B l o c k s
263C-----------------------------------------------
264#include "com04_c.inc"
265C-----------------------------------------------
266C D u m m y A r g u m e n t s
267C-----------------------------------------------
268 INTEGER ICNDS10(3,*)
269C REAL
270 my_real
271 . v(3,*),vnd(3,*),vmd(3,*)
272C-----------------------------------------------
273C L o c a l V a r i a b l e s
274C-----------------------------------------------
275 INTEGER I,ND,N1,N2
276C REAL
277 my_real
278 . vx,vy,vz
279C======================================================================|
280 vmd(1:3,1:numnod) = v(1:3,1:numnod)
281 DO i=1,ns10e
282 nd = iabs(icnds10(1,i))
283 vmd(1:3,nd) = v(1:3,nd)-vnd(1:3,i)
284 END DO
285C----6---------------------------------------------------------------7---------8
286 RETURN

◆ s10mvnd()

subroutine s10mvnd ( integer, dimension(3,*) icnds10,
vnd,
v )

Definition at line 79 of file s10cndv.F.

80C-----------------------------------------------
81C I m p l i c i t T y p e s
82C-----------------------------------------------
83#include "implicit_f.inc"
84C-----------------------------------------------
85C C o m m o n B l o c k s
86C-----------------------------------------------
87#include "com04_c.inc"
88C-----------------------------------------------
89C D u m m y A r g u m e n t s
90C-----------------------------------------------
91 INTEGER ICNDS10(3,*)
92 my_real
93 . v(3,*),vnd(3,*)
94C-----------------------------------------------
95C L o c a l V a r i a b l e s
96C-----------------------------------------------
97 INTEGER I,ND,N1,N2
98 my_real
99 . vx,vy,vz
100C======================================================================|
101C---relative velocity for KE compute----
102 DO i=1,ns10e
103 nd = icnds10(1,i)
104 IF (nd>0) THEN
105 v(1,nd) = v(1,nd)-vnd(1,i)
106 v(2,nd) = v(2,nd)-vnd(2,i)
107 v(3,nd) = v(3,nd)-vnd(3,i)
108 ELSEIF (nd<0) THEN
109 nd = -nd
110 v(1,nd) = zero
111 v(2,nd) = zero
112 v(3,nd) = zero
113 END IF
114 END DO
115C----6---------------------------------------------------------------7---------8
116 RETURN

◆ s10prtv()

subroutine s10prtv ( integer, dimension(3,*) icnds10,
v,
integer, dimension(*) itab )

Definition at line 121 of file s10cndv.F.

122C-----------------------------------------------
123C I m p l i c i t T y p e s
124C-----------------------------------------------
125#include "implicit_f.inc"
126C-----------------------------------------------
127C C o m m o n B l o c k s
128C-----------------------------------------------
129#include "com04_c.inc"
130#include "units_c.inc"
131C-----------------------------------------------
132C D u m m y A r g u m e n t s
133C-----------------------------------------------
134 INTEGER ICNDS10(3,*),itab(*)
135 my_real
136 . v(3,*)
137C-----------------------------------------------
138C L o c a l V a r i a b l e s
139C-----------------------------------------------
140 INTEGER I,ND,N1,N2
141C======================================================================|
142 DO i=1,ns10e
143 nd = iabs(icnds10(1,i))
144 n1 = icnds10(2,i)
145 n2 = icnds10(3,i)
146 write(iout,*)'ND,N1,N2,ND_id,I=',itab(nd),itab(n1),itab(n2),icnds10(1,i),i
147 write(iout,*)v(1,nd),v(1,n1) , v(1,n2)
148 write(iout,*)v(2,nd),v(2,n1) , v(2,n2)
149 write(iout,*)v(3,nd),v(3,n1) , v(3,n2)
150 END DO
151C----6---------------------------------------------------------------7---------8
152 RETURN