OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4cumu3.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "scr18_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine s4cumu3 (offg, e, nc1, nc2, nc3, nc4, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, them, fthe, condn, conde, nel, jthe, nodadt_therm)

Function/Subroutine Documentation

◆ s4cumu3()

subroutine s4cumu3 ( offg,
e,
integer, dimension(*) nc1,
integer, dimension(*) nc2,
integer, dimension(*) nc3,
integer, dimension(*) nc4,
stifn,
sti,
f11,
f21,
f31,
f12,
f22,
f32,
f13,
f23,
f33,
f14,
f24,
f34,
them,
fthe,
condn,
conde,
integer, intent(in) nel,
integer, intent(in) jthe,
integer, intent(in) nodadt_therm )

Definition at line 29 of file s4cumu3.F.

37C-----------------------------------------------
38C I m p l i c i t T y p e s
39C-----------------------------------------------
40#include "implicit_f.inc"
41C-----------------------------------------------
42C G l o b a l P a r a m e t e r s
43C-----------------------------------------------
44#include "mvsiz_p.inc"
45C-----------------------------------------------
46C D u m m y A r g u m e n t s
47C-----------------------------------------------
48 INTEGER, INTENT(IN) :: NEL
49 INTEGER, INTENT(IN) :: JTHE
50 INTEGER, INTENT(IN) :: NODADT_THERM
51 INTEGER NC1(*), NC2(*), NC3(*), NC4(*)
52C REAL
53 my_real
54 . offg(*),e(3,*),stifn(*),sti(*),
55 . f11(*),f21(*),f31(*),f12(*),f22(*),f32(*),
56 . f13(*),f23(*),f33(*),f14(*),f24(*),f34(*),
57 . them(mvsiz,4),fthe(*),condn(*),conde(*)
58CMasParINCLUDE 'scumu3.intmap.inc'
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "scr18_c.inc"
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER I,J
67 my_real
68 . off_l
69C-----------------------------------------------
70 off_l = 0.
71 DO i=1,nel
72 off_l = min(off_l,offg(i))
73 ENDDO
74 IF(off_l<zero)THEN
75 DO i=1,nel
76 IF(offg(i)<zero)THEN
77 f11(i)=zero
78 f21(i)=zero
79 f31(i)=zero
80 f12(i)=zero
81 f22(i)=zero
82 f32(i)=zero
83 f13(i)=zero
84 f23(i)=zero
85 f33(i)=zero
86 f14(i)=zero
87 f24(i)=zero
88 f34(i)=zero
89 sti(i)=zero
90 ENDIF
91 ENDDO
92 ENDIF
93 IF(jthe < 0 ) THEN
94 IF(off_l<=zero)THEN
95 DO j=1,4
96 DO i=1,nel
97 IF(offg(i)<=zero)THEN
98 them(i,j)=zero
99 ENDIF
100 ENDDO
101 ENDDO
102 ENDIF
103 IF(nodadt_therm == 1) THEN
104 IF(off_l<zero)THEN
105 DO i=1,nel
106 IF(offg(i)<zero)THEN
107 conde(i)=zero
108 ENDIF
109 ENDDO
110 ENDIF
111 ENDIF
112 ENDIF
113C-----------------------------------------------
114C
115C because 2*Mnodal = 2*Melement/4
116 DO i=1,nel
117 sti(i)=half*sti(i)
118 END DO
119 IF(nodadt_therm == 1 ) THEN
120 DO i=1,nel
121 conde(i)=fourth*conde(i)
122 END DO
123 ENDIF
124C
125 IF(jthe < 0 ) THEN
126 IF(nodadt_therm == 1) THEN
127 DO i=1,nel
128C
129 e(1,nc1(i))=e(1,nc1(i))+f11(i)
130 e(2,nc1(i))=e(2,nc1(i))+f21(i)
131 e(3,nc1(i))=e(3,nc1(i))+f31(i)
132 stifn(nc1(i))=stifn(nc1(i))+sti(i)
133 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
134 condn(nc1(i))= condn(nc1(i))+conde(i)
135C
136 e(1,nc2(i))=e(1,nc2(i))+f12(i)
137 e(2,nc2(i))=e(2,nc2(i))+f22(i)
138 e(3,nc2(i))=e(3,nc2(i))+f32(i)
139 stifn(nc2(i))=stifn(nc2(i))+sti(i)
140 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
141 condn(nc2(i))= condn(nc2(i))+ conde(i)
142C
143 e(1,nc3(i))=e(1,nc3(i))+f13(i)
144 e(2,nc3(i))=e(2,nc3(i))+f23(i)
145 e(3,nc3(i))=e(3,nc3(i))+f33(i)
146 stifn(nc3(i))=stifn(nc3(i))+sti(i)
147 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
148 condn(nc3(i))= condn(nc3(i))+ conde(i)
149C
150 e(1,nc4(i))=e(1,nc4(i))+f14(i)
151 e(2,nc4(i))=e(2,nc4(i))+f24(i)
152 e(3,nc4(i))=e(3,nc4(i))+f34(i)
153 stifn(nc4(i))=stifn(nc4(i))+sti(i)
154 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
155 condn(nc4(i))= condn(nc4(i))+ conde(i)
156C
157 ENDDO
158 ELSE
159 DO i=1,nel
160C
161 e(1,nc1(i))=e(1,nc1(i))+f11(i)
162 e(2,nc1(i))=e(2,nc1(i))+f21(i)
163 e(3,nc1(i))=e(3,nc1(i))+f31(i)
164 stifn(nc1(i))=stifn(nc1(i))+sti(i)
165 fthe(nc1(i)) = fthe(nc1(i)) + them(i,1)
166C
167 e(1,nc2(i))=e(1,nc2(i))+f12(i)
168 e(2,nc2(i))=e(2,nc2(i))+f22(i)
169 e(3,nc2(i))=e(3,nc2(i))+f32(i)
170 stifn(nc2(i))=stifn(nc2(i))+sti(i)
171 fthe(nc2(i)) = fthe(nc2(i)) + them(i,2)
172C
173 e(1,nc3(i))=e(1,nc3(i))+f13(i)
174 e(2,nc3(i))=e(2,nc3(i))+f23(i)
175 e(3,nc3(i))=e(3,nc3(i))+f33(i)
176 stifn(nc3(i))=stifn(nc3(i))+sti(i)
177 fthe(nc3(i)) = fthe(nc3(i)) + them(i,3)
178C
179 e(1,nc4(i))=e(1,nc4(i))+f14(i)
180 e(2,nc4(i))=e(2,nc4(i))+f24(i)
181 e(3,nc4(i))=e(3,nc4(i))+f34(i)
182 stifn(nc4(i))=stifn(nc4(i))+sti(i)
183 fthe(nc4(i)) = fthe(nc4(i)) + them(i,4)
184C
185 ENDDO
186 ENDIF
187 ELSE
188 DO 100 i=1,nel
189C
190 e(1,nc1(i))=e(1,nc1(i))+f11(i)
191 e(2,nc1(i))=e(2,nc1(i))+f21(i)
192 e(3,nc1(i))=e(3,nc1(i))+f31(i)
193 stifn(nc1(i))=stifn(nc1(i))+sti(i)
194C
195 e(1,nc2(i))=e(1,nc2(i))+f12(i)
196 e(2,nc2(i))=e(2,nc2(i))+f22(i)
197 e(3,nc2(i))=e(3,nc2(i))+f32(i)
198 stifn(nc2(i))=stifn(nc2(i))+sti(i)
199C
200 e(1,nc3(i))=e(1,nc3(i))+f13(i)
201 e(2,nc3(i))=e(2,nc3(i))+f23(i)
202 e(3,nc3(i))=e(3,nc3(i))+f33(i)
203 stifn(nc3(i))=stifn(nc3(i))+sti(i)
204C
205 e(1,nc4(i))=e(1,nc4(i))+f14(i)
206 e(2,nc4(i))=e(2,nc4(i))+f24(i)
207 e(3,nc4(i))=e(3,nc4(i))+f34(i)
208 stifn(nc4(i))=stifn(nc4(i))+sti(i)
209C
210 100 CONTINUE
211
212 ENDIF
213
214 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20