mcnp 4b 教程范本

1、MCNP4B教程教程第１頁(yè)2006年6月MCNP輸入的描述輸入的描述MCNP的輸入包括幾個(gè)文件，但主要一個(gè)是由用戶填寫(xiě)的INP(缺省文件名)文件，該文件包括描述問(wèn)題所必需的全部輸入信息。對(duì)于任何一個(gè)特定的問(wèn)題，只需用到INP的全部輸入卡片的一小部分“卡”這個(gè)詞描述的是一個(gè)最多可達(dá)80個(gè)字符的輸入行。有些MCNP的輸入項(xiàng)存在最大維的限制，用戶可以通過(guò)修改代碼來(lái)改變它們的最大值。MCNP的所有功能都應(yīng)謹(jǐn)慎使用并應(yīng)具備相應(yīng)的知識(shí)。尤其在探測(cè)器

2、的調(diào)試和降低方差的實(shí)現(xiàn)方面，因此，在運(yùn)行MCNP之前，建議閱讀第二章中相應(yīng)的部分。Ⅰ.INP文件文件INP文件有初始運(yùn)行及接續(xù)運(yùn)行的兩種形式，它們都可包括一個(gè)可選擇的信息塊，用以替換或補(bǔ)充MCNP的執(zhí)行命令行信息。A.信息塊信息塊用戶可以在INP文件中標(biāo)題卡的前面有選擇性的放一些信息卡，在沒(méi)有執(zhí)行行信息的計(jì)算環(huán)境下，只有信息塊能給出MCNP一個(gè)執(zhí)行信息。這是一個(gè)避免重復(fù)輸入一些信息的常用方法。信息塊用字符串“MESSAGE：”作為開(kāi)始，

3、并且限定在180列，字符可以是大寫(xiě)、小寫(xiě)或大小寫(xiě)混合，空行定界符前所有的卡都作為繼續(xù)卡。信息塊中$和&符號(hào)都是結(jié)束行標(biāo)志，在標(biāo)題卡之前用一個(gè)空行分隔符結(jié)束信息塊。信息塊上各部分的語(yǔ)法和在第一章所討論的執(zhí)行行信息一樣。信息塊上各部分的意義和執(zhí)行行信息是一樣的，但執(zhí)行行信息與信息塊中所指定的信息有沖突時(shí)，則執(zhí)行行信息優(yōu)先于信息塊上的同樣信息，特別地：a.在信息塊上，INP=文件名是不合法的，只能在MCNP的執(zhí)行行改變INP文件的名字。b.在

4、A=B（文件名替換）的情況下，如果A=這一結(jié)構(gòu)在執(zhí)行行信息上出現(xiàn)，也在信息塊上出現(xiàn)，則信息塊上這一項(xiàng)被忽略。c.如果在執(zhí)行行信息中有任意一個(gè)程序模塊執(zhí)行選項(xiàng)（如IP或IX），則信息塊的全部執(zhí)行選項(xiàng)被忽略。d.在執(zhí)行行上的任何關(guān)鍵詞項(xiàng)都將使得在信息塊上的相應(yīng)項(xiàng)被忽略。例如，在執(zhí)行行信息上的C7（指定接續(xù)第七次轉(zhuǎn)儲(chǔ)）將使信息塊上的C4被忽略。e.如果C或Cm在信息塊上出現(xiàn)，不在執(zhí)行行上，則這個(gè)運(yùn)行仍是一個(gè)接續(xù)運(yùn)行。B.初始運(yùn)行初始運(yùn)行用這種

5、形式建立一個(gè)蒙特卡羅問(wèn)題（描述問(wèn)題的幾何結(jié)構(gòu)、材料、記數(shù)要求等），如果需要，便可直接運(yùn)行。該初始運(yùn)行輸入文件的格式如下：信息塊可選空行分隔符可選標(biāo)題卡柵元卡┆空行分隔符曲面卡MCNP4B教程教程第３頁(yè)2006年6月如果上述項(xiàng)都沒(méi)有改變（并且計(jì)算環(huán)境允許執(zhí)行行信息），則不需要接續(xù)輸入文件；僅需要運(yùn)行文件RUNTPE及執(zhí)行行上的C選項(xiàng)。例如，如果一個(gè)作業(yè)運(yùn)行一分鐘左右，因時(shí)間限制中斷，但又想多運(yùn)行一些粒子，這時(shí)可簡(jiǎn)單地執(zhí)行命令行上的C或者C

6、N信息，作業(yè)將從中斷的地方重新開(kāi)始運(yùn)行，直到另一個(gè)新給的時(shí)間限制或者粒子截?cái)嗷蛘呷藶橥Ｖ?。這個(gè)例子假定來(lái)自初始運(yùn)行的一個(gè)名為RUNTPE覆蓋文件在當(dāng)前目錄下。完整的接續(xù)運(yùn)行執(zhí)行行選項(xiàng)為Cm或者CNm，m特指從轉(zhuǎn)儲(chǔ)文件RUNTPE接續(xù)運(yùn)行的覆蓋號(hào)，如果m沒(méi)有特指，最后覆蓋號(hào)被默認(rèn)使用。如果前次運(yùn)行是因處理完指定的總粒子數(shù)（NPS）而終止的，那么在接續(xù)運(yùn)行時(shí)，必須給出NPS卡，增大需處理的總粒子數(shù)NPS，NPS包括接續(xù)運(yùn)行和初始運(yùn)行的所有的

