Re: Joint calculation with source file

From: Ogawa Tatsuhiko <ogawa_at_n.t.u-tokyo.ac.jp>
Date: Mon, 11 Jan 2010 17:37:00 +0900

Dear FLUKA experst,

Please find the attached input file. This is the input file I used
actually. I also attached the source.f file again for convenience.

I would appreciate if you could point out what is wrong in this input.

Sincerely,
Ogawa

-- 
**************************************************
Department of Nuclear Engineering and Management
         Graduate School of Engineering
            The University of Tokyo
             Tatsuhiko Ogawa
      E-mail:  ogawa_at_n.t.u-tokyo.ac.jp
               togawa_at_cern.ch
         Tel:   090-6934-6732(Domestic call only)
               +81-3-5841-2905
Facsimile :   +81-3-3813-2010
Address: 2-11-16 Engineering Bld.12th Room421
              Yayoi Bunkyo-ku, Tokyo  113-0032
**************************************************
Stefan Roesler wrote:
> Dear Ogawa,
> 
> I have an account at SLAC and thus took the liberty to have a look at your 
> working directory
> 
> /afs/slac.stanford.edu/u/sa/komi/Himac
> 
> The source routine is working fine as far as 'score' and 'sample' are 
> concerned.
> 
> However, I am not sure which input file you are using. Could you please 
> tell me? I only see coll.inp however, without SOURCE card and heavy ion 
> primary. Thus, it cannot be this one..
> 
> Best regards
> Stefan
> 
> 
> 
> On Fri, 8 Jan 2010, Ogawa Tatsuhiko wrote:
> 
>> Dear FLUKA experts,
>>
>> Thank you very much for your suggestion, Dr. Fasso. I understood the 
>> variable "score" was used as a integer though it was defined as a real 
>> mumber.
>> But I saw the same message even after defining "score" explicitly as an 
>> integer or replacing "score" with "iscor"
>> I feel sorry for failing to understand what your suggestion meant but 
>> could you explain definition of what parameter must be modified?
>>
>> Sincerely,
>> Ogawa
>>
>>
>>
> 

TITLE
HMC CONC ACTIVATION
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
GLOBAL 9999.0 1.0
PHYSICS 3.0 EVAPORAT
PHYSICS 1.0 COALESCE
BEAM -0.40 0.0 0.0 0.0 0.0 1.0HEAVYION
HI-PROPE 6.0 12.0
BEAMPOS 0.0 0.0 -58.8 0.0 0.0 0.0
EVENTYPE 0.0 0.0 2.0 0.0 0.0 0.0DPMJET
SOURCE
*
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
* *
MATERIAL 83.0 208.98 9.78 26.0 0.0 0.0BISMUTH
MATERIAL 30.0 65.409 7.14 27.0 0.0 0.0ZINC
MATERIAL 49.0 114.818 7.31 28.0 0.0 0.0INDIUM
MATERIAL 19.0 39.0983 0.862 29.0 0.0 0.0POTASSIU
MATERIAL 16.0 32.07 1.96 30.0 0.0 0.0SULFUR
MATERIAL 25.0 54.93805 7.44 31.0 0.0 0.0MANGANES
* *
MATERIAL 0.0012048 32.0 0.0 0.0AIR
COMPOUND -.7555267 NITROGEN -0.231781 OXYGEN -0.012827 ARGONAIR
COMPOUND -0.000124 CARBON AIR
* *
MATERIAL 2.270 33.0 0.0 0.0CONC
COMPOUND -1.067E-2 HYDROGEN -0.48570 OXYGEN -2.204E-2 SODIUMCONC
COMPOUND -8.979E-3 MAGNESIU -6.02E-2 ALUMINUM -2.887E-1 SILICONCONC
COMPOUND -1.959E-2 POTASSIU -8.449E-2 CALCIUM -1.969E-2 IRONCONC
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
ASSIGNMA BLCKHOLE BLACKH
* outer vacuum
ASSIGNMA AIR AIROUT
*
* Shield 11 target
ASSIGNMA IRON TARIRON
* Au Foil
ASSIGNMA CONC CON1_1 D4_15
* *
GEOBEGIN 0.0 0.0 COMBNAME
    0 0 0.01
