--- /dev/null
+************************************************************************
+* Energy dependent factors for the conversion of fluence to *
+* effective dose and ambient dose equivalent for neutrons, protons, *
+* charged pions, muons, photons and electrons. *
+* *
+* The following sets are available: *
+* (for photons and electrons see note below) *
+* *
+* 1-3: Effective dose from ICRP74 and Pelliccioni data *
+* calculated with ICRP radiation weighting factors Wr *
+* 1: Anterior-Posterior irradiation *
+* 2: Rotational irradiation geometry *
+* 3: WORST possible geometry for the irradiation *
+* *
+* 4-6: Effective dose from ICRP74 and Pelliccioni data *
+* calculated with the Pelliccioni radiation weighting *
+* factors Wr *
+* 4: Anterior-Posterior irradiation *
+* 5: Rotational irradiation geometry *
+* 6: WORST possible geometry for the irradiation *
+* *
+* 7: Ambient dose equivalent from ICRP74 and Pelliccioni data *
+* *
+* 8: Ambient dose equivalent with old "GRS"-conversion factors *
+* *
+* The different sets are invoked by EXTRAWEI (WHAT(3)=1) and the *
+* first five characters of the SDUM input parameter of the USRBIN, *
+* USRTRACK or USRBDX cards (case-insensitive) : *
+* *
+* SDUM(1:5) = *
+* 'EAP74' = 1 'ERT74' = 2 'EWT74' = 3 *
+* 'EAPMP' = 4 'ERTMP' = 5 'EWTMP' = 6 *
+* 'AMB74' = 7 'AMBGS' = 8 *
+* *
+* The default (any other SDUM) is unit weight factor. *
+* *
+* Conversion factor sets for hadrons (n,p,pi) and muons can be *
+* enabled separately from those for photons and electrons by *
+* setting the 6th character of SDUM: *
+* *
+* SDUM(6:6) = 1 only hadrons and muons (zero factor returned for *
+* electrons and photons) *
+* = 2 only electrons/positrons and photons (zero factor *
+* returned for hadrons and muons) *
+* otherwise: all particles considered (default) *
+* *
+* Note: Photons and electrons: *
+* Only sets 1-3 and 7 are implemented for photons and electrons*
+* If sets 4-6 are requested sets 1-3 are used instead. *
+* For set 8 zero values are returned. *
+* *
+* Note 2: Stand-alone use: *
+* The function can be used "stand-alone" (i.e. independent *
+* from FLUKA). In this case different conversion factor sets *
+* are invoked by the first parameter (IIJ) *
+* *
+* IIJ = -(is*10000+id) where *
+* *
+* is = index of parameter set *
+* id = particle identity *
+* *
+* (example: -80008 = neutron with GRS-conversion factor set) *
+************************************************************************
+* *
+*=== fluscw ===========================================================*
+* *
+ DOUBLE PRECISION FUNCTION FLUSCW
+ & (IIJ,PLA,TXX,TYY,TZZ,WEE,XX,YY,ZZ,NREG,IOLREG,LLO,NSURF)
+
+ INCLUDE '(DBLPRC)'
+ INCLUDE '(DIMPAR)'
+ INCLUDE '(IOUNIT)'
+ SAVE
+
+ INCLUDE '(PAPROP)'
+ INCLUDE '(USRBIN)'
+ INCLUDE '(USRBDX)'
+ INCLUDE '(USRTRC)'
+ INCLUDE '(SCOHLP)'
+
+ LOGICAL LFIRST
+ CHARACTER*10 CSET
+ DIMENSION LFIRST(0:8)
+
+ PARAMETER (NBIN1N=81,NBIN2N=150,
+ & NBIN1P=64,NBIN2P=61,
+ & NBIN1I=71,
+ & NBIN1M=71,NBIN2M=16,
+ & NBIN1G=81,
+ & NBIN1E=41,NBIN2E=81)
+
+ DIMENSION ITT(40)
+* 1 - neutron, 2 - proton, 3 - pion+, 4 - pion-, 5 - muon+, 5 - muon-,
+* 7 - photon, 8 - electron, 99 - the rest
+ DATA ITT/ 2, 2, 8, 8,99,99, 7, 1, 1, 5,
+ & 5,99, 3, 4, 3, 4, 1, 1,99, 4,
+ & 3,99,99,99,99,99,99,99,99,99,
+ & 99,99,99,99,99,99,99,99,99,99/
+*
+*-----------------------------------------------------------------------
+*
+* Effective dose for neutrons from ICRP74 and Pelliccioni data
+* calculated with ICRP radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EBINN(NBIN1N),AEBINN(NBIN1N)
+ DATA EBINN /
+ 1 1.000D-12, 1.000D-11, 2.530D-11, 1.000D-10, 2.000D-10,
+ 2 5.000D-10, 1.000D-09, 2.000D-09, 5.000D-09, 1.000D-08,
+ 3 2.000D-08, 5.000D-08, 1.000D-07, 2.000D-07, 5.000D-07,
+ 4 1.000D-06, 2.000D-06, 5.000D-06, 1.000D-05, 2.000D-05,
+ 5 3.000D-05, 5.000D-05, 7.000D-05, 1.000D-04, 1.500D-04,
+ 6 2.000D-04, 3.000D-04, 5.000D-04, 7.000D-04, 9.000D-04,
+ 7 1.000D-03, 1.200D-03, 2.000D-03, 3.000D-03, 4.000D-03,
+ 8 5.000D-03, 6.000D-03, 7.000D-03, 8.000D-03, 9.000D-03,
+ 9 1.000D-02, 1.200D-02, 1.400D-02, 1.500D-02, 1.600D-02,
+ * 1.800D-02, 2.000D-02, 3.000D-02, 5.000D-02, 7.500D-02,
+ 1 1.000D-01, 1.500D-01, 2.000D-01, 3.000D-01, 5.000D-01,
+ 2 7.000D-01, 1.000D+00, 1.500D+00, 2.000D+00, 3.000D+00,
+ 3 5.000D+00, 7.000D+00, 1.000D+01, 1.500D+01, 2.000D+01,
+ 4 3.000D+01, 5.000D+01, 7.000D+01, 1.000D+02, 1.500D+02,
+ 5 2.000D+02, 3.000D+02, 5.000D+02, 7.000D+02, 1.000D+03,
+ 6 1.500D+03, 2.000D+03, 3.000D+03, 5.000D+03, 7.000D+03,
+ 7 1.000D+04/
+ DIMENSION AP74N(NBIN1N),RT74N(NBIN1N),WT74N(NBIN1N),
+ & AAP74N(NBIN1N),ART74N(NBIN1N),AWT74N(NBIN1N)
+* Anterior-Posterior irradiation
+ DATA AP74N /
+ 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
+ 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
+ 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
+ 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
+ 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
+ 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
+ 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
+ 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
+ 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
+ * 486.0D0, 480.0D0, 458.0D0, 437.0D0, 429.0D0,
+ 1 429.0D0, 431.0D0, 447.0D0, 491.0D0, 581.0D0,
+ 2 663.0D0, 767.0D0, 905.0D0, 1016.0D0, 1191.0D0,
+ 3 1443.0D0, 1628.0D0, 1840.0D0, 2097.0D0, 2291.0D0,
+ 4 2581.0D0, 2978.0D0, 3263.0D0, 3590.0D0, 3999.0D0,
+ 5 4315.0D0, 4801.0D0, 5484.0D0, 5980.0D0, 6550.0D0,
+ 6 7255.0D0, 7794.0D0, 8614.0D0, 9751.0D0, 10569.0D0,
+ 7 11500.0D0/
+* Rotational irradiation geometry
+ DATA RT74N /
+ 1 2.99D0, 3.72D0, 4.40D0, 5.75D0, 6.43D0,
+ 2 7.27D0, 7.84D0, 8.31D0, 8.72D0, 8.90D0,
+ 3 8.92D0, 8.82D0, 8.69D0, 8.56D0, 8.40D0,
+ 4 8.34D0, 8.39D0, 9.06D0, 10.6D0, 13.8D0,
+ 5 16.9D0, 22.7D0, 27.8D0, 34.8D0, 45.4D0,
+ 6 54.8D0, 71.6D0, 99.4D0, 123.0D0, 144.0D0,
+ 7 154.0D0, 173.0D0, 234.0D0, 283.0D0, 315.0D0,
+ 8 335.0D0, 348.0D0, 358.0D0, 366.0D0, 373.0D0,
+ 9 378.0D0, 385.0D0, 390.0D0, 391.0D0, 393.0D0,
+ * 394.0D0, 395.0D0, 395.0D0, 404.0D0, 422.0D0,
+ 1 443.0D0, 496.0D0, 535.0D0, 594.0D0, 681.0D0,
+ 2 754.0D0, 854.0D0, 1008.0D0, 1149.0D0, 1396.0D0,
+ 3 1792.0D0, 2102.0D0, 2460.0D0, 2884.0D0, 3192.0D0,
+ 4 3639.0D0, 4238.0D0, 4669.0D0, 5180.0D0, 5855.0D0,
+ 5 6406.0D0, 7296.0D0, 8634.0D0, 9664.0D0, 10900.0D0,
+ 6 12505.0D0, 13788.0D0, 15828.0D0, 18844.0D0, 21145.0D0,
+ 7 23900.0D0/
+* WORST possible geometry for the irradiation
+ DATA WT74N /
+ 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
+ 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
+ 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
+ 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
+ 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
+ 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
+ 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
+ 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
+ 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
+ * 486.0D0, 480.0D0, 458.0D0, 444.0D0, 459.0D0,
+ 1 477.0D0, 523.0D0, 559.0D0, 616.0D0, 707.0D0,
+ 2 785.0D0, 892.0D0, 1056.0D0, 1205.0D0, 1468.0D0,
+ 3 1898.0D0, 2244.0D0, 2660.0D0, 3184.0D0, 3586.0D0,
+ 4 4201.0D0, 5065.0D0, 5702.0D0, 6450.0D0, 7415.0D0,
+ 5 8184.0D0, 9409.0D0, 11224.0D0, 12618.0D0, 14300.0D0,
+ 6 16512.0D0, 18304.0D0, 21196.0D0, 25562.0D0, 28960.0D0,
+ 7 33100.0D0/
+*
+* Effective dose for protons from ICRP74 and Pelliccioni data
+* calculated with ICRP radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION E74P(NBIN1P),AE74P(NBIN1P)
+ DATA E74P /
+ 1 5.00D-03, 6.50D-03, 8.00D-03, 1.00D-02, 1.25D-02,
+ 2 1.50D-02, 2.00D-02, 2.50D-02, 3.00D-02, 4.00D-02,
+ 3 5.00D-02, 6.50D-02, 8.00D-02, 1.00D-01, 1.25D-01,
+ 4 1.50D-01, 2.00D-01, 2.50D-01, 3.00D-01, 4.00D-01,
+ 5 5.00D-01, 6.50D-01, 8.00D-01, 1.00D+00, 1.25D+00,
+ 6 1.50D+00, 2.00D+00, 2.50D+00, 3.00D+00, 4.00D+00,
+ 7 5.00D+00, 6.50D+00, 8.00D+00, 1.00D+01, 1.25D+01,
+ 8 1.50D+01, 2.00D+01, 2.50D+01, 3.00D+01, 4.00D+01,
+ 9 5.00D+01, 6.50D+01, 8.00D+01, 1.00D+02, 1.25D+02,
+ * 1.50D+02, 2.00D+02, 2.50D+02, 3.00D+02, 4.00D+02,
+ 1 5.00D+02, 6.50D+02, 8.00D+02, 1.00D+03, 1.25D+03,
+ 2 1.50D+03, 2.00D+03, 2.50D+03, 3.00D+03, 4.00D+03,
+ 3 5.00D+03, 6.50D+03, 8.00D+03, 1.00D+04/
+ DIMENSION AP74P(NBIN1P),RT74P(NBIN1P),WT74P(NBIN1P),
+ & AAP74P(NBIN1P),ART74P(NBIN1P),AWT74P(NBIN1P)
+* Anterior-Posterior irradiation
+ DATA AP74P /
+ 1 57.5D0, 65.4D0, 74.3D0, 89.5D0, 116.4D0,
+ 2 155.2D0, 289.0D0, 546.5D0, 982.1D0, 2540.0D0,
+ 3 4810.0D0, 6969.6D0, 7336.4D0, 6820.0D0, 6080.5D0,
+ 4 5467.0D0, 4570.0D0, 3977.7D0, 3573.9D0, 3089.0D0,
+ 5 2840.0D0, 2673.7D0, 2597.5D0, 2530.0D0, 2455.9D0,
+ 6 2410.1D0, 2420.0D0, 2534.3D0, 2688.1D0, 2996.8D0,
+ 7 3240.0D0, 3464.8D0, 3601.1D0, 3730.0D0, 3862.6D0,
+ 8 3975.3D0, 4155.7D0, 4292.5D0, 4398.5D0, 4547.2D0,
+ 9 4640.0D0, 4722.1D0, 4784.0D0, 4870.0D0, 4995.2D0,
+ * 5128.6D0, 5394.4D0, 5645.8D0, 5879.5D0, 6298.2D0,
+ 1 6662.4D0, 7129.9D0, 7525.5D0, 7970.0D0, 8428.5D0,
+ 2 8812.8D0, 9438.1D0, 9940.3D0, 10362.6D0, 11053.4D0,
+ 3 11611.9D0, 12297.0D0, 12863.8D0, 13500.0D0/
+* Rotational irradiation geometry
+ DATA RT74P /
+ 1 48.2D0, 54.5D0, 61.1D0, 71.4D0, 87.6D0,
+ 2 108.5D0, 169.7D0, 266.8D0, 409.7D0, 875.8D0,
+ 3 1582.0D0, 2499.1D0, 3087.1D0, 3642.5D0, 4289.0D0,
+ 4 4804.9D0, 5160.0D0, 4771.6D0, 4200.7D0, 3304.3D0,
+ 5 2832.5D0, 2593.8D0, 2546.2D0, 2540.0D0, 2518.5D0,
+ 6 2501.2D0, 2527.5D0, 2630.4D0, 2767.2D0, 3059.7D0,
+ 7 3330.0D0, 3664.4D0, 3928.8D0, 4200.0D0, 4446.4D0,
+ 8 4629.0D0, 4888.9D0, 5074.6D0, 5222.2D0, 5461.2D0,
+ 9 5665.0D0, 5941.3D0, 6192.5D0, 6497.5D0, 6842.1D0,
+ * 7154.8D0, 7708.2D0, 8190.8D0, 8621.5D0, 9370.8D0,
+ 1 10013.3D0, 10839.7D0, 11549.4D0, 12367.5D0, 13243.7D0,
+ 2 14005.1D0, 15295.5D0, 16377.1D0, 17317.0D0, 18910.3D0,
+ 3 20245.9D0, 21937.1D0, 23374.8D0, 25025.0D0/
+* WORST possible geometry for the irradiation
+ DATA WT74P /
+ 1 59.2D0, 66.8D0, 75.5D0, 90.6D0, 117.3D0,
+ 2 155.8D0, 289.0D0, 545.6D0, 980.1D0, 2540.0D0,
+ 3 4810.0D0, 6860.6D0, 7172.3D0, 6820.0D0, 6501.3D0,
+ 4 6255.9D0, 5640.0D0, 4889.3D0, 4237.9D0, 3379.9D0,
+ 5 2950.0D0, 2718.2D0, 2661.4D0, 2650.0D0, 2639.8D0,
+ 6 2636.5D0, 2670.0D0, 2751.5D0, 2859.1D0, 3109.6D0,
+ 7 3380.0D0, 3787.8D0, 4155.1D0, 4550.0D0, 4897.3D0,
+ 8 5139.5D0, 5459.3D0, 5672.7D0, 5839.6D0, 6122.3D0,
+ 9 6390.0D0, 6794.6D0, 7181.0D0, 7650.0D0, 8163.8D0,
+ * 8616.8D0, 9397.9D0, 10064.3D0, 10651.5D0, 11662.8D0,
+ 1 12524.9D0, 13633.2D0, 14589.0D0, 15700.0D0, 16904.4D0,
+ 2 17963.1D0, 19781.3D0, 21326.4D0, 22683.3D0, 25011.0D0,
+ 3 26986.3D0, 29515.1D0, 31686.5D0, 34200.0D0/
+*
+* Effective dose for charged pions from ICRP74 and Pelliccioni data
+* calculated with ICRP radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EBINI(NBIN1I),AEBINI(NBIN1I)
+ DATA EBINI/
+ 1 1.00D-03, 1.25D-03, 1.50D-03, 2.00D-03, 2.50D-03,
+ 2 3.00D-03, 4.00D-03, 5.00D-03, 6.50D-03, 8.00D-03,
+ 3 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
+ 4 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
+ 5 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
+ 6 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
+ 7 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
+ 8 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
+ 9 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
+ * 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
+ 1 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
+ 2 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
+ 3 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
+ 4 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
+ 5 1.00D+04/
+ DIMENSION AP74I(NBIN1I),RT74I(NBIN1I),WT74I(NBIN1I),
+ & AP74J(NBIN1I),RT74J(NBIN1I),WT74J(NBIN1I),
+ & AAP74I(NBIN1I),ART74I(NBIN1I),AWT74I(NBIN1I),
+ & AAP74J(NBIN1I),ART74J(NBIN1I),AWT74J(NBIN1I)
+* Anterior-Posterior irradiation
+* Positive pions
+ DATA AP74I /
+ 1 71.3D0, 76.3D0, 80.9D0, 89.3D0, 97.4D0,
+ 2 105.6D0, 122.3D0, 140.1D0, 168.8D0, 200.6D0,
+ 3 248.0D0, 315.4D0, 389.4D0, 547.5D0, 707.6D0,
+ 4 859.7D0, 1117.7D0, 1300.0D0, 1445.7D0, 1498.3D0,
+ 5 1498.6D0, 1454.3D0, 1396.7D0, 1290.0D0, 1209.1D0,
+ 6 1148.3D0, 1065.0D0, 1013.1D0, 968.2D0, 946.4D0,
+ 7 938.0D0, 946.1D0, 962.9D0, 1002.2D0, 1038.6D0,
+ 8 1068.4D0, 1107.9D0, 1125.0D0, 1124.7D0, 1115.3D0,
+ 9 1105.0D0, 1101.1D0, 1103.4D0, 1116.2D0, 1133.3D0,
+ * 1151.6D0, 1188.2D0, 1222.6D0, 1269.5D0, 1311.1D0,
+ 1 1360.0D0, 1412.9D0, 1459.0D0, 1537.5D0, 1603.3D0,
+ 2 1660.3D0, 1756.4D0, 1836.3D0, 1936.3D0, 2020.2D0,
+ 3 2115.0D0, 2214.8D0, 2300.1D0, 2442.1D0, 2558.6D0,
+ 4 2658.2D0, 2823.7D0, 2959.4D0, 3127.6D0, 3267.5D0,
+ 5 3425.0D0/
+* Negative pions
+ DATA AP74J /
+ 1 349.0D0, 359.7D0, 369.2D0, 387.2D0, 404.9D0,
+ 2 422.9D0, 460.0D0, 499.0D0, 560.9D0, 626.8D0,
+ 3 721.0D0, 847.1D0, 975.3D0, 1218.2D0, 1427.1D0,
+ 4 1594.6D0, 1806.2D0, 1880.0D0, 1830.5D0, 1712.3D0,
+ 5 1546.4D0, 1373.0D0, 1243.0D0, 1080.0D0, 996.1D0,
+ 6 950.4D0, 909.5D0, 898.1D0, 902.1D0, 915.7D0,
+ 7 938.0D0, 965.9D0, 991.4D0, 1033.7D0, 1065.3D0,
+ 8 1088.4D0, 1115.5D0, 1125.0D0, 1121.7D0, 1113.2D0,
+ 9 1105.0D0, 1102.9D0, 1106.1D0, 1119.7D0, 1137.0D0,
+ * 1155.2D0, 1191.3D0, 1225.1D0, 1271.0D0, 1311.9D0,
+ 1 1360.0D0, 1412.3D0, 1458.1D0, 1536.2D0, 1601.9D0,
+ 2 1658.9D0, 1755.2D0, 1835.3D0, 1935.6D0, 2019.8D0,
+ 3 2115.0D0, 2215.1D0, 2300.6D0, 2442.7D0, 2559.4D0,
+ 4 2659.0D0, 2824.4D0, 2960.0D0, 3128.0D0, 3267.8D0,
+ 5 3425.0D0/
+* Rotational irradiation geometry
+* Positive pions
+ DATA RT74I /
+ 1 37.0D0, 39.1D0, 40.9D0, 44.5D0, 48.0D0,
+ 2 51.7D0, 59.5D0, 68.0D0, 82.3D0, 98.6D0,
+ 3 123.5D0, 160.1D0, 201.6D0, 294.8D0, 395.2D0,
+ 4 496.8D0, 686.5D0, 842.0D0, 1002.6D0, 1099.3D0,
+ 5 1167.8D0, 1200.3D0, 1203.6D0, 1177.5D0, 1141.2D0,
+ 6 1106.0D0, 1047.8D0, 1005.6D0, 964.7D0, 942.3D0,
+ 7 931.0D0, 935.2D0, 949.0D0, 986.7D0, 1027.2D0,
+ 8 1065.8D0, 1133.5D0, 1188.0D0, 1249.8D0, 1295.9D0,
+ 9 1342.0D0, 1384.2D0, 1416.2D0, 1463.7D0, 1499.2D0,
+ * 1528.4D0, 1576.7D0, 1618.0D0, 1673.5D0, 1724.6D0,
+ 1 1789.0D0, 1865.3D0, 1936.9D0, 2067.0D0, 2182.4D0,
+ 2 2286.3D0, 2467.4D0, 2622.2D0, 2819.2D0, 2985.8D0,
+ 3 3174.0D0, 3370.4D0, 3537.1D0, 3811.8D0, 4035.6D0,
+ 4 4225.8D0, 4540.7D0, 4798.2D0, 5117.5D0, 5383.9D0,
+ 5 5685.0D0/
+* Negative pions
+ DATA RT74J /
+ 1 194.8D0, 196.1D0, 197.6D0, 201.6D0, 206.9D0,
+ 2 213.3D0, 228.3D0, 245.9D0, 275.9D0, 309.7D0,
+ 3 360.2D0, 431.1D0, 506.5D0, 659.1D0, 802.8D0,
+ 4 930.0D0, 1123.6D0, 1236.5D0, 1292.2D0, 1280.7D0,
+ 5 1229.6D0, 1156.7D0, 1093.0D0, 1004.5D0, 956.1D0,
+ 6 929.0D0, 904.6D0, 898.4D0, 902.7D0, 913.5D0,
+ 7 931.0D0, 953.5D0, 975.4D0, 1015.8D0, 1052.0D0,
+ 8 1084.5D0, 1140.7D0, 1188.0D0, 1246.7D0, 1293.6D0,
+ 9 1342.0D0, 1386.3D0, 1419.5D0, 1468.1D0, 1503.8D0,
+ * 1532.9D0, 1580.6D0, 1621.1D0, 1675.4D0, 1725.5D0,
+ 1 1789.0D0, 1864.5D0, 1935.7D0, 2065.3D0, 2180.6D0,
+ 2 2284.4D0, 2465.8D0, 2620.8D0, 2818.3D0, 2985.3D0,
+ 3 3174.0D0, 3370.8D0, 3537.7D0, 3812.8D0, 4036.7D0,
+ 4 4227.0D0, 4541.7D0, 4799.1D0, 5118.1D0, 5384.2D0,
+ 5 5685.0D0/
+* WORST possible geometry for the irradiation
+* Positive pions
+ DATA WT74I /
+ 1 71.3D0, 76.3D0, 80.8D0, 89.2D0, 97.3D0,
+ 2 105.5D0, 122.2D0, 139.9D0, 168.7D0, 200.5D0,
+ 3 248.0D0, 315.6D0, 389.9D0, 548.5D0, 709.1D0,
+ 4 861.5D0, 1119.2D0, 1300.0D0, 1442.3D0, 1491.8D0,
+ 5 1489.9D0, 1445.3D0, 1389.6D0, 1290.0D0, 1217.9D0,
+ 6 1165.1D0, 1093.7D0, 1049.0D0, 1008.9D0, 987.0D0,
+ 7 974.0D0, 972.3D0, 979.2D0, 1004.7D0, 1037.1D0,
+ 8 1072.1D0, 1144.0D0, 1215.0D0, 1315.6D0, 1402.7D0,
+ 9 1495.0D0, 1577.3D0, 1635.9D0, 1714.7D0, 1766.7D0,
+ * 1805.1D0, 1862.6D0, 1908.6D0, 1969.8D0, 2028.0D0,
+ 1 2105.0D0, 2201.5D0, 2295.6D0, 2472.3D0, 2633.3D0,
+ 2 2780.5D0, 3041.0D0, 3265.9D0, 3553.6D0, 3796.8D0,
+ 3 4070.0D0, 4351.8D0, 4588.0D0, 4972.1D0, 5280.3D0,
+ 4 5539.3D0, 5962.8D0, 6305.5D0, 6726.6D0, 7076.1D0,
+ 5 7470.0D0/
+* Negative pions
+ DATA WT74J /
+ 1 349.0D0, 359.6D0, 369.1D0, 386.9D0, 404.5D0,
+ 2 422.4D0, 459.5D0, 498.5D0, 560.4D0, 626.5D0,
+ 3 721.0D0, 847.7D0, 976.5D0, 1220.5D0, 1430.2D0,
+ 4 1597.9D0, 1808.6D0, 1880.0D0, 1826.1D0, 1704.9D0,
+ 5 1537.4D0, 1364.6D0, 1236.7D0, 1080.0D0, 1003.4D0,
+ 6 964.3D0, 934.0D0, 929.9D0, 939.9D0, 955.0D0,
+ 7 974.0D0, 992.6D0, 1008.2D0, 1036.3D0, 1063.9D0,
+ 8 1092.2D0, 1151.7D0, 1215.0D0, 1312.1D0, 1400.0D0,
+ 9 1495.0D0, 1579.7D0, 1639.9D0, 1720.2D0, 1772.5D0,
+ * 1810.8D0, 1867.4D0, 1912.4D0, 1972.1D0, 2029.2D0,
+ 1 2105.0D0, 2200.5D0, 2294.1D0, 2470.1D0, 2631.0D0,
+ 2 2778.2D0, 3038.8D0, 3264.1D0, 3552.4D0, 3796.2D0,
+ 3 4070.0D0, 4352.4D0, 4589.0D0, 4973.4D0, 5281.8D0,
+ 4 5540.9D0, 5964.3D0, 6306.8D0, 6727.5D0, 7076.6D0,
+ 5 7470.0D0/
+*
+* Effective dose for photons from ICRP74 and Pelliccioni data
+* calculated with ICRP radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION E74G(NBIN1G),AE74G(NBIN1G)
+ DATA E74G /
+ 1 0.10000D-05, 0.12589D-05, 0.15849D-05, 0.19953D-05, 0.25119D-05,
+ 2 0.31623D-05, 0.39811D-05, 0.50119D-05, 0.63096D-05, 0.79433D-05,
+ 3 0.10000D-04, 0.12589D-04, 0.15849D-04, 0.19953D-04, 0.25119D-04,
+ 4 0.31623D-04, 0.39811D-04, 0.50119D-04, 0.63096D-04, 0.79433D-04,
+ 5 0.10000D-03, 0.12589D-03, 0.15849D-03, 0.19953D-03, 0.25119D-03,
+ 6 0.31623D-03, 0.39811D-03, 0.50119D-03, 0.63096D-03, 0.79433D-03,
+ 7 0.10000D-02, 0.12589D-02, 0.15849D-02, 0.19953D-02, 0.25119D-02,
+ 8 0.31623D-02, 0.39811D-02, 0.50119D-02, 0.63096D-02, 0.79433D-02,
+ 9 0.10000D-01, 0.12589D-01, 0.15849D-01, 0.19953D-01, 0.25119D-01,
+ * 0.31623D-01, 0.39811D-01, 0.50119D-01, 0.63096D-01, 0.79433D-01,
+ 1 0.10000D+00, 0.12589D+00, 0.15849D+00, 0.19953D+00, 0.25119D+00,
+ 2 0.31623D+00, 0.39811D+00, 0.50119D+00, 0.63096D+00, 0.79433D+00,
+ 3 0.10000D+01, 0.12589D+01, 0.15849D+01, 0.19953D+01, 0.25119D+01,
+ 4 0.31623D+01, 0.39811D+01, 0.50119D+01, 0.63096D+01, 0.79433D+01,
+ 5 0.10000D+02, 0.12589D+02, 0.15849D+02, 0.19953D+02, 0.25119D+02,
+ 6 0.31623D+02, 0.39811D+02, 0.50119D+02, 0.63096D+02, 0.79433D+02,
+ 7 0.10000D+03/
+ DIMENSION AP74G(NBIN1G),RT74G(NBIN1G),WT74G(NBIN1G),
+ & AAP74G(NBIN1G),ART74G(NBIN1G),AWT74G(NBIN1G)
+* Anterior-Posterior irradiation
+ DATA AP74G /
+ 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
+ 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18658D+02, 0.11621D+02,
+ 3 0.71617D+01, 0.43328D+01, 0.26230D+01, 0.15912D+01, 0.97875D+00,
+ 4 0.62511D+00, 0.45300D+00, 0.41541D+00, 0.45697D+00, 0.50270D+00,
+ 5 0.55300D+00, 0.68970D+00, 0.86020D+00, 0.10728D+01, 0.13381D+01,
+ 6 0.16688D+01, 0.20814D+01, 0.25952D+01, 0.31566D+01, 0.38395D+01,
+ 7 0.46700D+01, 0.55604D+01, 0.65903D+01, 0.76977D+01, 0.89735D+01,
+ 8 0.10487D+02, 0.12358D+02, 0.13617D+02, 0.15456D+02, 0.17725D+02,
+ 9 0.21000D+02, 0.24813D+02, 0.29318D+02, 0.34640D+02, 0.41642D+02,
+ * 0.49614D+02, 0.57329D+02, 0.61958D+02, 0.67831D+02, 0.74261D+02,
+ 1 0.81300D+02, 0.85676D+02, 0.90288D+02, 0.95149D+02, 0.98906D+02,
+ 2 0.10280D+03, 0.10684D+03, 0.11104D+03, 0.11458D+03, 0.11823D+03,
+ 3 0.12200D+03, 0.12364D+03, 0.12530D+03, 0.12698D+03, 0.12798D+03,
+ 4 0.12898D+03, 0.12999D+03, 0.13101D+03, 0.13167D+03, 0.13233D+03,
+ 5 0.13300D+03, 0.13452D+03, 0.13606D+03, 0.13761D+03, 0.13919D+03,
+ 6 0.14078D+03, 0.14239D+03, 0.14401D+03, 0.14467D+03, 0.14533D+03,
+ 7 0.14600D+03/
+* Rotational irradiation geometry
+ DATA RT74G /
+ 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
+ 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18655D+02, 0.11608D+02,
+ 3 0.69154D+01, 0.37516D+01, 0.20643D+01, 0.11725D+01, 0.61014D+00,
+ 4 0.35636D+00, 0.27520D+00, 0.26661D+00, 0.30429D+00, 0.34763D+00,
+ 5 0.39750D+00, 0.50246D+00, 0.63525D+00, 0.80327D+00, 0.10159D+01,
+ 6 0.12851D+01, 0.16259D+01, 0.20570D+01, 0.25456D+01, 0.31506D+01,
+ 7 0.39000D+01, 0.46762D+01, 0.56061D+01, 0.67125D+01, 0.79141D+01,
+ 8 0.93241D+01, 0.10963D+02, 0.12696D+02, 0.14881D+02, 0.17439D+02,
+ 9 0.20700D+02, 0.24597D+02, 0.29229D+02, 0.34738D+02, 0.41889D+02,
+ * 0.50210D+02, 0.59027D+02, 0.68568D+02, 0.78333D+02, 0.89610D+02,
+ 1 0.10265D+03, 0.11038D+03, 0.11871D+03, 0.12770D+03, 0.13443D+03,
+ 2 0.14150D+03, 0.14896D+03, 0.15681D+03, 0.16272D+03, 0.16887D+03,
+ 3 0.17525D+03, 0.18162D+03, 0.18826D+03, 0.19518D+03, 0.19927D+03,
+ 4 0.20343D+03, 0.20768D+03, 0.21203D+03, 0.21473D+03, 0.21747D+03,
+ 5 0.22025D+03, 0.22378D+03, 0.22737D+03, 0.23102D+03, 0.23474D+03,
+ 6 0.23852D+03, 0.24237D+03, 0.24628D+03, 0.24898D+03, 0.25173D+03,
+ 7 0.25450D+03/
+* WORST possible geometry for the irradiation
+ DATA WT74G /
+ 1 0.54559D+03, 0.37003D+03, 0.24945D+03, 0.16496D+03, 0.10888D+03,
+ 2 0.71724D+02, 0.46941D+02, 0.29765D+02, 0.18658D+02, 0.11621D+02,
+ 3 0.71617D+01, 0.43328D+01, 0.26230D+01, 0.15912D+01, 0.97875D+00,
+ 4 0.62511D+00, 0.45300D+00, 0.41541D+00, 0.45697D+00, 0.50270D+00,
+ 5 0.55300D+00, 0.68970D+00, 0.86020D+00, 0.10728D+01, 0.13381D+01,
+ 6 0.16688D+01, 0.20814D+01, 0.25952D+01, 0.31566D+01, 0.38395D+01,
+ 7 0.46700D+01, 0.55604D+01, 0.65903D+01, 0.76977D+01, 0.89735D+01,
+ 8 0.10487D+02, 0.12358D+02, 0.13724D+02, 0.16218D+02, 0.18813D+02,
+ 9 0.23000D+02, 0.27007D+02, 0.31713D+02, 0.37238D+02, 0.45249D+02,
+ * 0.53784D+02, 0.59748D+02, 0.71072D+02, 0.83565D+02, 0.98879D+02,
+ 1 0.11700D+03, 0.12735D+03, 0.13861D+03, 0.15087D+03, 0.16030D+03,
+ 2 0.17027D+03, 0.18087D+03, 0.19208D+03, 0.20036D+03, 0.20899D+03,
+ 3 0.21800D+03, 0.22815D+03, 0.23877D+03, 0.24988D+03, 0.25669D+03,
+ 4 0.26363D+03, 0.27076D+03, 0.27805D+03, 0.28262D+03, 0.28897D+03,
+ 5 0.30000D+03, 0.30985D+03, 0.32002D+03, 0.33053D+03, 0.34138D+03,
+ 6 0.35259D+03, 0.36417D+03, 0.37616D+03, 0.39209D+03, 0.40869D+03,
+ 7 0.42600D+03/
+*
+* Effective dose for electrons from ICRP74 and Pelliccioni data
+* calculated with ICRP radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION E74E(NBIN1E),AE74E(NBIN1E)
+ DATA E74E /
+ 1 0.50000D-02, 0.64047D-02, 0.82039D-02, 0.10509D-01, 0.13461D-01,
+ 2 0.17242D-01, 0.22086D-01, 0.28291D-01, 0.36239D-01, 0.46420D-01,
+ 3 0.59460D-01, 0.76165D-01, 0.97562D-01, 0.12497D+00, 0.16008D+00,
+ 4 0.20505D+00, 0.26265D+00, 0.33644D+00, 0.43096D+00, 0.55203D+00,
+ 5 0.70711D+00, 0.90576D+00, 0.11602D+01, 0.14861D+01, 0.19037D+01,
+ 6 0.24385D+01, 0.31235D+01, 0.40010D+01, 0.51250D+01, 0.65647D+01,
+ 7 0.84090D+01, 0.10771D+02, 0.13797D+02, 0.17673D+02, 0.22638D+02,
+ 8 0.28998D+02, 0.37145D+02, 0.47580D+02, 0.60946D+02, 0.78068D+02,
+ 9 0.10000D+03/
+ DIMENSION AP74E(NBIN1E),RT74E(NBIN1E),WT74E(NBIN1E),
+ & AAP74E(NBIN1E),ART74E(NBIN1E),AWT74E(NBIN1E)
+* Anterior-Posterior irradiation
+ DATA AP74E /
+ 1 0.11700D+03, 0.14636D+03, 0.18310D+03, 0.22254D+03, 0.24108D+03,
+ 2 0.26117D+03, 0.28218D+03, 0.30367D+03, 0.32392D+03, 0.33732D+03,
+ 3 0.34585D+03, 0.35438D+03, 0.36312D+03, 0.36782D+03, 0.37210D+03,
+ 4 0.37679D+03, 0.38468D+03, 0.39275D+03, 0.40098D+03, 0.40825D+03,
+ 5 0.41392D+03, 0.41968D+03, 0.42412D+03, 0.42769D+03, 0.43128D+03,
+ 6 0.43436D+03, 0.43732D+03, 0.44030D+03, 0.44373D+03, 0.45113D+03,
+ 7 0.45866D+03, 0.46611D+03, 0.47323D+03, 0.48045D+03, 0.48778D+03,
+ 8 0.49522D+03, 0.50278D+03, 0.51045D+03, 0.51057D+03, 0.50878D+03,
+ 9 0.50700D+03/
+* Rotational irradiation geometry
+ DATA RT74E /
+ 1 0.46000D+02, 0.59898D+02, 0.78750D+02, 0.10262D+03, 0.12495D+03,
+ 2 0.15440D+03, 0.19064D+03, 0.23250D+03, 0.26904D+03, 0.29888D+03,
+ 3 0.31637D+03, 0.33039D+03, 0.34523D+03, 0.35974D+03, 0.37500D+03,
+ 4 0.39072D+03, 0.40301D+03, 0.41573D+03, 0.42889D+03, 0.44070D+03,
+ 5 0.45012D+03, 0.45980D+03, 0.47394D+03, 0.49154D+03, 0.50997D+03,
+ 6 0.52079D+03, 0.52974D+03, 0.53887D+03, 0.54941D+03, 0.57165D+03,
+ 7 0.59498D+03, 0.61689D+03, 0.63350D+03, 0.65061D+03, 0.66823D+03,
+ 8 0.68638D+03, 0.70507D+03, 0.72432D+03, 0.72548D+03, 0.72208D+03,
+ 9 0.71875D+03/
+* WORST possible geometry for the irradiation
+ DATA WT74E /
+ 1 0.11700D+03, 0.14636D+03, 0.18310D+03, 0.22254D+03, 0.24108D+03,
+ 2 0.26117D+03, 0.28218D+03, 0.30367D+03, 0.32392D+03, 0.34459D+03,
+ 3 0.35514D+03, 0.36261D+03, 0.37023D+03, 0.37763D+03, 0.38513D+03,
+ 4 0.39659D+03, 0.41277D+03, 0.42961D+03, 0.44714D+03, 0.46712D+03,
+ 5 0.48808D+03, 0.50997D+03, 0.52813D+03, 0.54373D+03, 0.57298D+03,
+ 6 0.58967D+03, 0.60557D+03, 0.63059D+03, 0.65754D+03, 0.69411D+03,
+ 7 0.73271D+03, 0.77101D+03, 0.80533D+03, 0.84117D+03, 0.87861D+03,
+ 8 0.91771D+03, 0.95856D+03, 0.10012D+04, 0.10295D+04, 0.10545D+04,
+ 9 0.10800D+04/
+*
+*-----------------------------------------------------------------------
+*
+* Effective dose for neutrons from ICRP74 and Pelliccioni data
+* calculated with the Pelliccioni radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION APPLN(NBIN1N),RTPLN(NBIN1N),WTPLN(NBIN1N),
+ & AAPPLN(NBIN1N),ARTPLN(NBIN1N),AWTPLN(NBIN1N)
+* Anterior-Posterior irradiation
+ DATA APPLN /
+ 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
+ 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
+ 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
+ 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
+ 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
+ 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
+ 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
+ 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
+ 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
+ * 486.0D0, 480.0D0, 458.0D0, 437.0D0, 429.0D0,
+ 1 409.0D0, 385.0D0, 378.0D0, 379.0D0, 399.0D0,
+ 2 422.0D0, 455.0D0, 503.0D0, 544.0D0, 610.0D0,
+ 3 707.0D0, 776.0D0, 852.0D0, 937.0D0, 996.0D0,
+ 4 1081.0D0, 1197.0D0, 1284.0D0, 1390.0D0, 1536.0D0,
+ 5 1658.0D0, 1856.0D0, 2150.0D0, 2368.0D0, 2620.0D0,
+ 6 2926.0D0, 3156.0D0, 3495.0D0, 3948.0D0, 4260.0D0,
+ 7 4600.0D0/
+* Rotational irradiation geometry
+ DATA RTPLN /
+ 1 2.99D0, 3.72D0, 4.40D0, 5.75D0, 6.43D0,
+ 2 7.27D0, 7.84D0, 8.31D0, 8.72D0, 8.90D0,
+ 3 8.92D0, 8.82D0, 8.69D0, 8.56D0, 8.40D0,
+ 4 8.34D0, 8.39D0, 9.06D0, 10.6D0, 13.8D0,
+ 5 16.9D0, 22.7D0, 27.8D0, 34.8D0, 45.4D0,
+ 6 54.8D0, 71.6D0, 99.4D0, 123.0D0, 144.0D0,
+ 7 154.0D0, 173.0D0, 234.0D0, 283.0D0, 315.0D0,
+ 8 335.0D0, 348.0D0, 358.0D0, 366.0D0, 373.0D0,
+ 9 378.0D0, 385.0D0, 390.0D0, 391.0D0, 393.0D0,
+ * 394.0D0, 395.0D0, 395.0D0, 404.0D0, 422.0D0,
+ 1 400.0D0, 417.0D0, 428.0D0, 441.0D0, 462.0D0,
+ 2 482.0D0, 513.0D0, 566.0D0, 615.0D0, 702.0D0,
+ 3 840.0D0, 947.0D0, 1070.0D0, 1216.0D0, 1324.0D0,
+ 4 1484.0D0, 1706.0D0, 1871.0D0, 2070.0D0, 2336.0D0,
+ 5 2555.0D0, 2910.0D0, 3447.0D0, 3862.0D0, 4360.0D0,
+ 6 5005.0D0, 5520.0D0, 6338.0D0, 7542.0D0, 8458.0D0,
+ 7 9550.0D0/
+* WORST possible geometry for the irradiation
+ DATA WTPLN /
+ 1 5.2D0, 6.6D0, 7.6D0, 10.0D0, 11.2D0,
+ 2 12.8D0, 13.8D0, 14.5D0, 15.0D0, 15.1D0,
+ 3 15.1D0, 14.8D0, 14.6D0, 14.4D0, 14.2D0,
+ 4 14.2D0, 14.4D0, 15.7D0, 18.3D0, 23.8D0,
+ 5 29.0D0, 38.5D0, 47.2D0, 59.8D0, 80.2D0,
+ 6 99.0D0, 133.0D0, 188.0D0, 231.0D0, 267.0D0,
+ 7 282.0D0, 310.0D0, 383.0D0, 432.0D0, 458.0D0,
+ 8 474.0D0, 483.0D0, 490.0D0, 494.0D0, 497.0D0,
+ 9 499.0D0, 499.0D0, 496.0D0, 494.0D0, 491.0D0,
+ * 486.0D0, 480.0D0, 458.0D0, 444.0D0, 459.0D0,
+ 1 460.0D0, 460.0D0, 461.0D0, 465.0D0, 481.0D0,
+ 2 501.0D0, 535.0D0, 595.0D0, 653.0D0, 757.0D0,
+ 3 927.0D0, 1062.0D0, 1220.0D0, 1413.0D0, 1557.0D0,
+ 4 1775.0D0, 2080.0D0, 2307.0D0, 2580.0D0, 2944.0D0,
+ 5 3242.0D0, 3726.0D0, 4460.0D0, 5030.0D0, 5720.0D0,
+ 6 6623.0D0, 7352.0D0, 8520.0D0, 10267.0D0, 11614.0D0,
+ 7 13240.0D0/
+*
+* Effective dose for protons from Pelliccioni data
+* calculated with the Pelliccioni radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EPLP(NBIN1P),AEPLP(NBIN1P)
+ DATA EPLP /
+ 1 5.00D-03, 6.50D-03, 8.00D-03, 1.00D-02, 1.25D-02,
+ 2 1.50D-02, 2.00D-02, 2.50D-02, 3.00D-02, 4.00D-02,
+ 3 5.00D-02, 6.50D-02, 8.00D-02, 1.00D-01, 1.25D-01,
+ 4 1.50D-01, 2.00D-01, 2.50D-01, 3.00D-01, 4.00D-01,
+ 5 5.00D-01, 6.50D-01, 8.00D-01, 1.00D+00, 1.25D+00,
+ 6 1.50D+00, 2.00D+00, 2.50D+00, 3.00D+00, 4.00D+00,
+ 7 5.00D+00, 6.50D+00, 8.00D+00, 1.00D+01, 1.25D+01,
+ 8 1.50D+01, 2.00D+01, 2.50D+01, 3.00D+01, 4.00D+01,
+ 9 5.00D+01, 6.50D+01, 8.00D+01, 1.00D+02, 1.25D+02,
+ * 1.50D+02, 2.00D+02, 2.50D+02, 3.00D+02, 4.00D+02,
+ 1 5.00D+02, 6.50D+02, 8.00D+02, 1.00D+03, 1.25D+03,
+ 2 1.50D+03, 2.00D+03, 2.50D+03, 3.00D+03, 4.00D+03,
+ 3 5.00D+03, 6.50D+03, 8.00D+03, 1.00D+04/
+ DIMENSION APPLP(NBIN1P),RTPLP(NBIN1P),WTPLP(NBIN1P),
+ & AAPPLP(NBIN1P),ARTPLP(NBIN1P),AWTPLP(NBIN1P)
+* Anterior-Posterior irradiation
+ DATA APPLP /
+ 1 23.0D0, 26.1D0, 29.7D0, 35.8D0, 46.6D0,
+ 2 62.1D0, 115.6D0, 218.6D0, 392.8D0, 1016.0D0,
+ 3 1924.0D0, 2787.8D0, 2934.5D0, 2728.0D0, 2432.2D0,
+ 4 2186.8D0, 1828.0D0, 1591.1D0, 1429.6D0, 1235.6D0,
+ 5 1136.0D0, 1069.5D0, 1039.0D0, 1012.0D0, 982.4D0,
+ 6 964.0D0, 968.0D0, 1013.7D0, 1075.3D0, 1198.7D0,
+ 7 1296.0D0, 1385.9D0, 1440.4D0, 1492.0D0, 1545.0D0,
+ 8 1590.1D0, 1662.3D0, 1717.0D0, 1759.4D0, 1818.9D0,
+ 9 1856.0D0, 1888.9D0, 1913.6D0, 1948.0D0, 1998.1D0,
+ * 2051.5D0, 2157.8D0, 2258.3D0, 2351.8D0, 2519.3D0,
+ 1 2664.9D0, 2852.0D0, 3010.2D0, 3188.0D0, 3371.4D0,
+ 2 3525.1D0, 3775.2D0, 3976.1D0, 4145.1D0, 4421.4D0,
+ 3 4644.7D0, 4918.8D0, 5145.5D0, 5400.0D0/
+* Rotational irradiation geometry
+ DATA RTPLP /
+ 1 19.3D0, 21.8D0, 24.4D0, 28.6D0, 35.0D0,
+ 2 43.4D0, 67.9D0, 106.7D0, 163.9D0, 350.3D0,
+ 3 632.8D0, 999.7D0, 1234.8D0, 1457.0D0, 1715.6D0,
+ 4 1921.9D0, 2064.0D0, 1908.6D0, 1680.3D0, 1321.7D0,
+ 5 1133.0D0, 1037.5D0, 1018.5D0, 1016.0D0, 1007.4D0,
+ 6 1000.5D0, 1011.0D0, 1052.2D0, 1106.9D0, 1223.9D0,
+ 7 1332.0D0, 1465.8D0, 1571.5D0, 1680.0D0, 1778.6D0,
+ 8 1851.6D0, 1955.6D0, 2029.8D0, 2088.9D0, 2184.5D0,
+ 9 2266.0D0, 2376.5D0, 2477.0D0, 2599.0D0, 2736.9D0,
+ * 2861.9D0, 3083.3D0, 3276.3D0, 3448.6D0, 3748.3D0,
+ 1 4005.3D0, 4335.9D0, 4619.7D0, 4947.0D0, 5297.5D0,
+ 2 5602.0D0, 6118.2D0, 6550.8D0, 6926.8D0, 7564.1D0,
+ 3 8098.4D0, 8774.8D0, 9349.9D0, 10010.0D0/
+* WORST possible geometry for the irradiation
+ DATA WTPLP /
+ 1 23.7D0, 26.7D0, 30.2D0, 36.2D0, 46.9D0,
+ 2 62.3D0, 115.6D0, 218.2D0, 392.1D0, 1016.0D0,
+ 3 1924.0D0, 2744.2D0, 2868.9D0, 2728.0D0, 2600.5D0,
+ 4 2502.4D0, 2256.0D0, 1955.7D0, 1695.1D0, 1351.9D0,
+ 5 1180.0D0, 1087.3D0, 1064.6D0, 1060.0D0, 1055.9D0,
+ 6 1054.6D0, 1068.0D0, 1100.6D0, 1143.6D0, 1243.8D0,
+ 7 1352.0D0, 1515.1D0, 1662.0D0, 1820.0D0, 1958.9D0,
+ 8 2055.8D0, 2183.7D0, 2269.1D0, 2335.9D0, 2448.9D0,
+ 9 2556.0D0, 2717.9D0, 2872.4D0, 3060.0D0, 3265.5D0,
+ * 3446.7D0, 3759.1D0, 4025.7D0, 4260.6D0, 4665.1D0,
+ 1 5009.9D0, 5453.3D0, 5835.6D0, 6280.0D0, 6761.8D0,
+ 2 7185.2D0, 7912.5D0, 8530.5D0, 9073.3D0, 10004.4D0,
+ 3 10794.5D0, 11806.0D0, 12674.6D0, 13680.0D0/
+*
+* Effective dose for charged pions from ICRP74 and Pelliccioni data
+* calculated with the Pelliccioni radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION APPLI(NBIN1I),RTPLI(NBIN1I),WTPLI(NBIN1I),
+ & APPLJ(NBIN1I),RTPLJ(NBIN1I),WTPLJ(NBIN1I),
+ & AAPPLI(NBIN1I),ARTPLI(NBIN1I),AWTPLI(NBIN1I),
+ & AAPPLJ(NBIN1I),ARTPLJ(NBIN1I),AWTPLJ(NBIN1I)
+* Anterior-Posterior irradiation
+* Positive pions
+ DATA APPLI /
+ 1 71.3D0, 77.4D0, 82.9D0, 92.9D0, 102.3D0,
+ 2 111.4D0, 129.4D0, 147.6D0, 175.8D0, 205.5D0,
+ 3 248.0D0, 305.6D0, 366.3D0, 491.1D0, 613.9D0,
+ 4 729.7D0, 930.0D0, 1083.3D0, 1233.9D0, 1320.6D0,
+ 5 1378.7D0, 1401.7D0, 1397.4D0, 1357.9D0, 1308.7D0,
+ 6 1262.3D0, 1186.8D0, 1133.1D0, 1082.0D0, 1054.7D0,
+ 7 1042.0D0, 1048.9D0, 1066.7D0, 1110.3D0, 1151.5D0,
+ 8 1185.4D0, 1230.6D0, 1250.0D0, 1249.3D0, 1238.7D0,
+ 9 1228.0D0, 1226.1D0, 1231.9D0, 1253.2D0, 1279.4D0,
+ * 1306.6D0, 1359.6D0, 1408.7D0, 1474.7D0, 1532.7D0,
+ 1 1600.0D0, 1671.9D0, 1734.0D0, 1838.2D0, 1924.4D0,
+ 2 1998.2D0, 2121.1D0, 2221.8D0, 2346.0D0, 2448.6D0,
+ 3 2563.0D0, 2681.5D0, 2781.4D0, 2945.1D0, 3077.4D0,
+ 4 3189.2D0, 3372.7D0, 3521.6D0, 3704.2D0, 3855.1D0,
+ 5 4024.0D0/
+* Negative pions
+ DATA APPLJ /
+ 1 471.6D0, 506.5D0, 537.4D0, 592.1D0, 641.1D0,
+ 2 686.9D0, 772.5D0, 853.8D0, 972.0D0, 1088.3D0,
+ 3 1243.1D0, 1435.8D0, 1618.2D0, 1929.9D0, 2159.7D0,
+ 4 2310.6D0, 2419.2D0, 2350.0D0, 2095.4D0, 1827.2D0,
+ 5 1537.2D0, 1284.1D0, 1118.2D0, 939.1D0, 867.2D0,
+ 6 841.0D0, 842.7D0, 871.4D0, 926.9D0, 981.2D0,
+ 7 1042.0D0, 1097.6D0, 1137.2D0, 1188.2D0, 1217.5D0,
+ 8 1234.5D0, 1249.0D0, 1250.0D0, 1242.0D0, 1233.6D0,
+ 9 1228.0D0, 1230.4D0, 1238.5D0, 1262.0D0, 1288.7D0,
+ * 1315.6D0, 1367.4D0, 1414.9D0, 1478.6D0, 1534.6D0,
+ 1 1600.0D0, 1670.4D0, 1731.5D0, 1834.7D0, 1920.6D0,
+ 2 1994.5D0, 2117.8D0, 2219.1D0, 2344.3D0, 2447.7D0,
+ 3 2563.0D0, 2682.2D0, 2782.6D0, 2946.8D0, 3079.3D0,
+ 4 3191.2D0, 3374.6D0, 3523.1D0, 3705.3D0, 3855.8D0,
+ 5 4024.0D0/
+* Rotational irradiation geometry
+* Positive pions
+ DATA RTPLI /
+ 1 37.0D0, 39.6D0, 42.0D0, 46.3D0, 50.4D0,
+ 2 54.5D0, 62.9D0, 71.7D0, 85.7D0, 101.0D0,
+ 3 123.5D0, 155.1D0, 189.7D0, 264.5D0, 342.9D0,
+ 4 421.7D0, 571.2D0, 701.7D0, 855.7D0, 969.0D0,
+ 5 1074.4D0, 1156.9D0, 1204.2D0, 1239.5D0, 1235.2D0,
+ 6 1215.9D0, 1167.8D0, 1124.8D0, 1078.1D0, 1050.2D0,
+ 7 1034.3D0, 1036.8D0, 1051.2D0, 1093.1D0, 1138.6D0,
+ 8 1182.3D0, 1258.8D0, 1320.0D0, 1388.6D0, 1439.8D0,
+ 9 1492.0D0, 1541.8D0, 1581.2D0, 1642.6D0, 1690.9D0,
+ * 1731.9D0, 1801.2D0, 1861.1D0, 1940.9D0, 2013.7D0,
+ 1 2104.0D0, 2209.2D0, 2306.8D0, 2482.1D0, 2636.1D0,
+ 2 2773.6D0, 3011.0D0, 3211.5D0, 3463.6D0, 3673.6D0,
+ 3 3907.0D0, 4145.7D0, 4344.4D0, 4664.8D0, 4919.7D0,
+ 4 5132.6D0, 5478.2D0, 5755.7D0, 6094.3D0, 6373.3D0,
+ 5 6686.0D0/
+* Negative pions
+ DATA RTPLJ /
+ 1 263.2D0, 276.2D0, 287.6D0, 308.3D0, 327.7D0,
+ 2 346.4D0, 383.5D0, 420.8D0, 478.2D0, 537.7D0,
+ 3 621.1D0, 730.7D0, 840.4D0, 1044.2D0, 1215.0D0,
+ 4 1347.7D0, 1504.9D0, 1545.6D0, 1479.2D0, 1366.5D0,
+ 5 1222.2D0, 1081.8D0, 983.3D0, 873.5D0, 832.4D0,
+ 6 822.1D0, 838.3D0, 871.8D0, 927.6D0, 978.9D0,
+ 7 1034.3D0, 1083.5D0, 1118.8D0, 1167.6D0, 1202.1D0,
+ 8 1229.9D0, 1277.1D0, 1320.0D0, 1380.7D0, 1434.0D0,
+ 9 1492.0D0, 1547.0D0, 1589.5D0, 1653.8D0, 1702.9D0,
+ * 1743.6D0, 1811.3D0, 1869.0D0, 1945.9D0, 2016.2D0,
+ 1 2104.0D0, 2207.2D0, 2303.5D0, 2477.5D0, 2631.1D0,
+ 2 2768.5D0, 3006.4D0, 3207.8D0, 3461.2D0, 3672.3D0,
+ 3 3907.0D0, 4146.8D0, 4346.2D0, 4667.4D0, 4922.7D0,
+ 4 5135.7D0, 5481.2D0, 5758.2D0, 6096.0D0, 6374.3D0,
+ 5 6686.0D0/
+* WORST possible geometry for the irradiation
+* Positive pions
+ DATA WTPLI /
+ 1 71.3D0, 77.4D0, 82.9D0, 92.9D0, 102.2D0,
+ 2 111.3D0, 129.2D0, 147.4D0, 175.6D0, 205.4D0,
+ 3 248.0D0, 305.8D0, 366.7D0, 492.1D0, 615.2D0,
+ 4 731.2D0, 931.3D0, 1083.3D0, 1231.0D0, 1315.0D0,
+ 5 1370.7D0, 1393.0D0, 1390.3D0, 1357.9D0, 1318.2D0,
+ 6 1280.7D0, 1218.8D0, 1173.2D0, 1127.4D0, 1100.0D0,
+ 7 1082.0D0, 1078.0D0, 1084.8D0, 1113.2D0, 1149.9D0,
+ 8 1189.6D0, 1270.7D0, 1350.0D0, 1461.3D0, 1557.6D0,
+ 9 1661.0D0, 1755.8D0, 1825.8D0, 1924.6D0, 1993.8D0,
+ * 2047.3D0, 2130.6D0, 2198.4D0, 2287.7D0, 2370.2D0,
+ 1 2476.0D0, 2604.6D0, 2727.7D0, 2955.3D0, 3160.1D0,
+ 2 3345.9D0, 3671.9D0, 3950.9D0, 4305.1D0, 4601.8D0,
+ 3 4932.0D0, 5269.1D0, 5548.6D0, 5997.6D0, 6353.0D0,
+ 4 6648.7D0, 7126.6D0, 7508.9D0, 7974.2D0, 8357.4D0,
+ 5 8787.0D0/
+* Negative pions
+ DATA WTPLJ /
+ 1 471.6D0, 506.4D0, 537.2D0, 591.6D0, 640.5D0,
+ 2 686.1D0, 771.6D0, 853.0D0, 971.2D0, 1087.8D0,
+ 3 1243.1D0, 1436.7D0, 1620.1D0, 1933.5D0, 2164.4D0,
+ 4 2315.4D0, 2422.4D0, 2350.0D0, 2090.4D0, 1819.3D0,
+ 5 1528.2D0, 1276.2D0, 1112.5D0, 939.1D0, 873.5D0,
+ 6 853.3D0, 865.4D0, 902.3D0, 965.7D0, 1023.3D0,
+ 7 1082.0D0, 1128.0D0, 1156.5D0, 1191.3D0, 1215.9D0,
+ 8 1238.9D0, 1289.7D0, 1350.0D0, 1452.7D0, 1551.2D0,
+ 9 1661.0D0, 1761.9D0, 1835.7D0, 1938.1D0, 2008.2D0,
+ * 2061.5D0, 2142.8D0, 2208.1D0, 2293.7D0, 2373.2D0,
+ 1 2476.0D0, 2602.2D0, 2723.8D0, 2949.7D0, 3153.9D0,
+ 2 3339.6D0, 3666.1D0, 3946.2D0, 4301.9D0, 4600.1D0,
+ 3 4932.0D0, 5270.5D0, 5551.0D0, 6001.1D0, 6357.0D0,
+ 4 6652.8D0, 7130.5D0, 7512.2D0, 7976.6D0, 8358.7D0,
+ 5 8787.0D0/
+*
+* Effective dose for muons from ICRP74 and Pelliccioni data
+* calculated with the ICRP/Pelliccioni radiation weighting factors Wr
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EBINM(NBIN1M),AEBINM(NBIN1M)
+ DATA EBINM/
+ 1 1.00D-03, 1.25D-03, 1.50D-03, 2.00D-03, 2.50D-03,
+ 2 3.00D-03, 4.00D-03, 5.00D-03, 6.50D-03, 8.00D-03,
+ 3 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
+ 4 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
+ 5 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
+ 6 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
+ 7 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
+ 8 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
+ 9 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
+ * 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
+ 1 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
+ 2 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
+ 3 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
+ 4 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
+ 5 1.00D+04/
+ DIMENSION APPLM(NBIN1M),RTPLM(NBIN1M),WTPLM(NBIN1M),
+ & AAPPLM(NBIN1M),ARTPLM(NBIN1M),AWTPLM(NBIN1M)
+* Anterior-Posterior irradiation
+ DATA APPLM /
+ 1 177.0D0, 178.0D0, 180.0D0, 182.0D0, 184.0D0,
+ 2 186.0D0, 188.0D0, 193.0D0, 200.0D0, 211.0D0,
+ 3 243.0D0, 337.7D0, 449.7D0, 672.5D0, 830.5D0,
+ 4 913.8D0, 921.6D0, 821.5D0, 635.1D0, 495.7D0,
+ 5 385.0D0, 353.0D0, 340.0D0, 339.0D0, 337.0D0,
+ 6 335.0D0, 334.0D0, 333.0D0, 332.0D0, 332.0D0,
+ 7 332.5D0, 333.0D0, 334.0D0, 335.0D0, 336.0D0,
+ 8 337.0D0, 338.0D0, 339.0D0, 340.0D0, 341.0D0,
+ 9 340.5D0, 337.0D0, 335.0D0, 333.3D0, 332.9D0,
+ * 333.1D0, 334.2D0, 335.5D0, 337.4D0, 338.9D0,
+ 1 340.5D0, 341.9D0, 343.0D0, 344.4D0, 345.3D0,
+ 2 346.0D0, 346.8D0, 347.4D0, 348.1D0, 348.5D0,
+ 3 349.0D0, 349.5D0, 349.9D0, 350.6D0, 351.2D0,
+ 4 351.7D0, 352.7D0, 353.7D0, 355.0D0, 356.2D0,
+ 5 357.7D0/
+* Rotational irradiation geometry
+ DATA RTPLM /
+ 1 94.0D0, 96.0D0, 98.0D0, 99.0D0, 100.0D0,
+ 2 101.0D0, 102.0D0, 105.0D0, 110.0D0, 116.0D0,
+ 3 124.2D0, 155.8D0, 192.5D0, 275.5D0, 363.2D0,
+ 4 445.4D0, 566.2D0, 614.1D0, 587.3D0, 529.0D0,
+ 5 463.8D0, 414.5D0, 385.4D0, 354.2D0, 339.2D0,
+ 6 331.4D0, 325.3D0, 324.5D0, 326.5D0, 329.4D0,
+ 7 332.8D0, 335.9D0, 338.1D0, 340.8D0, 342.4D0,
+ 8 343.5D0, 344.7D0, 345.4D0, 346.0D0, 346.3D0,
+ 9 346.6D0, 346.9D0, 347.1D0, 347.4D0, 347.7D0,
+ * 347.9D0, 348.2D0, 348.5D0, 348.8D0, 349.0D0,
+ 1 349.3D0, 349.6D0, 349.8D0, 350.1D0, 350.4D0,
+ 2 350.6D0, 351.0D0, 351.2D0, 351.6D0, 351.8D0,
+ 3 352.1D0, 352.4D0, 352.6D0, 353.0D0, 353.2D0,
+ 4 353.5D0, 353.8D0, 354.1D0, 354.4D0, 354.7D0,
+ 5 354.9D0/
+* WORST possible geometry for the irradiation
+ DATA WTPLM /
+ 1 177.0D0, 180.0D0, 182.0D0, 184.0D0, 186.0D0,
+ 2 188.0D0, 191.0D0, 195.0D0, 204.0D0, 215.0D0,
+ 3 243.0D0, 339.2D0, 452.5D0, 672.5D0, 818.9D0,
+ 4 892.6D0, 908.7D0, 847.5D0, 728.2D0, 628.6D0,
+ 5 537.5D0, 470.8D0, 430.9D0, 386.9D0, 364.4D0,
+ 6 351.6D0, 339.2D0, 334.4D0, 332.6D0, 333.1D0,
+ 7 334.5D0, 336.2D0, 337.6D0, 339.9D0, 341.6D0,
+ 8 343.0D0, 345.2D0, 346.8D0, 348.6D0, 349.8D0,
+ 9 351.0D0, 352.0D0, 352.7D0, 353.5D0, 354.0D0,
+ * 354.3D0, 354.7D0, 354.8D0, 355.0D0, 355.0D0,
+ 1 355.0D0, 355.0D0, 355.0D0, 354.9D0, 354.9D0,
+ 2 354.9D0, 354.8D0, 354.7D0, 354.6D0, 354.6D0,
+ 3 354.5D0, 354.4D0, 354.4D0, 354.3D0, 354.2D0,
+ 4 354.2D0, 354.1D0, 354.1D0, 354.0D0, 354.0D0,
+ 5 354.0D0/
+*
+*-----------------------------------------------------------------------
+*
+* Ambient dose equivalent for neutrons from ICRP74 and Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION AMBN(NBIN1N),AAMBN(NBIN1N)
+ DATA AMBN /
+ 1 6.6D0, 9.0D0, 10.6D0, 12.9D0, 13.5D0,
+ 2 13.6D0, 13.3D0, 12.9D0, 12.0D0, 11.3D0,
+ 3 10.6D0, 9.9D0, 9.4D0, 8.9D0, 8.3D0,
+ 4 7.9D0, 7.7D0, 8.0D0, 10.5D0, 16.6D0,
+ 5 23.7D0, 41.1D0, 60.0D0, 88.0D0, 132.0D0,
+ 6 170.0D0, 233.0D0, 322.0D0, 375.0D0, 400.0D0,
+ 7 416.0D0, 425.0D0, 420.0D0, 412.0D0, 408.0D0,
+ 8 405.0D0, 400.0D0, 405.0D0, 409.0D0, 420.0D0,
+ 9 440.0D0, 480.0D0, 520.0D0, 540.0D0, 555.0D0,
+ * 570.0D0, 600.0D0, 515.0D0, 400.0D0, 330.0D0,
+* high-energy extrapolation according to Pelliccioni
+ 1 285.0D0, 253.0D0, 247.0D0, 259.0D0, 298.0D0,
+ 2 335.0D0, 377.0D0, 421.0D0, 447.0D0, 476.0D0,
+ 3 499.0D0, 508.0D0, 511.0D0, 510.0D0, 508.0D0,
+ 4 504.0D0, 501.0D0, 504.0D0, 511.0D0, 528.0D0,
+ 5 545.0D0, 577.0D0, 629.0D0, 669.0D0, 717.0D0,
+ 6 777.0D0, 824.0D0, 896.0D0, 998.0D0, 1073.0D0,
+ 7 1160.0D0/
+* high-energy extrapolation according to Sannikov and Savitskaya
+C 1 285.0D0, 253.0D0, 285.0D0, 306.0D0, 420.0D0,
+C 2 500.0D0, 647.0D0, 733.0D0, 789.0D0, 862.0D0,
+C 3 951.0D0, 1000.0D0, 1050.0D0, 1050.0D0, 1046.0D0,
+C 4 1038.0D0, 1031.0D0, 1038.0D0, 1052.0D0, 1087.0D0,
+C 5 1122.0D0, 1188.0D0, 1295.0D0, 1377.0D0, 1476.0D0,
+C 6 1600.0D0, 1696.0D0, 1845.0D0, 2055.0D0, 2209.0D0,
+C 7 2388.0D0/
+* 1cm-Pb-modified Snoopy (A&B) response
+C DATA AMBN /
+CC & 0.13272E+01,0.35280E+01,0.44158E+01,0.61860E+01,0.76848E+01,
+CC & 0.91501E+01,0.10320E+02,0.11686E+02,0.13774E+02,0.15550E+02,
+CC & 0.17568E+02,0.20671E+02,0.23400E+02,0.26452E+02,0.31107E+02,
+CC & 0.35200E+02,0.40029E+02,0.47957E+02,0.55050E+02,0.63372E+02,
+CC & 0.68650E+02,0.76109E+02,0.83508E+02,0.98300E+02,0.12183E+03,
+CC & 0.14270E+03,0.18229E+03,0.25960E+03,0.33847E+03,0.39849E+03,
+CC & 0.41700E+03,0.44177E+03,0.46149E+03,0.45284E+03,0.44749E+03,
+CC & 0.44458E+03,0.43180E+03,0.40909E+03,0.38360E+03,0.36038E+03,
+CC & 0.34000E+03,0.33184E+03,0.32509E+03,0.32212E+03,0.31936E+03,
+CC & 0.31438E+03,0.31000E+03,0.34000E+03,0.32675E+03,0.32092E+03,
+CC & 0.33000E+03,0.33581E+03,0.34000E+03,0.35000E+03,0.42979E+03,
+CC & 0.47896E+03,0.53000E+03,0.10000E-12,0.10000E-12,0.10000E-12,
+CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
+CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
+CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
+CC & 0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,0.10000E-12,
+CC & 0.10000E-12/
+C & 0.13272E+01,0.35280E+01,0.44158E+01,0.61860E+01,0.76848E+01,
+C & 0.91501E+01,0.10320E+02,0.11686E+02,0.13774E+02,0.15550E+02,
+C & 0.17568E+02,0.20671E+02,0.23400E+02,0.26452E+02,0.31107E+02,
+C & 0.35200E+02,0.40029E+02,0.47957E+02,0.55050E+02,0.63372E+02,
+C & 0.68650E+02,0.76109E+02,0.83508E+02,0.98300E+02,0.12183E+03,
+C & 0.14270E+03,0.18229E+03,0.25960E+03,0.33847E+03,0.39849E+03,
+C & 0.41700E+03,0.44177E+03,0.46149E+03,0.45284E+03,0.44749E+03,
+C & 0.44458E+03,0.43180E+03,0.40909E+03,0.38360E+03,0.36038E+03,
+C & 0.34000E+03,0.33184E+03,0.32509E+03,0.32212E+03,0.31936E+03,
+C & 0.31438E+03,0.31000E+03,0.34000E+03,0.32675E+03,0.32092E+03,
+C & 0.33000E+03,0.33581E+03,0.34000E+03,0.35000E+03,0.42979E+03,
+C & 0.47896E+03,0.53000E+03,0.58919E+03,0.63516E+03,0.70610E+03,
+C & 0.80686E+03,0.88096E+03,0.96695E+03,0.10749E+04,0.11588E+04,
+C & 0.12882E+04,0.14721E+04,0.16072E+04,0.17641E+04,0.19612E+04,
+C & 0.21142E+04,0.23503E+04,0.26857E+04,0.29323E+04,0.32185E+04,
+C & 0.35780E+04,0.38571E+04,0.42879E+04,0.48998E+04,0.53498E+04,
+C & 0.58720E+04/
+*
+* Ambient dose equivalent for protons from Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EAMBP(NBIN2P),AEAMBP(NBIN2P)
+ DATA EAMBP /
+ 1 1.00D-02, 1.25D-02, 1.50D-02, 2.00D-02, 2.50D-02,
+ 2 3.00D-02, 4.00D-02, 5.00D-02, 6.50D-02, 8.00D-02,
+ 3 1.00D-01, 1.25D-01, 1.50D-01, 2.00D-01, 2.50D-01,
+ 4 3.00D-01, 4.00D-01, 5.00D-01, 6.50D-01, 8.00D-01,
+ 5 1.00D+00, 1.25D+00, 1.50D+00, 2.00D+00, 2.50D+00,
+ 6 3.00D+00, 4.00D+00, 5.00D+00, 6.50D+00, 8.00D+00,
+ 7 1.00D+01, 1.25D+01, 1.50D+01, 2.00D+01, 2.50D+01,
+ 8 3.00D+01, 4.00D+01, 5.00D+01, 6.50D+01, 8.00D+01,
+ 9 1.00D+02, 1.25D+02, 1.50D+02, 2.00D+02, 2.50D+02,
+ * 3.00D+02, 4.00D+02, 5.00D+02, 6.50D+02, 8.00D+02,
+ 1 1.00D+03, 1.25D+03, 1.50D+03, 2.00D+03, 2.50D+03,
+ 2 3.00D+03, 4.00D+03, 5.00D+03, 6.50D+03, 8.00D+03,
+ 3 1.00D+04/
+ DIMENSION AMBP(NBIN2P),AAMBP(NBIN2P)
+ DATA AMBP /
+ 1 .0D0, .0D0, .0D0, .0D0, .0D0,
+ 2 .0D0, 4289.5D0, 2790.0D0, 2000.8D0, 1710.4D0,
+ 3 1520.0D0, 1362.6D0, 1247.6D0, 1091.3D0, 991.0D0,
+ 4 921.9D0, 835.3D0, 786.0D0, 746.5D0, 726.7D0,
+ 5 714.4D0, 709.6D0, 710.3D0, 717.9D0, 727.7D0,
+ 6 737.2D0, 753.3D0, 765.0D0, 776.0D0, 782.7D0,
+ 7 788.1D0, 792.0D0, 794.4D0, 797.0D0, 798.6D0,
+ 8 799.8D0, 802.2D0, 804.9D0, 809.5D0, 814.6D0,
+ 9 822.0D0, 831.7D0, 841.2D0, 858.9D0, 874.6D0,
+ * 888.6D0, 912.4D0, 932.0D0, 955.9D0, 975.2D0,
+ 1 996.0D0, 1016.4D0, 1032.5D0, 1056.4D0, 1073.1D0,
+ 2 1085.2D0, 1100.8D0, 1109.3D0, 1114.7D0, 1115.1D0,
+ 3 1111.1D0/
+*
+* Ambient dose equivalent for ch. pions from ICRP74 and Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION AMBI(NBIN1I),AMBJ(NBIN1I),
+ & AAMBI(NBIN1I),AAMBJ(NBIN1I)
+* Positive pions
+ DATA AMBI /
+ 1 224.0D0, 245.3D0, 264.1D0, 296.5D0, 323.8D0,
+ 2 347.5D0, 387.2D0, 419.9D0, 459.9D0, 492.5D0,
+ 3 528.0D0, 563.6D0, 592.9D0, 640.3D0, 678.9D0,
+ 4 712.1D0, 769.2D0, 819.0D0, 884.5D0, 938.6D0,
+ 5 994.0D0, 1040.7D0, 1067.6D0, 1080.0D0, 1059.4D0,
+ 6 1026.9D0, 958.6D0, 900.6D0, 835.9D0, 792.4D0,
+ 7 757.0D0, 735.3D0, 726.0D0, 722.9D0, 727.1D0,
+ 8 733.1D0, 743.7D0, 750.0D0, 752.4D0, 751.2D0,
+ 9 748.0D0, 743.9D0, 740.2D0, 734.0D0, 729.2D0,
+ * 725.4D0, 720.1D0, 716.7D0, 714.0D0, 712.9D0,
+ 1 713.0D0, 714.7D0, 717.1D0, 722.7D0, 728.4D0,
+ 2 733.8D0, 743.6D0, 752.1D0, 762.9D0, 772.0D0,
+ 3 782.0D0, 792.2D0, 800.7D0, 814.4D0, 825.2D0,
+ 4 834.1D0, 848.4D0, 859.8D0, 873.3D0, 884.2D0,
+ 5 896.0D0/
+* Negative pions
+ DATA AMBJ /
+ 1 1130.0D0, 1245.2D0, 1345.3D0, 1507.9D0, 1631.0D0,
+ 2 1723.7D0, 1843.3D0, 1902.4D0, 1921.2D0, 1890.7D0,
+ 3 1810.0D0, 1684.7D0, 1563.7D0, 1368.4D0, 1231.1D0,
+ 4 1135.8D0, 1024.1D0, 975.0D0, 960.9D0, 974.3D0,
+ 5 1003.4D0, 1038.1D0, 1063.2D0, 1080.0D0, 1062.8D0,
+ 6 1031.8D0, 964.1D0, 905.4D0, 839.2D0, 794.2D0,
+ 7 757.0D0, 733.4D0, 722.6D0, 717.3D0, 720.2D0,
+ 8 725.5D0, 735.8D0, 743.0D0, 747.7D0, 748.8D0,
+ 9 748.0D0, 745.7D0, 742.9D0, 737.5D0, 732.8D0,
+ * 728.9D0, 722.9D0, 718.9D0, 715.2D0, 713.5D0,
+ 1 713.0D0, 714.2D0, 716.4D0, 721.7D0, 727.4D0,
+ 2 732.9D0, 742.8D0, 751.5D0, 762.6D0, 771.8D0,
+ 3 782.0D0, 792.4D0, 801.0D0, 814.7D0, 825.5D0,
+ 4 834.5D0, 848.8D0, 860.0D0, 873.5D0, 884.3D0,
+ 5 896.0D0/
+*
+* Ambient dose equivalent for muons from ICRP74 and Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION AMBM(NBIN1M),AAMBM(NBIN1M)
+ DATA AMBM /
+ 1 394.0D0, 411.7D0, 426.6D0, 450.7D0, 469.5D0,
+ 2 484.8D0, 508.3D0, 525.6D0, 544.3D0, 557.4D0,
+ 3 569.5D0, 578.9D0, 583.9D0, 585.1D0, 579.3D0,
+ 4 569.5D0, 543.9D0, 515.5D0, 474.1D0, 439.0D0,
+ 5 402.6D0, 369.8D0, 347.0D0, 320.0D0, 307.2D0,
+ 6 301.2D0, 297.8D0, 298.9D0, 302.8D0, 307.0D0,
+ 7 311.5D0, 315.1D0, 317.3D0, 319.7D0, 320.7D0,
+ 8 321.1D0, 321.2D0, 320.9D0, 320.5D0, 320.1D0,
+ 9 320.0D0, 320.2D0, 320.5D0, 321.3D0, 322.0D0,
+ * 322.7D0, 323.7D0, 324.4D0, 324.9D0, 325.1D0,
+ 1 325.0D0, 324.5D0, 323.9D0, 322.5D0, 321.3D0,
+ 2 320.2D0, 318.6D0, 317.4D0, 316.2D0, 315.4D0,
+ 3 315.0D0, 314.9D0, 315.1D0, 315.9D0, 316.8D0,
+ 4 317.7D0, 319.5D0, 321.0D0, 323.0D0, 324.7D0,
+ 5 326.5D0/
+*
+* Ambient dose equivalent for photons from ICRP74 and Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION AMBG(NBIN1G),AAMBG(NBIN1G)
+ DATA AMBG /
+ 1 0.29512E-08, 0.16406E-07, 0.91199E-07, 0.50699E-06, 0.28184E-05,
+ 2 0.15667E-04, 0.87096E-04, 0.48417E-03, 0.26915E-02, 0.14962E-01,
+ 3 0.83176E-01, 0.46238E+00, 0.93903E+00, 0.10505E+01, 0.93070E+00,
+ 4 0.76366E+00, 0.62173E+00, 0.52108E+00, 0.51322E+00, 0.55911E+00,
+ 5 0.61958E+00, 0.72946E+00, 0.92815E+00, 0.12264E+01, 0.15403E+01,
+ 6 0.19055E+01, 0.23480E+01, 0.27945E+01, 0.36249E+01, 0.42697E+01,
+ 7 0.51844E+01, 0.61688E+01, 0.71466E+01, 0.82509E+01, 0.95675E+01,
+ 8 0.10544E+02, 0.10715E+02, 0.10375E+02, 0.94276E+01, 0.91033E+01,
+ 9 0.87619E+01, 0.84957E+01, 0.83464E+01, 0.82813E+01, 0.82585E+01,
+ * 0.82149E+01, 0.82566E+01, 0.86437E+01, 0.88675E+01, 0.89207E+01,
+ 1 0.90012E+01, 0.92640E+01, 0.96872E+01, 0.10186E+02, 0.10715E+02,
+ 2 0.11169E+02, 0.11561E+02, 0.11803E+02, 0.11885E+02, 0.11830E+02,
+ 3 0.11695E+02, 0.11535E+02, 0.11429E+02, 0.11508E+02, 0.11830E+02,
+ 4 0.12331E+02, 0.12882E+02, 0.13305E+02, 0.13428E+02, 0.13092E+02,
+ 5 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02,
+ 6 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02, 0.12190E+02,
+ 7 0.12190E+02/
+*
+* Ambient dose equivalent for electrons from ICRP74 and Pelliccioni data
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EAMBE(NBIN2E),AEAMBE(NBIN2E)
+ DATA EAMBE /
+ 1 0.20000E-02, 0.22896E-02, 0.26212E-02, 0.30008E-02, 0.34354E-02,
+ 2 0.39329E-02, 0.45025E-02, 0.51545E-02, 0.59010E-02, 0.67556E-02,
+ 3 0.77339E-02, 0.88540E-02, 0.10136E-01, 0.11604E-01, 0.13285E-01,
+ 4 0.15209E-01, 0.17411E-01, 0.19932E-01, 0.22819E-01, 0.26124E-01,
+ 5 0.29907E-01, 0.34238E-01, 0.39196E-01, 0.44873E-01, 0.51371E-01,
+ 6 0.58811E-01, 0.67328E-01, 0.77078E-01, 0.88241E-01, 0.10102E+00,
+ 7 0.11565E+00, 0.13240E+00, 0.15157E+00, 0.17352E+00, 0.19865E+00,
+ 8 0.22742E+00, 0.26036E+00, 0.29806E+00, 0.34122E+00, 0.39064E+00,
+ 9 0.44721E+00, 0.51198E+00, 0.58612E+00, 0.67101E+00, 0.76818E+00,
+ * 0.87943E+00, 0.10068E+01, 0.11526E+01, 0.13195E+01, 0.15106E+01,
+ 1 0.17294E+01, 0.19798E+01, 0.22665E+01, 0.25948E+01, 0.29705E+01,
+ 2 0.34007E+01, 0.38932E+01, 0.44570E+01, 0.51025E+01, 0.58414E+01,
+ 3 0.66874E+01, 0.76559E+01, 0.87646E+01, 0.10034E+02, 0.11487E+02,
+ 4 0.13151E+02, 0.15055E+02, 0.17235E+02, 0.19731E+02, 0.22589E+02,
+ 5 0.25860E+02, 0.29605E+02, 0.33892E+02, 0.38801E+02, 0.44420E+02,
+ 6 0.50853E+02, 0.58217E+02, 0.66648E+02, 0.76300E+02, 0.87350E+02,
+ 7 0.10000E+03/
+*
+ DIMENSION AMBE(NBIN2E),AAMBE(NBIN2E)
+ DATA AMBE /
+ 1 1.00000E-12, 0.16082E+03, 0.19885E+03, 0.32509E+03, 0.41665E+03,
+ 2 0.45092E+03, 0.44206E+03, 0.41582E+03, 0.38863E+03, 0.36471E+03,
+ 3 0.34608E+03, 0.33323E+03, 0.32465E+03, 0.32115E+03, 0.32097E+03,
+ 4 0.32279E+03, 0.32481E+03, 0.32583E+03, 0.32546E+03, 0.32233E+03,
+ 5 0.31724E+03, 0.31179E+03, 0.30738E+03, 0.30701E+03, 0.31035E+03,
+ 6 0.31401E+03, 0.31705E+03, 0.31769E+03, 0.31635E+03, 0.31477E+03,
+ 7 0.31497E+03, 0.31646E+03, 0.31903E+03, 0.32163E+03, 0.32354E+03,
+ 8 0.32444E+03, 0.32485E+03, 0.32398E+03, 0.32261E+03, 0.32087E+03,
+ 9 0.31913E+03, 0.31669E+03, 0.31472E+03, 0.31256E+03, 0.31087E+03,
+ * 0.30919E+03, 0.30832E+03, 0.30848E+03, 0.30932E+03, 0.31016E+03,
+ 1 0.31154E+03, 0.31390E+03, 0.31573E+03, 0.31831E+03, 0.32056E+03,
+ 2 0.32277E+03, 0.32480E+03, 0.32730E+03, 0.32833E+03, 0.32986E+03,
+ 3 0.33037E+03, 0.33037E+03, 0.32972E+03, 0.32810E+03, 0.32810E+03,
+ 4 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
+ 5 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
+ 6 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03, 0.32810E+03,
+ 7 0.32810E+03/
+*
+*-----------------------------------------------------------------------
+*
+* Ambient dose equivalent for neutrons ("GRS"-conversion factors)
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION H10N(NBIN2N),H10PR(NBIN2N),H10PI(NBIN2N),H10MU(NBIN2M),
+ & AH10MU(NBIN2M)
+ DATA H10N/
+ 1 7.62D0, 7.84D0, 8.06D0, 8.28D0, 8.50D0,
+ 2 8.71D0, 8.92D0, 9.12D0, 9.32D0, 9.50D0,
+ 3 9.67D0, 9.83D0, 9.97D0, 10.10D0, 10.21D0,
+ 4 10.30D0, 10.38D0, 10.43D0, 10.46D0, 10.47D0,
+ 5 10.47D0, 10.43D0, 10.38D0, 10.31D0, 10.22D0,
+ 6 10.11D0, 9.99D0, 9.85D0, 9.69D0, 9.52D0,
+ 7 9.35D0, 9.16D0, 8.96D0, 8.76D0, 8.56D0,
+ 8 8.35D0, 8.15D0, 7.94D0, 7.74D0, 7.54D0,
+ 9 7.35D0, 7.16D0, 6.99D0, 6.82D0, 6.67D0,
+ * 6.53D0, 6.40D0, 6.30D0, 6.22D0, 6.16D0,
+ 1 6.13D0, 6.14D0, 6.19D0, 6.29D0, 6.45D0,
+ 2 6.69D0, 7.03D0, 7.49D0, 8.11D0, 8.93D0,
+ 3 10.04D0, 11.51D0, 13.47D0, 16.08D0, 19.58D0,
+ 4 24.24D0, 30.42D0, 38.53D0, 49.03D0, 62.35D0,
+ 5 78.83D0, 98.66D0, 121.76D0, 147.73D0, 175.91D0,
+ 6 205.40D0, 235.18D0, 264.24D0, 291.70D0, 316.84D0,
+ 7 339.20D0, 358.53D0, 374.78D0, 388.05D0, 398.56D0,
+ 8 406.61D0, 412.66D0, 417.53D0, 423.43D0, 436.87D0,
+ 9 472.90D0, 539.16D0, 598.43D0, 623.89D0, 595.0D0 ,
+ * 524.0D0 , 365.0D0 , 323.0D0 , 321.0D0 , 323.0D0 ,
+ 1 325.0D0 , 328.0D0 , 333.0D0 , 338.0D0 , 344.0D0 ,
+ 2 353.0D0 , 371.0D0 , 482.0D0 , 581.0D0 , 677.0D0 ,
+ 3 777.0D0 , 882.0D0 , 1005.0D0 , 1125.0D0 , 1265.0D0 ,
+ 4 1410.0D0 , 1570.0D0 , 1740.0D0 , 1930.0D0 , 2130.0D0 ,
+ 5 2360.0D0 , 2600.0D0 , 2880.0D0 , 3180.0D0 , 3500.0D0 ,
+ 6 3860.0D0 , 4250.0D0 , 4700.0D0 , 5180.0D0 , 5610.0D0 ,
+ 7 6290.0D0 , 6900.0D0 , 7570.0D0 , 8280.0D0 , 9000.0D0 ,
+ 8 9830.0D0 , 10700.0D0 , 11600.0D0 , 12650.0D0 , 13600.0D0 ,
+ 9 14750.0D0 , 15900.0D0 , 17100.0D0 , 18450.0D0 , 19700.0D0 ,
+ * 21100.0D0 , 22500.0D0 , 23900.0D0 , 25500.0D0 , 27000.0D0 /
+*
+* Ambient dose equivalent for protons ("GRS"-conversion factors)
+* Energy in GeV, dose in pSv.cm^2
+*
+ DATA H10PR/
+ 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ * 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ * 3130.0D0, 2530.0D0, 2130.0D0, 1825.0D0, 1585.0D0,
+ 1 1390.0D0, 1230.0D0, 1105.0D0, 1000.0D0, 923.0D0,
+ 2 864.0D0, 828.0D0, 898.0D0, 965.0D0, 1040.0D0,
+ 3 1125.0D0, 1220.0D0, 1335.0D0, 1455.0D0, 1595.0D0,
+ 4 1740.0D0, 1905.0D0, 2080.0D0, 2270.0D0, 2480.0D0,
+ 5 2710.0D0, 2960.0D0, 3240.0D0, 3550.0D0, 3870.0D0,
+ 6 4240.0D0, 4630.0D0, 5090.0D0, 5570.0D0, 6000.0D0,
+ 7 6690.0D0, 7300.0D0, 7980.0D0, 8690.0D0, 9410.0D0,
+ 8 10250.0D0, 11100.0D0, 12050.0D0, 13100.0D0, 14050.0D0,
+ 9 15200.0D0, 16350.0D0, 17550.0D0, 18900.0D0, 20100.0D0,
+ * 21600.0D0, 22300.0D0, 24400.0D0, 26000.0D0, 27500.0D0/
+*
+* Ambient dose equivalent for charged pions ("GRS"-conversion factors)
+* Energy in GeV, dose in pSv.cm^2
+*
+ DATA H10PI/
+ 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ * 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 1 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 2 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 3 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 4 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 5 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 6 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 7 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 8 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ 9 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0,
+ * 1640.0D0, 1560.0D0, 1495.0D0, 1445.0D0, 1405.0D0,
+ 1 1375.0D0, 1355.0D0, 1345.0D0, 1335.0D0, 1330.0D0,
+ 2 1330.0D0, 1330.0D0, 1330.0D0, 1335.0D0, 1340.0D0,
+ 3 1345.0D0, 1350.0D0, 1360.0D0, 1485.0D0, 1630.0D0,
+ 4 1780.0D0, 1945.0D0, 2120.0D0, 2310.0D0, 2520.0D0,
+ 5 2750.0D0, 3000.0D0, 3280.0D0, 3580.0D0, 3910.0D0,
+ 6 4270.0D0, 4670.0D0, 5120.0D0, 5600.0D0, 6040.0D0,
+ 7 6720.0D0, 7330.0D0, 8010.0D0, 8720.0D0, 9440.0D0,
+ 8 10300.0D0, 11150.0D0, 12050.0D0, 13100.0D0, 14050.0D0,
+ 9 15200.0D0, 16350.0D0, 17550.0D0, 18900.0D0, 20200.0D0,
+ * 21600.0D0, 23000.0D0, 24400.0D0, 26000.0D0, 27500.0D0/
+*
+* Ambient dose equivalent for muons ("GRS"-conversion factors)
+* Energy in GeV, dose in pSv.cm^2
+*
+ DIMENSION EAMBM(NBIN2M),AEAMBM(NBIN2M)
+ DATA EAMBM/
+ 1 1.0D-02, 2.0D-02, 5.0D-02, 1.0D-01, 2.0D-01,
+ 2 5.0D-01, 1.0D+00, 2.0D+00, 5.0D+00, 1.0D+01,
+ 3 2.0D+01, 5.0D+01, 1.0D+02, 2.0D+02, 5.0D+02,
+ 4 1.0D+03/
+ DATA H10MU /
+ 1 2600.0D0, 2600.0D0, 2600.0D0, 2600.0D0, 350.0D0,
+ 2 330.0D0, 350.0D0, 360.0D0, 380.0D0, 390.0D0,
+ 3 410.0D0, 430.0D0, 450.0D0, 480.0D0, 560.0D0,
+ 4 680.0D0/
+*
+*-----------------------------------------------------------------------
+*
+ DATA LFIRST /.TRUE.,.TRUE.,.TRUE.,.TRUE.,.TRUE.,
+ & .TRUE.,.TRUE.,.TRUE.,.TRUE./
+ DATA CSET /' '/
+
+* default is zero weighting
+ FLUSCW = 0.0D0
+ ITEST = 0
+ IJ = IIJ
+
+* calculate logarithms of energy-bins at the very first call
+ IF (LFIRST(0)) THEN
+ DO 1 I=1,NBIN1N
+ AEBINN(I) = LOG10(EBINN(I))
+ AAP74N(I) = LOG10(AP74N(I))
+ ART74N(I) = LOG10(RT74N(I))
+ AWT74N(I) = LOG10(WT74N(I))
+ AAPPLN(I) = LOG10(APPLN(I))
+ ARTPLN(I) = LOG10(RTPLN(I))
+ AWTPLN(I) = LOG10(WTPLN(I))
+ AAMBN(I) = LOG10(AMBN(I))
+ 1 CONTINUE
+ DO 2 I=1,NBIN1P
+ AE74P(I) = LOG10(E74P(I))
+ AEPLP(I) = LOG10(EPLP(I))
+ AAP74P(I) = LOG10(AP74P(I))
+ ART74P(I) = LOG10(RT74P(I))
+ AWT74P(I) = LOG10(WT74P(I))
+ AAPPLP(I) = LOG10(APPLP(I))
+ ARTPLP(I) = LOG10(RTPLP(I))
+ AWTPLP(I) = LOG10(WTPLP(I))
+ 2 CONTINUE
+ DO 3 I=1,NBIN1I
+ AEBINI(I) = LOG10(EBINI(I))
+ AAP74I(I) = LOG10(AP74I(I))
+ ART74I(I) = LOG10(RT74I(I))
+ AWT74I(I) = LOG10(WT74I(I))
+ AAP74J(I) = LOG10(AP74J(I))
+ ART74J(I) = LOG10(RT74J(I))
+ AWT74J(I) = LOG10(WT74J(I))
+ AAPPLI(I) = LOG10(APPLI(I))
+ ARTPLI(I) = LOG10(RTPLI(I))
+ AWTPLI(I) = LOG10(WTPLI(I))
+ AAPPLJ(I) = LOG10(APPLJ(I))
+ ARTPLJ(I) = LOG10(RTPLJ(I))
+ AWTPLJ(I) = LOG10(WTPLJ(I))
+ AAMBI(I) = LOG10(AMBI(I))
+ AAMBJ(I) = LOG10(AMBJ(I))
+ 3 CONTINUE
+ DO 4 I=1,NBIN1M
+ AEBINM(I) = LOG10(EBINM(I))
+ AAPPLM(I) = LOG10(APPLM(I))
+ ARTPLM(I) = LOG10(RTPLM(I))
+ AWTPLM(I) = LOG10(WTPLM(I))
+ AAMBM(I) = LOG10(AMBM(I))
+ 4 CONTINUE
+ DO 5 I=1,NBIN2P
+ AEAMBP(I) = LOG10(EAMBP(I))
+ IF (AMBP(I).GT.0.0D0) THEN
+ AAMBP(I) = LOG10(AMBP(I))
+ ELSE
+ AAMBP(I) = LOG10(ANGLGB)
+ ENDIF
+ 5 CONTINUE
+ DO 6 I=1,NBIN2M
+ AEAMBM(I) = LOG10(EAMBM(I))
+ AH10MU(I) = LOG10(H10MU(I))
+ 6 CONTINUE
+ DO 7 I=1,NBIN1G
+ AE74G(I) = LOG10(E74G(I))
+ AAP74G(I) = LOG10(AP74G(I))
+ ART74G(I) = LOG10(RT74G(I))
+ AWT74G(I) = LOG10(WT74G(I))
+ AAMBG(I) = LOG10(AMBG(I))
+ 7 CONTINUE
+ DO 8 I=1,NBIN1E
+ AE74E(I) = LOG10(E74E(I))
+ AAP74E(I) = LOG10(AP74E(I))
+ ART74E(I) = LOG10(RT74E(I))
+ AWT74E(I) = LOG10(WT74E(I))
+ 8 CONTINUE
+ DO 9 I=1,NBIN2E
+ AEAMBE(I) = LOG10(EAMBE(I))
+ AAMBE(I) = LOG10(AMBE(I))
+ 9 CONTINUE
+ LFIRST(0) = .FALSE.
+ ENDIF
+
+* test print
+ IF (IJ.LT.-10000) THEN
+* test print for set number itest (i.e. -20001 gives the proton
+* conversion factors for set 2)
+* return particle-id of conversion factor set used for particle IJ
+* via LLO
+ ITEST = -IJ/10000
+ IJ = -IJ-10000*ITEST
+ IF (ITT(IJ).EQ.1) THEN
+ LLO = 8
+ ELSEIF (ITT(IJ).EQ.2) THEN
+ LLO = 1
+ ELSEIF (ITT(IJ).EQ.3) THEN
+ LLO = 13
+ ELSEIF (ITT(IJ).EQ.4) THEN
+ LLO = 14
+ ELSEIF (ITT(IJ).EQ.5) THEN
+ LLO = 10
+ ELSEIF (ITT(IJ).EQ.6) THEN
+ LLO = 11
+ ELSEIF (ITT(IJ).EQ.7) THEN
+ LLO = 7
+ ELSEIF (ITT(IJ).EQ.8) THEN
+ LLO = 3
+ ELSEIF (ITT(IJ).EQ.99) THEN
+ LLO = 0
+ ELSE
+ LLO = -1
+ ENDIF
+ ELSE
+* ISCRNG = 1 for usrbdx
+* = 2 for usrbin
+* = 3 for usrtrack
+* return unit weight for all other detectors
+ IF ((ISCRNG.NE.1).AND.(ISCRNG.NE.2).AND.(ISCRNG.NE.3)) THEN
+ FLUSCW = 1.0D0
+ RETURN
+ ENDIF
+ ENDIF
+*
+* At this point the detector is either usrbdx or usrbin or usrtrack.
+*
+* Now check if any of the available conversion factor sets is requested
+* in which case the default conversion factor is zero. Otherwise the
+* default factor is one.
+ IF (ITEST.EQ.0) THEN
+* usrbdx
+ IF (ISCRNG.EQ.1) THEN
+ CSET = TITUSX(JSCRNG)
+* usrtrack
+ ELSEIF (ISCRNG.EQ.3) THEN
+ CSET = TITUTC(JSCRNG)
+* usrbin
+ ELSE
+ CSET = TITUSB(JSCRNG)
+ ENDIF
+ IF ((CSET(1:3).EQ.'EAP').OR.(CSET(1:3).EQ.'eap').OR.
+ & (CSET(1:3).EQ.'ERT').OR.(CSET(1:3).EQ.'ert').OR.
+ & (CSET(1:3).EQ.'EWT').OR.(CSET(1:3).EQ.'ewt').OR.
+ & (CSET(1:3).EQ.'AMB').OR.(CSET(1:3).EQ.'amb')) THEN
+ FLUSCW = 0.0D0
+ ELSE
+* skip the rest for all other sdum's
+ FLUSCW = 1.0D0
+ RETURN
+ ENDIF
+ ENDIF
+*
+* At this point one of the available conversion factor sets is
+* requested by sdum.
+*
+* check for particle type
+* return zero weight for all particles with ij < 1 or > 40 and for
+* all particles with flag=99 (see array ITT)
+ IF ((IJ.LT.1).OR.(IJ.GT.40)) RETURN
+ IPART = ITT(IJ)
+ IF (IPART.EQ.99) RETURN
+
+* kinetic energy
+ IF (PLA.LT.0.0D0) THEN
+ EKIN = ABS(PLA)
+ ELSEIF (PLA.GT.0.0D0) THEN
+ EKIN = SQRT(PLA**2+AM(IJ)**2)-AM(IJ)
+ ELSE
+ RETURN
+ ENDIF
+*
+*
+* conversion factor set
+ IF (ITEST.EQ.0) THEN
+* for muons there is only a Pelliccioni weighting factor set
+ IF (IPART.EQ.5) THEN
+ IF ((CSET(1:1).EQ.'E').AND.(CSET(4:5).EQ.'74'))
+ & CSET(4:5) = 'MP'
+ IF ((CSET(1:1).EQ.'e').AND.(CSET(4:5).EQ.'74'))
+ & CSET(4:5) = 'mp'
+ ENDIF
+* for photons and electrons there is only a ICRP weighting factor set
+ IF ((IPART.EQ.7).OR.(IPART.EQ.8)) THEN
+ IF (((CSET(1:1).EQ.'E' ).OR.(CSET(1:1).EQ.'e' )).AND.
+ & ((CSET(4:5).EQ.'MP').OR.(CSET(4:5).EQ.'mp'))) THEN
+ CSET(4:5) = '74'
+* for photons and electrons there is no GRS amb. dose equivalent data set
+ ELSEIF ((CSET(1:3).EQ.'AMB').OR.(CSET(1:3).EQ.'amb')) THEN
+ IF ((CSET(4:5).EQ.'GS').OR.(CSET(4:5).EQ.'gs')) RETURN
+ ENDIF
+* skip if only hadronic part is requested
+ IF (CSET(6:6).EQ.'1') RETURN
+ ELSE
+* skip if only electromagnetic part is requested
+ IF (CSET(6:6).EQ.'2') RETURN
+ ENDIF
+ ELSE
+ IF (IPART.EQ.5) THEN
+ IF (ITEST.EQ.1) ITEST = 4
+ IF (ITEST.EQ.2) ITEST = 5
+ IF (ITEST.EQ.3) ITEST = 6
+ ENDIF
+ IF ((IPART.EQ.7).OR.(IPART.EQ.8)) THEN
+ IF (ITEST.EQ.4) ITEST = 1
+ IF (ITEST.EQ.5) ITEST = 2
+ IF (ITEST.EQ.6) ITEST = 3
+ IF (ITEST.GE.9) RETURN
+ ENDIF
+ ENDIF
+*
+*
+* Effective dose (ICRP radiation weighting factors Wr)
+* Anterior-Posterior irradiation
+ IF ((CSET(1:5).EQ.'EAP74').OR.(CSET(1:5).EQ.'eap74')
+ & .OR.(ITEST.EQ.1)) THEN
+ IF (LFIRST(1).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1002)
+ WRITE(LUNOUT,1004)
+ LFIRST(1) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = AP74N(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = AP74N(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = AP74N(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAP74N(IBINHI)-AAP74N(IBINLO))+AAP74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 10 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 11
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 10
+ 11 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAP74N(IBINHI)-AAP74N(IBINLO))+AAP74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74P(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
+ FLUSCW = AP74P(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74P(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = AP74P(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AE74P(IBINLO))/
+ & (AE74P(IBINHI)-AE74P(IBINLO))*
+ & (AAP74P(IBINHI)-AAP74P(IBINLO))+AAP74P(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = AP74I(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = AP74I(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAP74I(IBINHI)-AAP74I(IBINLO))+AAP74I(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = AP74J(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = AP74J(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAP74J(IBINHI)-AAP74J(IBINLO))+AAP74J(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.7) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74G(1)) THEN
+ FLUSCW = AP74G(1)
+ ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
+ FLUSCW = AP74G(NBIN1G)
+ ELSE
+ XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
+ & *DBLE(NBIN1G-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74G(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1G) THEN
+ FLUSCW = AP74G(NBIN1G)
+ ELSE
+ FLUSCW =( AEKIN-AE74G(IBINLO))/
+ & (AE74G(IBINHI)-AE74G(IBINLO))*
+ & (AAP74G(IBINHI)-AAP74G(IBINLO))+AAP74G(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.8) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74E(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
+ FLUSCW = AP74E(NBIN1E)
+ ELSE
+ XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
+ & *DBLE(NBIN1E-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74E(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1E) THEN
+ FLUSCW = AP74E(NBIN1E)
+ ELSE
+ FLUSCW =( AEKIN-AE74E(IBINLO))/
+ & (AE74E(IBINHI)-AE74E(IBINLO))*
+ & (AAP74E(IBINHI)-AAP74E(IBINLO))+AAP74E(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Effective dose (ICRP radiation weighting factors Wr)
+* Rotational irradiation geometry
+ ELSEIF ((CSET(1:5).EQ.'ERT74').OR.(CSET(1:5).EQ.'ert74')
+ & .OR.(ITEST.EQ.2)) THEN
+ IF (LFIRST(2).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1002)
+ WRITE(LUNOUT,1005)
+ LFIRST(2) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = RT74N(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = RT74N(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = RT74N(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (ART74N(IBINHI)-ART74N(IBINLO))+ART74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 20 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 21
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 20
+ 21 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (ART74N(IBINHI)-ART74N(IBINLO))+ART74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74P(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
+ FLUSCW = RT74P(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74P(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = RT74P(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AE74P(IBINLO))/
+ & (AE74P(IBINHI)-AE74P(IBINLO))*
+ & (ART74P(IBINHI)-ART74P(IBINLO))+ART74P(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = RT74I(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = RT74I(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (ART74I(IBINHI)-ART74I(IBINLO))+ART74I(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = RT74J(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = RT74J(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (ART74J(IBINHI)-ART74J(IBINLO))+ART74J(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.7) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74G(1)) THEN
+ FLUSCW = RT74G(1)
+ ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
+ FLUSCW = RT74G(NBIN1G)
+ ELSE
+ XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
+ & *DBLE(NBIN1G-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74G(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1G) THEN
+ FLUSCW = RT74G(NBIN1G)
+ ELSE
+ FLUSCW =( AEKIN-AE74G(IBINLO))/
+ & (AE74G(IBINHI)-AE74G(IBINLO))*
+ & (ART74G(IBINHI)-ART74G(IBINLO))+ART74G(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.8) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74E(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
+ FLUSCW = RT74E(NBIN1E)
+ ELSE
+ XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
+ & *DBLE(NBIN1E-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74E(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1E) THEN
+ FLUSCW = RT74E(NBIN1E)
+ ELSE
+ FLUSCW =( AEKIN-AE74E(IBINLO))/
+ & (AE74E(IBINHI)-AE74E(IBINLO))*
+ & (ART74E(IBINHI)-ART74E(IBINLO))+ART74E(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Effective dose (ICRP radiation weighting factors Wr)
+* Worst possible geometry for the irradiation
+ ELSEIF ((CSET(1:5).EQ.'EWT74').OR.(CSET(1:5).EQ.'ewt74')
+ & .OR.(ITEST.EQ.3)) THEN
+ IF (LFIRST(3).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1002)
+ WRITE(LUNOUT,1006)
+ LFIRST(3) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = WT74N(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = WT74N(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = WT74N(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AWT74N(IBINHI)-AWT74N(IBINLO))+AWT74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 30 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 31
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 30
+ 31 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AWT74N(IBINHI)-AWT74N(IBINLO))+AWT74N(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74P(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74P(NBIN1P)) THEN
+ FLUSCW = WT74P(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AE74P(1))/(AE74P(NBIN1P)-AE74P(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74P(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74P(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = WT74P(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AE74P(IBINLO))/
+ & (AE74P(IBINHI)-AE74P(IBINLO))*
+ & (AWT74P(IBINHI)-AWT74P(IBINLO))+AWT74P(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = WT74I(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = WT74I(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AWT74I(IBINHI)-AWT74I(IBINLO))+AWT74I(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = WT74J(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = WT74J(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AWT74J(IBINHI)-AWT74J(IBINLO))+AWT74J(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.7) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74G(1)) THEN
+ FLUSCW = WT74G(1)
+ ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
+ FLUSCW = WT74G(NBIN1G)
+ ELSE
+ XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
+ & *DBLE(NBIN1G-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74G(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1G) THEN
+ FLUSCW = WT74G(NBIN1G)
+ ELSE
+ FLUSCW =( AEKIN-AE74G(IBINLO))/
+ & (AE74G(IBINHI)-AE74G(IBINLO))*
+ & (AWT74G(IBINHI)-AWT74G(IBINLO))+AWT74G(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.8) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74E(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AE74E(NBIN1E)) THEN
+ FLUSCW = WT74E(NBIN1E)
+ ELSE
+ XBIN = (AEKIN-AE74E(1))/(AE74E(NBIN1E)-AE74E(1))
+ & *DBLE(NBIN1E-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74E(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74E(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1E) THEN
+ FLUSCW = WT74E(NBIN1E)
+ ELSE
+ FLUSCW =( AEKIN-AE74E(IBINLO))/
+ & (AE74E(IBINHI)-AE74E(IBINLO))*
+ & (AWT74E(IBINHI)-AWT74E(IBINLO))+AWT74E(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+* Effective dose (Pelliccioni radiation weighting factors Wr)
+* Anterior-Posterior irradiation
+ ELSEIF ((CSET(1:5).EQ.'EAPMP').OR.(CSET(1:5).EQ.'eapmp')
+ & .OR.(ITEST.EQ.4)) THEN
+ IF (LFIRST(4).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1003)
+ WRITE(LUNOUT,1004)
+ LFIRST(4) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = APPLN(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = APPLN(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = APPLN(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAPPLN(IBINHI)-AAPPLN(IBINLO))+AAPPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 40 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 41
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 40
+ 41 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAPPLN(IBINHI)-AAPPLN(IBINLO))+AAPPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEPLP(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
+ FLUSCW = APPLP(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEPLP(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = APPLP(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AEPLP(IBINLO))/
+ & (AEPLP(IBINHI)-AEPLP(IBINLO))*
+ & (AAPPLP(IBINHI)-AAPPLP(IBINLO))+AAPPLP(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = APPLI(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = APPLI(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAPPLI(IBINHI)-AAPPLI(IBINLO))+AAPPLI(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = APPLJ(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = APPLJ(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAPPLJ(IBINHI)-AAPPLJ(IBINLO))+AAPPLJ(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.5) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINM(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
+ FLUSCW = APPLM(NBIN1M)
+ ELSE
+ XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
+ & *DBLE(NBIN1M-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINM(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1M) THEN
+ FLUSCW = APPLM(NBIN1M)
+ ELSE
+ FLUSCW =( AEKIN-AEBINM(IBINLO))/
+ & (AEBINM(IBINHI)-AEBINM(IBINLO))*
+ & (AAPPLM(IBINHI)-AAPPLM(IBINLO))+AAPPLM(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Effective dose (Pelliccioni radiation weighting factors Wr)
+* Rotational irradiation geometry
+ ELSEIF ((CSET(1:5).EQ.'ERTMP').OR.(CSET(1:5).EQ.'ertmp')
+ & .OR.(ITEST.EQ.5)) THEN
+ IF (LFIRST(5).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1003)
+ WRITE(LUNOUT,1005)
+ LFIRST(5) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = RTPLN(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = RTPLN(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = RTPLN(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (ARTPLN(IBINHI)-ARTPLN(IBINLO))+ARTPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 50 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 51
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 50
+ 51 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (ARTPLN(IBINHI)-ARTPLN(IBINLO))+ARTPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEPLP(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
+ FLUSCW = RTPLP(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEPLP(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = RTPLP(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AEPLP(IBINLO))/
+ & (AEPLP(IBINHI)-AEPLP(IBINLO))*
+ & (ARTPLP(IBINHI)-ARTPLP(IBINLO))+ARTPLP(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = RTPLI(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = RTPLI(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (ARTPLI(IBINHI)-ARTPLI(IBINLO))+ARTPLI(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = RTPLJ(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = RTPLJ(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (ARTPLJ(IBINHI)-ARTPLJ(IBINLO))+ARTPLJ(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.5) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINM(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
+ FLUSCW = RTPLM(NBIN1M)
+ ELSE
+ XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
+ & *DBLE(NBIN1M-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINM(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1M) THEN
+ FLUSCW = RTPLM(NBIN1M)
+ ELSE
+ FLUSCW =( AEKIN-AEBINM(IBINLO))/
+ & (AEBINM(IBINHI)-AEBINM(IBINLO))*
+ & (ARTPLM(IBINHI)-ARTPLM(IBINLO))+ARTPLM(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Effective dose (Pelliccioni radiation weighting factors Wr)
+* Worst possible geometry for the irradiation
+ ELSEIF ((CSET(1:5).EQ.'EWTMP').OR.(CSET(1:5).EQ.'ewtmp')
+ & .OR.(ITEST.EQ.6)) THEN
+ IF (LFIRST(6).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1000)
+ WRITE(LUNOUT,1003)
+ WRITE(LUNOUT,1006)
+ LFIRST(6) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = WTPLN(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = WTPLN(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = WTPLN(NBIN1N)
+ ELSE
+ FLUSCW =( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AWTPLN(IBINHI)-AWTPLN(IBINLO))+AWTPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 60 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 61
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 60
+ 61 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AWTPLN(IBINHI)-AWTPLN(IBINLO))+AWTPLN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEPLP(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEPLP(NBIN1P)) THEN
+ FLUSCW = WTPLP(NBIN1P)
+ ELSE
+ XBIN = (AEKIN-AEPLP(1))/(AEPLP(NBIN1P)-AEPLP(1))
+ & *DBLE(NBIN1P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEPLP(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEPLP(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1P) THEN
+ FLUSCW = WTPLP(NBIN1P)
+ ELSE
+ FLUSCW =( AEKIN-AEPLP(IBINLO))/
+ & (AEPLP(IBINHI)-AEPLP(IBINLO))*
+ & (AWTPLP(IBINHI)-AWTPLP(IBINLO))+AWTPLP(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = WTPLI(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = WTPLI(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AWTPLI(IBINHI)-AWTPLI(IBINLO))+AWTPLI(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = WTPLJ(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = WTPLJ(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AWTPLJ(IBINHI)-AWTPLJ(IBINLO))+AWTPLJ(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.5) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINM(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
+ FLUSCW = WTPLM(NBIN1M)
+ ELSE
+ XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
+ & *DBLE(NBIN1M-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINM(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1M) THEN
+ FLUSCW = WTPLM(NBIN1M)
+ ELSE
+ FLUSCW =( AEKIN-AEBINM(IBINLO))/
+ & (AEBINM(IBINHI)-AEBINM(IBINLO))*
+ & (AWTPLM(IBINHI)-AWTPLM(IBINLO))+AWTPLM(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Ambient dose equivalent (from ICRP74 and Pelliccioni data)
+ ELSEIF ((CSET(1:5).EQ.'AMB74').OR.(CSET(1:5).EQ.'amb74')
+ & .OR.(ITEST.EQ.7)) THEN
+ IF (LFIRST(7).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1001)
+ WRITE(LUNOUT,1008)
+ LFIRST(7) = .FALSE.
+ ENDIF
+ IF (IPART.EQ.1) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINN(1)) THEN
+ FLUSCW = AMBN(1)
+ ELSEIF (AEKIN.GT.AEBINN(NBIN1N)) THEN
+ FLUSCW = AMBN(NBIN1N)
+ ELSEIF ((AEKIN.GE.AEBINN(47)).AND.(AEKIN.LE.AEBINN(NBIN1N)))
+ & THEN
+ XBIN = (AEKIN-AEBINN(47))/(AEBINN(NBIN1N)-AEBINN(47))
+ & *DBLE(NBIN1N-47)+47.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINN(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINN(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1N) THEN
+ FLUSCW = AMBN(NBIN1N)
+ ELSE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAMBN(IBINHI)-AAMBN(IBINLO))+AAMBN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSE
+ IBINLO = 1
+ IBINHI = 47
+ 70 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 71
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEBINN(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 70
+ 71 CONTINUE
+ FLUSCW = ( AEKIN-AEBINN(IBINLO))/
+ & (AEBINN(IBINHI)-AEBINN(IBINLO))*
+ & (AAMBN(IBINHI)-AAMBN(IBINLO))+AAMBN(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ELSEIF (IPART.EQ.2) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEAMBP(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEAMBP(NBIN2P)) THEN
+ FLUSCW = AMBP(NBIN2P)
+ ELSE
+ XBIN = (AEKIN-AEAMBP(1))/(AEAMBP(NBIN2P)-AEAMBP(1))
+ & *DBLE(NBIN2P-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEAMBP(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEAMBP(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN2P) THEN
+ FLUSCW = AMBP(NBIN2P)
+ ELSE
+ FLUSCW = ( AEKIN-AEAMBP(IBINLO))/
+ & (AEAMBP(IBINHI)-AEAMBP(IBINLO))*
+ & (AAMBP(IBINHI)-AAMBP(IBINLO))+AAMBP(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.3) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = AMBI(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = AMBI(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAMBI(IBINHI)-AAMBI(IBINLO))+AAMBI(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.4) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINI(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINI(NBIN1I)) THEN
+ FLUSCW = AMBJ(NBIN1I)
+ ELSE
+ XBIN = (AEKIN-AEBINI(1))/(AEBINI(NBIN1I)-AEBINI(1))
+ & *DBLE(NBIN1I-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINI(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINI(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1I) THEN
+ FLUSCW = AMBJ(NBIN1I)
+ ELSE
+ FLUSCW =( AEKIN-AEBINI(IBINLO))/
+ & (AEBINI(IBINHI)-AEBINI(IBINLO))*
+ & (AAMBJ(IBINHI)-AAMBJ(IBINLO))+AAMBJ(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.5) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEBINM(1)) THEN
+ RETURN
+ ELSEIF (AEKIN.GT.AEBINM(NBIN1M)) THEN
+ FLUSCW = AMBM(NBIN1M)
+ ELSE
+ XBIN = (AEKIN-AEBINM(1))/(AEBINM(NBIN1M)-AEBINM(1))
+ & *DBLE(NBIN1M-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEBINM(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEBINM(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1M) THEN
+ FLUSCW = AMBM(NBIN1M)
+ ELSE
+ FLUSCW =( AEKIN-AEBINM(IBINLO))/
+ & (AEBINM(IBINHI)-AEBINM(IBINLO))*
+ & (AAMBM(IBINHI)-AAMBM(IBINLO))+AAMBM(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.7) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AE74G(1)) THEN
+ FLUSCW = AMBG(1)
+ ELSEIF (AEKIN.GT.AE74G(NBIN1G)) THEN
+ FLUSCW = AMBG(NBIN1G)
+ ELSE
+ XBIN = (AEKIN-AE74G(1))/(AE74G(NBIN1G)-AE74G(1))
+ & *DBLE(NBIN1G-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AE74G(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AE74G(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN1G) THEN
+ FLUSCW = AMBG(NBIN1G)
+ ELSE
+ FLUSCW =( AEKIN-AE74G(IBINLO))/
+ & (AE74G(IBINHI)-AE74G(IBINLO))*
+ & (AAMBG(IBINHI)-AAMBG(IBINLO))+AAMBG(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSEIF (IPART.EQ.8) THEN
+ AEKIN = LOG10(EKIN)
+ IF (AEKIN.LT.AEAMBE(1)) THEN
+ FLUSCW = AMBE(1)
+ ELSEIF (AEKIN.GT.AEAMBE(NBIN2E)) THEN
+ FLUSCW = AMBE(NBIN2E)
+ ELSE
+ XBIN = (AEKIN-AEAMBE(1))/(AEAMBE(NBIN2E)-AEAMBE(1))
+ & *DBLE(NBIN2E-1)+1.0D0
+ IBINLO = INT(XBIN)
+ IBINHI = IBINLO+1
+ IF (AEKIN.LT.AEAMBE(IBINLO)) THEN
+ IBINLO = IBINLO-1
+ IBINHI = IBINLO+1
+ ELSEIF (AEKIN.GT.AEAMBE(IBINHI)) THEN
+ IBINLO = IBINLO+1
+ IBINHI = IBINLO+1
+ ENDIF
+ IF (IBINHI.GT.NBIN2E) THEN
+ FLUSCW = AMBE(NBIN2E)
+ ELSE
+ FLUSCW =( AEKIN-AEAMBE(IBINLO))/
+ & (AEAMBE(IBINHI)-AEAMBE(IBINLO))*
+ & (AAMBE(IBINHI)-AAMBE(IBINLO))+AAMBE(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ENDIF
+*
+*
+* Ambient dose equivalent (old GRS conversion factors)
+ ELSEIF ((CSET(1:5).EQ.'AMBGS').OR.(CSET(1:5).EQ.'ambgs')
+ & .OR.(ITEST.EQ.8)) THEN
+ IF (LFIRST(8).AND.(ITEST.EQ.0)) THEN
+ WRITE(LUNOUT,1001)
+ WRITE(LUNOUT,1007)
+ LFIRST(8) = .FALSE.
+ ENDIF
+ AEKIN = LOG10(EKIN)
+ IF (IPART.EQ.1) THEN
+ XDUM = (11.0D0+AEKIN)*10.0D0
+ KENERG = INT(XDUM)+1
+ IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
+ FLUSCW = H10N(KENERG)
+ ELSEIF (IPART.EQ.2) THEN
+ XDUM = (11.0D0+AEKIN)*10.0D0
+ KENERG = INT(XDUM)+1
+ IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
+ FLUSCW = H10PR(KENERG)
+ ELSEIF ((IPART.EQ.3).OR.(IPART.EQ.4)) THEN
+ XDUM = (11.0D0+AEKIN)*10.0D0
+ KENERG = INT(XDUM)+1
+ IF ((KENERG.LT.1).OR.(KENERG.GT.NBIN2N)) RETURN
+ FLUSCW = H10PI(KENERG)
+ ELSEIF (IPART.EQ.5) THEN
+ IF (AEKIN.LT.AEAMBM(1)) THEN
+ FLUSCW = H10MU(1)
+ ELSEIF (AEKIN.GT.AEAMBM(NBIN2M)) THEN
+ FLUSCW = H10MU(NBIN2M)
+ ELSE
+ IBINLO = 1
+ IBINHI = NBIN2M
+ 80 CONTINUE
+ IF ((IBINHI-IBINLO).EQ.1) GOTO 81
+ KK = (IBINHI+IBINLO)/2
+ IF (AEKIN.LE.AEAMBM(KK)) THEN
+ IBINHI = KK
+ ELSE
+ IBINLO = KK
+ ENDIF
+ GOTO 80
+ 81 CONTINUE
+ FLUSCW = ( AEKIN-AEAMBM(IBINLO))/
+ & (AEAMBM(IBINHI)-AEAMBM(IBINLO))*
+ & (AH10MU(IBINHI)-AH10MU(IBINLO))+AH10MU(IBINLO)
+ FLUSCW = 10.0D0**FLUSCW
+ ENDIF
+ ENDIF
+ ELSE
+* this point should never be reached since the check for a valid
+* sdum has already been done above
+ FLUSCW = 0.0D0
+ ENDIF
+
+ LSCZER = .FALSE.
+
+ RETURN
+
+ 1000 FORMAT(1X,'FLUSCW: direct conversion of fluence to effective ',
+ & 'dose requested with')
+ 1001 FORMAT(1X,'FLUSCW: direct conversion of fluence to ambient ',
+ & 'dose eq. requested with')
+ 1002 FORMAT(12X,'- ICRP radiation weighting factors Wr')
+ 1003 FORMAT(12X,'- Pelliccioni radiation weighting factors Wr')
+ 1004 FORMAT(12X,'- anterior-posterior irradiation')
+ 1005 FORMAT(12X,'- rotational irradiation geometry')
+ 1006 FORMAT(12X,'- worst possible geometry for the irradiation')
+ 1007 FORMAT(12X,'- GRS conversion factors')
+ 1008 FORMAT(12X,'- ICRP74 and Pelliccioni conversion factors')
+
+ END
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff