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