C*********************************************************************
C*********************************************************************
C
      SUBROUTINE CLTURB
     O   ( SSTURB,
     I     TKENRG, CDTURB,
     B     TPOT0 , 
     G     ZGRID , FDZ   , FDZM  , DX    ,
     C     GRAV  , 
     W     STABZ , 
     D     NX    , NZ    , NXB   , NZB    )
C
C*******************************************************************
      IMPLICIT REAL*8 ( A-H, O-Z )
C*******************************************************************
C
      DIMENSION
     \  SSTURB(-2:NZB)
C
      DIMENSION
     \  TKENRG(-2:NZB),
     \  CDTURB(-2:NZB)
C
      DIMENSION
     \    ZGRID(-2:NZB), TPOT0(-2:NZB),
     \    FDZM (-2:NZB), FDZ  (-2:NZB)
C
      DIMENSION
     \    STABZ(-2:NZB)
C
C=================================================================
C     Static stability.
C-----------------------------------------------------------
C
      DO 21 IZ = 0, NZ-1
           STABZ(IZ)
     \  =  GRAV / DX / FDZM(IZ) / ( (TPOT0(IZ)+TPOT0(IZ-1))/2.0D0 )
     \   *( (  TPOT0(   IZ   ) - TPOT0(   IZ-1   ) ) )
   21 CONTINUE
C
      CALL BOUND( STABZ, NZ, NZB )
C
C=====================================================================
C
      CDIS = 0.2D0

      DO 31 IZ = 0, NZ

         MIXLNG = MIN( ZGRID(IZ), FDZ(IZ) )
         SSTURB(IZ) =
     \        - (  STABZ(IZ)    + STABZ(IZ+1) ) / 2.0D0
     \        - CDIS * TKENRG(IZ)**1.5D0 / MIXLNG

   31 CONTINUE
C
      CALL BOUND ( SSTURB, NZ, NZB )
C
C---------------------------------------------------------------
C
      RETURN
      E N D
