Libby-Williams gas 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
character(len=80) :: chaine
integer :: iscal, ivar, ii
double precision :: valmax, valmin
double precision, dimension(:), pointer :: cvar_yfm, cvar_fm, cvar_cyfp2m
double precision, dimension(:), pointer :: cvar_fp2m, cvar_coyfp
double precision, dimension(:), pointer :: cvar_scalt, cvar_scal
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:
do igg = 1, ngazgm
enddo
sommq = sommq +
qimp(izone)
enddo
if (abs(sommq).gt.
epzero)
then
fmelm = sommqf / sommq
tentm = sommqt / sommq
else
endif
coefg(1) = fmelm
coefg(2) = (1.d0-fmelm)
mode = -1
endif
cvar_yfm(iel) = 0.0d0*fmelm
cvar_fm(iel) = 0.d0*fmelm
endif
cvar_scalt(iel) = hinit
endif
enddo
valmax = max(valmax,cvar_scal(iel))
valmin = min(valmin,cvar_scal(iel))
enddo
endif
write(
nfecra,2010)chaine(1:8),valmin,valmax
enddo
9001 format( /,&
' cs_user_initialization: variables initialization by user' ,/,&
/)
2000 format( /,&
/,&
' -----------------------------------------------------------' ,/,&
/,&
' ** INITIALIZATION OF VARIABLES FOR Libby-Williams model' ,/,&
' --------------------------------------------------------' ,/,&
' ONLY ONE PASS' ,/,&
' ---------------------------------' ,/,&
' Variable Valeur min Valeur max' ,/,&
' ---------------------------------' )
2010 format( &
2x, a8, e12.4, e12.4 )
2020 format( &
' ---------------------------------' ,/)
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.