From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Thu Aug 04 2005 - 17:33:54 CEST
Have you tried the geometry debugger? This is the first thing
to do when starting a new job, even before geometry problems
show up.
In any case, if you wish people to look at your input, please send it
as an attachment. Some mailing systems chop lines at the 70th column,
making it very hard to read, and impossible to pass to the debugger.
Alberto
On Thu, 4 Aug 2005, Khnkanos N. Sanosyan wrote:
>
> Dear Fluka users,
>
> The global task is the same as in my previous message. in this time the
> run is stopped with core dump and with the following messages as output
> diagnostics:
>
> 1. in .err file:
> Region of particle in position 0.000000000E+00 0.000000000E+00
> -3.000000000E+00
> cannot be determined, error count: 1
> try again to establish the current region moving the particle of a
> 1.000000000E-06 long step.
> 2. in .log file:
> ENTER IN HIST INIT.
> HROPEN: ISTAT= 0
> , procedure georeg.f/georeg.
> Attempt to access the 3-th element of variable jomerr.
> 3. in .out file:
> 0****EXIT BEING CALLED FROM LOOKZ .
> ...
> Please, what are my mistakes in geometry establishment?
> 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 in advance.
> Best regards,
> Khnkanos
> .
> .
> .
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> BEAMPOS 0.0 0.0 -3.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 1.0E-04
> COMBINAT
> *
> * Geometry title (in special format):
> *
> 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.
> *
> * In our case we are interested in cascade longitudinal and lateral
> development
> * for the 300cm Cylinder of radii from 0.0cm to 30cm,
> * with axis coincide with the z axis.
> * An infinite cylinder of radius 30cm, with axis equal to the z axis
> * could be established as:
> *
> *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:
> *
> * external vacuum at the beginning
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> 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 -3.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.
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> 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 -5
> *
> * 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 -5
> *
> *
> END
> *
> *
> * N.B.
> * Region Volumes have to be inserted, inasmuch as we have defined: IVOPT=3
> !!!
> *
> *
> * Here is the stage of booking of the volume array VNOR(I=1,5) with the
> * format (7E10.5):
> *
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 12345678
> *
> 3.9250+11 4.4000+09 2.8260+04 8.4780+05 2.8000+09
> *
> *
> GEOEND
> *
> * {MATERIAL} card: defines a single-element material or a compound
> * (coupled to a COMPOUND card)
> * Here defines Materials and Compounds of Ordinary ConCreTe(CCT)
> ** atomic #Z at.wt #A g/cm^3 material # not used mass # name
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> *
> *
> * 1) Hydrogen (H)
> MATERIAL 1.0 1.0079 8.988E-5 3.0 0.0 1.
> HYDROGEN
> * 2) Carbon (C)
> MATERIAL 6.0 12.011 2.265 6.0 0.
> CARBON
> * 3) Oxygen (O)
> MATERIAL 8.0 15.9994 0.001429 8.0 0.
> OXYGEN
> * 4) Magnesium (MG)
> MATERIAL 12.0 24.305 1.738 9.0 0.
> MAGNESIU
> * 5) Aluminum (AL)
> MATERIAL 13.0 26.982 2.70 10.0 0.
> ALUMINUM
> * 6) Iron (FE)
> MATERIAL 26.0 55.847 7.877 11.0 0. IRON
> * 7) Calcium (CA)
> MATERIAL 20.0 40.08 1.54 12.0 0.
> CALCIUM
> * 8) Sodium (NA)
> MATERIAL 11.0 22.9898 0.975 13.0 0.
> SODIUM
> * 9) Silicon (SI)
> MATERIAL 14.0 28.0855 2.33 26.0 0.
> SILICON
> * 10) Potassium (K)
> MATERIAL 19.0 39.102 0.031165 27.0 0.
> POTASSIU
> *
> * 11) For Ordinary ConCreTe compound using atom relative contents
> *
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
> MATERIAL 0.0 0.0 2.30 28.0 0.0 0.
> CCTORDIN
> COMPOUND -1.0 3.0 -0.1 6.0 -52.9 8.
> CCTORDIN
> COMPOUND -0.2 9.0 -3.4 10.0 -1.4 11.
> CCTORDIN
> COMPOUND -4.4 12.0 -1.6 13.0 -33.7 26.
> CCTORDIN
> COMPOUND -1.3 27.0
> CCTORDIN
> *
> * {ASSIGNMAt} card: defines the correspondence between region indices and
> * material indices. It defines also regions with magnetic field.
> *
> * Here assigns materials to regions
> * Ordinary ConCreaT in region 3
> *
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> * mater.# lower
> * bound of
> * the region
> * indices
> * External Black Hole
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> ASSIGNMA 1.0 1.0
> * Internal Vacuum
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> ASSIGNMA 2.0 2.0
> *
> * Internal Vacuum bounded with Z<0.0cm half-space
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> ASSIGNMA 2.0 3.0
> *
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> ASSIGNMA 28.0 4.0 4.0 1.0 -1.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
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> 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)
>
> *23456789 123456789 123456789 123456789 123456789 123456789 123456789
> 123456789
> STEPSIZE 4.0E+01 3.0 5.0 1.0E-04
> *
> .
> .
> .
>
>
-- Alberto Fassò SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025 Phone: (1 650) 926 4762 Fax: (1 650) 926 3569 fasso@slac.stanford.edu
This archive was generated by hypermail 2.1.6 : Thu Aug 04 2005 - 23:04:59 CEST