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 - 17:36:42 CEST