Dear Fluka experts,
Regarding my previous question on the "not complete" 3D data exporting
via usbrea, I thought to share with you an example to clarify my problem
in case that my previous question wasn't clear!
Attached, you will find a very simple FLUKA code (usbrea.inp) that
scours Edep in a small Cu cylinder, with binning: x and y = -2.5 to 2.5
(1 cm width, i.e. 5 bins for each axis) and z = 0 to 5 (1 cm width. So,
5 bins for this axis as well). The runs were merged in "usbrea_98.bnn"
via flair. $FLUPRO/flutil/usbrea has been used to produce and export the
formatted Edep data values in 3D matrix and saved in
"usbrea_3D_DataExport.txt". A 2D plot (in x-y plan) at fixed z=0 was
generated, and fluka saved the results automatically in
"usbrea_z0_98_plot.dat" file.
- The correct number of binnings is:
x-axis: -2.5, -1.5, -0.5, 0.5, 1.5, 2.5 (6 points)
y-axis: -2.5, -1.5, -0.5, 0.5, 1.5, 2.5 (6 points)
z-axis: 0, 1, 2, 3 ,4 ,5 (6 points)
Now,
1- usbrea_3D_DataExport.txt gives/export 125 Edep data points. These
come from 5 (x-bins) * 5 (y-bins) * 5 (z-bins).
2- usbrea_z0_98_plot.dat gives/exports 36 Edep data points. These come
from 6 (x-bins) * 6 (y-bins) * 1 (for the first z bin value, z=0).
The two export methods are not compatible! The 3D data export method via
usbrea gives less data points than it must? Why is that? Do I use usbrea
in a wrong way?
It is critical for me to understand the exporting method of usbrea,
because I further use the exported 3D data points from fluka as
input/initial values in other program.
I really hope to get some feedback from the fluka experts on this issue!
With my best regards,
Tamer
On 13.03.2020 00:24, Tamer Tolba wrote:
>
> Dear all,
>
> I am scouring energy using USRBIN card, according to the following
> steps (copied from 3D data file exported with usbrea):
>
> Cartesian binning n. 1 "Edep_Core " , generalized particle n. 208
> X coordinate: from -2.0000E+02 to 2.0000E+02 cm, 200 bins (
> 2.0000E+00 cm wide)
> Y coordinate: from -2.0000E+02 to 2.0000E+02 cm, 200 bins (
> 2.0000E+00 cm wide)
> Z coordinate: from 3.0000E+03 to 3.3200E+03 cm, 160 bins (
> 2.0000E+00 cm wide)
> Data follow in a matrix A(ix,iy,iz), format (1(5x,1p,10(1x,e11.4)))
>
> I realized that for 3D data exporting via usbrea, I got 6,400,000 data
> points, which theoretically comes from 200*200*160 bins! However, this
> is not correct because from x (and y) coordinates that starts at x(y)=
> -200 to x(y)=200 (with 2 cm wide) I should have 201 bins not 200 bins!
> The same applies for the z coordinate I should have, for z=3000 to
> z=3320 (2 cm wide), 161 bins not 160. So, in total I should have
> 201*201*161= 6,504,561 data points!
>
> Since the 2D data exportation (via 2D *plot.dat) is the only way that
> fuka exports "coordinates+data_values", I decided to use it to check
> its results! For a big surprise I found that at fixed z (e.g. z=3000)
> I got, from the *plot.dat file 40,401 data points for the same x(y)
> range, i.e. x(y)= -200 to x(y)= 200 (with 2 cm wide), which indicates
> that it follows the correct binning of 201*201 (not 200*200).
>
> Why the 3D exportation produces less data points than the actual
> binning? Or I miss something here?
>
> Regards,
>
> Tamer
>
> --
> Dr. Tamer Tolba
> Forschungszentrum Jülich
> Institut für Kernphysik (IKP)
> 52425 Jülich, Germany
> Tel.: +49 2461 61 6860
> Fax.: +49 2461 61 96415
>
>
> ------------------------------------------------------------------------------------------------
> ------------------------------------------------------------------------------------------------
> 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,
> Prof. Dr. Sebastian M. Schmidt
> ------------------------------------------------------------------------------------------------
> ------------------------------------------------------------------------------------------------
>
--
Dr. Tamer Tolba
Forschungszentrum Jülich
Institut für Kernphysik (IKP)
52425 Jülich, Germany
Tel.: +49 2461 61 6860
Fax.: +49 2461 61 96415
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sun Mar 15 2020 - 20:44:48 CET