Dear Fluka experts,
I want to build a source with the following specifications:
It should read the spectrum from the input file.
Particles emitted should be isotropic.
The source should start within a region of the geometry.
I tried to edit the source_newgen.f subroutine with the help of FLUKA forum. Unfortunately, I ran into errors. The input file, source_newgen.f, particle spectrum file and the errors are uploaded for your kind help. n.txt (3.3 KB) input_source.flair (1.7 KB) input_source.inp (2.4 KB) error_output.txt (24.8 KB) source_newgen.f (21.8 KB)
You need to put the regionsource subroutine at the very end of the source routine file (or in a separate file), then you need to have a call to it from the source routine.
You need to provide the three coordinate variables, the number of the region, and the minimum and maximum values of the coordinates which encompasses the region in question (so it is not calculated in the regionsource routine).
Dear @horvathd,
Thank you for your kind reply. I have now corrected the source.f code as suggested by you and now it compiles successfully. When I now run the code together with my input file, I could see that the particles are starting from the origin (0,0,0), not from the region of my interest (here region no 3). Can you please help me where the things are getting wrong. I am attaching the files once again along with the USRBIN plot for your reference. n.txt (3.3 KB) source_newgen.f (22.1 KB) input_source.flair (2.0 KB) input_source.inp (2.3 KB)
Moreover, there are few more lines in the original code which I have not added. I am wrting them here. Please let me know if these lines are useful for my code as well.
The following line makes the starting region search much more
robust if particles are starting very close to a boundary:
CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
& NRGFLK(NPFLKA), IDISC )
Do not change these cards:
CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
NLATTC (NPFLKA) = MLATTC
CMPATH (NPFLKA) = ZERZER
CALL SOEVSV
RETURN
Dear @horvathd,
Thank you once again for your kind help in correcting the errors in the FORTRAN file. Just to make sure that the particles are starting within region 3 of the geometry, I changed the particle to alpha, which is having very small range and will be killed very near to its geometry. Moreover, I tried to print the starting coordinates of the source particle. I still find that there are many particles which are starting within the inside RPP (here region 4), whereas I want the particles to start within the region wall (here region 3).
I am enclosing the files once again. Please let me know if I am doing something wrong. source_newgen_regionsource.f (22.0 KB) input_source001.out (101.2 KB) input_source.inp (2.3 KB) input_source.flair (2.6 KB)
With regards,
Sujoy
if you limit the projection of the plotting not to include the front or back wall, or you use the Geoviewer to visualize the URSBIN result, you will be able to see that the particles ar indeed only starting from the wall.
Dear @horvathd ,
Thank you for explaining this. Now I have got the desired data plots which shows the particle starts only from region 3 of the geometry. The plots are attached here:
