Card 9: Calculation region details for tet mesh geometries
----------------------------------------------------------

.. _`param:imat-tet`:

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

.. list-table:: **IMAT** :phits-same:`PHITS: same name`
   :header-rows: 0

   * - Value
     - Description
   * -
     - The number of distinct materials present in the problem.

.. _`param:mtscore-tet`:

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

.. list-table:: **MTSCORE** :phits-same:`PHITS: same name`
   :header-rows: 0

   * - Value
     - Description
   * - (D\ :sub:`P`\ =0)
     - Select which materials will have activation assessed in them.
   * - = 0
     - Evaluate activation and decay in all materials.
   * - = N
     - Only evaluate activation and decay in the N\ :math:`^{th}` material.

**Card 9a: Material information listings**

This card must be repeated ``IMAT`` times; each entry is on a separate line.

.. _`param:mt-list-tet`:

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

.. list-table:: **MT-LIST** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Number of nuclides present in this material.

See :ref:`Card 9b <c9b>` from the region-wise entry. The nuclide composition of each material is specified by a Card 9b entry starting on the line immediately following each ``MT-LIST`` entry and should be ``MT-LIST`` lines long. :math:`\qquad` (Note that ``MT-LIST`` behaves like ``ITGNCLS`` for region-wise geometries.)

**Card 9c: Mesh information listings**

.. _`param:imesh-tet`:

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

.. list-table:: **IMESH** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Total number of mesh regions (number of tetrahedrons).

Following the ``IMESH`` entry line should be a single line containing the column header names listed below followed by ``IMESH`` lines with values in each column for each mesh region. Below are listed the column headers with descriptions for what information should be provided in the respective columns in each line for each mesh region.

.. _`param:non-tet`:

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

.. list-table:: **NON** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Index of tetrahedron region relative to all tetrahedron regions.

.. _`param:tetra`:

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

.. list-table:: **TETRA** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Tetrahedron number from ``*.ele`` elements file.

.. _`param:mt`:

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

.. list-table:: **MT** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Index of material (as ordered when listed in Card 9a) of the tetrahedron. (For instance, a ``MT`` of ``2`` refers to the second material listed in Card 9a.)

.. _`param:fluxs-tet`:

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

.. list-table:: **FLUXS** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0.0, D\ :sub:`P`\ = automatically calculated [#]_, requires ``JMODE``\ =1,2)
     - Total neutron flux below 20 MeV in [n/(cm\ :math:`^2`\ :math:`\cdot`\ sec)].

.. _`param:xcm`:

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

.. list-table:: **XCM** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0.0, D\ :sub:`P`\ = automatic)
     - x-coordinate of tetrahedron center of mass [cm].

.. _`param:ycm`:

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

.. list-table:: **YCM** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0.0, D\ :sub:`P`\ = automatic)
     - y-coordinate of tetrahedron center of mass [cm].

.. _`param:zcm`:

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

.. list-table:: **ZCM** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0.0, D\ :sub:`P`\ = automatic)
     - z-coordinate of tetrahedron center of mass [cm].

.. _`param:volume-tet`:

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

.. list-table:: **VOLUME** :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=1.0, D\ :sub:`P`\ =automatic)
     - Total volume of this mesh region in [cm\ :math:`^3`].

**Card 9d: Special output options for tetrahedral meshes**

.. _`param:foamout`:

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

.. list-table:: **FOAMOUT** :phits-same:`PHITS: same name`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0)
     - Select which additional special tet mesh files will be produced for each output time step ``ITOUT``.
   * - = 0
     - Do not output any additional files.
   * - = 1
     - Output OpenFOAM-formatted files (``foamout1_t``\ :math:`^*`\ ``.foam``).
   * - = 2
     - Output CSV-formatted files (``foamout2_t``\ :math:`^*`\ ``.csv``).
   * - = 3
     - Output both OpenFOAM and CSV-formatted files.

.. _`param:foamvals`:

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

.. list-table:: **FOAMVALS** :phits-same:`PHITS: same name`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0, requires ``FOAMOUT``\ ``=1`` or ``FOAMOUT``\ ``=3``)
     - Select what quantity from each tetrahedron will be written to the OpenFOAM-formatted files generated at each output time step ``ITOUT``.
   * - = 0,1
     - Total decay heat [W]
   * - = 2
     - Gamma decay heat [W]
   * - = 3
     - Beta decay heat [W]
   * - = 4
     - Alpha decay heat [W]
   * - = 5
     - Photon effective dose rate [:math:`\mu`\ Sv\ :math:`\cdot`\ m\ :math:`^2`/h] (calculated using the photon fluence-to-dose coefficients specified by ``IDOSECF``)
   * - = 6
     - Activity [Bq]
   * - = -N
     - Input any of the above values multiplied by -1 to have results normalized to volume (ex. [W/cm\ :math:`^3`])

.. [#]
   PHITS calculates ``FLUXS`` as ``amp`` (in source/sec, from [T-Dchain] input parameters) times the region’s total neutron flux (in n/(source\ \ :math:`\cdot`\ \ cm\ \ :math:`^2`) as determined by the automatically generated [T-Track] tally).