On 04.11.2016 19:03, adwasil wrote:
> Dear FLUKA experts.
> I've defined the source Cosines(tx, ty, tz) in the spherical coordinates in the file source.f. In the case of such a definition factor rNorm = 1.
> I've attached cross sections of the neutron flux in the XY plane (the result seems to be OK here) and in the XZ plane (I would expect a distribution similar to the XY one but the distribution I've obtained here is significantly different). What am I doing wrong?
> Best regards
> Adam Wasilewski
>
Dear Adam,
It looks like kind of isotropic distribution. The isotropic distribution
in spherical coordinates has two sample parameters: phi - azimuth angle,
cosine theta - cosine of polar angle. Sampling of the polar angle is
wrong in your case.
Try to use this part in your code:
***********************************************************************
************ Changes **********
***********************************************************************
C THETA - polar angle ( CTHETA - cosine THETA, STHETA - sine THETA)
C PHI - azimuthal angle
PHI = TWOTWO * PIPIPI * FLRNDM( XDUMMY )
CTHETA = TWOTWO * FLRNDM( XDUMMY ) - ONEONE
STHETA = SQRT( ONEONE - CTHETA**2 )
TXFLK (NPFLKA) = STHETA * COS(PHI)
TYFLK (NPFLKA) = STHETA * SIN(PHI)
TZFLK (NPFLKA) = CTHETA
***********************************************************************
Best regards,
Mikhail
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sat Nov 05 2016 - 11:30:38 CET