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.