Multiple cylindrical source as shown in figure

Thank you for quick reply.

Can I code like below:


random = FLRNDM(XDUMMY)

if (random .lt. beam1_probablility) then
call regionsource ( coordinate_x, coordinate_y, coordinate_z, [region number_1], [xmin], [xmax],[ymin], [ymax], [zmin], [zmax] )
else if (random .lt. beam1_probability + beam2_probability) then
call regionsource ( coordinate_x, coordinate_y, coordinate_z, [region number_2], [xmin], [xmax], [ymin], [ymax], [zmin], [zmax] )
else if (random .lt. beam1_probability + beam2_probability + beam3_probability) then
call regionsource ( coordinate_x, coordinate_y, coordinate_z, [region number_3], [xmin], [xmax], [ymin], [ymax], [zmin], [zmax] )
else
call regionsource ( coordinate_x, coordinate_y, coordinate_z, [region number_4], [xmin], [xmax], [ymin], [ymax], [zmin], [zmax] )
end if


Kindly correct me if the above code is wrong.

Whether you have some example source.f file using roto-translating the coordinates for such application?

Kindly give me some suggestions.