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
double precision, dimension(:), pointer :: cvar_k, cvar_ep, cvar_phi
double precision, dimension(:), pointer :: cvar_fb, cvar_omg
double precision, dimension(:), pointer :: cvar_r11, cvar_r22, cvar_r33
double precision, dimension(:), pointer :: cvar_r12, cvar_r13, cvar_r23
double precision, dimension(:), pointer :: cvar_yfol, cvar_ng, cvar_h2
double precision, dimension(:), pointer :: cvar_scalt
double precision, dimension(:), pointer :: cvar_fvap, cvar_f7m, cvar_fvp2m
double precision, dimension(:), pointer :: cvar_yco2, cvar_yhcn, cvar_yno
double precision, dimension(:), pointer :: cvar_hox
Before user initialization, work arrays lstelt must be allocated, like in basic example.
Before user initialization, work arrays lstelt must be allocated, like in basic example.
d2s3 = 2.d0/3.d0
xkent = 1.d-10
xeent = 1.d-10
cvar_k(iel) = xkent
cvar_ep(iel) = xeent
enddo
cvar_r11(iel) = d2s3*xkent
cvar_r22(iel) = d2s3*xkent
cvar_r33(iel) = d2s3*xkent
cvar_r12(iel) = 0.d0
cvar_r13(iel) = 0.d0
cvar_r23(iel) = 0.d0
cvar_ep(iel) = xeent
enddo
elseif (
iturb.eq.50)
then
cvar_k(iel) = xkent
cvar_ep(iel) = xeent
cvar_phi(iel) = d2s3
cvar_fb(iel) = 0.d0
enddo
elseif (
iturb.eq.60)
then
cvar_k(iel) = xkent
cvar_omg(iel) = xeent/
cmu/xkent
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
cvar_scalt(iel) = h1init
enddo
cvar_fvap(iel) = 0.d0
if ( ieqco2 .ge. 1 ) then
ioxy = 1
dmas = ( oxyo2(ioxy)*wmo2 +oxyn2(ioxy)*wmn2 &
+oxyh2o(ioxy)*wmh2o+oxyco2(ioxy)*wmco2 )
cvar_yco2(iel) = oxyco2(ioxy)*wmco2/dmas
endif
if ( ieqnox .eq. 1 ) then
cvar_hox(iel) = h1init
endif
enddo
endif
There is no work array in this subroutine, thus nothing to do.