注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

AlexYoung

做好大家都必须要做的事,发展自己感兴趣的事,结束圣人所不齿的事...

 
 
 

日志

 
 

MCNP5 用户手册卷II-附录B MCNP几何绘图与计数绘图  

2011-11-20 21:06:47|  分类: 专业mcnp |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

附录B MCNP几何绘图与计数绘图

译者:杨万奎(Alex

首先我想通过具体的操作命令来理清主要内容,可能更有利于程序的学习。几何绘图部分由于使用得比较多了,基本的命令都还算比较熟悉,故暂时不进行翻译。英文部分共有26页,经过理解简化,但是在不损意思表达的前提之下,可以分为以下几个方面。

0 X-Windows

MCNP5与以前的版本不同之处在于其需要单独的可视化平台,MCNP5 1.14版本支持的X-Windows平台包括:Reflection X (http://www.wrq.com/products/)Hummingbird Exceed_NT (http://www.hummingbird.com/products/nc/exceed/index.html)以及StarnetX-win32 (http://www.starnet.com/). 这三个版本的X-Windows Client 都是商用版本,正版是需要购买授权获得。也支持一个免费的X client——Cygwin支持下的X-Free86.这些Client都不需要开发者或专业版,因为它们所需的X11头文件和库文件已经包含在了MCNP5的发行版里。

另外,通过使用经验得知,通常使用HummingbirdExceed平台,因为这款Client软件已广泛应用于fluent流体计算软件中。由于Windows 7系统的普及,Exceed 13可以支持该系统,之前的版本都不支持Windows 7系统,在使用中应该针对自己的操作系统对版本进行选择。

在绘图之前需要的准备工作包括:首先,安装X-windows Client,之后设置为MCNP5设置好环境变量DISPLAY:  localhost:0,在绘图之前需要打开exceed程序。

1 The Geometry Plotter几何绘图

系统提供的可用命令还是很丰富的,但是一般就我现在的应用来看,用得最多的还是px py pz label center factor origin…但是随着X-windows Client的应用,这些命令都可以在绘图界面中输入。

plot>

?

available commands:

            term              extent            color

            file              px                shade

            viewport          py                ?

            &                pz                help

            end               label             pause

            return             scales            level

            mcplot            center            mbody

            options           factor            interact

            status            theta             mesh

            runtpe            cursor

            basis             restore

            origin            locate

 available colors:

          VioletRed         goldenrod         LightGoldenrod

          blue              plum              LightPink

          yellow            seashell          DodgerBlue

          green             sienna            LightSalmon

          cyan              thistle           LightSkyBlue

          magenta           tomato            LightYellow

          orange            turquoise         MediumOrchid

          pink              wheat             LightSteelBlue

          purple            salmon            MediumPurple

          brown             CadetBlue         OrangeRed

          SlateGray         DarkGoldenrod     PaleGreen

          azure             DarkOliveGreen    PaleTurquoise

          burlywood         SlateBlue         PaleVioletRed

          chartreuse        DarkOrange        LightCyan

          coral             DarkOrchid        RoyalBlue

          cornsilk          DarkSeaGreen      RosyBrown

          firebrick         DarkSlateGray     SkyBlue

          gold              DeepPink          SpringGreen

          honeydew          DeepSkyBlue       SteelBlue

          khaki             AntiqueWhite      red

          maroon            LavenderBlush

          orchid            LightBlue

2 The MCPLOT Tally and Cross-section Plotter计数绘图及截面绘图

mcplot>

?

 available commands:

            term              reset             hist

            file              title             plinear

            &                 below

            coplot            subtitle          bar

            end               xtitle            noerrbar

            return            ytitle            thick

            plot                                thin

            freq              label             legend

            iptal             free              contour

            options           fixed

            printal           set

            printpts          tfc               ?

            status            kcode             help

            runtpe            linlin            pert

            dump              linlog            xs

            wmctal            loglin            mt

            rmctal            loglog            par

            tally             xlims             wash

            nonorm            ylims

            factor            scales

MCPLOT可以用来绘制MCNP的计数结果以及MCNP计算用到的截面数据。可以绘制二维平面图、三维计数等高线图及色彩填充计数三维立体图等。提供诸多绘图选项,可以在同一个二维图中绘制多条曲线。该计数绘图可以来源于MCTAL文件或者是RUNTPE文件。但是如果用MCTAL文件绘图的话,并不能使用所有的绘图选项,因为该文件并未包含所有信息。

MCPLOT也可以绘制输入文件中用到的截面数据:既可以绘制单个核素,也可以绘制整个材料卡中的所有成分。绘制截面数据的目的在于反映MCNP程序对截面的调整,例如:能量截断、中子截面温度、1/v 处理、为提供总的光子截面的光子反映总和、对光子数据的简单物理处理、电子制止能及其它电子数据的生成等。需要注意的是,截面绘图不能从RUNTPE文件获得

该部分内容按以下顺序组织:执行行选项、绘图约定和绘图命令语法格式、根据功能分类的绘图命令、MCTAL文件。此处,MCPLOT绘图选项及关键词是以大写字母出现,但是为方便起见,通常用户可以键入小写字母实现相应功能。

最终的计数结果可以在粒子输运完成即程序执行完毕之后获得,但是也可以在输运进行时显示出一个或多个计数的临时状态(注意:在并行时,该功能不一定可以实现)。在粒子输运完成后,通过MCNP可执行行的Z选项可以调用MCPLOT,既可以作为单独的程序使用RUNTPE文件或MCTAL文件,也可以作为正式的不间断的MCNP运行的一部分进行。在程序运行过程中,有两种方法间歇地请求产生计数绘图:1在输入文件中使用MPLOT卡;2使用TTY中断功能(注意:在并行时,该功能不一定可以实现)。MPLOT输入卡的说明可以参考第三章第144页;TTY中断,<ctrl+c> m,使MCNP在当前历史输运结束时中断输运,从终端发出请求调用MCPLOT

实时计数绘图的限制:在实时绘图过程中,并不会向COMOUT文件写入输出。另外,以下命令无法使用:RMCTALRUNTPEDUMPENDRETURN命令用于推出MCPLOT和返回MCNP的输运模式。使用TTY中断或MPLOT输入卡之后,截面数据将不能显示

MCPLOT可以使用MCTAL文件在另外的计算机上获取计数绘图。在运行结束时要产生MCTAL文件,在输入文件中需要使用PRDMP文件,且其第三个参数非零。MCTAL文件包含了最后一个RUNTPE dump的所有计数数据信息,并且以ASCII码的形式记录,因此其可以转换或转移到另外的机器上。当MCTAL文件被用户定义为生成创建时,可以在运行程序时用以下命令行对该文件进行重命名:

mcnp  i=inpfile  mctal=filename

默认的命名是基于MCTAL的唯一标识名称。

A  Input for MCPLOT and Execution Line Options调用MCPLOT的执行行选项

绘制截面数据最常用的方式:mcnp inp=filename ixz options,在改命令模式下,问题所需的截面数据被读入,但并不进行粒子输运。以下命令可以在改模式下使用:BAR, CONTOUR, DUMP, FREQ, HIST, PLOT, RETURN, RMCTAL, RUNTPE, WASHWMCTAL

下列选项可以在执行行中输入:

NOTEK       禁止在终端产生绘图,而是发送所有的绘图到图形文件PLOTM.ps中。该命令用于批量生产或是终端没有图形处理能力的情形;

COM=aaaa 用文件aaaa作为绘图请求的来源。当文件末(EOF)被读入时,控制转向终端。在批量生产的情形下,在文件末写入END命令来禁止控制向终端转移。绝对不要让COM文件以空行结束。如果不用该命令,则将会以终端作为绘图请求来源。

RUNTPE=aaaa  读入文件aaaa作为MCNP计数数据的来源,默认的文件名是RUNTPE。如果默认的RUNTPE文件不存在,用户会被RMCTAL命令或RUNTPE命令提示。RUNTPE文件的名称可以在运行程序的执行行中使用命令mcnp runtpe=runtpe来指定文件名。

PLOTM=aaaa    将图形的图元文件命名为aaaa。默认的名称是PLOTM.ps。该文件是一个标准的PS文件(PostScript file)

COMOUT=aaaa       写入所有的绘图请求到文件aaaa中,默认的名称是COMOUTMCPLOT写入COMOUT文件的目的是用户以后可以进行同样的绘图,在下一次运行时可以使用部分或全部以前的COMOUT文件作为本次运行的COM文件。

对于output文件、PLOTM.ps文件和COMOUT文件,MCNP会选择唯一名称,以防对已有的文件覆盖写入。

绘图请求通常是通过终端的键盘输入的,但是用户也可以从一个文件输入。一个绘图请求是通过在命令提示符后的一系列绘图命令产生的。请求以ENTER键终止,但并不会立刻先于&COPLOT命令。绘图命令包括关键词,通常跟随着一些参数、空格或逗号。

几乎所有的命令都有缺省选项可用。例如:如果MCNPZ执行行运行,且RUNTPE文件中的首次计数具有多个能量箱,且在MCPLOT提示符后直接键入ENTER键,此时会出现一个计数与能量的lin-log直方图,包括误差棒和相应的标签。

B  Plot Conventions and Command Syntax绘图约定及绘图命令语法格式

1 二维绘图

坐标原点在图形左下角,水平轴叫x轴,它是一个独立变量轴,例如用户箱、栅元号或能量。竖直轴叫y轴,它是一个因变量轴,例如通量、流或剂量。每个轴都可以是线性或对数型式。

2 等高线绘图

坐标原点在图形左下角,水平轴叫x轴,它是两个独立变量的第一个变量轴。竖直轴叫y轴,它是第二个独立变量。等高线代表因变量的值,只允许线性型式的轴。每条等高线以不同的颜色标识,依据其值与极值的关系。

3 色彩填充绘图

该绘图选项与等高线图相似,只是与等高线图画z(x,y)的等高线不同,每个计数箱用一种颜色填充,颜色的选择表征该计数箱计数值的大小。该图的轴的约定与等高线图相同。该类用命令WASH调用。如果选择了两个自由变量模式(使用FREE命令),就会绘制色彩填充的图像。对射线照相计数而言,这是一个很有用的选项。采用从红到黄、绿、蓝的64色调色板。比色指数(color index)通过从Z值的最小值到最大值之间线性插值得到。

4 命令语法格式

每条命令包含一个关键字,大多数情况还跟着一些参数。关键字和参数通过空格隔开,每行不超过80个字符。逗号和等号都被解释成空格。绘图请求可以通过在上一行输入ENTER之前键入&符号延续到下一行继续输入,但是对其中的每条命令(包括关键字和参数)必须在同一行上。关键字可以简写至可以唯一标识的程度,当然必须要拼写正确。字符“current plot”意味着该图形是由当前输入的绘图命令定义的,这也许并不是显示在屏幕上的那个图形界面。只有下面C部分中标注*的命令可以用在首次以COPLOT命令作为绘图请求之后,因为其它的命令要么会影响图形的构架,要么只适用于等高线或三维绘图。

C  Plot Commands Grouped by Function根据功能分类的绘图命令

1 设备控制命令

通常情况下,MCPLOT在用户终端画图,而不在其它地方操作(译者注:指生成图元文件等)。通过以下命令,用户能够指定是要在终端画图还是发送至图像图元文件或PS文件,后者运行用户使用图形实用程序发送图像到其它的图形设备上进行后处理。

TERM n      该命令根据n值设置输出设备类型:

0        用于终端没有图形处理能力。不会在终端上画图,将所有图形发送至图形文件。TERM 0等同于在MCNP执行行中使用NOTEK

1        在下一次绘图请求时恢复可视化绘图窗口;

FILE aa       根据参数aa的值发送或者不发送图形到PS文件PLOTM.ps。图形文件直到第一次FILE命令输入之后才会被创建。并且FILE命令在NOTEKTERM 0情形下不起作用aa的运行值包括:

           即什么都不填,只是当前图形被发送至图形文件;

ALL              当前图形及后续图形被发送至PS文件(PLOTM.ps),直到下一个FILE命令输入;

NONE    当前图形及后续图形都不会被发送至PS文件(PLOTM.ps),直到下一个FILE命令输入;

2 通用命令(General Command

*  &                   从下一行连续读入命令,&符号必须是一行的最后一项;

*  COPLOT     根据目前为止输入的命令画一条曲线,并且为一条或多条曲线在同一张画布上绘制保留绘图界面;

FREQ  n           指定每隔n个历史调用MCPLOT的间隔。而在KCODE模式计算中,间隔是每n个循环。如果n是负数,则间隔是以CPU时间的分钟为单位。如果n=0,则MCNP在运行过程中不调用MCPLOT。默认情况就是n=0

RETURN           如果是正在运行的MCNP或是PLOT调用了MCPLOT(译者注:PLOTMCPLOT可以互相调用),该命令将控制返回到调用程序。其它情况下,该命令不起作用;

PLOT                 调用进入几何绘图,但是当从MCTAL文件进行的MCPLOT时,该命令不起作用;

PAUSE n           COM=aaaa配合使用。保持每个图片n秒钟。如果不提供n值,则每个图片保持到回车键键入;

*  END             终结MCPLOT

* = 可与COPLOT配合使用

3 查询命令

* OPTIONS or ? or HELP      显示MCPLOT命令的关键字列表;

* STATUS           显示画图参数的当前值;

* PRINTAL        显示当前RUNTPE文件或MCTAL文件中的计数数目;

* IPTAL              显示当前计数的IPTAL数组信息。该数组(参见附录E 31页)告诉我们当前计数(都是8维计数)每一维的元素的个数;

PRINTPTS          显示当前绘图中点的x-y坐标。该命令不适用于COPLOT、等高线和色彩填充绘图模式。

* = COPLOT配合使用

4 文件处理命令

* RUNTPE  aa  n   读取RUNTPE文件aa中的第ndump。如果参数n省略,则读取最后一个dump

* DUMP n                  读当前RUNTPE文件的第ndump

* WMCTAL aa          将当前RUMTPE文件中当前dump的计数数据写入MCTAL文件aa中;

* RMCTAL aa           读取MCTAL文件aa

* = 可与COPLOT配合使用

5 参数设置命令

针对一条曲线或者是一个图形的输入参数会对其后续的曲线或图形依然起作用,直到用RESET命令或者是用同样的命令的新参数覆盖,或者用互不相容的命令,亦或者是用FREE命令来重置许多参数。但需要注意的是有两个例外:FACTORLABEL只对当前曲线或图形起作用。对互不相容的命令,例如,BAR命令会关闭HISTPLINEAR两个命令。

a 一般命令

*  TALLY  n   定义计数n为当前计数。参数n指定了当前RUNTPE文件或MCTAL文件所描述的问题输入文件中Fn计数卡。默认的情况是第一个计数,这将是最小的中子计数,如果没有中子计数则显示最小光子计数,如果仍然没有计数则显示最小电子计数。

*  PERT  n      绘制一个关于计数的扰动。参数n表示在PERT n卡中的n值。如果是PERT 0则会重置PERT n

NONORM         禁止箱归一化。如果独立变量是方向角余弦、能量或者时间,默认的情况是,二维绘图要用计数值除以计数箱的宽度。但是,具体的请参考MCTAL文件的描述,位置在后面的Section II.D在等高线图和色彩填充图中是不进行箱归一化的

*  FACTOR  a  f  s    a轴数据乘以因子f再加上sa可以是xyzs是可选的,如果省略则被设置为0。对二维图形的初始曲线,会重置轴的限值到默认值,该命令只对当前曲线会图形有效。

*  RESET  aa 重置命令aa的参数到默认值。aa可以是参数设置命令,COPLOT或者是ALL。如果aaALL则所有的参数设置命令将会重置为默认值。在COPLOT命令后,只有COPLOTALL以及标注*的命令可以被重置,并且如果设置为COPLOTALL表示下一绘图为初始绘图。

* = 可与COPLOT配合使用

b 标题命令(其中的双引号是必须的)

TITLE  n  “aa”     使用字符aa作为位于图形顶部的主标题的第n行。n的允许值为12aa的最大长度为40个字符。默认值是FC卡的内容,如果有的话,否则,会命名为当前RUNTPE文件或MCTAL文件的名字+计数名称。KCODE图形有其自己独有的默认标题。

BELOW                   放置标题到图形之下,而不是在其上。(具体功能待实验验证

SUBTITLE  x  y  “aa”              在位置坐标(x,y)处插入副标题,名称为aa。可以放置在任何位置,包括轴和图形之间的边缘区域。

XTITLE  “aa”        定义aax轴的名称。默认是x轴变量的名称。

YTITLE  “aa”        定义aay轴的名称。默认是y轴变量的名称。

*  LABEL “aa”      定义aa为当前曲线的标签。放在图例的位置,默认参数是空格,表示用RUNTPE文件或MCTAL文件的名称作为图例标签。

* = 可与COPLOT配合使用

c 指定绘图内容的命令

MCNP的计数是按照8个不同的独立变量进行分箱计数的。由于只有其中的1个或2两个可以用于一个绘图,它们需要被指定为独立变量,剩余的就变成了固定变量,所有的固定变量(即箱的数目)可以自己定义或使用默认设置。每个固定变量的默认值是总的箱,如果存在的话;否则使用第一个计数箱的值。

 

FREE  xy           设置变量xy为空)或者是变量xy为一幅图中的独立变量。如果只指定了x,则绘制一个二维图形;如果xy都指定,则绘制等高线图或者是色彩填充图,这同时依赖于WASH卡是否生效。参见FIXED命令中的一系列可用于xy的符号。xy的缺省值是E,给出以E(能量)为自由变量的二维图形。FREE命令会重置XTITLE, YTITLE, XLIMS, YLIMS, HIST, BARPLINEAR命令的各参数为缺省值。

* FIXED q n      设置n为固定变量q的箱号。q的可用的符号及对应的箱的说明如下:

                           F      栅元、曲面或探测器;

                           D     总共(total)VS.直接(direct),或是flagged VS. unflagged

                           U     用户自定义

                           S      分段

                           M     乘子(multiplier)

                           C     方向角余弦(cosine)

                           E      能量

                           T      时间

SET   F  D  U  S  M  C  E  T

                           定义哪些变量是自由的,并且定义固定变量的箱的数目。该命令是命令FREE和若干FIXED命令的集合形式。每个参数的值可以是箱号(对应的变量就被定义为固定变量)或者是*(对应的变量就被定义为自由变量)。如果只有一个*,则表示绘制二维图。如果有两个*,则绘制等高线图或是色彩填充图。它对其它命令的覆盖影响与FREE命令相同。

TFC  x              绘制当前计数的计数波动表(tally fluctuation chart)。这里的独立变量是NPS(运行历史数)。x的允许值可以是以下的符号(符号前加?表示该项目可以从MCTAL文件绘制):

                           ?M  平均值(mean)

                           ?E  相对误差(relative error)

                           ?F  品质因子(figure of merit)

              L      201最大计数值 vs. x,如果有NONORM命令,则表示频率vs. x

                           N     累积号分额(cumulative number fraction)f(x) vs. x

                           P      概率f(x) vs. x,如果有NONORM命令,则表示频数 vs. x

                           S      高计数的斜率(SLOPE)作为NPS的函数

                           T      累积计数份额(cumulative tally fraction)f(x) vs. x

                           V     方差的方差(VOV)作为NPS的函数

                           1-8   表示f(x)1-8阶矩,即f(x)*x1 to 8 vs. x,如果有NONORM命令,则表示f(x)*x*x1 to 8 vs. x

                           < xmlnamespace prefix ="st1" ns ="urn:schemas-microsoft-com:office:smarttags" />1c-8c      表示1-8阶累积矩,即f(x)*x1 to 8 vs. x

KCODE  i        KCODE循环为独立变量,从第一个估计开始绘制单个估计量,但是平均估计col/abs/trk-len绘图是从第四个有效循环开始。根据i值绘制有效增殖系数keff或瞬时移出寿命(prompt remove lifetime)图形:

                           1      k(碰撞估计collision)

                           2      k(吸收估计absorption)

                           3      k(径迹长度估计track)

                           4      瞬时移出寿命(碰撞)

                           5      瞬时移出寿命(吸收)

                           11-15      对应于i-10的量在目前为止循环数的平均值

                           16    keff的平均估计col/abs/trk-len,且包含一个估计的标准偏差

                           17    keff的平均估计col/abs/trk-len,且包含一个通过跳过多个循环而估计的标准偏差,因此,少于10个有效循环将得不到绘图

                           18    keff的平均估计col/abs/trk-len的品质因子

                           19    keff的平均估计col/abs/trk-len的相对误差

* = 可与COPLOT配合使用

d 截面绘图命令

* XS  m            根据m值绘制截面(译者注:能量的函数):

                           Mn   输入文件中的材料卡中的序号。例如,执行命令XS  M15则绘制第15号材料的截面图。如果指定的材料号不存在,则会提示可用的材料号列表。

                           z      绘制某个核素ZAID的截面。例如,执行命令XS  92235.50C则绘制50C数据库中的核素92235的截面。如果指定的核素在输入文件中没有用到,则会提示可用的核素列表。

                           ?      打印出截面绘图的指导命令及说明。

* MT  n            绘制当前材料或核素的第n个反应截面。缺省值是总截面。可用的反应类型列于附录GB-1页(译者注:该页码只表示英文版的页码)中。如果输入了一个无效的反应号,则在数据文件中会列出可用的反应类型。

* PAR  p          绘制粒子类型p的截面,p值可以是N, P, E,分别表示材料Mn的中子、质子或电子截面。如果XS=Mn,则缺省绘制源粒子;如果XS=z,则缺省的粒子类型由数据库的类型决定。例如,92000.01g就定义了PAR=P。该命令必须在输入的第一行。

* = 可与COPLOT配合使用

e 指定二维绘图格式的命令

LINLIN             设定为线性x轴和线性y

LINLOG            设定为线性x轴和对数y轴,该设置是缺省值

LOGLIN            设定为对数x轴和线性y

LOGLOG          设定为双对数坐标轴

XLIMS  min  max  nsteps

YLIMS  min  max  nsteps

                           定义xy轴的下限和上限,及步长。nsteps选项对线性轴是可选的,对对数轴无效。当用户不指定时,两个坐标轴的上下限及步长是由MCPLOT的一种算法确定的(译者注:具体的算法没说)。

SCALES  n      根据n值显示或隐藏刻度和网格线:

                           0      没有刻度和网格线

                           1      在边缘显示刻度值,无网格线(缺省)

                           2      在边缘显示刻度值,在图形上显示网格线

*  HIST            绘制直方图。如果独立变量设定为方向角余弦C、能量E或时间T,则默认绘制直方图。

*  PLINEAR    绘制分段-线性图(piecewise-linear plot)。如果独立变量不是设定为方向角余弦C、能量E或时间T,则默认绘制分段-线性图。

*  BAR             绘制柱状图。

*  NOERRBAR       不绘制误差棒。默认情况是要绘制误差棒的。

*  THICK  x   设置图中曲线的线粗为x值。合法取值范围为0.01~0.10,默认值为0.02

*  THIN           设置图中曲线的线粗为允许的最小值,即0.01

LEGEND  x  y       根据xy的值确定是否显示图例:

                           x,无y     显示图例于通常位置(缺省)

                           x=0,无y     不显示图例

                           xy均赋值  则显示曲线标签的部分于坐标(x,y)处

* = 可与COPLOT配合使用

f 指定等高线图格式的命令

CONTOUR  cmin  cmax  cstep  %

                           定义等高线图形的最小值、最大值与步长分别为cmin, cmax, cstep。最后的百分号%是可选的,如果有该选项,则前面三个参数被解释为因变量最小值和最大值的百分数。默认值分别是5  95  10%

WASH  on/off   显示或隐藏等高线图形的色彩填充功能,默认是不显示的。

 

D.  MCTAL文件

 

该文件包含了RUNTPE文件中一个dump的计数数据信息。它的文件生成方式包括:MCNPMCRUN模式、MCPLOT模式、其它的自编程序或者甚至可以是手动地写入(为与MCNP的其它计数信息一起进行共同画图而发送至MCPLOT)。

当通过MCNP写入生成时,MCTAL文件具有如下的格式。但是只有那些对真实系统起决定性作用的信息是必要的。另外,数值项也并不需要完全按照其格式书写,只要顺序正确且用空格分开,且不要有多余的嵌入空格。例如给定某个量与能量的关系表,用户可以按以下格式简单写入MCTAL文件:

E 7 1

.2 .4 .7 1 3 8 12

VALS

4.00E-5 .022 5.78E-4 .054 3.70E-5 .079 1.22E-5 .122

7.60E-6 .187 2.20E-6 .245 9.10E-7 .307

如果需要多个自变量,则需要添加相应的行及数值。

MCNP生成的MCTAL文件的结构

kod, ver, probed, knod, nps, rnr                    (2A8, A19, 15, I11, I15)

kod  程序名称,即MCNP

ver  版本号,如:4A4B4C5…

probid  问题运行的时间和日期(如果机器标识符可用的话也会包括在内)

knod  dump

nps  运行历史数

rnr  用到的伪随机数号

问题标识行(即问题输入文件的第一行)1-79列之后有一个空格

以后再继续翻译

E  COPLOT举例

runtpe a  coplot  runtpe b

假设所有的参数设置命令在之前已经定义好。上述命令将会在同一张图上放置两条曲线。第一条曲线将显示RUNTPE文件a的计数数据,第二条曲线将显示RUNTPE文件b中相同计数号的计数数据。除非重置,否则MCPLOT将继续从RUNTPE文件b中读取数据。例如,下述命令改变x轴上下限,定义计数号11为当前计数,绘制第一条曲线来自于RUNTPE文件b,第二条曲线来自于MCTAL文件aux的计数号41,第三条曲线来自于RUNTPE文件a的计数号1,如果不重新设置,则后续绘图的数据就会来源于RUNTPE文件a中:

xlims min max  tally 11  coplot  rmctal aux  tally 41 coplot  runtpe a  tally1

下述命令将发送一个两条曲线的画面到一个postscript文件中:

tally 24  nonorm  file  coplot  tally 44


  评论这张
 
阅读(3426)| 评论(3)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018