1 ************************************************************************
2 * Energy dependent factors for the conversion of fluence to *
3 * effective dose and ambient dose equivalent for neutrons, protons, *
4 * charged pions, muons, photons and electrons. *
6 * The following sets are available: *
7 * (for photons and electrons see note below) *
9 * 1-3: Effective dose from ICRP74 and Pelliccioni data *
10 * calculated with ICRP radiation weighting factors Wr *
11 * 1: Anterior-Posterior irradiation *
12 * 2: Rotational irradiation geometry *
13 * 3: WORST possible geometry for the irradiation *
15 * 4-6: Effective dose from ICRP74 and Pelliccioni data *
16 * calculated with the Pelliccioni radiation weighting *
18 * 4: Anterior-Posterior irradiation *
19 * 5: Rotational irradiation geometry *
20 * 6: WORST possible geometry for the irradiation *
22 * 7: Ambient dose equivalent from ICRP74 and Pelliccioni data *
24 * 8: Ambient dose equivalent with old "GRS"-conversion factors *
26 * The different sets are invoked by EXTRAWEI (WHAT(3)=1) and the *
27 * first five characters of the SDUM input parameter of the USRBIN, *
28 * USRTRACK or USRBDX cards (case-insensitive) : *
31 * 'EAP74' = 1 'ERT74' = 2 'EWT74' = 3 *
32 * 'EAPMP' = 4 'ERTMP' = 5 'EWTMP' = 6 *
33 * 'AMB74' = 7 'AMBGS' = 8 *
35 * The default (any other SDUM) is unit weight factor. *
37 * Conversion factor sets for hadrons (n,p,pi) and muons can be *
38 * enabled separately from those for photons and electrons by *
39 * setting the 6th character of SDUM: *
41 * SDUM(6:6) = 1 only hadrons and muons (zero factor returned for *
42 * electrons and photons) *
43 * = 2 only electrons/positrons and photons (zero factor *
44 * returned for hadrons and muons) *
45 * otherwise: all particles considered (default) *
47 * Note: Photons and electrons: *
48 * Only sets 1-3 and 7 are implemented for photons and electrons*
49 * If sets 4-6 are requested sets 1-3 are used instead. *
50 * For set 8 zero values are returned. *
52 * Note 2: Stand-alone use: *
53 * The function can be used "stand-alone" (i.e. independent *
54 * from FLUKA). In this case different conversion factor sets *
55 * are invoked by the first parameter (IIJ) *
57 * IIJ = -(is*10000+id) where *
59 * is = index of parameter set *
60 * id = particle identity *
62 * (example: -80008 = neutron with GRS-conversion factor set) *
63 ************************************************************************
65 *=== fluscw ===========================================================*
67 DOUBLE PRECISION FUNCTION FLUSCW
68 & (IIJ,PLA,TXX,TYY,TZZ,WEE,XX,YY,ZZ,NREG,IOLREG,LLO,NSURF)
85 PARAMETER (NBIN1N=81,NBIN2N=150,
86 & NBIN1P=64,NBIN2P=61,
88 & NBIN1M=71,NBIN2M=16,
90 & NBIN1E=41,NBIN2E=81)
93 * 1 - neutron, 2 - proton, 3 - pion+, 4 - pion-, 5 - muon+, 5 - muon-,
94 * 7 - photon, 8 - electron, 99 - the rest
95 DATA ITT/ 2, 2, 8, 8,99,99, 7, 1, 1, 5,
96 & 5,99, 3, 4, 3, 4, 1, 1,99, 4,
97 & 3,99,99,99,99,99,99,99,99,99,
98 & 99,99,99,99,99,99,99,99,99,99/
100 *-----------------------------------------------------------------------
102 * Effective dose for neutrons from ICRP74 and Pelliccioni data
103 * calculated with ICRP radiation weighting factors Wr
104 * Energy in GeV, dose in pSv.cm^2
106 DIMENSION EBINN(NBIN1N),AEBINN(NBIN1N)
108 1 1.000D-12, 1.000D-11, 2.530D-11, 1.000D-10, 2.000D-10,
109 2 5.000D-10, 1.000D-09, 2.000D-09, 5.000D-09, 1.000D-08,
110 3 2.000D-08, 5.000D-08, 1.000D-07, 2.000D-07, 5.000D-07,
111 4 1.000D-06, 2.000D-06, 5.000D-06, 1.000D-05, 2.000D-05,
112 5 3.000D-05, 5.000D-05, 7.000D-05, 1.000D-04, 1.500D-04,
113 6 2.000D-04, 3.000D-04, 5.000D-04, 7.000D-04, 9.000D-04,
114 7 1.000D-03, 1.200D-03, 2.000D-03, 3.000D-03, 4.000D-03,
115 8 5.000D-03, 6.000D-03, 7.000D-03, 8.000D-03, 9.000D-03,
116 9 1.000D-02, 1.200D-02, 1.400D-02, 1.500D-02, 1.600D-02,
117 * 1.800D-02, 2.000D-02, 3.000D-02, 5.000D-02, 7.500D-02,
118 1 1.000D-01, 1.500D-01, 2.000D-01, 3.000D-01, 5.000D-01,
119 2 7.000D-01, 1.000D+00, 1.500D+00, 2.000D+00, 3.000D+00,
120 3 5.000D+00, 7.000D+00, 1.000D+01, 1.500D+01, 2.000D+01,
121 4 3.000D+01, 5.000D+01, 7.000D+01, 1.000D+02, 1.500D+02,
122 5 2.000D+02, 3.000D+02, 5.000D+02, 7.000D+02, 1.000D+03,
123 6 1.500D+03, 2.000D+03, 3.000D+03, 5.000D+03, 7.000D+03,
125 DIMENSION AP74N(NBIN1N),RT74N(NBIN1N),WT74N(NBIN1N),
126 & AAP74N(NBIN1N),ART74N(NBIN1N),AWT74N(NBIN1N)
127 * Anterior-Posterior irradiation
129 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
130 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
131 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
132 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
133 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
134 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
135 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
136 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
137 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
138 * 486.0D0, 480.0D0, 458.0D0, 437.0D0, 429.0D0,
139 1 429.0D0, 431.0D0, 447.0D0, 491.0D0, 581.0D0,
140 2 663.0D0, 767.0D0, 905.0D0, 1016.0D0, 1191.0D0,
141 3 1443.0D0, 1628.0D0, 1840.0D0, 2097.0D0, 2291.0D0,
142 4 2581.0D0, 2978.0D0, 3263.0D0, 3590.0D0, 3999.0D0,
143 5 4315.0D0, 4801.0D0, 5484.0D0, 5980.0D0, 6550.0D0,
144 6 7255.0D0, 7794.0D0, 8614.0D0, 9751.0D0, 10569.0D0,
146 * Rotational irradiation geometry
148 1 2.99D0, 3.72D0, 4.40D0, 5.75D0, 6.43D0,
149 2 7.27D0, 7.84D0, 8.31D0, 8.72D0, 8.90D0,
150 3 8.92D0, 8.82D0, 8.69D0, 8.56D0, 8.40D0,
151 4 8.34D0, 8.39D0, 9.06D0, 10.6D0, 13.8D0,
152 5 16.9D0, 22.7D0, 27.8D0, 34.8D0, 45.4D0,
153 6 54.8D0, 71.6D0, 99.4D0, 123.0D0, 144.0D0,
154 7 154.0D0, 173.0D0, 234.0D0, 283.0D0, 315.0D0,
155 8 335.0D0, 348.0D0, 358.0D0, 366.0D0, 373.0D0,
156 9 378.0D0, 385.0D0, 390.0D0, 391.0D0, 393.0D0,
157 * 394.0D0, 395.0D0, 395.0D0, 404.0D0, 422.0D0,
158 1 443.0D0, 496.0D0, 535.0D0, 594.0D0, 681.0D0,
159 2 754.0D0, 854.0D0, 1008.0D0, 1149.0D0, 1396.0D0,
160 3 1792.0D0, 2102.0D0, 2460.0D0, 2884.0D0, 3192.0D0,
161 4 3639.0D0, 4238.0D0, 4669.0D0, 5180.0D0, 5855.0D0,
162 5 6406.0D0, 7296.0D0, 8634.0D0, 9664.0D0, 10900.0D0,
163 6 12505.0D0, 13788.0D0, 15828.0D0, 18844.0D0, 21145.0D0,
165 * WORST possible geometry for the irradiation
167 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
168 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
169 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
170 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
171 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
172 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
173 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
174 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
175 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
176 * 486.0D0, 480.0D0, 458.0D0, 444.0D0, 459.0D0,
177 1 477.0D0, 523.0D0, 559.0D0, 616.0D0, 707.0D0,
178 2 785.0D0, 892.0D0, 1056.0D0, 1205.0D0, 1468.0D0,
179 3 1898.0D0, 2244.0D0, 2660.0D0, 3184.0D0, 3586.0D0,
180 4 4201.0D0, 5065.0D0, 5702.0D0, 6450.0D0, 7415.0D0,
181 5 8184.0D0, 9409.0D0, 11224.0D0, 12618.0D0, 14300.0D0,
182 6 16512.0D0, 18304.0D0, 21196.0D0, 25562.0D0, 28960.0D0,
185 * Effective dose for protons from ICRP74 and Pelliccioni data
186 * calculated with ICRP radiation weighting factors Wr
187 * Energy in GeV, dose in pSv.cm^2
189 DIMENSION E74P(NBIN1P),AE74P(NBIN1P)
191 1 5.00D-03, 6.50D-03, 8.00D-03, 1.00D-02, 1.25D-02,
192 2 1.50D-02, 2.00D-02, 2.50D-02, 3.00D-02, 4.00D-02,
193 3 5.00D-02, 6.50D-02, 8.00D-02, 1.00D-01, 1.25D-01,
194 4 1.50D-01, 2.00D-01, 2.50D-01, 3.00D-01, 4.00D-01,
195 5 5.00D-01, 6.50D-01, 8.00D-01, 1.00D+00, 1.25D+00,
196 6 1.50D+00, 2.00D+00, 2.50D+00, 3.00D+00, 4.00D+00,
197 7 5.00D+00, 6.50D+00, 8.00D+00, 1.00D+01, 1.25D+01,
198 8 1.50D+01, 2.00D+01, 2.50D+01, 3.00D+01, 4.00D+01,
199 9 5.00D+01, 6.50D+01, 8.00D+01, 1.00D+02, 1.25D+02,
200 * 1.50D+02, 2.00D+02, 2.50D+02, 3.00D+02, 4.00D+02,
201 1 5.00D+02, 6.50D+02, 8.00D+02, 1.00D+03, 1.25D+03,
202 2 1.50D+03, 2.00D+03, 2.50D+03, 3.00D+03, 4.00D+03,
203 3 5.00D+03, 6.50D+03, 8.00D+03, 1.00D+04/
204 DIMENSION AP74P(NBIN1P),RT74P(NBIN1P),WT74P(NBIN1P),
205 & AAP74P(NBIN1P),ART74P(NBIN1P),AWT74P(NBIN1P)
206 * Anterior-Posterior irradiation
208 1 57.5D0, 65.4D0, 74.3D0, 89.5D0, 116.4D0,
209 2 155.2D0, 289.0D0, 546.5D0, 982.1D0, 2540.0D0,
210 3 4810.0D0, 6969.6D0, 7336.4D0, 6820.0D0, 6080.5D0,
211 4 5467.0D0, 4570.0D0, 3977.7D0, 3573.9D0, 3089.0D0,
212 5 2840.0D0, 2673.7D0, 2597.5D0, 2530.0D0, 2455.9D0,
213 6 2410.1D0, 2420.0D0, 2534.3D0, 2688.1D0, 2996.8D0,
214 7 3240.0D0, 3464.8D0, 3601.1D0, 3730.0D0, 3862.6D0,
215 8 3975.3D0, 4155.7D0, 4292.5D0, 4398.5D0, 4547.2D0,
216 9 4640.0D0, 4722.1D0, 4784.0D0, 4870.0D0, 4995.2D0,
217 * 5128.6D0, 5394.4D0, 5645.8D0, 5879.5D0, 6298.2D0,
218 1 6662.4D0, 7129.9D0, 7525.5D0, 7970.0D0, 8428.5D0,
219 2 8812.8D0, 9438.1D0, 9940.3D0, 10362.6D0, 11053.4D0,
220 3 11611.9D0, 12297.0D0, 12863.8D0, 13500.0D0/
221 * Rotational irradiation geometry
223 1 48.2D0, 54.5D0, 61.1D0, 71.4D0, 87.6D0,
224 2 108.5D0, 169.7D0, 266.8D0, 409.7D0, 875.8D0,
225 3 1582.0D0, 2499.1D0, 3087.1D0, 3642.5D0, 4289.0D0,
226 4 4804.9D0, 5160.0D0, 4771.6D0, 4200.7D0, 3304.3D0,
227 5 2832.5D0, 2593.8D0, 2546.2D0, 2540.0D0, 2518.5D0,
228 6 2501.2D0, 2527.5D0, 2630.4D0, 2767.2D0, 3059.7D0,
229 7 3330.0D0, 3664.4D0, 3928.8D0, 4200.0D0, 4446.4D0,
230 8 4629.0D0, 4888.9D0, 5074.6D0, 5222.2D0, 5461.2D0,
231 9 5665.0D0, 5941.3D0, 6192.5D0, 6497.5D0, 6842.1D0,
232 * 7154.8D0, 7708.2D0, 8190.8D0, 8621.5D0, 9370.8D0,
233 1 10013.3D0, 10839.7D0, 11549.4D0, 12367.5D0, 13243.7D0,
234 2 14005.1D0, 15295.5D0, 16377.1D0, 17317.0D0, 18910.3D0,
235 3 20245.9D0, 21937.1D0, 23374.8D0, 25025.0D0/
236 * WORST possible geometry for the irradiation
238 1 59.2D0, 66.8D0, 75.5D0, 90.6D0, 117.3D0,
239 2 155.8D0, 289.0D0, 545.6D0, 980.1D0, 2540.0D0,
240 3 4810.0D0, 6860.6D0, 7172.3D0, 6820.0D0, 6501.3D0,
241 4 6255.9D0, 5640.0D0, 4889.3D0, 4237.9D0, 3379.9D0,
242 5 2950.0D0, 2718.2D0, 2661.4D0, 2650.0D0, 2639.8D0,
243 6 2636.5D0, 2670.0D0, 2751.5D0, 2859.1D0, 3109.6D0,
244 7 3380.0D0, 3787.8D0, 4155.1D0, 4550.0D0, 4897.3D0,
245 8 5139.5D0, 5459.3D0, 5672.7D0, 5839.6D0, 6122.3D0,
246 9 6390.0D0, 6794.6D0, 7181.0D0, 7650.0D0, 8163.8D0,
247 * 8616.8D0, 9397.9D0, 10064.3D0, 10651.5D0, 11662.8D0,
248 1 12524.9D0, 13633.2D0, 14589.0D0, 15700.0D0, 16904.4D0,
249 2 17963.1D0, 19781.3D0, 21326.4D0, 22683.3D0, 25011.0D0,
250 3 26986.3D0, 29515.1D0, 31686.5D0, 34200.0D0/
252 * Effective dose for charged pions from ICRP74 and Pelliccioni data
253 * calculated with ICRP radiation weighting factors Wr
254 * Energy in GeV, dose in pSv.cm^2
256 DIMENSION EBINI(NBIN1I),AEBINI(NBIN1I)
258 1 1.00D-03, 1.25D-03, 1.50D-03, 2.00D-03, 2.50D-03,
259 2 3.00D-03, 4.00D-03, 5.00D-03, 6.50D-03, 8.00D-03,
260 3 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
261 4 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
262 5 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
263 6 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
264 7 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
265 8 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
266 9 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
267 * 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
268 1 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
269 2 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
270 3 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
271 4 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
273 DIMENSION AP74I(NBIN1I),RT74I(NBIN1I),WT74I(NBIN1I),
274 & AP74J(NBIN1I),RT74J(NBIN1I),WT74J(NBIN1I),
275 & AAP74I(NBIN1I),ART74I(NBIN1I),AWT74I(NBIN1I),
276 & AAP74J(NBIN1I),ART74J(NBIN1I),AWT74J(NBIN1I)
277 * Anterior-Posterior irradiation
280 1 71.3D0, 76.3D0, 80.9D0, 89.3D0, 97.4D0,
281 2 105.6D0, 122.3D0, 140.1D0, 168.8D0, 200.6D0,
282 3 248.0D0, 315.4D0, 389.4D0, 547.5D0, 707.6D0,
283 4 859.7D0, 1117.7D0, 1300.0D0, 1445.7D0, 1498.3D0,
284 5 1498.6D0, 1454.3D0, 1396.7D0, 1290.0D0, 1209.1D0,
285 6 1148.3D0, 1065.0D0, 1013.1D0, 968.2D0, 946.4D0,
286 7 938.0D0, 946.1D0, 962.9D0, 1002.2D0, 1038.6D0,
287 8 1068.4D0, 1107.9D0, 1125.0D0, 1124.7D0, 1115.3D0,
288 9 1105.0D0, 1101.1D0, 1103.4D0, 1116.2D0, 1133.3D0,
289 * 1151.6D0, 1188.2D0, 1222.6D0, 1269.5D0, 1311.1D0,
290 1 1360.0D0, 1412.9D0, 1459.0D0, 1537.5D0, 1603.3D0,
291 2 1660.3D0, 1756.4D0, 1836.3D0, 1936.3D0, 2020.2D0,
292 3 2115.0D0, 2214.8D0, 2300.1D0, 2442.1D0, 2558.6D0,
293 4 2658.2D0, 2823.7D0, 2959.4D0, 3127.6D0, 3267.5D0,
297 1 349.0D0, 359.7D0, 369.2D0, 387.2D0, 404.9D0,
298 2 422.9D0, 460.0D0, 499.0D0, 560.9D0, 626.8D0,
299 3 721.0D0, 847.1D0, 975.3D0, 1218.2D0, 1427.1D0,
300 4 1594.6D0, 1806.2D0, 1880.0D0, 1830.5D0, 1712.3D0,
301 5 1546.4D0, 1373.0D0, 1243.0D0, 1080.0D0, 996.1D0,
302 6 950.4D0, 909.5D0, 898.1D0, 902.1D0, 915.7D0,
303 7 938.0D0, 965.9D0, 991.4D0, 1033.7D0, 1065.3D0,
304 8 1088.4D0, 1115.5D0, 1125.0D0, 1121.7D0, 1113.2D0,
305 9 1105.0D0, 1102.9D0, 1106.1D0, 1119.7D0, 1137.0D0,
306 * 1155.2D0, 1191.3D0, 1225.1D0, 1271.0D0, 1311.9D0,
307 1 1360.0D0, 1412.3D0, 1458.1D0, 1536.2D0, 1601.9D0,
308 2 1658.9D0, 1755.2D0, 1835.3D0, 1935.6D0, 2019.8D0,
309 3 2115.0D0, 2215.1D0, 2300.6D0, 2442.7D0, 2559.4D0,
310 4 2659.0D0, 2824.4D0, 2960.0D0, 3128.0D0, 3267.8D0,
312 * Rotational irradiation geometry
315 1 37.0D0, 39.1D0, 40.9D0, 44.5D0, 48.0D0,
316 2 51.7D0, 59.5D0, 68.0D0, 82.3D0, 98.6D0,
317 3 123.5D0, 160.1D0, 201.6D0, 294.8D0, 395.2D0,
318 4 496.8D0, 686.5D0, 842.0D0, 1002.6D0, 1099.3D0,
319 5 1167.8D0, 1200.3D0, 1203.6D0, 1177.5D0, 1141.2D0,
320 6 1106.0D0, 1047.8D0, 1005.6D0, 964.7D0, 942.3D0,
321 7 931.0D0, 935.2D0, 949.0D0, 986.7D0, 1027.2D0,
322 8 1065.8D0, 1133.5D0, 1188.0D0, 1249.8D0, 1295.9D0,
323 9 1342.0D0, 1384.2D0, 1416.2D0, 1463.7D0, 1499.2D0,
324 * 1528.4D0, 1576.7D0, 1618.0D0, 1673.5D0, 1724.6D0,
325 1 1789.0D0, 1865.3D0, 1936.9D0, 2067.0D0, 2182.4D0,
326 2 2286.3D0, 2467.4D0, 2622.2D0, 2819.2D0, 2985.8D0,
327 3 3174.0D0, 3370.4D0, 3537.1D0, 3811.8D0, 4035.6D0,
328 4 4225.8D0, 4540.7D0, 4798.2D0, 5117.5D0, 5383.9D0,
332 1 194.8D0, 196.1D0, 197.6D0, 201.6D0, 206.9D0,
333 2 213.3D0, 228.3D0, 245.9D0, 275.9D0, 309.7D0,
334 3 360.2D0, 431.1D0, 506.5D0, 659.1D0, 802.8D0,
335 4 930.0D0, 1123.6D0, 1236.5D0, 1292.2D0, 1280.7D0,
336 5 1229.6D0, 1156.7D0, 1093.0D0, 1004.5D0, 956.1D0,
337 6 929.0D0, 904.6D0, 898.4D0, 902.7D0, 913.5D0,
338 7 931.0D0, 953.5D0, 975.4D0, 1015.8D0, 1052.0D0,
339 8 1084.5D0, 1140.7D0, 1188.0D0, 1246.7D0, 1293.6D0,
340 9 1342.0D0, 1386.3D0, 1419.5D0, 1468.1D0, 1503.8D0,
341 * 1532.9D0, 1580.6D0, 1621.1D0, 1675.4D0, 1725.5D0,
342 1 1789.0D0, 1864.5D0, 1935.7D0, 2065.3D0, 2180.6D0,
343 2 2284.4D0, 2465.8D0, 2620.8D0, 2818.3D0, 2985.3D0,
344 3 3174.0D0, 3370.8D0, 3537.7D0, 3812.8D0, 4036.7D0,
345 4 4227.0D0, 4541.7D0, 4799.1D0, 5118.1D0, 5384.2D0,
347 * WORST possible geometry for the irradiation
350 1 71.3D0, 76.3D0, 80.8D0, 89.2D0, 97.3D0,
351 2 105.5D0, 122.2D0, 139.9D0, 168.7D0, 200.5D0,
352 3 248.0D0, 315.6D0, 389.9D0, 548.5D0, 709.1D0,
353 4 861.5D0, 1119.2D0, 1300.0D0, 1442.3D0, 1491.8D0,
354 5 1489.9D0, 1445.3D0, 1389.6D0, 1290.0D0, 1217.9D0,
355 6 1165.1D0, 1093.7D0, 1049.0D0, 1008.9D0, 987.0D0,
356 7 974.0D0, 972.3D0, 979.2D0, 1004.7D0, 1037.1D0,
357 8 1072.1D0, 1144.0D0, 1215.0D0, 1315.6D0, 1402.7D0,
358 9 1495.0D0, 1577.3D0, 1635.9D0, 1714.7D0, 1766.7D0,
359 * 1805.1D0, 1862.6D0, 1908.6D0, 1969.8D0, 2028.0D0,
360 1 2105.0D0, 2201.5D0, 2295.6D0, 2472.3D0, 2633.3D0,
361 2 2780.5D0, 3041.0D0, 3265.9D0, 3553.6D0, 3796.8D0,
362 3 4070.0D0, 4351.8D0, 4588.0D0, 4972.1D0, 5280.3D0,
363 4 5539.3D0, 5962.8D0, 6305.5D0, 6726.6D0, 7076.1D0,
367 1 349.0D0, 359.6D0, 369.1D0, 386.9D0, 404.5D0,
368 2 422.4D0, 459.5D0, 498.5D0, 560.4D0, 626.5D0,
369 3 721.0D0, 847.7D0, 976.5D0, 1220.5D0, 1430.2D0,
370 4 1597.9D0, 1808.6D0, 1880.0D0, 1826.1D0, 1704.9D0,
371 5 1537.4D0, 1364.6D0, 1236.7D0, 1080.0D0, 1003.4D0,
372 6 964.3D0, 934.0D0, 929.9D0, 939.9D0, 955.0D0,
373 7 974.0D0, 992.6D0, 1008.2D0, 1036.3D0, 1063.9D0,
374 8 1092.2D0, 1151.7D0, 1215.0D0, 1312.1D0, 1400.0D0,
375 9 1495.0D0, 1579.7D0, 1639.9D0, 1720.2D0, 1772.5D0,
376 * 1810.8D0, 1867.4D0, 1912.4D0, 1972.1D0, 2029.2D0,
377 1 2105.0D0, 2200.5D0, 2294.1D0, 2470.1D0, 2631.0D0,
378 2 2778.2D0, 3038.8D0, 3264.1D0, 3552.4D0, 3796.2D0,
379 3 4070.0D0, 4352.4D0, 4589.0D0, 4973.4D0, 5281.8D0,
380 4 5540.9D0, 5964.3D0, 6306.8D0, 6727.5D0, 7076.6D0,
383 * Effective dose for photons from ICRP74 and Pelliccioni data
384 * calculated with ICRP radiation weighting factors Wr
385 * Energy in GeV, dose in pSv.cm^2
387 DIMENSION E74G(NBIN1G),AE74G(NBIN1G)
389 1 0.10000D-05, 0.12589D-05, 0.15849D-05, 0.19953D-05, 0.25119D-05,
390 2 0.31623D-05, 0.39811D-05, 0.50119D-05, 0.63096D-05, 0.79433D-05,
391 3 0.10000D-04, 0.12589D-04, 0.15849D-04, 0.19953D-04, 0.25119D-04,
392 4 0.31623D-04, 0.39811D-04, 0.50119D-04, 0.63096D-04, 0.79433D-04,
393 5 0.10000D-03, 0.12589D-03, 0.15849D-03, 0.19953D-03, 0.25119D-03,
394 6 0.31623D-03, 0.39811D-03, 0.50119D-03, 0.63096D-03, 0.79433D-03,
395 7 0.10000D-02, 0.12589D-02, 0.15849D-02, 0.19953D-02, 0.25119D-02,
396 8 0.31623D-02, 0.39811D-02, 0.50119D-02, 0.63096D-02, 0.79433D-02,
397 9 0.10000D-01, 0.12589D-01, 0.15849D-01, 0.19953D-01, 0.25119D-01,
398 * 0.31623D-01, 0.39811D-01, 0.50119D-01, 0.63096D-01, 0.79433D-01,
399 1 0.10000D+00, 0.12589D+00, 0.15849D+00, 0.19953D+00, 0.25119D+00,
400 2 0.31623D+00, 0.39811D+00, 0.50119D+00, 0.63096D+00, 0.79433D+00,
401 3 0.10000D+01, 0.12589D+01, 0.15849D+01, 0.19953D+01, 0.25119D+01,
402 4 0.31623D+01, 0.39811D+01, 0.50119D+01, 0.63096D+01, 0.79433D+01,
403 5 0.10000D+02, 0.12589D+02, 0.15849D+02, 0.19953D+02, 0.25119D+02,
404 6 0.31623D+02, 0.39811D+02, 0.50119D+02, 0.63096D+02, 0.79433D+02,
406 DIMENSION AP74G(NBIN1G),RT74G(NBIN1G),WT74G(NBIN1G),
407 & AAP74G(NBIN1G),ART74G(NBIN1G),AWT74G(NBIN1G)
408 * Anterior-Posterior irradiation
410 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
411 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18658D+02, 0.11621D+02,
412 3 0.71617D+01, 0.43328D+01, 0.26230D+01, 0.15912D+01, 0.97875D+00,
413 4 0.62511D+00, 0.45300D+00, 0.41541D+00, 0.45697D+00, 0.50270D+00,
414 5 0.55300D+00, 0.68970D+00, 0.86020D+00, 0.10728D+01, 0.13381D+01,
415 6 0.16688D+01, 0.20814D+01, 0.25952D+01, 0.31566D+01, 0.38395D+01,
416 7 0.46700D+01, 0.55604D+01, 0.65903D+01, 0.76977D+01, 0.89735D+01,
417 8 0.10487D+02, 0.12358D+02, 0.13617D+02, 0.15456D+02, 0.17725D+02,
418 9 0.21000D+02, 0.24813D+02, 0.29318D+02, 0.34640D+02, 0.41642D+02,
419 * 0.49614D+02, 0.57329D+02, 0.61958D+02, 0.67831D+02, 0.74261D+02,
420 1 0.81300D+02, 0.85676D+02, 0.90288D+02, 0.95149D+02, 0.98906D+02,
421 2 0.10280D+03, 0.10684D+03, 0.11104D+03, 0.11458D+03, 0.11823D+03,
422 3 0.12200D+03, 0.12364D+03, 0.12530D+03, 0.12698D+03, 0.12798D+03,
423 4 0.12898D+03, 0.12999D+03, 0.13101D+03, 0.13167D+03, 0.13233D+03,
424 5 0.13300D+03, 0.13452D+03, 0.13606D+03, 0.13761D+03, 0.13919D+03,
425 6 0.14078D+03, 0.14239D+03, 0.14401D+03, 0.14467D+03, 0.14533D+03,
427 * Rotational irradiation geometry
429 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
430 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18655D+02, 0.11608D+02,
431 3 0.69154D+01, 0.37516D+01, 0.20643D+01, 0.11725D+01, 0.61014D+00,
432 4 0.35636D+00, 0.27520D+00, 0.26661D+00, 0.30429D+00, 0.34763D+00,
433 5 0.39750D+00, 0.50246D+00, 0.63525D+00, 0.80327D+00, 0.10159D+01,
434 6 0.12851D+01, 0.16259D+01, 0.20570D+01, 0.25456D+01, 0.31506D+01,
435 7 0.39000D+01, 0.46762D+01, 0.56061D+01, 0.67125D+01, 0.79141D+01,
436 8 0.93241D+01, 0.10963D+02, 0.12696D+02, 0.14881D+02, 0.17439D+02,
437 9 0.20700D+02, 0.24597D+02, 0.29229D+02, 0.34738D+02, 0.41889D+02,
438 * 0.50210D+02, 0.59027D+02, 0.68568D+02, 0.78333D+02, 0.89610D+02,
439 1 0.10265D+03, 0.11038D+03, 0.11871D+03, 0.12770D+03, 0.13443D+03,
440 2 0.14150D+03, 0.14896D+03, 0.15681D+03, 0.16272D+03, 0.16887D+03,
441 3 0.17525D+03, 0.18162D+03, 0.18826D+03, 0.19518D+03, 0.19927D+03,
442 4 0.20343D+03, 0.20768D+03, 0.21203D+03, 0.21473D+03, 0.21747D+03,
443 5 0.22025D+03, 0.22378D+03, 0.22737D+03, 0.23102D+03, 0.23474D+03,
444 6 0.23852D+03, 0.24237D+03, 0.24628D+03, 0.24898D+03, 0.25173D+03,
446 * WORST possible geometry for the irradiation
448 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
449 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18658D+02, 0.11621D+02,
450 3 0.71617D+01, 0.43328D+01, 0.26230D+01, 0.15912D+01, 0.97875D+00,
451 4 0.62511D+00, 0.45300D+00, 0.41541D+00, 0.45697D+00, 0.50270D+00,
452 5 0.55300D+00, 0.68970D+00, 0.86020D+00, 0.10728D+01, 0.13381D+01,
453 6 0.16688D+01, 0.20814D+01, 0.25952D+01, 0.31566D+01, 0.38395D+01,
454 7 0.46700D+01, 0.55604D+01, 0.65903D+01, 0.76977D+01, 0.89735D+01,
455 8 0.10487D+02, 0.12358D+02, 0.13724D+02, 0.16218D+02, 0.18813D+02,
456 9 0.23000D+02, 0.27007D+02, 0.31713D+02, 0.37238D+02, 0.45249D+02,
457 * 0.53784D+02, 0.59748D+02, 0.71072D+02, 0.83565D+02, 0.98879D+02,
458 1 0.11700D+03, 0.12735D+03, 0.13861D+03, 0.15087D+03, 0.16030D+03,
459 2 0.17027D+03, 0.18087D+03, 0.19208D+03, 0.20036D+03, 0.20899D+03,
460 3 0.21800D+03, 0.22815D+03, 0.23877D+03, 0.24988D+03, 0.25669D+03,
461 4 0.26363D+03, 0.27076D+03, 0.27805D+03, 0.28262D+03, 0.28897D+03,
462 5 0.30000D+03, 0.30985D+03, 0.32002D+03, 0.33053D+03, 0.34138D+03,
463 6 0.35259D+03, 0.36417D+03, 0.37616D+03, 0.39209D+03, 0.40869D+03,
466 * Effective dose for electrons from ICRP74 and Pelliccioni data
467 * calculated with ICRP radiation weighting factors Wr
468 * Energy in GeV, dose in pSv.cm^2
470 DIMENSION E74E(NBIN1E),AE74E(NBIN1E)
472 1 0.50000D-02, 0.64047D-02, 0.82039D-02, 0.10509D-01, 0.13461D-01,
473 2 0.17242D-01, 0.22086D-01, 0.28291D-01, 0.36239D-01, 0.46420D-01,
474 3 0.59460D-01, 0.76165D-01, 0.97562D-01, 0.12497D+00, 0.16008D+00,
475 4 0.20505D+00, 0.26265D+00, 0.33644D+00, 0.43096D+00, 0.55203D+00,
476 5 0.70711D+00, 0.90576D+00, 0.11602D+01, 0.14861D+01, 0.19037D+01,
477 6 0.24385D+01, 0.31235D+01, 0.40010D+01, 0.51250D+01, 0.65647D+01,
478 7 0.84090D+01, 0.10771D+02, 0.13797D+02, 0.17673D+02, 0.22638D+02,
479 8 0.28998D+02, 0.37145D+02, 0.47580D+02, 0.60946D+02, 0.78068D+02,
481 DIMENSION AP74E(NBIN1E),RT74E(NBIN1E),WT74E(NBIN1E),
482 & AAP74E(NBIN1E),ART74E(NBIN1E),AWT74E(NBIN1E)
483 * Anterior-Posterior irradiation
485 1 0.11700D+03, 0.14636D+03, 0.18310D+03, 0.22254D+03, 0.24108D+03,
486 2 0.26117D+03, 0.28218D+03, 0.30367D+03, 0.32392D+03, 0.33732D+03,
487 3 0.34585D+03, 0.35438D+03, 0.36312D+03, 0.36782D+03, 0.37210D+03,
488 4 0.37679D+03, 0.38468D+03, 0.39275D+03, 0.40098D+03, 0.40825D+03,
489 5 0.41392D+03, 0.41968D+03, 0.42412D+03, 0.42769D+03, 0.43128D+03,
490 6 0.43436D+03, 0.43732D+03, 0.44030D+03, 0.44373D+03, 0.45113D+03,
491 7 0.45866D+03, 0.46611D+03, 0.47323D+03, 0.48045D+03, 0.48778D+03,
492 8 0.49522D+03, 0.50278D+03, 0.51045D+03, 0.51057D+03, 0.50878D+03,
494 * Rotational irradiation geometry
496 1 0.46000D+02, 0.59898D+02, 0.78750D+02, 0.10262D+03, 0.12495D+03,
497 2 0.15440D+03, 0.19064D+03, 0.23250D+03, 0.26904D+03, 0.29888D+03,
498 3 0.31637D+03, 0.33039D+03, 0.34523D+03, 0.35974D+03, 0.37500D+03,
499 4 0.39072D+03, 0.40301D+03, 0.41573D+03, 0.42889D+03, 0.44070D+03,
500 5 0.45012D+03, 0.45980D+03, 0.47394D+03, 0.49154D+03, 0.50997D+03,
501 6 0.52079D+03, 0.52974D+03, 0.53887D+03, 0.54941D+03, 0.57165D+03,
502 7 0.59498D+03, 0.61689D+03, 0.63350D+03, 0.65061D+03, 0.66823D+03,
503 8 0.68638D+03, 0.70507D+03, 0.72432D+03, 0.72548D+03, 0.72208D+03,
505 * WORST possible geometry for the irradiation
507 1 0.11700D+03, 0.14636D+03, 0.18310D+03, 0.22254D+03, 0.24108D+03,
508 2 0.26117D+03, 0.28218D+03, 0.30367D+03, 0.32392D+03, 0.34459D+03,
509 3 0.35514D+03, 0.36261D+03, 0.37023D+03, 0.37763D+03, 0.38513D+03,
510 4 0.39659D+03, 0.41277D+03, 0.42961D+03, 0.44714D+03, 0.46712D+03,
511 5 0.48808D+03, 0.50997D+03, 0.52813D+03, 0.54373D+03, 0.57298D+03,
512 6 0.58967D+03, 0.60557D+03, 0.63059D+03, 0.65754D+03, 0.69411D+03,
513 7 0.73271D+03, 0.77101D+03, 0.80533D+03, 0.84117D+03, 0.87861D+03,
514 8 0.91771D+03, 0.95856D+03, 0.10012D+04, 0.10295D+04, 0.10545D+04,
517 *-----------------------------------------------------------------------
519 * Effective dose for neutrons from ICRP74 and Pelliccioni data
520 * calculated with the Pelliccioni radiation weighting factors Wr
521 * Energy in GeV, dose in pSv.cm^2
523 DIMENSION APPLN(NBIN1N),RTPLN(NBIN1N),WTPLN(NBIN1N),
524 & AAPPLN(NBIN1N),ARTPLN(NBIN1N),AWTPLN(NBIN1N)
525 * Anterior-Posterior irradiation
527 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
528 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
529 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
530 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
531 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
532 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
533 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
534 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
535 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
536 * 486.0D0, 480.0D0, 458.0D0, 437.0D0, 429.0D0,
537 1 409.0D0, 385.0D0, 378.0D0, 379.0D0, 399.0D0,
538 2 422.0D0, 455.0D0, 503.0D0, 544.0D0, 610.0D0,
539 3 707.0D0, 776.0D0, 852.0D0, 937.0D0, 996.0D0,
540 4 1081.0D0, 1197.0D0, 1284.0D0, 1390.0D0, 1536.0D0,
541 5 1658.0D0, 1856.0D0, 2150.0D0, 2368.0D0, 2620.0D0,
542 6 2926.0D0, 3156.0D0, 3495.0D0, 3948.0D0, 4260.0D0,
544 * Rotational irradiation geometry
546 1 2.99D0, 3.72D0, 4.40D0, 5.75D0, 6.43D0,
547 2 7.27D0, 7.84D0, 8.31D0, 8.72D0, 8.90D0,
548 3 8.92D0, 8.82D0, 8.69D0, 8.56D0, 8.40D0,
549 4 8.34D0, 8.39D0, 9.06D0, 10.6D0, 13.8D0,
550 5 16.9D0, 22.7D0, 27.8D0, 34.8D0, 45.4D0,
551 6 54.8D0, 71.6D0, 99.4D0, 123.0D0, 144.0D0,
552 7 154.0D0, 173.0D0, 234.0D0, 283.0D0, 315.0D0,
553 8 335.0D0, 348.0D0, 358.0D0, 366.0D0, 373.0D0,
554 9 378.0D0, 385.0D0, 390.0D0, 391.0D0, 393.0D0,
555 * 394.0D0, 395.0D0, 395.0D0, 404.0D0, 422.0D0,
556 1 400.0D0, 417.0D0, 428.0D0, 441.0D0, 462.0D0,
557 2 482.0D0, 513.0D0, 566.0D0, 615.0D0, 702.0D0,
558 3 840.0D0, 947.0D0, 1070.0D0, 1216.0D0, 1324.0D0,
559 4 1484.0D0, 1706.0D0, 1871.0D0, 2070.0D0, 2336.0D0,
560 5 2555.0D0, 2910.0D0, 3447.0D0, 3862.0D0, 4360.0D0,
561 6 5005.0D0, 5520.0D0, 6338.0D0, 7542.0D0, 8458.0D0,
563 * WORST possible geometry for the irradiation
565 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
566 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
567 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
568 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
569 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
570 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
571 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
572 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
573 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
574 * 486.0D0, 480.0D0, 458.0D0, 444.0D0, 459.0D0,
575 1 460.0D0, 460.0D0, 461.0D0, 465.0D0, 481.0D0,
576 2 501.0D0, 535.0D0, 595.0D0, 653.0D0, 757.0D0,
577 3 927.0D0, 1062.0D0, 1220.0D0, 1413.0D0, 1557.0D0,
578 4 1775.0D0, 2080.0D0, 2307.0D0, 2580.0D0, 2944.0D0,
579 5 3242.0D0, 3726.0D0, 4460.0D0, 5030.0D0, 5720.0D0,
580 6 6623.0D0, 7352.0D0, 8520.0D0, 10267.0D0, 11614.0D0,
583 * Effective dose for protons from Pelliccioni data
584 * calculated with the Pelliccioni radiation weighting factors Wr
585 * Energy in GeV, dose in pSv.cm^2
587 DIMENSION EPLP(NBIN1P),AEPLP(NBIN1P)
589 1 5.00D-03, 6.50D-03, 8.00D-03, 1.00D-02, 1.25D-02,
590 2 1.50D-02, 2.00D-02, 2.50D-02, 3.00D-02, 4.00D-02,
591 3 5.00D-02, 6.50D-02, 8.00D-02, 1.00D-01, 1.25D-01,
592 4 1.50D-01, 2.00D-01, 2.50D-01, 3.00D-01, 4.00D-01,
593 5 5.00D-01, 6.50D-01, 8.00D-01, 1.00D+00, 1.25D+00,
594 6 1.50D+00, 2.00D+00, 2.50D+00, 3.00D+00, 4.00D+00,
595 7 5.00D+00, 6.50D+00, 8.00D+00, 1.00D+01, 1.25D+01,
596 8 1.50D+01, 2.00D+01, 2.50D+01, 3.00D+01, 4.00D+01,
597 9 5.00D+01, 6.50D+01, 8.00D+01, 1.00D+02, 1.25D+02,
598 * 1.50D+02, 2.00D+02, 2.50D+02, 3.00D+02, 4.00D+02,
599 1 5.00D+02, 6.50D+02, 8.00D+02, 1.00D+03, 1.25D+03,
600 2 1.50D+03, 2.00D+03, 2.50D+03, 3.00D+03, 4.00D+03,
601 3 5.00D+03, 6.50D+03, 8.00D+03, 1.00D+04/
602 DIMENSION APPLP(NBIN1P),RTPLP(NBIN1P),WTPLP(NBIN1P),
603 & AAPPLP(NBIN1P),ARTPLP(NBIN1P),AWTPLP(NBIN1P)
604 * Anterior-Posterior irradiation
606 1 23.0D0, 26.1D0, 29.7D0, 35.8D0, 46.6D0,
607 2 62.1D0, 115.6D0, 218.6D0, 392.8D0, 1016.0D0,
608 3 1924.0D0, 2787.8D0, 2934.5D0, 2728.0D0, 2432.2D0,
609 4 2186.8D0, 1828.0D0, 1591.1D0, 1429.6D0, 1235.6D0,
610 5 1136.0D0, 1069.5D0, 1039.0D0, 1012.0D0, 982.4D0,
611 6 964.0D0, 968.0D0, 1013.7D0, 1075.3D0, 1198.7D0,
612 7 1296.0D0, 1385.9D0, 1440.4D0, 1492.0D0, 1545.0D0,
613 8 1590.1D0, 1662.3D0, 1717.0D0, 1759.4D0, 1818.9D0,
614 9 1856.0D0, 1888.9D0, 1913.6D0, 1948.0D0, 1998.1D0,
615 * 2051.5D0, 2157.8D0, 2258.3D0, 2351.8D0, 2519.3D0,
616 1 2664.9D0, 2852.0D0, 3010.2D0, 3188.0D0, 3371.4D0,
617 2 3525.1D0, 3775.2D0, 3976.1D0, 4145.1D0, 4421.4D0,
618 3 4644.7D0, 4918.8D0, 5145.5D0, 5400.0D0/
619 * Rotational irradiation geometry
621 1 19.3D0, 21.8D0, 24.4D0, 28.6D0, 35.0D0,
622 2 43.4D0, 67.9D0, 106.7D0, 163.9D0, 350.3D0,
623 3 632.8D0, 999.7D0, 1234.8D0, 1457.0D0, 1715.6D0,
624 4 1921.9D0, 2064.0D0, 1908.6D0, 1680.3D0, 1321.7D0,
625 5 1133.0D0, 1037.5D0, 1018.5D0, 1016.0D0, 1007.4D0,
626 6 1000.5D0, 1011.0D0, 1052.2D0, 1106.9D0, 1223.9D0,
627 7 1332.0D0, 1465.8D0, 1571.5D0, 1680.0D0, 1778.6D0,
628 8 1851.6D0, 1955.6D0, 2029.8D0, 2088.9D0, 2184.5D0,
629 9 2266.0D0, 2376.5D0, 2477.0D0, 2599.0D0, 2736.9D0,
630 * 2861.9D0, 3083.3D0, 3276.3D0, 3448.6D0, 3748.3D0,
631 1 4005.3D0, 4335.9D0, 4619.7D0, 4947.0D0, 5297.5D0,
632 2 5602.0D0, 6118.2D0, 6550.8D0, 6926.8D0, 7564.1D0,
633 3 8098.4D0, 8774.8D0, 9349.9D0, 10010.0D0/
634 * WORST possible geometry for the irradiation
636 1 23.7D0, 26.7D0, 30.2D0, 36.2D0, 46.9D0,
637 2 62.3D0, 115.6D0, 218.2D0, 392.1D0, 1016.0D0,
638 3 1924.0D0, 2744.2D0, 2868.9D0, 2728.0D0, 2600.5D0,
639 4 2502.4D0, 2256.0D0, 1955.7D0, 1695.1D0, 1351.9D0,
640 5 1180.0D0, 1087.3D0, 1064.6D0, 1060.0D0, 1055.9D0,
641 6 1054.6D0, 1068.0D0, 1100.6D0, 1143.6D0, 1243.8D0,
642 7 1352.0D0, 1515.1D0, 1662.0D0, 1820.0D0, 1958.9D0,
643 8 2055.8D0, 2183.7D0, 2269.1D0, 2335.9D0, 2448.9D0,
644 9 2556.0D0, 2717.9D0, 2872.4D0, 3060.0D0, 3265.5D0,
645 * 3446.7D0, 3759.1D0, 4025.7D0, 4260.6D0, 4665.1D0,
646 1 5009.9D0, 5453.3D0, 5835.6D0, 6280.0D0, 6761.8D0,
647 2 7185.2D0, 7912.5D0, 8530.5D0, 9073.3D0, 10004.4D0,
648 3 10794.5D0, 11806.0D0, 12674.6D0, 13680.0D0/
650 * Effective dose for charged pions from ICRP74 and Pelliccioni data
651 * calculated with the Pelliccioni radiation weighting factors Wr
652 * Energy in GeV, dose in pSv.cm^2
654 DIMENSION APPLI(NBIN1I),RTPLI(NBIN1I),WTPLI(NBIN1I),
655 & APPLJ(NBIN1I),RTPLJ(NBIN1I),WTPLJ(NBIN1I),
656 & AAPPLI(NBIN1I),ARTPLI(NBIN1I),AWTPLI(NBIN1I),
657 & AAPPLJ(NBIN1I),ARTPLJ(NBIN1I),AWTPLJ(NBIN1I)
658 * Anterior-Posterior irradiation
661 1 71.3D0, 77.4D0, 82.9D0, 92.9D0, 102.3D0,
662 2 111.4D0, 129.4D0, 147.6D0, 175.8D0, 205.5D0,
663 3 248.0D0, 305.6D0, 366.3D0, 491.1D0, 613.9D0,
664 4 729.7D0, 930.0D0, 1083.3D0, 1233.9D0, 1320.6D0,
665 5 1378.7D0, 1401.7D0, 1397.4D0, 1357.9D0, 1308.7D0,
666 6 1262.3D0, 1186.8D0, 1133.1D0, 1082.0D0, 1054.7D0,
667 7 1042.0D0, 1048.9D0, 1066.7D0, 1110.3D0, 1151.5D0,
668 8 1185.4D0, 1230.6D0, 1250.0D0, 1249.3D0, 1238.7D0,
669 9 1228.0D0, 1226.1D0, 1231.9D0, 1253.2D0, 1279.4D0,
670 * 1306.6D0, 1359.6D0, 1408.7D0, 1474.7D0, 1532.7D0,
671 1 1600.0D0, 1671.9D0, 1734.0D0, 1838.2D0, 1924.4D0,
672 2 1998.2D0, 2121.1D0, 2221.8D0, 2346.0D0, 2448.6D0,
673 3 2563.0D0, 2681.5D0, 2781.4D0, 2945.1D0, 3077.4D0,
674 4 3189.2D0, 3372.7D0, 3521.6D0, 3704.2D0, 3855.1D0,
678 1 471.6D0, 506.5D0, 537.4D0, 592.1D0, 641.1D0,
679 2 686.9D0, 772.5D0, 853.8D0, 972.0D0, 1088.3D0,
680 3 1243.1D0, 1435.8D0, 1618.2D0, 1929.9D0, 2159.7D0,
681 4 2310.6D0, 2419.2D0, 2350.0D0, 2095.4D0, 1827.2D0,
682 5 1537.2D0, 1284.1D0, 1118.2D0, 939.1D0, 867.2D0,
683 6 841.0D0, 842.7D0, 871.4D0, 926.9D0, 981.2D0,
684 7 1042.0D0, 1097.6D0, 1137.2D0, 1188.2D0, 1217.5D0,
685 8 1234.5D0, 1249.0D0, 1250.0D0, 1242.0D0, 1233.6D0,
686 9 1228.0D0, 1230.4D0, 1238.5D0, 1262.0D0, 1288.7D0,
687 * 1315.6D0, 1367.4D0, 1414.9D0, 1478.6D0, 1534.6D0,
688 1 1600.0D0, 1670.4D0, 1731.5D0, 1834.7D0, 1920.6D0,
689 2 1994.5D0, 2117.8D0, 2219.1D0, 2344.3D0, 2447.7D0,
690 3 2563.0D0, 2682.2D0, 2782.6D0, 2946.8D0, 3079.3D0,
691 4 3191.2D0, 3374.6D0, 3523.1D0, 3705.3D0, 3855.8D0,
693 * Rotational irradiation geometry
696 1 37.0D0, 39.6D0, 42.0D0, 46.3D0, 50.4D0,
697 2 54.5D0, 62.9D0, 71.7D0, 85.7D0, 101.0D0,
698 3 123.5D0, 155.1D0, 189.7D0, 264.5D0, 342.9D0,
699 4 421.7D0, 571.2D0, 701.7D0, 855.7D0, 969.0D0,
700 5 1074.4D0, 1156.9D0, 1204.2D0, 1239.5D0, 1235.2D0,
701 6 1215.9D0, 1167.8D0, 1124.8D0, 1078.1D0, 1050.2D0,
702 7 1034.3D0, 1036.8D0, 1051.2D0, 1093.1D0, 1138.6D0,
703 8 1182.3D0, 1258.8D0, 1320.0D0, 1388.6D0, 1439.8D0,
704 9 1492.0D0, 1541.8D0, 1581.2D0, 1642.6D0, 1690.9D0,
705 * 1731.9D0, 1801.2D0, 1861.1D0, 1940.9D0, 2013.7D0,
706 1 2104.0D0, 2209.2D0, 2306.8D0, 2482.1D0, 2636.1D0,
707 2 2773.6D0, 3011.0D0, 3211.5D0, 3463.6D0, 3673.6D0,
708 3 3907.0D0, 4145.7D0, 4344.4D0, 4664.8D0, 4919.7D0,
709 4 5132.6D0, 5478.2D0, 5755.7D0, 6094.3D0, 6373.3D0,
713 1 263.2D0, 276.2D0, 287.6D0, 308.3D0, 327.7D0,
714 2 346.4D0, 383.5D0, 420.8D0, 478.2D0, 537.7D0,
715 3 621.1D0, 730.7D0, 840.4D0, 1044.2D0, 1215.0D0,
716 4 1347.7D0, 1504.9D0, 1545.6D0, 1479.2D0, 1366.5D0,
717 5 1222.2D0, 1081.8D0, 983.3D0, 873.5D0, 832.4D0,
718 6 822.1D0, 838.3D0, 871.8D0, 927.6D0, 978.9D0,
719 7 1034.3D0, 1083.5D0, 1118.8D0, 1167.6D0, 1202.1D0,
720 8 1229.9D0, 1277.1D0, 1320.0D0, 1380.7D0, 1434.0D0,
721 9 1492.0D0, 1547.0D0, 1589.5D0, 1653.8D0, 1702.9D0,
722 * 1743.6D0, 1811.3D0, 1869.0D0, 1945.9D0, 2016.2D0,
723 1 2104.0D0, 2207.2D0, 2303.5D0, 2477.5D0, 2631.1D0,
724 2 2768.5D0, 3006.4D0, 3207.8D0, 3461.2D0, 3672.3D0,
725 3 3907.0D0, 4146.8D0, 4346.2D0, 4667.4D0, 4922.7D0,
726 4 5135.7D0, 5481.2D0, 5758.2D0, 6096.0D0, 6374.3D0,
728 * WORST possible geometry for the irradiation
731 1 71.3D0, 77.4D0, 82.9D0, 92.9D0, 102.2D0,
732 2 111.3D0, 129.2D0, 147.4D0, 175.6D0, 205.4D0,
733 3 248.0D0, 305.8D0, 366.7D0, 492.1D0, 615.2D0,
734 4 731.2D0, 931.3D0, 1083.3D0, 1231.0D0, 1315.0D0,
735 5 1370.7D0, 1393.0D0, 1390.3D0, 1357.9D0, 1318.2D0,
736 6 1280.7D0, 1218.8D0, 1173.2D0, 1127.4D0, 1100.0D0,
737 7 1082.0D0, 1078.0D0, 1084.8D0, 1113.2D0, 1149.9D0,
738 8 1189.6D0, 1270.7D0, 1350.0D0, 1461.3D0, 1557.6D0,
739 9 1661.0D0, 1755.8D0, 1825.8D0, 1924.6D0, 1993.8D0,
740 * 2047.3D0, 2130.6D0, 2198.4D0, 2287.7D0, 2370.2D0,
741 1 2476.0D0, 2604.6D0, 2727.7D0, 2955.3D0, 3160.1D0,
742 2 3345.9D0, 3671.9D0, 3950.9D0, 4305.1D0, 4601.8D0,
743 3 4932.0D0, 5269.1D0, 5548.6D0, 5997.6D0, 6353.0D0,
744 4 6648.7D0, 7126.6D0, 7508.9D0, 7974.2D0, 8357.4D0,
748 1 471.6D0, 506.4D0, 537.2D0, 591.6D0, 640.5D0,
749 2 686.1D0, 771.6D0, 853.0D0, 971.2D0, 1087.8D0,
750 3 1243.1D0, 1436.7D0, 1620.1D0, 1933.5D0, 2164.4D0,
751 4 2315.4D0, 2422.4D0, 2350.0D0, 2090.4D0, 1819.3D0,
752 5 1528.2D0, 1276.2D0, 1112.5D0, 939.1D0, 873.5D0,
753 6 853.3D0, 865.4D0, 902.3D0, 965.7D0, 1023.3D0,
754 7 1082.0D0, 1128.0D0, 1156.5D0, 1191.3D0, 1215.9D0,
755 8 1238.9D0, 1289.7D0, 1350.0D0, 1452.7D0, 1551.2D0,
756 9 1661.0D0, 1761.9D0, 1835.7D0, 1938.1D0, 2008.2D0,
757 * 2061.5D0, 2142.8D0, 2208.1D0, 2293.7D0, 2373.2D0,
758 1 2476.0D0, 2602.2D0, 2723.8D0, 2949.7D0, 3153.9D0,
759 2 3339.6D0, 3666.1D0, 3946.2D0, 4301.9D0, 4600.1D0,
760 3 4932.0D0, 5270.5D0, 5551.0D0, 6001.1D0, 6357.0D0,
761 4 6652.8D0, 7130.5D0, 7512.2D0, 7976.6D0, 8358.7D0,
764 * Effective dose for muons from ICRP74 and Pelliccioni data
765 * calculated with the ICRP/Pelliccioni radiation weighting factors Wr
766 * Energy in GeV, dose in pSv.cm^2
768 DIMENSION EBINM(NBIN1M),AEBINM(NBIN1M)
770 1 1.00D-03, 1.25D-03, 1.50D-03, 2.00D-03, 2.50D-03,
771 2 3.00D-03, 4.00D-03, 5.00D-03, 6.50D-03, 8.00D-03,
772 3 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
773 4 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
774 5 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
775 6 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
776 7 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
777 8 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
778 9 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
779 * 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
780 1 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
781 2 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
782 3 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
783 4 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
785 DIMENSION APPLM(NBIN1M),RTPLM(NBIN1M),WTPLM(NBIN1M),
786 & AAPPLM(NBIN1M),ARTPLM(NBIN1M),AWTPLM(NBIN1M)
787 * Anterior-Posterior irradiation
789 1 177.0D0, 178.0D0, 180.0D0, 182.0D0, 184.0D0,
790 2 186.0D0, 188.0D0, 193.0D0, 200.0D0, 211.0D0,
791 3 243.0D0, 337.7D0, 449.7D0, 672.5D0, 830.5D0,
792 4 913.8D0, 921.6D0, 821.5D0, 635.1D0, 495.7D0,
793 5 385.0D0, 353.0D0, 340.0D0, 339.0D0, 337.0D0,
794 6 335.0D0, 334.0D0, 333.0D0, 332.0D0, 332.0D0,
795 7 332.5D0, 333.0D0, 334.0D0, 335.0D0, 336.0D0,
796 8 337.0D0, 338.0D0, 339.0D0, 340.0D0, 341.0D0,
797 9 340.5D0, 337.0D0, 335.0D0, 333.3D0, 332.9D0,
798 * 333.1D0, 334.2D0, 335.5D0, 337.4D0, 338.9D0,
799 1 340.5D0, 341.9D0, 343.0D0, 344.4D0, 345.3D0,
800 2 346.0D0, 346.8D0, 347.4D0, 348.1D0, 348.5D0,
801 3 349.0D0, 349.5D0, 349.9D0, 350.6D0, 351.2D0,
802 4 351.7D0, 352.7D0, 353.7D0, 355.0D0, 356.2D0,
804 * Rotational irradiation geometry
806 1 94.0D0, 96.0D0, 98.0D0, 99.0D0, 100.0D0,
807 2 101.0D0, 102.0D0, 105.0D0, 110.0D0, 116.0D0,
808 3 124.2D0, 155.8D0, 192.5D0, 275.5D0, 363.2D0,
809 4 445.4D0, 566.2D0, 614.1D0, 587.3D0, 529.0D0,
810 5 463.8D0, 414.5D0, 385.4D0, 354.2D0, 339.2D0,
811 6 331.4D0, 325.3D0, 324.5D0, 326.5D0, 329.4D0,
812 7 332.8D0, 335.9D0, 338.1D0, 340.8D0, 342.4D0,
813 8 343.5D0, 344.7D0, 345.4D0, 346.0D0, 346.3D0,
814 9 346.6D0, 346.9D0, 347.1D0, 347.4D0, 347.7D0,
815 * 347.9D0, 348.2D0, 348.5D0, 348.8D0, 349.0D0,
816 1 349.3D0, 349.6D0, 349.8D0, 350.1D0, 350.4D0,
817 2 350.6D0, 351.0D0, 351.2D0, 351.6D0, 351.8D0,
818 3 352.1D0, 352.4D0, 352.6D0, 353.0D0, 353.2D0,
819 4 353.5D0, 353.8D0, 354.1D0, 354.4D0, 354.7D0,
821 * WORST possible geometry for the irradiation
823 1 177.0D0, 180.0D0, 182.0D0, 184.0D0, 186.0D0,
824 2 188.0D0, 191.0D0, 195.0D0, 204.0D0, 215.0D0,
825 3 243.0D0, 339.2D0, 452.5D0, 672.5D0, 818.9D0,
826 4 892.6D0, 908.7D0, 847.5D0, 728.2D0, 628.6D0,
827 5 537.5D0, 470.8D0, 430.9D0, 386.9D0, 364.4D0,
828 6 351.6D0, 339.2D0, 334.4D0, 332.6D0, 333.1D0,
829 7 334.5D0, 336.2D0, 337.6D0, 339.9D0, 341.6D0,
830 8 343.0D0, 345.2D0, 346.8D0, 348.6D0, 349.8D0,
831 9 351.0D0, 352.0D0, 352.7D0, 353.5D0, 354.0D0,
832 * 354.3D0, 354.7D0, 354.8D0, 355.0D0, 355.0D0,
833 1 355.0D0, 355.0D0, 355.0D0, 354.9D0, 354.9D0,
834 2 354.9D0, 354.8D0, 354.7D0, 354.6D0, 354.6D0,
835 3 354.5D0, 354.4D0, 354.4D0, 354.3D0, 354.2D0,
836 4 354.2D0, 354.1D0, 354.1D0, 354.0D0, 354.0D0,
839 *-----------------------------------------------------------------------
841 * Ambient dose equivalent for neutrons from ICRP74 and Pelliccioni data
842 * Energy in GeV, dose in pSv.cm^2
844 DIMENSION AMBN(NBIN1N),AAMBN(NBIN1N)
846 1 6.6D0, 9.0D0, 10.6D0, 12.9D0, 13.5D0,
847 2 13.6D0, 13.3D0, 12.9D0, 12.0D0, 11.3D0,
848 3 10.6D0, 9.9D0, 9.4D0, 8.9D0, 8.3D0,
849 4 7.9D0, 7.7D0, 8.0D0, 10.5D0, 16.6D0,
850 5 23.7D0, 41.1D0, 60.0D0, 88.0D0, 132.0D0,
851 6 170.0D0, 233.0D0, 322.0D0, 375.0D0, 400.0D0,
852 7 416.0D0, 425.0D0, 420.0D0, 412.0D0, 408.0D0,
853 8 405.0D0, 400.0D0, 405.0D0, 409.0D0, 420.0D0,
854 9 440.0D0, 480.0D0, 520.0D0, 540.0D0, 555.0D0,
855 * 570.0D0, 600.0D0, 515.0D0, 400.0D0, 330.0D0,
856 * high-energy extrapolation according to Pelliccioni
857 1 285.0D0, 253.0D0, 247.0D0, 259.0D0, 298.0D0,
858 2 335.0D0, 377.0D0, 421.0D0, 447.0D0, 476.0D0,
859 3 499.0D0, 508.0D0, 511.0D0, 510.0D0, 508.0D0,
860 4 504.0D0, 501.0D0, 504.0D0, 511.0D0, 528.0D0,
861 5 545.0D0, 577.0D0, 629.0D0, 669.0D0, 717.0D0,
862 6 777.0D0, 824.0D0, 896.0D0, 998.0D0, 1073.0D0,
864 * high-energy extrapolation according to Sannikov and Savitskaya
865 C 1 285.0D0, 253.0D0, 285.0D0, 306.0D0, 420.0D0,
866 C 2 500.0D0, 647.0D0, 733.0D0, 789.0D0, 862.0D0,
867 C 3 951.0D0, 1000.0D0, 1050.0D0, 1050.0D0, 1046.0D0,
868 C 4 1038.0D0, 1031.0D0, 1038.0D0, 1052.0D0, 1087.0D0,
869 C 5 1122.0D0, 1188.0D0, 1295.0D0, 1377.0D0, 1476.0D0,
870 C 6 1600.0D0, 1696.0D0, 1845.0D0, 2055.0D0, 2209.0D0,
872 * 1cm-Pb-modified Snoopy (A&B) response
874 CC & 0.13272E+01,0.35280E+01,0.44158E+01,0.61860E+01,0.76848E+01,
875 CC & 0.91501E+01,0.10320E+02,0.11686E+02,0.13774E+02,0.15550E+02,
876 CC & 0.17568E+02,0.20671E+02,0.23400E+02,0.26452E+02,0.31107E+02,
877 CC & 0.35200E+02,0.40029E+02,0.47957E+02,0.55050E+02,0.63372E+02,
878 CC & 0.68650E+02,0.76109E+02,0.83508E+02,0.98300E+02,0.12183E+03,
879 CC & 0.14270E+03,0.18229E+03,0.25960E+03,0.33847E+03,0.39849E+03,
880 CC & 0.41700E+03,0.44177E+03,0.46149E+03,0.45284E+03,0.44749E+03,
881 CC & 0.44458E+03,0.43180E+03,0.40909E+03,0.38360E+03,0.36038E+03,
882 CC & 0.34000E+03,0.33184E+03,0.32509E+03,0.32212E+03,0.31936E+03,
883 CC & 0.31438E+03,0.31000E+03,0.34000E+03,0.32675E+03,0.32092E+03,
884 CC & 0.33000E+03,0.33581E+03,0.34000E+03,0.35000E+03,0.42979E+03,
885 CC & 0.47896E+03,0.53000E+03,0.10000E-12,0.10000E-12,0.10000E-12,
886 CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
887 CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
888 CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
889 CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
891 C & 0.13272E+01,0.35280E+01,0.44158E+01,0.61860E+01,0.76848E+01,
892 C & 0.91501E+01,0.10320E+02,0.11686E+02,0.13774E+02,0.15550E+02,
893 C & 0.17568E+02,0.20671E+02,0.23400E+02,0.26452E+02,0.31107E+02,
894 C & 0.35200E+02,0.40029E+02,0.47957E+02,0.55050E+02,0.63372E+02,
895 C & 0.68650E+02,0.76109E+02,0.83508E+02,0.98300E+02,0.12183E+03,
896 C & 0.14270E+03,0.18229E+03,0.25960E+03,0.33847E+03,0.39849E+03,
897 C & 0.41700E+03,0.44177E+03,0.46149E+03,0.45284E+03,0.44749E+03,
898 C & 0.44458E+03,0.43180E+03,0.40909E+03,0.38360E+03,0.36038E+03,
899 C & 0.34000E+03,0.33184E+03,0.32509E+03,0.32212E+03,0.31936E+03,
900 C & 0.31438E+03,0.31000E+03,0.34000E+03,0.32675E+03,0.32092E+03,
901 C & 0.33000E+03,0.33581E+03,0.34000E+03,0.35000E+03,0.42979E+03,
902 C & 0.47896E+03,0.53000E+03,0.58919E+03,0.63516E+03,0.70610E+03,
903 C & 0.80686E+03,0.88096E+03,0.96695E+03,0.10749E+04,0.11588E+04,
904 C & 0.12882E+04,0.14721E+04,0.16072E+04,0.17641E+04,0.19612E+04,
905 C & 0.21142E+04,0.23503E+04,0.26857E+04,0.29323E+04,0.32185E+04,
906 C & 0.35780E+04,0.38571E+04,0.42879E+04,0.48998E+04,0.53498E+04,
909 * Ambient dose equivalent for protons from Pelliccioni data
910 * Energy in GeV, dose in pSv.cm^2
912 DIMENSION EAMBP(NBIN2P),AEAMBP(NBIN2P)
914 1 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
915 2 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
916 3 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
917 4 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
918 5 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
919 6 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
920 7 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
921 8 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
922 9 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
923 * 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
924 1 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
925 2 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
927 DIMENSION AMBP(NBIN2P),AAMBP(NBIN2P)
929 1 .0D0, .0D0, .0D0, .0D0, .0D0,
930 2 .0D0, 4289.5D0, 2790.0D0, 2000.8D0, 1710.4D0,
931 3 1520.0D0, 1362.6D0, 1247.6D0, 1091.3D0, 991.0D0,
932 4 921.9D0, 835.3D0, 786.0D0, 746.5D0, 726.7D0,
933 5 714.4D0, 709.6D0, 710.3D0, 717.9D0, 727.7D0,
934 6 737.2D0, 753.3D0, 765.0D0, 776.0D0, 782.7D0,
935 7 788.1D0, 792.0D0, 794.4D0, 797.0D0, 798.6D0,
936 8 799.8D0, 802.2D0, 804.9D0, 809.5D0, 814.6D0,
937 9 822.0D0, 831.7D0, 841.2D0, 858.9D0, 874.6D0,
938 * 888.6D0, 912.4D0, 932.0D0, 955.9D0, 975.2D0,
939 1 996.0D0, 1016.4D0, 1032.5D0, 1056.4D0, 1073.1D0,
940 2 1085.2D0, 1100.8D0, 1109.3D0, 1114.7D0, 1115.1D0,
943 * Ambient dose equivalent for ch. pions from ICRP74 and Pelliccioni data
944 * Energy in GeV, dose in pSv.cm^2
946 DIMENSION AMBI(NBIN1I),AMBJ(NBIN1I),
947 & AAMBI(NBIN1I),AAMBJ(NBIN1I)
950 1 224.0D0, 245.3D0, 264.1D0, 296.5D0, 323.8D0,
951 2 347.5D0, 387.2D0, 419.9D0, 459.9D0, 492.5D0,
952 3 528.0D0, 563.6D0, 592.9D0, 640.3D0, 678.9D0,
953 4 712.1D0, 769.2D0, 819.0D0, 884.5D0, 938.6D0,
954 5 994.0D0, 1040.7D0, 1067.6D0, 1080.0D0, 1059.4D0,
955 6 1026.9D0, 958.6D0, 900.6D0, 835.9D0, 792.4D0,
956 7 757.0D0, 735.3D0, 726.0D0, 722.9D0, 727.1D0,
957 8 733.1D0, 743.7D0, 750.0D0, 752.4D0, 751.2D0,
958 9 748.0D0, 743.9D0, 740.2D0, 734.0D0, 729.2D0,
959 * 725.4D0, 720.1D0, 716.7D0, 714.0D0, 712.9D0,
960 1 713.0D0, 714.7D0, 717.1D0, 722.7D0, 728.4D0,
961 2 733.8D0, 743.6D0, 752.1D0, 762.9D0, 772.0D0,
962 3 782.0D0, 792.2D0, 800.7D0, 814.4D0, 825.2D0,
963 4 834.1D0, 848.4D0, 859.8D0, 873.3D0, 884.2D0,
967 1 1130.0D0, 1245.2D0, 1345.3D0, 1507.9D0, 1631.0D0,
968 2 1723.7D0, 1843.3D0, 1902.4D0, 1921.2D0, 1890.7D0,
969 3 1810.0D0, 1684.7D0, 1563.7D0, 1368.4D0, 1231.1D0,
970 4 1135.8D0, 1024.1D0, 975.0D0, 960.9D0, 974.3D0,
971 5 1003.4D0, 1038.1D0, 1063.2D0, 1080.0D0, 1062.8D0,
972 6 1031.8D0, 964.1D0, 905.4D0, 839.2D0, 794.2D0,
973 7 757.0D0, 733.4D0, 722.6D0, 717.3D0, 720.2D0,
974 8 725.5D0, 735.8D0, 743.0D0, 747.7D0, 748.8D0,
975 9 748.0D0, 745.7D0, 742.9D0, 737.5D0, 732.8D0,
976 * 728.9D0, 722.9D0, 718.9D0, 715.2D0, 713.5D0,
977 1 713.0D0, 714.2D0, 716.4D0, 721.7D0, 727.4D0,
978 2 732.9D0, 742.8D0, 751.5D0, 762.6D0, 771.8D0,
979 3 782.0D0, 792.4D0, 801.0D0, 814.7D0, 825.5D0,
980 4 834.5D0, 848.8D0, 860.0D0, 873.5D0, 884.3D0,
983 * Ambient dose equivalent for muons from ICRP74 and Pelliccioni data
984 * Energy in GeV, dose in pSv.cm^2
986 DIMENSION AMBM(NBIN1M),AAMBM(NBIN1M)
988 1 394.0D0, 411.7D0, 426.6D0, 450.7D0, 469.5D0,
989 2 484.8D0, 508.3D0, 525.6D0, 544.3D0, 557.4D0,
990 3 569.5D0, 578.9D0, 583.9D0, 585.1D0, 579.3D0,
991 4 569.5D0, 543.9D0, 515.5D0, 474.1D0, 439.0D0,
992 5 402.6D0, 369.8D0, 347.0D0, 320.0D0, 307.2D0,
993 6 301.2D0, 297.8D0, 298.9D0, 302.8D0, 307.0D0,
994 7 311.5D0, 315.1D0, 317.3D0, 319.7D0, 320.7D0,
995 8 321.1D0, 321.2D0, 320.9D0, 320.5D0, 320.1D0,
996 9 320.0D0, 320.2D0, 320.5D0, 321.3D0, 322.0D0,
997 * 322.7D0, 323.7D0, 324.4D0, 324.9D0, 325.1D0,
998 1 325.0D0, 324.5D0, 323.9D0, 322.5D0, 321.3D0,
999 2 320.2D0, 318.6D0, 317.4D0, 316.2D0, 315.4D0,
1000 3 315.0D0, 314.9D0, 315.1D0, 315.9D0, 316.8D0,
1001 4 317.7D0, 319.5D0, 321.0D0, 323.0D0, 324.7D0,
1004 * Ambient dose equivalent for photons from ICRP74 and Pelliccioni data
1005 * Energy in GeV, dose in pSv.cm^2
1007 DIMENSION AMBG(NBIN1G),AAMBG(NBIN1G)
1009 1 0.29512E-08, 0.16406E-07, 0.91199E-07, 0.50699E-06, 0.28184E-05,
1010 2 0.15667E-04, 0.87096E-04, 0.48417E-03, 0.26915E-02, 0.14962E-01,
1011 3 0.83176E-01, 0.46238E+00, 0.93903E+00, 0.10505E+01, 0.93070E+00,
1012 4 0.76366E+00, 0.62173E+00, 0.52108E+00, 0.51322E+00, 0.55911E+00,
1013 5 0.61958E+00, 0.72946E+00, 0.92815E+00, 0.12264E+01, 0.15403E+01,
1014 6 0.19055E+01, 0.23480E+01, 0.27945E+01, 0.36249E+01, 0.42697E+01,
1015 7 0.51844E+01, 0.61688E+01, 0.71466E+01, 0.82509E+01, 0.95675E+01,
1016 8 0.10544E+02, 0.10715E+02, 0.10375E+02, 0.94276E+01, 0.91033E+01,
1017 9 0.87619E+01, 0.84957E+01, 0.83464E+01, 0.82813E+01, 0.82585E+01,
1018 * 0.82149E+01, 0.82566E+01, 0.86437E+01, 0.88675E+01, 0.89207E+01,
1019 1 0.90012E+01, 0.92640E+01, 0.96872E+01, 0.10186E+02, 0.10715E+02,
1020 2 0.11169E+02, 0.11561E+02, 0.11803E+02, 0.11885E+02, 0.11830E+02,
1021 3 0.11695E+02, 0.11535E+02, 0.11429E+02, 0.11508E+02, 0.11830E+02,
1022 4 0.12331E+02, 0.12882E+02, 0.13305E+02, 0.13428E+02, 0.13092E+02,
1023 5 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02,
1024 6 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02,
1027 * Ambient dose equivalent for electrons from ICRP74 and Pelliccioni data
1028 * Energy in GeV, dose in pSv.cm^2
1030 DIMENSION EAMBE(NBIN2E),AEAMBE(NBIN2E)
1032 1 0.20000E-02, 0.22896E-02, 0.26212E-02, 0.30008E-02, 0.34354E-02,
1033 2 0.39329E-02, 0.45025E-02, 0.51545E-02, 0.59010E-02, 0.67556E-02,
1034 3 0.77339E-02, 0.88540E-02, 0.10136E-01, 0.11604E-01, 0.13285E-01,
1035 4 0.15209E-01, 0.17411E-01, 0.19932E-01, 0.22819E-01, 0.26124E-01,
1036 5 0.29907E-01, 0.34238E-01, 0.39196E-01, 0.44873E-01, 0.51371E-01,
1037 6 0.58811E-01, 0.67328E-01, 0.77078E-01, 0.88241E-01, 0.10102E+00,
1038 7 0.11565E+00, 0.13240E+00, 0.15157E+00, 0.17352E+00, 0.19865E+00,
1039 8 0.22742E+00, 0.26036E+00, 0.29806E+00, 0.34122E+00, 0.39064E+00,
1040 9 0.44721E+00, 0.51198E+00, 0.58612E+00, 0.67101E+00, 0.76818E+00,
1041 * 0.87943E+00, 0.10068E+01, 0.11526E+01, 0.13195E+01, 0.15106E+01,
1042 1 0.17294E+01, 0.19798E+01, 0.22665E+01, 0.25948E+01, 0.29705E+01,
1043 2 0.34007E+01, 0.38932E+01, 0.44570E+01, 0.51025E+01, 0.58414E+01,
1044 3 0.66874E+01, 0.76559E+01, 0.87646E+01, 0.10034E+02, 0.11487E+02,
1045 4 0.13151E+02, 0.15055E+02, 0.17235E+02, 0.19731E+02, 0.22589E+02,
1046 5 0.25860E+02, 0.29605E+02, 0.33892E+02, 0.38801E+02, 0.44420E+02,
1047 6 0.50853E+02, 0.58217E+02, 0.66648E+02, 0.76300E+02, 0.87350E+02,
1050 DIMENSION AMBE(NBIN2E),AAMBE(NBIN2E)
1052 1 1.00000E-12, 0.16082E+03, 0.19885E+03, 0.32509E+03, 0.41665E+03,
1053 2 0.45092E+03, 0.44206E+03, 0.41582E+03, 0.38863E+03, 0.36471E+03,
1054 3 0.34608E+03, 0.33323E+03, 0.32465E+03, 0.32115E+03, 0.32097E+03,
1055 4 0.32279E+03, 0.32481E+03, 0.32583E+03, 0.32546E+03, 0.32233E+03,
1056 5 0.31724E+03, 0.31179E+03, 0.30738E+03, 0.30701E+03, 0.31035E+03,
1057 6 0.31401E+03, 0.31705E+03, 0.31769E+03, 0.31635E+03, 0.31477E+03,
1058 7 0.31497E+03, 0.31646E+03, 0.31903E+03, 0.32163E+03, 0.32354E+03,
1059 8 0.32444E+03, 0.32485E+03, 0.32398E+03, 0.32261E+03, 0.32087E+03,
1060 9 0.31913E+03, 0.31669E+03, 0.31472E+03, 0.31256E+03, 0.31087E+03,
1061 * 0.30919E+03, 0.30832E+03, 0.30848E+03, 0.30932E+03, 0.31016E+03,
1062 1 0.31154E+03, 0.31390E+03, 0.31573E+03, 0.31831E+03, 0.32056E+03,
1063 2 0.32277E+03, 0.32480E+03, 0.32730E+03, 0.32833E+03, 0.32986E+03,
1064 3 0.33037E+03, 0.33037E+03, 0.32972E+03, 0.32810E+03, 0.32810E+03,
1065 4 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
1066 5 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
1067 6 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
1070 *-----------------------------------------------------------------------
1072 * Ambient dose equivalent for neutrons ("GRS"-conversion factors)
1073 * Energy in GeV, dose in pSv.cm^2
1075 DIMENSION H10N(NBIN2N),H10PR(NBIN2N),H10PI(NBIN2N),H10MU(NBIN2M),
1078 1 7.62D0, 7.84D0, 8.06D0, 8.28D0, 8.50D0,
1079 2 8.71D0, 8.92D0, 9.12D0, 9.32D0, 9.50D0,
1080 3 9.67D0, 9.83D0, 9.97D0, 10.10D0, 10.21D0,
1081 4 10.30D0, 10.38D0, 10.43D0, 10.46D0, 10.47D0,
1082 5 10.47D0, 10.43D0, 10.38D0, 10.31D0, 10.22D0,
1083 6 10.11D0, 9.99D0, 9.85D0, 9.69D0, 9.52D0,
1084 7 9.35D0, 9.16D0, 8.96D0, 8.76D0, 8.56D0,
1085 8 8.35D0, 8.15D0, 7.94D0, 7.74D0, 7.54D0,
1086 9 7.35D0, 7.16D0, 6.99D0, 6.82D0, 6.67D0,
1087 * 6.53D0, 6.40D0, 6.30D0, 6.22D0, 6.16D0,
1088 1 6.13D0, 6.14D0, 6.19D0, 6.29D0, 6.45D0,
1089 2 6.69D0, 7.03D0, 7.49D0, 8.11D0, 8.93D0,
1090 3 10.04D0, 11.51D0, 13.47D0, 16.08D0, 19.58D0,
1091 4 24.24D0, 30.42D0, 38.53D0, 49.03D0, 62.35D0,
1092 5 78.83D0, 98.66D0, 121.76D0, 147.73D0, 175.91D0,
1093 6 205.40D0, 235.18D0, 264.24D0, 291.70D0, 316.84D0,
1094 7 339.20D0, 358.53D0, 374.78D0, 388.05D0, 398.56D0,
1095 8 406.61D0, 412.66D0, 417.53D0, 423.43D0, 436.87D0,
1096 9 472.90D0, 539.16D0, 598.43D0, 623.89D0, 595.0D0 ,
1097 * 524.0D0 , 365.0D0 , 323.0D0 , 321.0D0 , 323.0D0 ,
1098 1 325.0D0 , 328.0D0 , 333.0D0 , 338.0D0 , 344.0D0 ,
1099 2 353.0D0 , 371.0D0 , 482.0D0 , 581.0D0 , 677.0D0 ,
1100 3 777.0D0 , 882.0D0 , 1005.0D0 , 1125.0D0 , 1265.0D0 ,
1101 4 1410.0D0 , 1570.0D0 , 1740.0D0 , 1930.0D0 , 2130.0D0 ,
1102 5 2360.0D0 , 2600.0D0 , 2880.0D0 , 3180.0D0 , 3500.0D0 ,
1103 6 3860.0D0 , 4250.0D0 , 4700.0D0 , 5180.0D0 , 5610.0D0 ,
1104 7 6290.0D0 , 6900.0D0 , 7570.0D0 , 8280.0D0 , 9000.0D0 ,
1105 8 9830.0D0 , 10700.0D0 , 11600.0D0 , 12650.0D0 , 13600.0D0 ,
1106 9 14750.0D0 , 15900.0D0 , 17100.0D0 , 18450.0D0 , 19700.0D0 ,
1107 * 21100.0D0 , 22500.0D0 , 23900.0D0 , 25500.0D0 , 27000.0D0 /
1109 * Ambient dose equivalent for protons ("GRS"-conversion factors)
1110 * Energy in GeV, dose in pSv.cm^2
1113 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1114 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1115 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1116 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1117 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1118 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1119 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1120 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1121 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1122 * 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1123 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1124 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1125 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1126 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1127 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1128 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1129 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1130 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1131 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1132 * 3130.0D0, 2530.0D0, 2130.0D0, 1825.0D0, 1585.0D0,
1133 1 1390.0D0, 1230.0D0, 1105.0D0, 1000.0D0, 923.0D0,
1134 2 864.0D0, 828.0D0, 898.0D0, 965.0D0, 1040.0D0,
1135 3 1125.0D0, 1220.0D0, 1335.0D0, 1455.0D0, 1595.0D0,
1136 4 1740.0D0, 1905.0D0, 2080.0D0, 2270.0D0, 2480.0D0,
1137 5 2710.0D0, 2960.0D0, 3240.0D0, 3550.0D0, 3870.0D0,
1138 6 4240.0D0, 4630.0D0, 5090.0D0, 5570.0D0, 6000.0D0,
1139 7 6690.0D0, 7300.0D0, 7980.0D0, 8690.0D0, 9410.0D0,
1140 8 10250.0D0, 11100.0D0, 12050.0D0, 13100.0D0, 14050.0D0,
1141 9 15200.0D0, 16350.0D0, 17550.0D0, 18900.0D0, 20100.0D0,
1142 * 21600.0D0, 22300.0D0, 24400.0D0, 26000.0D0, 27500.0D0/
1144 * Ambient dose equivalent for charged pions ("GRS"-conversion factors)
1145 * Energy in GeV, dose in pSv.cm^2
1148 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1149 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1150 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1151 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1152 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1153 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1154 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1155 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1156 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1157 * 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1158 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1159 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1160 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1161 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1162 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1163 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1164 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1165 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1166 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
1167 * 1640.0D0, 1560.0D0, 1495.0D0, 1445.0D0, 1405.0D0,
1168 1 1375.0D0, 1355.0D0, 1345.0D0, 1335.0D0, 1330.0D0,
1169 2 1330.0D0, 1330.0D0, 1330.0D0, 1335.0D0, 1340.0D0,
1170 3 1345.0D0, 1350.0D0, 1360.0D0, 1485.0D0, 1630.0D0,
1171 4 1780.0D0, 1945.0D0, 2120.0D0, 2310.0D0, 2520.0D0,
1172 5 2750.0D0, 3000.0D0, 3280.0D0, 3580.0D0, 3910.0D0,
1173 6 4270.0D0, 4670.0D0, 5120.0D0, 5600.0D0, 6040.0D0,
1174 7 6720.0D0, 7330.0D0, 8010.0D0, 8720.0D0, 9440.0D0,
1175 8 10300.0D0, 11150.0D0, 12050.0D0, 13100.0D0, 14050.0D0,
1176 9 15200.0D0, 16350.0D0, 17550.0D0, 18900.0D0, 20200.0D0,
1177 * 21600.0D0, 23000.0D0, 24400.0D0, 26000.0D0, 27500.0D0/
1179 * Ambient dose equivalent for muons ("GRS"-conversion factors)
1180 * Energy in GeV, dose in pSv.cm^2
1182 DIMENSION EAMBM(NBIN2M),AEAMBM(NBIN2M)
1184 1 1.0D-02, 2.0D-02, 5.0D-02, 1.0D-01, 2.0D-01,
1185 2 5.0D-01, 1.0D+00, 2.0D+00, 5.0D+00, 1.0D+01,
1186 3 2.0D+01, 5.0D+01, 1.0D+02, 2.0D+02, 5.0D+02,
1189 1 2600.0D0, 2600.0D0, 2600.0D0, 2600.0D0, 350.0D0,
1190 2 330.0D0, 350.0D0, 360.0D0, 380.0D0, 390.0D0,
1191 3 410.0D0, 430.0D0, 450.0D0, 480.0D0, 560.0D0,
1194 *-----------------------------------------------------------------------
1196 DATA LFIRST /.TRUE.,.TRUE.,.TRUE.,.TRUE.,.TRUE.,
1197 & .TRUE.,.TRUE.,.TRUE.,.TRUE./
1200 * default is zero weighting
1205 * calculate logarithms of energy-bins at the very first call
1208 AEBINN(I) = LOG10(EBINN(I))
1209 AAP74N(I) = LOG10(AP74N(I))
1210 ART74N(I) = LOG10(RT74N(I))
1211 AWT74N(I) = LOG10(WT74N(I))
1212 AAPPLN(I) = LOG10(APPLN(I))
1213 ARTPLN(I) = LOG10(RTPLN(I))
1214 AWTPLN(I) = LOG10(WTPLN(I))
1215 AAMBN(I) = LOG10(AMBN(I))
1218 AE74P(I) = LOG10(E74P(I))
1219 AEPLP(I) = LOG10(EPLP(I))
1220 AAP74P(I) = LOG10(AP74P(I))
1221 ART74P(I) = LOG10(RT74P(I))
1222 AWT74P(I) = LOG10(WT74P(I))
1223 AAPPLP(I) = LOG10(APPLP(I))
1224 ARTPLP(I) = LOG10(RTPLP(I))
1225 AWTPLP(I) = LOG10(WTPLP(I))
1228 AEBINI(I) = LOG10(EBINI(I))
1229 AAP74I(I) = LOG10(AP74I(I))
1230 ART74I(I) = LOG10(RT74I(I))
1231 AWT74I(I) = LOG10(WT74I(I))
1232 AAP74J(I) = LOG10(AP74J(I))
1233 ART74J(I) = LOG10(RT74J(I))
1234 AWT74J(I) = LOG10(WT74J(I))
1235 AAPPLI(I) = LOG10(APPLI(I))
1236 ARTPLI(I) = LOG10(RTPLI(I))
1237 AWTPLI(I) = LOG10(WTPLI(I))
1238 AAPPLJ(I) = LOG10(APPLJ(I))
1239 ARTPLJ(I) = LOG10(RTPLJ(I))
1240 AWTPLJ(I) = LOG10(WTPLJ(I))
1241 AAMBI(I) = LOG10(AMBI(I))
1242 AAMBJ(I) = LOG10(AMBJ(I))
1245 AEBINM(I) = LOG10(EBINM(I))
1246 AAPPLM(I) = LOG10(APPLM(I))
1247 ARTPLM(I) = LOG10(RTPLM(I))
1248 AWTPLM(I) = LOG10(WTPLM(I))
1249 AAMBM(I) = LOG10(AMBM(I))
1252 AEAMBP(I) = LOG10(EAMBP(I))
1253 IF (AMBP(I).GT.0.0D0) THEN
1254 AAMBP(I) = LOG10(AMBP(I))
1256 AAMBP(I) = LOG10(ANGLGB)
1260 AEAMBM(I) = LOG10(EAMBM(I))
1261 AH10MU(I) = LOG10(H10MU(I))
1264 AE74G(I) = LOG10(E74G(I))
1265 AAP74G(I) = LOG10(AP74G(I))
1266 ART74G(I) = LOG10(RT74G(I))
1267 AWT74G(I) = LOG10(WT74G(I))
1268 AAMBG(I) = LOG10(AMBG(I))
1271 AE74E(I) = LOG10(E74E(I))
1272 AAP74E(I) = LOG10(AP74E(I))
1273 ART74E(I) = LOG10(RT74E(I))
1274 AWT74E(I) = LOG10(WT74E(I))
1277 AEAMBE(I) = LOG10(EAMBE(I))
1278 AAMBE(I) = LOG10(AMBE(I))
1284 IF (IJ.LT.-10000) THEN
1285 * test print for set number itest (i.e. -20001 gives the proton
1286 * conversion factors for set 2)
1287 * return particle-id of conversion factor set used for particle IJ
1290 IJ = -IJ-10000*ITEST
1291 IF (ITT(IJ).EQ.1) THEN
1293 ELSEIF (ITT(IJ).EQ.2) THEN
1295 ELSEIF (ITT(IJ).EQ.3) THEN
1297 ELSEIF (ITT(IJ).EQ.4) THEN
1299 ELSEIF (ITT(IJ).EQ.5) THEN
1301 ELSEIF (ITT(IJ).EQ.6) THEN
1303 ELSEIF (ITT(IJ).EQ.7) THEN
1305 ELSEIF (ITT(IJ).EQ.8) THEN
1307 ELSEIF (ITT(IJ).EQ.99) THEN
1313 * ISCRNG = 1 for usrbdx
1316 * return unit weight for all other detectors
1317 IF ((ISCRNG.NE.1).AND.(ISCRNG.NE.2).AND.(ISCRNG.NE.3)) THEN
1323 * At this point the detector is either usrbdx or usrbin or usrtrack.
1325 * Now check if any of the available conversion factor sets is requested
1326 * in which case the default conversion factor is zero. Otherwise the
1327 * default factor is one.
1328 IF (ITEST.EQ.0) THEN
1330 IF (ISCRNG.EQ.1) THEN
1331 CSET = TITUSX(JSCRNG)
1333 ELSEIF (ISCRNG.EQ.3) THEN
1334 CSET = TITUTC(JSCRNG)
1337 CSET = TITUSB(JSCRNG)
1339 IF ((CSET(1:3).EQ.'EAP').OR.(CSET(1:3).EQ.'eap').OR.
1340 & (CSET(1:3).EQ.'ERT').OR.(CSET(1:3).EQ.'ert').OR.
1341 & (CSET(1:3).EQ.'EWT').OR.(CSET(1:3).EQ.'ewt').OR.
1342 & (CSET(1:3).EQ.'AMB').OR.(CSET(1:3).EQ.'amb')) THEN
1345 * skip the rest for all other sdum's
1351 * At this point one of the available conversion factor sets is
1352 * requested by sdum.
1354 * check for particle type
1355 * return zero weight for all particles with ij < 1 or > 40 and for
1356 * all particles with flag=99 (see array ITT)
1357 IF ((IJ.LT.1).OR.(IJ.GT.40)) RETURN
1359 IF (IPART.EQ.99) RETURN
1362 IF (PLA.LT.0.0D0) THEN
1364 ELSEIF (PLA.GT.0.0D0) THEN
1365 EKIN = SQRT(PLA**2+AM(IJ)**2)-AM(IJ)
1371 * conversion factor set
1372 IF (ITEST.EQ.0) THEN
1373 * for muons there is only a Pelliccioni weighting factor set
1374 IF (IPART.EQ.5) THEN
1375 IF ((CSET(1:1).EQ.'E').AND.(CSET(4:5).EQ.'74'))
1377 IF ((CSET(1:1).EQ.'e').AND.(CSET(4:5).EQ.'74'))
1380 * for photons and electrons there is only a ICRP weighting factor set
1381 IF ((IPART.EQ.7).OR.(IPART.EQ.8)) THEN
1382 IF (((CSET(1:1).EQ.'E' ).OR.(CSET(1:1).EQ.'e' )).AND.
1383 & ((CSET(4:5).EQ.'MP').OR.(CSET(4:5).EQ.'mp'))) THEN
1385 * for photons and electrons there is no GRS amb. dose equivalent data set
1386 ELSEIF ((CSET(1:3).EQ.'AMB').OR.(CSET(1:3).EQ.'amb')) THEN
1387 IF ((CSET(4:5).EQ.'GS').OR.(CSET(4:5).EQ.'gs')) RETURN
1389 * skip if only hadronic part is requested
1390 IF (CSET(6:6).EQ.'1') RETURN
1392 * skip if only electromagnetic part is requested
1393 IF (CSET(6:6).EQ.'2') RETURN
1396 IF (IPART.EQ.5) THEN
1397 IF (ITEST.EQ.1) ITEST = 4
1398 IF (ITEST.EQ.2) ITEST = 5
1399 IF (ITEST.EQ.3) ITEST = 6
1401 IF ((IPART.EQ.7).OR.(IPART.EQ.8)) THEN
1402 IF (ITEST.EQ.4) ITEST = 1
1403 IF (ITEST.EQ.5) ITEST = 2
1404 IF (ITEST.EQ.6) ITEST = 3
1405 IF (ITEST.GE.9) RETURN
1410 * Effective dose (ICRP radiation weighting factors Wr)
1411 * Anterior-Posterior irradiation
1412 IF ((CSET(1:5).EQ.'EAP74').OR.(CSET(1:5).EQ.'eap74')
1413 & .OR.(ITEST.EQ.1)) THEN
1414 IF (LFIRST(1).AND.(ITEST.EQ.0)) THEN
1420 IF (IPART.EQ.1) THEN
1422 IF (AEKIN.LT.AEBINN(1)) THEN
1424 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
1425 FLUSCW = AP74N(NBIN1N)
1426 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
1428 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
1429 & *DBLE(NBIN1N-47)+47.0D0
1432 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
1435 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
1439 IF (IBINHI.GT.NBIN1N) THEN
1440 FLUSCW = AP74N(NBIN1N)
1442 FLUSCW =( AEKIN-AEBINN(IBINLO))/
1443 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1444 & (AAP74N(IBINHI)-AAP74N(IBINLO))+AAP74N(IBINLO)
1445 FLUSCW = 10.0D0**FLUSCW
1451 IF ((IBINHI-IBINLO).EQ.1) GOTO 11
1452 KK = (IBINHI+IBINLO)/2
1453 IF (AEKIN.LE.AEBINN(KK)) THEN
1460 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
1461 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1462 & (AAP74N(IBINHI)-AAP74N(IBINLO))+AAP74N(IBINLO)
1463 FLUSCW = 10.0D0**FLUSCW
1465 ELSEIF (IPART.EQ.2) THEN
1467 IF (AEKIN.LT.AE74P(1)) THEN
1469 ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
1470 FLUSCW = AP74P(NBIN1P)
1472 XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
1473 & *DBLE(NBIN1P-1)+1.0D0
1476 IF (AEKIN.LT.AE74P(IBINLO)) THEN
1479 ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
1483 IF (IBINHI.GT.NBIN1P) THEN
1484 FLUSCW = AP74P(NBIN1P)
1486 FLUSCW =( AEKIN-AE74P(IBINLO))/
1487 & (AE74P(IBINHI)-AE74P(IBINLO))*
1488 & (AAP74P(IBINHI)-AAP74P(IBINLO))+AAP74P(IBINLO)
1489 FLUSCW = 10.0D0**FLUSCW
1492 ELSEIF (IPART.EQ.3) THEN
1494 IF (AEKIN.LT.AEBINI(1)) THEN
1496 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1497 FLUSCW = AP74I(NBIN1I)
1499 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1500 & *DBLE(NBIN1I-1)+1.0D0
1503 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1506 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1510 IF (IBINHI.GT.NBIN1I) THEN
1511 FLUSCW = AP74I(NBIN1I)
1513 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1514 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1515 & (AAP74I(IBINHI)-AAP74I(IBINLO))+AAP74I(IBINLO)
1516 FLUSCW = 10.0D0**FLUSCW
1519 ELSEIF (IPART.EQ.4) THEN
1521 IF (AEKIN.LT.AEBINI(1)) THEN
1523 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1524 FLUSCW = AP74J(NBIN1I)
1526 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1527 & *DBLE(NBIN1I-1)+1.0D0
1530 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1533 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1537 IF (IBINHI.GT.NBIN1I) THEN
1538 FLUSCW = AP74J(NBIN1I)
1540 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1541 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1542 & (AAP74J(IBINHI)-AAP74J(IBINLO))+AAP74J(IBINLO)
1543 FLUSCW = 10.0D0**FLUSCW
1546 ELSEIF (IPART.EQ.7) THEN
1548 IF (AEKIN.LT.AE74G(1)) THEN
1550 ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
1551 FLUSCW = AP74G(NBIN1G)
1553 XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
1554 & *DBLE(NBIN1G-1)+1.0D0
1557 IF (AEKIN.LT.AE74G(IBINLO)) THEN
1560 ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
1564 IF (IBINHI.GT.NBIN1G) THEN
1565 FLUSCW = AP74G(NBIN1G)
1567 FLUSCW =( AEKIN-AE74G(IBINLO))/
1568 & (AE74G(IBINHI)-AE74G(IBINLO))*
1569 & (AAP74G(IBINHI)-AAP74G(IBINLO))+AAP74G(IBINLO)
1570 FLUSCW = 10.0D0**FLUSCW
1573 ELSEIF (IPART.EQ.8) THEN
1575 IF (AEKIN.LT.AE74E(1)) THEN
1577 ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
1578 FLUSCW = AP74E(NBIN1E)
1580 XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
1581 & *DBLE(NBIN1E-1)+1.0D0
1584 IF (AEKIN.LT.AE74E(IBINLO)) THEN
1587 ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
1591 IF (IBINHI.GT.NBIN1E) THEN
1592 FLUSCW = AP74E(NBIN1E)
1594 FLUSCW =( AEKIN-AE74E(IBINLO))/
1595 & (AE74E(IBINHI)-AE74E(IBINLO))*
1596 & (AAP74E(IBINHI)-AAP74E(IBINLO))+AAP74E(IBINLO)
1597 FLUSCW = 10.0D0**FLUSCW
1603 * Effective dose (ICRP radiation weighting factors Wr)
1604 * Rotational irradiation geometry
1605 ELSEIF ((CSET(1:5).EQ.'ERT74').OR.(CSET(1:5).EQ.'ert74')
1606 & .OR.(ITEST.EQ.2)) THEN
1607 IF (LFIRST(2).AND.(ITEST.EQ.0)) THEN
1613 IF (IPART.EQ.1) THEN
1615 IF (AEKIN.LT.AEBINN(1)) THEN
1617 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
1618 FLUSCW = RT74N(NBIN1N)
1619 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
1621 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
1622 & *DBLE(NBIN1N-47)+47.0D0
1625 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
1628 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
1632 IF (IBINHI.GT.NBIN1N) THEN
1633 FLUSCW = RT74N(NBIN1N)
1635 FLUSCW =( AEKIN-AEBINN(IBINLO))/
1636 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1637 & (ART74N(IBINHI)-ART74N(IBINLO))+ART74N(IBINLO)
1638 FLUSCW = 10.0D0**FLUSCW
1644 IF ((IBINHI-IBINLO).EQ.1) GOTO 21
1645 KK = (IBINHI+IBINLO)/2
1646 IF (AEKIN.LE.AEBINN(KK)) THEN
1653 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
1654 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1655 & (ART74N(IBINHI)-ART74N(IBINLO))+ART74N(IBINLO)
1656 FLUSCW = 10.0D0**FLUSCW
1658 ELSEIF (IPART.EQ.2) THEN
1660 IF (AEKIN.LT.AE74P(1)) THEN
1662 ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
1663 FLUSCW = RT74P(NBIN1P)
1665 XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
1666 & *DBLE(NBIN1P-1)+1.0D0
1669 IF (AEKIN.LT.AE74P(IBINLO)) THEN
1672 ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
1676 IF (IBINHI.GT.NBIN1P) THEN
1677 FLUSCW = RT74P(NBIN1P)
1679 FLUSCW =( AEKIN-AE74P(IBINLO))/
1680 & (AE74P(IBINHI)-AE74P(IBINLO))*
1681 & (ART74P(IBINHI)-ART74P(IBINLO))+ART74P(IBINLO)
1682 FLUSCW = 10.0D0**FLUSCW
1685 ELSEIF (IPART.EQ.3) THEN
1687 IF (AEKIN.LT.AEBINI(1)) THEN
1689 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1690 FLUSCW = RT74I(NBIN1I)
1692 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1693 & *DBLE(NBIN1I-1)+1.0D0
1696 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1699 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1703 IF (IBINHI.GT.NBIN1I) THEN
1704 FLUSCW = RT74I(NBIN1I)
1706 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1707 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1708 & (ART74I(IBINHI)-ART74I(IBINLO))+ART74I(IBINLO)
1709 FLUSCW = 10.0D0**FLUSCW
1712 ELSEIF (IPART.EQ.4) THEN
1714 IF (AEKIN.LT.AEBINI(1)) THEN
1716 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1717 FLUSCW = RT74J(NBIN1I)
1719 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1720 & *DBLE(NBIN1I-1)+1.0D0
1723 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1726 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1730 IF (IBINHI.GT.NBIN1I) THEN
1731 FLUSCW = RT74J(NBIN1I)
1733 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1734 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1735 & (ART74J(IBINHI)-ART74J(IBINLO))+ART74J(IBINLO)
1736 FLUSCW = 10.0D0**FLUSCW
1739 ELSEIF (IPART.EQ.7) THEN
1741 IF (AEKIN.LT.AE74G(1)) THEN
1743 ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
1744 FLUSCW = RT74G(NBIN1G)
1746 XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
1747 & *DBLE(NBIN1G-1)+1.0D0
1750 IF (AEKIN.LT.AE74G(IBINLO)) THEN
1753 ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
1757 IF (IBINHI.GT.NBIN1G) THEN
1758 FLUSCW = RT74G(NBIN1G)
1760 FLUSCW =( AEKIN-AE74G(IBINLO))/
1761 & (AE74G(IBINHI)-AE74G(IBINLO))*
1762 & (ART74G(IBINHI)-ART74G(IBINLO))+ART74G(IBINLO)
1763 FLUSCW = 10.0D0**FLUSCW
1766 ELSEIF (IPART.EQ.8) THEN
1768 IF (AEKIN.LT.AE74E(1)) THEN
1770 ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
1771 FLUSCW = RT74E(NBIN1E)
1773 XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
1774 & *DBLE(NBIN1E-1)+1.0D0
1777 IF (AEKIN.LT.AE74E(IBINLO)) THEN
1780 ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
1784 IF (IBINHI.GT.NBIN1E) THEN
1785 FLUSCW = RT74E(NBIN1E)
1787 FLUSCW =( AEKIN-AE74E(IBINLO))/
1788 & (AE74E(IBINHI)-AE74E(IBINLO))*
1789 & (ART74E(IBINHI)-ART74E(IBINLO))+ART74E(IBINLO)
1790 FLUSCW = 10.0D0**FLUSCW
1796 * Effective dose (ICRP radiation weighting factors Wr)
1797 * Worst possible geometry for the irradiation
1798 ELSEIF ((CSET(1:5).EQ.'EWT74').OR.(CSET(1:5).EQ.'ewt74')
1799 & .OR.(ITEST.EQ.3)) THEN
1800 IF (LFIRST(3).AND.(ITEST.EQ.0)) THEN
1806 IF (IPART.EQ.1) THEN
1808 IF (AEKIN.LT.AEBINN(1)) THEN
1810 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
1811 FLUSCW = WT74N(NBIN1N)
1812 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
1814 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
1815 & *DBLE(NBIN1N-47)+47.0D0
1818 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
1821 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
1825 IF (IBINHI.GT.NBIN1N) THEN
1826 FLUSCW = WT74N(NBIN1N)
1828 FLUSCW =( AEKIN-AEBINN(IBINLO))/
1829 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1830 & (AWT74N(IBINHI)-AWT74N(IBINLO))+AWT74N(IBINLO)
1831 FLUSCW = 10.0D0**FLUSCW
1837 IF ((IBINHI-IBINLO).EQ.1) GOTO 31
1838 KK = (IBINHI+IBINLO)/2
1839 IF (AEKIN.LE.AEBINN(KK)) THEN
1846 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
1847 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
1848 & (AWT74N(IBINHI)-AWT74N(IBINLO))+AWT74N(IBINLO)
1849 FLUSCW = 10.0D0**FLUSCW
1851 ELSEIF (IPART.EQ.2) THEN
1853 IF (AEKIN.LT.AE74P(1)) THEN
1855 ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
1856 FLUSCW = WT74P(NBIN1P)
1858 XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
1859 & *DBLE(NBIN1P-1)+1.0D0
1862 IF (AEKIN.LT.AE74P(IBINLO)) THEN
1865 ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
1869 IF (IBINHI.GT.NBIN1P) THEN
1870 FLUSCW = WT74P(NBIN1P)
1872 FLUSCW =( AEKIN-AE74P(IBINLO))/
1873 & (AE74P(IBINHI)-AE74P(IBINLO))*
1874 & (AWT74P(IBINHI)-AWT74P(IBINLO))+AWT74P(IBINLO)
1875 FLUSCW = 10.0D0**FLUSCW
1878 ELSEIF (IPART.EQ.3) THEN
1880 IF (AEKIN.LT.AEBINI(1)) THEN
1882 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1883 FLUSCW = WT74I(NBIN1I)
1885 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1886 & *DBLE(NBIN1I-1)+1.0D0
1889 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1892 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1896 IF (IBINHI.GT.NBIN1I) THEN
1897 FLUSCW = WT74I(NBIN1I)
1899 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1900 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1901 & (AWT74I(IBINHI)-AWT74I(IBINLO))+AWT74I(IBINLO)
1902 FLUSCW = 10.0D0**FLUSCW
1905 ELSEIF (IPART.EQ.4) THEN
1907 IF (AEKIN.LT.AEBINI(1)) THEN
1909 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
1910 FLUSCW = WT74J(NBIN1I)
1912 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
1913 & *DBLE(NBIN1I-1)+1.0D0
1916 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
1919 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
1923 IF (IBINHI.GT.NBIN1I) THEN
1924 FLUSCW = WT74J(NBIN1I)
1926 FLUSCW =( AEKIN-AEBINI(IBINLO))/
1927 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
1928 & (AWT74J(IBINHI)-AWT74J(IBINLO))+AWT74J(IBINLO)
1929 FLUSCW = 10.0D0**FLUSCW
1932 ELSEIF (IPART.EQ.7) THEN
1934 IF (AEKIN.LT.AE74G(1)) THEN
1936 ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
1937 FLUSCW = WT74G(NBIN1G)
1939 XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
1940 & *DBLE(NBIN1G-1)+1.0D0
1943 IF (AEKIN.LT.AE74G(IBINLO)) THEN
1946 ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
1950 IF (IBINHI.GT.NBIN1G) THEN
1951 FLUSCW = WT74G(NBIN1G)
1953 FLUSCW =( AEKIN-AE74G(IBINLO))/
1954 & (AE74G(IBINHI)-AE74G(IBINLO))*
1955 & (AWT74G(IBINHI)-AWT74G(IBINLO))+AWT74G(IBINLO)
1956 FLUSCW = 10.0D0**FLUSCW
1959 ELSEIF (IPART.EQ.8) THEN
1961 IF (AEKIN.LT.AE74E(1)) THEN
1963 ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
1964 FLUSCW = WT74E(NBIN1E)
1966 XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
1967 & *DBLE(NBIN1E-1)+1.0D0
1970 IF (AEKIN.LT.AE74E(IBINLO)) THEN
1973 ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
1977 IF (IBINHI.GT.NBIN1E) THEN
1978 FLUSCW = WT74E(NBIN1E)
1980 FLUSCW =( AEKIN-AE74E(IBINLO))/
1981 & (AE74E(IBINHI)-AE74E(IBINLO))*
1982 & (AWT74E(IBINHI)-AWT74E(IBINLO))+AWT74E(IBINLO)
1983 FLUSCW = 10.0D0**FLUSCW
1988 * Effective dose (Pelliccioni radiation weighting factors Wr)
1989 * Anterior-Posterior irradiation
1990 ELSEIF ((CSET(1:5).EQ.'EAPMP').OR.(CSET(1:5).EQ.'eapmp')
1991 & .OR.(ITEST.EQ.4)) THEN
1992 IF (LFIRST(4).AND.(ITEST.EQ.0)) THEN
1998 IF (IPART.EQ.1) THEN
2000 IF (AEKIN.LT.AEBINN(1)) THEN
2002 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
2003 FLUSCW = APPLN(NBIN1N)
2004 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
2006 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
2007 & *DBLE(NBIN1N-47)+47.0D0
2010 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
2013 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
2017 IF (IBINHI.GT.NBIN1N) THEN
2018 FLUSCW = APPLN(NBIN1N)
2020 FLUSCW =( AEKIN-AEBINN(IBINLO))/
2021 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2022 & (AAPPLN(IBINHI)-AAPPLN(IBINLO))+AAPPLN(IBINLO)
2023 FLUSCW = 10.0D0**FLUSCW
2029 IF ((IBINHI-IBINLO).EQ.1) GOTO 41
2030 KK = (IBINHI+IBINLO)/2
2031 IF (AEKIN.LE.AEBINN(KK)) THEN
2038 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
2039 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2040 & (AAPPLN(IBINHI)-AAPPLN(IBINLO))+AAPPLN(IBINLO)
2041 FLUSCW = 10.0D0**FLUSCW
2043 ELSEIF (IPART.EQ.2) THEN
2045 IF (AEKIN.LT.AEPLP(1)) THEN
2047 ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
2048 FLUSCW = APPLP(NBIN1P)
2050 XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
2051 & *DBLE(NBIN1P-1)+1.0D0
2054 IF (AEKIN.LT.AEPLP(IBINLO)) THEN
2057 ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
2061 IF (IBINHI.GT.NBIN1P) THEN
2062 FLUSCW = APPLP(NBIN1P)
2064 FLUSCW =( AEKIN-AEPLP(IBINLO))/
2065 & (AEPLP(IBINHI)-AEPLP(IBINLO))*
2066 & (AAPPLP(IBINHI)-AAPPLP(IBINLO))+AAPPLP(IBINLO)
2067 FLUSCW = 10.0D0**FLUSCW
2070 ELSEIF (IPART.EQ.3) THEN
2072 IF (AEKIN.LT.AEBINI(1)) THEN
2074 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2075 FLUSCW = APPLI(NBIN1I)
2077 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2078 & *DBLE(NBIN1I-1)+1.0D0
2081 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2084 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2088 IF (IBINHI.GT.NBIN1I) THEN
2089 FLUSCW = APPLI(NBIN1I)
2091 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2092 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2093 & (AAPPLI(IBINHI)-AAPPLI(IBINLO))+AAPPLI(IBINLO)
2094 FLUSCW = 10.0D0**FLUSCW
2097 ELSEIF (IPART.EQ.4) THEN
2099 IF (AEKIN.LT.AEBINI(1)) THEN
2101 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2102 FLUSCW = APPLJ(NBIN1I)
2104 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2105 & *DBLE(NBIN1I-1)+1.0D0
2108 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2111 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2115 IF (IBINHI.GT.NBIN1I) THEN
2116 FLUSCW = APPLJ(NBIN1I)
2118 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2119 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2120 & (AAPPLJ(IBINHI)-AAPPLJ(IBINLO))+AAPPLJ(IBINLO)
2121 FLUSCW = 10.0D0**FLUSCW
2124 ELSEIF (IPART.EQ.5) THEN
2126 IF (AEKIN.LT.AEBINM(1)) THEN
2128 ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
2129 FLUSCW = APPLM(NBIN1M)
2131 XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
2132 & *DBLE(NBIN1M-1)+1.0D0
2135 IF (AEKIN.LT.AEBINM(IBINLO)) THEN
2138 ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
2142 IF (IBINHI.GT.NBIN1M) THEN
2143 FLUSCW = APPLM(NBIN1M)
2145 FLUSCW =( AEKIN-AEBINM(IBINLO))/
2146 & (AEBINM(IBINHI)-AEBINM(IBINLO))*
2147 & (AAPPLM(IBINHI)-AAPPLM(IBINLO))+AAPPLM(IBINLO)
2148 FLUSCW = 10.0D0**FLUSCW
2154 * Effective dose (Pelliccioni radiation weighting factors Wr)
2155 * Rotational irradiation geometry
2156 ELSEIF ((CSET(1:5).EQ.'ERTMP').OR.(CSET(1:5).EQ.'ertmp')
2157 & .OR.(ITEST.EQ.5)) THEN
2158 IF (LFIRST(5).AND.(ITEST.EQ.0)) THEN
2164 IF (IPART.EQ.1) THEN
2166 IF (AEKIN.LT.AEBINN(1)) THEN
2168 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
2169 FLUSCW = RTPLN(NBIN1N)
2170 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
2172 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
2173 & *DBLE(NBIN1N-47)+47.0D0
2176 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
2179 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
2183 IF (IBINHI.GT.NBIN1N) THEN
2184 FLUSCW = RTPLN(NBIN1N)
2186 FLUSCW =( AEKIN-AEBINN(IBINLO))/
2187 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2188 & (ARTPLN(IBINHI)-ARTPLN(IBINLO))+ARTPLN(IBINLO)
2189 FLUSCW = 10.0D0**FLUSCW
2195 IF ((IBINHI-IBINLO).EQ.1) GOTO 51
2196 KK = (IBINHI+IBINLO)/2
2197 IF (AEKIN.LE.AEBINN(KK)) THEN
2204 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
2205 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2206 & (ARTPLN(IBINHI)-ARTPLN(IBINLO))+ARTPLN(IBINLO)
2207 FLUSCW = 10.0D0**FLUSCW
2209 ELSEIF (IPART.EQ.2) THEN
2211 IF (AEKIN.LT.AEPLP(1)) THEN
2213 ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
2214 FLUSCW = RTPLP(NBIN1P)
2216 XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
2217 & *DBLE(NBIN1P-1)+1.0D0
2220 IF (AEKIN.LT.AEPLP(IBINLO)) THEN
2223 ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
2227 IF (IBINHI.GT.NBIN1P) THEN
2228 FLUSCW = RTPLP(NBIN1P)
2230 FLUSCW =( AEKIN-AEPLP(IBINLO))/
2231 & (AEPLP(IBINHI)-AEPLP(IBINLO))*
2232 & (ARTPLP(IBINHI)-ARTPLP(IBINLO))+ARTPLP(IBINLO)
2233 FLUSCW = 10.0D0**FLUSCW
2236 ELSEIF (IPART.EQ.3) THEN
2238 IF (AEKIN.LT.AEBINI(1)) THEN
2240 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2241 FLUSCW = RTPLI(NBIN1I)
2243 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2244 & *DBLE(NBIN1I-1)+1.0D0
2247 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2250 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2254 IF (IBINHI.GT.NBIN1I) THEN
2255 FLUSCW = RTPLI(NBIN1I)
2257 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2258 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2259 & (ARTPLI(IBINHI)-ARTPLI(IBINLO))+ARTPLI(IBINLO)
2260 FLUSCW = 10.0D0**FLUSCW
2263 ELSEIF (IPART.EQ.4) THEN
2265 IF (AEKIN.LT.AEBINI(1)) THEN
2267 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2268 FLUSCW = RTPLJ(NBIN1I)
2270 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2271 & *DBLE(NBIN1I-1)+1.0D0
2274 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2277 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2281 IF (IBINHI.GT.NBIN1I) THEN
2282 FLUSCW = RTPLJ(NBIN1I)
2284 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2285 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2286 & (ARTPLJ(IBINHI)-ARTPLJ(IBINLO))+ARTPLJ(IBINLO)
2287 FLUSCW = 10.0D0**FLUSCW
2290 ELSEIF (IPART.EQ.5) THEN
2292 IF (AEKIN.LT.AEBINM(1)) THEN
2294 ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
2295 FLUSCW = RTPLM(NBIN1M)
2297 XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
2298 & *DBLE(NBIN1M-1)+1.0D0
2301 IF (AEKIN.LT.AEBINM(IBINLO)) THEN
2304 ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
2308 IF (IBINHI.GT.NBIN1M) THEN
2309 FLUSCW = RTPLM(NBIN1M)
2311 FLUSCW =( AEKIN-AEBINM(IBINLO))/
2312 & (AEBINM(IBINHI)-AEBINM(IBINLO))*
2313 & (ARTPLM(IBINHI)-ARTPLM(IBINLO))+ARTPLM(IBINLO)
2314 FLUSCW = 10.0D0**FLUSCW
2320 * Effective dose (Pelliccioni radiation weighting factors Wr)
2321 * Worst possible geometry for the irradiation
2322 ELSEIF ((CSET(1:5).EQ.'EWTMP').OR.(CSET(1:5).EQ.'ewtmp')
2323 & .OR.(ITEST.EQ.6)) THEN
2324 IF (LFIRST(6).AND.(ITEST.EQ.0)) THEN
2330 IF (IPART.EQ.1) THEN
2332 IF (AEKIN.LT.AEBINN(1)) THEN
2334 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
2335 FLUSCW = WTPLN(NBIN1N)
2336 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
2338 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
2339 & *DBLE(NBIN1N-47)+47.0D0
2342 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
2345 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
2349 IF (IBINHI.GT.NBIN1N) THEN
2350 FLUSCW = WTPLN(NBIN1N)
2352 FLUSCW =( AEKIN-AEBINN(IBINLO))/
2353 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2354 & (AWTPLN(IBINHI)-AWTPLN(IBINLO))+AWTPLN(IBINLO)
2355 FLUSCW = 10.0D0**FLUSCW
2361 IF ((IBINHI-IBINLO).EQ.1) GOTO 61
2362 KK = (IBINHI+IBINLO)/2
2363 IF (AEKIN.LE.AEBINN(KK)) THEN
2370 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
2371 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2372 & (AWTPLN(IBINHI)-AWTPLN(IBINLO))+AWTPLN(IBINLO)
2373 FLUSCW = 10.0D0**FLUSCW
2375 ELSEIF (IPART.EQ.2) THEN
2377 IF (AEKIN.LT.AEPLP(1)) THEN
2379 ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
2380 FLUSCW = WTPLP(NBIN1P)
2382 XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
2383 & *DBLE(NBIN1P-1)+1.0D0
2386 IF (AEKIN.LT.AEPLP(IBINLO)) THEN
2389 ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
2393 IF (IBINHI.GT.NBIN1P) THEN
2394 FLUSCW = WTPLP(NBIN1P)
2396 FLUSCW =( AEKIN-AEPLP(IBINLO))/
2397 & (AEPLP(IBINHI)-AEPLP(IBINLO))*
2398 & (AWTPLP(IBINHI)-AWTPLP(IBINLO))+AWTPLP(IBINLO)
2399 FLUSCW = 10.0D0**FLUSCW
2402 ELSEIF (IPART.EQ.3) THEN
2404 IF (AEKIN.LT.AEBINI(1)) THEN
2406 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2407 FLUSCW = WTPLI(NBIN1I)
2409 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2410 & *DBLE(NBIN1I-1)+1.0D0
2413 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2416 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2420 IF (IBINHI.GT.NBIN1I) THEN
2421 FLUSCW = WTPLI(NBIN1I)
2423 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2424 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2425 & (AWTPLI(IBINHI)-AWTPLI(IBINLO))+AWTPLI(IBINLO)
2426 FLUSCW = 10.0D0**FLUSCW
2429 ELSEIF (IPART.EQ.4) THEN
2431 IF (AEKIN.LT.AEBINI(1)) THEN
2433 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2434 FLUSCW = WTPLJ(NBIN1I)
2436 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2437 & *DBLE(NBIN1I-1)+1.0D0
2440 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2443 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2447 IF (IBINHI.GT.NBIN1I) THEN
2448 FLUSCW = WTPLJ(NBIN1I)
2450 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2451 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2452 & (AWTPLJ(IBINHI)-AWTPLJ(IBINLO))+AWTPLJ(IBINLO)
2453 FLUSCW = 10.0D0**FLUSCW
2456 ELSEIF (IPART.EQ.5) THEN
2458 IF (AEKIN.LT.AEBINM(1)) THEN
2460 ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
2461 FLUSCW = WTPLM(NBIN1M)
2463 XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
2464 & *DBLE(NBIN1M-1)+1.0D0
2467 IF (AEKIN.LT.AEBINM(IBINLO)) THEN
2470 ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
2474 IF (IBINHI.GT.NBIN1M) THEN
2475 FLUSCW = WTPLM(NBIN1M)
2477 FLUSCW =( AEKIN-AEBINM(IBINLO))/
2478 & (AEBINM(IBINHI)-AEBINM(IBINLO))*
2479 & (AWTPLM(IBINHI)-AWTPLM(IBINLO))+AWTPLM(IBINLO)
2480 FLUSCW = 10.0D0**FLUSCW
2486 * Ambient dose equivalent (from ICRP74 and Pelliccioni data)
2487 ELSEIF ((CSET(1:5).EQ.'AMB74').OR.(CSET(1:5).EQ.'amb74')
2488 & .OR.(ITEST.EQ.7)) THEN
2489 IF (LFIRST(7).AND.(ITEST.EQ.0)) THEN
2494 IF (IPART.EQ.1) THEN
2496 IF (AEKIN.LT.AEBINN(1)) THEN
2498 ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
2499 FLUSCW = AMBN(NBIN1N)
2500 ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
2502 XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
2503 & *DBLE(NBIN1N-47)+47.0D0
2506 IF (AEKIN.LT.AEBINN(IBINLO)) THEN
2509 ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
2513 IF (IBINHI.GT.NBIN1N) THEN
2514 FLUSCW = AMBN(NBIN1N)
2516 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
2517 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2518 & (AAMBN(IBINHI)-AAMBN(IBINLO))+AAMBN(IBINLO)
2519 FLUSCW = 10.0D0**FLUSCW
2525 IF ((IBINHI-IBINLO).EQ.1) GOTO 71
2526 KK = (IBINHI+IBINLO)/2
2527 IF (AEKIN.LE.AEBINN(KK)) THEN
2534 FLUSCW = ( AEKIN-AEBINN(IBINLO))/
2535 & (AEBINN(IBINHI)-AEBINN(IBINLO))*
2536 & (AAMBN(IBINHI)-AAMBN(IBINLO))+AAMBN(IBINLO)
2537 FLUSCW = 10.0D0**FLUSCW
2539 ELSEIF (IPART.EQ.2) THEN
2541 IF (AEKIN.LT.AEAMBP(1)) THEN
2543 ELSEIF (AEKIN.GT.AEAMBP(NBIN2P)) THEN
2544 FLUSCW = AMBP(NBIN2P)
2546 XBIN = (AEKIN-AEAMBP(1))/(AEAMBP(NBIN2P)-AEAMBP(1))
2547 & *DBLE(NBIN2P-1)+1.0D0
2550 IF (AEKIN.LT.AEAMBP(IBINLO)) THEN
2553 ELSEIF (AEKIN.GT.AEAMBP(IBINHI)) THEN
2557 IF (IBINHI.GT.NBIN2P) THEN
2558 FLUSCW = AMBP(NBIN2P)
2560 FLUSCW = ( AEKIN-AEAMBP(IBINLO))/
2561 & (AEAMBP(IBINHI)-AEAMBP(IBINLO))*
2562 & (AAMBP(IBINHI)-AAMBP(IBINLO))+AAMBP(IBINLO)
2563 FLUSCW = 10.0D0**FLUSCW
2566 ELSEIF (IPART.EQ.3) THEN
2568 IF (AEKIN.LT.AEBINI(1)) THEN
2570 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2571 FLUSCW = AMBI(NBIN1I)
2573 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2574 & *DBLE(NBIN1I-1)+1.0D0
2577 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2580 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2584 IF (IBINHI.GT.NBIN1I) THEN
2585 FLUSCW = AMBI(NBIN1I)
2587 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2588 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2589 & (AAMBI(IBINHI)-AAMBI(IBINLO))+AAMBI(IBINLO)
2590 FLUSCW = 10.0D0**FLUSCW
2593 ELSEIF (IPART.EQ.4) THEN
2595 IF (AEKIN.LT.AEBINI(1)) THEN
2597 ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
2598 FLUSCW = AMBJ(NBIN1I)
2600 XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
2601 & *DBLE(NBIN1I-1)+1.0D0
2604 IF (AEKIN.LT.AEBINI(IBINLO)) THEN
2607 ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
2611 IF (IBINHI.GT.NBIN1I) THEN
2612 FLUSCW = AMBJ(NBIN1I)
2614 FLUSCW =( AEKIN-AEBINI(IBINLO))/
2615 & (AEBINI(IBINHI)-AEBINI(IBINLO))*
2616 & (AAMBJ(IBINHI)-AAMBJ(IBINLO))+AAMBJ(IBINLO)
2617 FLUSCW = 10.0D0**FLUSCW
2620 ELSEIF (IPART.EQ.5) THEN
2622 IF (AEKIN.LT.AEBINM(1)) THEN
2624 ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
2625 FLUSCW = AMBM(NBIN1M)
2627 XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
2628 & *DBLE(NBIN1M-1)+1.0D0
2631 IF (AEKIN.LT.AEBINM(IBINLO)) THEN
2634 ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
2638 IF (IBINHI.GT.NBIN1M) THEN
2639 FLUSCW = AMBM(NBIN1M)
2641 FLUSCW =( AEKIN-AEBINM(IBINLO))/
2642 & (AEBINM(IBINHI)-AEBINM(IBINLO))*
2643 & (AAMBM(IBINHI)-AAMBM(IBINLO))+AAMBM(IBINLO)
2644 FLUSCW = 10.0D0**FLUSCW
2647 ELSEIF (IPART.EQ.7) THEN
2649 IF (AEKIN.LT.AE74G(1)) THEN
2651 ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
2652 FLUSCW = AMBG(NBIN1G)
2654 XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
2655 & *DBLE(NBIN1G-1)+1.0D0
2658 IF (AEKIN.LT.AE74G(IBINLO)) THEN
2661 ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
2665 IF (IBINHI.GT.NBIN1G) THEN
2666 FLUSCW = AMBG(NBIN1G)
2668 FLUSCW =( AEKIN-AE74G(IBINLO))/
2669 & (AE74G(IBINHI)-AE74G(IBINLO))*
2670 & (AAMBG(IBINHI)-AAMBG(IBINLO))+AAMBG(IBINLO)
2671 FLUSCW = 10.0D0**FLUSCW
2674 ELSEIF (IPART.EQ.8) THEN
2676 IF (AEKIN.LT.AEAMBE(1)) THEN
2678 ELSEIF (AEKIN.GT.AEAMBE(NBIN2E)) THEN
2679 FLUSCW = AMBE(NBIN2E)
2681 XBIN = (AEKIN-AEAMBE(1))/(AEAMBE(NBIN2E)-AEAMBE(1))
2682 & *DBLE(NBIN2E-1)+1.0D0
2685 IF (AEKIN.LT.AEAMBE(IBINLO)) THEN
2688 ELSEIF (AEKIN.GT.AEAMBE(IBINHI)) THEN
2692 IF (IBINHI.GT.NBIN2E) THEN
2693 FLUSCW = AMBE(NBIN2E)
2695 FLUSCW =( AEKIN-AEAMBE(IBINLO))/
2696 & (AEAMBE(IBINHI)-AEAMBE(IBINLO))*
2697 & (AAMBE(IBINHI)-AAMBE(IBINLO))+AAMBE(IBINLO)
2698 FLUSCW = 10.0D0**FLUSCW
2704 * Ambient dose equivalent (old GRS conversion factors)
2705 ELSEIF ((CSET(1:5).EQ.'AMBGS').OR.(CSET(1:5).EQ.'ambgs')
2706 & .OR.(ITEST.EQ.8)) THEN
2707 IF (LFIRST(8).AND.(ITEST.EQ.0)) THEN
2713 IF (IPART.EQ.1) THEN
2714 XDUM = (11.0D0+AEKIN)*10.0D0
2715 KENERG = INT(XDUM)+1
2716 IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
2717 FLUSCW = H10N(KENERG)
2718 ELSEIF (IPART.EQ.2) THEN
2719 XDUM = (11.0D0+AEKIN)*10.0D0
2720 KENERG = INT(XDUM)+1
2721 IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
2722 FLUSCW = H10PR(KENERG)
2723 ELSEIF ((IPART.EQ.3).OR.(IPART.EQ.4)) THEN
2724 XDUM = (11.0D0+AEKIN)*10.0D0
2725 KENERG = INT(XDUM)+1
2726 IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
2727 FLUSCW = H10PI(KENERG)
2728 ELSEIF (IPART.EQ.5) THEN
2729 IF (AEKIN.LT.AEAMBM(1)) THEN
2731 ELSEIF (AEKIN.GT.AEAMBM(NBIN2M)) THEN
2732 FLUSCW = H10MU(NBIN2M)
2737 IF ((IBINHI-IBINLO).EQ.1) GOTO 81
2738 KK = (IBINHI+IBINLO)/2
2739 IF (AEKIN.LE.AEAMBM(KK)) THEN
2746 FLUSCW = ( AEKIN-AEAMBM(IBINLO))/
2747 & (AEAMBM(IBINHI)-AEAMBM(IBINLO))*
2748 & (AH10MU(IBINHI)-AH10MU(IBINLO))+AH10MU(IBINLO)
2749 FLUSCW = 10.0D0**FLUSCW
2753 * this point should never be reached since the check for a valid
2754 * sdum has already been done above
2762 1000 FORMAT(1X,'FLUSCW: direct conversion of fluence to effective ',
2763 & 'dose requested with')
2764 1001 FORMAT(1X,'FLUSCW: direct conversion of fluence to ambient ',
2765 & 'dose eq. requested with')
2766 1002 FORMAT(12X,'- ICRP radiation weighting factors Wr')
2767 1003 FORMAT(12X,'- Pelliccioni radiation weighting factors Wr')
2768 1004 FORMAT(12X,'- anterior-posterior irradiation')
2769 1005 FORMAT(12X,'- rotational irradiation geometry')
2770 1006 FORMAT(12X,'- worst possible geometry for the irradiation')
2771 1007 FORMAT(12X,'- GRS conversion factors')
2772 1008 FORMAT(12X,'- ICRP74 and Pelliccioni conversion factors')
git://git.uio.no / u / mrichter / AliRoot.git / blob