Atmospheric example
Local variables to be added
The following local variables need to be defined for the examples in this section:
integer iel
double precision d2s3
double precision zent,xuent,xvent,xkent,xeent,tpent
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
(nbmetd, nbmetm, &
zdmet, tmmet, umet , zent ,
ttcabs, xuent )
(nbmetd, nbmetm, &
zdmet, tmmet, vmet , zent ,
ttcabs, xvent )
(nbmetd, nbmetm, &
zdmet, tmmet, ekmet, zent ,
ttcabs, xkent )
(nbmetd, nbmetm, &
zdmet, tmmet, epmet, zent ,
ttcabs, xeent )
rtp(iel,
ir11) = d2s3*xkent
rtp(iel,
ir22) = d2s3*xkent
rtp(iel,
ir33) = d2s3*xkent
elseif (
iturb.eq.50)
then
elseif (
iturb.eq.60)
then
elseif (
iturb.eq.70)
then
endif
(nbmett, nbmetm, &
ztmet, tmmet, tpmet, zent ,
ttcabs, tpent )
endif
enddo
endif
Finalization
At the end of the subroutine, it is recommended to deallocate the work array lstelt, like in basic example.