Advance Multiple- complex Radioactive source

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:
image

I have also attached the source file.
source_newgen.f (18.3 KB)

Have a nice day,
Mihaela

4 Likes