You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
READ(1,20) IZE(I),INE(I),EL<br />
20 FORMAT(8X,I3,6X,I3,3X,D11.5)<br />
IF (IZE(I).EQ.0) IZE(I)=IZE(I-1)<br />
ELAS(IZE(I),INE(I))=EL<br />
END IF<br />
25 CONTINUE<br />
C<br />
C Read residual nuclei (with elastic collisions)<br />
C from 'resid' (HTAPE3X output).<br />
C Then elastic collisions are subtracted.<br />
C 'all z' points the end of the file.<br />
C TOTH: sum of residual nuclei per source particle<br />
C NUCLH: total number of nuclides<br />
C<br />
30 OPEN(2,FILE='resid',STATUS='OLD',FORM='FORMATTED')<br />
DO 35, I=1,IS<br />
READ(2,15) C,D<br />
IF (C.EQ.'all z') GOTO 40<br />
IF (D.EQ.'n =') THEN<br />
BACKSPACE 2<br />
READ(2,20) IZ(I),IN(I),HIS<br />
IF (IZ(I).EQ.0) IZ(I)=IZ(I-1)<br />
IA(I)=IN(I)+IZ(I)<br />
HISOT(IZ(I),IN(I))=HIS-ELAS(IZ(I),IN(I))<br />
TOTH=TOTH+HISOT(IZ(I),IN(I))<br />
IF (HISOT(IZ(I),IN(I)).GT.0) NUCLH=NUCLH+1<br />
ENDIF<br />
35 CONTINUE<br />
C<br />
C Read metastable nuclides from 'resid'<br />
C NMET: Number of metastable nuclides<br />
C which are not only metastable<br />
C<br />
40 READ(2,42) B<br />
42 FORMAT(11X,A6)<br />
IF (B.EQ.'metast') THEN<br />
READ(2,'(/)')<br />
DO 48, I=1,IS<br />
READ(2,45) IZM(I),IAM(I),FRACM<br />
45 FORMAT(2X,I3,2X,I3,27X,D11.5)<br />
IF (IZM(I).EQ.0) GOTO 50<br />
IF (.NOT.FRACM.EQ.1) NMET=NMET+1<br />
INM(I)=IAM(I)-IZM(I)<br />
HISOTM(IZM(I),INM(I))=HISOT(IZM(I),INM(I))*FRACM<br />
19