Dear @horvathd,
Thank you very much for your reply and all your help! Maybe this is worth another post but how would you do an SOBP in the “source-newgen.f” subroutine? I have been trying to write several files, each with its own histogram as suggested in point 3 of Region of particle cannot be determined in extreme position - #3 by horvathd, but I cannot make it work, any suggestion on how to proceed?
On another note, is it possible to use both “source.f” subroutines: the old and new ones, where the old one would be in charge of sampling the SOBP and the new one in charge of sampling the circular distribution?
Best wishes,
Daniel