Hello,
I strongly suggest using the source_newgen.f file instead of the source.f one.
For example, if you want a beam with 50-50 Co-60 and Ba-133:
double precision, save :: choose_isotope
choose_isotope = FLRNDM (xdummy)
radioactive_isotope = .true.
if (choose_isotope < 0.5) then
heavyion_atomic_number = 27
heavyion_mass_number = 60
else
heavyion_atomic_number = 56
heavyion_mass_number = 133
end if
The photon spectrum will look like this:
I have also attached the source file.
source_newgen.f (18.3 KB)
Have a nice day,
Mihaela