isca(1) is the number related to the first user-defined scalar variable. rtp(iel,isca(1)) is the value of this variable in cell number iel. ONLY done if there is no restart computation
Unified combustion 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 ( noxyd .ge. 2 ) then
endif
if ( noxyd .ge. 3 ) then
endif
endif
if ( ihtco2 .eq. 1 ) then
endif
if ( ihth2o .eq. 1 ) 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(ihox )) = h1init
endif
enddo
endif
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.