*
* REFERENCE FRAME
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*
* /|\ _ Z[cm]
* X[cm]| /|
* | /
* | /
* | /
* | /
* |/
* |----------> Y[cm]
*
* *******************************************************************************
* OUTER BODIES:
RCC InfCyl 0.0 0.0 -999.9 0.0 0.0 9999.0 9999.9
RCC BlaCyl 0.0 0.0 -200.0 0.0 0.0 520.0 1000.0
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* TARGET
RPP TARIRO -8.0 8.0 -8.0 8.0 0.0 6.85
RPP CONP -125.0 55.0 -60.0 60.0 56.85 146.85
XYP CONP1 62.85
XYP CONP2 68.85
XYP CONP3 74.85
XYP CONP4 80.85
XYP CONP5 86.85
XYP CONP6 92.85
XYP CONP7 98.85
XYP CONP8 104.85
XYP CONP9 110.85
XYP CONP10 116.85
XYP CONP11 122.85
XYP CONP12 128.85
XYP CONP13 134.85
XYP CONP14 140.85
*
YZP CONP0 -35.0
*
*DETECTOR«
*
RPP D0_0 -2.0 2.0 -2.0 2.0 56.851 56.852
RPP D0_6 -2.0 2.0 -2.0 2.0 62.851 62.852
RPP D0_12 -2.0 2.0 -2.0 2.0 68.851 68.852
RPP D0_18 -2.0 2.0 -2.0 2.0 74.851 74.852
RPP D0_24 -2.0 2.0 -2.0 2.0 80.851 80.852
RPP D0_30 -2.0 2.0 -2.0 2.0 86.851 86.852
RPP D0_36 -2.0 2.0 -2.0 2.0 92.851 92.852
RPP D0_42 -2.0 2.0 -2.0 2.0 98.851 98.852
RPP D0_48 -2.0 2.0 -2.0 2.0 104.851 104.852
RPP D0_54 -2.0 2.0 -2.0 2.0 110.851 110.852
RPP D0_60 -2.0 2.0 -2.0 2.0 116.851 116.852
RPP D0_66 -2.0 2.0 -2.0 2.0 122.851 122.852
RPP D0_72 -2.0 2.0 -2.0 2.0 128.851 128.852
RPP D0_78 -2.0 2.0 -2.0 2.0 134.851 134.852
RPP D0_84 -2.0 2.0 -2.0 2.0 140.851 140.852
RPP D0_90 -2.0 2.0 -2.0 2.0 146.848 146.849
*
RPP D15_0 -17.0 -13.0 -2.0 2.0 56.851 56.852
RPP D15_6 -17.0 -13.0 -2.0 2.0 62.851 62.852
RPP D15_12 -17.0 -13.0 -2.0 2.0 68.851 68.852
RPP D15_18 -17.0 -13.0 -2.0 2.0 74.851 74.852
RPP D15_24 -17.0 -13.0 -2.0 2.0 80.851 80.852
RPP D15_30 -17.0 -13.0 -2.0 2.0 86.851 86.852
RPP D15_36 -17.0 -13.0 -2.0 2.0 92.851 92.852
RPP D15_42 -17.0 -13.0 -2.0 2.0 98.851 98.852
RPP D15_48 -17.0 -13.0 -2.0 2.0 104.851 104.852
RPP D15_54 -17.0 -13.0 -2.0 2.0 110.851 110.852
RPP D15_60 -17.0 -13.0 -2.0 2.0 116.851 116.852
RPP D15_66 -17.0 -13.0 -2.0 2.0 122.851 122.852
RPP D15_72 -17.0 -13.0 -2.0 2.0 128.851 128.852
RPP D15_78 -17.0 -13.0 -2.0 2.0 134.851 134.852
RPP D15_84 -17.0 -13.0 -2.0 2.0 140.851 140.852
RPP D15_90 -17.0 -13.0 -2.0 2.0 146.848 146.849
*
RPP D30_0 -32.0 -28.0 -2.0 2.0 56.851 56.852
RPP D30_6 -32.0 -28.0 -2.0 2.0 62.851 62.852
RPP D30_12 -32.0 -28.0 -2.0 2.0 68.851 68.852
RPP D30_18 -32.0 -28.0 -2.0 2.0 74.851 74.852
RPP D30_24 -32.0 -28.0 -2.0 2.0 80.851 80.852
RPP D30_30 -32.0 -28.0 -2.0 2.0 86.851 86.852
RPP D30_36 -32.0 -28.0 -2.0 2.0 92.851 92.852
RPP D30_42 -32.0 -28.0 -2.0 2.0 98.851 98.852
RPP D30_48 -32.0 -28.0 -2.0 2.0 104.851 104.852
RPP D30_54 -32.0 -28.0 -2.0 2.0 110.851 110.852
RPP D30_60 -32.0 -28.0 -2.0 2.0 116.851 116.852
RPP D30_66 -32.0 -28.0 -2.0 2.0 122.851 122.852
RPP D30_72 -32.0 -28.0 -2.0 2.0 128.851 128.852
RPP D30_78 -32.0 -28.0 -2.0 2.0 134.851 134.852
RPP D30_84 -32.0 -28.0 -2.0 2.0 140.851 140.852
RPP D30_90 -32.0 -28.0 -2.0 2.0 146.848 146.849
*
RPP D45_0 -47.0 -43.0 -2.0 2.0 56.851 56.852
RPP D45_6 -47.0 -43.0 -2.0 2.0 62.851 62.852
RPP D45_12 -47.0 -43.0 -2.0 2.0 68.851 68.852
RPP D45_18 -47.0 -43.0 -2.0 2.0 74.851 74.852
RPP D45_24 -47.0 -43.0 -2.0 2.0 80.851 80.852
RPP D45_30 -47.0 -43.0 -2.0 2.0 86.851 86.852
RPP D45_36 -47.0 -43.0 -2.0 2.0 92.851 92.852
RPP D45_42 -47.0 -43.0 -2.0 2.0 98.851 98.852
RPP D45_48 -47.0 -43.0 -2.0 2.0 104.851 104.852
RPP D45_54 -47.0 -43.0 -2.0 2.0 110.851 110.852
RPP D45_60 -47.0 -43.0 -2.0 2.0 116.851 116.852
RPP D45_66 -47.0 -43.0 -2.0 2.0 122.851 122.852
RPP D45_72 -47.0 -43.0 -2.0 2.0 128.851 128.852
RPP D45_78 -47.0 -43.0 -2.0 2.0 134.851 134.852
RPP D45_84 -47.0 -43.0 -2.0 2.0 140.851 140.852
RPP D45_90 -47.0 -43.0 -2.0 2.0 146.848 146.849
END
* ** REGIONS
TARIRON 25 | +TARIRO
AIROUT 25 | +BlaCyl -TARIRO -CONP
* outer volumes
BLACKH 25 | -BlaCyl +InfCyl
*
CON1_1 25 | +CONP +CONP1 -CONP0 -D0_0 -D15_0 -D30_0
CON1_2 25 | +CONP +CONP2 -CONP1 -CONP0 -D0_6 -D15_6 -D30_6
CON1_3 25 | +CONP +CONP3 -CONP2 -CONP0 -D0_12 -D15_12 -D30_12
CON1_4 25 | +CONP +CONP4 -CONP3 -CONP0 -D0_18 -D15_18 -D30_18
CON1_5 25 | +CONP +CONP5 -CONP4 -CONP0 -D0_24 -D15_24 -D30_24
CON1_6 25 | +CONP +CONP6 -CONP5 -CONP0 -D0_30 -D15_30 -D30_30
CON1_7 25 | +CONP +CONP7 -CONP6 -CONP0 -D0_36 -D15_36 -D30_36
CON1_8 25 | +CONP +CONP8 -CONP7 -CONP0 -D0_42 -D15_42 -D30_42
CON1_9 25 | +CONP +CONP9 -CONP8 -CONP0 -D0_48 -D15_48 -D30_48
CON1_10 25 | +CONP +CONP10 -CONP9 -CONP0 -D0_54 -D15_54 -D30_54
CON1_11 25 | +CONP +CONP11 -CONP10 -CONP0 -D0_60 -D15_60 -D30_60
CON1_12 25 | +CONP +CONP12 -CONP11 -CONP0 -D0_66 -D15_66 -D30_66
CON1_13 25 | +CONP +CONP13 -CONP12 -CONP0 -D0_72 -D15_72 -D30_72
CON1_14 25 | +CONP +CONP14 -CONP13 -CONP0 -D0_78 -D15_78 -D30_78
CON1_15 25 | +CONP -CONP14 -CONP0 -D0_84 -D15_84 -D30_84
    -D0_90 -D15_90 -D30_90
