pulverized_coal example
Local variables to be added
The following local variables need to be defined for the examples in this section:
integer iel, ige, mode, icla, icha
integer ioxy
double precision t1init, h1init, coefe(
ngazem)
double precision t2init, h2init
double precision xkent, xeent, d2s3
double precision wmh2o,wmco2,wmn2,wmo2,dmas
integer, allocatable, dimension(:) :: lstelt
Allocation
Before user initialization, work arrays lstelt must be allocated, like in basic example.
Initialization
The following initialization block needs to be added for the following examples:
d2s3 = 2.d0/3.d0
xkent = 1.d-10
xeent = 1.d-10
enddo
rtp(iel,
ir11) = d2s3*xkent
rtp(iel,
ir22) = d2s3*xkent
rtp(iel,
ir33) = d2s3*xkent
enddo
elseif (
iturb.eq.50)
then
enddo
elseif (
iturb.eq.60)
then
enddo
elseif (
iturb.eq.70)
then
enddo
endif
do icla = 1, nclacp
icha = ichcor(icla)
enddo
enddo
enddo
coefe(io2) =
wmole(io2)*oxyo2(1) &
coefe(ih2o) =
wmole(ih2o)*oxyh2o(1) &
coefe(ico2) =
wmole(ico2)*oxyco2(1) &
coefe(in2) = 1.d0-coefe(io2)-coefe(ih2o)-coefe(ico2)
enddo
mode = -1
enddo
do icha = 1, ncharb
enddo
enddo
if ( ihtco2 .eq. 1 ) then
endif
if ( noxyd .ge. 2 ) then
endif
if ( noxyd .eq. 3 ) then
endif
if ( ieqco2.ge.1 ) then
ioxy = 1
dmas = ( oxyo2(ioxy)*wmo2 +oxyn2(ioxy)*wmn2 &
+oxyh2o(ioxy)*wmh2o+oxyco2(ioxy)*wmco2 )
xco2 = oxyco2(ioxy)*wmco2/dmas
rtp(iel,
isca(iyco2)) = oxyco2(ioxy)*wmco2/dmas
endif
if ( ieqnox .eq. 1 ) then
rtp(iel,
isca(iyhcn)) = 0.d0
rtp(iel,
isca(iyno )) = 0.d0
rtp(iel,
isca(itaire)) = 293.d0
endif
enddo
endif
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.