c23456789012345678901234567890123456789012345678901234567890123456789012


      REAL*8 FUNCTION therm_U(i,j,k) 

c     ******************************************************************
c
c     This routine returns the vertical shear in transport component U 
c
c     ******************************************************************

c     -------------
c     Declarations
c     -------------

c     ''''''''''''''''''''''
      implicit none
      include 'commonvar.f'      
c     ''''''''''''''''''''''

c---- Argument variables

      integer	i,j,k




c     ----------------------
c     Transport component U
c     ----------------------

      therm_U = -sv*g*dy(j) / (2*omega*rho_0)

      therm_U =
     >therm_U * (r(i,j+1,k)-
     >	         r(i,j  ,k)) / (cos(lat(j+1))-
     >				cos(lat(j  )))

      therm_U =
     >therm_U * (dz(k-1)+dz(k))/2.




c     --------------------
c     End of the function
c     --------------------

      return

      end

c     ******************************************************************
