《mcnp从入门到精通讲座》由会员分享，可在线阅读，更多相关《mcnp从入门到精通讲座（88页珍藏版）》请在金锄头文库上搜索。
1、-第一讲： 如何运行MP程序2-Ej4I 执行命令语句 Tw;3_Lj *Y %&n g m0 )r mp i=(输入文件名) o=(输出文件名) 选项 sEGO2*eI L4mL, 选项 i|p|*|r|z SfaQvstN i=process input file 过程输入文件 dE4L=sTEsy p=plot geometry 几何画图 U z=plot tally results,plot cross-sections 画记数结果,画截面。 2O&NdPZk 3M|m后敲入， IDB+%*l*S C:mpbinmp i|p|*|r|z i=(输入文件名) o=(输出文件名)=E.t*
2、= 其中为回车键。 V9O V% 缺省文件名： :7!0OVQla inpascii input file ascii码输入文件 ,+meTvn outpascii output file ascii码输出文件 7BK0s*O runtpebinary restart file 二进制转储文件 Z,qs 举例如下： 1DEmp i*r i=greatwall name=greatwall.1 O4cr*MCb5 其中greatwall是你自己编写的源程序输入文件。 19Mu61 这个命令行语句，将读入输入文件greatwall，调用*S数据库，并且运行mp。 ktFhc3);! .1是你自己选
3、择的输出文件名，每次输出必须更改输出文件名，或者不更改但删除上一次的输出文件名，否则程序运行出错。你也可以改成.2,.8,等等，任你选择。 p) ;eJtV/ 这样运行的结果，就会有如下输出： /H( 2IHC ,9h greatwall.1.oascii码输出文件 h5C|M| greatwall.1.r转储文件 v/68*,z 其中，转储文件也可以作为下一次运行的输入文件，进展接续运行。关于接续运行，将在以后的讲座中详细介绍。 LcB*dsa( b7zJUz 你自己定义的文件名，将使你能够跟踪上述的辅助文件。 OWjk=u2Lz sl|E 练习： a ueL/dgN 众坛友，请运行下面的g
4、reatwall文件，暂时不要计较文件里的容，以后的讲座，会陆续予以讲解。程序的目的是，画出在一个立方体当中有一个球体。 75f2; message: datapath=C:mpbin R7$:g title -greatwall: U sphere in Al shpere in water cube- ein4ojLRb|7Ts 30 3 -1 +2 -3 imp:n=1 $ water uw 5 c LR.Hh 40 0 +3 imp:n=0 $ outside world 2 mJ+v c no fUD. S4vpY DeN 1 sph 0 0 0 5 $ Sphere at orig
5、in radius 2 cm 7ePqmB. c pTpK ) 2 sph 0 0 0 5.02 $ outter Al spherical shell 0.5 cm thick qUjmB sB c jG*O:s O 3 bo* +20 -20 -20 0 40 0 0 0 40 -40 0 0 $ surrounding H2O ap&rTu J)-+7$v c * Source Card* bu:S: c (d(hR0HKE sdef $ default source ZAv,*5& c n*qu, c * Material Cards * z9 )VIe c S*ua36 m3 100
6、1 +1 8016 +2 $ water zUKm*y c BRlT7grgq c * Print first 50 histories * =atVp/y print 110 &Qk7L ! c db 17j 1 j :*D0 prdmp 2j 1 1 f W - c $lLz 3YS fc4 flu* integrated over the Pu/Al interface mcEl f4:n 10 p(Y*3$* f14:n 20 0URji|* c * time and nps cut cards * hS)aFV nps 10 b*lk2oMa/ ctme 1 J4WzA%i .7l&
7、1C)i End of File EJsM(iGM 你可以用Wordpad或者Notepad翻开任何一个你已经有的INP文件，删去其所有容，拷贝上述程序，另存为greatwall文件。注意，不要敲Tab键，否则程序出错。 M72Gp S8*t 尝试下面的命令，看看你都得到了什么输出结果，以及输出的文件都在哪里？ /emp i=greatwall Bs O+NP 其中为回车键。 Zo g= (2)再次执行一遍上述命令 kzt(i Y_6 C:mpbinmp i=greatwall a=h *J1O (3)执行 8 0Gn%1A9 C:mpbinmp i=greatwall name=greatw
8、all.1 =Urw&H8 (4)执行 qI*lW C:mpbinmp i=greatwall name=greatwall.2 *a*rM6P (5)执行 NhaImp i*r i=greatwall name=greatwall.2 L IZRoG8 YKO)f5 下面是mp程序运行的构造框图，从中你可以看到，用户输入文件被mp.e*e执行后，要根据执行命令语句的设置，来调用*S数据库，然后计算输运过程，把输出文件放到所执行的子目录里。 %*J k*W: W9D863Y =j*8.K4 我不知道，国怎样翻译*S数据库，索性就把*S放到这里，不再忽发奇想，闭门造车，以免毁人不倦。如果有哪位朋
9、友知道的话，请跟帖加以说明，非常感。 M3fmAP u%$Zqee 这一讲，就说到这里。欲知后事如何，请听下一讲分解：mp输入文件。*_*第二讲，MP输入文件 LtH j(1) MP中物理量的默认单位KTt+-vP g4=C1 长度：cm ysZ(*K n( 通量：MeV |m 5(S 时间：刹 shake (10-8秒) )Nq*=ms(! 能量：MeV Zo-,TKgY 温度：MeV (kT) 0ETT/)z 原子密度：1024个原子/cm3 *E_KW1RV 质量密度：g/cm3u:jw=f 截面：巴 barns (10-24/cm2) M*bsA/Z 加热量：MeV/collision
10、 u;(4sa cLp_ 此外， 原子质量按照中子质量为1.0计算，这种单位下阿佛伽德罗常数是0.59703109; 程序运行时间以分钟为单位。 50dN(;p MP的源代码是用FORTRAN语言编写的。 *v *ci W (1q(6! (2) 输入INP文件的根本形式5*|$q kI /s%-c!o INP文件由一些被空行分隔的由一些被空行分隔的输入块组成，主要的输入块是信息块、标题和栅元块、曲面块和数据块等。输入块又由一些被称为卡的输入行组成。论坛的朋友可以参照比照一下，上一讲中的greatwall程序容。 V*5PF INP文件的格式如下： =+iN*Sc 0fBPHWr 信息块 (可选
11、项) !C 空行分隔符 (可选项) aCH;l+U 标题卡 7TmAllIk%_:CC 空行分隔符 Fl O,| 曲面卡 + r*_ !5h-$; 空行分隔符 g18zoLZ 数据卡 =bPcC=3b (pM& eow 空行分隔符 -DJ ,f*$ 其它 (可选项) /&G|.C* 6Y6t.j0vN. 说明如下： _*g .U=u 信息块的第一行，必须在它的18列填写MESSAGE:，后面跟着用空格分隔的参数项。其后的980列和空行分隔符之前各行的180列都看作信息块容。可用A=B参数项更改输出文件名，如OUTP = MYOUT。信息块是可选的。信息块提供应MP一个执行信息的方法。 4_CL
12、1g 在信息块之后的第一行是问题的标题卡，它仅限于一行，占用180列，可以是任何信息，将作为OUTP文件中各个输出表的标题被复制。 gaNe 用户在栅元块和曲面块中描述问题的几何。栅元由栅元卡描述。空间必须由彼此相邻的栅元填满，栅元之间不能重叠，也不能出现无栅元的空区，否则会出现错误。构建栅元的曲面由曲面卡定义，曲面卡在曲面块中给出。曲面卡和栅元卡的填写方法，将在以后的讲座中予以介绍。 uit.r8l 曲面块之后是数据块，在数据块中用户描述源、记数方式、材料等。数据卡在以后的讲座中予以详细介绍。 *ws u8 数据卡后不管有无空行分隔符均可以运行，不同之处是，如果数据卡后面有空行分隔符，则MP将不再读后面的附加行(如果附加行存在的话)。这对*些用户喜欢保存附加的容是很方便的，如改变这个问题的说明或者与这个输入文件本身有关的信息，则，这个分隔符就能防止读入这些附加信息。用户可以把希望保存的一些附加信息写在数据块的空行之后，MP会将它们复制到OUTP文件末尾。 (RrC5 Ktt(l-e+ (3)卡片格式S
