Dear Arghya,
what about choosing randomly among the 5 sphere centers?
Let us say with a code like:
REAL*4 CENTERS(5), X0
DATA CENTERS/-10.,-5.,0.,5.,10./
X0 = CENTERS(FLOOR(6*FLRNDM(XDUMMY)))
where X0 is the center of your current sphere and then:
CALL RACO(TXX, TYY, TZZ)
RADIUS = R*FLRNDM(XDUMMY)
XFLK (NPFLKA) = X0 + RADIUS*TXX
YFLK (NPFLKA) = RADIUS*TYY
ZFLK (NPFLKA) = RADIUS*TZZ
You have to generate 5 times histories, but the result should be
the one you are looking for. I hope this helps.
Kind regards,
Andrea
Il 15/09/2017 14:15, arghyac_at_barc.gov.in ha scritto:
> Dear Experts,
>
> I need to simulate 5 spherical sources with different spatial
> co-ordinate in a single source.f card.
>
> I model one sphere (centre at (0,0,0) in source.f card as follows
>
> *particle corodinate
> CALL RACO(TXX, TYY, TZZ)
> RADIUS = R*FLRNDM(XDUMMY)
> XFLK (NPFLKA) = RADIUS*TXX
> YFLK (NPFLKA) = RADIUS*TYY
> ZFLK (NPFLKA) = RADIUS*TZZ
>
> I guess I need to edit XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA) to
> sample 5 spherical sources in single subroutine.
> Co-ordinate of other 4 sphere's centre are (-10,0,0), (-5,0,0),(5,0,0)
> and (10,0,0). radius of each = 1.5 cm.
>
> Now how to sample other spheres here or how to edit it to accommodate
> all the 5 sources at the same time.
>
> Kindly give some suggestions.
>
> Thanking you in advance.
>
> With regards,
> Arghya
>
> __________________________________________________________________________
>
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id=acc_info
>
--
========================================================================
Dr. Andrea Fontana tel: +39 0382 987991
Istituto Nazionale fax: +39 0382 423241
di Fisica Nucleare
Sezione di Pavia e-mail: andrea.fontana_at_pv.infn.it
Via Bassi 6 web : www.pv.infn.it/~fontana
27100 PAVIA, Italy
========================================================================
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Sep 15 2017 - 17:47:45 CEST