*
CON2_1 25 | +CONP +CONP1 +CONP0 -D45_0
CON2_2 25 | +CONP +CONP2 -CONP1 +CONP0 -D45_6
CON2_3 25 | +CONP +CONP3 -CONP2 +CONP0 -D45_12
CON2_4 25 | +CONP +CONP4 -CONP3 +CONP0 -D45_18
CON2_5 25 | +CONP +CONP5 -CONP4 +CONP0 -D45_24
CON2_6 25 | +CONP +CONP6 -CONP5 +CONP0 -D45_30
CON2_7 25 | +CONP +CONP7 -CONP6 +CONP0 -D45_36
CON2_8 25 | +CONP +CONP8 -CONP7 +CONP0 -D45_42
CON2_9 25 | +CONP +CONP9 -CONP8 +CONP0 -D45_48
CON2_10 25 | +CONP +CONP10 -CONP9 +CONP0 -D45_54
CON2_11 25 | +CONP +CONP11 -CONP10 +CONP0 -D45_60
CON2_12 25 | +CONP +CONP12 -CONP11 +CONP0 -D45_66
CON2_13 25 | +CONP +CONP13 -CONP12 +CONP0 -D45_72
CON2_14 25 | +CONP +CONP14 -CONP13 +CONP0 -D45_78
CON2_15 25 | +CONP -CONP14 +CONP0 -D45_84 -D45_90
*
D1_0 25 | +D0_0
D1_1 25 | +D0_6
D1_2 25 | +D0_12
D1_3 25 | +D0_18
D1_4 25 | +D0_24
D1_5 25 | +D0_30
D1_6 25 | +D0_36
D1_7 25 | +D0_42
D1_8 25 | +D0_48
D1_9 25 | +D0_54
D1_10 25 | +D0_60
D1_11 25 | +D0_66
D1_12 25 | +D0_72
D1_13 25 | +D0_78
D1_14 25 | +D0_84
D1_15 25 | +D0_90
*
D2_0 25 | +D15_0
D2_1 25 | +D15_6
D2_2 25 | +D15_12
D2_3 25 | +D15_18
D2_4 25 | +D15_24
D2_5 25 | +D15_30
D2_6 25 | +D15_36
D2_7 25 | +D15_42
D2_8 25 | +D15_48
D2_9 25 | +D15_54
D2_10 25 | +D15_60
D2_11 25 | +D15_66
D2_12 25 | +D15_72
D2_13 25 | +D15_78
D2_14 25 | +D15_84
D2_15 25 | +D15_90
*
D3_0 25 | +D30_0
D3_1 25 | +D30_6
D3_2 25 | +D30_12
D3_3 25 | +D30_18
D3_4 25 | +D30_24
D3_5 25 | +D30_30
D3_6 25 | +D30_36
D3_7 25 | +D30_42
D3_8 25 | +D30_48
D3_9 25 | +D30_54
D3_10 25 | +D30_60
D3_11 25 | +D30_66
D3_12 25 | +D30_72
D3_13 25 | +D30_78
D3_14 25 | +D30_84
D3_15 25 | +D30_90
*
D4_0 25 | +D45_0
D4_1 25 | +D45_6
D4_2 25 | +D45_12
D4_3 25 | +D45_18
D4_4 25 | +D45_24
D4_5 25 | +D45_30
D4_6 25 | +D45_36
D4_7 25 | +D45_42
D4_8 25 | +D45_48
D4_9 25 | +D45_54
D4_10 25 | +D45_60
D4_11 25 | +D45_66
D4_12 25 | +D45_72
D4_13 25 | +D45_78
D4_14 25 | +D45_84
D4_15 25 | +D45_90
END
GEOEND
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
* BIASING 3.0 1.0 1.51 CON1_1
* BIASING 3.0 1.0 1.12 CON1_2
* BIASING 3.0 1.0 1.0 CON1_3
* BIASING 3.0 1.0 1.04 CON1_4
* BIASING 3.0 1.0 1.14 CON1_5
* BIASING 3.0 1.0 1.38 CON1_6
* BIASING 3.0 1.0 1.72 CON1_7
* BIASING 3.0 1.0 2.38 CON1_8
* BIASING 3.0 1.0 4.04 CON1_9
* BIASING 3.0 1.0 17.8 CON1_10
* BIASING 3.0 1.0 22.4 CON1_11
* BIASING 3.0 1.0 45.3 CON1_12
* BIASING 3.0 1.0 88.4 CON1_13
* BIASING 3.0 1.0 164.0 CON1_14
* BIASING 3.0 1.0 328.0 CON1_15
* *
* BIASING 3.0 1.0 3.02 CON2_1
* BIASING 3.0 1.0 2.24 CON2_2
* BIASING 3.0 1.0 2.00 CON2_3
* BIASING 3.0 1.0 2.08 CON2_4
* BIASING 3.0 1.0 2.28 CON2_5
* BIASING 3.0 1.0 2.77 CON2_6
* BIASING 3.0 1.0 3.44 CON2_7
* BIASING 3.0 1.0 4.76 CON2_8
* BIASING 3.0 1.0 8.08 CON2_9
* BIASING 3.0 1.0 35.6 CON2_10
* BIASING 3.0 1.0 44.8 CON2_11
* BIASING 3.0 1.0 90.7 CON2_12
* BIASING 3.0 1.0 177.0 CON2_13
* BIASING 3.0 1.0 329.0 CON2_14
* BIASING 3.0 1.0 658.0 CON2_15
* ---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
DISCARD KAONS E+&E- MUONS 7.0 NEUTRIM
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* ?X?R?A?????O?Í ONE-WAY,FLUENCE?ōs?¤?×?«(-101.0)B?ʐÏ?Í?Ê?Ì?Ð?Ê?Ì?Ý?Ì?lB?~???Ì?ê?Ê?©A?¼?û?Ì?Ì?å?«?¢?Ê
USRBDX 99.0 8.0 -61.0 4.0 34.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 5.0 35.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 6.0 36.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 7.0 37.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 8.0 38.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 9.0 39.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 10.0 40.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 11.0 41.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 12.0 42.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 13.0 43.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 14.0 44.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 15.0 45.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 16.0 46.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 17.0 47.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 48.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 49.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
*
USRBDX 99.0 8.0 -61.0 4.0 50.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 5.0 51.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 6.0 52.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 7.0 53.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 8.0 54.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 9.0 55.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 10.0 56.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 11.0 57.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 12.0 58.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 13.0 59.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 14.0 60.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 15.0 61.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 16.0 62.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 17.0 63.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 64.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 65.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
*
USRBDX 99.0 8.0 -61.0 4.0 66.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 5.0 67.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 6.0 68.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 7.0 69.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 8.0 70.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 9.0 71.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 10.0 72.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 11.0 73.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 12.0 74.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 13.0 75.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 14.0 76.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 15.0 77.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 16.0 78.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 17.0 79.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 80.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 18.0 81.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
*
USRBDX 99.0 8.0 -61.0 19.0 82.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 20.0 83.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 21.0 84.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 22.0 85.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 23.0 86.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 24.0 87.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 25.0 88.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 26.0 89.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 27.0 90.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 28.0 91.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 29.0 92.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 30.0 93.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 31.0 94.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 32.0 95.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 33.0 96.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
USRBDX 99.0 8.0 -61.0 33.0 97.0 16.0F10
USRBDX 1.0 1.0E-14 200.0 &
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
USRBDX 99.0 1.0 -62.0 4.0 34.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 5.0 35.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 6.0 36.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 7.0 37.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 8.0 38.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 9.0 39.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 10.0 40.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 11.0 41.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 12.0 42.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 13.0 43.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 14.0 44.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 15.0 45.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 16.0 46.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 17.0 47.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 48.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 49.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 1.0 -62.0 4.0 50.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 5.0 51.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 6.0 52.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 7.0 53.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 8.0 54.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 9.0 55.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 10.0 56.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 11.0 57.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 12.0 58.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 13.0 59.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 14.0 60.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 15.0 61.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 16.0 62.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 17.0 63.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 64.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 65.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 1.0 -62.0 4.0 66.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 5.0 67.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 6.0 68.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 7.0 69.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 8.0 70.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 9.0 71.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 10.0 72.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 11.0 73.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 12.0 74.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 13.0 75.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 14.0 76.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 15.0 77.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 16.0 78.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 17.0 79.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 80.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 18.0 81.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 1.0 -62.0 19.0 82.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 20.0 83.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 21.0 84.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 22.0 85.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 23.0 86.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 24.0 87.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 25.0 88.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 26.0 89.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 27.0 90.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 28.0 91.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 29.0 92.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 30.0 93.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 31.0 94.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 32.0 95.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 33.0 96.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 1.0 -62.0 33.0 97.0 16.0F10
USRBDX 1.0 0.001 125.0 &
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
USRBDX 99.0 209.0 -63.0 4.0 34.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 5.0 35.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 6.0 36.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 7.0 37.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 8.0 38.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 9.0 39.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 10.0 40.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 11.0 41.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 12.0 42.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 13.0 43.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 14.0 44.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 15.0 45.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 16.0 46.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 17.0 47.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 48.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 49.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 209.0 -63.0 4.0 50.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 5.0 51.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 6.0 52.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 7.0 53.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 8.0 54.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 9.0 55.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 10.0 56.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 11.0 57.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 12.0 58.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 13.0 59.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 14.0 60.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 15.0 61.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 16.0 62.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 17.0 63.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 64.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 65.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 209.0 -63.0 4.0 66.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 5.0 67.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 6.0 68.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 7.0 69.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 8.0 70.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 9.0 71.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 10.0 72.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 11.0 73.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 12.0 74.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 13.0 75.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 14.0 76.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 15.0 77.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 16.0 78.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 17.0 79.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 80.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 18.0 81.0 16.0F10
USRBDX 1.0 0.001 125.0 &
*
USRBDX 99.0 209.0 -63.0 19.0 82.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 20.0 83.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 21.0 84.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 22.0 85.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 23.0 86.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 24.0 87.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 25.0 88.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 26.0 89.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 27.0 90.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 28.0 91.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 29.0 92.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 30.0 93.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 31.0 94.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 32.0 95.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 33.0 96.0 16.0F10
USRBDX 1.0 0.001 125.0 &
USRBDX 99.0 209.0 -63.0 33.0 97.0 16.0F10
USRBDX 1.0 0.001 125.0 &
* ...+....1....+....2....+....3....+....4....+....5....+....6....+....7...
RANDOMIZ 1.0 1.0
START 3E5
STOP

