copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
***s/r vrtstc
*     

      integer function vrtstc (refip1,datestp,nkfin,iun,verbose) 4,2
      implicit none
*     
      integer datestp,nkfin,iun
      logical verbose
      integer refip1(nkfin,3)
**
*
#include "yomdyn.cdk"
*     
**
      character*12 mesglvl
      integer k,nk_tt,nk_qx
      real dum
*
*----------------------------------------------------------------------
*
      vrtstc=-1
*
      do k=1,nkfin
         refip1(k,1) = -1
         refip1(k,2) = -1
         refip1(k,3) = -1
      end do 
*
      call getlvls (refip1(1,2),nk_tt,nkfin,datestp,"TT",iun,verbose)
      if (nk_tt.lt.1) 
     $call getlvls (refip1(1,2),nk_tt,nkfin,datestp,"VT",iun,verbose)
      if (nk_tt.lt.1) 
     $call getlvls (refip1(1,2),nk_tt,nkfin,datestp,"BUOY",iun,verbose)
      if (nk_tt.lt.1) return
*
      if (gngalsig.eq.1) then
         call getlvls (refip1,nk_qx,nkfin,datestp,"PREG",iun,verbose)
         if (nk_qx.lt.1) return
         refip1(1    ,3) = refip1(1,1)
         refip1(nk_tt,3) = refip1(nk_qx,1)
         do k=2,nk_tt-1
            refip1(k,3) = refip1(k,2)
         end do
      else
         do k=1,nk_tt
            refip1(k,1) = refip1(k,2)
            refip1(k,3) = refip1(k,2)
         end do
      endif
*
      if (nk_tt.lt.1) return
      vrtstc = nk_tt
*
      if (gngalsig.ne.1) then  ! Invert column
         do k=1,nk_tt/2
            dum               = refip1 (k,1)
            refip1(k,1)         = refip1 (nk_tt-k+1,1)
            refip1(nk_tt-k+1,1) = dum
            dum               = refip1 (k,2)
            refip1(k,2)         = refip1 (nk_tt-k+1,2)
            refip1(nk_tt-k+1,2) = dum
            dum               = refip1 (k,3)
            refip1(k,3)         = refip1 (nk_tt-k+1,3)
            refip1(nk_tt-k+1,3) = dum
         end do
      endif
*     
      if (.not.verbose) then
         mesglvl='Pressure'
         if (gngalsig.eq.1) mesglvl='Gal-Chen'
         if (gngalsig.eq.2) mesglvl='Sigma'
         write (6,400) mesglvl
      endif
*
 400  format (' (VRTSTC) Vertical structure of this dataset is: ',a9/)
*--------------------------------------------------------------------
*
      return
      end