Gas EBU example
Local variables to be added
The following local variables need to be defined for the examples in this section:
integer iel, mode, igg, izone
double precision hinit, coefg(ngazgm)
double precision sommqf, sommqt, sommq, tentm, fmelm
integer, allocatable, dimension(:) :: lstelt
double precision, dimension(:), pointer :: cvar_ygfm, cvar_fm, cvar_scalt
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:
allocate(lstelt(ncel))
write(nfecra,9001)
do igg = 1, ngazgm
coefg(igg) = zero
enddo
if ( isuite.eq.0 ) then
sommqf = zero
sommq = zero
sommqt = zero
do izone = 1, nozapm
sommqf = sommqf + qimp(izone)*fment(izone)
sommqt = sommqt + qimp(izone)*tkent(izone)
sommq = sommq + qimp(izone)
enddo
if (abs(sommq).gt.epzero) then
fmelm = sommqf / sommq
tentm = sommqt / sommq
else
fmelm = zero
tentm = t0
endif
if ( ippmod(icoebu).eq.1 .or. ippmod(icoebu).eq.3 ) then
coefg(1) = fmelm
coefg(2) = (1.d0-fmelm)
coefg(3) = zero
mode = -1
( mode , ngazg , ngazgm , coefg , &
npo , npot , th , ehgazg , &
hinit , tentm )
endif
do iel = 1, ncel
cvar_ygfm(iel) = 5.d-1
if ( ippmod(icoebu).eq.2 .or. ippmod(icoebu).eq.3 ) then
cvar_fm(iel) = fmelm
endif
if ( ippmod(icoebu).eq.1 .or. ippmod(icoebu).eq.3 ) then
cvar_scalt(iel) = hinit
endif
enddo
endif
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.