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