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.