.. _[t-volume]:

[ T-Volume ] セクション
==================================================

各タリーにおいて **mesh=reg** とした場合、特定の **unit** の結果を正しく求めるためには、各セルの体積が必要となります。
**[t-volume]** セクションは各セルの体積を自動で計算するためのタリーです。
この体積自動計算機能では、各セルの物質を真空 (void) とし、特別な線源領域を設定したモンテカルロ積分を実行します。

本機能を使用する際は、**[t-volume]** セクションを設定し、**[parameters]** セクションにおいて **icntl=14** としてください。
本タリーは、1 つのインプットファイルにおいて 1 つのみ指定できます。

**[t-volume]** の入力書式は、以下の通りです。

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

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

   * - value
     - explanation
   * - **reg**
     - メッシュ型です。
       **reg** のみ指定可能です。

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

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

   * - value
     - explanation
   * - ``file name``
     - ファイル名を定義します。

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

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

   * - value
     - explanation
   * - (optional, D=**file**)
     - 再開始計算時の過去タリーファイル名です。

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

.. list-table:: **r-out**
   :header-rows: 0

   * - value
     - explanation
   * - (optional, D=0.0)
     - 外部ボイドの半径 [cm] です。

.. include:: ./commontally/title.rst
.. rst-class:: no-caption-number

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

   * - value
     - explanation
   * - (optional, D=0)
     - モンテカルロ積分に関するオプションです。
   * - **0**
     - 粒子軌跡による積分です。
   * - **1**
     - 線源発生点による積分です。

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

.. list-table:: **s-type**
   :header-rows: 0

   * - value
     - explanation
   * - **1, 2**
     - 体積計算の線源の種類です。
   * - **1**
     - 球体系線源です。
       原点 **(x0, y0, z0)** と半径 **r0** が必要です。
   * - **2**
     - 直方体線源です。
       **x0, x1, y0, y1, z0, z1** が必要です。

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

.. list-table:: **x0, y0, z0**
   :header-rows: 0

   * - value
     - explanation
   * -
     - 線源の形状を定義する座標です。

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

.. list-table:: **x1, y1, z1**
   :header-rows: 0

   * - value
     - explanation
   * -
     - 線源の形状を定義する座標です。

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

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

   * - value
     - explanation
   * -
     - 球体系線源の半径です。

.. include:: ./commontally/stdcut.rst

**s-type** により、体積計算時の線源の種類を指定できます。
**s-type=1** では、中心を **(x0,y0,z0)** とし半径を **r0** とする球の表面から、法線内向きに :math:`\cos` 分布となるよう線源を発生させます。
**s-type=2** では、**x=x0** と **x=x1**、**y=y0** と **y=y1**、**z=z0** と **z=z1** の面により形作られる直方体の 6 面から、法線内向きに一様な線源を発生させます。
いずれの種類の場合でも、対象とするセルが線源領域に十分な余裕をもって囲まれるようにしてください。
また、これらの領域に外部ボイドである物質番号 -1 の領域が含まれないように注意してください。

本機能を用いて PHITS を実行した場合、**file(6)** において出力される **[source]** には、体積計算で使用した線源の情報が書き出されます。
インプットファイルにある **[source]** の内容は出力されませんので注意してください。

.. _ex-icntl14:

.. code-block:: text
   :caption: **[t-volume]** の例題

        1:   [ T - V o l u m e ]
        2:     mesh =  reg            # mesh type is region-wise
        3:      reg = 101  102  103  104  105
        4:     file = volume.out      # file name of output for [volume]
        5:   s-type =    1            # 1: Sphere source, 2: Rectangular source
        6:       x0 =  0.0            # (D=0.0) x of sphere center
        7:       y0 =  0.0            # (D=0.0) y of sphere center
        8:       z0 =  0.0            # (D=0.0) z of sphere center
        9:       r0 =  50.0           # radius of sphere

**[t-volume]** において **mesh=reg** とし、体積を求めたいセルの番号を **reg=** にて指定します。
**s-type** の選択とこれに関連した **x0** などのパラメータを設定することで、体積計算用の特別な線源領域を定義します。
基本的に、対象とするセルを囲むように線源領域を決定してください。
対象とするセルが線源領域を越えていると正しく計算できません。
また、逆に線源領域を広く取り過ぎると、計算時間が長くなってしまいますので注意してください。
求めた体積は、**file=** で指定したファイルに出力されます。
この出力ファイルの形式は以下の通りです。

.. code-block:: text
   :caption: **[t-volume]** の出力形式

   [ T - V o l u m e ] off
        mesh =  reg            # mesh type is region-wise
             ....   ....     ....   ....
             ....   ....     ....   ....

   [ V o l u m e ]
     non      reg      vol      non
       1      101   5.0370E+02  0.2909
       2      102   2.4727E+03  0.1634
             ....   ....     ....   ....
             ....   ....     ....   ....

   # Information for Restart Calculation
             ....   ....     ....   ....
             ....   ....     ....   ....

**[volume]** セクションの形式で出力されるので、**infl** コマンド等により、**icntl=0** とした本計算で使用できます。
ここで、**[volume]** の最後の列にある数字は計算した体積の統計誤差で、相対誤差です。
再開始計算のための情報も出力されていますので、統計誤差が大きい場合は、**[parameters]** に **istdev=-1 or -2** を加えて、再開始計算を行ってください。