Direction Cosines are 0 via mgdraw.f

Advanced Features and User Routines
1

Dear FLUKA experts,

I’ve written an mgdraw.f routine meant to print out the direction cosines for neutrons found in a particular region of interest. While I’m able to print neutron energies, event positions, etc, the direction cosines of each neutron are printed as 0. Is there any indication as to what the problem could be? Here is the code in the mgdraw.f file:

IF (MREG .EQ. 3 .AND. JTRACK .EQ. 8 .AND. LLOUSE .EQ. 0) THEN ! Neutron in the ROIK'
            WRITE(99,*) NCASE, ETRACK, LTRACK, WTRACK,
     &       XSCO, YSCO, ZSCO, CXTRACK, CYTRACK, CZTRACK
            ! Tagging the particle for killing by setting LLOUSE = 1
            LLOUSE = 1
      END IF

Thank you for your help!

2

It’s just due to a typo, the concerned variables are CXTRCK, CYTRCK, CZTRCK, without any A in their name.

3

Thank you Francesco.

I am humbled.

1 Like