*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
      SUBROUTINE SOURCE ( NOMORE )

      INCLUDE '(DBLPRC)'
      INCLUDE '(DIMPAR)'
      INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2006 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* New source for FLUKA9x-FLUKA200x: *
* *
* Created on 07 january 1990 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 03-mar-06 by Alfredo Ferrari *
* *
* This is just an example of a possible user written source routine. *
* note that the beam card still has some meaning - in the scoring the *
* maximum momentum used in deciding the binning is taken from the *
* beam momentum. Other beam card parameters are obsolete. *
* *
*----------------------------------------------------------------------*
*
      INCLUDE '(BEAMCM)'
      INCLUDE '(FHEAVY)'
      INCLUDE '(FLKSTK)'
      INCLUDE '(IOIOCM)'
      INCLUDE '(LTCLCM)'
      INCLUDE '(PAPROP)'
      INCLUDE '(SOURCM)'
      INCLUDE '(SUMCOU)'
*
      LOGICAL LFIRST
*
      integer*4 sample
      integer*4 score
      parameter(score = 900000)
      dimension xcoo(score), ycoo(score), zcoo(score), weigh(score)
      dimension cosx(score), cosy(score), cosz(score), jpart(score)
      dimension enrgy(score)
      SAVE LFIRST, kount
      DATA LFIRST / .TRUE. /
