Thank you very much for the help Mikhail! Your modified source file did indeed work and I'm finally able to run simulations without an issue. Thank you for the prompt reply and help!
Regards,
Michel Georges Najarian
________________________________
From: Mikhail Polkovnikov <Михаил Полковников> <pmk_at_ihep.ru>
Sent: Tuesday, April 9, 2019 10:25 AM
To: Michel Georges Najarian; fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: Custom source.f file error while filtering by particle type
Dear Michel,
I have checked your source file, and there problems are (in my opinion):
1. You must sample the number N from 1 to NNN, but in your file N=N+1.
2. Statistical weight is 0 for negative muons (IJ=11) and some other particles, the simulation just has crushed on my computer.
You can try a changed source file from attachment.
Best regards,
Mikhail
On 08.04.2019 18:09, Michel Georges Najarian wrote:
Dear Fluka experts,
(TLDR at the bottom)
I'm currently using the two step method to simulate GCR hitting a tepc detector situated at 86m above sea level, in Ottawa. Currently I have been successful in the first step, using a custom mgdraw.f file to save the results in a temporary file and then, using a small piece of fortran code that I wrote, manipulate the data to my specifications. In the end, the source.f file reads the information from tempfileDOWN which contains the readings for all particle types. the second argument that I'm passing to the SOURCE subroutine through the input file, is the particle type that I would like to score or, 201 if I want to score everything. When selecting 201 as the argument, source.f works without an issue and I have results. However, when selecting a particle type, for example 11 which is the negative muons, the program fails. I have added multiple write statements to look at the log and see where the issue happens but every time, the issue happens after source.f finishes running. I am by no means an expert in fluka thus I am unsure on how to proceed. I have been trying to fix this issue for multiple days now but have failed.
Additionally, I used the rdsource.solution.f file from the 4th advanced course and workshop as my template and modified it. The main piece of code that I have added is that after reading the value of each required element (particle type, x, y, z coordinates, etc), I verify if the particle type of the element is the same as the type of element we want. If not, then go the the next line of the document ( in source.f: GOTO 310). If this is indeed the issue, I do not see why. any help will be greatly appreciated. I have attached the input file, the tempfile as well as my source.f file. Thank you very much.
TLDR: When trying to filter the particles in the source.f file, it fails but when selecting all types of particles, the program completes succesfully. Not Sure where the error is.
Regards,
Michel Georges Najarian
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Apr 09 2019 - 18:32:20 CEST