Fuel 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
integer ioxy
double precision t1init, h1init, coefe(
ngazem)
double precision t2init, h2init
double precision xkent, xeent, d2s3
double precision dmas , wmco2 , wmh2o , wmn2 , wmo2
Allocation
Before user initialization, work arrays lstelt must be allocated, like in basic example.
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
endif
t1init = 1000.d0
t2init = 1000.d0
h2init = h02fol + cp2fol*(t2init-
trefth)
do icla = 1, nclafu
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)
mode = -1
enddo
if ( ieqco2 .ge. 1 ) then
ioxy = 1
dmas = ( oxyo2(ioxy)*wmo2 +oxyn2(ioxy)*wmn2 &
+oxyh2o(ioxy)*wmh2o+oxyco2(ioxy)*wmco2 )
rtp(iel,
isca(iyco2)) = oxyco2(ioxy)*wmco2/dmas
endif
if ( ieqnox .eq. 1 ) then
rtp(iel,
isca(ihox)) = h1init
endif
enddo
endif
Finalization
There is no work array in this subroutine, thus nothing to do.