*======================================================================*
* *
* BASIC VERSION *
* *
*======================================================================*
      NOMORE = 0
* +-------------------------------------------------------------------*
* | First call initializations:
      IF ( LFIRST ) THEN
* | *** The following 3 cards are mandatory ***
         TKESUM = ZERZER
         LFIRST = .FALSE.
         LUSSRC = .TRUE.
* | *** User initialization ***
        write(lunout,*)
        write(lunout,'(a,132a)') ("*",i=1,132)
        write(lunout,*)
        write(lunout,*)
        write(lunout,'(a)') "@ Source: collision file around shield11"
          open(unit=28,file=
     $ "/afs/slac.stanford.edu/u/sa/komi/Himac/Cion/secondaries.dat",
     $ status="old")
        open(unit=29,file="../sampledso",status="unknown")
        write(lunout,*)
        write(lunout,*)'@ Beam momentum [GeV/c]:',Pbeam
        write(lunout,'(a,132a)') ("*",i=1,132)
        write(lunout,*)
        do i = 1, (score - 1)
           read(28,*) jpart(i), xcoo(i), ycoo(i), zcoo(i), cosx(i),
     $ cosy(i), cosz(i), enrgy(i), weigh(i)
        end do
        kount = 0
      END IF
* |
* +-------------------------------------------------------------------*
      kount = kount + 1
