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