Hi Jenny,
Sorry I tried, but I don’t understand your explanation, nor what image you are referring at and what is the relationship of that image with the energy distribution.
In any case, I suggest that you add print statements in the source.f card so that you have some way to know what you are doing.
For example:
write(lunout,*)’This is my printout ‘,ZFLK(NPFLKA)
(Then you can parse “This is my printout” in the output file, *.out).
From that point on, it is up to you to know what type of physics you want to implement.
-M
On Mar 13, 2019, at 8:18 PM, JENNY MARGOT RAMOS LAZARO <jenny.ramos_at_unmsm.edu.pe<mailto:jenny.ramos_at_unmsm.edu.pe>> wrote:
Dear Santana ,Mario
es, it seems that I am wrong in the order, I am changing the initial energy distribution and I made this change by placing the already calculated values
of the parameters of my distribution so that I did not have to initialize them and inserted this distribution .
* Kinetic energy of the particle (GeV).
PMOFLK (NPFLKA) = PBEAM
TKEFLK (NPFLKA) = (1.D0 / (PBEAM ** 2)) * 0.25
But when I draw the images I get a straight line like the previous one, maybe There is anything incorrect in the distribution.In the example, using JBEAM,
is it important to use it ? Do I have to add more values to my distribution such as AM or JBEAM?
Here the example of the manual:
PMOFLK (NPFLKA) = PBEAM
Therefore, the kinetic energy (in GeV) must be derived:
TKEFLK (NPFLKA) = SQRT (PBEAM ** 2 + AM (IJBEAM) ** 2) - AM (IJBEAM)
(where AM is the rest mass, in COMMON PAPROP, and IJBEAM is the
type of particle, in COMMON BEAMCM)
Thanks for the help.
Best Regards
regards
El mié., 13 mar. 2019 a las 21:08, Santana, Mario (<msantana_at_slac.stanford.edu<mailto:msantana_at_slac.stanford.edu>>) escribió:
Where do you initialize the values for Z2 and Z1? It seems like you assign Z1 to Whasou(1) and Z2 to Whasou(2), while I believe that you wanted to do the reverse assignment, i.e. assign Whasou(1) to Z1, ….
Right?
On Mar 13, 2019, at 1:02 PM, JENNY MARGOT RAMOS LAZARO <jenny.ramos_at_unmsm.edu.pe<mailto:jenny.ramos_at_unmsm.edu.pe>> wrote:
Dear Santana, Mario
Thank you, I managed to compile as you recommended, but when plotting my results in FLAIR, my graphics are shown in Y in 0 as seen in the file that I mention and maybe the reason is in something that I modified in source.f.
In the source file modify the following:
*** User initialization ***
WHASOU(1)=Z1
WHASOU(2)=Z2
END IF
* Particle coordinates
XFLK (NPFLKA) = XBEAM
YFLK (NPFLKA) = YBEAM
ZFLK (NPFLKA) = 10.D0 + (Z2-Z1)*FLRNDM(XXX)
To place a one-report distribution along the z-axis, I followed the instructions in the manual.
I wanted to start making this distribution as a test, but my goal is to place an initial distribution of energy, power law type, which leads me to ask, the final part that I am modifying in FUENTE.f It is the correct one to implement. this type of distribution? .
Thanks for the help.
Best Regards
regards
El mar., 12 mar. 2019 a las 21:07, Santana, Mario (<msantana_at_slac.stanford.edu<mailto:msantana_at_slac.stanford.edu>>) escribió:
How did you compile the source routine? It seems like FLUKA is told that there is a source routine through the SOURCE card in the input file, but you haven’t one.
Try this:
$FLUPRO/flutil/lfluka -m fluka -o ./flukahps source.f
In the directory where you have the ‘source.f’ routine (I assume that it is also the directory where you have the input file?).
Then re-run using flukahps as your new executable.
-M
On Mar 12, 2019, at 4:54 PM, JENNY MARGOT RAMOS LAZARO <jenny.ramos_at_unmsm.edu.pe<mailto:jenny.ramos_at_unmsm.edu.pe>> wrote:
Dear FLUKA experts:
I'm using the source.f file, following the instructions in the manual to insert a uniform distribution along the z axis, but when I run the script, the .out file says that source.f does not exist or is not executable.
These are the lines that I have added in the source.f file
* Particle coordinates
XFLK (NPFLKA) = XBEAM
YFLK (NPFLKA) = YBEAM
Z1 = 10.D0
Z2 = 80.D0
ZFLK (NPFLKA) = 10.D0 + (Z2-Z1)*FLRNDM(XXX)
Inside input.inp I did not declare any (WHAT (1,2,3 ..)) in the source card because I declared Z1 and Z2 in source.f.
I attach the .inp, .out and source.f files.
Thanks for the help.
Best Regards
<exampl.out><exampl.inp><source.f>
<source.f><exampl_48_plot.jpg><exampl.inp>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Mar 14 2019 - 06:05:49 CET