7、初始粒子歷史，接續(xù)運(yùn)行中的CTME是指這個(gè)運(yùn)行所花的時(shí)間，而不是總時(shí)間。為了運(yùn)行更多的KCODE循環(huán)，只有第四個(gè)輸入KCT起作用，和NPS卡類似，KCT包括以前的循環(huán)在內(nèi)的所有循環(huán)。在接續(xù)運(yùn)行文件中，NPS卡上給出是負(fù)值時(shí)，將只是對(duì)前次轉(zhuǎn)儲(chǔ)的中間結(jié)果進(jìn)行編輯輸出，并不再多計(jì)算歷史。當(dāng)原有的輸出文件已經(jīng)丟失?；蛘呦胪ㄟ^(guò)PRINT卡或FQ卡改變輸出內(nèi)容時(shí)，使用這一方式是十分方便的。但在初始運(yùn)行使用FILES卡時(shí)要謹(jǐn)慎小心。見(jiàn)3119頁(yè)。D.

8、卡片格式卡片格式INP輸入文件的每一行（稱之為一個(gè)卡）都限于使用180列。英文字符可以用大寫(xiě)、小寫(xiě)或大小寫(xiě)混合。大部分輸入卡按行填寫(xiě)；然而，對(duì)數(shù)據(jù)卡允許按列格式填寫(xiě)。注釋可以加入到任意一個(gè)輸入卡中。$符號(hào)（美元符號(hào)）為它所在那行數(shù)據(jù)的結(jié)束符，在$符號(hào)后面的內(nèi)容作為注釋?？招锌梢员挥米鞣纸绶徒K止符，數(shù)據(jù)項(xiàng)用一個(gè)或多個(gè)空格分開(kāi)。輸入文件中，在標(biāo)題卡之后及最后的空行結(jié)束卡之前，任何地方都可插入注釋卡。注釋卡的第1—5列必須有字母“C”，且后

9、面至少要有1個(gè)空格。注釋卡僅在輸入內(nèi)容的原形輸出部分打印出，不再在輸出文件的其它任何部分出現(xiàn)。FCn卡是作為用戶注釋用，它作為記數(shù)類型n的表頭文字，作為記數(shù)的標(biāo)題。SCn卡也是作為注釋用，它作為源概率分布n的表頭文字。標(biāo)題卡被限定在1—80列的范圍內(nèi)使用，可以包括任何用戶所期望的信息甚至也可以是空白。它通常描述的是特別問(wèn)題的信息，但要注意其它地方使用空格或空行表示分界符和終止符。1.行輸入格式行輸入格式柵元卡，曲面卡，及數(shù)據(jù)卡都必須從第

10、1～5列開(kāi)始，數(shù)據(jù)輸入緊跟在卡名或卡號(hào)和指定粒子之后，而這些數(shù)據(jù)項(xiàng)又被一個(gè)或多個(gè)空格分開(kāi)。如果1～5列為空白，則表示它是前一張卡的繼續(xù)卡。而至少一個(gè)空格以后的符號(hào)&表示隨后的卡是繼續(xù)卡，繼續(xù)卡數(shù)據(jù)限定在180列。完全空白卡為在兩個(gè)輸入塊之間的定界符。一個(gè)數(shù)據(jù)項(xiàng)必須在一張卡上寫(xiě)完，不得跨到下一張卡上。對(duì)任何給定的帶有粒子標(biāo)識(shí)符的類型卡只有一張。需要整數(shù)輸入的地方必須填寫(xiě)整數(shù)。其它數(shù)據(jù)可填寫(xiě)為整數(shù)或浮點(diǎn)數(shù)，MCNP能正確地讀入。實(shí)際上，非整

11、數(shù)數(shù)據(jù)能用任何FTRAN語(yǔ)言接受的E格式描述。為書(shū)寫(xiě)方便，可以使用四項(xiàng)書(shū)寫(xiě)功能：AnR功能功能，表示將它前面的一個(gè)數(shù)據(jù)再重復(fù)n遍（R表示重復(fù)）。例如：24R和22222是一樣的。BnI功能功能，表示在與其前、后相鄰的兩個(gè)數(shù)之間插進(jìn)n個(gè)線性插值點(diǎn)（I表示插值）。例如：1.52I3.0就相當(dāng)于給出1.52.02.53。在XnIY這種結(jié)構(gòu)中，如果X和Y是整數(shù)，并且YX剛好是n1的倍數(shù)，則產(chǎn)生標(biāo)準(zhǔn)的整數(shù)插值。否則產(chǎn)生實(shí)數(shù)插值。而y被保存到所有的