The uspt1d subroutine is used to set the 1D-wall thermal module parameters.
Arguments of uspt1d
integer nvar , nscal , nfpt1d
integer iappel
integer ifpt1d(nfpt1d), nppt1d(nfpt1d), iclt1d(nfpt1d)
integer izft1d(nfabor)
double precision dt(ncelet)
double precision eppt1d(nfpt1d) , rgpt1d(nfpt1d) , tppt1d(nfpt1d)
double precision tept1d(nfpt1d) , hept1d(nfpt1d) , fept1d(nfpt1d)
double precision xlmt1d(nfpt1d) , rcpt1d(nfpt1d) , dtpt1d(nfpt1d)
Local variables declaration
integer ifbt1d , ii , ifac
integer ilelt, nlelt
integer izone
integer, allocatable, dimension(:) :: lstelt
Allocation
Rereading of the restart file
isuit1 = isuite
izone = 0
ifbt1d = 0
iappel = 1 or 2
if (iappel.eq.1.or.iappel.eq.2) then
call getfbr(
'3', nlelt, lstelt)
izone = izone + 1
do ilelt = 1, nlelt
ifac = lstelt(ilelt)
izft1d(ifac) = izone
ifbt1d = ifbt1d + 1
if (iappel.eq.2) ifpt1d(ifbt1d) = ifac
enddo
endif
if (iappel.eq.1) then
nfpt1d = ifbt1d
endif
iappel = 2
if (iappel.eq.2) then
do ii = 1, nfpt1d
ifac = ifpt1d(ii)
nppt1d(ii) = 8
eppt1d(ii) = 0.01144d0
rgpt1d(ii) = 1.d0
tppt1d(ii) = 25.d0 + tkelvi
enddo
endif
iappel = 3
if (iappel.eq.3) then
do ii = 1, nfpt1d
iclt1d(ii) = 1
ifac = ifpt1d(ii)
if (cdgfbo(2,ifac).le.0.025d0) then
iclt1d(ii) = 3
fept1d(ii) = -1.d4
else
iclt1d(ii) = 3
fept1d(ii) = 1.d4
endif
xlmt1d(ii) = 31.5d0
rcpt1d(ii) = 3.5d6
dtpt1d(ii) = 0.3d0
enddo
endif
Deallocation