c----------------------------------------------------------------------
      subroutine TRACE(a,n)
c  produces the trace of an NxN matrix
c
c
c  author: Bill Terilla 6.7.99
c      in:
c     out:
c   calls:
      parameter(q=100)
      integer a(q,q)