1!
2!     CalculiX - A 3-dimensional finite element program
3!              Copyright (C) 1998-2021 Guido Dhondt
4!
5!     This program is free software; you can redistribute it and/or
6!     modify it under the terms of the GNU General Public License as
7!     published by the Free Software Foundation(version 2);
8!
9!
10!     This program is distributed in the hope that it will be useful,
11!     but WITHOUT ANY WARRANTY; without even the implied warranty of
12!     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13!     GNU General Public License for more details.
14!
15!     You should have received a copy of the GNU General Public License
16!     along with this program; if not, write to the Free Software
17!     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
18!
19      subroutine printoutnodefluid(prlab,vold,xturb,physcon,ii,node,
20     &  trab,inotr,ntrans,co,mi,xkappa,xmach)
21!
22!     stores results in the .dat file
23!
24      implicit none
25!
26      character*1 local
27      character*6 prlab(*)
28!
29      integer node,ii,j,inotr(2,*),ntrans,mi(*)
30!
31      real*8 trab(7,*),xkappa(*),xmach(*),
32     &  co(3,*),a(3,3),xturb(2,*),physcon(*),vold(0:mi(2),*)
33!
34      local='L'
35!
36      if(prlab(ii)(1:4).eq.'VF  ') then
37         if((ntrans.eq.0).or.(prlab(ii)(6:6).eq.'G')) then
38            write(5,'(i10,1p,3(1x,e13.6))') node,
39     &           (vold(j,node),j=1,3)
40         elseif(inotr(1,node).eq.0) then
41            write(5,'(i10,1p,3(1x,e13.6))') node,
42     &           (vold(j,node),j=1,3)
43         else
44            call transformatrix(trab(1,inotr(1,node)),co(1,node),a)
45            write(5,'(i10,1p,3(1x,e13.6),1x,a1)') node,
46     &      vold(1,node)*a(1,1)+vold(2,node)*a(2,1)+vold(3,node)*a(3,1),
47     &      vold(1,node)*a(1,2)+vold(2,node)*a(2,2)+vold(3,node)*a(3,2),
48     &      vold(1,node)*a(1,3)+vold(2,node)*a(2,3)+vold(3,node)*a(3,3),
49     &      local
50         endif
51      elseif(prlab(ii)(1:4).eq.'PSF ') then
52         write(5,'(i10,1x,1p,e13.6)') node,
53     &           vold(4,node)
54      elseif(prlab(ii)(1:4).eq.'TSF ') then
55         write(5,'(i10,1x,1p,e13.6)') node,
56     &           vold(0,node)
57      elseif(prlab(ii)(1:4).eq.'PTF ') then
58         write(5,'(i10,1x,1p,e13.6)') node,vold(4,node)*
59     &       (1.d0+(xkappa(node)-1.d0)/2*xmach(node)**2)**(xkappa(node)/
60     &       (xkappa(node)-1.d0))
61      elseif(prlab(ii)(1:4).eq.'TTF ') then
62         write(5,'(i10,1x,1p,e13.6)') node,
63     &     vold(0,node)*(1.d0+(xkappa(node)-1.d0)/2*xmach(node)**2)
64      elseif(prlab(ii)(1:4).eq.'CP  ') then
65         write(5,'(i10,1x,1p,e13.6)') node,
66     &            (vold(4,node)-physcon(6))*2.d0/
67     &            (physcon(7)*physcon(5)**2)
68      elseif(prlab(ii)(1:4).eq.'TURB') then
69         write(5,'(i10,1x,1p,e13.6,1p,e13.6)') node,
70     &            xturb(1,node),xturb(2,node)
71      elseif(prlab(ii)(1:4).eq.'MACH') then
72         write(5,'(i10,1x,1p,e13.6)') node,
73     &            xmach(node)
74      endif
75!
76      flush(5)
77!
78      return
79      end
80
81
82
83
84
85
86