5.2.26. Dumpall オプション

dumpall=1, -1 を指定して計算すると、再計算に必要な輸送の情報を全てファイルに書き出します。このファイルを用いてicntl=12 を指定すると、同じdumpallの値で計算した時にダンプしたデータを基に再計算します。輸送に関しては、最初に計算した全ての情報を再現します。再計算の時には、最初に計算したインプットファイルが必要になります。タリーを除いて同じものを使わないと矛盾が生じます。maxcas, maxbch は、ダンプファイルに指定されているものを使いますので変更できません。最初に計算したときと異なるタリーを入れて再計算したいときなど有効です。ただし、dumpallの機能を利用するときは、巨大なデータファイルが作成されることがありますから注意してください。

バイナリ形式で書き出したデータは、マシン間での互換性がありませんので注意してください。以下に書き出されるデータの種類と順序を説明します。

(1). NCOL

NCOLは、プログラムの中の変数で、粒子輸送のタイミングの認識変数です。

NCOL


1
:
start of calculation
2
:
end of calculation
3
:
end of a batch
4
:
source
5
:
detection of geometry error
6
:
recovery of geometry error
7
:
termination by geometry error
8
:
termination by weight cut-off
9
:
termination by time cut-off
10
:
geometry boundary crossing
11
:
termination by energy cut-off
12
:
termination by escape or leakage
13
:
(n,x) reaction
14
:
(n,n'x) reaction
15
:
sequential transport only for tally
16
:
surface cross for WW of xyz mesh

NCOL=1, 2, 3 の時は、書き出しはNCOLだけです。従って、以下の書き出しは、NCOL \(\ge\) 4の時です。

(2). NOCAS, NOBCH, RCASC, RSOUIN

これら4つの変数は、NCOL=4の時にだけ書き出されます。それぞれの意味は、

NOCAS
:
current event number in this batch
NOBCH
:
current batch number
RCASC
:
real number of NOCAS+maxcas*(NOBCH-1)
RSOUIN
:
sum of the weight of source particle

(3). NO, MAT, ITYP, KTYP, JTYP, MTYP, RTYP, OLDWT

これらの変数の意味は、

NO
:
cascade id in this event
MAT
:
material id
ITYP
:
particle type
KTYP
:
particle kf-code
JTYP
:
charge number of the particle
MTYP
:
baryon number of the particle
RTYP
:
rest mass of the particle (MeV)
OLDWT
:
weight of the particle at (x,y,z)

(a). QS

この変数は、ITYP=12, 13 の時、すなわち電子、陽電子の時のみ出力されます。QS は、電子の \(dE/dx\) です。

QS
:
\(dE/dx\) of electron at (x,y,z)

(4). IBLZ1, IBLZ2, ILEV1, ILEV2

これらの変数の意味は、

IBLZ1
:
cell id at (x,y,z)
IBLZ2
:
cell id after crossing
ILEV1
:
level structure id of the cell at (x,y,z)
ILEV2
:
level structure id of the cell after crossing

(a). ILAT1

階層構造をもつCELLの変数です。書き出しは、ILEV1 \(>\) 0 の時に

write(io) ( ( ILAT1(i,j), i=1,5 ), j=1,ILEV1 )

で行われます

(b). ILAT2

階層構造をもつCELLの変数です。書き出しは、ILEV2 \(>\) 0 の時に

write(io) ( ( ILAT2(i,j), i=1,5 ), j=1,ILEV2 )

で行われます

(5). COSTH, UANG(1), UANG(2), UANG(3), NSURF

各変数の意味は以下の通りです。これらはバージョン2.30から出力されるようになりました。

COSTH
:
cosine of an angle of incidence in a surface crossing
UANG(1,2,3)
:
x,y,z component of a normal vector of its surface, respectively
NSURF
:
internal number of the surface
Note that this is different from the surface number defined in the [surface] section

(6). NAME, NCNT(1), NCNT(2), NCNT(3)

これらの変数の意味は、

NAME
:
collision number of the particle
NCNT(1,2,3)
:
values of counter 1, 2, and 3

(7). WT, U, V, W

これらの変数の意味は、

WT
:
weight of the particle at (xc,yc,zc)
U, V, W
:
unit vector of momentum of the particle

(8). E, T, X, Y, Z

これらの変数の意味は、

E
:
energy of the particle at (x,y,z) (MeV)
T
:
time of the particle at (x,y,z) (nsec)
X, Y, Z
:
position coordinate of the preceding event point (cm)

(9). EC, TC, XC, YC, ZC

これらの変数の意味は、

EC
:
energy of the particle at (xc,yc,zc) (MeV)
TC
:
time of the particle at (xc,yc,zc) (nsec)
XC, YC, ZC
:
position coordinate of the particle (cm)

(10). SPX, SPY, SPZ

これらの変数の意味は、

SPX, SPY, SPZ
:
unit vector of spin direction of the particle

(11). NZST

この変数は粒子がもつ電荷(charge state)を表します。

(12). NCLSTS

この変数は、Collisionが起こった時に、すなわち、NCOL=13, 14 の時にのみ書き出されます。意味は、衝突もしくは崩壊による生成粒子の個数です。原子核も含みます。NCLSTS \(>\) 0 の時に、次の変数が書き出されます。

(a). MATHZ, MATHN, JCOLL, KCOLL

これらの変数の意味は、

MATHZ
:
Z number of the mother nucleus
MATHN
:
N number of the mother nucleus
JCOLL
:
reaction type id1
KCOLL
:
reaction type id2

JCOLL, KCOLL はそれぞれ次のような意味を持ちます。

JCOLL


0
:
nothing happen
1
:
Hydrogen collisions
2
:
Particle Decays
3
:
Elastic collisions
4
:
High Energy Nuclear collisions
5
:
Heavy Ion reactions
6
:
Neutron reactions by data
7
:
Photon reactions by data
8
:
Electron reactions by data
9
:
P,d,a, and photo-nuclear reactions by data
10
:
Neutron event mode
11
:
Delta Ray production
12
:
Muon atomic interaction
13
:
Photon reactions by EGS5
14
:
Electron reactions by EGS5
15
:
Photon photonuclear interaction
16
:
Negative muon captured by nucleon
17
:
Muon photonuclear interaction
18
:
Electron recoil by track structure mode
19
:
Muon pair production (photon -> mu+ mu-)
20
:
User defined interaction
KCOLL


0
:
normal
1
:
high energy fission
2
:
high energy absorption
3
:
low energy n elastic
4
:
low energy n non-elastic
5
:
low energy n fission
6
:
low energy n absorption

(b). ICLUSTS, JCLUSTS, QCLUSTS, JCOUNT

これらの変数は、各生成粒子の情報を記述するもので、NCLSTSだけ繰り返し、また、それぞれの変数も以下のように配列を持ちます。書き出しは以下のように行われます。

do i = 1, NCLSTS
write(io) ICLUSTS(i)
write(io) ( JCLUSTS(j,i), j=0,8)
write(io) ( QCLUSTS(j,i), j=0,12)
write(io) ( JCOUNT(j,i),  j=1,3)
end do

これらの変数の意味は、

ICLUSTS

kind of particle
0
:
nucleus
1
:
proton
2
:
neutron
3
:
pion
4
:
photon
5
:
kaon
6
:
muon
7
:
others
JCLUSTS(i)


i = 0
:
angular momentum
= 1
:
proton number
= 2
:
neutron number
= 3
:
ityp
= 4
:
status of the particle 0: real, \(<\) 0 : dead
= 5
:
charge number
= 6
:
baryon number
= 7
:
kf code
= 8
:
isomer level (0: ground state, 1: first isomer, 2: second isomer)
QCLUSTS(i)


i = 0
:
impact parameter
= 1
:
x-component of unit vector of momentum
= 2
:
y-component of unit vector of momentum
= 3
:
z-component of unit vector of momentum
= 4
:
\(e_{tot} = \sqrt{ p^2 + m^2 }\) (GeV)
= 5
:
rest mass (GeV)
= 6
:
excitation energy (MeV)
= 7
:
kinetic energy (MeV)
= 8
:
weight
= 9
:
time (nsec)
= 10
:
x coordinate (cm)
= 11
:
y coordinate (cm)
= 12
:
z coordinate (cm)