.. _standalone-example:

例題：PHITSを使用しないDCHAIN単体計算
--------------------------------------

DCHAINは主としてPHITSの[T-Dchain]と組み合わせて使用するが、PHITSの[T-Track]・[T-Yield]出力を必要としない単純な崩壊計算であれば、DCHAIN単体でも実行できる。この場合は、DCHAINの既定値からいくつかのパラメータを変更し、DCHAINデータフォルダーのパスを1行で記載した\ ``dch_link.dat``\ を計算フォルダーに置く必要がある。既定の内容を以下に示す。

::

   c:/phits/dchain-sp/data/

以下は、\ :math:`^{183}`\ Hf、\ :math:`^{160m}`\ Tm、\ :math:`^{169}`\ Ta、\ :math:`^{164}`\ Hfをそれぞれ初期放射能1 MBqとし、10日間の崩壊を複数の出力時刻で評価する例である。単体崩壊計算で通常必要となる主要なパラメータだけを示しており、必要に応じて他の出力パラメータも追加できる。変更時に重要な値は入力内のコメントで説明している。

原則として、\ ``jmode``\ ``= -1``、\ ``inmtcf``\ ``= 0``、\ ``hnfluxs``\ ``= no``\ の設定が必要である。ただし、\ ``jmode``\ ``= -1``\ を指定すれば、DCHAINが他の2つを強制的に上書きするため、実際には\ ``jmode``\ だけでも単体崩壊計算を実行できる。

::

   htitle = Standalone DCHAIN decay calculation example

     imode =       2        ! output inventory, decay heats, and photon spectra
     jmode =      -1        ! decay only, no external source
   idcylib =       5        ! default decay library using data from JENDL and ENDF

     acmin =  1.0000E-20    ! threshold activity [Bq] to appear in *.act output file
    istabl =       0        ! toggle output of stable nuclides
   iwrtchn =       1        ! decay chain scheme settings for *.dcs file
   chrlvth =  1.0000E-04    !   threshold for determining if chains printed to *.dcs file
   iwrchdt =       1        !   add extra information for contribution by reaction

   ! --- irradiation time ---
    itstep =       1        ! a decay calculation is just a single 'cooling' time step
      10.0   d  0.0         ! this time step must be >= to the final requested output time

   ! --- output time ---
     itout =        9       ! request times when inventory information is outputted
        1  m
        5  m
        10 m
        30 m
        3  h
        12 h
        1  d
        5  d
        10 d

   ! --- irradiation condition ---
    inmtcf =       0       ! set to 0 so DCHAIN won't look for [T-Yield] file

   ! --- calculation region details ---
    iregon =      1
   !1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
   DUMMY001  -1  4  0.0  no  1.0E+0 ! 4 nuclides to be tracked; set HNFLUXS to 'no' so
   Hf-183       b1.00E+06           !  DCHAIN won't look for [T-Track] file; vol = 1 cm^3.
   Tm-160m      b1.00E+06           ! Inventory is normally in 1E24 atoms/cm^3, but putting
   Ta-169       b1.00E+06           ! 'b' before the value changes it to an activity in Bq.
   Hf-164       b1.00E+06           !  Here each nuclide has a starting activity of 1 MBq.

   end