*$debug*$declare        SUBROUTINE MAT2	include	'stpwtbl.i'	integer	j,k	real	zasum	REAL ZDIFI(8),DIFI(8)	REAL INDPOT(8)        real		DIADJ        COMMON	/GEN/	DIADJ(10) *	REAL IADJ,INDPOT(8),NATM(10)*        REAL INCIN*        INTEGER FRGN *        CHARACTER MAT*12,BEAM*12,MATD*12,MATFRD*12,FRNAME*12*        CHARACTER INNAME*8,TGNAME*8,NAMEIN*8,ENY*        COMMON /AZSUM/ASUM,ZSUM	DATA ZDIFI/1.,2.,3.,4.,7.,8.,10.,12./	DATA DIFI/16.5,42.,38.,60.,85.,96.,131.,156./*	INCLUDE 'HE1.I'*	include	'HE10.I'	zasum=0	DO 10 J=1,M1            ZASUM=ZASUM+SZ(J)*FM(J)	    INDPOT(J)=0.	    DO 20 K=1,8 20	        IF(ABS(SZ(J)-ZDIFI(K)).LT..1)INDPOT(J)=DIFI(K)	    IF(DIADJ(J).GT.1.)GOTO 10		IF(SZ(J).le.13.)then			DIADJ(J)=12.*SZ(J)+7.		else			DIADJ(J)=9.76*SZ(J)+58.8/SZ(J)**.19 		endif 10	CONTINUE	aiadj=0	DO 50 J=1,M1 50	    AIADJ = AIADJ + FM(J)*SZ(J)*ALOG(DIADJ(J))        AIADJ = EXP(AIADJ/ZASUM)	RETURN	END