* Sample a line from the file
      xi = flrndm(dummy)
      sample = int(xi * dble(score)) + 1
      if (sample.eq.score) then
        sample=score-1
      end if
      if(kount .le. 50) then
         write(29,*) sample,xcoo(sample),ycoo(sample),zcoo(sample),
     $ enrgy(sample),weigh(sample)
      end if
* Push one source particle to the stack. Note that you could as well
* push many but this way we reserve a maximum amount of space in the
* stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
      NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
      WTFLK (NPFLKA) = weigh(sample)
      WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
* +-------------------------------------------------------------------*
* | (Radioactive) isotope:
      IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
         IARES = IPROA
         IZRES = IPROZ
         IISRES = IPROM
         CALL STISBM ( IARES, IZRES, IISRES )
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
* |
* +-------------------------------------------------------------------*
* | Heavy ion:
      ELSE IF ( IJBEAM .EQ. -2 ) THEN
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
         ILOFLK (NPFLKA) = IJHION
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
* |
* +-------------------------------------------------------------------*
* | Normal hadron:
      ELSE
         IONID = IJBEAM
* ILOFLK (NPFLKA) = IJBEAM
* Warning!!!! in the input file BEAM is of electrons!!!!!
         ILOFLK (NPFLKA) = jpart(sample)
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
      END IF
