内容发布更新时间 : 2024/11/8 17:47:55星期一 下面是文章的全部内容请认真阅读。
(1) MCNP中物理量的默认单位
长度：cm 通量：MeV
时间：刹 shake (10-8秒) 能量：MeV
温度：MeV (kT)
原子密度：1024个原子/cm3 质量密度：g/cm3
截面：巴 barns (10-24/cm2) 加热量：MeV/collision
此外， 原子质量按照中子质量为1.0计算，这种单位下阿佛伽德罗常数是0.59703109; 程序运行时间以分钟为单位。
MCNP的源代码是用FORTRAN语言编写的。
(2) 输入INP文件的基本形式
INP文件由一些被空行分隔的由一些被空行分隔的输入块组成，主要的输入块是信息块、标题和栅元块、曲面块和数据块等。输入块又由一些被称为卡的输入行组成。论坛的朋友可以参照对比一下，上一讲中的greatwall程序内容。 INP文件的格式如下：
信息块 (可选项) 空行分隔符 (可选项) 标题卡 栅元卡 ……
空行分隔符 曲面卡 ……
空行分隔符 数据卡 ……
空行分隔符 其它 (可选项)
说明如下：
? 信息块的第一行，必须在它的1～8列填写MESSAGE:，后面跟着用空格分隔的参数项。其后的9～80列和空行分隔符之前各行的1～80列都看作信息块内容。可用A=B参数项更改输出文件名，如OUTP = MYOUT。信息块是可选的。信息块提供给MCNP一个执行信息的方法。
? 在信息块之后的第一行是问题的标题卡，它仅限于一行，占用1～80列，可以是任何信息，将作为OUTP文件中各个输出表的标题被复制。
? 用户在栅元块和曲面块中描述问题的几何。栅元由栅元卡描述。空间必须由彼此相邻的栅元填满，栅元之间不能重叠，也不能出现无栅元的空区，否则会出现错误。构建栅元的曲面由曲面卡定义，曲面卡在曲面块中给出。曲面卡和栅元卡的填写方法，将在以后的讲座中予以介绍。
? 曲面块之后是数据块，在数据块中用户描述源、记数方式、材料等。数据卡在以后的讲座中予以详细介绍。
? 数据卡后不管有无空行分隔符均可以运行，不同之处是，如果数据卡后面有空行分隔符，则MCNP将不再读后面的附加行(如果附加行存在的话)。这对某些用户喜欢保留附加的内容是很方便的，如改变这个问题的说明或者与这个输入文件本身有关的信息，那么，这个分隔符就能防止读入这些附加信息。用户可以把希望保留的一些附加信息写在数据块的空行之后，MCNP会将它们复制到OUTP文件末尾。
(3)卡片格式
INP文件的每一行(称之为一张卡片)都限于使用1～80列。大多数卡片都可以按行填写，然而数据卡也允许按列填写。
a) 行输入格式
通常卡片的1～5列用于填写栅元号、曲面号或数据卡的助记名，6～72列填写卡片参数，73～80列为注释，$符号之后也为注释。序号或卡片助记名可以写在1～5列的任何地方。带有粒子标识符助记名可能需要5列以上，但冒号必须写在6列以内。如果1～5列空白，表示本行为上一行的接续行。6列之后可以写数据项，多个数据项之间用空格分隔。一个数据项必须在一行上写完，不得接续到下一行。相同编号的卡片只能有一张。
需要注意的地方： 每行最多80个字符?
不含控制字符，比如：Tab键，你是看不见的。它导致程序出错。? 注释行：?
—标题卡之后的任何位置都可插入?C?; —第一列是字母“C”，且随后四个空格; —从输入数据之后的$符号后开始的内容 (以上三种情况可以单独或同时存在)
? $符号为它所在的那行数据的结束符，在$符号后面的内容作为注释，它可从$符号后面的任一列开始。
填写卡片参数时可以使用以下输入简写功能：
1>. nR功能，表示它前面的数据再重复n遍。例如：2 4R => 2 2 2 2 2 2>. nI功能，表示与其前后相临两个数之间插进n个线性插值点。 例如：1 5I 7 => 1 2 3 4 5 6 7
3>. xM功能，表示数值等于它前面数据的x倍。
例如：5 4M => 4 20
4>. nJ功能，表示从它所在位置跳过n项不指定的数据而使用缺省值。 这四项输入简写功能可以综合运用。
另外，FCn卡是作为注释用的，但它将作为记数类型n的表头文字印出，比如可作为记数的标题。
SCn卡也是作为注释用的，它作为源概率分布n的表头文字印出。具体使用方法请见手册。
b) 列输入格式
列输入只能用于数据块中，对栅元参数和源的描述比较有用。按行输入的栅元重要性、体积、权窗等数据项可读性较差，而且增加或删除栅元时要在行输入卡上仔细寻找相应项。列输入的可读性有很大提高，删除或增加与某一栅元相对应的数据项时也比较方便。
列输入格式的第一行以#开始，#可以放在1～5列的任意位置，卡片助记名逐个放在该行6列以后，在这些助记名之下按列给出数据项。同一个列输入格式块中的卡片必须是同一类卡片，比如都是栅元参数卡、都是曲面参数卡或都是源参数卡等，在#号下面的1～5列放置栅元号、曲面号或源分布号。
c) 粒子标识符
粒子标识符是卡片助记名的一部分，:N表示中子，:P表示光子，有时也能遇到表示中子-光子联合输运的:N P。下面一些数据卡需要粒子标识符：IMP、EXT、FCL、WWN、WWE、WWP、WWCE、DXT、DXC、Fna、PHYS、ESPLT和CUT。
d) 缺省值
MCNP许多输入卡的参数项有缺省值，用户不必每次都给出这些参数，如果卡片输入项有固定顺序，可以使用nJ功能跳过n个输入项。如果卡片上所有数据项都想缺省，只给出卡片助记名即可。有些卡片不给出也有缺省值，如MODE N卡就可以省略。
练习：
请找出上一讲中的greatwall文件里的信息卡MESSAGE:，标题卡Title,栅元卡cell Card(从cell descriptions开始)，曲面卡Surface card(从1 sph 0 0 0 5 $ Sphere at origin radius 2 cm 开始的三行),数据卡Data card，包括源卡Source Card，材料卡Material Cards 。论坛朋友们不必彻底搞懂它们的作用，以后讲座会详细讲解；只需体会书写的格式，看看是否符合上述的规则。