From: Jason Crawford (jaycrawf@triumf.ca)
Date: Wed Oct 17 2007 - 23:44:28 CEST
Hello, I am having troubles using the LOW-NEUT/PART-THR/LOW-MAT
cards.When I disable these cards my input runs fine, but when enabled,
these cards are responsible for an output message, which reads:
**** Fluka stopped in Lowset: undefined atomic density for the low en.
xsec mat. corresponding to the Fluka comp. mat. n. 10 ****
I have looked through past messages on the fluka discussion board to
find out that this may be caused by my LOW-MAT card for Nitrogen,
suggesting it should read:
LOW-MAT NITROGEN 1.0 -2.0 ...
(referencing the new code for Nitrogen, -2.0)
After making sure my :LOW-MAT card for Nitrogen was correct my
simulation still gave me the same error. Can this be explained? Can
Someone help me find the solution?
I thank you in advance for any help you can give me.
Jason Crawford
I have enclosed my input file, but be warned it is a working copy, and I
use the FLAIR interface.
My input file:
****************************
TITLE
TNF Neutron Channel Spectrum
GEOBEGIN 0.0 0.0 0.0
COMBNAME
0 0 TNF NEUTRON CHANNEL GEOMETRY
*
* X-axis up, Y-axis horizontal, Z-axis along beam
* Origin at centre of beam dump
*
* ========= Body Definitions ==============================================
*
* A large box for the blackhole
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
RPP B1 -999999.9 999999.9 -999999.9 999999.9 -999999.9 999999.9
*
* A cylinder for the target + proton beam channel and the water tank
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
ZCC B2 0.0 0.0 10.0
XCC B3 0.0 0.0 36.5
*
* Top and bottom of shielding
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
YZP B4 100.0
YZP B5 -137.16
*
* Top and bottom of neutron channel
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
YZP B6 -26.67
YZP B7 -34.29
*
* Other horizontal planes
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
YZP B8 45.72
YZP B9 0.0
YZP B10 -30.48
YZP B11 -45.72
YZP B12 -52.3875
*
* Cylinders for inner and outer shielding
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
XCC B13 0.0 0.0 260.0
XCC B14 0.0 0.0 400.0
*
* Target boundaries
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
XYP B15 25.0
XYP B16 -25.0
*
* Planes perpendicular to neutron channel
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
PLA B17 0.0 0.8660254 0.5 0.0 0.0 0.0
PLA B18 0.0 0.8660254 0.5 0.0 225.16661 130.0
PLA B19 0.0 0.8660254 0.5 0.0 226.26646 130.635
PLA B20 0.0 0.8660254 0.5 0.0 277.95951 160.48
PLA B21 0.0 0.8660254 0.5 0.0 311.76915 180.0
PLA B22 0.0 0.8660254 0.5 0.0 330.75242 190.96
*
* Planes parallel to neutron channel
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
PLA B23 0.0 -0.5 0.8660254 0.0 -4.445 7.698966
PLA B24 0.0 -0.5 0.8660254 0.0 4.445 -7.698966
PLA B25 0.0 -0.5 0.8660254 0.0 -7.62 13.19823
PLA B26 0.0 -0.5 0.8660254 0.0 7.62 -13.19823
PLA B27 0.0 -0.5 0.8660254 0.0 -15.24 26.39645
PLA B28 0.0 -0.5 0.8660254 0.0 15.24 -26.39645
*
* An additional cylinder for biasing
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
XCC B29 0.0 0.0 150.0
END
*
* ========= Region Definitions ===========================================
*
* Blackhole
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
BH1 5 | +B1 -B4 | +B1 +B5 | +B1 +B4 -B5 -B14
*
* Proton Channel
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
PC2 5 +B2 +B3 +B16
PC3 5 +B2 +B29 -B3 +B17
PC4 5 +B2 +B13 -B29 +B17
PC5 5 +B2 +B14 -B13 +B17
*
* Target
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
TA6 5 +B2 +B15 -B16
*
* Water tank
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
WT7 5 | +B3 -B2 +B4 -B12 -B17 | +B3 +B2 -B15
WT8 5 +B3 -B2 +B4 -B12 +B17
*
* Neutron Channel Part 1
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
NC9 5 +B6 -B7 +B23 -B24 -B17 +B29 -B3
N10 5 +B6 -B7 +B23 -B24 -B17 -B29 +B18
*
* Steel Box in Neutron Channel
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
S11 5 +B6 -B7 +B23 -B24 +B19 -B18
*
* Neutron Channel Part 2
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
N12 5 +B6 -B7 +B23 -B24 +B20 -B19
*
* Steel Box - 1ft by 1ft with neutron channel in the middle
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
S13 5 | +B20 -B18 +B9 -B6 +B25 -B26 | +B20 -B18 +B7
-B10 +B25 -B26 | +B20 -B18 +B6 -B7 +B25 -B23
| +B20 -B18 +B6 -B7 +B24 -B26
*
* Neutron Channel Part 3
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
N14 5 +B6 -B7 +B23 -B24 +B21 -B20
*
* Detector Box
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
D15 5 +B6 -B7 +B23 -B24 +B22 -B21
*
* Air Box
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
A16 5 | +B22 -B20 +B8 -B6 +B27 -B28 | +B22 -B20 +B7
-B11 +B27 -B28 | +B22 -B20 +B6 -B7 +B27 -B23
| +B22 -B20 +B6 -B7 +B24 -B28
*
* Concrete Boxes
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
C17 5 | +B20 -B18 +B8 -B9 +B27 -B28 | +B20 -B18 +B10
-B5 +B27 -B28 | +B20 -B18 +B9 -B10 +B27 -B25
| +B20 -B18 +B9 -B10 +B26 -B28
C18 5 +B22 -B20 +B27 -B28 +B11 -B5
*
* Rest of Concrete
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
C19 5 | +B14 -B13 +B8 -B5 +B17 -B2 | +B14 -B13 +B8
-B5 -B17 -B27 | +B14 -B13 +B8 -B5 -B17 +B28
| +B14 -B22 +B8 -B5 +B27 -B28 | +B18 -B13 -B17
+B8 -B5 +B27 -B23 | +B18 +B14 -B13 -B17 +B8
-B5 +B24 -B28 | +B18 +B14 -B13 +B8 -B6 +B23
-B24 | +B18 +B14 -B13 +B7 -B5 +B23 -B24
*
* Top Steel
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
S20 5 +B14 -B3 +B4 -B8
S21 5 +B13 -B29 +B8 -B6 -B2 +B17
S22 5 +B13 -B29 +B8 -B6 -B17
S23 5 +B29 -B3 +B8 -B6 -B2 +B17
S24 5 +B29 -B3 +B8 -B6 -B17
*
* Bottom Steel
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
S25 5 +B13 -B29 +B12 -B5 +B17
S26 5 +B13 -B29 +B12 -B5 -B17
S27 5 +B29 -B3 +B12 -B5 +B17
S28 5 +B29 -B3 +B12 -B5 -B17
S29 5 +B3 +B12 -B5 +B17
S30 5 +B3 +B12 -B5 -B17
S31 5 +B13 -B29 +B7 -B12 +B17
S32 5 +B13 -B29 +B7 -B12 -B17
S33 5 +B29 -B3 +B7 -B12 +B17
S34 5 +B29 -B3 +B7 -B12 -B17
*
* Rest of Steel
* *REG.....OR....1OR....2OR....3OR....4OR....5OR....6OR....7OR....8OR....9
S35 5 +B13 -B29 +B6 -B7 +B17
S36 5 +B29 -B3 +B6 -B7 +B17
S37 5 +B13 -B29 +B6 -B7 -B17 -B23
S38 5 +B29 -B3 +B6 -B7 -B17 -B23
S39 5 +B13 -B29 +B6 -B7 -B17 +B24
S40 5 +B29 -B3 +B6 -B7 -B17 +B24
END
#if 0
GEOEND 200.0 400.0 400.0 -200.0 -400.0 -400.0DEBUG
GEOEND 200.0 400.0 400.0 &
#endif
GEOEND
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* MATERIAL DEFINITIONS ******************************************
MATERIAL 1.0 1.0079 8.9D-05 3.0
1.0HYDROGEN
MATERIAL 7.0 14.0067 0.0012606 4.0
0.0NITROGEN
MATERIAL 8.0 15.9994 0.001429 5.0 0.0OXYGEN
MATERIAL 13.0 26.982 2.6989 6.0
0.0ALUMINUM
MATERIAL 14.0 28.09 2.33 7.0
0.0SILICON
MATERIAL 20.0 40.08 1.55 8.0
0.0CALCIUM
MATERIAL 26.0 55.847 7.8658 9.0 0.0IRON
MATERIAL 0.0 0.0 1.0 10.0 0.0WATER
COMPOUND 2.0 HYDROGEN 1.0 OXYGEN 0.0 WATER
MATERIAL 0.001205 11.0 AIR
COMPOUND 4.195D-05 NITROGEN 1.128D-05 OXYGEN 0.0 AIR
MAT-PROP 85.7 AIR
STERNHEI 10.5961 1.7418 4.2759 0.10914 3.3994 AIR
MATERIAL 0.0 0.0 2.35 12.0
0.0CONCRETE
COMPOUND 0.02015 SILICON 0.04587 OXYGEN 0.001743
ALUMINUMCONCRETE
COMPOUND 0.01375 HYDROGEN 0.00266 CALCIUM 0.0
CONCRETE
MATERIAL 0.0 0.0 2.654 13.0 0.0TARGET
COMPOUND 0.08 SILICON 0.92 OXYGEN 0.0 TARGET
#if 0
MATERIAL 6.0 12.011 1.8 CARBON
#endif
#if 0
MATERIAL 12.0 24.305 1.738
MAGNESIU
#endif
#if 0
* Concrete
* Concrete has a wide variation in density and composition. The above
* description is for poured structural concrete with 10% moisture
* content. Concrete block will have a density of about 2.05 g/cm3.
* Ranges of concrete composition are : C (8-25%), O (38-60%), Si (8-18%).
* Concrete composition can be analyzed cheaply by commercial laboratories.
MATERIAL 2.34
Concrete
#endif
#if 0
COMPOUND 23.0 CARBON 40.0 OXYGEN 12.0
SILICONConcrete
COMPOUND 12.0 CALCIUM 10.0 HYDROGEN 2.0
MAGNESIUConcrete
#endif
*
* ========= Region Materials ===============================================
ASSIGNMA BLCKHOLE BH1
BLACKHOLE
ASSIGNMA VACUUM PC2 PC5 1.0 VACUUM
ASSIGNMA TARGET TA6 TARGET
ASSIGNMA WATER WT7 WT8 1.0 WATER
ASSIGNMA VACUUM NC9 N10 1.0 VACUUM
ASSIGNMA IRON S11 S13 2.0 IRON
#if 0
LOW-MAT WATER 1.0 -2.0 293.0
HYDROGEN
#endif
ASSIGNMA AIR N12 A16 2.0 AIR
ASSIGNMA AIR D15 AIR
ASSIGNMA CONCRETE C17 C19 1.0
CONCRETE
ASSIGNMA IRON S20 S40 1.0 IRON
#if 0
PLOTGEOM 1.0 1.0 0.0 -5000.0 1.0
NOFORMAT
*Vertical Section along neutron channel
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
#endif
#if 0
* Old card tag was: 0
error
#endif
LOW-MAT BLCKHOLE 0.0 0.0 0.0
LOW-MAT WATER 0.0 0.0 0.0
LOW-MAT HYDROGEN 1.0 -2.0 293.0
HYDROGEN
LOW-MAT NITROGEN 1.0 -2.0 293.0
HYDROGEN
LOW-MAT OXYGEN 8.0 16.0 293.0 OXYGEN
LOW-MAT ALUMINUM 13.0 -2.0 293.0
LOW-MAT SILICON 14.0 -2.0 293.0
LOW-MAT IRON 26.0 -2.0 293.0
LOW-MAT CALCIUM 20.0 -2.0 293.0
#if 0
LOW-MAT MAGNESIU 25.0 55.0 293.0
MANGANES
#endif
#if 0
LOW-MAT CARBON 6.0 -2.0 293.0 CARBON
#endif
#if 0
STOP
#endif
#if 0
PLOTGEOM 1.0 1.0 0.0 -5000.0 1.0 FORMAT
*Vertical Section along proton beam line
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
#endif
#if 0
* Old card tag was: 0
error
#endif
#if 0
STOP
#endif
#if 0
PLOTGEOM 1.0 1.0 0.0 -5000.0 1.0
NOFORMAT
*Horizontal Section along proton beam line
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
#endif
#if 0
* Old card tag was: 0
error
#endif
#if 0
STOP
#endif
#if 0
PLOTGEOM 1.0 1.0 0.0 -5000.0 1.0
NOFORMAT
*Horizontal Section along neutron channel
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 0.0
#endif
#if 0
* Old card tag was: 0
error
#endif
#if 0
STOP
#endif
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* BEAM SPECIFICATION *********************************************
*
* This is a 400 MeV Proton beam travelling down the z-axis.
BEAM -0.37 0.0 0.0 -2.0 -2.0 1.0PROTON
BEAMPOS 0.0 0.0 -350.0 0.0 0.0
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* EVENT GENERATOR, PARTICLES AND FIELDS **************************
*
* Switch on the generator with default level densities
EVENTYPE 0.0 0.0 1.0 0.0 0.0 0.0EVAP
*
* Discard neutrinos, antineutrinos, pi-zeros, electrons, photons, and muons
DISCARD NEUTRIE ANEUTRIE PIZERO ELECTRON
DISCARD POSITRON PHOTON MUON+ MUON-
*
* Turn off electromagnetic cascade
EMF
EMF-OFF
*
* Set threshold to 19.6 MeV for neutrons
PART-THR -0.0196 NEUTRON
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* LOW ENERGY NEUTRONS ********************************************
*
* Switch on low-energy neutron transport
LOW-NEUT 72.0 22.0 0.0196 0.0
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* OUTPUT AND SCORING *********************************************
*
* ========= Biasing ========================================================
BIASING 3.0 1.0 0.0001 PC5 PRINT
BIASING 3.0 1.0 0.001 PC4 PRINT
BIASING 3.0 1.0 0.001 S21 S25 4.0PRINT
BIASING 3.0 1.0 0.001 S31 S35 4.0PRINT
BIASING 3.0 1.0 0.01 PC3 PRINT
BIASING 3.0 1.0 0.01 S23 S27 4.0PRINT
BIASING 3.0 1.0 0.01 S33 S36 3.0PRINT
BIASING 3.0 1.0 0.1 PC2 WT8 6.0PRINT
BIASING 3.0 1.0 0.1 S29 PRINT
BIASING 3.0 1.0 1.0 TA6 PRINT
BIASING 3.0 1.0 4.6 WT7 S30 23.0PRINT
BIASING 3.0 1.0 21.5 NC9 PRINT
BIASING 3.0 1.0 21.5 S24 S28 4.0PRINT
BIASING 3.0 1.0 21.5 S34 PRINT
BIASING 3.0 1.0 21.5 S38 S40 2.0PRINT
BIASING 3.0 1.0 100.0 N10 PRINT
BIASING 3.0 1.0 100.0 C17 C18 1.0PRINT
BIASING 3.0 1.0 100.0 S22 S26 4.0PRINT
BIASING 3.0 1.0 100.0 S32 PRINT
BIASING 3.0 1.0 100.0 S37 S39 2.0PRINT
BIASING 3.0 1.0 460.0 S11 S13 1.0PRINT
BIASING 3.0 1.0 2150.0 N14 A16 2.0PRINT
BIASING 3.0 1.0 10000.0 D15 PRINT
*
* Score energy and all particles
SCORE ENERGY ALL-PART
#if 0
* Ignore this card (for test only)
USRBDX NEUTRON -94.0 WT7 NC9 136.86
USRBDX 0.371 1D-07 200.0 &
#endif
*
* Score neutron fluence and energy
USRBIN 10.0 NEUTRON -77.0 -11.0 400.0 400.0YZ
USRBIN -19.0 -400.0 -400.0 1.0 200.0 200.0&
USRBIN 10.0 NEUTRON -78.0 200.0 50.0 400.0XZ
USRBIN -200.0 -50.0 -400.0 100.0 1.0 200.0&
#if 0
USRBDX 99.0 NEUTRON -72.0 WT7 NC9
136.86DETECT1
USRBDX 0.37 0.0196 10.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -72.0 WT7 NC9
136.86DETECT1
USRBDX 0.0196 1.2277D-4 42.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -73.0 NC9 N10
135.563DETECT2
USRBDX 0.37 0.0196 10.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -73.0 NC9 N10
135.563DETECT2
USRBDX 0.0196 1.2277D-4 42.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -74.0 N10 S11
135.4836DETECT3
USRBDX 0.37 0.0196 10.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -74.0 N10 S11
135.4836DETECT3
USRBDX 0.0196 1.2277D-4 42.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -75.0 N12 N14
135.4836DETECT4
USRBDX 0.37 0.0196 10.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -75.0 N12 N14
135.4836DETECT4
USRBDX 0.0196 1.2277D-4 42.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -76.0 N14 D15
135.4836DETECT5
USRBDX 0.37 0.0196 10.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX 99.0 NEUTRON -76.0 N14 D15
135.4836DETECT5
USRBDX 0.0196 1.2277D-4 42.0 0.0 0.0 0.0&
#endif
#if 0
USRBDX -1.0 NEUTRON 62.0 NC9 N10
135.563DETECT2
USRBDX 0.4 0.0196 1.0 0.25 0.0 100.0&
#endif
*
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
* ********* START SIMULATION ***********************************************
RANDOMIZ 1.0
* Start the Calculation
START 10000.0
STOP
This archive was generated by hypermail 2.1.6 : Thu Oct 18 2007 - 10:08:23 CEST