Hi Stuart
I modified and compiled the source file like you said but I got error
which was something along the lines of "floating point error". Also, if i
keep the original line in the source file , I get a reasonable spectrum
(please see attached. X axis is in eV).
saurabh
On Mon, Jun 25, 2018 at 8:01 PM, Stuart Ansell <stuart.ansell_at_maxiv.lu.se>
wrote:
> Dear Saurabuh,
>
>
> Not sure that it will fix anything but you have a typo error in your
> source1.f file (I think). You write:
>
> ESAMPLE=XYZ*(ENEPRO(I)-ENEPRO(I-1))*(ENEPRO(I)+ENEPRO(I-1))
> ESAMPLE=SQRT(ESAMPLE+ENEPRO(I-1)**2)-ENEPRO(I-1)
> ESAMPLE=ESAMPLE/(ENEPRO(I)-ENEPRO(I-1))
>
>
> attempting to get a linear energy based on XYZ being between 0.0 and 1.0.
>
>
> and I believe that you intended that middle line should be
>
> ESAMPLE=SQRT(ESAMPLE**2+ENEPRO(I-1)**2)-ENEPRO(I-1)
>
>
> It leaves rather large chunks of the energy spectrum unsampled as it
> currently stands.
>
>
>
> Regards,
>
>
> Stuart Ansell
>
>
>
>
>
>
> ------------------------------
> *From:* owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it>
> on behalf of Saurabh Mukherjee <sxm3816_at_gmail.com>
> *Sent:* Sunday, June 24, 2018 7:43:29 PM
> *To:* fluka-discuss
> *Subject:* [fluka-discuss]: thermal neutron induced radioactivity in
> natural tungsten
>
> Dear experts,
> I am trying to simulate residual activity due to thermal neutron incident
> on natural tungsten. The input file is attached, along with the source and
> spectrum file. The isotopes generated by FLUKA are - W-187, W-185 and
> Re-187.
> But when I compare it to online calculator (https://www.ncnr.nist.gov/
> resources/activation/) (the inputs to the online calculator are shown in
> attachment 2 and the results in attachment3), the results are different. I
> see that in online calculator, the W-187, W-185 isotopes are present but
> there are some extra isotopes present as well-W-181, W-188 and Ta-183.
> These are not present in FLUKA results. Is there some problem with my input
> file? Can some one please take a look?
> sincerely
> saurabh
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Jun 25 2018 - 18:31:15 CEST