************************************************************************ * 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) * * * * * * Version for FLUKA2000 by S. Roesler and G. R. Stevenson, 25.11.2000 * ************************************************************************ * * *=== 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 WRITE(LUNOUT,*) 'FLUSCW for deq conversion activated ' 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