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

.. _`param:imat-xyz`:

.. 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-xyz`:

.. 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-xyz`:

.. 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>` on the previous page. 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-xyz`:

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

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

   * - Value
     - Description
   * -
     - Total number of mesh regions (nx\ :math:`\times`\ ny\ :math:`\times`\ nz). [#]_

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-xyz`:

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

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

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

.. _`param:ix`:

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

.. list-table:: **IX** :phits-other:`PHITS: other means`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Index of mesh region in the x-axis mesh.

.. _`param:iy`:

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

.. list-table:: **IY** :phits-other:`PHITS: other means`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Index of mesh region in the y-axis mesh.

.. _`param:iz`:

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

.. list-table:: **IZ** :phits-other:`PHITS: other means`
   :header-rows: 0

   * - Value
     - Description
   * -
     - Index of mesh region in the z-axis mesh.

.. _`param:fluxs-xyz`:

.. 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:volume-xyz`:

.. 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`].

.. _`param:volmat`:

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

.. list-table:: **VOLMAT**\ ``NNN`` :phits-not-settable:`PHITS: not settable`
   :header-rows: 0

   * - Value
     - Description
   * - (D=1.0, D\ :sub:`P`\ =automatic)
     - This column should be repeated ``IMAT`` times with ``NNN`` starting as ``001`` and incrementing for each material. The value listed here is the volume in [cm\ :math:`^3`] of this mesh region which is filled with each material; the sum of these ``VOLMAT``\ ``NNN`` values should equal ``VOLUME``, minus whatever portion of the mesh region, if any, is in a void region.

**Card 9d: 2D output plotting options for xyz mesh geometries**

Requires all x, y, and z axes to have mesh type 2, 3, 4, or 5 (evenly spaced meshes; breaks if any meshes are type 1). In [T-Dchain], these variables are called ``IPLTMODE`` and ``IPLTAXIS`` (without the “o”).

.. _`param:iplotmode`:

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

.. list-table:: **IPLOTMODE** :phits-other:`PHITS: other means`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0)
     - Select which axes combinations will have 2D plots made.

.. _`param:iplotaxis`:

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

.. list-table:: **IPLOTAXIS** :phits-other:`PHITS: other means`
   :header-rows: 0

   * - Value
     - Description
   * - (D=0)
     - Select axis units. If ``IPLOTMODE`` is specified, this defaults to 1 instead of 0.
   * - = 0
     - Grid index number.
   * - = 1
     - Distance in [cm].

.. [#]
   Note: large meshes in [T-Dchain] can require high amounts of memory in PHITS (:math:`\sim`\ \ 4 MB per voxel).

.. [#]
   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).