materialdata_cp.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	materialdata_cp (imat, ntmat_, t1l, shcon, nshcon, cp)

Function/Subroutine Documentation

materialdata_cp()

subroutine materialdata_cp	(	integer	imat,
		integer	ntmat_,
		real*8	t1l,
		real*8, dimension(0:3,ntmat_,*)	shcon,
		integer, dimension(*)	nshcon,
		real*8	cp
	)

!
      implicit none
!
!     determines the specific heat
!
      integer imat,ntmat_,id,nshcon(*),four
!
      real*8 t1l,shcon(0:3,ntmat_,*),cp
!
      four=4
!     
!     calculating the specific heat
!
      call ident2(shcon(0,1,imat),t1l,nshcon(imat),four,id)
      if(nshcon(imat).eq.0) then
         continue
      elseif(nshcon(imat).eq.1) then
         cp=shcon(1,1,imat)
      elseif(id.eq.0) then
         cp=shcon(1,1,imat)
      elseif(id.eq.nshcon(imat)) then
         cp=shcon(1,id,imat)
      else
         cp=shcon(1,id,imat)+
     &        (shcon(1,id+1,imat)-shcon(1,id,imat))*
     &        (t1l-shcon(0,id,imat))/
     &        (shcon(0,id+1,imat)-shcon(0,id,imat))
      endif
!
      return