Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
GOTO 75<br />
ENDIF<br />
C Same for metastable nuclides<br />
75 IF (HISOTM(IZ(I),IN(I)).EQ.0) GOTO 80<br />
IF (IDFM(IZ(I),IN(I)).EQ.0) THEN<br />
WRITE(6,77) IZ(I),IA(I)<br />
77 FORMAT (8X,I4,1X,I4,'m')<br />
GOTO 80<br />
ELSE<br />
NUCLF=NUCLF+1<br />
MCP2=MCP2+1<br />
TOTF=TOTF+HISOTM(IZ(I),IN(I))<br />
NP2(MCP2)=IDTARG<br />
ND2(MCP2)=IDFM(IZ(I),IN(I))<br />
ZSEKT(MCP2)=HISOTM(IZ(I),IN(I))*PROTSC*CURR/(ATTARG*CPFLUX)<br />
ENDIF<br />
80 CONTINUE<br />
C<br />
C Adds metastable nuclides (NMET) to the number of nuclides (NUCLH).<br />
C Writes out how many nuclides have been expunged by FISPACT.<br />
C<br />
NUCLH=NUCLH+NMET<br />
NUCLF=NUCLH-NUCLF<br />
TOTF=TOTH-TOTF<br />
PERCN=100*(REAL(NUCLF)/REAL(NUCLH))<br />
PERCR=100*(TOTF/TOTH)<br />
WRITE(6,90) NUCLF,NUCLH,PERCN<br />
90 FORMAT(/,1X,'FISPACT2001 has not accepted ',I4,1X,<br />
*'nuclides out of ',I4,' nuclides produced by HTAPE3X (',<br />
*0PF7.3,'%)')<br />
WRITE(6,95) TOTF,TOTH,PERCR<br />
95 FORMAT(1X,'These corre<strong>sp</strong>ond to ',1PE10.4,<br />
*' residual nuclei per source particle out of ',<br />
*1PE10.4,' (',0PF7.3,'%)',/)<br />
c<br />
RETURN<br />
END<br />
21
Change language