1      subroutine opprint(symops,rotoop,maxops,nops,itype)
2C$Id$
3      implicit none
4      integer maxops, nops, itype
5      character*2 rotoop(maxops)
6      double precision symops(maxops*3,4)
7      integer i, j, k, indx
8c***********************************************************************
9c
10c     print the matrix representations of the operators
11c
12c***********************************************************************
13        write(*,9)
14        write(*,12)
15        write(*,13) nops
16        do 100 i=1,nops
17          indx=(i-1)*3
18c
19c         Add next block to avoid printing out bad rotoop data for pt. grps
20c
21          if(itype.eq.0) then
22            rotoop(i)='pt'
23            write(*,8) rotoop(i)
24          elseif(itype.eq.3) then
25            write(*,8) rotoop(i)
26          endif
27          do 150 j=1,3
28            write(*,10) (symops((indx+j),k), k=1,4)
29  150     continue
30  100   continue
31    8 format(/,25x,a2,' fold Rotoinversion Operator')
32    9 format(//,20x,'---------- FULL LISTING OF GROUP ----------')
33   10 format(19x,4(f10.6))
34   12 format(/,18x,'MATRIX REPRESENTATIONS OF THE GROUP OPERATORS',/)
35   13 format(5x,'EXCLUDING THE IDENTITY OPERATOR THERE ARE',I3,' OPERATO
36     &RS IN THIS GROUP')
37      return
38      end
39