* |
* +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
      LOFLK (NPFLKA) = 1
* User dependent flag:
      LOUSE (NPFLKA) = 0
* User dependent spare variables:
      DO 100 ISPR = 1, MKBMX1
         SPAREK (ISPR,NPFLKA) = ZERZER
 100 CONTINUE
* User dependent spare flags:
      DO 200 ISPR = 1, MKBMX2
         ISPARK (ISPR,NPFLKA) = 0
 200 CONTINUE
* Save the track number of the stack particle:
      ISPARK (MKBMX2,NPFLKA) = NPFLKA
      NPARMA = NPARMA + 1
      NUMPAR (NPFLKA) = NPARMA
      NEVENT (NPFLKA) = 0
      DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
      AGESTK (NPFLKA) = +ZERZER
      AKNSHR (NPFLKA) = -TWOTWO
* Group number for "low" energy neutrons, set to 0 anyway
      IGROUP (NPFLKA) = 0
* Kinetic energy of the particle (GeV)
      TKEFLK (NPFLKA) = enrgy(sample)
* Particle momentum
      PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
     & + TWOTWO * AM (ILOFLK(NPFLKA)) ) )
* Cosines (tx,ty,tz)
      TXFLK (NPFLKA) = cosx(sample)
      TYFLK (NPFLKA) = cosy(sample)
      TZFLK (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
     & - TYFLK (NPFLKA)**2 )
* Polarization cosines:
      TXPOL (NPFLKA) = -TWOTWO
      TYPOL (NPFLKA) = +ZERZER
      TZPOL (NPFLKA) = +ZERZER
* Particle coordinates
      XFLK (NPFLKA) = xcoo(sample)
      YFLK (NPFLKA) = ycoo(sample)
      ZFLK (NPFLKA) = zcoo(sample)
* Calculate the total kinetic energy of the primaries: don't change
          IF ( ILOFLK (NPFLKA) .EQ. -2 .OR.
     & ILOFLK (NPFLKA) .GT. 2000000 )
     & THEN

         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
         TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) )
     & * WTFLK (NPFLKA)
      ELSE
         TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
      END IF
      RADDLY (NPFLKA) = ZERZER
* Here we ask for the region number of the hitting point.
* NREG (NPFLKA) = ...
* The following line makes the starting region search much more
* robust if particles are starting very close to a boundary:
      CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
      CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
     & NRGFLK(NPFLKA), IDISC )
* Do not change these cards:
      CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
      NLATTC (NPFLKA) = MLATTC
      CMPATH (NPFLKA) = ZERZER
      CALL SOEVSV
      RETURN
*=== End of subroutine Source =========================================*
      END
Received on Mon Jan 11 2010 - 10:22:08 CET

This archive was generated by hypermail 2.2.0 : Mon Jan 11 2010 - 10:22:13 CET