Hi Qi,
Here is a mgdraw routine with :
do 20 ip = 1, NP
IF (KPART(ip) .EQ. 8) THEN
WRITE(IODRAW,*) XSCO, YSCO, ZSCO,
& TKI(ip), CXR(ip), CXR(ip), CXR(ip), WEI(IP)
Thus, for only neutron secondaries (id = 8), you record position of reaction, kinetic energy of created neutron, direction and weight.
But you can change what you want. If you want (p,n) reactions only, change the IF with IF (JTRACK .EQ. 1 .AND. KPART(ip) .EQ. 8)...
Hope this helps.
Be carreful with neutrons and multigroup theory!
Julien
De: "Ding, Qi" <q.ding_at_fz-juelich.de>
À: "fluka-discuss" <fluka-discuss_at_fluka.org>
Envoyé: Lundi 4 Janvier 2021 18:34:59
Objet: [fluka-discuss]: (p.n) reaction position-MGDRAW
Dear FLUKA experts,
My simulation is a proton beam incident on a metal target(there exist water channels inside target).
I want to get the information about the position of (p,n) reaction, and I have read manual and the forum.
it seems that I need a customized mgdraw (usdraw entry) user routine.
I have never used mgdraw.f.
for this purpose, could you please help me have a look if the code is correctly modified?
in the attachment you can find my input file and mgdraw.f.
any help will be appreciated.
Best wishes and a happy new year.
Qi Ding
E-mail: q.ding_at_fz-juelich.de
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich
Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzender des Aufsichtsrats: MinDir Volker Rieke
Geschaeftsfuehrung: Prof. Dr.-Ing. Wolfgang Marquardt (Vorsitzender),
Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Jan 05 2021 - 11:44:17 CET