READ(iocins,'(2i6)') nnzab, nnzae
A typical line is
1 99125
In 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.