.. _sec-dumpall:

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

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

(1). NCOL

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

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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 :math:`\ge` 4の時です。

(2). NOCAS, NOBCH, RCASC, RSOUIN

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

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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 は、電子の :math:`dE/dx`  です。

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

     * -  | QS
       -  | :
       -  | :math:`dE/dx`  of electron at (x,y,z)

(4). IBLZ1, IBLZ2, ILEV1, ILEV2

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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 :math:`>` 0 の時に

  ::

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

  で行われます

 (b). ILAT2

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

  ::

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

  で行われます

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

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

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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)

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

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

(7). WT, U, V, W

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

    * -  | 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

 これらの変数の意味は、

 .. rst-class:: no-caption-number

 .. list-table::
    :header-rows: 0

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

(11). NZST

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

(12). NCLSTS

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

 (a). MATHZ, MATHN, JCOLL, KCOLL

  これらの変数の意味は、

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

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

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

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

     * -  | 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

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

     * -  | 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

  これらの変数の意味は、

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

     * -  | ICLUSTS
       -  | 
       -  | kind of particle
     * -  | 0
       -  | :
       -  | nucleus
     * -  | 1
       -  | :
       -  | proton
     * -  | 2
       -  | :
       -  | neutron
     * -  | 3
       -  | :
       -  | pion
     * -  | 4
       -  | :
       -  | photon
     * -  | 5
       -  | :
       -  | kaon
     * -  | 6
       -  | :
       -  | muon
     * -  | 7
       -  | :
       -  | others

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

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

  .. rst-class:: no-caption-number

  .. list-table::
     :header-rows: 0

     * -  | 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
       -  | :
       -  | :math:`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)