extrapolatefluid.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	extrapolatefluid (nk, iponofa, inofa, inum, vfa, v, ielfa, ithermal, imach, ikappa, xmach, xkappa, shcon, nshcon, ntmat_, ielmat, physcon, mi, iturb, xturb)

Function/Subroutine Documentation

extrapolatefluid()

subroutine extrapolatefluid	(	integer	nk,
		integer, dimension(*)	iponofa,
		integer, dimension(2,*)	inofa,
		integer, dimension(*)	inum,
		real*8, dimension(0:7,*)	vfa,
		real*8, dimension(0:4,*)	v,
		integer, dimension(4,*)	ielfa,
		integer	ithermal,
		integer	imach,
		integer	ikappa,
		real*8, dimension(*)	xmach,
		real*8, dimension(*)	xkappa,
		real*8, dimension(0:3,ntmat_,*)	shcon,
		integer, dimension(*)	nshcon,
		integer	ntmat_,
		integer, dimension(mi(3),*)	ielmat,
		real*8, dimension(*)	physcon,
		integer, dimension(*)	mi,
		integer	iturb,
		real*8, dimension(2,*)	xturb
	)

!
!     extrapolates the field values at the center of the faces to
!     the nodes
!
      implicit none
!
      integer nk,iponofa(*),inofa(2,*),inum(*),ielfa(4,*),i,l,indexf,
     &  iface,ithermal,imach,ikappa,imat,nshcon(*),ntmat_,mi(*),
     &  ielmat(mi(3),*),iturb
!
      real*8 vfa(0:7,*),v(0:4,*),cp,r,xk,xmach(*),xkappa(*),t1l,
     &  shcon(0:3,ntmat_,*),physcon(*),xturb(2,*)
!
      do i=1,nk
         if(ithermal.eq.0) then
            do l=1,4
               v(l,i)=0.d0
            enddo
            inum(i)=0
            indexf=iponofa(i)
            do
               if(indexf.eq.0) exit
               iface=inofa(1,indexf)
               do l=1,4
                  v(l,i)=v(l,i)+vfa(l,iface)
               enddo
               inum(i)=inum(i)+1
               indexf=inofa(2,indexf)
            enddo
            if(inum(i).gt.0) then
               do l=1,4
                  v(l,i)=v(l,i)/inum(i)
               enddo
            endif
         else
            do l=0,4
               v(l,i)=0.d0
            enddo
            inum(i)=0
            indexf=iponofa(i)
            do
               if(indexf.eq.0) exit
               iface=inofa(1,indexf)
               do l=0,4
                  v(l,i)=v(l,i)+vfa(l,iface)
               enddo
               if(imach.eq.1) then
                  t1l=vfa(0,iface)
                  imat=ielmat(1,ielfa(1,iface))
                  r=shcon(3,1,imat)
                  call materialdata_cp_sec(imat,ntmat_,t1l,
     &                 shcon,nshcon,cp,physcon)
                  xk=cp/(cp-r)
                  xmach(i)=xmach(i)+dsqrt((vfa(1,iface)**2+
     &               vfa(2,iface)**2+vfa(3,iface)**2)/(xk*r*t1l))
               endif
               if(ikappa.eq.1) then
                  xkappa(i)=xkappa(i)+xk
               endif
!
               inum(i)=inum(i)+1
               indexf=inofa(2,indexf)
            enddo
            if(inum(i).gt.0) then
               do l=0,4
                  v(l,i)=v(l,i)/inum(i)
               enddo
               if(imach.eq.1) xmach(i)=xmach(i)/inum(i)
               if(ikappa.eq.1) xkappa(i)=xkappa(i)/inum(i)
            endif
         endif
!
!        turbulence output
!
         if(iturb.ne.0) then
            do l=1,2
               xturb(l,i)=0.d0
            enddo
            indexf=iponofa(i)
            do
               if(indexf.eq.0) exit
               iface=inofa(1,indexf)
               do l=1,2
                  xturb(l,i)=xturb(l,i)+vfa(l+5,iface)
               enddo
               indexf=inofa(2,indexf)
            enddo
            if(inum(i).gt.0) then
               do l=1,2
                  xturb(l,i)=xturb(l,i)/inum(i)
               enddo
            endif
         endif
      enddo
!  
      return