mcnp4b快速入门一个例子 下载本文

451①2③367②8

例3:LB0.inp

#闪烁体无屏模拟单能

1 1 -5.08 -1 -5 6 IMP:P=1 IMP:E=1 $ 某种闪烁晶体,栅元号1,材料号1,该材料质量密度5.08g/cm3 2 2 -0.001293 (-2 -5 7) #1 IMP:P=1 IMP:E=1 $空气 3 3 -2.70 (-3 -4 8) #1 #2 IMP:P=1 IMP:E=1 $铝 4 4 -0.001293 -9 #1 #2 #3 IMP:P=1 IMP:E=1 $空气 5 0 9 IMP:P=0 IMP:E=0 $真空 1 CZ 1.27 $曲面1定义,单位是cm

2 CZ 1.47 $曲面2平行于Z轴,铝壳内径 3 CZ 1.52 $曲面3,铝壳外径

4

PZ 2.66 $曲面4垂直于Z轴,铝壳外径 5 PZ 2.61 6 PZ 0.07 7 PZ 0.05 8 PZ 0 9 SO 150

MODE P E $源定义

SDEF ERG=0.662 POS=0 0 -1 PAR=2

M1 57139 -0.3484 35000 -0.6016 58140 -0.05 $材料定义卡 M2 8016 -0.24377 7014 -0.75623 $空气包括氧和氮, M3 13027 1

M4 8016 -0.24377 7014 -0.75623 $空气包括氧和氮 F8:P 1 $记数卡

E8 0 1022I 2.784 $记数能量卡

1

FT8 GEB 0 0.0252 $按照FWHM= 0.662*3.1% NPS 800000

PRINT $打印输出结果

错例1:一个错误的排出,下列源代码运行是没有错误提示,但是NPS设置很小,例如为800,但是运行的时间很长还没有结束,正常运行时瞬间即结束。为此,调试了一个下午,最终找出,问题在不该有的一个空行,在源的定义卡中,见下面的F8:P 1的上一行,即两行红色直接,不应该有空行。

程序运行正常结束的提示如下:

2

错例2:程序运行非正常结束一例,原因何在?其中一条正确应该为5 3 -0.001293 -13 #1 #2 #3 #4 IMP:P=1 $外罩内是空气 如果写成5 3 -0.001293 -13 $外罩内是空气,就是错误的,运行会出现如下错误。

错例3:编辑输入文件*.inp时,最好在文本文件中进行,如果在word中编辑输入文件,在英文状态下编辑也可以,但是,在中文状态下,输入空格,在运行时就会出错。

3