c
c ____________________________
c
      subroutine vecpro(k,t1)
c  based on LApp routine
c     *****************
      real t1(3,3)
c
c *** calculates axis x(k) as vector product of x(k1),x(k2) taken cyclically
c
      k1=k+1
      if(k1.gt.3) k1=k1-3
      k2=k+2
      if(k2.gt.3) k2=k2-3
      t1(k,1)=t1(k1,2)*t1(k2,3)-t1(k1,3)*t1(k2,2)
      t1(k,2)=t1(k1,3)*t1(k2,1)-t1(k1,1)*t1(k2,3)
      t1(k,3)=t1(k1,1)*t1(k2,2)-t1(k1,2)*t1(k2,1)
      return
c  look for vector product in the kth column of t1
      end
c
c -----------------------------------------
c