[en] Double differential distributions of neutrons produced by 100, 150, 200 and 250 MeV protons stopped in a thick iron target were calculated with the FLUKA Monte Carlo code at four emission angles: forward, 45 deg., transverse and 135 deg. backwards. The attenuation in thick iron shields of the dose equivalent due to neutrons, protons, photons and electrons was also calculated. The contribution to the total ambient dose equivalent from photons and protons is limited to a few percent at maximum. Source terms and attenuation lengths are given as a function of energy and emission angle, along with fits to the Monte Carlo data, for shallow depth and deep penetration in the shield. A brief discussion of simulations performed with composite iron/concrete shields is also given, showing the need for further investigations