OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8sforc3.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!|| s8sforc3 ../engine/source/elements/solid/solide8s/s8sforc3.F
25!||--- called by ------------------------------------------------------
26!|| forint ../engine/source/elements/forint.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| basis8 ../engine/source/elements/solid/solide8/basis8.F
30!|| csmall3 ../engine/source/elements/solid/solide/csmall3.F
31!|| degenes8 ../engine/source/elements/solid/solide/degenes8.F
32!|| mmain ../engine/source/materials/mat_share/mmain.F90
33!|| s8edefoc3 ../engine/source/elements/solid/solide8e/s8edefoc3.F
34!|| s8edefw3 ../engine/source/elements/solid/solide8e/s8edefw3.F
35!|| s8ederi_2 ../engine/source/elements/solid/solide8e/s8ederi_2.F
36!|| s8ederic3 ../engine/source/elements/solid/solide8e/s8ederic3.F
37!|| s8edericm3 ../engine/source/elements/solid/solide8e/s8edericm3.F
38!|| s8ederig3 ../engine/source/elements/solid/solide8e/s8ederig3.F
39!|| s8ederipr3 ../engine/source/elements/solid/solide8e/s8ederipr3.F
40!|| s8efmoy3 ../engine/source/elements/solid/solide8e/s8efmoy3.F
41!|| s8ejacip3 ../engine/source/elements/solid/solide8e/s8ejacip3.F
42!|| s8eprst_ini ../engine/source/elements/solid/solide8e/s8eprst_ini.f
43!|| s8sbdefo3 ../engine/source/elements/solid/solide8s/s8sbdefo3.F
44!|| s8sderi3 ../engine/source/elements/solid/solide8s/s8sderi3.F
45!|| s8sfint3 ../engine/source/elements/solid/solide8s/s8sfint3.F
46!|| s8sfint3_crimp ../engine/source/elements/solid/solide8s/s8sfint3_crimp.F
47!|| s8xref_imp ../engine/source/elements/solid/solide8s/s8xref_imp.F
48!|| s8zderims3 ../engine/source/elements/solid/solide8z/s8zderims3.F
49!|| s8zfint_reg ../engine/source/elements/solid/solide8z/s8zfint_reg.F
50!|| s8zfintp3 ../engine/source/elements/solid/solide8z/s8zfintp3.F
51!|| s8zzero3 ../engine/source/elements/solid/solide8z/s8zzero3.f
52!|| sbilan ../engine/source/elements/solid/solide/sbilan.F
53!|| scumu3 ../engine/source/elements/solid/solide/scumu3.F
54!|| scumu3p ../engine/source/elements/solid/solide/scumu3p.F
55!|| sdlen_dege ../engine/source/elements/solid/solide/sdlen_dege.F
56!|| sfillopt ../engine/source/elements/solid/solide/sfillopt.F
57!|| smallb3 ../engine/source/elements/solid/solide/smallb3.F
58!|| smallg3 ../engine/source/elements/solid/solide/smallg3.F
59!|| srbilan ../engine/source/elements/solid/solide/srbilan.F
60!|| srcoor3_imp ../engine/source/elements/solid/solide8s/srcoor3_imp.F
61!|| srho3 ../engine/source/elements/solid/solide/srho3.F
62!|| sstra3 ../engine/source/elements/solid/solide/sstra3.F
63!|| startime ../engine/source/system/timer_mod.F90
64!|| stoptime ../engine/source/system/timer_mod.F90
65!||--- uses -----------------------------------------------------
66!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
67!|| dt_mod ../engine/source/modules/dt_mod.F
68!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
69!|| element_mod ../common_source/modules/elements/element_mod.F90
70!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
71!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
72!|| message_mod ../engine/share/message_module/message_mod.F
73!|| mmain_mod ../engine/source/materials/mat_share/mmain.F90
74!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.F
75!|| output_mod ../common_source/modules/output/output_mod.F90
76!|| sensor_mod ../common_source/modules/sensor_mod.F90
77!|| table_mod ../engine/share/modules/table_mod.F
78!|| timer_mod ../engine/source/system/timer_mod.F90
79!||====================================================================
80 SUBROUTINE s8sforc3(TIMERS, OUTPUT, ELBUF_TAB,NG ,
81 1 PM ,GEO ,IXS ,X ,
82 2 A ,V ,MS ,W ,FLUX ,
83 3 FLU1 ,VEUL ,FV ,ALE_CONNECT ,IPARG ,
84 4 TF ,NPF ,BUFMAT ,PARTSAV ,NLOC_DMG,
85 5 DT2T ,NELTST ,ITYPTST ,STIFN ,FSKY ,
86 6 IADS ,OFFSET ,EANI ,IPARTS ,ICP ,
87 7 F11 ,F21 ,F31 ,F12 ,F22 ,
88 8 F32 ,F13 ,F23 ,F33 ,F14 ,
89 9 F24 ,F34 ,F15 ,F25 ,F35 ,
90 A F16 ,F26 ,F36 ,F17 ,F27 ,
91 B F37 ,F18 ,F28 ,F38 ,NEL ,
92 C NVC ,IPM ,ITASK ,ISTRAIN ,
93 D TEMP ,FTHE ,FTHESKY ,IEXPAN ,GRESAV ,
94 E GRTH ,IGRTH ,MSSA ,DMELS ,TABLE ,
95 F IGEO ,XDP ,VOLN ,CONDN ,CONDNSKY,
96 G D ,IOUTPRT,MAT_ELEM,H3D_STRAIN,SNPC,STF,
97 H SBUFMAT ,SVIS,NSVOIS,IDTMINS,IRESP,IDTMIN,MAXFUNC,
98 I USERL_AVAIL,DT ,GLOB_THERM,SENSORS)
99C-----------------------------------------------
100C M o d u l e s
101C-----------------------------------------------
102 USE timer_mod
103 USE output_mod, only : output_
104 USE mmain_mod
105 USE table_mod
106 USE mat_elem_mod
107 USE message_mod
108 USE nlocal_reg_mod
109 USE ale_connectivity_mod
110 USE elbufdef_mod
111 USE dt_mod
112 use glob_therm_mod
113 USE sensor_mod
114 use element_mod , only : nixs
115C-----------------------------------------------
116C I m p l i c i t T y p e s
117C-----------------------------------------------
118#include "implicit_f.inc"
119C-----------------------------------------------
120C G l o b a l P a r a m e t e r s
121C-----------------------------------------------
122#include "mvsiz_p.inc"
123C-----------------------------------------------
124C C o m m o n B l o c k s
125C-----------------------------------------------
126#include "com01_c.inc"
127#include "com04_c.inc"
128#include "com08_c.inc"
129#include "scr07_c.inc"
130#include "vect01_c.inc"
131#include "parit_c.inc"
132#include "param_c.inc"
133#include "timeri_c.inc"
134#include "impl1_c.inc"
135#include "scr17_c.inc"
136C-----------------------------------------------
137C D u m m y A r g u m e n t s
138C-----------------------------------------------
139 TYPE(timer_) , INTENT(INOUT) :: TIMERS
140 TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
141 INTEGER, INTENT(IN) :: SNPC
142 INTEGER, INTENT(IN) :: STF
143 INTEGER, INTENT(IN) :: SBUFMAT
144 INTEGER, INTENT(IN) :: NSVOIS
145 INTEGER, INTENT(IN) :: IDTMINS
146 INTEGER ,INTENT(IN) :: IRESP
147 integer,dimension(102) :: IDTMIN
148 INTEGER ,INTENT(IN) :: MAXFUNC
149 INTEGER, INTENT(IN) :: USERL_AVAIL
150
151 INTEGER IXS(NIXS,*), IPARG(NPARG,NGROUP),NPF(*),
152 . IADS(8,*),IPARTS(*), IPM(NPROPMI,*),GRTH(*),IGRTH(*),IGEO(*),
153 . IOUTPRT
154C
155 INTEGER NELTST,ITYPTST,OFFSET,NEL,ICP,
156 . NVC ,ITASK, ISTRAIN, IEXPAN,NG,H3D_STRAIN
157
158 double precision
159 . xdp(3,*)
160 my_real
161 . dt2t
162C
163 my_real
164 . pm(npropm,*), geo(*), x(*), a(*), v(*), ms(*), w(*),
165 . flux(6,*),flu1(*), veul(*), fv(*), tf(*),
166 . partsav(*),stifn(*), fsky(*),eani(*),bufmat(*),
167 . f11(mvsiz),f21(mvsiz),f31(mvsiz),
168 . f12(mvsiz),f22(mvsiz),f32(mvsiz),
169 . f13(mvsiz),f23(mvsiz),f33(mvsiz),
170 . f14(mvsiz),f24(mvsiz),f34(mvsiz),
171 . f15(mvsiz),f25(mvsiz),f35(mvsiz),
172 . f16(mvsiz),f26(mvsiz),f36(mvsiz),
173 . f17(mvsiz),f27(mvsiz),f37(mvsiz),
174 . f18(mvsiz),f28(mvsiz),f38(mvsiz),
175 . temp(*),fthe(*),fthesky(*),gresav(*), mssa(*), dmels(*), voln(mvsiz),
176 . condn(*),condnsky(*),d(*)
177 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
178 TYPE(ttable) TABLE(*)
179 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
180 TYPE (NLOCAL_STR_) , TARGET :: NLOC_DMG
181 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
182 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
183 TYPE (DT_), INTENT(IN) :: DT
184 type (glob_therm_) ,intent(inout) :: glob_therm
185 type (sensors_),INTENT(INOUT) :: SENSORS
186C-----------------------------------------------
187C L o c a l V a r i a b l e s
188C-----------------------------------------------
189 INTEGER NF1, IFLAG, I,
190 . ILAY,IP,IR,IS,IT,NPTR,NPTS,NPTT,j
191 INTEGER IBID, IBIDV(1), IFVM22
192 INTEGER MXT(MVSIZ),NGL(MVSIZ),NGEO(MVSIZ)
193 my_real
194 . VD2(MVSIZ) , DVOL(MVSIZ),DELTAX(MVSIZ),
195 . vis(mvsiz) , qvis(mvsiz), cxx(mvsiz) ,
196 . s1(mvsiz) , s2(mvsiz) , s3(mvsiz) ,
197 . s4(mvsiz) , s5(mvsiz) , s6(mvsiz) ,
198 . dxx(mvsiz) , dyy(mvsiz) , dzz(mvsiz) ,
199 . d4(mvsiz) , d5(mvsiz) , d6(mvsiz) ,
200 . ajc1(mvsiz) , ajc2(mvsiz) , ajc3(mvsiz) ,
201 . ajc4(mvsiz) , ajc5(mvsiz) , ajc6(mvsiz) ,
202 . ajc7(mvsiz) , ajc8(mvsiz) , ajc9(mvsiz) ,
203 . aj1(mvsiz) , aj2(mvsiz) , aj3(mvsiz) ,
204 . aj4(mvsiz) , aj5(mvsiz) , aj6(mvsiz) ,
205 . ajp1(mvsiz,8) , ajp2(mvsiz,8) , ajp3(mvsiz,8) ,
206 . ajp4(mvsiz,8) , ajp5(mvsiz,8) , ajp6(mvsiz,8) ,
207 . ajp7(mvsiz,8) , ajp8(mvsiz,8) , ajp9(mvsiz,8) ,
208 . vdx(mvsiz) , vdy(mvsiz) , vdz(mvsiz),ssp_eq(mvsiz),aire(mvsiz),
209 . fac(mvsiz) ,pr(8,8),ps(8,8),pt(8,8)
210C-----
211C variables used as arguments by the materials.
212 my_real
213 . sti(mvsiz),stin(mvsiz),gama(mvsiz,6),
214 . wxx(mvsiz) , wyy(mvsiz) , wzz(mvsiz),
215 . wxx0(mvsiz) , wyy0(mvsiz) , wzz0(mvsiz),
216 . conde(mvsiz), conden(mvsiz),divde(mvsiz)
217C Variables used in argument by materials if SPH only.
218 my_real
219 . muvoid(mvsiz)
220 INTEGER IOFFS
221 my_real
222 . OFFS(MVSIZ),DSV(MVSIZ),SDV(MVSIZ)
223 INTEGER SZ_R1_FREE
224C-----
225C variables used only in solid routines (as arguments).
226 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ),
227 . NC5(MVSIZ), NC6(MVSIZ), NC7(MVSIZ), NC8(MVSIZ)
228
229 DOUBLE PRECISION
230 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
231 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
232 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
233 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
234 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
235 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),
236 . x0(mvsiz,8),y0(mvsiz,8),z0(mvsiz,8)
237
238 my_real
239 . off(mvsiz) ,offg(mvsiz) , rhoo(mvsiz),
240 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
241 . x5(mvsiz), x6(mvsiz), x7(mvsiz), x8(mvsiz),
242 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
243 . y5(mvsiz), y6(mvsiz), y7(mvsiz), y8(mvsiz),
244 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
245 . z5(mvsiz), z6(mvsiz), z7(mvsiz), z8(mvsiz),
246 . vx1(mvsiz),vx2(mvsiz),vx3(mvsiz),vx4(mvsiz),
247 . vx5(mvsiz),vx6(mvsiz),vx7(mvsiz),vx8(mvsiz),
248 . vy1(mvsiz),vy2(mvsiz),vy3(mvsiz),vy4(mvsiz),
249 . vy5(mvsiz),vy6(mvsiz),vy7(mvsiz),vy8(mvsiz),
250 . vz1(mvsiz),vz2(mvsiz),vz3(mvsiz),vz4(mvsiz),
251 . vz5(mvsiz),vz6(mvsiz),vz7(mvsiz),vz8(mvsiz),
252 . pxc1(mvsiz),pxc2(mvsiz),pxc3(mvsiz),pxc4(mvsiz),
253 . pyc1(mvsiz),pyc2(mvsiz),pyc3(mvsiz),pyc4(mvsiz),
254 . pzc1(mvsiz),pzc2(mvsiz),pzc3(mvsiz),pzc4(mvsiz),
255 .
256 .
257 .
258 .
259 .
260 .
261 . vgxa(mvsiz),vgya(mvsiz),vgza(mvsiz), vga2(mvsiz),
262 . hx(mvsiz,4), hy(mvsiz,4), hz(mvsiz,4),
263 . xgxa(mvsiz),xgya(mvsiz),xgza(mvsiz),
264 . xgxya(mvsiz),xgyza(mvsiz),xgzxa(mvsiz),
265 . xgxa2(mvsiz),xgya2(mvsiz),xgza2(mvsiz)
266 my_real
267 . dxy(mvsiz),dyx(mvsiz),
268 . dyz(mvsiz),dzy(mvsiz),
269 . dzx(mvsiz),dxz(mvsiz)
270 my_real
271 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
272 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
273 . r31(mvsiz),r32(mvsiz),r33(mvsiz)
274 my_real
275 . wi,smax(mvsiz),volg(mvsiz),pp(mvsiz),bid(mvsiz)
276 my_real
277 . sigy(mvsiz), et(mvsiz), nu(mvsiz),
278 . r1_free(mvsiz),r3_free(mvsiz)
279 my_real
280 .
281 . tempel(mvsiz),them(mvsiz,8),die(mvsiz)
282 INTEGER NNPT,IDEG(MVSIZ)
283 PARAMETER (NNPT = 8)
284 my_real
285 . mfxx(mvsiz,nnpt),mfxy(mvsiz,nnpt),mfyx(mvsiz,nnpt),
286 . mfyy(mvsiz,nnpt),mfyz(mvsiz,nnpt),mfzy(mvsiz,nnpt),
287 . mfzz(mvsiz,nnpt),mfzx(mvsiz,nnpt),mfxz(mvsiz,nnpt),
288 .
289 .
290 . jfac(mvsiz,nnpt)
291C-----
292 TYPE(g_bufel_) ,POINTER :: GBUF
293 TYPE(L_BUFEL_) ,POINTER :: LBUF
294c-----------------------------------------------------
295 my_real
296 . bxy1(mvsiz),bxy2(mvsiz),bxy3(mvsiz),bxy4(mvsiz),
297 . bxy5(mvsiz),bxy6(mvsiz),bxy7(mvsiz),bxy8(mvsiz),
298 . byx1(mvsiz),byx2(mvsiz),byx3(mvsiz),byx4(mvsiz),
299 . byx5(mvsiz),byx6(mvsiz),byx7(mvsiz),byx8(mvsiz),
300 . bxz1(mvsiz),bxz2(mvsiz),bxz3(mvsiz),bxz4(mvsiz),
301 . bxz5(mvsiz),bxz6(mvsiz),bxz7(mvsiz),bxz8(mvsiz),
302 . bzx1(mvsiz),bzx2(mvsiz),bzx3(mvsiz),bzx4(mvsiz),
303 . bzx5(mvsiz),bzx6(mvsiz),bzx7(mvsiz),bzx8(mvsiz),
304 . byz1(mvsiz),byz2(mvsiz),byz3(mvsiz),byz4(mvsiz),
305 . byz5(mvsiz),byz6(mvsiz),byz7(mvsiz),byz8(mvsiz),
306 . bzy1(mvsiz),bzy2(mvsiz),bzy3(mvsiz),bzy4(mvsiz),
307 . bzy5(mvsiz),bzy6(mvsiz),bzy7(mvsiz),bzy8(mvsiz),
308 . bxx1(mvsiz),bxx2(mvsiz),bxx3(mvsiz),bxx4(mvsiz),
309 . bxx5(mvsiz),bxx6(mvsiz),bxx7(mvsiz),bxx8(mvsiz),
310 . byy1(mvsiz),byy2(mvsiz),byy3(mvsiz),byy4(mvsiz),
311 . byy5(mvsiz),byy6(mvsiz),byy7(mvsiz),byy8(mvsiz),
312 . bzz1(mvsiz),bzz2(mvsiz),bzz3(mvsiz),bzz4(mvsiz),
313 . bzz5(mvsiz),bzz6(mvsiz),bzz7(mvsiz),bzz8(mvsiz)
314C
315 INTEGER NNEGA, INDEX(MVSIZ), ISELECT
316 my_real
317 .
318 .
319 .
320 .
321 .
322 .
323 .
324 .
325 .
326 .
327 .
328 .
329 . AJI1(MVSIZ,NNPT) , AJI2(MVSIZ,NNPT) , AJI3(MVSIZ,NNPT) ,
330 . AJI4(MVSIZ,NNPT) , AJI5(MVSIZ,NNPT) , AJI6(MVSIZ,NNPT) ,
331 . AJI7(MVSIZ,NNPT) , AJI8(MVSIZ,NNPT) , AJI9(MVSIZ,NNPT),
332 . PX1(MVSIZ,NNPT),PX2(MVSIZ,NNPT),PX3(MVSIZ,NNPT),PX4(MVSIZ,NNPT),
333 . px5(mvsiz,nnpt),px6(mvsiz,nnpt),px7(mvsiz,nnpt),px8(mvsiz,nnpt),
334 . py1(mvsiz,nnpt),py2(mvsiz,nnpt),py3(mvsiz,nnpt),py4(mvsiz,nnpt),
335 . py5(mvsiz,nnpt),py6(mvsiz,nnpt),py7(mvsiz,nnpt),py8(mvsiz,nnpt),
336 . pz1(mvsiz,nnpt),pz2(mvsiz,nnpt),pz3(mvsiz,nnpt),pz4(mvsiz,nnpt),
337 . pz5(mvsiz,nnpt),pz6(mvsiz,nnpt),pz7(mvsiz,nnpt),pz8(mvsiz,nnpt),
338 .
339 .
340 .
341 .
342 .
343 .
344 .
345 .
346 .
347 .
348 .
349 .
350 . px1h1(mvsiz), px1h2(mvsiz), px1h3(mvsiz), px1h4(mvsiz),
351 . px2h1(mvsiz), px2h2(mvsiz), px2h3(mvsiz), px2h4(mvsiz),
352 . px3h1(mvsiz), px3h2(mvsiz), px3h3(mvsiz), px3h4(mvsiz),
353 . px4h1(mvsiz), px4h2(mvsiz), px4h3(mvsiz), px4h4(mvsiz),
354 . jr_1(mvsiz),js_1(mvsiz),jt_1(mvsiz),
355 . bb(6,24,mvsiz),amu(mvsiz)
356 double precision
357 . volp(mvsiz,nnpt)
358C----- Variables used for iterative power calculation
359C-----
360C variables used for the non-local
361 my_real,
362 . DIMENSION(:,:), ALLOCATABLE :: var_reg
363 INTEGER :: INLOC, L_NLOC, INOD(8), IPOS(8), IMAT
364 my_real,
365 . DIMENSION(:), POINTER :: DNL
366 my_real
367 . H(8),PS2(8),PR2(8),PT2(8), ZR,ZS,ZT
368C
369 DOUBLE PRECISION
370 . ULX1(MVSIZ), ULX2(MVSIZ), ULX3(MVSIZ), ULX4(MVSIZ),
371 . ULX5(MVSIZ), ULX6(MVSIZ), ULX7(MVSIZ), ULX8(MVSIZ),
372 . ULY1(MVSIZ), ULY2(MVSIZ), ULY3(MVSIZ), ULY4(MVSIZ),
373 . ULY5(MVSIZ), ULY6(MVSIZ), ULY7(MVSIZ), ULY8(MVSIZ),
374 . ulz1(mvsiz), ulz2(mvsiz), ulz3(mvsiz), ulz4(mvsiz),
375 . ulz5(mvsiz), ulz6(mvsiz), ulz7(mvsiz), ulz8(mvsiz),
376 . dn_x(mvsiz,8),dn_y(mvsiz,8),dn_z(mvsiz,8) ,
377 . trm(nel,24,24),dn_r(8),dn_s(8),dn_t(8),invj(9,mvsiz)
378 my_real
379 . a11(mvsiz), a12(mvsiz), a13(mvsiz),
380 . a21(mvsiz), a22(mvsiz), a23(mvsiz),
381 . a31(mvsiz), a32(mvsiz), a33(mvsiz)
382 INTEGER SZ_IX
383C
384 my_real
385 . w_gauss(9,9),a_gauss(9,9)
386 DATA w_gauss /
387 1 2. ,0. ,0. ,
388 1 0. ,0. ,0. ,
389 1 0. ,0. ,0. ,
390 2 1. ,1. ,0. ,
391 2 0. ,0. ,0. ,
392 2 0. ,0. ,0. ,
393 3 0.555555555555556,0.888888888888889,0.555555555555556,
394 3 0. ,0. ,0. ,
395 3 0. ,0. ,0. ,
396 4 0.347854845137454,0.652145154862546,0.652145154862546,
397 4 0.347854845137454,0. ,0. ,
398 4 0. ,0. ,0. ,
399 5 0.236926885056189,0.478628670499366,0.568888888888889,
400 5 0.478628670499366,0.236926885056189,0. ,
401 5 0. ,0. ,0. ,
402 6 0.171324492379170,0.360761573048139,0.467913934572691,
403 6 0.467913934572691,0.360761573048139,0.171324492379170,
404 6 0. ,0. ,0. ,
405 7 0.129484966168870,0.279705391489277,0.381830050505119,
406 7 0.417959183673469,0.381830050505119,0.279705391489277,
407 7 0.129484966168870,0. ,0. ,
408 8 0.101228536290376,0.222381034453374,0.313706645877887,
409 8 0.362683783378362,0.362683783378362,0.313706645877887,
410 8 0.222381034453374,0.101228536290376,0. ,
411 9 0.081274388361574,0.180648160694857,0.260610696402935,
412 9 0.312347077040003,0.330239355001260,0.312347077040003,
413 9 0.260610696402935,0.180648160694857,0.081274388361574/
414 DATA a_gauss /
415 1 0. ,0. ,0. ,
416 1 0. ,0. ,0. ,
417 1 0. ,0. ,0. ,
418 2 -.577350269189626,0.577350269189626,0. ,
419 2 0. ,0. ,0. ,
420 2 0. ,0. ,0. ,
421 3 -.774596669241483,0. ,0.774596669241483,
422 3 0. ,0. ,0. ,
423 3 0. ,0. ,0. ,
424 4 -.861136311594053,-.339981043584856,0.339981043584856,
425 4 0.861136311594053,0. ,0. ,
426 4 0. ,0. ,0. ,
427 5 -.906179845938664,-.538469310105683,0. ,
428 5 0.538469310105683,0.906179845938664,0. ,
429 5 0. ,0. ,0. ,
430 6 -.932469514203152,-.661209386466265,-.238619186083197,
431 6 0.238619186083197,0.661209386466265,0.932469514203152,
432 6 0. ,0. ,0. ,
433 7 -.949107912342759,-.741531185599394,-.405845151377397,
434 7 0. ,0.405845151377397,0.741531185599394,
435 7 0.949107912342759,0. ,0. ,
436 8 -.960289856497536,-.796666477413627,-.525532409916329,
437 8 -.183434642495650,0.183434642495650,0.525532409916329,
438 8 0.796666477413627,0.960289856497536,0. ,
439 9 -.968160239507626,-.836031107326636,-.613371432700590,
440 9 -.324253423403809,0. ,0.324253423403809,
441 9 0.613371432700590,0.836031107326636,0.968160239507626/
442C
443c Flag Bolt Preloading
444 INTEGER IBOLTP,NBPRELD
445 my_real,
446 . DIMENSION(:), POINTER :: BPRELD
447C-----------------------------------------------
448C S o u r c e L i n e s
449C-----------------------------------------------
450 NPTR = elbuf_tab(ng)%NPTR
451 npts = elbuf_tab(ng)%NPTS
452 nptt = elbuf_tab(ng)%NPTT
453 gbuf => elbuf_tab(ng)%GBUF
454C
455 iboltp = iparg(72,ng)
456 inloc = iparg(78,ng)
457 ALLOCATE(var_reg(nel,nptr*npts*nptt))
458 nbpreld = gbuf%G_BPRELD
459 bpreld =>gbuf%BPRELD(1:nbpreld*nel)
460 tempel(:) = zero
461 bid(:) = zero
462 sz_r1_free=mvsiz
463 sz_ix=numelq+numels+nsvois
464c
465 ilay = 1
466 ibid = 0
467 ibidv = 0
468 nf1=nft+1
469C------------------------shear locking is alleviated with ANS for shear
470 iselect=0
471C-----------
472C GATHER NODAL VARIABLES AND COMPUTE INTRINSIC ROTATION.
473C-----------
474C GATHER NODAL VARIABLES AND COMPUTE INTRINSIC ROTATION.
475 CALL srcoor3_imp(x,ixs(1,nf1),v,w,gbuf%GAMA,gama,
476 . x1, x2, x3, x4, x5, x6, x7, x8,
477 . y1, y2, y3, y4, y5, y6, y7, y8,
478 . z1, z2, z3, z4, z5, z6, z7, z8,
479 . vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
480 . vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
481 . vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
482 . vd2,vis,gbuf%OFF,offg,gbuf%SMSTR,gbuf%RHO,rhoo, gbuf%COR_FR,
483 . nc1,nc2,nc3,nc4,nc5,nc6,nc7,nc8,
484 . ngl,mxt,ngeo,ioutprt, vgxa, vgya, vgza, vga2,
485 . xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8,
486 . yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8,
487 . zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8,
488 . xdp, x0 , y0 , z0 , nel, trm, gbuf%COR_XR,
489 . ulx1 ,ulx2 ,ulx3 ,ulx4 ,ulx5 ,ulx6 ,ulx7 ,ulx8 ,
490 . uly1 ,uly2 ,uly3 ,uly4 ,uly5 ,uly6 ,uly7 ,uly8 ,
491 . ulz1 ,ulz2 ,ulz3 ,ulz4 ,ulz5 ,ulz6 ,ulz7 ,ulz8 ,
492 . xgxa ,xgya ,xgza, xgxa2,xgya2,xgza2,xgxya,xgyza,xgzxa,iparg(1,ng))
493c---
494C-----------
495C GATHER NODAL VARIABLES FOR TOTAL STRAIN CASE.
496C-----------
497C-----------
498 ioffs=0
499 DO i=1,nel
500 offs(i)=ep20
501 deltax(i)=ep30
502 ideg(i)=0
503 ENDDO
504 IF(jthe < 0) them(1:nel,1:8) =zero
505
506 CALL s8eprst_ini(pr ,ps ,pt )
507C-------now [B](Pij) is always in global system w/ ISMSTR10
508C-----------------------------
509C----------JACOBIEN-inverse- in case of Vol<=0-------
510 IF (ismstr/=11) THEN
511 CALL s8ederic3(
512 1 offg, volg, ngl, xd1,
513 2 xd2, xd3, xd4, xd5,
514 3 xd6, xd7, xd8, yd1,
515 4 yd2, yd3, yd4, yd5,
516 5 yd6, yd7, yd8, zd1,
517 6 zd2, zd3, zd4, zd5,
518 7 zd6, zd7, zd8, pxc1,
519 8 pxc2, pxc3, pxc4, pyc1,
520 9 pyc2, pyc3, pyc4, pzc1,
521 a pzc2, pzc3, pzc4, px1h1,
522 b px1h2, px1h3, px1h4, px2h1,
523 c px2h2, px2h3, px2h4, px3h1,
524 d px3h2, px3h3, px3h4, px4h1,
525 e px4h2, px4h3, px4h4, hx,
526 f hy, hz, jr_1, js_1,
527 g jt_1, ajc1, ajc2, ajc3,
528 h ajc4, ajc5, ajc6, ajc7,
529 i ajc8, ajc9, smax, gbuf%SMSTR,
530 j gbuf%OFF, nel, ismstr, jlag)
531 CALL s8ejacip3(
532 1 hx, hy, hz, ajc1,
533 2 ajc2, ajc3, ajc4, ajc5,
534 3 ajc6, ajc7, ajc8, ajc9,
535 4 ajp1, ajp2, ajp3, ajp4,
536 5 ajp5, ajp6, ajp7, ajp8,
537 6 ajp9, nel)
538 IF (idts6>0) THEN
539 CALL sdlen_dege(
540 1 volg, deltax, x1, x2,
541 2 x3, x4, x5, x6,
542 3 x7, x8, y1, y2,
543 4 y3, y4, y5, y6,
544 5 y7, y8, z1, z2,
545 6 z3, z4, z5, z6,
546 7 z7, z8, ixs(1,nf1),ideg,
547 8 nel)
548 END IF
549 END IF !(ISMSTR/=11) THEN
550C----------JACOBIEN-inverse- in case of Vol<=0-------
551 nnega = 0
552C-----------Begin integrating points-----
553c
554 DO ir=1,nptr
555 DO is=1,npts
556 DO it=1,nptt
557 ip = ir + ( (is-1) + (it-1)*npts )*nptr
558 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
559C
560 CALL s8ederipr3(
561 1 gbuf%OFF, volp(1,ip),ngl, wi,
562 2 ajp1(1,ip),ajp2(1,ip),ajp3(1,ip),ajp4(1,ip),
563 3 ajp5(1,ip),ajp6(1,ip),ajp7(1,ip),ajp8(1,ip),
564 4 ajp9(1,ip),aji1(1,ip),aji2(1,ip),aji3(1,ip),
565 5 aji4(1,ip),aji5(1,ip),aji6(1,ip),aji7(1,ip),
566 6 aji8(1,ip),aji9(1,ip),nnega, index,
567 7 ip, nel)
568 ENDDO ! IT=1,NPTT
569 ENDDO ! IS=1,NPTS
570 ENDDO ! IR=1,NPTR
571C
572 IF (nnega>0) THEN
573 CALL s8edericm3(
574 1 volg, ngl, xd1, xd2,
575 2 xd3, xd4, xd5, xd6,
576 3 xd7, xd8, yd1, yd2,
577 4 yd3, yd4, yd5, yd6,
578 5 yd7, yd8, zd1, zd2,
579 6 zd3, zd4, zd5, zd6,
580 7 zd7, zd8, pxc1, pxc2,
581 8 pxc3, pxc4, pyc1, pyc2,
582 9 pyc3, pyc4, pzc1, pzc2,
583 a pzc3, pzc4, hx, hy,
584 b hz, ajc1, ajc2, ajc3,
585 c ajc4, ajc5, ajc6, ajc7,
586 d ajc8, ajc9, smax, gbuf%SMSTR,
587 e gbuf%OFF, nnega, index, nel,
588 f ismstr)
589 DO ir=1,nptr
590 DO is=1,npts
591 DO it=1,nptt
592 ip = ir + ( (is-1) + (it-1)*npts )*nptr
593 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
594C
595 CALL s8zderims3(volp(1,ip),
596 . a_gauss(ir,nptr),a_gauss(is,npts),a_gauss(it,nptt),wi,
597 . hx, hy, hz,
598 . ajc1,ajc2,ajc3,
599 . ajc4,ajc5,ajc6,
600 . ajc7,ajc8,ajc9,
601 . ajp1(1,ip),ajp2(1,ip),ajp3(1,ip),
602 . ajp4(1,ip),ajp5(1,ip),ajp6(1,ip),
603 . ajp7(1,ip),ajp8(1,ip),ajp9(1,ip),
604 . aji1(1,ip),aji2(1,ip),aji3(1,ip),
605 . aji4(1,ip),aji5(1,ip),aji6(1,ip),
606 . aji7(1,ip),aji8(1,ip),aji9(1,ip),nnega,index)
607 ENDDO ! IT=1,NPTT
608 ENDDO ! IS=1,NPTS
609 ENDDO ! IR=1,NPTR
610 IF (ineg_v ==0) THEN
611C--- /NEGAVOL/STOP
612 CALL ancmsg(msgid=280,anmode=aninfo)
613 mstop = 1
614 END IF !(INEG_V /=0) THEN
615 END IF
616C --------------------------
617C --- UPDATE REF CONFIGURATION (possible future change to small strain option)
618C --- ! Total strain option doesn't change the Ref CONF.
619C --------------------------
620C
621 IF (icp==1.OR.icp==2) THEN
622 CALL s8edefoc3(
623 1 pxc1, pxc2, pxc3, pxc4,
624 2 pyc1, pyc2, pyc3, pyc4,
625 3 pzc1, pzc2, pzc3, pzc4,
626 4 vx1, vx2, vx3, vx4,
627 5 vx5, vx6, vx7, vx8,
628 6 vy1, vy2, vy3, vy4,
629 7 vy5, vy6, vy7, vy8,
630 8 vz1, vz2, vz3, vz4,
631 9 vz5, vz6, vz7, vz8,
632 a dsv, nel)
633C-----don't do it w/ degenerated elm, but limit to law28 after QA
634 IF (idts6==0) CALL degenes8(
635 1 ixs(1,nf1),ideg, nel)
636 IF (icp==1.AND.mtn==28) THEN
637 DO i=1,nel
638 IF (ideg(i)>0) ideg(i) = ideg(i) + 10
639 ENDDO
640 END IF !(ICP==1.AND.MTN==28) THEN
641 ENDIF
642C-----------
643C INITIALIZATION---Before increment-----
644C-----------
645 CALL s8zzero3(
646 1 f11, f21, f31, f12,
647 2 f22, f32, f13, f23,
648 3 f33, f14, f24, f34,
649 4 f15, f25, f35, f16,
650 5 f26, f36, f17, f27,
651 6 f37, f18, f28, f38,
652 7 gbuf%SIG, gbuf%EINT, gbuf%RHO, gbuf%QVIS,
653 8 gbuf%PLA, gbuf%EPSD, stin, pp,
654 9 gbuf%G_PLA, gbuf%G_EPSD,iexpan, gbuf%EINTTH,
655 a nel, conden)
656C-------------------------------------------
657C COMPUTE AVERAGE TEMPERATURE IN ELEMENT (for output)
658C-------------------------------------------
659 IF(jthe < 0 ) THEN
660 DO i=1,nel
661 IF(gbuf%OFF(i) == zero) cycle
662 tempel(i) = one_over_8 *( temp(nc1(i)) + temp(nc2(i))
663 . + temp(nc3(i)) + temp(nc4(i))
664 . + temp(nc5(i)) + temp(nc6(i))
665 . + temp(nc7(i)) + temp(nc8(i)))
666 gbuf%TEMP(i) = tempel(i)
667 ENDDO
668 ENDIF
669c-------------------------------------------
670c COMPUTE Regularized non local variable in Gauss point
671c-------------------------------------------
672 IF (inloc > 0) THEN
673 l_nloc = nloc_dmg%L_NLOC
674 dnl => nloc_dmg%DNL(1:l_nloc) ! DNL = non local variable increment
675 DO ir=1,nptr
676 DO is=1,npts
677 DO it=1,nptt
678 zr = a_gauss(ir,nptr)
679 zs = a_gauss(is,npts)
680 zt = a_gauss(it,nptt)
681 ip = ir + ( (is-1) + (it-1)*npts )*nptr
682 CALL basis8 (zr,zs,zt,h,pr2,ps2,pt2)
683 DO i=1,nel
684 inod(1) = nloc_dmg%IDXI(nc1(i))
685 inod(2) = nloc_dmg%IDXI(nc2(i))
686 inod(3) = nloc_dmg%IDXI(nc3(i))
687 inod(4) = nloc_dmg%IDXI(nc4(i))
688 inod(5) = nloc_dmg%IDXI(nc5(i))
689 inod(6) = nloc_dmg%IDXI(nc6(i))
690 inod(7) = nloc_dmg%IDXI(nc7(i))
691 inod(8) = nloc_dmg%IDXI(nc8(i))
692 DO j = 1, 8
693 ipos(j) = nloc_dmg%POSI(inod(j))
694 ENDDO
695 var_reg(i,ip) = h(1)*dnl(ipos(1)) + h(2)*dnl(ipos(2)) + h(3)*dnl(ipos(3))
696 . + h(4)*dnl(ipos(4)) + h(5)*dnl(ipos(5)) + h(6)*dnl(ipos(6))
697 . + h(7)*dnl(ipos(7)) + h(8)*dnl(ipos(8))
698 ENDDO
699 ENDDO
700 ENDDO
701 ENDDO
702 ENDIF
703C---------[Bj] first---
704 DO ir=1,nptr
705 DO is=1,npts
706 DO it=1,nptt
707C-----------
708 ip = ir + ( (is-1) + (it-1)*npts )*nptr
709C
710 CALL s8ederig3(
711 1 px1(1,ip), px2(1,ip), px3(1,ip), px4(1,ip),
712 2 px5(1,ip), px6(1,ip), px7(1,ip), px8(1,ip),
713 3 py1(1,ip), py2(1,ip), py3(1,ip), py4(1,ip),
714 4 py5(1,ip), py6(1,ip), py7(1,ip), py8(1,ip),
715 5 pz1(1,ip), pz2(1,ip), pz3(1,ip), pz4(1,ip),
716 6 pz5(1,ip), pz6(1,ip), pz7(1,ip), pz8(1,ip),
717 7 aji1(1,ip),aji2(1,ip),aji3(1,ip),aji4(1,ip),
718 8 aji5(1,ip),aji6(1,ip),aji7(1,ip),aji8(1,ip),
719 9 aji9(1,ip),pr(1,ip), ps(1,ip), pt(1,ip),
720 a nel)
721 ENDDO ! IT=1,NPTT
722 ENDDO ! IS=1,NPTS
723 ENDDO ! IR=1,NPTR
724C-----------Begin integrating points-----
725 DO is=1,npts
726 DO ir=1,nptr
727 DO it=1,nptt
728 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
729C-----avoid multi-print
730 IF (ioffs == 1)THEN
731 DO i=1,nel
732 IF (offs(i)<=two) lbuf%OFF(i)=offs(i)
733 ENDDO
734 END IF
735C-----------
736 ip = ir + ( (is-1) + (it-1)*npts )*nptr
737 wi = w_gauss(ir,nptr)*w_gauss(is,npts)*w_gauss(it,nptt)
738C
739 CALL s8ederi_2(
740 1 offg, off, voln, a_gauss(ir,nptr),
741 2 a_gauss(is,npts),a_gauss(it,nptt),wi, px1(1,ip),
742 3 px2(1,ip), px3(1,ip), px4(1,ip), px5(1,ip),
743 4 px6(1,ip), px7(1,ip), px8(1,ip), py1(1,ip),
744 5 py2(1,ip), py3(1,ip), py4(1,ip), py5(1,ip),
745 6 py6(1,ip), py7(1,ip), py8(1,ip), pz1(1,ip),
746 7 pz2(1,ip), pz3(1,ip), pz4(1,ip), pz5(1,ip),
747 8 pz6(1,ip), pz7(1,ip), pz8(1,ip), pxc1,
748 9 pxc2, pxc3, pxc4, pyc1,
749 a pyc2, pyc3, pyc4, pzc1,
750 b pzc2, pzc3, pzc4, bxy1,
751 c bxy2, bxy3, bxy4, bxy5,
752 d bxy6, bxy7, bxy8, byx1,
753 e byx2, byx3, byx4, byx5,
754 f byx6, byx7, byx8, bxz1,
755 g bxz2, bxz3, bxz4, bxz5,
756 h bxz6, bxz7, bxz8, bzx1,
757 i bzx2, bzx3, bzx4, bzx5,
758 j bzx6, bzx7, bzx8, byz1,
759 k byz2, byz3, byz4, byz5,
760 l byz6, byz7, byz8, bzy1,
761 m bzy2, bzy3, bzy4, bzy5,
762 n bzy6, bzy7, bzy8, bxx1,
763 o bxx2, bxx3, bxx4, bxx5,
764 p bxx6, bxx7, bxx8, byy1,
765 q byy2, byy3, byy4, byy5,
766 r byy6, byy7, byy8, bzz1,
767 s bzz2, bzz3, bzz4, bzz5,
768 t bzz6, bzz7, bzz8, ajp4(1,ip),
769 u ajp5(1,ip), ajp6(1,ip), ajp7(1,ip), ajp8(1,ip),
770 v ajp9(1,ip), aj1, aj2, aj3,
771 w aj4, aj5, aj6, smax,
772 x deltax, icp, ideg, nu,
773 y volp(1,ip), nel)
774C
775 CALL s8sderi3(
776 1 offg, off, volp(1,ip), ngl,
777 2 a_gauss(it,nptt),a_gauss(ir,nptr),a_gauss(is,npts),wi,
778 3 xd1, xd2, xd3, xd4,
779 4 xd5, xd6, xd7, xd8,
780 5 yd1, yd2, yd3, yd4,
781 6 yd5, yd6, yd7, yd8,
782 7 zd1, zd2, zd3, zd4,
783 8 zd5, zd6, zd7, zd8,
784 9 a11, a12, a13, a21,
785 a a22, a23, a31, a32,
786 b a33, dn_r, dn_s, dn_t,
787 c invj, dn_x, dn_y, dn_z,
788 d voln, nel)
789 CALL s8sbdefo3(
790 1 ulx1, ulx2, ulx3, ulx4,
791 2 ulx5, ulx6, ulx7, ulx8,
792 3 uly1, uly2, uly3, uly4,
793 4 uly5, uly6, uly7, uly8,
794 5 ulz1, ulz2, ulz3, ulz4,
795 6 ulz5, ulz6, ulz7, ulz8,
796 7 xd1, xd2, xd3, xd4,
797 8 xd5, xd6, xd7, xd8,
798 9 yd1, yd2, yd3, yd4,
799 a yd5, yd6, yd7, yd8,
800 b zd1, zd2, zd3, zd4,
801 c zd5, zd6, zd7, zd8,
802 d invj, a_gauss(it,nptt),a_gauss(ir,nptr),a_gauss(is,npts),
803 e a11, a12, a13, a21,
804 f a22, a23, a31, a32,
805 g a33, dn_r, dn_s, dn_t,
806 h bb, dxx, dxy, dxz,
807 i dyx, dyy, dyz, dzx,
808 j dzy, dzz, d4, d5,
809 k d6, wxx, wyy, wzz,
810 l lbuf%VOL, off, lbuf%EINT, gbuf%OFF,
811 m dsv, icp, fac, sdv,
812 n iselect, ideg, nel, ismstr)
813c
814 DO i=1,nel
815 rhoo(i) = lbuf%RHO(i)
816 ENDDO
817 divde(1:nel) = dt1*(dxx(1:nel)+ dyy(1:nel)+ dzz(1:nel))+sdv(1:nel)
818 CALL srho3(
819 1 pm, lbuf%VOL, lbuf%RHO, lbuf%EINT,
820 2 divde, flux(1,nf1),flu1(nf1), voln,
821 3 dvol, ngl, mxt, off,
822 4 0, gbuf%TAG22, volp(1,ip), lbuf%VOL0DP,
823 5 amu, gbuf%OFF, nel, mtn,
824 6 jale, ismstr, jeul, jlag)
825C
826C-----------------------------
827C EXTRACT STRESSES + SMALL STRAIN
828C-----------------------------
829 CALL csmall3(lbuf%SIG,s1,s2,s3,s4,s5,s6,
830 . gbuf%OFF,off,nel)
831C
832C for heat transfert
833C
834C------------------------------------------------------
835C CALCUL DES CONTRAINTES SUIVANT LOIS CONSTITUTIVES
836C TIME STEP EST CALCULE ICI DT2T
837C------------------------------------------------------
838 IF ((itask==0).AND.(imon_mat==1)) CALL startime(timers,35)
839c
840c
841 CALL mmain(timers,output,
842 1 elbuf_tab, ng, pm, geo,
843 2 ale_connect, ixs, iparg,
844 3 v, tf, npf, bufmat,
845 4 sti, x, dt2t, neltst,
846 5 ityptst, offset, nel, w,
847 6 off, ngeo, mxt, ngl,
848 7 voln, vd2, dvol, deltax,
849 8 vis, qvis, cxx, s1,
850 9 s2, s3, s4, s5,
851 a s6, dxx, dyy, dzz,
852 b d4, d5, d6, wxx,
853 c wyy, wzz, aj1, aj2,
854 d aj3, aj4, aj5, aj6,
855 e vdx, vdy, vdz, muvoid,
856 f ssp_eq, aire, sigy, et,
857 g r1_free, lbuf%PLA, r3_free, amu,
858 h mfxx(1,ip), mfxy(1,ip), mfxz(1,ip), mfyx(1,ip),
859 i mfyy(1,ip), mfyz(1,ip), mfzx(1,ip), mfzy(1,ip),
860 j mfzz(1,ip), ipm, gama, bid,
861 k bid, bid, bid, bid,
862 l bid, bid, istrain, tempel,
863 m die, iexpan, ilay, mssa,
864 n dmels, ir, is, it,
865 o table, bid, bid, bid,
866 p bid, iparg(1,ng), igeo, conde,
867 q itask, nloc_dmg, var_reg(1,ip),mat_elem,
868 r h3d_strain, jplasol, jsph, sz_r1_free,
869 s snpc, stf, sbufmat, glob_therm,
870 * svis, sz_ix, iresp,
871 t n2d, th_strain, ngroup, tt,
872 . dt1, ntable, numelq, nummat,
873 . numgeo, numnod, numels,
874 . idel7nok, idtmin, maxfunc,
875 . imon_mat, userl_avail, impl_s,
876 . idyna, dt, bid ,sensors)
877C
878 IF (istrain == 1)THEN
879 CALL sstra3(
880 1 dxx, dyy, dzz, d4,
881 2 d5, d6, lbuf%STRA,wxx,
882 3 wyy, wzz, off, nel,
883 4 jcvt)
884 ENDIF
885 IF ((itask==0).AND.(imon_mat==1)) CALL stoptime(timers,35)
886C----------------------------
887C INTERNAL FORCES
888C----------------------------
889 CALL s8sfint3(
890 1 lbuf%SIG, f11, f21, f31,
891 2 f12, f22, f32, f13,
892 3 f23, f33, f14, f24,
893 4 f34, f15, f25, f35,
894 5 f16, f26, f36, f17,
895 6 f27, f37, f18, f28,
896 7 f38, dn_x, dn_y, dn_z,
897 8 bb, voln, qvis, icp,
898 9 jfac(1,ip),nel, iselect, ideg,
899 a ismstr, svis)
900 CALL s8efmoy3(
901 1 lbuf%SIG, voln, qvis, pp,
902 2 lbuf%EINT, lbuf%RHO, lbuf%QVIS, lbuf%PLA,
903 3 lbuf%EPSD, gbuf%EPSD, gbuf%SIG, gbuf%EINT,
904 4 gbuf%RHO, gbuf%QVIS, gbuf%PLA, volg,
905 5 sti, stin, icp, off,
906 6 lbuf%VOL, gbuf%VOL, gbuf%G_PLA, gbuf%G_EPSD,
907 7 lbuf%EINTTH,gbuf%EINTTH,iexpan, nel,
908 8 conde, conden,svis,glob_therm%NODADT_THERM,
909 9 gbuf%WPLA, lbuf%WPLA, gbuf%G_WPLA)
910c-------------------------
911c Virtual internal forces of regularized non local ddl
912c--------------------------
913 IF (inloc > 0) THEN
914 imat = mxt(1)
915 zr = a_gauss(ir,nptr)
916 zs = a_gauss(is,npts)
917 zt = a_gauss(it,nptt)
918 CALL basis8 (zr,zs,zt,h,pr2,ps2,pt2)
919 CALL s8zfint_reg(
920 1 nloc_dmg, var_reg(1,ip),nel, gbuf%OFF,
921 2 voln, nc1, nc2, nc3,
922 3 nc4, nc5, nc6, nc7,
923 4 nc8, px1, px2, px3,
924 5 px4, px5, px6, px7,
925 6 px8, py1, py2, py3,
926 7 py4, py5, py6, py7,
927 8 py8, pz1, pz2, pz3,
928 9 pz4, pz5, pz6, pz7,
929 a pz8, imat, h, wi,
930 b ip, itask, dt2t, gbuf%VOL,
931 c nft)
932 ENDIF
933 DO i=1,nel
934 offg(i)=min(offg(i),off(i))
935 IF (lbuf%OFF(i)>one .AND. gbuf%OFF(i) == one) THEN
936 offs(i)=min(lbuf%OFF(i),offs(i))
937 ioffs =1
938 END IF
939 ENDDO
940C-------------------------
941c finite element heat transfert
942C--------------------------
943C
944 ENDDO ! IT=1,NPTT
945 ENDDO ! IS=1,NPTS
946 ENDDO ! IR=1,NPTR
947C-----------End integrating points-----
948c
949 IF (ioffs == 1) THEN
950 DO i=1,nel
951 IF (offs(i)<=two) gbuf%OFF(i)=offs(i)
952 END DO
953c-------------------
954 DO ir=1,nptr
955 DO is=1,npts
956 DO it=1,nptt
957 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
958 ip = ir + ( (is-1) + (it-1)*npts )*nptr
959 DO i=1,nel
960 IF (gbuf%OFF(i) > one) lbuf%OFF(i)=gbuf%OFF(i)
961 END DO
962 END DO
963 END DO
964 END DO
965 END IF ! IOFFS == 1
966 IF(icp==1.AND.ismstr/=10.AND.ismstr/=12) THEN
967 CALL s8zfintp3(
968 1 pxc1, pxc2, pxc3, pxc4,
969 2 pyc1, pyc2, pyc3, pyc4,
970 3 pzc1, pzc2, pzc3, pzc4,
971 4 f11, f21, f31, f12,
972 5 f22, f32, f13, f23,
973 6 f33, f14, f24, f34,
974 7 f15, f25, f35, f16,
975 8 f26, f36, f17, f27,
976 9 f37, f18, f28, f38,
977 a volg, pp, ideg, nel)
978 ENDIF
979C-----------------------------
980C SMALL STRAIN
981C-----------------------------
982 CALL smallb3(
983 1 gbuf%OFF,offg, nel, ismstr)
984 IF (ismstr==11.OR.(jcvt == 0 .AND. ismstr>0)) THEN
985 CALL s8edefw3(
986 1 pxc1, pxc2, pxc3, pxc4,
987 2 pyc1, pyc2, pyc3, pyc4,
988 3 pzc1, pzc2, pzc3, pzc4,
989 4 vx1, vx2, vx3, vx4,
990 5 vx5, vx6, vx7, vx8,
991 6 vy1, vy2, vy3, vy4,
992 7 vy5, vy6, vy7, vy8,
993 8 vz1, vz2, vz3, vz4,
994 9 vz5, vz6, vz7, vz8,
995 a wxx0, wyy0, wzz0, nel)
996 CALL smallg3(
997 1 gbuf%SMSTR,offg, wxx0, wyy0,
998 2 wzz0, r11, r12, r13,
999 3 r21, r22, r23, r31,
1000 4 r32, r33, nel, ismstr,
1001 5 jcvt)
1002 END IF
1003 IF (inconv==1.AND.tt>zero ) THEN
1004 CALL s8xref_imp(
1005 1 gbuf%OFF, gbuf%COR_XR,x1, x2,
1006 2 x3, x4, x5, x6,
1007 3 x7, x8, y1, y2,
1008 4 y3, y4, y5, y6,
1009 5 y7, y8, z1, z2,
1010 6 z3, z4, z5, z6,
1011 7 z7, z8, gbuf%COR_FR,nel)
1012 END IF
1013
1014C--------------------------
1015C balances by material
1016C--------------------------
1017 iflag=mod(ncycle,ncpri)
1018 IF(ioutprt>0)THEN
1019 IF(jcvt == 0) THEN
1020 ifvm22 = 0
1021 CALL sbilan(partsav,gbuf%EINT,gbuf%RHO,gbuf%RK,gbuf%VOL,
1022 . vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
1023 . vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
1024 . vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
1025 . volg,iparts,gresav,grth,igrth,gbuf%OFF,
1026 . iexpan,gbuf%EINTTH,gbuf%FILL,gbuf%MOM,nel,ifvm22,
1027 . x1, x2, x3, x4, x5, x6, x7, x8,
1028 . y1, y2, y3, y4, y5, y6, y7, y8,
1029 . z1, z2, z3, z4, z5, z6, z7, z8,
1030 . itask,iparg(1,ng),sensors,gbuf%G_WPLA,gbuf%WPLA)
1031 ELSE
1032 CALL srbilan(partsav,gbuf%EINT,gbuf%RHO,gbuf%RK,gbuf%VOL,
1033 . vgxa, vgya, vgza, vga2, volg,iparts,
1034 . gresav,grth,igrth,gbuf%OFF,iexpan,gbuf%EINTTH,
1035 . gbuf%FILL, xgxa, xgya, xgza,
1036 . xgxa2,xgya2,xgza2,xgxya,xgyza,xgzxa,itask,iparg(1,ng),sensors,
1037 . nel,gbuf%G_WPLA,gbuf%WPLA)
1038 ENDIF
1039 ENDIF
1040C----------------------------
1041C CONVECTE --> GLOBAL.
1042C----------------------------
1043 CALL s8sfint3_crimp(trm,gbuf%COR_NF,gbuf%COR_FR,
1044 . f11,f21,f31,f12,f22,f32,f13,f23,f33,f14,f24,f34,
1045 . f15,f25,f35,f16,f26,f36,f17,f27,f37,f18,f28,f38,
1046 . nel )
1047C----------------------------
1048 IF(nfilsol/=0) CALL sfillopt(
1049 1 gbuf%FILL,stin, f11, f21,
1050 2 f31, f12, f22, f32,
1051 3 f13, f23, f33, f14,
1052 4 f24, f34, f15, f25,
1053 5 f35, f16, f26, f36,
1054 6 f17, f27, f37, f18,
1055 7 f28, f38, nel)
1056C----------------------------
1057 IF(iparit == 0)THEN
1058 CALL scumu3(
1059 1 gbuf%OFF,a, nc1, nc2,
1060 2 nc3, nc4, nc5, nc6,
1061 3 nc7, nc8, stifn, stin,
1062 4 f11, f21, f31, f12,
1063 5 f22, f32, f13, f23,
1064 6 f33, f14, f24, f34,
1065 7 f15, f25, f35, f16,
1066 8 f26, f36, f17, f27,
1067 9 f37, f18, f28, f38,
1068 a nvc, bid, bid, bid,
1069 b bid, bid, bid, bid,
1070 c bid, bid, bid, bid,
1071 d bid, bid, bid, bid,
1072 e bid, bid, bid, bid,
1073 f bid, bid, bid, bid,
1074 g bid, bid, bid, bid,
1075 h them, fthe, condn, conden,
1076 i nel, jthe, isrot, ipartsph,glob_therm%NODADT_THERM)
1077 ELSE
1078 CALL scumu3p(
1079 1 gbuf%OFF,stin, fsky, fsky,
1080 2 iads, f11, f21, f31,
1081 3 f12, f22, f32, f13,
1082 4 f23, f33, f14, f24,
1083 5 f34, f15, f25, f35,
1084 6 f16, f26, f36, f17,
1085 7 f27, f37, f18, f28,
1086 8 f38, nc1, nc2, nc3,
1087 9 nc4, nc5, nc6, nc7,
1088 a nc8, bid, bid, bid,
1089 b bid, bid, bid, bid,
1090 c bid, bid, bid, bid,
1091 d bid, bid, bid, bid,
1092 e bid, bid, bid, bid,
1093 f bid, bid, bid, bid,
1094 g bid, bid, bid, bid,
1095 h them, fthesky, condnsky,conden,
1096 i nel, nft, jthe, isrot,
1097 j ipartsph,glob_therm%NODADT_THERM)
1098 ENDIF
1099C-----------
1100 IF (ALLOCATED(var_reg)) DEALLOCATE(var_reg)
1101 RETURN
1102 END
1103
my_real
#define my_real
Definition cppsort.cpp:32
csmall3
subroutine csmall3(sig, s1, s2, s3, s4, s5, s6, offg, off, nel)
Definition csmall3.F:36
degenes8
subroutine degenes8(ixs, idege, nel)
Definition degenes8.F:39
min
#define min(a, b)
Definition macros.h:20
mmain
subroutine mmain(pm, elbuf_str, ix, nix, x, geo, iparg, nel, skew, bufmat, ipart, ipartel, nummat, matparam, imat, ipm, ngl, pid, npf, tf, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, rx, ry, rz, sx, sy, sz, gama, voln, dvol, s1, s2, s3, s4, s5, s6, dxx, dyy, dzz, d4, d5, d6, wxx, wyy, wzz)
Definition mmain.F:43
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
dt_mod
Definition dt_mod.F:80
message_mod
Definition message_mod.F:1257
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
table_mod
Definition table_mod.F:113
s8edefoc3
subroutine s8edefoc3(pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dvca, nel)
Definition s8edefoc3.F:40
s8edefw3
subroutine s8edefw3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, wxx, wyy, wzz, nel)
Definition s8edefw3.F:41
s8ederi_2
subroutine s8ederi_2(offg, off, vol, ksi, eta, zeta, wi, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, pxc1, pxc2, pxc3, pxc4, pyc1, pyc2, pyc3, pyc4, pzc1, pzc2, pzc3, pzc4, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, bxx1, bxx2, bxx3, bxx4, bxx5, bxx6, bxx7, bxx8, byy1, byy2, byy3, byy4, byy5, byy6, byy7, byy8, bzz1, bzz2, bzz3, bzz4, bzz5, bzz6, bzz7, bzz8, aj4, aj5, aj6, aj7, aj8, aj9, rx, ry, rz, sx, sy, sz, smax, deltax, icp, ideg, nu, volp, nel)
Definition s8ederi_2.F:66
s8ederic3
subroutine s8ederic3(off, det, ngl, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px1h1, px1h2, px1h3, px1h4, px2h1, px2h2, px2h3, px2h4, px3h1, px3h2, px3h3, px3h4, px4h1, px4h2, px4h3, px4h4, hx, hy, hz, jr_1, js_1, jt_1, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, smax, sav, offg, nel, ismstr, jlag)
Definition s8ederic3.F:55
s8edericm3
subroutine s8edericm3(det, ngl, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, hx, hy, hz, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, smax, sav, offg, nnega, index, nel, ismstr)
Definition s8edericm3.F:49
s8ederig3
subroutine s8ederig3(px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, aji1, aji2, aji3, aji4, aji5, aji6, aji7, aji8, aji9, pr, ps, pt, nel)
Definition s8ederig3.F:40
s8ederipr3
subroutine s8ederipr3(offg, voldp, ngl, wi, aj1, aj2, aj3, aj4, aj5, aj6, aj7, aj8, aj9, aji1, aji2, aji3, aji4, aji5, aji6, aji7, aji8, aji9, nnega, index, ipt, nel)
Definition s8ederipr3.F:39
s8efmoy3
subroutine s8efmoy3(sigor, vol, qvis, pp, eint, rho, q, defp, epsd, epsdm, sigm, eintm, rhom, qm, defpm, volg, sti, stin, icp, off, vol0, vol0g, g_pla, g_epsd, eintth, eintthm, iexpan, nel, conde, conden, svis, nodadt_therm, g_wpla, l_wpla, g_wpla_flag)
Definition s8efmoy3.F:40
s8eprst_ini
subroutine s8eprst_ini(pr, ps, pt)
Definition s8eprst_ini.F:30
s8sbdefo3
subroutine s8sbdefo3(ulx1, ulx2, ulx3, ulx4, ulx5, ulx6, ulx7, ulx8, uly1, uly2, uly3, uly4, uly5, uly6, uly7, uly8, ulz1, ulz2, ulz3, ulz4, ulz5, ulz6, ulz7, ulz8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, invj, ksi, eta, zeta, a11, a12, a13, a21, a22, a23, a31, a32, a33, dn_r, dn_s, dn_t, bb, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz, volo, off, eint, offs, dsv, icp, fac, sdv, i_sh, idege, nel, ismstr)
Definition s8sbdefo3.F:52
s8sderi3
subroutine s8sderi3(offg, off, voldp, ngl, ksi, eta, zeta, wi, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, a11, a12, a13, a21, a22, a23, a31, a32, a33, dn_r, dn_s, dn_t, invj, dn_x, dn_y, dn_z, voln, nel)
Definition s8sderi3.F:43
s8sfint3
subroutine s8sfint3(sig, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, dn_x, dn_y, dn_z, bb, vol, qvis, icp, jfac, nel, i_sh, idege, ismstr, svis)
Definition s8sfint3.F:39
s8sfint3_crimp
subroutine s8sfint3_crimp(trm, qf, r, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel)
Definition s8sfint3_crimp.F:32
s8sforc3
subroutine s8sforc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, nloc_dmg, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, icp, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, nvc, ipm, itask, istrain, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, igeo, xdp, voln, condn, condnsky, d, ioutprt, mat_elem, h3d_strain, snpc, stf, sbufmat, svis, nsvois, idtmins, iresp, idtmin, maxfunc, userl_avail, dt, glob_therm, sensors)
Definition s8sforc3.F:99
s8xref_imp
subroutine s8xref_imp(offg, xref, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, r, nel)
Definition s8xref_imp.F:38
s8zderims3
subroutine s8zderims3(voldp, ksi, eta, zeta, wi, hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9, jaci1, jaci2, jaci3, jaci4, jaci5, jaci6, jaci7, jaci8, jaci9, nnega, index)
Definition s8zderims3.F:42
s8zfint_reg
subroutine s8zfint_reg(nloc_dmg, var_reg, nel, offg, vol, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, imat, h, wi, ip, itask, dt2t, vol0, nft)
Definition s8zfint_reg.F:45
s8zfintp3
subroutine s8zfintp3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, vol, pp, idege, nel)
Definition s8zfintp3.F:42
s8zzero3
subroutine s8zzero3(fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, fx5, fy5, fz5, fx6, fy6, fz6, fx7, fy7, fz7, fx8, fy8, fz8, sigm, eintm, rhom, qm, defpm, epsdm, stin, pp, g_pla, g_epsd, iexpan, eintthm, nel, conden)
Definition s8zzero3.F:42
sbilan
subroutine sbilan(partsav, eint, rho, rk, vol, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, vnew, iparts, gresav, grth, igrth, off, iexpan, eintth, fill, mom, nel, ifvm22, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, itask, iparg, sensors, g_wpla, wpla)
Definition sbilan.F:47
scumu3
subroutine scumu3(offg, e, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nvc, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthe, condn, conde, nel, jthe, isrot, ipartsph, nodadt_therm)
Definition scumu3.F:55
scumu3p
subroutine scumu3p(offg, sti, fsky, fskyv, iads, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthesky, condnsky, conde, nel, nft, jthe, isrot, ipartsph, nodadt_therm)
Definition scumu3p.F:56
sdlen_dege
subroutine sdlen_dege(volg, lat, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, ixs, idege, nel)
Definition sdlen_dege.F:47
sfillopt
subroutine sfillopt(fill, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel)
Definition sfillopt.F:43
smallb3
subroutine smallb3(offg, off, nel, ismstr)
Definition smallb3.F:45
smallg3
subroutine smallg3(sav, offg, wxx, wyy, wzz, r11, r12, r13, r21, r22, r23, r31, r32, r33, nel, ismstr, jcvt)
Definition smallg3.F:36
srbilan
subroutine srbilan(partsav, eint, rho, rk, vol, vxa, vya, vza, va2, vnew, iparts, gresav, grth, igrth, off, iexpan, eintth, fill, xx, yy, zz, xx2, yy2, zz2, xy, yz, zx, itask, iparg, sensors, nel, g_wpla, wpla)
Definition srbilan.F:45
srcoor3_imp
subroutine srcoor3_imp(x, ixs, v, w, gama0, gama, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, vd2, vis, offg, off, sav, rho, rhoo, r, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ngl, mxt, ngeo, ioutprt, vgax, vgay, vgaz, vga2, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, xdp, x0, y0, z0, nel, trm, xref, ulx1, ulx2, ulx3, ulx4, ulx5, ulx6, ulx7, ulx8, uly1, uly2, uly3, uly4, uly5, uly6, uly7, uly8, ulz1, ulz2, ulz3, ulz4, ulz5, ulz6, ulz7, ulz8, xgax, xgay, xgaz, xgxa2, xgya2, xgza2, xgxya, xgyza, xgzxa, iparg)
Definition srcoor3_imp.F:55
sstra3
subroutine sstra3(dxx, dyy, dzz, d4, d5, d6, strain, wxx, wyy, wzz, off, nel, jcvt)
Definition sstra3.F:47
basis8
subroutine basis8(r, s, t, h, pr, ps, pt)
Definition basis8.F:30
s8ejacip3
subroutine s8ejacip3(hx, hy, hz, cj1, cj2, cj3, cj4, cj5, cj6, cj7, cj8, cj9, jac1, jac2, jac3, jac4, jac5, jac6, jac7, jac8, jac9)
Definition s8zderi3.F:2180
srho3
subroutine srho3(pm, volo, rhon, eint, dxx, dyy, dzz, voln, dvol, mat)
Definition srho3.F:31
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:256
table_mod::ttable
Definition table_mod.F:126