OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
tensors.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!|| tensors ../starter/source/output/anim/tensors.F
25!||--- called by ------------------------------------------------------
26!|| genani1 ../starter/source/output/anim/genani1.f
27!||--- calls -----------------------------------------------------
28!|| initbuf ../starter/source/output/anim/initbuf.F
29!|| srota6 ../starter/source/output/anim/srota6.F
30!||--- uses -----------------------------------------------------
31!||====================================================================
32 SUBROUTINE tensors(ELBUF_TAB,IPARG ,ITENS ,IXS ,PM ,
33 2 EL2FA ,NBF ,TENS ,
34 3 X ,IPART ,IPARTSP ,ISPH3D ,IPM )
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE elbufdef_mod
39 use element_mod , only : nixs
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "vect01_c.inc"
48#include "mvsiz_p.inc"
49#include "com01_c.inc"
50#include "com04_c.inc"
51#include "param_c.inc"
52#include "scr17_c.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56C REAL
57 my_real
58 . tens(6,*),pm(npropm,*),x(3,*)
59 INTEGER IPARG(NPARG,*),ITENS, ISPH3D,NBF,
60 . IXS(NIXS,*),EL2FA(*),IPM(NPROPMI,*),IPART(LIPART1,*),IPARTSP(*)
61 REAL WA(6*NBF)
62 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
63C-----------------------------------------------
64C REAL
65 my_real
66 . evar(6,mvsiz),off, fac, a1, a2, a3, thk, gama(6)
67 REAL R4(18)
68 INTEGER I, II,NG, NEL, IPT,MT1,NLAY,IL,
69 . n, j, mlw,ialel,nn1,nn2,nn3,iprt, kcvt,isolnod,
70 . nptr, npts, nptt,nptg, is, ir, it,
71 . jhbe, jivf, jclose, jplasol, irep, igtyp,
72 . icsen, isorthg, ifailure, iint,jj(6)
73 TYPE(g_bufel_) ,POINTER :: GBUF
74 TYPE(L_BUFEL_) ,POINTER :: LBUF
75C=======================================================================
76 DO j=1,18
77 r4(j) = zero
78 ENDDO
79 nn1 = 1
80 nn2 = 1
81 nn3 = nn2 + numels
82C
83 DO 490 ng=1,ngroup
84 isolnod = iparg(28,ng)
85 CALL initbuf (iparg ,ng ,
86 2 mlw ,nel ,nft ,iad ,ity ,
87 3 npt ,jale ,ismstr ,jeul ,jtur ,
88 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
89 5 jpor ,kcvt ,jclose ,jplasol ,
90 6 irep ,iint ,igtyp ,israt ,isrot ,
91 7 icsen ,isorth ,isorthg ,ifailure)
92 lft=1
93 llt=nel
94!
95 DO j=1,6 ! length max of GBUF%L_STRA = 6
96 jj(j) = (j-1)*nel
97 ENDDO
98!
99C-----------------------------------------------
100C SOLID 8N
101C-----------------------------------------------
102 IF (ity == 1) THEN
103 gbuf => elbuf_tab(ng)%GBUF
104 nlay = elbuf_tab(ng)%NLAY
105 nptr = elbuf_tab(ng)%NPTR
106 npts = elbuf_tab(ng)%NPTS
107 nptt = elbuf_tab(ng)%NPTT
108 nptg = nptt*npts*nptr
109 npt = nptg*nlay
110 mt1 = ixs(1,1 + nft)
111 IF (kcvt==1 .AND. isorth/=0) kcvt=2
112C-----------------------------------------------
113 IF (itens == 1)THEN
114C STRESS
115C-----------------------------------------------
116 DO i=lft,llt
117 n = i + nft
118 evar(1,i) = gbuf%SIG(jj(1) + i)
119 evar(2,i) = gbuf%SIG(jj(2) + i)
120 evar(3,i) = gbuf%SIG(jj(3) + i)
121 evar(4,i) = gbuf%SIG(jj(4) + i)
122 evar(5,i) = gbuf%SIG(jj(5) + i)
123 evar(6,i) = gbuf%SIG(jj(6) + i)
124 ENDDO
125 IF (kcvt/=0) THEN
126C STRESS TENSOR IN GLOBAL SYSTEM
127 DO i=lft,llt
128 n = i + nft
129 IF(el2fa(nn2+n)/=0)THEN
130C For jhbe = 14, average values is in rep.Corota.
131 IF (kcvt==2.AND.jhbe/=14) THEN
132 gama(1)=gbuf%GAMA(jj(1) + i)
133 gama(2)=gbuf%GAMA(jj(2) + i)
134 gama(3)=gbuf%GAMA(jj(3) + i)
135 gama(4)=gbuf%GAMA(jj(4) + i)
136 gama(5)=gbuf%GAMA(jj(5) + i)
137 gama(6)=gbuf%GAMA(jj(6) + i)
138 ELSE
139 gama(1)=one
140 gama(2)=zero
141 gama(3)=zero
142 gama(4)=zero
143 gama(5)=one
144 gama(6)=zero
145 END IF
146 CALL srota6(x,ixs(1:nixs,n),kcvt,evar(1:6,i),gama)
147 ENDIF
148 ENDDO
149 ENDIF
150C-----------------------------------------------
151 ELSEIF(itens == 2)THEN
152C STRAIN
153C-----------------------------------------------
154 DO i=lft,llt
155 DO j=1,6
156 evar(j,i) = zero
157 ENDDO
158 ENDDO
159c
160 IF (isolnod == 8.AND.npt == 8.AND.jhbe/=14.AND.
161 . jhbe/=24.AND.jhbe/=15 )THEN
162 IF (mlw >= 28) THEN
163 DO i=lft,llt
164 n = i + nft
165 DO j=1,8
166 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,j)
167 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1)+i)*one_over_8
168 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2)+i)*one_over_8
169 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3)+i)*one_over_8
170 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4)+i)*one_over_8
171 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5)+i)*one_over_8
172 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6)+i)*one_over_8
173 ENDDO
174 ENDDO
175 ENDIF
176c
177 ELSEIF(isolnod == 8 .AND. npt == 1 .AND.
178 . jhbe/=14.AND.jhbe/=15.AND.jhbe/=24)THEN
179 IF (mlw>=28) THEN
180 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
181 DO i=lft,llt
182 n = i + nft
183 evar(1,i) = lbuf%STRA(jj(1)+i)*one_over_8
184 evar(2,i) = lbuf%STRA(jj(2)+i)*one_over_8
185 evar(3,i) = lbuf%STRA(jj(3)+i)*one_over_8
186 evar(4,i) = lbuf%STRA(jj(4)+i)*one_over_8
187 evar(5,i) = lbuf%STRA(jj(5)+i)*one_over_8
188 evar(6,i) = lbuf%STRA(jj(6)+i)*one_over_8
189 ENDDO
190 ELSEIF(mlw == 14)THEN
191 DO i=lft,llt
192 n = i + nft
193 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1)+i)
194 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2)+i)
195 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3)+i)
196 ENDDO
197 ELSEIF(mlw == 24)THEN
198 DO i=lft,llt
199 n = i + nft
200 evar(1,i) = lbuf%STRA(jj(1)+i)
201 evar(2,i) = lbuf%STRA(jj(2)+i)
202 evar(3,i) = lbuf%STRA(jj(3)+i)
203 evar(4,i) = lbuf%STRA(jj(4)+i)*half
204 evar(5,i) = lbuf%STRA(jj(5)+i)*half
205 evar(6,i) = lbuf%STRA(jj(6)+i)*half
206 ENDDO
207 ENDIF
208c
209 ELSEIF (isolnod == 16 .OR. isolnod == 20 .OR.
210 . (isolnod == 8 .AND. jhbe == 14)) THEN
211 IF (mlw>=28) THEN
212 DO i=lft,llt
213 DO il=1,nlay
214 DO is=1,npts
215 DO it=1,nptt
216 DO ir=1,nptr
217 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
218 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1)+i)/nptg
219 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2)+i)/nptg
220 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3)+i)/nptg
221 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4)+i)/nptg
222 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5)+i)/nptg
223 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6)+i)/nptg
224 ENDDO
225 ENDDO
226 ENDDO
227 ENDDO
228 ENDDO
229 ELSEIF (mlw == 14) THEN
230 DO i=lft,llt
231 n = i + nft
232 DO il=1,nlay
233 DO is=1,npts
234 DO it=1,nptt
235 DO ir=1,nptr
236 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
237 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1)+i)/nptg
238 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2)+i)/nptg
239 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3)+i)/nptg
240 ENDDO
241 ENDDO
242 ENDDO
243 ENDDO
244 ENDDO
245 ELSEIF(mlw == 24)THEN
246 DO i=lft,llt
247 n = i + nft
248 DO il=1,nlay
249 DO is=1,npts
250 DO it=1,nptt
251 DO ir=1,nptr
252 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
253 IF (elbuf_tab(ng)%BUFLY(il)%L_STRA > 0) THEN
254 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1)+i)/nptg
255 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2)+i)/nptg
256 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3)+i)/nptg
257 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4)+i)/nptg
258 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5)+i)/nptg
259 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6)+i)/nptg
260 ENDIF
261 ENDDO
262 ENDDO
263 ENDDO
264 ENDDO
265 ENDDO
266 ENDIF
267c
268 ELSEIF (isolnod == 10) THEN
269 IF(mlw>=28)THEN
270 DO i=lft,llt
271 n = i + nft
272 DO ipt=1,npt
273 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
274 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1)+i)/npt
275 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2)+i)/npt
276 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3)+i)/npt
277 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4)+i)/npt
278 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5)+i)/npt
279 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6)+i)/npt
280 ENDDO
281 ENDDO
282 ELSEIF(mlw == 14)THEN
283 DO i=lft,llt
284 n = i + nft
285 DO ipt=1,npt
286 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
287 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1)+i)/npt
288 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2)+i)/npt
289 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3)+i)/npt
290 ENDDO
291 ENDDO
292 ELSEIF (mlw == 24) THEN
293 DO i=lft,llt
294 n = i + nft
295 DO ipt=1,npt
296 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
297 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1)+i)/npt
298 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2)+i)/npt
299 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3)+i)/npt
300 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4)+i)/npt
301 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5)+i)/npt
302 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6)+i)/npt
303 ENDDO
304 ENDDO
305 ENDIF
306 ELSEIF ((isolnod == 6.OR.isolnod == 8) .AND. jhbe == 15) THEN
307 IF (mlw>=28) THEN
308 DO i=lft,llt
309 n = i + nft
310 DO il= 1,nlay
311 DO ipt=1,nptg
312 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
313 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1)+i)/nptg
314 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2)+i)/nptg
315 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3)+i)/nptg
316 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4)+i)/nptg
317 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5)+i)/nptg
318 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6)+i)/nptg
319 ENDDO
320 ENDDO
321 ENDDO
322 ELSEIF(mlw == 14)THEN
323 DO i=lft,llt
324 n = i + nft
325 DO il= 1,nlay
326 DO ipt=1,nptg
327 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
328 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1)+i)/nptg
329 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2)+i)/nptg
330 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3)+i)/nptg
331 ENDDO
332 ENDDO
333 ENDDO
334 ELSEIF (mlw == 24) THEN
335 DO i=lft,llt
336 n = i + nft
337 DO il= 1,nlay
338 DO ipt=1,nptg
339 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
340 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1)+i)/nptg
341 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2)+i)/nptg
342 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3)+i)/nptg
343 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4)+i)/nptg
344 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5)+i)/nptg
345 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6)+i)/nptg
346 ENDDO
347 ENDDO
348 ENDDO
349 ENDIF
350 ENDIF
351C-----------------------------------------------
352 ELSEIF (itens == 4.AND.mlw == 24.AND.
353 . nint(pm(56,mt1)) == 1) THEN
354C CRACKS
355C-----------------------------------------------
356 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
357 DO i=lft,llt
358 evar(1,i) = lbuf%DGLO(jj(1)+i)
359 evar(2,i) = lbuf%DGLO(jj(2)+i)
360 evar(3,i) = lbuf%DGLO(jj(3)+i)
361 evar(4,i) = lbuf%DGLO(jj(4)+i)
362 evar(5,i) = lbuf%DGLO(jj(5)+i)
363 evar(6,i) = lbuf%DGLO(jj(6)+i)
364 ENDDO
365 IF (kcvt/=0) THEN
366C DAMAGE IN GLOBAL SYSTEM
367 DO i=lft,llt
368 n = i + nft
369 IF(el2fa(nn2+n)/=0)THEN
370 IF(kcvt==2)THEN
371 gama(1)= gbuf%GAMA(jj(1) + i)
372 gama(2)= gbuf%GAMA(jj(2) + i)
373 gama(3)= zero
374 gama(4)=-gama(2)
375 gama(5)= gama(1)
376 gama(6)= zero
377 ELSE
378 gama(1)=one
379 gama(2)=zero
380 gama(3)=zero
381 gama(4)=zero
382 gama(5)=one
383 gama(6)=zero
384 END IF
385 CALL srota6(x,ixs(1:nixs,n),kcvt,evar(1:6,i),gama)
386 ENDIF
387 ENDDO
388 ENDIF
389C--------
390 ELSE
391C--------
392 DO i=lft,llt
393 n = i + nft
394 evar(1,i) = zero
395 evar(2,i) = zero
396 evar(3,i) = zero
397 evar(4,i) = zero
398 evar(5,i) = zero
399 evar(6,i) = zero
400 ENDDO
401 ENDIF
402C-----------------------------------------------
403 IF(isolnod == 16)THEN
404 DO i=lft,llt
405 n = i + nft
406 IF(el2fa(nn2+n)/=0)THEN
407 tens(1,el2fa(nn2+n)) = evar(1,i)
408 tens(2,el2fa(nn2+n)) = evar(2,i)
409 tens(3,el2fa(nn2+n)) = evar(3,i)
410 tens(4,el2fa(nn2+n)) = evar(4,i)
411 tens(5,el2fa(nn2+n)) = evar(5,i)
412 tens(6,el2fa(nn2+n)) = evar(6,i)
413 tens(1,el2fa(nn2+n)+1) = evar(1,i)
414 tens(2,el2fa(nn2+n)+1) = evar(2,i)
415 tens(3,el2fa(nn2+n)+1) = evar(3,i)
416 tens(4,el2fa(nn2+n)+1) = evar(4,i)
417 tens(5,el2fa(nn2+n)+1) = evar(5,i)
418 tens(6,el2fa(nn2+n)+1) = evar(6,i)
419 tens(1,el2fa(nn2+n)+2) = evar(1,i)
420 tens(2,el2fa(nn2+n)+2) = evar(2,i)
421 tens(3,el2fa(nn2+n)+2) = evar(3,i)
422 tens(4,el2fa(nn2+n)+2) = evar(4,i)
423 tens(5,el2fa(nn2+n)+2) = evar(5,i)
424 tens(6,el2fa(nn2+n)+2) = evar(6,i)
425 tens(1,el2fa(nn2+n)+3) = evar(1,i)
426 tens(2,el2fa(nn2+n)+3) = evar(2,i)
427 tens(3,el2fa(nn2+n)+3) = evar(3,i)
428 tens(4,el2fa(nn2+n)+3) = evar(4,i)
429 tens(5,el2fa(nn2+n)+3) = evar(5,i)
430 tens(6,el2fa(nn2+n)+3) = evar(6,i)
431 ENDIF
432 ENDDO
433 ELSE
434 DO i=lft,llt
435 n = i + nft
436 IF(el2fa(nn2+n)/=0)THEN
437 tens(1,el2fa(nn2+n)) = evar(1,i)
438 tens(2,el2fa(nn2+n)) = evar(2,i)
439 tens(3,el2fa(nn2+n)) = evar(3,i)
440 tens(4,el2fa(nn2+n)) = evar(4,i)
441 tens(5,el2fa(nn2+n)) = evar(5,i)
442 tens(6,el2fa(nn2+n)) = evar(6,i)
443 ENDIF
444 ENDDO
445 ENDIF
446C-----------------------------------------------
447 ELSEIF(isph3d == 1.AND.ity == 51)THEN
448C-----------------------------------------------
449C TETRAS SPH.
450C-----------------------------------------------
451 iprt=ipartsp(1 + nft)
452 mt1 =ipart(1,iprt)
453C-----------------------------------------------
454 IF (itens == 1) THEN
455C STRESS
456C-----------------------------------------------
457 DO i=lft,llt
458 n = i + nft
459 IF(el2fa(nn3+n)/=0)THEN
460 tens(1,el2fa(nn3+n)) = lbuf%SIG(jj(1) + i)
461 tens(2,el2fa(nn3+n)) = lbuf%SIG(jj(2) + i)
462 tens(3,el2fa(nn3+n)) = lbuf%SIG(jj(3) + i)
463 tens(4,el2fa(nn3+n)) = lbuf%SIG(jj(4) + i)
464 tens(5,el2fa(nn3+n)) = lbuf%SIG(jj(5) + i)
465 tens(6,el2fa(nn3+n)) = lbuf%SIG(jj(6) + i)
466 ENDIF
467 ENDDO
468C-----------------------------------------------
469 ELSEIF(itens == 4.AND.mlw == 24.
470 . and.nint(pm(56,mt1)) == 1)THEN
471C CRACKS
472C-----------------------------------------------
473 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
474 DO i=lft,llt
475 n = i + nft
476 IF (el2fa(nn3+n)/=0) THEN
477 evar(1,i) = lbuf%DGLO(jj(1)+i)
478 evar(2,i) = lbuf%DGLO(jj(2)+i)
479 evar(3,i) = lbuf%DGLO(jj(3)+i)
480 evar(4,i) = lbuf%DGLO(jj(4)+i)
481 evar(5,i) = lbuf%DGLO(jj(5)+i)
482 evar(6,i) = lbuf%DGLO(jj(6)+i)
483 ENDIF
484 ENDDO
485 ELSE
486 DO i=lft,llt
487 n = i + nft
488 IF(el2fa(nn3+n)/=0)THEN
489 tens(1,el2fa(nn3+n)) = zero
490 tens(2,el2fa(nn3+n)) = zero
491 tens(3,el2fa(nn3+n)) = zero
492 tens(4,el2fa(nn3+n)) = zero
493 tens(5,el2fa(nn3+n)) = zero
494 tens(6,el2fa(nn3+n)) = zero
495 ENDIF
496 ENDDO
497 ENDIF
498C-----------------------------------------------
499 ENDIF ! ITY
500C-----------------------------------------------
501 490 CONTINUE
502C-----------------------------------------------
503 DO n=1,nbf
504 r4(1) = tens(1,n)
505 r4(2) = tens(2,n)
506 r4(3) = tens(3,n)
507 r4(4) = tens(4,n)
508 r4(5) = tens(5,n)
509 r4(6) = tens(6,n)
510 CALL write_r_c(r4,6)
511 ENDDO
512C-----------
513 RETURN
514 END
my_real
#define my_real
Definition cppsort.cpp:32
genani1
subroutine genani1(x, bufel, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, swaft, iparg, pm, geo, skew, itab, lpby, npby, nstrf, rwbuf, nprw, ipart, iparts, ipartq, ipartc, ipartt, ipartp, ipartr, iparttg, rby, swa4, igrsurf, bufsf, ipartx, kxsp, ixsp, ipartsp, spbuf, ixs10, ixs20, ixs16, ipm, igeo, smater, sel2fa, snfacptx, sixedge, soffx1, snumx1, sxnorm, sinvert, sfunc1, siad, nmanim, d, smas, ms, fxani, mbufel, mdepl, nlevel, elsub, dsanim, nelem, cep, cepsp, nom_opt, ptr_nopt_rwall, ptr_nopt_sect, elbuf_tab, sph2sol, subset)
Definition genani1.F:98
initbuf
subroutine initbuf(iparg, ng, mtn, llt, nft, iad, ity, npt, jale, ismstr, jeul, jtur, jthe, jlag, jmult, jhbe, jivf, jpor, jcvt, jclose, jpla, irep, iint, igtyp, israt, isrot, icsen, isorth, isorthg, ifailure)
Definition initbuf.F:38
srota6
subroutine srota6(x, ixs, kcvt, tens, gama)
Definition srota6.F:33
starter
program starter
Definition starter.F:39
tensors
subroutine tensors(elbuf_tab, iparg, itens, ixs, pm, el2fa, nbf, tens, x, ipart, ipartsp, isph3d, ipm)
Definition tensors.F:35
write_r_c
void write_r_c(float *w, int *len)
Definition write_routines.c:519