.. _app-comprehensive-example-dchain-input:

:ref:`4.12節 <comprehensive-example>` でPHITSが生成するDCHAIN入力ファイル
----------------------------------------------------------------------------------------

::

       htitle = 250MeV protons incident on Be-encased W/H2O/Fe column
    
   ! --- control parameters ---
        imode =       2
        jmode =       2
    
   ! --- calculation parameters ---
        idivs =      50
       ichain =     100
       itdecs =       1
       itdecn =       1
       isomtr =       2
       ifisyd =       0
       ifisye =       0
      ixsrall =       1
    
   ! --- data library parameters ---
      inxslib =     100
      idcylib =       5
   !   hnxslib = h-JP-US-EU                                                                      
   !   hdcylib = J15-EB8-ES                                                                      
    
   ! --- output parameters ---
        iyild =       2
        iggrp =       3
       ibetap =       1
        acmin =  1.0000E-20
       istabl =       0
       igsdef =       1
       igsorg =       1
      iwrtchn =       1
      chrlvth = -1.0000E+00
      iwrchdt =       0
      iwrchss =       0
     iwrchnuc =       0
      idosecf =       1
     idosunit =     333
      imtcard =       1
      imtcnum =       0
     imtcmeta =       2
      thmatnd =  1.0000E-06
    
   ! --- Proton beam current and neutron flux ---
          amp =  1.0000E-04
        ebeam =  3.0000E+00
       prodnp =  1.0000E+00
    
   ! --- irradiation time ---
       itstep =       3
         6.0000E+00 m  1.0000E+00
         4.0000E+00 m  5.0000E-01
         5.0000E+01 m  0.0000E+00
    
   ! --- output time ---
        itout =       8
         1.0000E+00 m
         3.0000E+00 m
         6.0000E+00 m
         1.0000E+01 m
         1.5000E+01 m
         2.0000E+01 m
         3.0000E+01 m
         4.0000E+01 m
    
   ! --- irradiation condition ---
       inmtcf =       1
   !  ffnmtcf =    1.0
     itdchout =       1
      hhnmtcf = spallation.dyld
     iertdcho =       1
      henmtcf = spallation_err.dyld
     iredufmt =       1
      irdonce =       1
    
   ! --- calculation region details ---
    
       iregon =      4
   !1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
      DUMMY001        1    5  2.3797E+10   spallation.dtrk  1.5080E+02
      W-180     7.5882E-05
      W-182     1.6756E-02
      W-183     9.0514E-03
      W-184     1.9376E-02
      W-186     1.7975E-02
   !1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
      DUMMY002        2    3  7.3851E+09   spallation.dtrk  5.0265E+02
      H-1       6.6880E-02
      H-2       1.0033E-05
      O-16      3.3446E-02
   !1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
      DUMMY003        3    4  9.4207E+08   spallation.dtrk  2.5133E+02
     Fe-54      4.9642E-03
     Fe-56      7.7928E-02
     Fe-57      1.7997E-03
     Fe-58      2.3951E-04
   !1)HRGCMM 2)IREGS 3)ITGNCLS 4)FLUXS 5)HNFLUXS 6)VOLUMES
      DUMMY004        4    1  1.6216E+09   spallation.dtrk  3.7196E+03
     Be-9       1.2352E-01
    
   ! --- output option for ANGEL ---
     phitsout =  1
     angelout =  1
     iangbpwr =  1
   !   angelout_act        = 0
   !   angelout_heat       = 0
   !   angelout_heat_alpha = 0
   !   angelout_heat_beta  = 0
   !   angelout_heat_gamma = 0
   !   angelout_dose       = 0
   !   angelout_region     = all, sum
   !   angelout_nuclides   = 0
    
   end