From: Khnkanos N. Sanosyan (khnko@asls.candle.am)
Date: Thu Aug 04 2005 - 10:26:12 CEST
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
*
.
.
.
-- 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 : Thu Aug 04 2005 - 10:55:55 CEST