【正文】 須沿著棱柱的一個(gè)方向或者另一個(gè)方向增加。然后按照順序的規(guī)定在柵元卡上填寫界定（0,0,0）元素的曲面。對(duì)一個(gè)六面體柵格柵元，在第一個(gè)列的曲面的外面是（1,0,0）元素，在第二個(gè)列的曲面的外面是（1,0,0）元素；然后按照那個(gè)順序依次為（0,1,0）元素，（0,1,0）元素，（0,0,l）元素，及（０,０,1）元素，如圖23(a)所示。這個(gè)方法提供柵格排列的順序，因此當(dāng)你指定元素（7,9,3）時(shí)，程序就會(huì)知道這個(gè)元素是哪一個(gè)。對(duì)六棱柱柵格柵元，在第一個(gè)列的曲面外面是（l,0,0）元素，在第二個(gè)列的曲面外面是（1,0,0）元素；然后依次是（0,l,0），（0,1,0）,（l,l,0）,（l,1,0），（0,0,1）,（0,0,l）,如圖23(b)所示。描述（0,0,0）元素的最后兩個(gè)曲面必須是棱柱的上、下底面。(a)六面體柵元排列順序 (b)六棱柱柵元排列順序圖23柵元排列順序六面體不一定是直角的，六棱柱也不一定是規(guī)則的，但是由它們構(gòu)成的一些柵格必須正好填滿空間，這意味著柵格的對(duì)邊必須是相等并平行的。六面體柵格柵元可以在一個(gè)或兩個(gè)方向上是無(wú)限的。六棱柱柵格柵元可以沿著棱柱的長(zhǎng)度方向無(wú)限。橫截面必須是凸的。這個(gè)柵格是左旋或右旋都可以。一個(gè)柵格必須在它的Universe里。一個(gè)問(wèn)題的實(shí)際領(lǐng)域本身可以是一個(gè)柵格。如果一個(gè)粒子離開實(shí)際領(lǐng)域的最后一個(gè)柵元（限制柵格的范圍），這個(gè)粒子死亡。4） FILL卡FILL卡上的非零項(xiàng)指明用來(lái)填充相應(yīng)柵元的Universe號(hào)。在U卡上相同號(hào)標(biāo)記組成這個(gè)Universe的一些柵元，F(xiàn)ILL項(xiàng)后面的括號(hào)里可放有變換號(hào)或變換本身。這個(gè)變換是被填充柵元的坐標(biāo)系和在輔助坐標(biāo)系下用來(lái)填充柵元的Universe之間的變換。如果沒(méi)有指定變換，則這個(gè)Universe繼承被填充柵元的變換。如果旋轉(zhuǎn)角是以度為單位而不是余弦，則使用*FILL。在輸入文件的數(shù)據(jù)卡部分不能同時(shí)有FILL卡，又有*FILL卡。如果要填寫一些以度為單位的角度而又要按余弦來(lái)填寫一些角度，那么所有FILL及*FILL數(shù)據(jù)都必須放在輸人文件INP的柵元卡上。如果被填充的柵元是一個(gè)柵陣，F(xiàn)ILL說(shuō)明或者是一項(xiàng)或者是一個(gè)數(shù)組。如果是一項(xiàng)，這個(gè)柵陣的每一個(gè)柵元都是用同一個(gè)Universe填充。如果是一個(gè)數(shù)組，填充由這個(gè)數(shù)組所覆蓋柵陣的那部分；柵陣的其余部分不存在?？捎貌煌腢niverse去填充這個(gè)柵陣的各種元素。用一個(gè)柵格填充一個(gè)柵元的數(shù)組說(shuō)明是：三個(gè)尺寸，緊接著是這個(gè)數(shù)組值的本身。這三個(gè)尺寸說(shuō)明定義了三個(gè)柵格指標(biāo)的范圍。它與 FORTRAN的 DO循環(huán)一樣，下界和上界必須用正數(shù)、負(fù)數(shù)或者零表示，并用：號(hào)隔開。每個(gè)柵格元素的指標(biāo)是由它相對(duì)于（0,0,0）元素的位置確定的。在LAT卡部分已強(qiáng)調(diào)，柵格柵元曲面說(shuō)明的順序就定義了柵陣元素的順序。在該柵陣柵元卡上前兩個(gè)曲面定義第一個(gè)柵格指標(biāo)的方向。指標(biāo)的數(shù)值范圍取決于這個(gè)柵格的（0,0,0）反素的位置在哪。例如，5:5， 0:10，及10:6都是定義11個(gè)元素，但（0,0,0）元素的位置的是不同的。在柵格柵元卡上第三和第四列的曲面定義第二個(gè)柵格指標(biāo)的方向。數(shù)組大小遵守尺寸說(shuō)明，在這個(gè)數(shù)組中的每個(gè)元素。實(shí)際上只有與這個(gè)數(shù)組元素對(duì)應(yīng)的那些柵陣柵元存在。每一個(gè)數(shù)組元素的值是填充相應(yīng)柵格的Universe號(hào)。在這個(gè)數(shù)組里可以使用兩個(gè)有特殊意義的值：零意味著柵格元素不存在，這樣使用零就可以指定一個(gè)非直角陣列；如果數(shù)組值與這個(gè)柵格的Universe號(hào)一樣，就意味著與這個(gè)數(shù)組元素對(duì)應(yīng)的那個(gè)柵陣元素不用任何Universe去填充。如果零Universe（問(wèn)題的實(shí)際領(lǐng)域）是一個(gè)柵格，這個(gè)特殊值不起作用。跟帶有單項(xiàng)的FILL說(shuō)明一樣，在數(shù)組的任何一個(gè)值后面的括號(hào)里可跟有一個(gè)變換號(hào)或者變換本身。例如： fill＝0:2 1:2 0:1 4 4 2 0 4 0 4 3 3 0 4 0這個(gè)柵陣只有8個(gè)元素被填充（存在8個(gè)元素）。元素（0,l,0），（l,l,0），（l,2,0），（0,1,1）和（l,2,1）都是用Universe 4填充；元素（2,l,0）用 Universe 2填充；元素（l,l,1）和（2,l,l）用 Universe 3填充。 TRn 坐標(biāo)變換卡格式：trn O1 O2 O3 B1 B2 B3 B4 B5 B6 B7 B8 B9 Mn：變換號(hào)， 1n999。O1 ，O2 ，O3：變換的位移矢量；B1～B9：變換的旋轉(zhuǎn)矩陣；M：M=l意味著位移矢量是輔助坐標(biāo)系的原點(diǎn)在基本坐標(biāo)系定義的位置。M=l意味著位移矢量是基本坐標(biāo)系的原點(diǎn)在輔助坐標(biāo)系定義的位置。缺?。?trn 0 0 0 1 0 0 0 1 0 0 0 1 1使用：可選擇卡MCNP所有類型的曲面都有某種對(duì)稱性，當(dāng)對(duì)稱軸平行于某些坐標(biāo)軸時(shí)，即使較為復(fù)雜的曲面都是不難給出描述參數(shù)的。但是，當(dāng)曲面對(duì)稱軸不平行于坐標(biāo)軸時(shí)，對(duì)曲面給出描述參數(shù)難度就大多了。在這種情況下，用戶可以在一個(gè)較為方便的輔助坐標(biāo)系來(lái)描述該曲面，然后使用如下變換功能，將該曲面的描述轉(zhuǎn)換到問(wèn)題的基本坐標(biāo)系中。這個(gè)變換是通過(guò)坐標(biāo)變換卡TRn卡實(shí)現(xiàn)的。TRn卡指明輔助坐標(biāo)系與基本坐標(biāo)系的關(guān)系，在輔助坐標(biāo)系描述曲面是容易的，但這時(shí)在該曲面描述卡上要有變換號(hào)n，在相應(yīng)的TRn卡指定這個(gè)變換。環(huán)面（Torus）不能做坐標(biāo)變換。在MCNP的輸入文件中，基本上是在輔助坐標(biāo)系下給出較為簡(jiǎn)單的曲面描述卡。當(dāng)輔助坐標(biāo)系并不相應(yīng)于基本坐標(biāo)系時(shí)，便需使用TRn卡，用以指明輔助坐標(biāo)系與基本坐標(biāo)系之間的關(guān)系。一個(gè)問(wèn)題中最多可有10個(gè)TRn卡。B矩陣說(shuō)明兩個(gè)坐標(biāo)系軸的方向之間的關(guān)系。Bi是在基本坐標(biāo)系（X Y Z）軸與輔助坐標(biāo)系（x’y’z’）軸之間的夾角余弦（若在TRn前加了“*”，如*TRn，則由指明余弦值改為指明角度，0o＜θ＜180o＝，其排列規(guī)則為：元素 B1 B2 B3 B4 B5 B6 B7 B8 B9軸 Bi的意義與M無(wú)關(guān)。用戶不必給出所有Bi值，可用kJ跳過(guò)不指定的項(xiàng)，Bi可按以下任一種方式指定：i）Bi全部給出（9個(gè)值）；ii）給出矩陣中列或行方向上三個(gè)矢量中的兩個(gè)（6個(gè)值）。程序?qū)⒂伤鼈兊牟娉水a(chǎn)生第三個(gè)矢量；iii）給出矩陣中列與行方向上各一個(gè)矢量，其中的公共分量必須小于1（共5個(gè)值）。程序?qū)⒂脷W拉角方法填完該矩陣；iv）給出矩陣中列或行方向上的一個(gè)矢量（3個(gè)值），程序?qū)茨撤N隨機(jī)方法生成其他兩個(gè)矢量；v）沒(méi)有輸入，程序?qū)⒔⑾嗤仃?。一個(gè)矢量由三個(gè)元素組成，這些元素可以按矩陣的橫向或豎向排列，但不能走對(duì)角錢。在所有情況下，MCNP解決任何小的非正交性，并歸一這個(gè)矩陣。在歸一過(guò)程，形似（l，0，0）的矢量將保持原樣。，則給出一個(gè)警告信息。當(dāng)變換只是一個(gè)平移時(shí)，使用方式v是合適的。當(dāng)一些斜置曲面都是關(guān)于某個(gè)斜軸的旋轉(zhuǎn)面時(shí)，使用方式iv）是合適的。曲面不是旋轉(zhuǎn)面（一般情況）時(shí)，方式ii）與iii）大體是同等適用的。如果輔助坐標(biāo)系和基本坐標(biāo)系中一個(gè)是右手坐標(biāo)系，而另一個(gè)是左手坐標(biāo)系，則使用方式i）是合適的。使用坐標(biāo)變換以簡(jiǎn)化曲面的幾何描述，以及建立曲面源問(wèn)題的坐標(biāo)系統(tǒng)與寫這個(gè)曲面源文件中這個(gè)問(wèn)題的坐標(biāo)系統(tǒng)之間的關(guān)系。通常一個(gè)柵元組分將有同一個(gè)坐標(biāo)系統(tǒng)。在這個(gè)坐標(biāo)系統(tǒng)下描述這個(gè)組分的曲面，然后用TRn卡指定這個(gè)系統(tǒng)與基本坐標(biāo)系統(tǒng)的關(guān)系。例子： 17 4 px 5 tr4 7 .9 0 1 0 0 0 1 1 0 0 降低方差卡下面是定義降低方差卡的參數(shù)的卡片。助記名卡片內(nèi)容助記名卡片內(nèi)容IMP柵元重要性卡ESPLT 能量分裂與輪盤賭卡TSPLT 時(shí)間分裂或輪盤賭卡PWT光子權(quán)重卡EXT指數(shù)變換參量VECT向量輸入FCL強(qiáng)迫碰撞參量WWE粒子權(quán)窗能量WWN粒子權(quán)窗邊界WWP粒子權(quán)窗參量WWG粒子權(quán)窗產(chǎn)生器WWGE粒子權(quán)窗產(chǎn)生能量MESH網(wǎng)格權(quán)窗產(chǎn)生器的網(wǎng)格層次重要性PDnN記數(shù)的探測(cè)器貢獻(xiàn)DXCDXTRAN（球）柵元貢獻(xiàn)BBREM軔致輻射偏倚在方差降低技巧中大部分卡片都是選擇性的，但I(xiàn)MP卡和WWN卡必有其一。一個(gè)柵元或者曲面參數(shù)卡的輸入要與先前出現(xiàn)在INP文件的柵元、曲面卡相符合。為了在卡片上取到特定的柵元或曲面，需要提供適當(dāng)?shù)娜笔≈诞?dāng)作空格間隔（可使用nR重復(fù)或者nJ跳躍）。當(dāng)IMP卡上的數(shù)據(jù)項(xiàng)數(shù)不等于問(wèn)題中的柵元數(shù)時(shí)，MCNP將指出為致命性錯(cuò)誤。這些卡既需要蒙特卡羅方法的知識(shí)，也需要一定的降低方差技巧。 IMP 柵元重要性卡格式：IMP：n x1 x2…xi…xIn=對(duì)于中子是N；對(duì)于光子是 P；對(duì)于電子是E。如果不同粒子的重要性相同，也允許出現(xiàn)N，P或N，P，E 或P，E。xi=柵元i的重要性I=問(wèn)題中的柵元數(shù)缺?。喝绻侵凶又匾詥?wèn)題可缺省。如果中子重要性不是0，所有的光子?xùn)旁闹匾砸O(shè)成1。如果中子重要性不是0，光子重要性也是0。使用：對(duì)于一個(gè)柵元卡IMP卡和WWN卡必須有其中一個(gè)。柵元的重要性有兩個(gè)用途：a）當(dāng)重要性為0時(shí)，用它來(lái)終止粒子的歷史；b）用于幾何分裂與輪盤賭技巧，促迸粒子軌跡朝向幾何上較為重要的區(qū)域運(yùn)動(dòng)。在中子光子耦合輸運(yùn)問(wèn)題中，用戶可能希望對(duì)光子的分裂與中子有所不同，所以可能需要用到兩組重要性（中子和光子）。如果在運(yùn)行方式為n p的問(wèn)題中，未給出 imp:p卡，那么所有柵元的光子重要性均置為1。當(dāng)IMP卡上的數(shù)據(jù)項(xiàng)數(shù)不等于問(wèn)題中的柵元數(shù)時(shí)，MCNP將指出為致命性錯(cuò)誤。對(duì)IMP卡，nR及xM書寫功能是很有用的，例如給出一張IMP卡為： imp:n l 2 2m 0 1 20r則使柵元 1的重要性為 l，2為2，3為 4，4為 0，而 525，則都是 l。此時(shí)若有一條軌跡從塊2走進(jìn)塊3，則將被分裂成兩條軌跡，每條新軌跡的權(quán)是分裂前原軌跡權(quán)的一半。相反若軌跡是從塊3走進(jìn)塊2，則將以50％的概率終止這條軌跡；如果未被終止，則該軌跡的權(quán)將加倍。應(yīng)記住，凡進(jìn)入“0”重要性柵元的軌跡都將被殺死。對(duì)探測(cè)器進(jìn)行記數(shù)或向DXTRAN球上生成假粒子時(shí)，如果軌跡經(jīng)過(guò)一個(gè)“0”重要性柵元，則也將被殺死。當(dāng)粒子軌跡進(jìn)入真空柵元時(shí)，即使重要性增大，也不進(jìn)行分裂；但對(duì)離開真空塊的軌跡，若進(jìn)人了較低重要性的柵元，則仍將常規(guī)實(shí)行輪盤賭。 ESPLT 能量分裂與輪盤賭卡格式：esplt:n R1E1 …R20E20 n： n表示中子，P表示光子；缺?。喝绻÷粤诉@個(gè)卡片表示這些粒子不發(fā)生能量分裂。使用：可選擇卡；不能在多組計(jì)算中使用。Ri的意義取決于對(duì)于粒子類型n有沒(méi)有權(quán)窗，這將在下面進(jìn)行討論。（注意與TSPLT卡相反，ESPLT 中的Ri表示下降能量的分裂或輪盤賭。）對(duì)于大于最大能量Ei的能量，絕對(duì)重要性作用設(shè)成相同，并附加一個(gè)R0=1,代碼中的分裂與輪盤賭都要改變。規(guī)定ESPLT卡和與時(shí)間有關(guān)的權(quán)窗一起使用。對(duì)于與能量有關(guān)的權(quán)窗不推薦使用ESPLT卡，因?yàn)檫@兩個(gè)卡可能相互關(guān)聯(lián)?？梢杂靡粋€(gè)能量有關(guān)的權(quán)窗卡代替一個(gè)單區(qū)間權(quán)窗加一個(gè)ESPLT卡。沒(méi)有權(quán)窗Ei：粒子實(shí)行分裂的能量限（Mev）。對(duì)于能量下降過(guò)程中經(jīng)過(guò)Ei的粒子：Ri=（Ri 1）粒子要分裂成的軌跡數(shù)；(0Ri1)俄羅斯輪盤賭的概率。對(duì)于能量上升過(guò)程中（例如裂變）經(jīng)過(guò)Ei的粒子（Ri0參看“其他”）：1/ Ri=（1/ Ri 1）粒子要分裂成的軌跡數(shù)；(01/ Ri 1)輪盤賭的概率。其他：如果第一輸入項(xiàng)是負(fù)的，也就是Ri0，能量上升時(shí)不能做分裂或輪盤賭游戲。缺省：省略這張卡的粒子沒(méi)有能量重要性。對(duì)于問(wèn)題范圍的粒子能量分裂與輪盤賭可以使用ESPLT卡，而對(duì)于幾何作用的分裂與輪盤賭可以使用IMP卡。在MODE P E F8加上高度記數(shù)和多組問(wèn)題不能使用ESPLT卡。ESPLT卡引起的粒子權(quán)重的改變會(huì)造成中止能量和權(quán)窗能量的能量修正。卡上的輸入數(shù)據(jù)包括最多20對(duì)的能量重要性比率參數(shù)RiEi。如果能量Ei不是單調(diào)的，會(huì)發(fā)出警告信息。Ri可以不是整數(shù)也可以是0和1之間的數(shù)。如果能量下降到Ei以下，Ri 1，粒子會(huì)發(fā)生分裂。0Ri1，R i代表輪盤賭游戲中的存活概率。當(dāng)能量上升Ri0時(shí)，若1/ Ri 1，粒子會(huì)分裂。若01/ Ri 1，R i代表輪盤賭游戲中的存活概率。如果粒子能量下降到Ei以下，抽樣進(jìn)行指定的分裂或賭。如果在一個(gè)粒子軌跡中經(jīng)過(guò)多于一個(gè)能量分界，用Ri得到的值決定輸出。如果能量下降到Ei以下，一定會(huì)發(fā)生特定的分裂或賭。如果能量上升到Ei以上，通常要進(jìn)行相反的游戲。例如，假定在lev處以 50％的概率指定了輪盤賭，如果粒子能量升到lev以上時(shí)，它又將被一分為二。中子能量上升可能是因?yàn)榱炎兓蛘呤菬嶂凶拥南蛏仙⑸洹Ｔ谟械那闆r下，在能量上升時(shí)，可能并不希望做相反的分裂或輪盤賭游戲（尤其在裂變占優(yōu)勢(shì)的問(wèn)題中）。為此設(shè)置了一種選擇。當(dāng)分裂參數(shù)Ri給成負(fù)值時(shí)，只對(duì)能量下降時(shí)做分裂或輪盤賭游戲，而在能量上升時(shí)，則不做相反的游戲。例如卡片： esplt:n 2 .1 2 .01 .25 .001，粒子將被一分為二；，將又一次一分為二；，將以25％的殘存概率進(jìn)行輪盤賭。有權(quán)窗如果指定了權(quán)窗，能量的分裂或賭就由權(quán)窗獨(dú)自完成。這時(shí)的Ri是權(quán)窗的能量修正。如果能量下降到Ei以下，權(quán)窗除以Ri劃分。相反的，如果能量上升到Ei以上，權(quán)窗乘以Ri。如果經(jīng)過(guò)了多個(gè)能量分界，權(quán)窗的修正由Ri得到的值決定。例如卡片：ESPLT:N 2 .1 2 .01 .25 .001WWP:N 5 3 5 0 0 0 J J，權(quán)窗將被2除；，將又一次被2除；。如果WWP卡的第八項(xiàng)輸入是1（缺省值為0），除了上一段的