c23456789012345678901234567890123456789012345678901234567890123456789012


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

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

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

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

c---- Argument variables

      integer	i,j,k




c     ---------------------
c     Transport component V
c     ---------------------

      therm_V = -sv*g / (2*omega*rho_0)

      therm_V =
     >therm_V * (r(i+1,j,k)-
     >	         r(i  ,j,k)) / sin(lat(j))

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




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

      return

      end

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