From: Khnkanos N. Sanosyan (khnko@asls.candle.am)
Date: Fri Jul 29 2005 - 12:39:17 CEST
Dear FLUKA users,
I am calculating 3 GeV e- beam induced radiation dose rate in
Ordinary Concrete by following to the electron propagation in Concrete
absorber.
I receive the following messages as output diagnostics:
1. in .err file:
Abort called from CDWPRS reason LNUMRC = .TRUE. Run stopped!
STOP LNUMRC = .TRUE.
2. in .log file:
STOP STOP: FLUKA ABORTED statement executed
3. in .out file:
Abort called from CDWPRS reason LNUMRC = .TRUE. Run stopped!
STOP LNUMRC = .TRUE.
And I can not fix the reason/meaning of such kind error diagnostic!
What kind of connection one should keep during the establishment of the
"WHAT(1-5)" parameters in {STEPSIZE} card depended on the accuracy
parameter in {GEOBEGIN} input card?
Please, I am appending to my message some lines from my task's
steering input file,
and I am asking you if you have any suggestions that might help me to
overcome my problem?
Thank you for your help and time in advance!
Sincerely yours,
Khnkanos
...
.
.
.
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
BEAMPOS 0.0 0.0 -1.0E+00
*
* Here starts the geometry descripton (see 9} Combinatorial Geometry):
* -----------------------
* a large cylinder of concrete delimited by vacuum layers at the top
* and bottom. Here depth is fixed at 3m and the radius 30cm.
* Electrons are injected in the initial vacuum.
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
* geometry the Comb.Geom
* errors are accuracy is used
* printed parameter
GEOBEGIN 0.0 0.1 COMBINAT
*
3 10 Ordinary ConCreTe(CCT)
*
* Body data:
*******************************************************************************
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
* SPHere with coordinates of the center and R radius:
SPH 001 0.0 0.0 0.0 +5000.0
* code #/# V_x V_y V_z Radius
*
*******************************************************************************
* Rectangular Parallelepiped. Code: RPP
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
RPP 002 -1000.0 +1000.0 -1000.0 +1000.0 -100.0 +1000.0
*
*******************************************************************************
* Infinite Circular Cylinder parallel to the axis z.
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
* code #/# Z_Axis_x Z_Axis_y R_Radius
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
*
ZCC 003 0.0 0.0 30.0
*
******************************************************************************
* Infinite half-space delimited by a coordinate plane. Code: XYP.
* The half-space "inside the body" is the locus of points for which z < V_z.
*
XYP 004 0.0
*
*******************************************************************************
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
*
XYP 005 300.0
*
*
******************************************************************************
*
END
*
* Region data:
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
* External vacuum at the beginning
1 5 1 -2
*
* Internal vacuum at the beginning created as a buffer vacuum
* region containing the whole geometry, which must itself be
* contained within the external (mandatory) blackhole region.
* It is then suggested that the beam impact point be chosen
* in vacuum, slightly upstream the actual one on a
* material boundary.
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
2 5 2 3 4 5 -2
*
* Internal vacuum in which the region of starting point of particle
* is established by moving the particle of a -1.0E+00 cm long step:
* it is a Vacuum Cylinder with coordinates r = 30 cm and -100 <= z < 0.0 cm
* before the Ordinary Concrete layer.
*
3 5 2 -3 4 5 -3
*
* Ordinary concrete layer: a Cylinder with r = 30 cm, and 0.0 <= z < 300
cm
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
4 5 2 3 -4 5 -4
*
* Internal vacuum at the end: a Rectangular Parallelepiped with
* -1000 <= x <= 1000 cm, -1000 <= y <= 1000 cm, 300 <= z <= 1000 cm
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
5 5 2 -3 -4
*
*
END
*
*
* N.B.
* Region Volumes have to be inserted, inasmuch as we have defined: IVOPT=3 !!!
*
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 12345678
*
3.9250+11 4.4000+09 2.8260+04 8.4780+05 2.8000+09
*
GEOEND
.
.
.
* {ASSIGNMAt} card: defines the correspondence between region indices and
* material indices. It defines also regions with magnetic field.
*
* Here is assigned the materials to regions
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
*
* External Black Hole
ASSIGNMA 1.0 1.0
* Internal Vacuum
ASSIGNMA 2.0 2.0
*
* Internal Vacuum bounded with Z<0.0cm half-space
ASSIGNMA 2.0 3.0
*
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
ASSIGNMA 28.0 4.0 5.0 1.0 -1.0 0.0
* CCTORDIN from region WHAT(2) to region WHAT(3) in steps of WHAT(4);
* WHAT(5)< 0.0 resets the default (no magnetic field)
* in the region(s) defined by WHAT(2),(3), and (4);
* WHAT(6), SDUM : not used.
*
* Internal Vacuum
ASSIGNMA 2.0 5.0
*
* {STEPSIZE} card: sets the minimum and maximum step size on a region-by-region
* basis for transport of all charged particles
* (hadrons, muons & electrons).
*23456789 123456789 123456789 123456789 123456789 123456789 123456789 123456789
* min max lowerbound upper_ in steps
* step size stp size from bound: to of 1 cm
* in cm in cm region region (default)
STEPSIZE 1.0-05 1.0-01 1.0 5.0 0.001
*
.
.
.
-- Khnkanos N. SANOSYAN ASLS (Armenian Synchrotron Light Source) CANDLE (Center for the Advancement of Atjaryan St. 31, Avan, Natural Discoveries using Light Emission) Yerevan 375040 / ARMENIA Tel.: (+374-1) 612569 Private: (+374-1) 645865 Fax: (+374-1) 629806 Mail: khnko@asls.candle.am http: //www.candle.am/
This archive was generated by hypermail 2.1.6 : Fri Jul 29 2005 - 13:40:31 CEST