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.