38 SUBROUTINE tensors(ELBUF_TAB,IPARG ,ITENS ,IXS ,PM ,
39 2 EL2FA ,NBF ,TENS ,EPSDOT ,
40 3 NBPART ,X ,IADG ,IPART ,
41 4 IPARTSP ,ISPH3D ,IPM ,IGEO )
51#include
"implicit_f.inc"
55#include
"vect01_c.inc"
67 my_real tens(6,*),epsdot(6,*),pm(npropm,*),x(3,*)
68 INTEGER IPARG(NPARG,*),ITENS,
69 . IXS(NIXS,*),EL2FA(*),IADG(NSPMD,*),IPM(NPROPMI,*),
70 . nbf,nbpart,ipart(lipart1,*),ipartsp(*),
71 . isph3d,igeo(npropgi,*)
72 TYPE (ELBUF_STRUCT_),
DIMENSION(NGROUP),
TARGET :: ELBUF_TAB
76 my_real off, fac, a1, a2, a3, thk, gama
78 INTEGER I,I1,I2,II,N,J,NG,,IPT,MT1,MLW, ISTRAIN,TSHELL,
79 . ipid, ns1, ns2 ,ialel, istre,iprt,pti,iadpg,pid,
80 . nn1,nn2,nn3,nn4,icsig,ior_tsh,nuvar,buf,l_pla,l_stra,khbe,
81 . kcvt,isolnod,nlay,nptr,npts,nptt,nptg,il,is,ir,it,ivisc,iok,
83 REAL,
DIMENSION(:),
ALLOCATABLE :: WA
84 TYPE(G_BUFEL_) ,
POINTER :: GBUF
85 TYPE(L_BUFEL_) ,
POINTER :: LBUF
88 CALL my_alloc(wa,6*nbf)
95 nn4 = nn3 + isph3d*(numsph+maxpjet)
98 gbuf => elbuf_tab(ng)%GBUF
99 istrain = iparg(44,ng)
100 isolnod = iparg(28,ng)
103 2 mlw ,nel ,nft ,iad ,ity ,
104 3 npt ,jale ,ismstr ,jeul ,jtur ,
105 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
106 5 nvaux ,jpor ,kcvt ,jclose ,jplasol ,
107 6 irep ,iint ,igtyp ,israt ,isrot ,
108 7 icsen ,isorth ,isorthg ,ifailure,jsms )
123 IF (igtyp==20 .OR. igtyp==21 .OR. igtyp==22) tshell = 1
124 IF (igtyp == 21.OR.igtyp == 22) ior_tsh = 1
125 nlay = elbuf_tab(ng)%NLAY
126 nptr = elbuf_tab(ng)%NPTR
127 npts = elbuf_tab(ng)%NPTS
128 nptt = elbuf_tab(ng)%NPTT
129 nptg = nptt*npts*nptr
134 IF (kcvt==1.AND.isorth/=0) kcvt=2
136 IF (igtyp /= 22)
THEN
137 IF (isorth > 0) isorthg = 0
142 tens(1,el2fa(nn2+n)) = zero
143 tens(2,el2fa(nn2+n)) = zero
144 tens(3,el2fa(nn2+n)) = zero
145 tens(4,el2fa(nn2+n)) = zero
146 tens(5,el2fa(nn2+n)) = zero
147 tens(6,el2fa(nn2+n)) = zero
151 evar(1:6,lft:llt)=zero
158 evar(1,i) = gbuf%SIG(jj(1) + i)
159 evar(2,i) = gbuf%SIG(jj(2) + i)
160 evar(3,i) = gbuf%SIG(jj(3) + i)
161 evar(4,i) = gbuf%SIG(jj(4) + i)
162 evar(5,i) = gbuf%SIG(jj(5) + i)
163 evar(6,i) = gbuf%SIG(jj(6) + i)
166 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
168 evar(1,i) =evar(1,i)+ lbuf%VISC(jj(1) + i)
169 evar(2,i) =evar(2,i)+ lbuf%VISC(jj(2) + i)
170 evar(3,i) =evar(3,i)+ lbuf%VISC(jj(3) + i)
171 evar(4,i) =evar(4,i)+ lbuf%VISC(jj(4) + i)
172 evar(5,i) =evar(5,i)+ lbuf%VISC(jj(5) + i)
173 evar(6,i) =evar(6,i)+ lbuf%VISC(jj(6) + i)
177 IF( nfilsol /= 0 .AND. gbuf%G_FILL /= 0 )
THEN
179 evar(1,i) = evar(1,i) * gbuf%FILL(i)
180 evar(2,i) = evar(2,i) * gbuf%FILL(i)
181 evar(3,i) = evar(3,i) * gbuf%FILL(i)
182 evar(4,i) = evar(4,i) * gbuf%FILL(i)
183 evar(5,i) = evar(5,i) * gbuf%FILL(i)
184 evar(6,i) = evar(6,i) * gbuf%FILL(i)
188 IF (jhbe == 17 .AND. iint ==3)
THEN
189! stress tensor in global system
192 IF(el2fa(nn2+n) /= 0)
THEN
194 IF(kcvt==2.AND.jhbe/=14.AND.jhbe/=15)
THEN
195 gama(1)=gbuf%GAMA(jj(1) + i)
196 gama(2)=gbuf%GAMA(jj(2) + i)
197 gama(3)=gbuf%GAMA(jj(3) + i)
198 gama(4)=gbuf%GAMA(jj(4) + i)
199 gama(5)=gbuf%GAMA(jj(5) + i)
200 gama(6)=gbuf%GAMA(jj(6) + i)
210 2 igtyp, gbuf%COR_FR(9*(i-1)+1),iint, isorth)
213 ELSE IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
217 IF(el2fa(nn2+n) /= 0)
THEN
219 IF(kcvt==2.AND.jhbe/=14.AND.jhbe/=15)
THEN
220 gama(1)=gbuf%GAMA(jj(1) + i)
221 gama(2)=gbuf%GAMA(jj(2) + i)
222 gama(3)=gbuf%GAMA(jj(3) + i)
223 gama(4)=gbuf%GAMA(jj(4) + i)
224 gama(5)=gbuf%GAMA(jj(5) + i)
225 gama(6)=gbuf%GAMA(jj(6) + i)
234 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
239 ELSEIF (itens == 2)
THEN
244 IF (isolnod == 8 .AND. igtyp == 43)
THEN
247 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
248 evar(3,i) = evar(3,i) + lbuf%EPE(jj(1) + i)/npt
249 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)/npt
250 evar(1,i) = evar(1,i) + lbuf%EPE(jj(3) + i)/npt
255 IF(el2fa(nn2+n) /= 0)
THEN
256 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
260 ELSEIF (isolnod == 8 .AND. npt == 8 .AND. jhbe /= 14 .AND. jhbe /= 24 .AND. jhbe /= 15 .AND. jhbe /= 17 )
THEN
266 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,j)
267 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)*one_over_8
268 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)*one_over_8
269 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)*one_over_8
270 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*one_over_8
271 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*one_over_8
272 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*one_over_8
277 ELSEIF ((isolnod==8 .OR. (isolnod == 4 .AND. isrot==0)) .AND. npt==1 .AND. jhbe /= 14 .AND. jhbe /= 15)
THEN
278 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
279 IF (isorth > 0) isorthg = 1
280 IF (mlw>=28.AND.mlw /= 49)
THEN
283 evar(1,i) = lbuf%STRA(jj(1) + i)
284 evar(2,i) = lbuf%STRA(jj(2) + i)
285 evar(3,i) = lbuf%STRA(jj(3) + i)
286 evar(4,i) = lbuf%STRA(jj(4) + i)*half
287 evar(5,i) = lbuf%STRA(jj(5) + i)*half
288 evar(6,i) = lbuf%STRA(jj(6) + i)*half
290 IF (isorth > 0) kcvt = 2
291 ELSEIF (mlw == 12 .OR. mlw == 14)
THEN
294 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)
295 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)
296 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)
298 IF (isorth > 0) kcvt = 2
299 ELSEIF (mlw == 24 .OR. mlw == 25)
THEN
302 evar(1,i) = lbuf%STRA(jj(1) + i)
303 evar(2,i) = lbuf%STRA(jj(2) + i)
304 evar(3,i) = lbuf%STRA(jj(3) + i)
305 evar(4,i) = lbuf%STRA(jj(4) + i)*half
306 evar(5,i) = lbuf%STRA(jj(5) + i)*half
307 evar(6,i) = lbuf%STRA(jj(6) + i)*half
309 IF (isorth > 0) kcvt = 2
310 ELSEIF (istrain > 0)
THEN
311 IF (mlw /= 14 .AND. mlw /= 24 .AND. mlw<28 .OR. mlw == 49)
THEN
314 evar(1,i) = lbuf%STRA(jj(1) + i)
315 evar(2,i) = lbuf%STRA(jj(2) + i)
316 evar(3,i) = lbuf%STRA(jj(3) + i)
317 evar(4,i) = lbuf%STRA(jj(4) + i)*half
318 evar(5,i) = lbuf%STRA(jj(5) + i)*half
319 evar(6,i) = lbuf%STRA(jj(6) + i)*half
336 IF(el2fa(nn2+n) /= 0)
THEN
338 gama(1)=gbuf%GAMA(jj(1) + i)
339 gama(2)=gbuf%GAMA(jj(2) + i)
340 gama(3)=gbuf%GAMA(jj(3) + i)
341 gama(4)=gbuf%GAMA(jj(4) + i)
342 gama(5)=gbuf%GAMA(jj(5) + i)
343 gama(6)=gbuf%GAMA(jj(6) + i)
352 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
357 ELSEIF(isolnod == 16 .OR. isolnod == 20 .OR. (isolnod == 8 .AND. (jhbe == 14 .OR. jhbe == 17)))
THEN
359 IF (mlw>=28.AND.mlw /= 49)
THEN
366 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
367 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)/npt
368 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)/npt
369 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3)
370 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*half/npt
371 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*half/npt
372 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*half/npt
378 ELSEIF (mlw == 12 .OR. mlw == 14)
THEN
385 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
386 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)/npt
387 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)/npt
388 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)/npt
394 ELSEIF(mlw == 24 .OR. mlw == 25)
THEN
401 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
402 IF (elbuf_tab(ng)%BUFLY(il)%L_STRA > 0)
THEN
403 evar_tmp(1) = lbuf%STRA(jj
404 evar_tmp(2) = lbuf%STRA(jj(2) + i)/npt
405 evar_tmp(3) = lbuf%STRA(jj(3) + i)/npt
406 evar_tmp(4) = lbuf%STRA(jj(4) + i)*half/npt
407 evar_tmp(5) = lbuf%STRA(jj(5) + i)*half/npt
408 evar_tmp(6) = lbuf%STRA(jj(6) + i)*half/npt
410 IF (kcvt /= 0 .AND.icsig > 0)
THEN
411 IF (igtyp == 21)
THEN
412! strain tensor in global system
416 IF(el2fa(nn2+n) /= 0)
THEN
419 gama(2)= gbuf%GAMA(jj(1) + i)
420 gama(3)= gbuf%GAMA(jj(2) + i)
432 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
435 IF(el2fa(nn2+n) /= 0)
THEN
437 gama(1)= gbuf%GAMA(jj(1) + i)
438 gama(2)= gbuf%GAMA(jj(2) + i)
451 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
454 IF(el2fa(nn2+n) /= 0)
THEN
456 gama(1)= gbuf%GAMA(jj(2) + i)
458 gama(3)= gbuf%GAMA(jj(1) + i)
470 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
479 IF(el2fa(nn2+n) /= 0)
THEN
482 gama(2)= lbuf%GAMA(jj(1) + i)
483 gama(3)= lbuf%GAMA(jj(2) + i)
495 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
498 IF(el2fa(nn2+n) /= 0)
THEN
500 gama(1)= lbuf%GAMA(jj(1) + i)
501 gama(2)= lbuf%GAMA(jj(2) + i)
514 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth
517 IF(el2fa(nn2+n) /= 0)
THEN
519 gama(1)= lbuf%GAMA(jj(2) + i)
521 gama(3)= lbuf%GAMA(jj(1) + i)
533 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
539 evar(1,i) = evar(1,i)+evar_tmp(1)
540 evar(2,i) = evar(2,i)+evar_tmp(2)
541 evar(3,i) = evar(3,i)+evar_tmp(3)
542 evar(4,i) = evar(4,i)+evar_tmp(4)
543 evar(5,i) = evar(5,i)+evar_tmp(5)
544 evar(6,i) = evar(6,i)+evar_tmp(6)
551 ELSEIF(istrain > 0)
THEN
552 IF (mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
560 evar_tmp(1) = lbuf%STRA(jj(1) + i
561 evar_tmp(2) = lbuf%STRA(jj
562 evar_tmp(3) = lbuf%STRA(jj(3) + i)/npt
563 evar_tmp(4) = lbuf%STRA(jj(4) + i)*half/npt
564 evar_tmp(5) = lbuf%STRA(jj(5) + i)*half/npt
565 evar_tmp(6) = lbuf%STRA(jj(6) +
567 IF (kcvt /= 0 .AND.icsig > 0)
THEN
572 IF(el2fa(nn2+n) /= 0)
THEN
575 gama(2)= lbuf%GAMA(jj(1) + i)
576 gama(3)= lbuf%GAMA(jj(2) + i)
588 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
591 IF(el2fa(nn2+n) /= 0)
THEN
593 gama(1)= lbuf%GAMA(jj(1) + i)
594 gama(2)= lbuf%GAMA(jj(2) + i)
607 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
610 IF(el2fa(nn2+n) /= 0)
THEN
612 gama(1)= lbuf%GAMA(jj(2) + i)
614 gama(3)= lbuf%GAMA(jj(1) + i)
626 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
631 evar(1,i) = evar(1,i)+evar_tmp(1)
632 evar(2,i) = evar(2,i)+evar_tmp(2)
633 evar(3,i) = evar(3,i)+evar_tmp(3)
634 evar(4,i) = evar(4,i)+evar_tmp(4)
635 evar(5,i) = evar(5,i)+evar_tmp(5)
636 evar(6,i) = evar(6,i)+evar_tmp(6)
655 IF (mlw == 12 .OR. mlw == 14 .OR. mlw == 24 .OR. mlw == 25 .OR. (mlw >= 28 .AND. mlw /= 49))
THEN
656 IF (isorth > 0) kcvt = 2
659 IF (kcvt /= 0 .AND.icsig == 0 .AND. jhbe /= 16)
THEN
663 IF(el2fa(nn2+n) /= 0)
THEN
665 gama(1)=gbuf%GAMA(jj(1) + i)
666 gama(2)=gbuf%GAMA(jj(2) + i)
667 gama(3)=gbuf%GAMA(jj(3) + i)
668 gama(4)=gbuf%GAMA(jj(4) + i)
669 gama(5)=gbuf%GAMA(jj(5) + i)
670 gama(6)=gbuf%GAMA(jj(6) + i)
679 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
684 ELSEIF (isolnod==10 .OR. (isolnod==4 .AND. isrot==1))
THEN
686 IF (mlw>=28.AND.mlw /= 49)
THEN
690 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
691 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/npt
692 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/npt
693 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/npt
694 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4) + i)*half/npt
695 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/npt
696 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/npt
699 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
703 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
704 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)/npt
705 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)/npt
706 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)/npt
709 ELSEIF ((mlw == 24 .OR. mlw == 25) .and. istrain > 0)
THEN
713 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
714 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/npt
715 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/npt
716 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/npt
717 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4) + i)*half/npt
718 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/npt
719 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/npt
722 ELSEIF(istrain > 0)
THEN
723 IF (mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
727 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
728 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/npt
729 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/npt
730 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/npt
731 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4) + i)*half/npt
732 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/npt
733 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/npt
751 IF (el2fa(nn2+n) /= 0)
THEN
753 gama(1)=gbuf%GAMA(jj(1) + i)
754 gama(2)=gbuf%GAMA(jj(2) + i)
755 gama(3)=gbuf%GAMA(jj(3) + i)
756 gama(4)=gbuf%GAMA(jj(4) + i)
757 gama(5)=gbuf%GAMA(jj(5) + i)
758 gama(6)=gbuf%GAMA(jj(6) + i)
767 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp
772 ELSEIF((isolnod == 6.OR.isolnod == 8).AND.jhbe == 15)
THEN
774 IF (mlw>=28.AND.mlw /= 49.AND.istrain > 0)
THEN
779 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
780 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/(nptg*nlay)
781 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/(nptg*nlay)
782 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/(nptg*nlay)
784 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/(nptg*nlay)
785 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/(nptg*nlay)
789 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
793 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
794 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)/(nptg*nlay)
795 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)/(nptg*nlay)
796 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)/(nptg*nlay)
800 ELSEIF ((mlw == 24 .OR. mlw == 25) .and. istrain > 0)
THEN
805 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
806 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/(nptg*nlay)
807 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/(nptg*nlay)
808 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/(nptg*nlay)
809 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4) + i)*half/(nptg*nlay)
810 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/(nptg*nlay)
811 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/(nptg*nlay)
815 ELSEIF (istrain > 0)
THEN
816 IF(mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
821 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
822 evar(1,i) = evar(1,i)+lbuf%STRA(jj(1) + i)/(nptg*nlay)
823 evar(2,i) = evar(2,i)+lbuf%STRA(jj(2) + i)/(nptg*nlay)
824 evar(3,i) = evar(3,i)+lbuf%STRA(jj(3) + i)/(nptg*nlay)
825 evar(4,i) = evar(4,i)+lbuf%STRA(jj(4) + i)*half/(nptg*nlay)
826 evar(5,i) = evar(5,i)+lbuf%STRA(jj(5) + i)*half/(nptg*nlay)
827 evar(6,i) = evar(6,i)+lbuf%STRA(jj(6) + i)*half/(nptg*nlay)
846 IF (el2fa(nn2+n) /= 0)
THEN
848 gama(1)= gbuf%GAMA(jj(1) + i)
849 gama(2)= gbuf%GAMA(jj(2) + i)
862 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
871 ELSEIF (itens == 4 .AND. mlw == 24 .AND. nint(pm(56,mt1)) == 1)
THEN
882 IF (isolnod == 8 .AND.(jhbe == 14 .OR. jhbe == 15))
THEN
885 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
887 evar(1,i) = evar(1,i)+lbuf%DGLO(jj(1) + i)
888 evar(2,i) = evar(2,i)+lbuf%DGLO(jj(2) + i)
889 evar(3,i) = evar(3,i)+lbuf%DGLO(jj(3) + i)
890 evar(4,i) = evar(4,i)+lbuf%DGLO(jj(4) + i)
891 evar(5,i) = evar(5,i)+lbuf%DGLO(jj(5) + i)
892 evar(6,i) = evar(6,i)+lbuf%DGLO(jj(6) + i)
899 IF(el2fa(nn2+n) /= 0)
THEN
901 gama(1)= gbuf%GAMA(jj(1) + i)
902 gama(2)= gbuf%GAMA(jj(2) + i)
915 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
920 ELSEIF (itens == 5)
THEN
933 IF ((isolnod == 8 .OR. (isolnod == 4 .AND. isrot == 0)) .AND. npt == 1 .AND. jhbe /= 14 .AND. jhbe /= 15)
THEN
935 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
936 IF (isorth > 0) isorthg = 1
940 evar(1,i) = lbuf%PLA(jj(1) + i + nel)
941 evar(2,i) = lbuf%PLA(jj(2) + i + nel)
942 evar(3,i) = lbuf%PLA(jj(3) + i + nel)
943 evar(4,i) = lbuf%PLA(jj(4) + i + nel)*half
944 evar(5,i) = lbuf%PLA(jj(5) + i + nel)*half
945 evar(6,i) = lbuf%PLA(jj(6) + i + nel)*half
953 IF (el2fa(nn2+n) /= 0)
THEN
955 gama(1) = gbuf%GAMA(jj(1) + i)
956 gama(2) = gbuf%GAMA(jj(2) + i)
957 gama(3) = gbuf%GAMA(jj(3) + i)
958 gama(4) = gbuf%GAMA(jj(4) + i)
959 gama(5) = gbuf%GAMA(jj(5) + i)
960 gama(6) = gbuf%GAMA(jj(6) + i)
969 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe
974 ELSEIF (isolnod == 16 .OR. isolnod == 20 .OR. (isolnod == 8 .AND. (jhbe == 14 .OR. jhbe == 17)))
THEN
983 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ir,is,it)
984 IF (elbuf_tab(ng)%BUFLY(il)%L_PLA > 0)
THEN
986 evar_tmp(1) = lbuf%PLA(jj(1) + i + nel)/npt
987 evar_tmp(2) = lbuf%PLA(jj(2) + i + nel)/npt
988 evar_tmp(3) = lbuf%PLA(jj(3) + i + nel
989 evar_tmp(4) = lbuf%PLA(jj(4) + i + nel)*half/npt
990 evar_tmp(5) = lbuf%PLA(jj(5) + i + nel)*half/npt
991 evar_tmp(6) = lbuf%PLA(jj(6) + i + nel)*half/npt
994 IF (kcvt /= 0 .AND.icsig > 0)
THEN
995 IF (igtyp == 21)
THEN
996! plastic strain tensor in global system
1000 IF (el2fa(nn2+n) /= 0)
THEN
1003 gama(2) = gbuf%GAMA(jj(1) + i)
1004 gama(3) = gbuf%GAMA(jj(2) + i)
1016 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
1019 IF (el2fa(nn2+n) /= 0)
THEN
1021 gama(1) = gbuf%GAMA(jj(1) + i)
1022 gama(2) = gbuf%GAMA(jj(2) + i)
1035 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
1038 IF (el2fa(nn2+n) /= 0)
THEN
1040 gama(1) = gbuf%GAMA(jj(2) + i)
1042 gama(3) = gbuf%GAMA(jj(1) + i)
1054 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
1060 IF (jhbe == 14)
THEN
1063 IF (el2fa(nn2+n) /= 0)
THEN
1066 gama(2) = lbuf%GAMA(jj(1) + i)
1067 gama(3) = lbuf%GAMA(jj(2) + i)
1079 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
1082 IF (el2fa(nn2+n) /= 0)
THEN
1084 gama(1) = lbuf%GAMA(jj(1) + i)
1085 gama(2) = lbuf%GAMA(jj(2) + i)
1098 CALL srota6(x, ixs(1,n),kcvt, evar_tmp, gama, jhbe, igtyp, isorth)
1101 IF (el2fa(nn2+n) /= 0)
THEN
1103 gama(1) = lbuf%GAMA(jj(2) + i)
1117 CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
1123 evar(1,i) = evar(1,i)+evar_tmp
1124 evar(2,i) = evar(2,i)+evar_tmp(2)
1125 evar(3,i) = evar(3,i)+evar_tmp(3)
1126 evar(4,i) = evar(4,i)+evar_tmp(4)
1127 evar(5,i) = evar(5,i)+evar_tmp(5)
1128 evar(6,i) = evar(6,i)+evar_tmp(6)
1137 icsig = iparg(17,ng)
1138 IF (kcvt /= 0 .AND. icsig == 0 .AND. jhbe /= 16)
THEN
1142 IF (el2fa(nn2+n) /= 0)
THEN
1144 gama(1) = gbuf%GAMA(jj(1) + i)
1145 gama(2) = gbuf%GAMA(jj(2) + i)
1146 gama(3) = gbuf%GAMA(jj(3) + i)
1147 gama(4) = gbuf%GAMA(jj(4) + i)
1148 gama(5) = gbuf%GAMA(jj(5) + i)
1149 gama(6) = gbuf%GAMA(jj(6) + i)
1158 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1163 ELSEIF (isolnod == 10 .OR. (isolnod == 4 .AND. isrot == 1))
THEN
1165 IF (mlw == 24 .AND. istrain > 0)
THEN
1169 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
1170 evar(1,i) = evar(1,i) + lbuf%PLA(jj(1) + i + nel)/npt
1171 evar(2,i) = evar(2,i) + lbuf%PLA(jj(2) + i + nel)/npt
1172 evar(3,i) = evar(3,i) + lbuf%PLA(jj(3) + i + nel)/npt
1173 evar(4,i) = evar(4,i) + lbuf%PLA(jj(4) + i + nel)*half/npt
1174 evar(5,i) = evar(5,i) + lbuf%PLA(jj(5) + i + nel)*half/npt
1175 evar(6,i) = evar(6,i) + lbuf%PLA(jj(6) + i + nel)*half/npt
1184 IF (el2fa(nn2+n) /= 0)
THEN
1186 gama(1) = gbuf%GAMA(jj(1) + i)
1187 gama(2) = gbuf%GAMA(jj(2) + i)
1188 gama(3) = gbuf%GAMA(jj(3) + i)
1189 gama(4) = gbuf%GAMA(jj(4) + i)
1190 gama(5) = gbuf%GAMA(jj(5) + i)
1191 gama(6) = gbuf%GAMA(jj(6) + i)
1200 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1205 ELSEIF ((isolnod == 6 .OR. isolnod == 8) .AND. jhbe == 15)
THEN
1207 IF (mlw == 24 .AND. istrain > 0)
THEN
1212 lbuf => elbuf_tab(ng)%BUFLY(il)%LBUF(ipt,1,1)
1213 evar(1,i) = evar(1,i)+lbuf%PLA(jj(1) + i + nel)/(nptg*nlay)
1214 evar(2,i) = evar(2,i)+lbuf%PLA(jj(2) + i + nel)/(nptg*nlay)
1215 evar(3,i) = evar(3,i)+lbuf%PLA(jj(3) + i + nel)/(nptg*nlay)
1216 evar(4,i) = evar(4,i)+lbuf%PLA(jj(4) + i + nel)*half/(nptg*nlay)
1217 evar(5,i) = evar(5,i)+lbuf%PLA(jj(5) + i + nel)*half/(nptg*nlay)
1218 evar(6,i) = evar(6,i)+lbuf%PLA(jj(6) + i + nel)*half/(nptg*nlay)
1228 IF (el2fa(nn2+n) /= 0)
THEN
1230 gama(1) = gbuf%GAMA(jj(1) + i)
1231 gama(2) = gbuf%GAMA(jj(2) + i)
1244 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1245 ENDIF !
IF (el2fa(nn2+n) /= 0)
1255 ELSEIF (itens>=10.AND.itens<=10
THEN
1258 IF (isolnod == 8 .AND. igtyp == 43)
THEN
1263 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
1264 evar(3,i) = evar(3,i) + lbuf%SIG(jj(3) + i)/nptr
1265 evar(2,i) = evar(2,i) + lbuf%SIG(jj(5) + i)/nptr
1266 evar(1,i) = evar(1,i) + lbuf%SIG(jj(6) + i)/nptr
1272 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
1273 evar(3,i)= evar(3,i)+ lbuf%SIG(jj(3) + i)/nptr+ lbuf%VISC(jj(3) + i)/nptr
1274 evar(2,i)= evar(2,i)+ lbuf%SIG
1275 evar(1,i)= evar(1,i)+ lbuf%SIG(jj(6) + i)/nptr+ lbuf%VISC(jj(6) + i)/nptr
1281 IF(el2fa(nn2+n) /= 0)
THEN
1282 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1286 ELSEIF (isolnod == 8 .AND. npt == 8.AND. jhbe /= 14 .AND. jhbe /= 24 .AND. jhbe /= 15)
THEN
1289 is = mod(abs(pti)/10,10)
1290 it = mod(abs(pti),10)
1291 IF (ir == 0 .AND. it =
THEN
1293 ELSEIF(ir <= nptr .AND. is <= npts .AND. it <= nptt)
THEN
1294 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1295 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
1298 evar(1,i) = lbuf%SIG(jj(1) + i)
1299 evar(2,i) = lbuf%SIG(jj(2) + i)
1300 evar(3,i) = lbuf%SIG
1301 evar(4,i) = lbuf%SIG(jj(4) + i)
1302 evar(5,i) = lbuf%SIG(jj(5) + i)
1303 evar(6,i) = lbuf%SIG(jj(6) + i)
1307 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1308 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1309 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1310 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1311 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1312 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1320 IF(el2fa(nn2+n) /= 0)
THEN
1322 gama(1)= gbuf%GAMA(jj(1) + i)
1323 gama(2)= gbuf%GAMA(jj(2) + i)
1324 gama(3)= gbuf%GAMA(jj(3) + i)
1325 gama(4)= gbuf%GAMA(jj(4) + i)
1326 gama(5)= gbuf%GAMA(jj(5) + i)
1327 gama(6)= gbuf%GAMA(jj(6) + i)
1336 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1351 ELSEIF((isolnod == 8.OR.npt == 1) .AND. jhbe /= 14.AND.jhbe /= 15.AND.jhbe /= 17)
THEN
1358 it = mod(abs(pti),10)
1359 IF (ir == 0 .AND. it == 0)
THEN
1361 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1363 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
1365 evar(1,i) = lbuf%SIG(jj(1) + i)
1366 evar(2,i) = lbuf%SIG(jj(2) + i)
1367 evar(3,i) = lbuf%SIG(jj(3) + i)
1368 evar(4,i) = lbuf%SIG(jj(4) + i)
1369 evar(5,i) = lbuf%SIG(jj(5) + i)
1370 evar(6,i) = lbuf%SIG(jj(6) + i)
1374 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1376 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1377 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1378 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1379 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1387 IF(el2fa(nn2+n) /= 0)
THEN
1389 gama(1)=gbuf%GAMA(jj(1) + i)
1390 gama(2)=gbuf%GAMA(jj(2) + i)
1391 gama(3)=gbuf%GAMA(jj(3) + i)
1392 gama(4)=gbuf%GAMA(jj(4) + i)
1393 gama(5)=gbuf%GAMA(jj(5) + i)
1394 gama(6)=gbuf%GAMA(jj(6) + i)
1403 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1409 ELSEIF (isolnod == 20 .OR. isolnod == 16 )
THEN
1412 is=mod(abs(pti)/10,10)
1414 IF (ir == 0 .OR. is == 0.OR. it == 0) cycle
1415 IF (tshell == 1 .AND. is <= nlay )
THEN
1416 lbuf => elbuf_tab(ng)%BUFLY(is)%LBUF(ir,1,it)
1418 ELSEIF(ir <= nptr .AND. is <= npts .AND. it <= nptt)
THEN
1419 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
1422 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1425 evar(1,i) = lbuf%SIG(jj(1) + i)
1426 evar(2,i) = lbuf%SIG(jj(2) + i)
1427 evar(3,i) = lbuf%SIG(jj(3) + i)
1428 evar(4,i) = lbuf%SIG(jj(4) + i)
1429 evar(5,i) = lbuf%SIG(jj(5) + i)
1430 evar(6,i) = lbuf%SIG(jj(6) + i)
1434 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1435 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1436 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1437 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1438 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1439 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1443 IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
1447 IF(el2fa(nn2+n) /= 0)
THEN
1449 gama(1)=gbuf%GAMA(jj(1) + i)
1450 gama(2)=gbuf%GAMA(jj(2) + i)
1451 gama(3)=gbuf%GAMA(jj(3) + i)
1452 gama(4)=gbuf%GAMA(jj(4) + i
1453 gama(5)=gbuf%GAMA(jj(5) + i)
1454 gama(6)=gbuf%GAMA(jj(6) + i)
1463 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1468 ELSEIF (isolnod == 8 .AND. jhbe == 14 )
THEN
1470 icsig = iparg(17,ng)
1471 nptg = nptr * npts * nptt * nlay
1473 is0=mod(abs(pti)/10,10)
1474 it0=mod(abs(pti),10)
1475 ipid = ixs(10,1 + nft)
1476 IF (ir0==0.OR.is0==0.OR.it0==0) cycle
1480 IF (tshell == 1)
THEN
1481 IF (icsig==100)
THEN
1485 ELSEIF (icsig==10)
THEN
1496 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1498 IF (tshell == 1 .AND. it <= nlay )
THEN
1499 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
1501 ELSEIF(ir0 <= nptr .AND. is0 <= npts .AND. it0 <= nptt)
THEN
1502 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
1507 evar(1,i) = lbuf%SIG(jj(1) + i)
1508 evar(2,i) = lbuf%SIG(jj(2) + i)
1509 evar(3,i) = lbuf%SIG(jj(3) + i)
1510 evar(4,i) = lbuf%SIG(jj(4) + i)
1511 evar(5,i) = lbuf%SIG(jj(5) + i)
1512 evar(6,i) = lbuf%SIG(jj(6) + i)
1516 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1517 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1518 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1519 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1520 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1521 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1529 IF (igtyp == 21)
THEN
1534 IF(el2fa(nn2+n) /= 0)
THEN
1537 gama(2)=gbuf%GAMA(jj(1) + i)
1538 gama(3)=gbuf%GAMA(jj(2) + i)
1550 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1556 IF(el2fa(nn2+n) /= 0)
THEN
1558 gama(1)=gbuf%GAMA(jj(1) + i)
1559 gama(2)=gbuf%GAMA(jj(2) + i)
1572 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1578 IF(el2fa(nn2+n) /= 0)
THEN
1580 gama(1)=gbuf%GAMA(jj(2) + i)
1582 gama(3)=gbuf%GAMA(jj(1) + i)
1594 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1603 IF(el2fa(nn2+n) /= 0)
THEN
1606 gama(2)=lbuf%GAMA(jj(1) + i)
1607 gama(3)=lbuf%GAMA(jj(2) + i)
1619 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1625 IF(el2fa(nn2+n) /= 0)
THEN
1627 gama(1)=lbuf%GAMA(jj(1) + i)
1628 gama(2)=lbuf%GAMA(jj(2) + i)
1641 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1647 IF(el2fa(nn2+n) /= 0)
THEN
1649 gama(1)=lbuf%GAMA(jj(2) + i)
1651 gama(3)=lbuf%GAMA(jj(1) + i)
1663 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1671 IF(el2fa(nn2+n) /= 0)
THEN
1673 gama(1)=gbuf%GAMA(jj(1) + i)
1674 gama(2)=gbuf%GAMA(jj(2) + i)
1675 gama(3)=gbuf%GAMA(jj(3) + i)
1676 gama(4)=gbuf%GAMA(jj(4) + i)
1677 gama(5)=gbuf%GAMA(jj(5) + i)
1678 gama(6)=gbuf%GAMA(jj(6) + i)
1687 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1693 ELSEIF(isolnod == 10.OR.(isolnod == 4 .AND. isrot == 1))
THEN
1696 is=mod(abs(pti)/10,10)
1698 IF (ir == 0 .AND. it == 0)
THEN
1701 IF (ir == 1 .AND. is == 1 .AND. it == 1) ipt = 1
1702 IF (ir == 2 .AND. is == 1 .AND. it == 1) ipt = 2
1703 IF (ir == 1 .AND. is == 2 .AND. it == 1) ipt = 3
1704 IF (ir == 1 .AND. is == 1 .AND. it == 2) ipt = 4
1706 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
1708 evar(1,i) = lbuf%SIG(jj(1) + i)
1709 evar(2,i) = lbuf%SIG(jj(2) + i)
1710 evar(3,i) = lbuf%SIG(jj(3) + i)
1711 evar(4,i) = lbuf%SIG(jj(4) + i)
1712 evar(5,i) = lbuf%SIG(jj(5) + i)
1713 evar(6,i) = lbuf%SIG(jj(6) + i)
1717 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1718 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1719 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1720 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1721 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1722 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1730 IF(el2fa(nn2+n) /= 0)
THEN
1732 gama(1)=gbuf%GAMA(jj(1) + i)
1733 gama(2)=gbuf%GAMA(jj(2) + i)
1734 gama(3)=gbuf%GAMA(jj(3) + i)
1735 gama(4)=gbuf%GAMA(jj(4) + i)
1736 gama(5)=gbuf%GAMA(jj(5) + i)
1737 gama(6)=gbuf%GAMA(jj(6) + i)
1746 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1752 ELSEIF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
1753 ipt = mod(abs(pti)/10,10)
1754 IF ( ipt > 0 .AND. ipt<=nlay)
THEN
1756 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
1758 evar(1,i) = lbuf%SIG(jj(1) + i)
1759 evar(2,i) = lbuf%SIG(jj(2) + i)
1760 evar(3,i) = lbuf%SIG(jj(3) + i)
1761 evar(4,i) = lbuf%SIG(jj(4) + i)
1762 evar(5,i) = lbuf%SIG(jj(5) + i)
1763 evar(6,i) = lbuf%SIG(jj(6) + i)
1767 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1768 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1769 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1770 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1771 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1772 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1779 IF(el2fa(nn2+n) /= 0)
THEN
1780 gama(1)= gbuf%GAMA(jj(1) + i)
1786 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1793 IF( nfilsol /= 0 .AND. gbuf%G_FILL /= 0 )
THEN
1795 evar(1,i) = evar(1,i) * gbuf%FILL(i)
1796 evar(2,i) = evar(2,i) * gbuf%FILL(i)
1797 evar(3,i) = evar(3,i) * gbuf%FILL(i)
1798 evar(4,i) = evar(4,i) * gbuf%FILL(i)
1799 evar(5,i) = evar(5,i) * gbuf%FILL(i)
1800 evar(6,i) = evar(6,i) * gbuf%FILL(i)
1805 ELSEIF(itens>=2010.AND.itens<=22109)
THEN
1808 IF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
1809 ipt = mod(abs(pti)/10,201)
1810 IF ( ipt > 0 .AND. ipt<=nlay .AND.nlay>9)
THEN
1812 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
1814 evar(1,i) = lbuf%SIG(jj(1) + i)
1815 evar(2,i) = lbuf%SIG(jj(2) + i)
1816 evar(3,i) = lbuf%SIG(jj(3) + i)
1817 evar(4,i) = lbuf%SIG(jj(4) + i)
1818 evar(5,i) = lbuf%SIG(jj(5) + i)
1819 evar(6,i) = lbuf%SIG(jj(6) + i)
1823 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1824 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1825 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1826 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1827 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1828 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1835 IF(el2fa(nn2+n) /= 0)
THEN
1836 gama(1)= lbuf%GAMA(jj(1) + i)
1837 gama(2)= lbuf%GAMA(jj(2) + i)
1842 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1848 ELSEIF (isolnod == 16.OR.(isolnod == 8 .AND.jhbe == 14))
THEN
1850 icsig = iparg(17,ng)
1852 is0=mod(abs(pti)/10,201)
1853 it0=mod(abs(pti),10)
1854 IF (ir0==0.OR.is0==0.OR.it0==0.OR.nlay<10) cycle
1858 IF (tshell == 1)
THEN
1859 IF (icsig==100)
THEN
1863 ELSEIF (icsig==10)
THEN
1873 IF (ir>nptr.OR.is>npts.OR.it>nlay) cycle
1874 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1875 IF ( ipt <= npt )
THEN
1876 IF (tshell == 1)
THEN
1877 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
1879 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
1882 evar(1,i) = lbuf%SIG(jj(1) + i)
1883 evar(2,i) = lbuf%SIG(jj(2) + i)
1884 evar(3,i) = lbuf%SIG(jj(3) + i)
1885 evar(4,i) = lbuf%SIG(jj(4) + i)
1886 evar(5,i) = lbuf%SIG(jj(5) + i)
1887 evar(6,i) = lbuf%SIG(jj(6) + i)
1891 evar(1,i) =evar(1,i)+lbuf%VISC(jj(1) + i)
1892 evar(2,i) =evar(2,i)+lbuf%VISC(jj(2) + i)
1893 evar(3,i) =evar(3,i)+lbuf%VISC(jj(3) + i)
1894 evar(4,i) =evar(4,i)+lbuf%VISC(jj(4) + i)
1895 evar(5,i) =evar(5,i)+lbuf%VISC(jj(5) + i)
1896 evar(6,i) =evar(6,i)+lbuf%VISC(jj(6) + i)
1900 IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
1907 IF(el2fa(nn2+n) /= 0)
THEN
1910 gama(2)=lbuf%GAMA(jj(1) + i)
1911 gama(3)=lbuf%GAMA(jj(2) + i)
1923 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1929 IF(el2fa(nn2+n) /= 0)
THEN
1931 gama(1)=lbuf%GAMA(jj(1) + i)
1932 gama(2)=lbuf%GAMA(jj(2) + i)
1945 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1951 IF(el2fa(nn2+n) /= 0)
THEN
1953 gama(1)=lbuf%GAMA(jj(2) + i)
1955 gama(3)=lbuf%GAMA(jj(1) + i)
1967 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
1975 IF( nfilsol /= 0 .AND. gbuf%G_FILL /= 0 )
THEN
1977 evar(1,i) = evar(1,i) * gbuf%FILL(i)
1978 evar(2,i) = evar(2,i) * gbuf%FILL(i)
1979 evar(3,i) = evar(3,i) * gbuf%FILL(i)
1980 evar(4,i) = evar(4,i) * gbuf%FILL(i)
1981 evar(5,i) = evar(5,i) * gbuf%FILL(i)
1982 evar(6,i) = evar(6,i) * gbuf%FILL(i)
1987 ELSEIF (itens>=1010.AND.itens<=2009)
THEN
1990 IF (isolnod == 8.AND.npt == 8 .AND. jhbe /= 14 .AND. jhbe /= 24 .AND. jhbe /= 15 .AND. jhbe /= 17)
THEN
1992 is=mod(abs(pti)/10,10)
1994 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
1996 IF(ir <= nptr .AND. is <= npts .AND. it <= nptt)
THEN
1997 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
2000 evar(1,i) = evar(1,i) + lbuf%STRA(jj
2001 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)
2002 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)
2003 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)
2004 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)
2005 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)
2023 IF(el2fa(nn2+n) /= 0)
THEN
2025 gama(1)=gbuf%GAMA(jj(1) + i)
2026 gama(2)=gbuf%GAMA(jj(2) + i)
2027 gama(3)=gbuf%GAMA(jj(3) + i)
2028 gama(4)=gbuf%GAMA(jj(4) + i)
2029 gama(5)=gbuf%GAMA(jj(5) + i)
2030 gama(6)=gbuf%GAMA(jj(6) + i)
2039 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2044 ELSEIF ((isolnod == 8 .OR. npt == 1 .OR. (isolnod == 4 .AND. isrot == 0)) .AND.
2045 . jhbe /= 14 .AND. jhbe /= 15 .AND. jhbe /= 17)
THEN
2048 is=mod(abs(pti)/10,10)
2050 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
2052 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
2053 IF (mlw>=28.AND.mlw /= 49 .OR. mlw == 24)
THEN
2055 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)
2056 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)
2057 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)
2058 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*half
2059 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*half
2060 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*half
2062 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
2064 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)
2065 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)
2066 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)
2068 ELSEIF (istrain > 0)
THEN
2069 IF (mlw /= 14.AND.mlw /= 24.AND.mlw<28.OR. mlw == 49)
THEN
2071 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)
2072 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)
2073 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)
2074 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*half
2075 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*half
2076 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*half
2086 IF(el2fa(nn2+n) /= 0)
THEN
2088 gama(1)=gbuf%GAMA(jj(1) + i)
2089 gama(2)=gbuf%GAMA(jj(2) + i)
2090 gama(3)=gbuf%GAMA(jj(3) + i)
2091 gama(4)=gbuf%GAMA(jj(4) + i)
2092 gama(5)=gbuf%GAMA(jj(5) + i)
2093 gama(6)=gbuf%GAMA(jj(6) + i)
2102 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2107 ELSEIF (isolnod == 16.OR.isolnod == 20.OR.(isolnod == 8.AND. (jhbe == 14 .OR. jhbe == 17)))
THEN
2109 icsig = iparg(17,ng)
2111 is0=mod(abs(pti)/10,10)
2112 it0=mod(abs(pti),10)
2113 IF (ir0==0.OR.is0==0.OR.it0==0) cycle
2117 IF (tshell == 1)
THEN
2118 IF (icsig==100)
THEN
2122 ELSEIF (icsig==10)
THEN
2132 IF (ir>nptr.OR.is>npts) cycle
2133 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
2135 IF (tshell == 1)
THEN
2136 IF (isolnod == 16.AND. is0 <= nlay)
THEN
2137 lbuf => elbuf_tab(ng)%BUFLY(is0)%LBUF(ir,1,it)
2139 ELSEIF (it <= nlay)
THEN
2140 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
2144 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
2148 IF(mlw>=28.AND.mlw /= 49)
THEN
2150! 3*9*3 points d
'integration (r*s*t)
2151 EVAR(1,I) = LBUF%STRA(JJ(1) + I)
2152 EVAR(2,I) = LBUF%STRA(JJ(2) + I)
2153 EVAR(3,I) = LBUF%STRA(JJ(3) + I)
2154 EVAR(4,I) = LBUF%STRA(JJ(4) + I)*HALF
2155 EVAR(5,I) = LBUF%STRA(JJ(5) + I)*HALF
2156 EVAR(6,I) = LBUF%STRA(JJ(6) + I)*HALF
2158.OR.
ELSEIF(MLW == 12 MLW == 14)THEN
2160! 3*9*3 points d'integration(r*s*t)
2161 evar(1,i) = lbuf%EPE(jj(1) + i)
2162 evar(2,i) = lbuf%EPE(jj(2) + i)
2163 evar(3,i) = lbuf%EPE(jj(3) + i)
2165 ELSEIF(mlw == 24 .OR. mlw == 25)
THEN
2168 evar(1,i) = lbuf%STRA(jj(1) + i)
2169 evar(2,i) = lbuf%STRA(jj(2) + i)
2170 evar(3,i) = lbuf%STRA(jj(3) + i)
2171 evar(4,i) = lbuf%STRA(jj(4) + i) * half
2172 evar(5,i) = lbuf%STRA(jj(5) + i) * half
2173 evar(6,i) = lbuf%STRA(jj(6) + i) * half
2175 ELSEIF (mlw == 25)
THEN
2177 evar(1,i) = lbuf%STRA(jj(1) + i)
2178 evar(2,i) = lbuf%STRA(jj(2) + i)
2179 evar(3,i) = lbuf%STRA(jj(3) + i)
2180 evar(4,i) = lbuf%STRA(jj(4) + i) * half
2181 evar(5,i) = lbuf%STRA(jj(5) + i) * half
2182 evar(6,i) = lbuf%STRA(jj(6) + i) * half
2184 ELSEIF(istrain > 0)
THEN
2185 IF(mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
2188 evar(1,i) = lbuf%STRA(jj(1) + i)
2189 evar(2,i) = lbuf%STRA(jj(2) + i)
2190 evar(3,i) = lbuf%STRA(jj(3) + i)
2191 evar(4,i) = lbuf%STRA(jj(4) + i) * half
2192 evar(5,i) = lbuf%STRA(jj(5) + i) * half
2193 evar(6,i) = lbuf%STRA(jj(6) + i) * half
2198 IF (kcvt >1 .AND. jhbe /= 16)
THEN
2201 IF (jhbe == 14.AND.icsig > 0)
THEN
2202 IF (igtyp == 21)
THEN
2207 IF(el2fa(nn2+n) /= 0)
THEN
2210 gama(2)=gbuf%GAMA(jj(1) + i)
2211 gama(3)=gbuf%GAMA(jj(2) + i)
2223 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2229 IF(el2fa(nn2+n) /= 0)
THEN
2231 gama(1)=gbuf%GAMA(jj(1) + i)
2232 gama(2)=gbuf%GAMA(jj(2) + i)
2245 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2251 IF(el2fa(nn2+n) /= 0)
THEN
2253 gama(1)=gbuf%GAMA(jj(2) + i)
2255 gama(3)=gbuf%GAMA(jj(1) + i)
2267 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2276 IF(el2fa(nn2+n) /= 0)
THEN
2279 gama(2)=lbuf%GAMA(jj(1) + i)
2292 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2298 IF(el2fa(nn2+n) /= 0)
THEN
2300 gama(1)=lbuf%GAMA(jj(1) + i)
2301 gama(2)=lbuf%GAMA(jj(2) + i)
2314 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2320 IF(el2fa(nn2+n) /= 0)
THEN
2322 gama(1)=lbuf%GAMA(jj(2) + i)
2324 gama(3)=lbuf%GAMA(jj(1) + i)
2336 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2344 IF(el2fa(nn2+n) /= 0)
THEN
2346 gama(1)=gbuf%GAMA(jj(1) + i)
2347 gama(2)=gbuf%GAMA(jj(2) + i)
2348 gama(3)=gbuf%GAMA(jj(3) + i)
2349 gama(4)=gbuf%GAMA(jj(4) + i)
2350 gama(5)=gbuf%GAMA(jj(5) + i)
2351 gama(6)=gbuf%GAMA(jj(6) + i)
2360 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2367 ELSEIF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
2368 ipt = mod(abs(pti)/10,10)
2369 IF ( ipt > 0 .AND. ipt<=nlay )
THEN
2370 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
2371 IF(mlw>=28.AND.mlw /= 49)
THEN
2373 evar(1,i) = lbuf%STRA(jj(1) + i)
2374 evar(2,i) = lbuf%STRA(jj(2) + i)
2375 evar(3,i) = lbuf%STRA(jj(3) + i)
2376 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2377 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2378 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2380 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
2382 evar(1,i) = lbuf%EPE(jj(1) + i)
2383 evar(2,i) = lbuf%EPE(jj(2) + i)
2384 evar(3,i) = lbuf%EPE(jj(3) + i)
2388 evar(1,i) = lbuf%STRA(jj(1) + i)
2389 evar(2,i) = lbuf%STRA(jj(2) + i)
2390 evar(3,i) = lbuf%STRA(jj(3) + i)
2391 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2392 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2393 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2396 IF (kcvt /= 0 )
THEN
2400 IF(el2fa(nn2+n) /= 0)
THEN
2402 gama(1)=gbuf%GAMA(jj(1) + i)
2403 gama(2)=gbuf%GAMA(jj(2) + i)
2416 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2422 ELSEIF (isolnod==10 .OR. (isolnod==4 .AND. isrot==1))
THEN
2425 is=mod(abs(pti)/10,10)
2428 IF (ir == 1 .AND. is == 1 .AND. it == 1) ipt = 1
2429 IF (ir == 2 .AND. is == 1 .AND. it == 1) ipt = 2
2430 IF (ir == 1 .AND. is == 2 .AND. it == 1) ipt = 3
2431 IF (ir == 1 .AND. is == 1 .AND. it == 2) ipt = 4
2433 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
2434 IF (mlw>=28.AND.mlw /= 49 .OR. mlw == 24)
THEN
2436 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)
2437 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)
2438 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)
2439 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*half
2440 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*half
2441 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*half
2443 ELSEIF (mlw == 12 .OR. mlw == 14)
THEN
2445 evar(1,i) = evar(1,i) + lbuf%EPE(jj(1) + i)
2446 evar(2,i) = evar(2,i) + lbuf%EPE(jj(2) + i)
2447 evar(3,i) = evar(3,i) + lbuf%EPE(jj(3) + i)
2449 ELSEIF (istrain > 0)
THEN
2450 IF (mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
2452 evar(1,i) = evar(1,i) + lbuf%STRA(jj(1) + i)
2453 evar(2,i) = evar(2,i) + lbuf%STRA(jj(2) + i)
2454 evar(3,i) = evar(3,i) + lbuf%STRA(jj(3) + i)
2455 evar(4,i) = evar(4,i) + lbuf%STRA(jj(4) + i)*half
2456 evar(5,i) = evar(5,i) + lbuf%STRA(jj(5) + i)*half
2457 evar(6,i) = evar(6,i) + lbuf%STRA(jj(6) + i)*half
2466 IF(el2fa(nn2+n) /= 0)
THEN
2468 gama(1)=gbuf%GAMA(jj(1) + i)
2469 gama(2)=gbuf%GAMA(jj(2) + i)
2470 gama(3)=gbuf%GAMA(jj(3) + i)
2471 gama(4)=gbuf%GAMA(jj(4) + i)
2472 gama(5)=gbuf%GAMA(jj(5) + i)
2473 gama(6)=gbuf%GAMA(jj(6) + i)
2482 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2490 ELSEIF (itens>=22110.AND.itens<=42209)
THEN
2493 IF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
2494 ipt = mod(abs(pti)/10,201)
2495 IF ( ipt > 0 .AND. ipt<=nlay.AND.nlay>9)
THEN
2496 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
2497 IF(mlw>=28.AND.mlw /= 49)
THEN
2499 evar(1,i) = lbuf%STRA(jj(1) + i)
2500 evar(2,i) = lbuf%STRA(jj(2) + i)
2501 evar(3,i) = lbuf%STRA(jj(3) + i)
2502 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2503 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2504 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2506 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
2508 evar(1,i) = lbuf%EPE(jj(1) + i)
2509 evar(2,i) = lbuf%EPE(jj(2) + i)
2510 evar(3,i) = lbuf%EPE(jj(3) + i)
2514 evar(1,i) = lbuf%STRA(jj(1) + i)
2515 evar(2,i) = lbuf%STRA(jj(2) + i)
2516 evar(3,i) = lbuf%STRA(jj(3) + i)
2517 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2518 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2519 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2522 IF (kcvt /= 0 )
THEN
2526 IF(el2fa(nn2+n) /= 0)
THEN
2528 gama(1)=lbuf%GAMA(jj(1) + i)
2529 gama(2)=lbuf%GAMA(jj(2) + i)
2542 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2548 ELSEIF (isolnod==16.OR.(isolnod==8.AND.jhbe==14))
THEN
2550 icsig = iparg(17,ng)
2552 is0=mod(abs(pti)/10,201)
2553 it0=mod(abs(pti),10)
2554 IF (ir0==0.OR.is0==0.OR.it0==0.OR.nlay<10) cycle
2558 IF (tshell == 1)
THEN
2559 IF (icsig==100)
THEN
2563 ELSEIF (icsig==10)
THEN
2573 IF (ir>nptr.OR.is>npts.OR.it>nlay) cycle
2574 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
2575 IF (ipt <= npt )
THEN
2576 IF (tshell == 1)
THEN
2577 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
2579 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
2581 IF(mlw>=28.AND.mlw /= 49)
THEN
2583 evar(1,i) = lbuf%STRA(jj(1) + i)
2584 evar(2,i) = lbuf%STRA(jj(2) + i)
2585 evar(3,i) = lbuf%STRA(jj(3) + i)
2586 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2587 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2588 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2590 ELSEIF(mlw == 12 .OR. mlw == 14)
THEN
2592 evar(1,i) = lbuf%EPE(jj(1) + i)
2593 evar(2,i) = lbuf%EPE(jj(2) + i)
2594 evar(3,i) = lbuf%EPE(jj(3) + i)
2596 ELSEIF(mlw == 24 .OR. mlw == 25)
THEN
2598 evar(1,i) = lbuf%STRA(jj(1) + i)
2599 evar(2,i) = lbuf%STRA(jj(2) + i)
2600 evar(3,i) = lbuf%STRA(jj(3) + i)
2601 evar(4,i) = lbuf%STRA(jj(4) + i)*half
2602 evar(5,i) = lbuf%STRA(jj(5) + i)*half
2603 evar(6,i) = lbuf%STRA(jj(6) + i)*half
2605 ELSEIF (mlw == 25)
THEN
2607 evar(1,i) = lbuf%STRA(jj(1) + i)
2608 evar(2,i) = lbuf%STRA(jj(2) + i)
2609 evar(3,i) = lbuf%STRA(jj(3) + i)
2610 evar(4,i) = lbuf%STRA(jj(4) + i) * half
2611 evar(5,i) = lbuf%STRA(jj(5) + i) * half
2612 evar(6,i) = lbuf%STRA(jj(6) + i) * half
2614 ELSEIF(istrain > 0)
THEN
2615 IF(mlw /= 14.AND.mlw /= 24.AND.mlw<28)
THEN
2617 evar(1,i) = lbuf%STRA(jj(1) + i)
2618 evar(2,i) = lbuf%STRA(jj(2) + i)
2619 evar(3,i) = lbuf%STRA(jj(3) + i)
2620 evar(4,i) = lbuf%STRA(jj(4) + i) * half
2621 evar(5,i) = lbuf%STRA(jj(5) + i) * half
2622 evar(6,i) = lbuf%STRA(jj(6) + i) * half
2628 IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
2631 IF (jhbe == 14.AND.icsig > 0)
THEN
2636 IF(el2fa(nn2+n) /= 0)
THEN
2639 gama(2)=lbuf%GAMA(jj(1) + i)
2640 gama(3)=lbuf%GAMA(jj(2) + i)
2652 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2658 IF(el2fa(nn2+n) /= 0)
THEN
2660 gama(1)=lbuf%GAMA(jj(1) + i)
2661 gama(2)=lbuf%GAMA(jj(2) + i)
2674 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
2680 IF(el2fa(nn2+n) /= 0)
THEN
2682 gama(1)=lbuf%GAMA(jj(2) + i)
2684 gama(3)=lbuf%GAMA(jj(1) + i)
2696 CALL srota6(x, ixs(1,n), kcvt
2705 ELSEIF (itens >= 42210 .AND. itens <= 43209)
THEN
2709 IF (isolnod == 16.OR.isolnod == 20.OR.(isolnod == 8.AND. (jhbe == 14 .OR. jhbe == 17)))
THEN
2711 icsig = iparg(17,ng)
2713 is0=mod(abs(pti)/10,10)
2714 it0=mod(abs(pti),10)
2715 ipid = ixs(10,1 + nft)
2716 IF (ir0==0.OR.is0==0.OR.it0==0) cycle
2720 IF (tshell == 1)
THEN
2721 IF (icsig==100)
THEN
2725 ELSEIF (icsig==10)
THEN
2735 IF (ir>nptr.OR.is>npts) cycle
2736 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
2738 IF (tshell == 1)
THEN
2739 IF (isolnod == 16.AND. is0 <= nlay)
THEN
2740 lbuf => elbuf_tab(ng)%BUFLY(is0)%LBUF(ir,1,it)
2742 ELSEIF (it <= nlay)
THEN
2743 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
2747 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ir,is,it)
2755 evar(1,i) = lbuf%PLA(jj(1) + i + nel)
2756 evar(2,i) = lbuf%PLA(jj(2) + i + nel)
2757 evar(3,i) = lbuf%PLA(jj(3) + i + nel)
2758 evar(4,i) = lbuf%PLA(jj(4) + i + nel) * half
2759 evar(5,i) = lbuf%PLA(jj(5) + i + nel) * half
2760 evar(6,i) = lbuf%PLA(jj(6) + i + nel) * half
2764 IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
2767 IF (jhbe == 14 .AND. icsig > 0)
THEN
2768 IF (igtyp == 21)
THEN
2773 IF (el2fa(nn2+n) /= 0)
THEN
2776 gama(2) = gbuf%GAMA(jj(1) + i)
2777 gama(3) = gbuf%GAMA(jj(2) + i)
2789 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2795 IF (el2fa(nn2+n) /= 0)
THEN
2797 gama(1) = gbuf%GAMA(jj(1) + i)
2798 gama(2) = gbuf%GAMA(jj(2) + i)
2811 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2817 IF (el2fa(nn2+n) /= 0)
THEN
2819 gama(1) = gbuf%GAMA(jj(2) + i)
2821 gama(3) = gbuf%GAMA(jj(1) + i)
2833 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2842 IF (el2fa(nn2+n) /= 0)
THEN
2845 gama(2) = lbuf%GAMA(jj(1) + i)
2846 gama(3) = lbuf%GAMA(jj(2) + i)
2858 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2864 IF (el2fa(nn2+n) /= 0)
THEN
2866 gama(1) = lbuf%GAMA(jj(1) + i)
2867 gama(2) = lbuf%GAMA(jj(2) + i)
2880 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2886 IF (el2fa(nn2+n) /= 0)
THEN
2888 gama(1) = lbuf%GAMA(jj(2) + i)
2890 gama(3) = lbuf%GAMA(jj(1) + i)
2902 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2910 IF (el2fa(nn2+n) /= 0)
THEN
2912 gama(1) = lbuf%GAMA(jj(1) + i)
2913 gama(2) = lbuf%GAMA(jj(2) + i)
2914 gama(3) = lbuf%GAMA(jj(3) + i)
2915 gama(4) = lbuf%GAMA(jj(4) + i)
2916 gama(5) = lbuf%GAMA(jj(5) + i)
2917 gama(6) = lbuf%GAMA(jj(6) + i)
2926 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
2933 ELSEIF (isolnod == 10 .OR. (isolnod==4 .AND. isrot==1))
THEN
2936 is = mod(abs(pti)/10,10)
2937 it = mod(abs(pti),10)
2939 IF (ir == 1 .AND. is == 1 .AND. it == 1) ipt = 1
2940 IF (ir == 2 .AND. is == 1 .AND. it == 1) ipt = 2
2941 IF (ir == 1 .AND. is == 2 .AND. it == 1) ipt = 3
2942 IF (ir == 1 .AND. is == 1 .AND. it == 2) ipt = 4
2944 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(ipt,1,1)
2947 evar(1,i) = evar(1,i) + lbuf%PLA(jj(1) + i + nel)
2948 evar(2,i) = evar(2,i) + lbuf%PLA(jj(2) + i + nel)
2949 evar(3,i) = evar(3,i) + lbuf%PLA(jj(3) + i + nel)
2950 evar(4,i) = evar(4,i) + lbuf%PLA(jj(4) + i + nel)*half
2951 evar(5,i) = evar(5,i) + lbuf%PLA(jj(5) + i + nel)*half
2952 evar(6,i) = evar(6,i) + lbuf%PLA(jj(6) + i + nel)*half
2960 IF (el2fa(nn2+n) /= 0)
THEN
2962 gama(1) = gbuf%GAMA(jj(1) + i)
2963 gama(2) = gbuf%GAMA(jj(2) + i)
2964 gama(3) = gbuf%GAMA(jj(3) + i)
2965 gama(4) = gbuf%GAMA(jj(4) + i)
2966 gama(5) = gbuf%GAMA(jj(5) + i)
2967 gama(6) = gbuf%GAMA(jj(6) + i)
2976 CALL srota6(x, ixs(1,n), kcvt, evar(1,i),gama, jhbe, igtyp, isorth)
2981 ELSEIF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
2983 ipt = mod(abs(pti)/10,10)
2984 IF ( ipt > 0 .AND. ipt<=nlay)
THEN
2985 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
2988 evar(1,i) = lbuf%PLA(jj(1) + i + nel)
2989 evar(2,i) = lbuf%PLA(jj(2) + i + nel)
2990 evar(3,i) = lbuf%PLA(jj(3) + i + nel)
2991 evar(4,i) = lbuf%PLA(jj(4) + i + nel)*half
2992 evar(5,i) = lbuf%PLA(jj(5) + i + nel)*half
2993 evar(6,i) = lbuf%PLA(jj(6) + i + nel)*half
2998 IF (el2fa(nn2+n) /= 0)
THEN
3000 gama(1)=gbuf%GAMA(jj(1) + i)
3001 gama(2)=gbuf%GAMA(jj(2) + i)
3014 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3021 ELSEIF (itens >= 43210 .AND. itens <= 63309)
THEN
3025 IF ((isolnod == 6 .OR. isolnod == 8).AND.jhbe == 15)
THEN
3027 ipt = mod(abs(pti)/10,201)
3028 IF ( ipt > 0 .AND. ipt<=nlay .AND. nlay>9)
THEN
3029 lbuf => elbuf_tab(ng)%BUFLY(ipt)%LBUF(1,1,1)
3032 evar(1,i) = lbuf%PLA(jj(1) + i + nel)
3033 evar(2,i) = lbuf%PLA(jj(2) + i + nel)
3034 evar(3,i) = lbuf%PLA(jj(3) + i + nel)
3035 evar(4,i) = lbuf%PLA(jj(4) + i + nel)*half
3036 evar(5,i) = lbuf%PLA(jj(5) + i + nel)*half
3037 evar(6,i) = lbuf%PLA(jj(6) + i + nel)*half
3042 IF (el2fa(nn2+n) /= 0)
THEN
3044 gama(1)=gbuf%GAMA(jj(1) + i)
3045 gama(2)=gbuf%GAMA(jj(2) + i)
3058 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3062 ELSEIF ((isolnod == 16.OR.(isolnod == 8 .AND.jhbe == 14)))
THEN
3064 icsig = iparg(17,ng)
3066 is0=mod(abs(pti)/10,201)
3067 it0=mod(abs(pti),10)
3068 IF (ir0==0.OR.is0==0.OR.it0==0.OR.nlay<10) cycle
3072 IF (tshell == 1)
THEN
3073 IF (icsig==100)
THEN
3077 ELSEIF (icsig==10)
THEN
3087 IF (ir>nptr.OR.is>npts) cycle
3088 ipt = ir + ( (is-1) + (it-1)*npts )*nptr
3089 IF (ipt <= npt )
THEN
3090 IF (isolnod == 16)
THEN
3091 lbuf => elbuf_tab(ng)%BUFLY(is0)%LBUF(ir,1,it)
3093 lbuf => elbuf_tab(ng)%BUFLY(it)%LBUF(ir,is,1)
3097 evar(1,i) = lbuf%PLA(jj(1) + i + nel)
3098 evar(2,i) = lbuf%PLA(jj(2) + i + nel)
3099 evar(3,i) = lbuf%PLA(jj(3) + i + nel)
3100 evar(4,i) = lbuf%PLA(jj(4) + i + nel)*half
3101 evar(5,i) = lbuf%PLA(jj(5) + i + nel)*half
3102 evar(6,i) = lbuf%PLA(jj(6) + i + nel)*half
3106 IF (kcvt /= 0 .AND. jhbe /= 16)
THEN
3109 IF (jhbe == 14 .AND. icsig > 0)
THEN
3114 IF (el2fa(nn2+n) /= 0)
THEN
3117 gama(2) = lbuf%GAMA(jj(1) + i)
3118 gama(3) = lbuf%GAMA(jj(2) + i)
3130 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3136 IF (el2fa(nn2+n) /= 0)
THEN
3138 gama(1) = lbuf%GAMA(jj(1) + i)
3139 gama(2) = lbuf%GAMA(jj(2) + i)
3152 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3158 IF (el2fa(nn2+n) /= 0)
THEN
3160 gama(1) = lbuf%GAMA(jj(2) + i)
3162 gama(3) = lbuf%GAMA(jj(1) + i)
3174 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3181 IF (el2fa(nn2+n) /= 0)
THEN
3183 gama(1) = lbuf%GAMA(jj(1) + i)
3184 gama(2) = lbuf%GAMA(jj(2) + i)
3185 gama(3) = lbuf%GAMA(jj(3) + i)
3186 gama(4) = lbuf%GAMA(jj(4) + i)
3187 gama(5) = lbuf%GAMA(jj(5) + i)
3188 gama(6) = lbuf%GAMA(jj(6) + i)
3197 CALL srota6(x, ixs(1,n), kcvt, evar(1,i), gama, jhbe, igtyp, isorth)
3214 IF (isolnod == 16)
THEN
3218 IF(el2fa(nn2+n) /= 0)
THEN
3219 tens(1,el2fa(nn2+n)) = evar(1,i)
3220 tens(2,el2fa(nn2+n)) = evar(2,i)
3221 tens(3,el2fa(nn2+n)) = evar(3,i)
3222 tens(4,el2fa(nn2+n)) = evar(4,i)
3223 tens(5,el2fa(nn2+n)) = evar(5,i)
3224 tens(6,el2fa(nn2+n)) = evar(6,i)
3225 tens(1,el2fa(nn2+n)+1) = evar(1,i)
3226 tens(2,el2fa(nn2+n)+1) = evar(2,i)
3227 tens(3,el2fa(nn2+n)+1) = evar(3,i)
3228 tens(4,el2fa(nn2+n)+1) = evar(4,i)
3229 tens(5,el2fa(nn2+n)+1) = evar(5,i)
3230 tens(6,el2fa(nn2+n)+1) = evar(6,i)
3231 tens(1,el2fa(nn2+n)+2) = evar(1,i)
3232 tens(2,el2fa(nn2+n)+2) = evar(2,i)
3233 tens(3,el2fa(nn2+n)+2) = evar(3,i)
3234 tens(4,el2fa(nn2+n)+2) = evar(4,i)
3235 tens(5,el2fa(nn2+n)+2) = evar(5,i)
3236 tens(6,el2fa(nn2+n)+2) = evar(6,i)
3237 tens(1,el2fa(nn2+n)+3) = evar(1,i)
3238 tens(2,el2fa(nn2+n)+3) = evar(2,i)
3239 tens(3,el2fa(nn2+n)+3) = evar(3,i)
3240 tens(4,el2fa(nn2+n)+3) = evar(4,i)
3241 tens(5,el2fa(nn2+n)+3) = evar(5,i)
3242 tens(6,el2fa(nn2+n)+3) = evar(6,i)
3248 IF(el2fa(nn2+n) /= 0)
THEN
3249 tens(1,el2fa(nn2+n)) = evar(1,i)
3250 tens(2,el2fa(nn2+n)) = evar(2,i)
3251 tens(3,el2fa(nn2+n)) = evar(3,i)
3252 tens(4,el2fa(nn2+n)) = evar(4,i)
3253 tens(5,el2fa(nn2+n)) = evar(5,i)
3254 tens(6,el2fa(nn2+n)) = evar(6,i)
3260 ELSEIF (isph3d == 1.AND.ity == 51)
THEN
3264 iprt=ipartsp(1 + nft)
3266 gbuf => elbuf_tab(ng)%GBUF
3267 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
3270 IF (itens == 1)
THEN
3275 IF (el2fa(nn3+n) /= 0)
THEN
3276 tens(1,el2fa(nn3+n)) = lbuf%SIG(jj(1) + i)
3277 tens(2,el2fa(nn3+n)) = lbuf%SIG(jj(2) + i)
3278 tens(3,el2fa(nn3+n)) = lbuf%SIG(jj(3) + i)
3279 tens(4,el2fa(nn3+n)) = lbuf%SIG(jj(4) + i)
3280 tens(5,el2fa(nn3+n)) = lbuf%SIG(jj(5) + i)
3281 tens(6,el2fa(nn3+n)) = lbuf%SIG(jj(6) + i)
3287 IF (el2fa(nn3+n) /= 0)
THEN
3288 tens(1,el2fa(nn3+n)) = lbuf%SIG(jj(1)+i) + lbuf%VISC(jj(1)+i)
3289 tens(2,el2fa(nn3+n)) = lbuf%SIG(jj(2)+i) + lbuf%VISC(jj(2)+i)
3290 tens(3,el2fa(nn3+n)) = lbuf%SIG(jj(3)+i) + lbuf%VISC(jj(3)+i)
3291 tens(4,el2fa(nn3+n)) = lbuf%SIG(jj(4)+i) + lbuf%VISC(jj(4)+i)
3292 tens(5,el2fa(nn3+n)) = lbuf%SIG(jj(5)+i) + lbuf%VISC(jj(5)+i)
3293 tens(6,el2fa(nn3+n)) = lbuf%SIG(jj(6)+i) + lbuf%VISC(jj(6)+i)
3300 ELSEIF(itens == 4.AND.mlw == 24 .AND. nint(pm(56,mt1)) == 1)
THEN
3304 IF(el2fa(nn3+n) /= 0)
THEN
3305 tens(1,el2fa(nn3+n)) = lbuf%DGLO(jj(1) + i)
3306 tens(2,el2fa(nn3+n)) = lbuf%DGLO(jj(2) + i)
3307 tens(3,el2fa(nn3+n)) = lbuf%DGLO(jj(3) + i)
3308 tens(4,el2fa(nn3+n)) = lbuf%DGLO(jj(4) + i)
3309 tens(5,el2fa(nn3+n)) = lbuf%DGLO(jj(5) + i)
3310 tens(6,el2fa(nn3+n)) = lbuf%DGLO(jj(6) + i)
3318 IF (el2fa(nn3+n) /= 0)
THEN
3319 tens(1,el2fa(nn3+n)) = zero
3320 tens(2,el2fa(nn3+n)) = zero
3321 tens(3,el2fa(nn3+n)) = zero
3322 tens(4,el2fa(nn3+n)) = zero
3323 tens(5,el2fa(nn3+n)) = zero
3324 tens(6,el2fa(nn3+n)) = zero
3329 ELSEIF (ity==101)
THEN
3344 IF (el2fa(nn4+n) /= 0)
THEN
3346 tens(1,el2fa(nn4+n)+j-1) = evar(1,i)
3347 tens(2,el2fa(nn4+n)+j-1) = evar(1,i)
3348 tens(3,el2fa(nn4+n)+j-1) = evar(1,i)
3349 tens(4,el2fa(nn4+n)+j-1) = evar(1,i)
3350 tens(5,el2fa(nn4+n)+j-1) = evar(1,i)
3351 tens(6,el2fa(nn4+n)+j-1) = evar(1,i)
3373 wa(6*n-5) = tens(1,n)
3374 wa(6*n-4) = tens(2,n)
3375 wa(6*n-3) = tens(3,n)
3376 wa(6*n-2) = tens(4,n)
3377 wa(6*n-1) = tens(5,n)
3378 wa(6*n ) = tens(6,n)
3381 buf = numelsg*6 + numels16g*18+numsphg*6
3906 SUBROUTINE tensgps3(ELBUF_TAB,FUNC1 ,FUNC2 ,IPARG ,GEO ,
3907 . IXS ,IXS10 ,IXS16 ,IXS20 ,IXQ ,
3908 . IXC ,IXTG ,IXT ,IXP ,IXR ,
3920#include "implicit_f.inc"
3924#include "vect01_c.inc"
3925#include "mvsiz_p.inc"
3926#include "com01_c.inc"
3927#include "com04_c.inc"
3928#include "param_c.inc"
3934 . func1(3,*),func2(3,*),geo(npropg,*),x(3,*),
3936 INTEGER IPARG(NPARG,*),
3937 . IXS(NIXS,*),IXQ(NIXQ,*),IXC(NIXC,*),IXTG(NIXTG,*),
3938 . IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),
3939 . IXS10(6,*) ,IXS16(8,*) ,IXS20(12,*) ,ITAGPS(*)
3940 TYPE (ELBUF_STRUCT_),
DIMENSION(NGROUP),
TARGET :: ELBUF_TAB
3946 . off, p, vonm2, vonm, s1, s2, s12, s3,
VALUE,
3947 . a1,b1,b2,b3,yeq,f1,m1,m2,m3,
for,
area(mvsiz),
3948 . a_gauss_r,a_gauss_s,a_gauss_t,n1,
3949 . a_gauss_r1,a_gauss_s1,a_gauss_t1,
3950 . a_gauss_p_r,a_gauss_p_s,a_gauss_p_t,
3952 my_real,
ALLOCATABLE,
DIMENSION(:,:) :: evar
3953 INTEGER I,II, NG, NEL, ISS, ISC,NBGAMA,KCVT,
3955 . ISTRAIN,NN, JTURB,MT, IMID, IALEL,IPID,
3956 . NN1,NF,OFFSET,K,INC,KK, IUS, NUVAR,
3957 . inod, isolnod, iprt, liad, nptr, npts, nptt, ipt,
3958 . is, ir, it, nptg,nc(20,mvsiz),nnod,iexpan,ihbe,mpt,ilay,
3959 . icsig,dir,ivisc,jj(6),mat(mvsiz)
3961 TYPE(G_BUFEL_) ,
POINTER :: GBUF
3962 TYPE(L_BUFEL_) ,
POINTER :: LBUF
3964 . a_gauss(9,9),evar_tmp(6),
alpha,beta,alpha_1,beta_1,
3965 . jr0(mvsiz),js0(mvsiz),jt0(mvsiz),nu(mvsiz),sig_hour(mvsiz,6),
3966 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz), xd5(mvsiz),
3967 . xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
3968 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz), yd5(mvsiz),
3969 . yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
3970 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz), zd5(mvsiz),
3971 . zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),
3972 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
3973 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
3974 . r31(mvsiz),r32(mvsiz),r33(mvsiz),
3975 . rx(mvsiz),ry(mvsiz),rz(mvsiz),sx(mvsiz),sy(mvsiz),sz(mvsiz),
3976 . tx(mvsiz),ty(mvsiz),tz(mvsiz),
3977 . xdl(mvsiz), ydl(mvsiz), zdl(mvsiz),evar_t10(6,10),a_heph(3,8)
3978 INTEGER SOL_NODE(3,8), IPERM1(10),IPERM2(10),NN2,ITSH
3979 DATA IPERM1/0,0,0,0,1,2,3,1,2,3/
3980 DATA IPERM2/0,0,0,0,2,3,1,4,4,4/
3986 2 -.577350269189626,0.577350269189626,0. ,
3989 3 -.774596669241483,0. ,0.774596669241483,
3992 4 -.861136311594053,-.339981043584856,0.339981043584856,
3993 4 0.861136311594053,0. ,0. ,
3995 5 -.906179845938664,-.538469310105683,0. ,
3996 5 0.538469310105683,0.906179845938664,0. ,
3998 6 -.932469514203152,-.661209386466265,-.238619186083197,
3999 6 0.238619186083197,0.661209386466265,0.932469514203152,
4001 7 -.949107912342759,-.741531185599394,-.405845151377397,
4002 7 0. ,0.405845151377397,0.741531185599394,
4003 7 0.949107912342759,0. ,0. ,
4004 8 -.960289856497536,-.796666477413627,-.525532409916329,
4005 8 -.183434642495650,0.183434642495650,0.525532409916329,
4006 8 0.796666477413627,0.960289856497536,0. ,
4007 9 -.968160239507626,-.836031107326636,-.613371432700590,
4008 9 -.324253423403809,0. ,0.324253423403809,
4009 9 0.613371432700590,0.836031107326636,0.968160239507626/
4032 CALL my_alloc(evar,6,numnod)
4043 ivisc = iparg(61,ng)
4044 gbuf => elbuf_tab(ng)%GBUF
4046 2 mlw ,nel ,nft ,iad ,ity ,
4047 3 npt ,jale ,ismstr ,jeul ,jtur ,
4048 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
4049 5 nvaux ,jpor ,kcvt ,jclose ,jplasol ,
4050 6 irep ,iint ,igtyp ,israt ,isrot ,
4051 7 icsen ,isorth ,isorthg ,ifailure,jsms )
4053 IF (iparg(8,ng)==1.OR.mlw==0.OR.mlw==13) cycle
4055 isolnod = iparg(28,ng)
4068 gbuf => elbuf_tab(ng)%GBUF
4069 IF (kcvt==1.AND.isorth/=0) kcvt=2
4073 IF(isolnod == 8)
THEN
4075 nc(j,i) = ixs(j+1,n)
4077 ELSEIF(isolnod == 4)
THEN
4082 ELSEIF(isolnod == 6)
THEN
4089 ELSEIF(isolnod == 10)
THEN
4096 nc(j+4,i) = ixs10(j,nn1)
4098 ELSEIF(isolnod == 16)
THEN
4100 nc(j,i) = ixs(j+1,n)
4102 nn1 = n - (numels8+numels10+numels20)
4104 nc(j+8,i) = ixs16(j,nn1)
4106 ELSEIF(isolnod == 20)
THEN
4108 nc(j,i) = ixs(j+1,n)
4110 nn1 = n - (numels8+numels10)
4112 nc(j+8,i) = ixs20(j,nn1)
4117 nptr = elbuf_tab(ng)%NPTR
4118 npts = elbuf_tab(ng)%NPTS
4119 nptt = elbuf_tab(ng)%NPTT
4120 nlay = elbuf_tab(ng)%NLAY
4121 npt = nptr*npts*nptt
4124 IF (jhbe == 24)
THEN
4125 CALL pre_heph(x,ixs,jr0,js0,jt0,pm,mat,nu,nft,nel)
4127 IF(igtyp == 20 .OR. igtyp ==21 .OR. igtyp == 22)
THEN
4133 IF (isolnod == 8.AND. jhbe<9)
THEN
4139 gama(1:6) = gbuf%GAMA(jj(1:6) + i)
4151 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
4152 evar_tmp(1:6) = gbuf%SIG(jj(1:6) + i)
4154 evar_tmp(1:6) =evar_tmp(1:6)+ lbuf%VISC(jj(1:6) + i)
4156 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
4158 evar(1:6,nc(j,i)) = evar(1:6,nc(j,i)) + evar_tmp(1:6)
4161 ELSEIF(isolnod == 6 .OR. isolnod == 8 .OR. isolnod == 16 .OR. isolnod == 20)
THEN
4164 IF(itsh>0 .AND. jhbe /= 14)
THEN
4169 gama(1) = gbuf%GAMA(jj(1) + i)
4170 gama(2) = gbuf%GAMA(jj(2) + i)
4171 gama(3) = gbuf%GAMA(jj(3) + i)
4172 gama(4) = gbuf%GAMA(jj(4) + i)
4173 gama(5) = gbuf%GAMA(jj(5) + i)
4174 gama(6) = gbuf%GAMA(jj(6) + i)
4186 DO j=1,
min(8,isolnod)
4187 DO k=1,
min(8,isolnod)
4188 IF(sol_node(2,k) == sol_node(2,j))
THEN
4190 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == -1)
4192 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == 1)
4193 . ir =
max(1,nptr-1)
4194 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == 1)
4196 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == -1)
4198 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == -1)
4200 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == 1)
4201 . is =
max(1,npts-1)
4202 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == 1)
4204 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == -1)
4206 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == -1)
4208 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == 1)
4209 . it =
max(1,nptt-1)
4210 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == 1)
4212 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == -1)
4221 ELSEIF (sol_node(1,j) == -1 )
THEN
4222 a_gauss_r = a_gauss(1,nptr)
4223 a_gauss_r1 = a_gauss(2,nptr)
4225 . (-one-half*(a_gauss_r1+a_gauss_r))/
4226 . (half*(a_gauss_r1-a_gauss_r))
4227 ELSEIF(sol_node(1,j) == 1 )
THEN
4228 a_gauss_r = a_gauss(nptr-1,nptr)
4229 a_gauss_r1 = a_gauss(nptr,nptr)
4231 . (one+half*(a_gauss_r1+a_gauss_r))/
4232 . (half*(a_gauss_r1-a_gauss_r))
4237 ELSEIF (sol_node(2,j) == -1 )
THEN
4238 a_gauss_s = a_gauss(1,npts)
4239 a_gauss_s1 = a_gauss(2,npts)
4241 . (-one-half*(a_gauss_s1+a_gauss_s))/
4242 . (half*(a_gauss_s1-a_gauss_s))
4243 ELSEIF(sol_node(2,j) == 1 )
THEN
4244 a_gauss_s = a_gauss(npts-1,npts)
4245 a_gauss_s1 = a_gauss(npts,npts)
4247 . (one+half*(a_gauss_s1+a_gauss_s))/
4248 . (half*(a_gauss_s1-a_gauss_s))
4253 ELSEIF (sol_node(3,j) == -1 )
THEN
4254 a_gauss_t = a_gauss(1,nptt)
4255 a_gauss_t1 = a_gauss(2,nptt)
4257 . (-one-half*(a_gauss_t1+a_gauss_t))/
4258 . (half*(a_gauss_t1-a_gauss_t))
4259 ELSEIF(sol_node(3,j) == 1 )
THEN
4260 a_gauss_t = a_gauss(nptt-1,nptt)
4261 a_gauss_t1 = a_gauss(nptt,nptt)
4263 . (one+half*(a_gauss_t1+a_gauss_t))/
4264 . (half*(a_gauss_t1-a_gauss_t))
4267 IF (jhbe == 15 .OR. jhbe == 16)
THEN
4271 . (one+sol_node(1,k) * a_gauss_p_r) *
4272 . (one+sol_node(3,k) * a_gauss_p_t) )
4275 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
4276 evar_tmp(1) = lbuf%SIG(jj(1) + i)
4277 evar_tmp(2) = lbuf%SIG(jj(2) + i)
4278 evar_tmp(3) = lbuf%SIG(jj(3) + i)
4279 evar_tmp(4) = lbuf%SIG(jj(4) + i)
4280 evar_tmp(5) = lbuf%SIG(jj(5) + i)
4281 evar_tmp(6) = lbuf%SIG(jj(6) + i)
4283 evar_tmp(1) = evar_tmp(1) + lbuf%VISC(jj(1) + i)
4284 evar_tmp(2) = evar_tmp(2) + lbuf%VISC(jj(2) + i)
4285 evar_tmp(3) = evar_tmp(3) + lbuf%VISC(jj(3) + i)
4286 evar_tmp(4) = evar_tmp(4) + lbuf%VISC(jj(4) + i)
4287 evar_tmp(5) = evar_tmp(5) + lbuf%VISC(jj(5) + i)
4288 evar_tmp(6) = evar_tmp(6) + lbuf%VISC(jj(6) + i)
4290 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
4291 evar(1,nc(j,i)) = evar(1,nc(j,i)) + n1 * evar_tmp(1)
4292 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
4293 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
4294 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
4295 evar(5,nc(j,i)) = evar(5,nc(j,i)) + n1 * evar_tmp(5)
4296 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
4301 ELSEIF (jhbe == 24)
THEN
4306 gama(1:6) = gbuf%GAMA(jj(1:6) + i)
4323 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
4325 CALL szsigpara(jr0 ,js0 ,jt0 ,gbuf%HOURG ,gbuf%SIG ,
4326 . sig_hour ,ksi ,eta ,zeta ,nu ,nel , i)
4327 evar_tmp(1:6) = sig_hour(i,1:6)
4329 evar_tmp(1:6) =evar_tmp(1:6)+ lbuf%VISC(jj(1:6) + i)
4331 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
4332 evar(1:6,nc(j,i)) = evar(1:6,nc(j,i)) + evar_tmp(1:6)
4340 gama(1) = gbuf%GAMA(jj(1) + i)
4341 gama(2) = gbuf%GAMA(jj(2) + i)
4342 gama(3) = gbuf%GAMA(jj(3) + i)
4343 gama(4) = gbuf%GAMA(jj(4) + i)
4344 gama(5) = gbuf%GAMA(jj(5) + i)
4345 gama(6) = gbuf%GAMA(jj(6) + i)
4355 IF(itsh>0) nptt = nlay
4356 DO j=1,
min(8,isolnod)
4357 DO k=1,
min(8,isolnod)
4358 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == -1)
4360 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == 1)
4361 . is =
max(1,npts-1)
4362 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == 1)
4364 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == -1)
4366 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == -1)
4368 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == 1)
4369 . it =
max(1,nptt-1)
4370 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == 1)
4372 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == -1)
4374 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == -1)
4376 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == 1)
4377 . ir =
max(1,nptr-1)
4378 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == 1)
4380 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == -1)
4389 ELSEIF (sol_node(1,j) == -1 )
THEN
4390 a_gauss_r = a_gauss(1,nptr)
4391 a_gauss_r1 = a_gauss(2,nptr)
4393 . (-one-half*(a_gauss_r1+a_gauss_r))/
4394 . (half*(a_gauss_r1-a_gauss_r))
4395 ELSEIF(sol_node(1,j) == 1 )
THEN
4396 a_gauss_r = a_gauss(nptr-1,nptr)
4397 a_gauss_r1 = a_gauss(nptr,nptr)
4399 . (one+half*(a_gauss_r1+a_gauss_r))/
4400 . (half*(a_gauss_r1-a_gauss_r))
4405 ELSEIF (sol_node(2,j) == -1 )
THEN
4406 a_gauss_s = a_gauss(1,npts)
4407 a_gauss_s1 = a_gauss(2,npts)
4409 . (-one-half*(a_gauss_s1+a_gauss_s))/
4410 . (half*(a_gauss_s1-a_gauss_s))
4411 ELSEIF(sol_node(2,j) == 1 )
THEN
4412 a_gauss_s = a_gauss(npts-1,npts)
4413 a_gauss_s1 = a_gauss(npts,npts)
4415 . (one+half*(a_gauss_s1+a_gauss_s))/
4416 . (half*(a_gauss_s1-a_gauss_s))
4421 ELSEIF (sol_node(3,j) == -1 )
THEN
4422 a_gauss_t = a_gauss(1,nptt)
4423 a_gauss_t1 = a_gauss(2,nptt)
4425 . (-one-half*(a_gauss_t1+a_gauss_t))/
4426 . (half*(a_gauss_t1-a_gauss_t))
4427 ELSEIF(sol_node(3,j) == 1 )
THEN
4428 a_gauss_t = a_gauss(nptt-1,nptt)
4429 a_gauss_t1 = a_gauss(nptt,nptt)
4431 . (one+half*(a_gauss_t1+a_gauss_t))/
4432 . (half*(a_gauss_t1-a_gauss_t))
4436 . (one+sol_node(1,k) * a_gauss_p_r) *
4437 . (one+sol_node(2,k) * a_gauss_p_s) *
4438 . (one+sol_node(3,k) * a_gauss_p_t) )
4440 IF (igtyp == 20 .OR. igtyp ==21 .OR. igtyp == 22)
THEN
4449 zeta = a_gauss(it,2)
4451 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
4453 evar_tmp(1) = lbuf%SIG(jj(1) + i)
4454 evar_tmp(2) = lbuf%SIG(jj(2) + i)
4455 evar_tmp(3) = lbuf%SIG(jj(3) + i)
4456 evar_tmp(4) = lbuf%SIG(jj(4) + i)
4457 evar_tmp(5) = lbuf%SIG(jj(5) + i)
4458 evar_tmp(6) = lbuf%SIG(jj(6) + i)
4461 evar_tmp(1) =evar_tmp(1)+ lbuf%VISC(jj(1) + i)
4462 evar_tmp(2) =evar_tmp(2)+ lbuf%VISC(jj(2) + i)
4463 evar_tmp(3) =evar_tmp(3)+ lbuf%VISC(jj(3) + i)
4464 evar_tmp(4) =evar_tmp(4)+ lbuf%VISC(jj(4) + i)
4465 evar_tmp(5) =evar_tmp(5)+ lbuf%VISC(jj(5) + i)
4466 evar_tmp(6) =evar_tmp(6)+ lbuf%VISC(jj(6) + i)
4468 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
4469 evar(1,nc(j,i)) = evar(1,nc(j,i)) + n1 * evar_tmp(1)
4470 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
4471 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
4472 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
4473 evar(5,nc(j,i)) = evar(5,nc(j,i)) + n1 * evar_tmp(5)
4474 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
4481 ELSEIF(isolnod == 4 .AND. isrot/=1)
THEN
4487 gama(1) = gbuf%GAMA(jj(1) + i)
4488 gama(2) = gbuf%GAMA(jj(2) + i)
4489 gama(3) = gbuf%GAMA(jj(3) + i)
4490 gama(4) = gbuf%GAMA(jj(4) + i)
4491 gama(5) = gbuf%GAMA(jj(5) + i)
4492 gama(6) = gbuf%GAMA(jj(6) + i)
4504 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
4505 evar_tmp(1) = lbuf%SIG(jj(1) + i)
4506 evar_tmp(2) = lbuf%SIG(jj(2) + i)
4507 evar_tmp(3) = lbuf%SIG(jj(3) + i)
4508 evar_tmp(4) = lbuf%SIG(jj(4) + i)
4509 evar_tmp(5) = lbuf%SIG(jj(5) + i)
4510 evar_tmp(6) = lbuf%SIG(jj(6) + i)
4512 evar_tmp(1) =evar_tmp(1)+ lbuf%VISC(jj(1) + i)
4513 evar_tmp(2) =evar_tmp(2)+ lbuf%VISC(jj(2) + i)
4514 evar_tmp(3) =evar_tmp(3)+ lbuf%VISC(jj(3) + i)
4515 evar_tmp(4) =evar_tmp(4)+ lbuf%VISC(jj(4) + i)
4516 evar_tmp(5) =evar_tmp(5)+ lbuf%VISC(jj(5) + i)
4517 evar_tmp(6) =evar_tmp(6)+ lbuf%VISC(jj(6) + i)
4519 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
4521 evar(1,nc(j,i)) = evar(1,nc(j,i)) + n1 * evar_tmp(1)
4522 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
4523 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
4524 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
4525 evar(5,nc(j,i)) = evar(5,nc(j,i)) + n1 * evar_tmp(5)
4526 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
4529 ELSEIF(isolnod == 10 .OR. (isolnod == 4 .AND. isrot==1))
THEN
4537 gama(1) = gbuf%GAMA(jj(1) + i)
4538 gama(2) = gbuf%GAMA(jj(2) + i)
4539 gama(3) = gbuf%GAMA(jj(3) + i)
4540 gama(4) = gbuf%GAMA(jj(4) + i)
4541 gama(5) = gbuf%GAMA(jj(5) + i)
4542 gama(6) = gbuf%GAMA(jj(6) + i)
4553 evar_t10(1:6,j)=zero
4565 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
4566 evar_t10(1,j) = evar_t10(1,j)+ n1 *lbuf%SIG(jj(1) + i)
4567 evar_t10(2,j) = evar_t10(2,j)+ n1 *lbuf%SIG(jj(2) + i)
4568 evar_t10(3,j) = evar_t10(3,j)+ n1 *lbuf%SIG(jj(3) + i)
4569 evar_t10(4,j) = evar_t10(4,j)+ n1 *lbuf%SIG(jj(4) + i)
4570 evar_t10(5,j) = evar_t10(5,j)+ n1 *lbuf%SIG(jj(5) + i)
4571 evar_t10(6,j) = evar_t10(6,j)+ n1 *lbuf%SIG(jj(6) + i)
4573 evar_t10(1,j) =evar_t10(1,j)+ n1 *lbuf%VISC(jj(1) + i)
4574 evar_t10(2,j) =evar_t10(2,j)+ n1 *lbuf%VISC(jj(2) + i)
4575 evar_t10(3,j) =evar_t10(3,j)+ n1 *lbuf%VISC(jj(3) + i)
4576 evar_t10(4,j) =evar_t10(4,j)+ n1 *lbuf%VISC(jj(4) + i)
4577 evar_t10(5,j) =evar_t10(5,j)+ n1 *lbuf%VISC(jj(5) + i)
4578 evar_t10(6,j) =evar_t10(6,j)+ n1 *lbuf%VISC(jj(6) + i)
4581 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n), kcvt, evar_t10(1,j),gama, jhbe, igtyp, isorth)
4584 evar(1,nc(j,i)) = evar(1,nc(j,i)) + evar_t10(1,j)
4585 evar(2,nc(j,i)) = evar(2,nc(j,i)) + evar_t10(2,j)
4586 evar(3,nc(j,i)) = evar(3,nc(j,i)) + evar_t10(3,j)
4587 evar(4,nc(j,i)) = evar(4,nc(j,i)) + evar_t10(4,j)
4588 evar(5,nc(j,i)) = evar(5,nc(j,i)) + evar_t10(5,j)
4589 evar(6,nc(j,i)) = evar(6,nc(j,i)) + evar_t10(6,j)
4591 IF(isolnod == 10 )
THEN
4595 evar_t10(1:6,j) = half*(evar_t10(1:6,nn1)+evar_t10(1:6,nn2))
4598 evar(1,nc(j,i)) = evar(1,nc(j,i)) + evar_t10(1,j)
4599 evar(2,nc(j,i)) = evar(2,nc(j,i)) + evar_t10(2,j)
4600 evar(3,nc(j,i)) = evar(3,nc(j,i
4601 evar(4,nc(j,i)) = evar(4,nc(j,i)) + evar_t10(4,j)
4602 evar(5,nc(j,i)) = evar(5,nc(j,i)) + evar_t10(5,j)
4603 evar(6,nc(j,i)) = evar(6,nc(j,i)) + evar_t10(6,j)
4613 func1(k,n) = evar(k,n)
4614 func2(k,n) = evar(k+3,n)
4616 itagps(n) = itagps(n)+1
4749 . IXS ,IXS10 ,IXS16 ,IXS20 ,X ,
4750 . ITAGPS ,PM ,TAG_SKIN_ND )
4760#include "implicit_f.inc"
4764#include "vect01_c.inc"
4765#include "mvsiz_p.inc"
4766#include "com01_c.inc"
4767#include "com04_c.inc"
4768#include "param_c.inc"
4774 . func1(3,*),func2(3,*),x(3,*),
4776 INTEGER IPARG(NPARG,*),TAG_SKIN_ND(*),
4777 . IXS(NIXS,*),IXS10(6,*) ,IXS16(8,*) ,IXS20(12,*) ,ITAGPS(*)
4778 TYPE (ELBUF_STRUCT_),
DIMENSION(NGROUP),
TARGET :: ELBUF_TAB
4783 . OFF, P, VONM2, VONM, S1, S2, S12, S3,
VALUE,
4784 . A1,B1,B2,B3,YEQ,F1,M1,M2,M3,FOR,AREA(MVSIZ),
4785 . A_GAUSS_R,A_GAUSS_S,A_GAUSS_T,N1,
4786 . A_GAUSS_R1,A_GAUSS_S1,A_GAUSS_T1,
4787 . A_GAUSS_P_R,A_GAUSS_P_S,A_GAUSS_P_T,
4789 my_real,
ALLOCATABLE,
DIMENSION(:,:) :: evar
4790 INTEGER I,II, NG, NEL, ISS, ISC,NBGAMA,KCVT,
4792 . istrain,nn, jturb,mt, imid, ialel,ipid,
4793 . nn1,nf,offset,k,inc,kk, ius, nuvar,
4794 . inod, isolnod, iprt, liad, nptr, npts, nptt, ipt,
4795 . is, ir, it, nptg,nc(20,mvsiz),nnod,iexpan,ihbe,mpt,ilay,
4796 . icsig,dir,ivisc,jj(6),mat(mvsiz),iskin(mvsiz)
4798 TYPE(G_BUFEL_) ,
POINTER :: GBUF
4799 TYPE(L_BUFEL_) ,
POINTER :: LBUF
4801 . A_GAUSS(9,9),EVAR_TMP(6),ALPHA,BETA,ALPHA_1,BETA_1,
4802 . JR0(MVSIZ),JS0(MVSIZ),JT0(MVSIZ),NU(MVSIZ),SIG_HOUR(MVSIZ,6),
4803 . XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ), XD5(MVSIZ),
4804 . xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
4805 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz), yd5(mvsiz),
4806 . yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
4807 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz), zd5(mvsiz),
4808 . zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),
4809 . r11(mvsiz),r12(mvsiz),r13(mvsiz),
4810 . r21(mvsiz),r22(mvsiz),r23(mvsiz),
4811 . r31(mvsiz),r32(mvsiz),r33(mvsiz),
4812 . rx(mvsiz),ry(mvsiz),rz(mvsiz),sx(mvsiz),sy(mvsiz),sz(mvsiz),
4813 . tx(mvsiz),ty(mvsiz),tz(mvsiz),
4814 . xdl(mvsiz), ydl(mvsiz), zdl(mvsiz),evar_t10(6,10),a_heph(3,8)
4816 . SOL_NODE(3,8), IPERM1(10),IPERM2(10),NN2
4817 DATA IPERM1/0,0,0,0,1,2,3,1,2,3/
4818 DATA IPERM2/0,0,0,0,2,3,1,4,4,4/
4824 2 -.577350269189626,0.577350269189626,0. ,
4827 3 -.774596669241483,0. ,0.774596669241483,
4830 4 -.861136311594053,-.339981043584856,0.339981043584856,
4831 4 0.861136311594053,0. ,0. ,
4833 5 -.906179845938664,-.538469310105683,0. ,
4834 5 0.538469310105683,0.906179845938664,0. ,
4836 6 -.932469514203152,-.661209386466265,-.238619186083197,
4837 6 0.238619186083197,0.661209386466265,0.932469514203152,
4839 7 -.949107912342759,-.741531185599394,-.405845151377397,
4840 7 0. ,0.405845151377397,0.741531185599394,
4841 7 0.949107912342759,0. ,0. ,
4842 8 -.960289856497536,-.796666477413627,-.525532409916329,
4843 8 -.183434642495650,0.183434642495650,0.525532409916329,
4844 8 0.796666477413627,0.960289856497536,0. ,
4845 9 -.968160239507626,-.836031107326636,-.613371432700590,
4846 9 -.324253423403809,0. ,0.324253423403809,
4847 9 0.613371432700590,0.836031107326636,0.968160239507626/
4870 CALL my_alloc(evar,6,numnod)
4880 ivisc = iparg(61,ng)
4881 gbuf => elbuf_tab(ng)%GBUF
4883 2 mlw ,nel ,nft ,iad ,ity ,
4884 3 npt ,jale ,ismstr ,jeul ,jtur ,
4885 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
4886 5 nvaux ,jpor ,kcvt ,jclose ,jplasol ,
4887 6 irep ,iint ,igtyp ,israt ,isrot ,
4888 7 icsen ,isorth ,isorthg ,ifailure,jsms )
4891 isolnod = iparg(28,ng)
4903 IF (ity == 1.AND.(igtyp==14.OR.igtyp==6))
THEN
4904 gbuf => elbuf_tab(ng)%GBUF
4905 IF (kcvt==1.AND.isorth/=0) kcvt=2
4910 IF(isolnod == 8)
THEN
4912 nc(j,i) = ixs(j+1,n)
4915 iskin(i) = iskin(i) + tag_skin_nd(nc(j,i))
4917 ELSEIF(isolnod == 4)
THEN
4923 iskin(i) = iskin(i) + tag_skin_nd(nc(j,i))
4925 ELSEIF(isolnod == 6)
THEN
4932 ELSEIF(isolnod == 10)
THEN
4939 nc(j+4,i) = ixs10(j,nn1)
4942 iskin(i) = iskin(i) + tag_skin_nd(nc(j,i))
4944 ELSEIF(isolnod == 16)
THEN
4946 nc(j,i) = ixs(j+1,n)
4948 nn1 = n - (numels8+numels10+numels20)
4950 nc(j+8,i) = ixs16(j,nn1)
4952 ELSEIF(isolnod == 20)
THEN
4954 nc(j,i) = ixs(j+1,n)
4956 nn1 = n - (numels8+numels10)
4958 nc(j+8,i) = ixs20(j,nn1)
4961 iskin(i) = iskin(i) + tag_skin_nd(nc(j,i))
4966 nptr = elbuf_tab(ng)%NPTR
4967 npts = elbuf_tab(ng)%NPTS
4968 nptt = elbuf_tab(ng)%NPTT
4969 nlay = elbuf_tab(ng)%NLAY
4970 npt = nptr*npts*nptt
4973 IF (jhbe == 24)
THEN
4974 CALL pre_heph(x,ixs,jr0,js0,jt0,pm,mat,nu,nft,nel)
4977 IF(isolnod == 6 .OR. isolnod == 8 .OR.
4978 . isolnod == 16 .OR. isolnod == 20)
THEN
4981 IF(nlay > 1 .AND. jhbe /= 14)
THEN
4983 IF (iskin(i)==0) cycle
4987 gama(1) = gbuf%GAMA(jj(1) + i)
4988 gama(2) = gbuf%GAMA(jj(2) + i)
4989 gama(3) = gbuf%GAMA(jj(3) + i)
4990 gama(4) = gbuf%GAMA(jj(4) + i)
4991 gama(5) = gbuf%GAMA(jj(5) + i)
4992 gama(6) = gbuf%GAMA(jj(6) + i)
5004 DO j=1,
min(8,isolnod)
5005 DO k=1,
min(8,isolnod)
5006 IF(sol_node(2,k) == sol_node(2,j))
THEN
5008 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == -1)
5010 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == 1)
5011 . ir =
max(1,nptr-1)
5012 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == 1)
5014 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == -1)
5016 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == -1)
5018 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == 1)
5019 . is =
max(1,npts-1)
5020 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == 1)
5022 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == -1)
5024 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == -1)
5026 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == 1)
5027 . it =
max(1,nptt-1)
5028 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == 1)
5030 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == -1)
5039 ELSEIF (sol_node(1,j) == -1 )
THEN
5040 a_gauss_r = a_gauss(1,nptr)
5041 a_gauss_r1 = a_gauss(2,nptr)
5043 . (-one-half*(a_gauss_r1+a_gauss_r))/
5044 . (half*(a_gauss_r1-a_gauss_r))
5045 ELSEIF(sol_node(1,j) == 1 )
THEN
5046 a_gauss_r = a_gauss(nptr-1,nptr)
5047 a_gauss_r1 = a_gauss(nptr,nptr)
5049 . (one+half*(a_gauss_r1+a_gauss_r))/
5050 . (half*(a_gauss_r1-a_gauss_r))
5055 ELSEIF (sol_node(2,j) == -1 )
THEN
5056 a_gauss_s = a_gauss(1,npts)
5057 a_gauss_s1 = a_gauss(2,npts)
5059 . (-one-half*(a_gauss_s1+a_gauss_s))/
5060 . (half*(a_gauss_s1-a_gauss_s))
5061 ELSEIF(sol_node(2,j) == 1 )
THEN
5062 a_gauss_s = a_gauss(npts-1,npts)
5063 a_gauss_s1 = a_gauss(npts,npts)
5065 . (one+half*(a_gauss_s1+a_gauss_s))/
5066 . (half*(a_gauss_s1-a_gauss_s))
5071 ELSEIF (sol_node(3,j) == -1 )
THEN
5072 a_gauss_t = a_gauss(1,nptt)
5073 a_gauss_t1 = a_gauss(2,nptt)
5075 . (-one-half*(a_gauss_t1+a_gauss_t))/
5076 . (half*(a_gauss_t1-a_gauss_t))
5077 ELSEIF(sol_node(3,j) == 1 )
THEN
5078 a_gauss_t = a_gauss(nptt-1,nptt)
5079 a_gauss_t1 = a_gauss(nptt,nptt)
5081 . (one+half*(a_gauss_t1+a_gauss_t))/
5082 . (half*(a_gauss_t1-a_gauss_t))
5085 IF (jhbe == 15 .OR. jhbe == 16)
THEN
5089 . (one+sol_node(1,k) * a_gauss_p_r) *
5090 . (one+sol_node(3,k) * a_gauss_p_t) )
5093 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
5094 evar_tmp(1) = lbuf%SIG(jj(1) + i)
5095 evar_tmp(2) = lbuf%SIG(jj(2) + i)
5096 evar_tmp(3) = lbuf%SIG(jj(3) + i)
5097 evar_tmp(4) = lbuf%SIG(jj(4) + i)
5098 evar_tmp(5) = lbuf%SIG(jj(5) + i)
5099 evar_tmp(6) = lbuf%SIG(jj(6) + i)
5101 evar_tmp(1) = evar_tmp(1) + lbuf%VISC(jj(1) + i)
5102 evar_tmp(2) = evar_tmp(2) + lbuf%VISC(jj(2) + i)
5103 evar_tmp(3) = evar_tmp(3) + lbuf%VISC(jj(3) + i)
5104 evar_tmp(4) = evar_tmp(4) + lbuf%VISC(jj(4) + i)
5105 evar_tmp(5) = evar_tmp(5) + lbuf%VISC(jj(5) + i)
5106 evar_tmp(6) = evar_tmp(6) + lbuf%VISC(jj(6) + i)
5108 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
5109 evar(1,nc(j,i)) = evar(1,nc(j,i)) + n1 * evar_tmp(1)
5110 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
5111 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
5112 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
5113 evar(5,nc(j,i)) = evar(5,nc(j,i)) + n1 * evar_tmp(5)
5114 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
5119 ELSEIF (jhbe == 24)
THEN
5121 IF (iskin(i)==0) cycle
5125 gama(1:6) = gbuf%GAMA(jj(1:6) + i)
5142 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
5144 CALL szsigpara(jr0 ,js0 ,jt0 ,gbuf%HOURG ,gbuf%SIG ,
5145 . sig_hour ,ksi ,eta ,zeta ,nu ,nel , i)
5146 evar_tmp(1:6) = sig_hour(i,1:6)
5148 evar_tmp(1:6) =evar_tmp(1:6)+ lbuf%VISC(jj(1:6) + i)
5150 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
5151 evar(1:6,nc(j,i)) = evar(1:6,nc(j,i)) + evar_tmp(1:6)
5156 IF (iskin(i)==0) cycle
5160 gama(1) = gbuf%GAMA(jj(1) + i)
5161 gama(2) = gbuf%GAMA(jj(2) + i)
5162 gama(3) = gbuf%GAMA(jj(3) + i)
5163 gama(4) = gbuf%GAMA(jj(4) + i)
5164 gama(5) = gbuf%GAMA(jj(5) + i)
5165 gama(6) = gbuf%GAMA(jj(6) + i)
5175 IF(igtyp == 20 .OR. igtyp ==21 .OR. igtyp == 22)
THEN
5178 DO j=1,
min(8,isolnod)
5179 DO k=1,
min(8,isolnod)
5180 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == -1)
5182 IF (sol_node(1,k) == -1 .AND. sol_node(1,j) == 1)
5183 . is =
max(1,npts-1)
5184 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == 1)
5186 IF (sol_node(1,k) == 1 .AND. sol_node(1,j) == -1)
5188 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == -1)
5190 IF (sol_node(2,k) == -1 .AND. sol_node(2,j) == 1)
5191 . it =
max(1,nptt-1)
5192 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == 1)
5194 IF (sol_node(2,k) == 1 .AND. sol_node(2,j) == -1)
5196 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == -1)
5198 IF (sol_node(3,k) == -1 .AND. sol_node(3,j) == 1)
5199 . ir =
max(1,nptr-1)
5200 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == 1)
5202 IF (sol_node(3,k) == 1 .AND. sol_node(3,j) == -1)
5211 ELSEIF (sol_node(1,j) == -1 )
THEN
5212 a_gauss_r = a_gauss(1,nptr)
5213 a_gauss_r1 = a_gauss(2,nptr)
5215 . (-one-half*(a_gauss_r1+a_gauss_r))/
5217 ELSEIF(sol_node(1,j) == 1 )
THEN
5218 a_gauss_r = a_gauss(nptr-1,nptr)
5219 a_gauss_r1 = a_gauss(nptr,nptr)
5221 . (one+half*(a_gauss_r1+a_gauss_r))/
5222 . (half*(a_gauss_r1-a_gauss_r))
5227 ELSEIF (sol_node(2,j) == -1 )
THEN
5228 a_gauss_s = a_gauss(1,npts)
5229 a_gauss_s1 = a_gauss(2,npts)
5231 . (-one-half*(a_gauss_s1+a_gauss_s))/
5232 . (half*(a_gauss_s1-a_gauss_s))
5233 ELSEIF(sol_node(2,j) == 1 )
THEN
5234 a_gauss_s = a_gauss(npts-1,npts)
5235 a_gauss_s1 = a_gauss(npts,npts)
5237 . (one+half*(a_gauss_s1+a_gauss_s))/
5238 . (half*(a_gauss_s1-a_gauss_s))
5243 ELSEIF (sol_node(3,j) == -1 )
THEN
5244 a_gauss_t = a_gauss(1,nptt)
5245 a_gauss_t1 = a_gauss(2,nptt)
5247 . (-one-half*(a_gauss_t1+a_gauss_t))/
5248 . (half*(a_gauss_t1-a_gauss_t))
5249 ELSEIF(sol_node(3,j) == 1 )
THEN
5250 a_gauss_t = a_gauss(nptt-1,nptt)
5251 a_gauss_t1 = a_gauss(nptt,nptt)
5253 . (one+half*(a_gauss_t1+a_gauss_t))/
5254 . (half*(a_gauss_t1-a_gauss_t))
5258 . (one+sol_node(1,k) * a_gauss_p_r) *
5259 . (one+sol_node(2,k) * a_gauss_p_s) *
5260 . (one+sol_node(3,k) * a_gauss_p_t) )
5262 IF (igtyp == 20 .OR. igtyp ==21 .OR. igtyp == 22)
THEN
5271 zeta = a_gauss(it,2)
5273 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
5275 evar_tmp(1) = lbuf%SIG(jj(1) + i)
5276 evar_tmp(2) = lbuf%SIG(jj(2) + i)
5277 evar_tmp(3) = lbuf%SIG(jj(3) + i)
5278 evar_tmp(4) = lbuf%SIG(jj(4) + i)
5279 evar_tmp(5) = lbuf%SIG(jj(5) + i)
5280 evar_tmp(6) = lbuf%SIG(jj(6) + i)
5282 evar_tmp(1) =evar_tmp(1)+ lbuf%VISC(jj(1) + i)
5283 evar_tmp(2) =evar_tmp(2)+ lbuf%VISC(jj(2) + i)
5284 evar_tmp(3) =evar_tmp(3)+ lbuf%VISC(jj(3) + i)
5285 evar_tmp(4) =evar_tmp(4)+ lbuf%VISC(jj(4) + i)
5286 evar_tmp(5) =evar_tmp(5)+ lbuf%VISC(jj(5) + i)
5287 evar_tmp(6) =evar_tmp(6)+ lbuf%VISC(jj(6) + i)
5289 IF (kcvt /= 0)
CALL srota6(x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
5290 evar(1,nc(j,i)) = evar(1,nc(j,i)) + n1 * evar_tmp(1)
5291 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
5292 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
5293 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
5295 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
5301 ELSEIF(isolnod == 4 )
THEN
5304 IF (iskin(i)==0) cycle
5308 gama(1) = gbuf%GAMA(jj(1) + i)
5309 gama(2) = gbuf%GAMA(jj(2) + i)
5310 gama(3) = gbuf%GAMA(jj(3) + i)
5311 gama(4) = gbuf%GAMA(jj(4) + i)
5312 gama(5) = gbuf%GAMA(jj(5) + i)
5313 gama(6) = gbuf%GAMA(jj(6) + i)
5325 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(1,1,1)
5326 evar_tmp(1) = lbuf%SIG(jj(1) + i)
5327 evar_tmp(2) = lbuf%SIG(jj(2) + i)
5328 evar_tmp(3) = lbuf%SIG(jj(3) + i)
5329 evar_tmp(4) = lbuf%SIG(jj(4) + i)
5330 evar_tmp(5) = lbuf%SIG(jj(5) + i)
5331 evar_tmp(6) = lbuf%SIG(jj(6) + i)
5333 evar_tmp(1) =evar_tmp(1)+ lbuf%VISC(jj(1) + i)
5334 evar_tmp(2) =evar_tmp(2)+ lbuf%VISC(jj(2) + i)
5335 evar_tmp(3) =evar_tmp(3)+ lbuf%VISC(jj(3) + i)
5336 evar_tmp(4) =evar_tmp(4)+ lbuf%VISC(jj(4) + i)
5337 evar_tmp(5) =evar_tmp(5)+ lbuf%VISC(jj(5) + i)
5338 evar_tmp(6) =evar_tmp(6)+ lbuf%VISC(jj(6) + i)
5340 IF (kcvt /= 0)
CALL srota6( x, ixs(1,n),kcvt, evar_tmp,gama, jhbe, igtyp, isorth)
5342 evar(1,nc(j,i)) = evar(1,nc(j
5343 evar(2,nc(j,i)) = evar(2,nc(j,i)) + n1 * evar_tmp(2)
5344 evar(3,nc(j,i)) = evar(3,nc(j,i)) + n1 * evar_tmp(3)
5345 evar(4,nc(j,i)) = evar(4,nc(j,i)) + n1 * evar_tmp(4)
5346 evar(5,nc(j,i)) = evar
5347 evar(6,nc(j,i)) = evar(6,nc(j,i)) + n1 * evar_tmp(6)
5350 ELSEIF(isolnod == 10)
THEN
5352 alpha_1 = -alpha/(beta-alpha
5353 beta_1 = (one-alpha)/(beta-alpha)
5355 IF (iskin(i)==0) cycle
5359 gama(1) = gbuf%GAMA(jj(1) + i)
5360 gama(2) = gbuf%GAMA(jj(2) + i)
5361 gama(3) = gbuf%GAMA(jj(3) + i)
5362 gama(4) = gbuf%GAMA(jj(4) + i)
5363 gama(5) = gbuf%GAMA(jj(5) + i)
5375 evar_t10(1:6,j)=zero
5387 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(ir,is,it)
5388 evar_t10(1,j) = evar_t10(1,j)+ n1 *lbuf%SIG(jj(1) + i)
5389 evar_t10(2,j) = evar_t10(2,j)+ n1 *lbuf%SIG(jj(2) + i)
5390 evar_t10(3,j) = evar_t10(3,j)+ n1 *lbuf%SIG(jj(3) + i)
5391 evar_t10(4,j) = evar_t10(4,j)+ n1 *lbuf%SIG
5392 evar_t10(5,j) = evar_t10(5,j)+ n1 *lbuf%SIG(jj(5) + i)
5393 evar_t10(6,j) = evar_t10(6,j)+ n1 *lbuf%SIG(jj(6) + i)
5395 evar_t10(1,j) =evar_t10(1,j)+ n1 *lbuf%VISC(jj(1) + i)
5396 evar_t10(2,j) =evar_t10(2,j)+ n1 *lbuf%VISC(jj(2) + i)
5397 evar_t10(3,j) =evar_t10(3,j)+ n1 *lbuf%VISC(jj(3) + i)
5398 evar_t10(4,j) =evar_t10(4,j)+ n1 *lbuf%VISC(jj(4) + i)
5399 evar_t10(5,j) =evar_t10(5,j)+ n1 *lbuf%VISC(jj(5) + i)
5400 evar_t10(6,j) =evar_t10(6,j)+ n1 *lbuf%VISC(jj(6) + i)
5403 IF (kcvt /= 0)
CALL srota6( x, ixs(1,n), kcvt, evar_t10(1,j), gama, jhbe, igtyp, isorth)
5408 evar_t10(1:6,j) = half*(evar_t10(1:6,nn1)+evar_t10(1:6,nn2))
5411 evar(1,nc(j,i)) = evar(1,nc(j,i)) + evar_t10(1,j)
5412 evar(2,nc(j,i)) = evar(2,nc(j,i)) + evar_t10(2,j)
5413 evar(3,nc(j,i)) = evar(3,nc(j,i)) + evar_t10(3,j)
5414 evar(4,nc(j,i)) = evar(4,nc(j,i)) + evar_t10(4,j)
5415 evar(5,nc(j,i)) = evar(5,nc(j,i)) + evar_t10(5,j)
5416 evar(6,nc(j,i)) = evar(6,nc(j,i)) + evar_t10(6,j)
5421 IF (iskin(i)==0) cycle
5426 func1(k,n) = evar(k,n)
5427 func2(k,n) = evar(k+3,n)
5429 itagps(n) = itagps(n)+1