29 integer ithermal,node,ii,j,inotr(2,*),ntrans,mi(*)
31 real*8 v(0:mi(2),*),t1(*),fn(0:mi(2),*),rftot(0:3),trab(7,*),
32 & co(3,*),a(3,3),veold(0:mi(2),*)
36 if(prlab(ii)(1:4).eq.
'U ')
then 37 if((ntrans.eq.0).or.(prlab(ii)(6:6).eq.
'G'))
then 38 write(5,
'(i10,1p,6(1x,e13.6))') node,
39 & (v(j,node),j=1,mi(2))
40 elseif(inotr(1,node).eq.0)
then 41 write(5,
'(i10,1p,6(1x,e13.6))') node,
42 & (v(j,node),j=1,mi(2))
43 elseif(mi(2).eq.3)
then 45 write(5,
'(i10,1p,3(1x,e13.6),1x,a1)') node,
46 & v(1,node)*a(1,1)+v(2,node)*a(2,1)+v(3,node)*a(3,1),
47 & v(1,node)*a(1,2)+v(2,node)*a(2,2)+v(3,node)*a(3,2),
48 & v(1,node)*a(1,3)+v(2,node)*a(2,3)+v(3,node)*a(3,3),
51 write(*,*)
'*WARNING in printoutnode:' 52 write(*,*)
' for output purposes only 4, 5 or 6' 53 write(*,*)
' degrees of freedom are allowed' 54 write(*,*)
' for generalized vectors;' 55 write(*,*)
' actual degrees of freedom = ',mi(2)
56 write(*,*)
' output request ist not performed;' 58 elseif(prlab(ii)(1:4).eq.
'V ')
then 59 if((ntrans.eq.0).or.(prlab(ii)(6:6).eq.
'G'))
then 60 write(5,
'(i10,1p,3(1x,e13.6))') node,
61 & (veold(j,node),j=1,3)
62 elseif(inotr(1,node).eq.0)
then 63 write(5,
'(i10,1p,3(1x,e13.6))') node,
64 & (veold(j,node),j=1,3)
67 write(5,
'(i10,1p,3(1x,e13.6),1x,a1)') node,
68 & veold(1,node)*a(1,1)+veold(2,node)*a(2,1)+
69 & veold(3,node)*a(3,1),
70 & veold(1,node)*a(1,2)+veold(2,node)*a(2,2)+
71 & veold(3,node)*a(3,2),
72 & veold(1,node)*a(1,3)+veold(2,node)*a(2,3)+
73 & veold(3,node)*a(3,3),
76 elseif((prlab(ii)(1:4).eq.
'NT ').or.
77 & (prlab(ii)(1:4).eq.
'TS '))
then 78 if(ithermal.le.1)
then 79 write(5,
'(i10,1x,1p,e13.6)') node,
82 write(5,
'(i10,1x,1p,e13.6)') node,
85 elseif(prlab(ii)(1:4).eq.
'PS ')
then 86 write(5,
'(i10,1x,1p,e13.6)') node,
88 elseif(prlab(ii)(1:4).eq.
'PN ')
then 89 write(5,
'(i10,1x,1p,e13.6)') node,
91 elseif(prlab(ii)(1:4).eq.
'MF ')
then 92 write(5,
'(i10,1x,1p,e13.6)') node,
94 elseif(prlab(ii)(1:4).eq.
'RF ')
then 96 rftot(j)=rftot(j)+fn(j,node)
98 if(prlab(ii)(5:5).ne.
'O')
then 99 if((ntrans.eq.0).or.(prlab(ii)(6:6).eq.
'G'))
then 100 write(5,
'(i10,1p,3(1x,e13.6))') node,
102 elseif(inotr(1,node).eq.0)
then 103 write(5,
'(i10,1p,3(1x,e13.6))') node,
107 write(5,
'(i10,1p,3(1x,e13.6),1x,a1)') node,
108 & fn(1,node)*a(1,1)+fn(2,node)*a(2,1)+fn(3,node)*a(3,1),
109 & fn(1,node)*a(1,2)+fn(2,node)*a(2,2)+fn(3,node)*a(3,2),
110 & fn(1,node)*a(1,3)+fn(2,node)*a(2,3)+fn(3,node)*a(3,3),
114 elseif(prlab(ii)(1:4).eq.
'RFL ')
then 115 rftot(0)=rftot(0)+fn(0,node)
116 if(prlab(ii)(5:5).ne.
'O')
then 117 write(5,
'(i10,1p,3(1x,e13.6))') node,