IMPLICIT REAL*8 (A-H,O-Z)
CHARACTER*5 SCHNAM
CHARACTER*10 HMNAME
COMPLEX*16 V
CHARACTER*2 NP
DIMENSION KB(3),KKM(3,1000),ND(12)
DIMENSION NDES(12),JTRS(12),IPAS(12)
DIMENSION KO(2000),IND(4,500),V(2000),BB(500),D(500)
C OPEN(3,FILE='generator',iostat=ISO,STATUS='OLD')
C WRITE(6,*) ISO
C IF(ISO.NE.0) STOP
READ(1,*) NUMBER
IF(NUMBER.EQ.0) STOP
CALL TSNTNM(NUMBER,3,NC,SCHNAM,HMNAME)
DO 20 INC=1,NC
NNC=INC
CALL TSPNGE(NUMBER,NNC,3)
IF(NNC.EQ.1) READ(1,*) A,B,C
IF(NNC.EQ.1)READ(1,*) CA,CB,CC
CALL TSLATC(A,B,C,CA,CB,CC)
IF(NNC.EQ.1)READ(1,*) NX,NY,NZ
IF(NNC.EQ.1)READ(1,*) AM
CALL TSKPGN(NX,NY,NZ,KKM,ICC,NK)
DO 10 K=1,NK
KX=KKM(1,K)
KY=KKM(2,K)
KZ=KKM(3,K)
CALL TSNMKP(KX,KY,KZ,ICC,NP)
CALL TSKFBZ(KX,KY,KZ,ICC,INDZB)
KB(1)=KKM(1,K)
KB(2)=KKM(2,K)
KB(3)=KKM(3,K)
WRITE(6,600) NP,KX,KY,KZ,ICC,INDZB
600 FORMAT(//' NAME OF K-POINT ',A2,3I3,'/',I3
& ,' IN OR ON B.Z',I2)
CALL TSIREP(KB,ICC,0)
CALL TSPKDS
CALL TSIRDS
CALL DGTRST(JDUB,NRS,MMG,NSTR,NDES,JTRS,IPAS)
WRITE(6,607) MMG,NSTR,NRS
607 FORMAT(' ORDER OF K-POINT GROUP=',I3,' NUMBER OF STAR=',I3
& /' NUMBER OF REPRESENTATION=',I3)
WRITE(6,608) (I,I=1,NRS)
WRITE(6,604) (NDES(I),I=1,NRS)
WRITE(6,605) (JTRS(I),I=1,NRS)
WRITE(6,606) (IPAS(I),I=1,NRS)
608 FORMAT(' NO ',12I4)
604 FORMAT(' DEGENERACY ',12I4)
605 FORMAT(' HERRING SUM',12I4)
606 FORMAT(' PARTNER NO ',12I4)
CALL TSIRNR(NR,NH,ND)
CALL ZZZY41(KB,ICC,AM)
DO 13 I=1,NR
IR=I
CALL TSPWA(KB,ICC,IR,NDI,KO,V,BB,D,IND,ND1,ND2,NNWW)
CALL DSPWA(KB,ICC,IR,NDI,KO,V,BB,D,IND,ND1,ND2,NNWW)
IF(ND(IR).GT.1) THEN
CALL TSPWAM(KB,ICC,IR,NDI,KO,V,BB,D,IND,ND1,ND2,NNWW)
CALL DSPWA(KB,ICC,IR,NDI,KO,V,BB,D,IND,ND1,ND2,NNWW)
END IF
13 CONTINUE
10 CONTINUE
20 CONTINUE
STOP