mat2.f

来自「计算能量损失,重离子经过介质后的能量沉积LET,射程」· F 代码 · 共 54 行

*$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