?? stopp.f
字號:
*$declare*$debug SUBROUTINE STOPP(NT) integer NT integer i,j external RANGE,DEDX real RANGE,DEDX**************** real sum********************** real z1,a1 real DIADJ(10) COMMON /GEN/DIADJ include 'stpwtbl.i'****** SUBROUTINE FOR CALCULATING ENERGY-RANGE AND ENERGY-DEDX ******************** TABLES FOR THE PRIMARY PARTICLES ************************************** Original code L. Sihver* 1995-10-31 H. Tomura AIADJ=0.0 MAT=TGNAME(NT) M1=TG1(1,NT) RHO=TG1(2,NT) AIADJ=TG1(3,NT)*?????? NELECT=0.0************** SUM=0.0 DO 10 I=1,M1 SZ(I)=TG2(1,I,NT) SA(I)=TG2(2,I,NT) FM(I)=TG2(3,I,NT) DIADJ(I)=TG2(4,I,NT)*************** SUM=SUM+FM(I)*SA(I) 10 CONTINUE************************* Z1=PZ************************** A1=PA CALL MATPAR(M1,SZ,SA,FM) IF(DIADJ(1).GT.0.0) CALL MAT2 ER1(22)=10. ER1(23)=10.4 I=23 DO 408 J=12,50,2 I=I+1 408 ER1(I)=J DO 409 J=55,200,5 I=I+1 409 ER1(I)=J DO 410 J=210,1000,10 I=I+1 410 ER1(I)=J NR=I ER1(1)=.01 ER1(2)=.015 DO 519 I=2,9 519 ER1(I+1)=.01*FLOAT(I) ER1(11)=.1 ER1(12)=.15 J=12 DO 524 I=2,9 J=J+1 524 ER1(J)=0.1*FLOAT(I) J=20 DO 525 I=2,20 J=J+1 525 ER1(J)=0.5*FLOAT(I) ER1(40)=10.4 I=40 DO 504 J=12,50,2 I=I+1 504 ER1(I)=FLOAT(J) DO 505 J=55,200,5 I=I+1 505 ER1(I)=FLOAT(J) DO 506 J=210,1000,10 I=I+1 506 ER1(I)=FLOAT(J) NR=I CALL STNEW DO 30 I=1,170*********** write(8,*) i,er1(i) RR(I)=RANGE(pz,pa,ER1(I))*********** write(8,*) rr(i) DEDX1R(I)=DEDX(pZ,pA,ER1(I))*********** write(8,*) dedx1r(i) 30 CONTINUE RETURN END?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -