calcumfa.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calcumfa (nface, vfa, shcon, nshcon, ielmat, ntmat_, ithermal, mi, ielfa, umfa)

Function/Subroutine Documentation

calcumfa()

subroutine calcumfa	(	integer	nface,
		real*8, dimension(0:7,*)	vfa,
		real*8, dimension(0:3,ntmat_,*)	shcon,
		integer, dimension(*)	nshcon,
		integer, dimension(mi(3),*)	ielmat,
		integer	ntmat_,
		integer	ithermal,
		integer, dimension(*)	mi,
		integer, dimension(4,*)	ielfa,
		real*8, dimension(*)	umfa
	)

!
!     calculation of the dynamic viscosity at the face centers
!
      implicit none
!
      integer nface,i,nshcon(*),imat,ithermal,ntmat_,mi(*),
     &  ielmat(mi(3),*),ielfa(4,*)
!
      real*8 t1l,vfa(0:7,*),dvi,shcon(0:3,ntmat_,*),umfa(*)
!     
c$omp parallel default(none)
c$omp& shared(nface,vfa,ielmat,ielfa,shcon,nshcon,ntmat_,ithermal,umfa)
c$omp& private(i,t1l,imat,dvi)
c$omp do
      do i=1,nface
         t1l=vfa(0,i)
!
!        take the material of the first adjacent element
!
         imat=ielmat(1,ielfa(1,i))
         call materialdata_dvi(shcon,nshcon,imat,dvi,t1l,ntmat_,
     &            ithermal)
         umfa(i)=dvi
c         write(*,*) 'calcumfa ',umfa(i)
      enddo
c$omp end do
c$omp end parallel
!            
      return