READ(iocins,'(2i6)') nnzab, nnzaeA typical line is
1 99125In this case the target ranges from a neutron ZA = 1 to a fission fragment ZA = 99125. It is very common to use this option because the standard first line of the `mcfgen.input' file permits the specification of only a single target.
The remaining options are usually not used, and they are given here
only for the sake of completeness. These options all depend on the
values of flags as described above in the discussion of the first line
of the `mcfgen.input' file.
If = 1 or = 3, the next line of the `mcfgen.input' file gives the range of excitation levels in the format:
READ(iocins,'(2e11.4)') elvb, elve
If = 2 or = 3, the next line of the `mcfgen.input' file gives the range of temperatures in the format:
READ(iocins,'(2e11.4)') tempb, tempe
If = 0, we get the energy bin boundaries from the `mcfgen.input' file. But we first read the number of bin boundaries.
READ (iocins, '(i4)') nbounds READ (iocins, '(6e12.5)') (bounds(i), i = 1, nbounds)
If = 0, we get the flux from the `mcfgen.input' file. We first read the order of the Legendre expansion lmax, which must be zero. We then read nflux, the number of (energy, flux) entries. Finally, we read the (energy, flux) pairs.
READ (iocins, '(i4)') lmax READ (iocins, '(i4)') nflux READ (iocins, '(6e12.5)') (eflux(k), k = 1, nflux)
If = 0, we get the temperatures from the `mcfgen.input' file.
READ (iocins, '(i4, 68x, i1)') numtemp Read (iocins, '(6e12.5)') (temture(i), i = 1, numtemp)Otherwise, mcfgen reads this information from the `bdfls' file. Note again that numtemp must be 1 because Doppler broadening is not implemented here.