FLUKA run with EMF option

From: Khnkanos N. Sanosyan (khnko@asls.candle.am)
Date: Fri Jul 29 2005 - 12:39:17 CEST

  • Next message: Alfredo Ferrari: "Re: FLUKA run with EMF option"

    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/
    

  • Next message: Alfredo Ferrari: "Re: FLUKA run with EMF option"

    This archive was generated by hypermail 2.1.6 : Fri Jul 29 2005 - 13:40:31 CEST