calcrhoel.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calcrhoel (nef, vel, rhcon, nrhcon, ielmat, ntmat_, ithermal, mi)

Function/Subroutine Documentation

calcrhoel()

subroutine calcrhoel	(	integer	nef,
		real*8, dimension(nef,0:7)	vel,
		real*8, dimension(0:1,ntmat_,*)	rhcon,
		integer, dimension(*)	nrhcon,
		integer, dimension(mi(3),*)	ielmat,
		integer	ntmat_,
		integer	ithermal,
		integer, dimension(*)	mi
	)

!
!     calculation of rho in the element centers (incompressible
!     fluids)
!
      implicit none
!
      integer nef,i,nrhcon(*),imat,ithermal,ntmat_,mi(*),
     &  ielmat(mi(3),*)
!
      real*8 t1l,vel(nef,0:7),rho,rhcon(0:1,ntmat_,*)
!     
      do i=1,nef
         t1l=vel(i,0)
         imat=ielmat(1,i)
         call materialdata_rho(rhcon,nrhcon,imat,rho,t1l,ntmat_,
     &            ithermal)
         vel(i,5)=rho
c         write(*,*) 'calcrhoel rho',i,rho
      enddo
c      write(*,*)
!            
      return