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.