Normalisation in USRBDX and USRYIELD plus discrepancies in

From: Tomasz Cybulski <T.Cybulski_at_liv.ac.uk>
Date: Mon, 04 Jul 2011 13:34:41 +0100

estimators
Content-Type: multipart/mixed;
  boundary="------------050008020502040502000001"
X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.2.7
(smtp1.mi.infn.it [192.84.138.69]); Mon, 04 Jul 2011 14:36:04 +0200 (CEST)
Sender: owner-fluka-discuss_at_mi.infn.it

This is a multi-part message in MIME format.
--------------050008020502040502000001
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit

Dear Fluka Users,

I have started using fluka some time ago. I started my simulations but
in the meanwhile I compared the results I get from usrbdx and usryield.
I have problem with the normalisation of the results.

I am shooting with a 60MeV proton beam on a copper block (input file
attached) - the first approximation of the experiment before further
simulations. In the same run of the Fluka I use both usrbdx and usryield
to score the double differential proton fluence from the surface of the
Cu block (surface irradiated with the beam) to vacuum (in the opposite
direction of the beam direction). I set up 20 runs - 1e6 protons in each
cycle. The usrbdx settings: fluence one way, 200 energy bins from 0.0 -
0.062 GeV, one angular bin, the surface area 16cm2 ; the usryield
settings: polar angle 90-180 degrees in 1 bin, energy 0.0 - 0.062GeV in
200 bins.

According to the manual, the the results of usrbdx are normalized to
particles/GeV/sr/cm2 per primary. The usryield results are normalized to
particles/ GeV /sr. When I multiply all the results of usrbdx by the bin
energy width, 2pi solid angle and 16cm2 the result differ from the
Usryield, multiplied by bin energy width and 2pi solid angle,
approximately by the factor 16. I wanted to see the total particle fluence.

Is the difference, visible in plots + output files, my fault or this is
the issue of the normalization in program?

The second matter is concerning the estimators itself. When I multiply
the usryield by the factor of 16 (plot 'ProtYield_energy16.gif') and
have 'the same' particle flux like in usrbdx (plot
'Proton_fluence_usr_bdx.gif') with there is quite important
inconsistency between results. Namely, some of the results agree within
the estimated uncertainty but most of them does not. Once again the
estimators were used simultaneously during the same run of the program.

Is it just the matter of the estimators' statistical nature or there is
any other - mine or software - reason?

Thank you for your response in advance.

With best regards,
Tomasz Cybulski

P.S. Please, find .inp, .tab.lis and .gif files attached.

-- 
------------------------------------
Tomasz Cybulski
QUASAR GROUP
Marie Curie Fellow
The Cockcroft Institute
4 Keckwick Lane
Daresbury
Warrington
WA4 4AD
Department of Physics
The University of Liverpool
Liverpool
L69 7ZE
United Kingdom
phone: +44 1925 864057
e-mail: t.cybulski_at_liv.ac.uk
------------------------------------
--------------050008020502040502000001
Content-Type: text/plain;
  name="Prot_60_usryield_91_tab.lis"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
  filename="Prot_60_usryield_91_tab.lis"
  # Detector n:  1  YilProt
  # N. of x1 intervals 200
    0.000000      3.1000000E-04  1.3348479E-03   13.78303
   3.1000000E-04  6.2000001E-04  1.1808270E-03   13.04348
   6.2000001E-04  9.2999998E-04  1.9766018E-03   10.86980
   9.2999998E-04  1.2400000E-03  1.8995912E-03   10.20719
   1.2400000E-03  1.5500000E-03  2.7980467E-03   9.734241
   1.5500000E-03  1.8600000E-03  2.7723764E-03   8.325208
   1.8600000E-03  2.1700000E-03  3.3371197E-03   7.276195
   2.1700000E-03  2.4800000E-03  2.8750570E-03   10.94673
   2.4800000E-03  2.7900001E-03  3.5424810E-03   9.801307
   2.7900001E-03  3.1000001E-03  3.7221720E-03   8.424259
   3.1000001E-03  3.4100001E-03  4.0045436E-03   6.997402
   3.4100001E-03  3.7199999E-03  4.0815542E-03   7.980237
   3.7199999E-03  4.0300000E-03  3.8761930E-03   8.347911
   4.0300000E-03  4.3400000E-03  4.3125856E-03   7.887274
   4.3400000E-03  4.6500000E-03  3.0804183E-03   8.289357
   4.6500000E-03  4.9600000E-03  4.0815542E-03   8.236984
   4.9600000E-03  5.2700001E-03  3.6708317E-03   9.369750
   5.2700001E-03  5.5800001E-03  3.3884600E-03   6.447012
   5.5800001E-03  5.8900001E-03  3.5424810E-03   7.352390
   5.8900001E-03  6.2000002E-03  3.2857794E-03   8.736480
   6.2000002E-03  6.5100002E-03  3.4141303E-03   10.01500
   6.5100002E-03  6.8200002E-03  2.8493868E-03   6.466832
   6.8200002E-03  7.1300003E-03  3.0034077E-03   9.864689
   7.1300003E-03  7.4399998E-03  2.9520674E-03   8.545616
   7.4399998E-03  7.7499999E-03  2.4129944E-03   9.144898
   7.7499999E-03  8.0599999E-03  2.8493870E-03   7.449117
   8.0599999E-03  8.3699999E-03  2.5670151E-03   9.403247
   8.3699999E-03  8.6800000E-03  2.2589734E-03   12.49184
   8.6800000E-03  8.9900000E-03  2.5670151E-03   9.514532
   8.9900000E-03  9.3000000E-03  2.3103138E-03   12.43567
   9.3000000E-03  9.6100001E-03  3.0034077E-03   8.966322
   9.6100001E-03  9.9200001E-03  2.5156748E-03   7.534892
   9.9200001E-03  1.0230000E-02  1.9509315E-03   9.846467
   1.0230000E-02  1.0540000E-02  2.1306225E-03   10.95553
   1.0540000E-02  1.0850000E-02  2.2333032E-03   8.384545
   1.0850000E-02  1.1160000E-02  1.6428897E-03   9.511074
   1.1160000E-02  1.1470000E-02  1.7712405E-03   13.70474
   1.1470000E-02  1.1780000E-02  2.0022718E-03   11.89653
   1.1780000E-02  1.2090000E-02  1.8225808E-03   11.09486
   1.2090000E-02  1.2400000E-02  1.8482510E-03   11.47079
   1.2400000E-02  1.2710000E-02  1.6172195E-03   9.001259
   1.2710000E-02  1.3020000E-02  1.5658792E-03   14.16722
   1.3020000E-02  1.3330000E-02  1.5145390E-03   9.356437
   1.3330000E-02  1.3640000E-02  1.4118584E-03   12.89028
   1.3640000E-02  1.3950001E-02  1.5145390E-03   13.57604
   1.3950001E-02  1.4260001E-02  1.9509316E-03   10.90055
   1.4260001E-02  1.4570001E-02  1.5915495E-03   11.67780
   1.4570001E-02  1.4880000E-02  1.6685599E-03   10.67629
   1.4880000E-02  1.5190000E-02  1.3861882E-03   12.63152
   1.5190000E-02  1.5500000E-02  1.4118584E-03   14.41931
   1.5500000E-02  1.5810000E-02  1.3861882E-03   17.22584
   1.5810000E-02  1.6120000E-02  1.3861882E-03   13.98763
   1.6120000E-02  1.6430000E-02  1.3605180E-03   10.34345
   1.6430000E-02  1.6740000E-02  1.4375285E-03   17.30345
   1.6740000E-02  1.7050000E-02  1.3605180E-03   12.62735
   1.7050000E-02  1.7360000E-02  1.3091777E-03   14.91255
   1.7360000E-02  1.7670000E-02  9.2412549E-04   17.38230
   1.7670000E-02  1.7980000E-02  1.1551569E-03   16.40140
   1.7980000E-02  1.8290000E-02  1.1038166E-03   17.62248
   1.8290000E-02  1.8600000E-02  9.7546575E-04   18.66674
   1.8600000E-02  1.8910000E-02  1.2835076E-03   16.02629
   1.8910000E-02  1.9220000E-02  1.0524762E-03   12.49941
   1.9220000E-02  1.9530000E-02  1.0524762E-03   16.01332
   1.9530000E-02  1.9840000E-02  8.9845533E-04   21.88732
   1.9840000E-02  2.0150000E-02  1.2064972E-03   14.23935
   2.0150000E-02  2.0460000E-02  9.4979565E-04   22.62093
   2.0460000E-02  2.0770000E-02  1.5915495E-03   12.58086
   2.0770000E-02  2.1080000E-02  9.4979565E-04   12.57021
   2.1080000E-02  2.1390000E-02  9.4979565E-04   13.73922
   2.1390000E-02  2.1700000E-02  6.1608368E-04   20.59069
   2.1700000E-02  2.2010000E-02  6.6742394E-04   20.19833
   2.2010000E-02  2.2320000E-02  6.4175378E-04   17.29009
   2.2320000E-02  2.2630000E-02  7.7010458E-04   20.23257
   2.2630000E-02  2.2940001E-02  7.7010458E-04   20.80261
   2.2940001E-02  2.3250001E-02  9.7546575E-04   15.69687
   2.3250001E-02  2.3560001E-02  1.0524762E-03   14.79374
   2.3560001E-02  2.3870001E-02  8.7278517E-04   15.44578
   2.3870001E-02  2.4180001E-02  7.1876426E-04   20.33553
   2.4180001E-02  2.4490001E-02  6.9309410E-04   22.97302
   2.4490001E-02  2.4800001E-02  7.9577474E-04   16.53148
   2.4800001E-02  2.5110001E-02  6.6742394E-04   20.19833
   2.5110001E-02  2.5420001E-02  8.4711501E-04   11.01403
   2.5420001E-02  2.5730001E-02  6.9309410E-04   16.36614
   2.5730001E-02  2.6040001E-02  4.8773288E-04   20.87878
   2.6040001E-02  2.6350001E-02  4.6206274E-04   24.04782
   2.6350001E-02  2.6660001E-02  4.8773288E-04   14.23564
   2.6660001E-02  2.6970001E-02  4.8773288E-04   28.03264
   2.6970001E-02  2.7280001E-02  5.9041352E-04   20.22164
   2.7280001E-02  2.7590001E-02  7.7010458E-04   14.90712
   2.7590001E-02  2.7900001E-02  4.8773288E-04   22.23153
   2.7900001E-02  2.8210001E-02  6.6742394E-04   19.41210
   2.8210001E-02  2.8520001E-02  4.8773288E-04   28.03264
   2.8520001E-02  2.8830001E-02  6.4175378E-04   14.06751
   2.8830001E-02  2.9140001E-02  5.3907320E-04   21.26785
   2.9140001E-02  2.9450001E-02  6.1608362E-04   17.72928
   2.9450001E-02  2.9759999E-02  4.8773288E-04   23.50657
   2.9759999E-02  3.0069999E-02  5.6474336E-04   19.67549
   3.0069999E-02  3.0379999E-02  3.0804181E-04   30.58876
   3.0379999E-02  3.0689999E-02  3.5938213E-04   25.59706
   3.0689999E-02  3.0999999E-02  3.8505229E-04   25.36286
   3.0999999E-02  3.1310000E-02  4.8773288E-04   17.86865
   3.1310000E-02  3.1620000E-02  4.3639258E-04   23.02082
   3.1620000E-02  3.1930000E-02  1.5402091E-04   48.96592
   3.1930000E-02  3.2240000E-02  3.0804181E-04   32.89181
   3.2240000E-02  3.2550000E-02  3.5938213E-04   27.61559
   3.2550000E-02  3.2860000E-02  1.2835076E-04   39.73597
   3.2860000E-02  3.3170000E-02  2.8237168E-04   24.58883
   3.3170000E-02  3.3480000E-02  3.8505229E-04   25.36286
   3.3480000E-02  3.3790000E-02  2.5670152E-04   30.77935
   3.3790000E-02  3.4100000E-02  3.8505229E-04   28.81682
   3.4100000E-02  3.4410000E-02  2.3103137E-04   30.05302
   3.4410000E-02  3.4720000E-02  2.3103137E-04   30.05302
   3.4720000E-02  3.5030000E-02  3.0804181E-04   22.29520
   3.5030000E-02  3.5340000E-02  3.5938213E-04   27.61559
   3.5340000E-02  3.5650000E-02  2.5670152E-04   39.73597
   3.5650000E-02  3.5960000E-02  2.0536121E-04   42.14636
   3.5960000E-02  3.6270000E-02  2.8237168E-04   40.60225
   3.6270000E-02  3.6580000E-02  2.0536121E-04   28.09757
   3.6580000E-02  3.6890000E-02  2.0536121E-04   45.88315
   3.6890000E-02  3.7200000E-02  3.0804181E-04   25.36286
   3.7200000E-02  3.7510000E-02  3.3371197E-04   23.07692
   3.7510000E-02  3.7820000E-02  2.8237168E-04   20.75143
   3.7820000E-02  3.8130000E-02  4.8773288E-04   17.86865
   3.8130000E-02  3.8440000E-02  2.0536121E-04   33.44280
   3.8440000E-02  3.8750000E-02  4.1072242E-04   26.59393
   3.8750000E-02  3.9060000E-02  2.0536121E-04   28.09757
   3.9060000E-02  3.9370000E-02  1.2835076E-04   64.07233
   3.9370000E-02  3.9680000E-02  3.3371197E-04   25.63428
   3.9680000E-02  3.9990000E-02  2.5670152E-04   30.77935
   3.9990000E-02  4.0300000E-02  2.3103137E-04   30.05302
   4.0300000E-02  4.0610000E-02  3.0804181E-04   25.36286
   4.0610000E-02  4.0920001E-02  4.1072242E-04   26.59393
   4.0920001E-02  4.1230001E-02  2.8237168E-04   27.90339
   4.1230001E-02  4.1540001E-02  4.8773288E-04   23.50657
   4.1540001E-02  4.1850001E-02  2.8237168E-04   30.86403
   4.1850001E-02  4.2160001E-02  2.8237168E-04   30.86403
   4.2160001E-02  4.2470001E-02  3.3371197E-04   32.10656
   4.2470001E-02  4.2780001E-02  3.8505229E-04   21.35744
   4.2780001E-02  4.3090001E-02  3.0804181E-04   32.89181
   4.3090001E-02  4.3400001E-02  3.3371197E-04   30.10415
   4.3400001E-02  4.3710001E-02  3.3371197E-04   30.10415
   4.3710001E-02  4.4020001E-02  2.5670152E-04   34.02785
   4.4020001E-02  4.4330001E-02  2.0536121E-04   28.09757
   4.4330001E-02  4.4640001E-02  1.5402091E-04   42.57776
   4.4640001E-02  4.4950001E-02  2.8237168E-04   40.60225
   4.4950001E-02  4.5260001E-02  3.3371197E-04   25.63428
   4.5260001E-02  4.5570001E-02  2.3103137E-04   34.10415
   4.5570001E-02  4.5880001E-02  4.3639258E-04   17.64706
   4.5880001E-02  4.6190001E-02  2.5670152E-04   34.02785
   4.6190001E-02  4.6500001E-02  3.0804181E-04   30.58876
   4.6500001E-02  4.6810001E-02  2.3103137E-04   25.36286
   4.6810001E-02  4.7120001E-02  2.5670152E-04   36.99218
   4.7120001E-02  4.7430001E-02  2.0536121E-04   33.44280
   4.7430001E-02  4.7740001E-02  3.0804181E-04   22.29520
   4.7740001E-02  4.8050001E-02  1.7969107E-04   31.26410
   4.8050001E-02  4.8360001E-02  1.7969107E-04   37.51119
   4.8360001E-02  4.8670001E-02  2.0536121E-04   38.04430
   4.8670001E-02  4.8980001E-02  1.0268061E-04   45.88315
   4.8980001E-02  4.9290001E-02  1.7969107E-04   31.26410
   4.9290001E-02  4.9600001E-02  3.3371197E-04   23.07692
   4.9600001E-02  4.9910001E-02  2.3103137E-04   30.05302
   4.9910001E-02  5.0220001E-02  4.8773288E-04   23.50657
   5.0220001E-02  5.0530002E-02  3.3371197E-04   30.10415
   5.0530002E-02  5.0840002E-02  1.5402091E-04   35.04383
   5.0840002E-02  5.1150002E-02  1.0268061E-04   45.88315
   5.1150002E-02  5.1460002E-02  1.7969107E-04   37.51119
   5.1460002E-02  5.1770002E-02  2.3103137E-04   37.72270
   5.1770002E-02  5.2080002E-02  1.7969107E-04   37.51119
   5.2080002E-02  5.2390002E-02  4.1072242E-04   30.88604
   5.2390002E-02  5.2700002E-02  2.3103137E-04   34.10415
   5.2700002E-02  5.3010002E-02  2.3103137E-04   30.05302
   5.3010002E-02  5.3320002E-02  3.5938213E-04   25.59706
   5.3320002E-02  5.3630002E-02  3.8505229E-04   25.36286
   5.3630002E-02  5.3940002E-02  2.8237168E-04   24.58883
   5.3940002E-02  5.4250002E-02  2.3103137E-04   30.05302
   5.4250002E-02  5.4560002E-02  3.3371197E-04   25.63428
   5.4560002E-02  5.4870002E-02  3.0804181E-04   30.58876
   5.4870002E-02  5.5180002E-02  2.0536121E-04   33.44280
   5.5180002E-02  5.5490002E-02  7.7010453E-05   54.61187
   5.5490002E-02  5.5800002E-02  2.0536121E-04   38.04430
   5.5800002E-02  5.6110002E-02  2.5670152E-04   30.77935
   5.6110002E-02  5.6420002E-02  1.5402091E-04   42.57776
   5.6420002E-02  5.6730002E-02  2.3103137E-04   34.10415
   5.6730002E-02  5.7040002E-02  3.3371197E-04   30.10415
   5.7040002E-02  5.7350002E-02  4.3639258E-04   27.35868
   5.7350002E-02  5.7660002E-02  3.0804181E-04   28.09757
   5.7660002E-02  5.7970002E-02  4.1072242E-04   19.44967
   5.7970002E-02  5.8280002E-02  3.8505229E-04   30.39698
   5.8280002E-02  5.8590002E-02  2.3103137E-04   37.72270
   5.8590002E-02  5.8900002E-02  5.1340303E-05   68.82472
   5.8900002E-02  5.9210002E-02   0.000000       0.000000
   5.9210002E-02  5.9519999E-02   0.000000       0.000000
   5.9519999E-02  5.9829999E-02   0.000000       0.000000
   5.9829999E-02  6.0139999E-02   0.000000       0.000000
   6.0139999E-02  6.0449999E-02   0.000000       0.000000
   6.0449999E-02  6.0759999E-02   0.000000       0.000000
   6.0759999E-02  6.1069999E-02   0.000000       0.000000
   6.1069999E-02  6.1379999E-02   0.000000       0.000000
   6.1379999E-02  6.1689999E-02   0.000000       0.000000
   6.1689999E-02  6.1999999E-02   0.000000       0.000000
--------------050008020502040502000001
Content-Type: text/plain;
  name="Prot_60_usrbdx_90_tab.lis"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
  filename="Prot_60_usrbdx_90_tab.lis"
  # Detector n:  1  OneABin   (integrated over solid angle)
  # N. of energy intervals 200
   0.000E+00  3.150E-04  1.185E-03  1.969E+01
   3.150E-04  6.300E-04  2.134E-03  6.014E+01
   6.300E-04  9.450E-04  1.996E-03  1.763E+01
   9.450E-04  1.260E-03  1.614E-03  1.439E+01
   1.260E-03  1.575E-03  2.176E-03  1.235E+01
   1.575E-03  1.890E-03  2.372E-03  9.758E+00
   1.890E-03  2.205E-03  2.429E-03  9.030E+00
   2.205E-03  2.520E-03  2.015E-03  1.228E+01
   2.520E-03  2.835E-03  2.623E-03  1.404E+01
   2.835E-03  3.150E-03  2.835E-03  1.003E+01
   3.150E-03  3.465E-03  4.083E-03  2.403E+01
   3.465E-03  3.780E-03  3.483E-03  9.203E+00
   3.780E-03  4.095E-03  2.690E-03  1.024E+01
   4.095E-03  4.410E-03  3.392E-03  8.824E+00
   4.410E-03  4.725E-03  2.507E-03  1.143E+01
   4.725E-03  5.040E-03  3.685E-03  8.255E+00
   5.040E-03  5.355E-03  4.204E-03  2.035E+01
   5.355E-03  5.670E-03  2.843E-03  1.834E+01
   5.670E-03  5.985E-03  2.804E-03  1.040E+01
   5.985E-03  6.300E-03  3.003E-03  9.662E+00
   6.300E-03  6.615E-03  2.584E-03  1.278E+01
   6.615E-03  6.930E-03  2.811E-03  9.657E+00
   6.930E-03  7.245E-03  2.747E-03  1.149E+01
   7.245E-03  7.560E-03  2.413E-03  1.522E+01
   7.560E-03  7.875E-03  1.967E-03  1.480E+01
   7.875E-03  8.190E-03  2.256E-03  8.376E+00
   8.190E-03  8.505E-03  2.234E-03  1.395E+01
   8.505E-03  8.820E-03  2.049E-03  1.369E+01
   8.820E-03  9.135E-03  2.057E-03  1.048E+01
   9.135E-03  9.450E-03  2.438E-03  1.932E+01
   9.450E-03  9.765E-03  2.361E-03  1.264E+01
   9.765E-03  1.008E-02  1.665E-03  1.578E+01
   1.008E-02  1.040E-02  1.882E-03  1.509E+01
   1.040E-02  1.071E-02  1.832E-03  1.398E+01
   1.071E-02  1.103E-02  1.442E-03  1.788E+01
   1.103E-02  1.134E-02  2.015E-03  3.167E+01
   1.134E-02  1.166E-02  1.572E-03  1.161E+01
   1.166E-02  1.197E-02  1.402E-03  1.365E+01
   1.197E-02  1.229E-02  2.012E-03  1.476E+01
   1.229E-02  1.260E-02  1.397E-03  1.296E+01
   1.260E-02  1.292E-02  1.720E-03  2.262E+01
   1.292E-02  1.323E-02  1.224E-03  2.576E+01
   1.323E-02  1.355E-02  1.292E-03  1.253E+01
   1.355E-02  1.386E-02  1.154E-03  2.193E+01
   1.386E-02  1.418E-02  2.020E-03  1.864E+01
   1.418E-02  1.449E-02  2.203E-03  4.437E+01
   1.449E-02  1.481E-02  2.411E-03  3.464E+01
   1.481E-02  1.512E-02  1.524E-03  1.759E+01
   1.512E-02  1.544E-02  1.327E-03  2.222E+01
   1.544E-02  1.575E-02  2.185E-03  3.155E+01
   1.575E-02  1.607E-02  1.188E-03  2.177E+01
   1.607E-02  1.638E-02  1.229E-03  1.731E+01
   1.638E-02  1.670E-02  1.363E-03  1.802E+01
   1.670E-02  1.701E-02  1.259E-03  2.121E+01
   1.701E-02  1.733E-02  1.470E-03  1.959E+01
   1.733E-02  1.764E-02  1.243E-03  2.203E+01
   1.764E-02  1.796E-02  1.337E-03  2.556E+01
   1.796E-02  1.827E-02  8.795E-04  1.609E+01
   1.827E-02  1.859E-02  1.160E-03  2.790E+01
   1.859E-02  1.890E-02  1.411E-03  2.023E+01
   1.890E-02  1.922E-02  1.277E-03  2.373E+01
   1.922E-02  1.953E-02  7.235E-04  1.711E+01
   1.953E-02  1.985E-02  8.775E-04  2.585E+01
   1.985E-02  2.016E-02  1.008E-03  2.564E+01
   2.016E-02  2.048E-02  1.165E-03  2.755E+01
   2.048E-02  2.079E-02  2.202E-03  3.038E+01
   2.079E-02  2.111E-02  1.165E-03  3.859E+01
   2.111E-02  2.142E-02  1.107E-03  2.732E+01
   2.142E-02  2.174E-02  8.506E-04  3.070E+01
   2.174E-02  2.205E-02  5.917E-04  2.587E+01
   2.205E-02  2.237E-02  4.294E-04  2.135E+01
   2.237E-02  2.268E-02  8.937E-04  2.101E+01
   2.268E-02  2.300E-02  9.498E-04  2.626E+01
   2.300E-02  2.331E-02  7.094E-04  2.438E+01
   2.331E-02  2.363E-02  1.180E-03  1.972E+01
   2.363E-02  2.394E-02  7.131E-04  2.487E+01
   2.394E-02  2.426E-02  7.712E-04  2.972E+01
   2.426E-02  2.457E-02  8.793E-04  3.085E+01
   2.457E-02  2.489E-02  5.747E-04  2.575E+01
   2.489E-02  2.520E-02  8.607E-04  3.476E+01
   2.520E-02  2.552E-02  9.217E-04  3.108E+01
   2.552E-02  2.583E-02  5.841E-04  1.735E+01
   2.583E-02  2.615E-02  5.110E-04  3.182E+01
   2.615E-02  2.646E-02  4.269E-04  3.425E+01
   2.646E-02  2.678E-02  4.243E-04  1.912E+01
   2.678E-02  2.709E-02  3.828E-04  2.818E+01
   2.709E-02  2.741E-02  6.035E-04  2.607E+01
   2.741E-02  2.772E-02  9.383E-04  3.572E+01
   2.772E-02  2.804E-02  5.847E-04  3.349E+01
   2.804E-02  2.835E-02  1.351E-03  4.564E+01
   2.835E-02  2.867E-02  5.081E-04  3.007E+01
   2.867E-02  2.898E-02  6.550E-04  2.388E+01
   2.898E-02  2.930E-02  5.966E-04  3.216E+01
   2.930E-02  2.961E-02  6.123E-04  3.209E+01
   2.961E-02  2.993E-02  7.379E-04  3.130E+01
   2.993E-02  3.024E-02  4.947E-04  2.966E+01
   3.024E-02  3.056E-02  4.197E-04  2.553E+01
   3.056E-02  3.087E-02  7.920E-04  7.032E+01
   3.087E-02  3.119E-02  5.940E-04  3.969E+01
   3.119E-02  3.150E-02  6.655E-04  3.016E+01
   3.150E-02  3.181E-02  5.969E-04  5.695E+01
   3.181E-02  3.213E-02  3.598E-04  3.184E+01
   3.213E-02  3.245E-02  3.486E-04  3.768E+01
   3.245E-02  3.276E-02  1.521E-04  4.009E+01
   3.276E-02  3.308E-02  2.564E-04  2.824E+01
   3.308E-02  3.339E-02  6.495E-04  3.803E+01
   3.339E-02  3.370E-02  2.869E-04  3.439E+01
   3.370E-02  3.402E-02  7.375E-04  4.965E+01
   3.402E-02  3.434E-02  3.109E-04  3.356E+01
   3.434E-02  3.465E-02  2.590E-04  3.176E+01
   3.465E-02  3.497E-02  2.134E-04  3.696E+01
   3.497E-02  3.528E-02  5.334E-04  3.471E+01
   3.528E-02  3.559E-02  2.405E-04  3.812E+01
   3.559E-02  3.591E-02  2.819E-04  3.641E+01
   3.591E-02  3.623E-02  2.249E-04  4.873E+01
   3.623E-02  3.654E-02  1.899E-04  3.579E+01
   3.654E-02  3.686E-02  2.743E-04  4.018E+01
   3.686E-02  3.717E-02  5.295E-04  3.518E+01
   3.717E-02  3.748E-02  4.049E-04  4.092E+01
   3.748E-02  3.780E-02  5.981E-04  4.383E+01
   3.780E-02  3.812E-02  6.304E-04  3.595E+01
   3.812E-02  3.843E-02  2.318E-04  4.959E+01
   3.843E-02  3.875E-02  6.063E-04  3.379E+01
   3.875E-02  3.906E-02  1.734E-04  3.434E+01
   3.906E-02  3.937E-02  9.438E-05  6.994E+01
   3.937E-02  3.969E-02  3.466E-04  2.752E+01
   3.969E-02  4.001E-02  3.627E-04  5.150E+01
   4.001E-02  4.032E-02  3.066E-04  4.232E+01
   4.032E-02  4.064E-02  5.097E-04  3.199E+01
   4.064E-02  4.095E-02  4.585E-04  2.779E+01
   4.095E-02  4.127E-02  2.253E-04  2.980E+01
   4.127E-02  4.158E-02  6.788E-04  2.829E+01
   4.158E-02  4.190E-02  3.397E-04  3.651E+01
   4.190E-02  4.221E-02  3.448E-04  2.554E+01
   4.221E-02  4.253E-02  4.515E-04  4.586E+01
   4.253E-02  4.284E-02  3.671E-04  3.249E+01
   4.284E-02  4.316E-02  4.431E-04  3.215E+01
   4.316E-02  4.347E-02  3.512E-04  4.715E+01
   4.347E-02  4.379E-02  3.523E-04  3.106E+01
   4.379E-02  4.410E-02  1.722E-04  4.997E+01
   4.410E-02  4.442E-02  4.650E-04  6.052E+01
   4.442E-02  4.473E-02  3.006E-04  4.409E+01
   4.473E-02  4.505E-02  4.639E-04  3.429E+01
   4.505E-02  4.536E-02  4.040E-04  4.225E+01
   4.536E-02  4.568E-02  2.489E-04  3.930E+01
   4.568E-02  4.599E-02  2.838E-04  2.568E+01
   4.599E-02  4.631E-02  3.332E-04  2.924E+01
   4.631E-02  4.662E-02  3.656E-04  3.140E+01
   4.662E-02  4.694E-02  2.619E-04  3.456E+01
   4.694E-02  4.725E-02  3.757E-04  3.086E+01
   4.725E-02  4.757E-02  4.677E-04  5.674E+01
   4.757E-02  4.788E-02  2.792E-04  2.986E+01
   4.788E-02  4.820E-02  2.070E-04  2.987E+01
   4.820E-02  4.851E-02  4.288E-04  4.767E+01
   4.851E-02  4.883E-02  1.676E-04  4.970E+01
   4.883E-02  4.914E-02  1.583E-04  6.110E+01
   4.914E-02  4.946E-02  4.267E-04  3.770E+01
   4.946E-02  4.977E-02  4.462E-04  2.896E+01
   4.977E-02  5.009E-02  4.153E-04  3.049E+01
   5.009E-02  5.040E-02  8.574E-04  3.604E+01
   5.040E-02  5.072E-02  2.692E-04  5.261E+01
   5.072E-02  5.103E-02  1.580E-04  4.190E+01
   5.103E-02  5.135E-02  2.897E-04  4.811E+01
   5.135E-02  5.166E-02  2.601E-04  4.399E+01
   5.166E-02  5.198E-02  1.248E-04  5.996E+01
   5.198E-02  5.229E-02  5.286E-04  3.411E+01
   5.229E-02  5.261E-02  2.857E-04  5.274E+01
   5.261E-02  5.292E-02  3.639E-04  4.534E+01
   5.292E-02  5.324E-02  9.628E-04  5.255E+01
   5.324E-02  5.355E-02  4.037E-04  3.564E+01
   5.355E-02  5.387E-02  2.920E-04  3.895E+01
   5.387E-02  5.418E-02  1.588E-04  4.206E+01
   5.418E-02  5.450E-02  3.066E-04  3.634E+01
   5.450E-02  5.481E-02  5.737E-04  3.010E+01
   5.481E-02  5.513E-02  2.064E-04  4.423E+01
   5.513E-02  5.544E-02  2.217E-04  8.931E+01
   5.544E-02  5.576E-02  3.659E-04  4.645E+01
   5.576E-02  5.607E-02  3.957E-04  4.617E+01
   5.607E-02  5.639E-02  2.189E-04  4.062E+01
   5.639E-02  5.670E-02  2.862E-04  3.392E+01
   5.670E-02  5.702E-02  3.149E-04  3.571E+01
   5.702E-02  5.733E-02  5.289E-04  3.250E+01
   5.733E-02  5.765E-02  3.724E-04  2.701E+01
   5.765E-02  5.796E-02  3.433E-04  2.647E+01
   5.796E-02  5.828E-02  2.241E-04  3.104E+01
   5.828E-02  5.859E-02  1.101E-04  3.917E+01
   5.859E-02  5.891E-02  2.156E-05  6.883E+01
   5.891E-02  5.922E-02  0.000E+00  0.000E+00
   5.922E-02  5.954E-02  0.000E+00  0.000E+00
   5.954E-02  5.985E-02  0.000E+00  0.000E+00
   5.985E-02  6.017E-02  0.000E+00  0.000E+00
   6.017E-02  6.048E-02  0.000E+00  0.000E+00
   6.048E-02  6.080E-02  0.000E+00  0.000E+00
   6.080E-02  6.111E-02  0.000E+00  0.000E+00
   6.111E-02  6.143E-02  0.000E+00  0.000E+00
   6.143E-02  6.174E-02  0.000E+00  0.000E+00
   6.174E-02  6.206E-02  0.000E+00  0.000E+00
   6.206E-02  6.237E-02  0.000E+00  0.000E+00
   6.237E-02  6.269E-02  0.000E+00  0.000E+00
   6.269E-02  6.300E-02  0.000E+00  0.000E+00
--------------050008020502040502000001
Content-Type: image/gif;
  name="Proton_Fluence_usr_bdx.gif"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
  filename="Proton_Fluence_usr_bdx.gif"
R0lGODlhgALgAcYAAP///wAAAKCgoP8AAADAAACA/8AA/wDu7sBAAO7uACAgwP/AIACAQKCA
/4BAAP+A/wDAYADAwABggMBggACAAED/gDBggIBgAEBAQECAAAAAgIBgEIBgYIBggAAAwAAA
/wBgAOOwwEDAgGCgwGDAAGDAoIAAAIAAgGAggGBgYCAgICBAQCBAgGCAIGCAYGCAgICAQCCA
IICAgKCgoKDQ4MAgIACAgMBgAIDA4MBgwMCAAMCAYP9AAP9AQIDA//+AYP+AgMCgAMDAwMD/
wP8AAP8A//+AoMDAoP9gYAD/AP+AAP+gAIDg4KDg4KD/IMAAAMAAwKAgIKAg/4AgAIAgIIBA
IIBAgIBgwIBg/4CAAMDAAP+AQP+gQP+gYP+gcP/AwP//AP//gP//wAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAACH5BAEAAAAALAAAAACAAuABAAf+gACCg4SFhoeIiYqLjI2O
j5CRkpOUlZaXmJmam5ydnp+goaKjpKWmp6ipqqusra6vsLGys7S1tre4ubq7vL2+v8DBwsPE
xcbHyMnKy8zNzs/Q0dLT1NXW19jZ2tvc3d7f4OHi4+Tl5ufo6err7O3u7/Dx8vP09fb3+Pn6
+/z9/v8AAwocSLCgwYMHAyhEyLChQ08KI0ZaiCmiREEUH2rcyJERxYyNQFYSiTFAx5MoT340
adFkyZUXAVh8SXJQzYwradJ0KTPmzJRAg45rWbJozp6Ech41RNQmT6UuF0I1GpWn0KtYr4nc
WrWoU6pemVr96lUq0alJm2Zdy/YZ10L+R0GirYlULFmzcLu+bcu377K9d/WyFBw27yGchJHO
HVzYr+PHuwA7lYs47liyky82/Wm58mXIoEOLHk26tOnTqFOrXs26tevXsGPLnk27tu3buHPr
3s27t+/fwIMLH068uPHjyJOHo+xZufPnaZ8mpgu9OnDGZbFTt869N3PBn3u2HE++vPnz6NOr
X8++vfv38OPLn0+/vv37+PPr38+/P/9fiFGFHVPOhFeMgQc2gyAxCwrToIMASgdWXQQqaOGF
fxWIYYa+BKhYV4dpyMyDwZAIoIgcjtjLeDq1GOKGyZh4ooowIiNjh9nc2N2O6+h4jI88BnnY
dtEAmaCQsqn+5QiRS6rFZCwDdGIkg0gmCWKTmdw05SVROtJlI18qsiWEVcIGlU+a/SReVTHZ
5WZdMLHZnIuGhLmInYrg+aJWZZrJ2XRXwqldeEp6iNeHEx5K1wCMNuroo5BG+mhIOfb5mmSC
FrYYIloGtqaiA+4Unp6IkHqIqUlVamlrmFqWVqKEfmZoqJlieuojqBaSK0aqrrpaqyDitemQ
xGZXK6CZcillr76mZqua4mVmbLGZxSnntGpCO8uYJTbrLTncAhPut+QaMy6O5aa7zbkrqusu
NuzyEu+79Noyry731qsvLPni0u++AKvyr70BF2wjswYnLC7CCjccGcMOR0wwvBL+V+wvxBZn
3MrAtXCs8ceLeLwtyCSzIrIsJ5escsr8quyyKSy/EvPLGc+8Mc04g2KzyTn3zMnOqwDts8FC
p1L00AAffYrSSNfLdClPN+1u1KNQLXW5VoeS9dXebm2YtHTCxfXYHnU4a7IVkq122vhSuBdd
Xq9dZdyYbTZYrG166Wiee8udrpKRSfghrTZhsqsgh/u9Kt1hDcs2JYknrnifjLuN7J6VRD45
vXSzuCbYmEPu5ebvVq7JmJqTrq7pWRo+uup/YwxJ6rCTy3pFroNZO9ay4/r67l33/vudwH97
+yWoD1+8pcdbkrzuy/va/Ei5MyJ59NBNT8nz1mO/uPD+0HfvPeXgi2/++EhqPwn3xKM/d/nt
n+8+j+pLwn6e86cPP/7h579j/ROpXvz81x0AQuJ+ibgeAYdjwEcgsFQL/N/+Eqi8CD6ngUsS
IP8sWB0MUkpZ/eOgczxYNhDKT4TJIWHINEhBFF5wghAMoQuPo0IxsTCGM0whDG8lwxwSp4aJ
eCAPfUhDUpiHZyYcIBF/aMTyINEStFsiE6sWuqXdcIhSZKARqwizK9Ypi8UBIqe8qCswTpFP
SdygGa+DiqjcLI0tXCMbYXa3J2augnLcjcjcODMhfjGPv9mjTHg1EVlda4xwxCEgebPHQzmQ
JBJZitgSicVF6uZkbjwghXT+IslUUfKPltRjhF4lLLhpy3c9DGVsAFc1lqxPcGYhHCGhiEdV
0qaRsgyZdGbSycJ9soy2zI0gM3lAQ6JtkrRMZTBXaUVihgRwcQoiGQmhwGWeBpd5s+IvqWnN
22Aym9pM5gm7aaYdglKJ5Cxn0OwoOmWm81fsDOcd3flOZzWRPPGURBTr6ZqBHXGd2xxENfkJ
GX9yEWrTFChB1SmwfEZinwtVjRgPOrtaRtQ0mJxlQ7cZqYtKdGmBamNCEedR1gjSlxsVJzpL
WpqTalSkASUpS+1Jx1N2MaYAGOhMs/LNJ+lspDndKWom+rhJQFSooemjQytKT6T2RakAVaka
nZr+1CZ+CqpGhdRKqeoYf94NnAjFqU65ihKXDjKq8xwnWZ+6VCqKda1Vpdhb4VpQtBoNqGOl
q0Zc6lOIeFGreu2qPLGa1hjmNbAN4ZgT7WrU9h0WsQhRrOcYq0/HQpatN5XZNLv02MsWhKjI
LCwPO+vZgWRUsxyVaWl5CtKWpTaoq2VtZl3hR5mSNrb/MOsbpfrH2+K2H9js67Je69vf7iO4
OvoO6IpaWf4V17j5YByaYIVI3gLzudC9R0+FOyDHhRZyfbtudoGCXPsZM5acsqnuTMXZ8XKE
lVo7a3LHkqZcJo+9qnXve5v50me+ylMU9R1+Yavf/XZRvWISHICZ+1D+Ave2wB1hmXCvSuHP
VfeOA8YuhOEBWk/SMsMb3khPaes6EIf4IcHdLYYR17f2ntghfFWx6MIUpY6+OLHhtNl9ualQ
B9/YIMidcOuURWMea/jH6QiyhEvMY9siObLmLGORe3zkJ5ujwyhdcY+dbOXPRpmbU+Zyl037
ZYWGGbZVHrM4sNzfyp7ZxWoOCJvP+uEmoznOApnzjrd8ZzwDRM9M5jOc/QzcsMZiz/kdNKGP
a1WUBTrR+V20Pgza5rtCUatTTrOkteFVfFK2wW+O9Ka1CzVPp1TLkPbxqO1xWhIT2c6KXnU9
AP1qQYta1vPY7sbCC2pY3xrX8SjvOvNa45z+YvrXwH6HWZVKbAdzdlLJnvVgedbsMwuU19Fu
h7AbWm07G7tR2ZaHrqldZz4DM9zKrlS3ze1tdPdI3eVGtqbdnSJrmGTd8mYxuHsBbXqLAtD4
VrWtgzFvfwM83gJPNTAKTu+DozrhzhYGw93t8Bm3e+AL97eh7Q3x9V5c4b+YOLpbfbOADxjZ
uRB5uBvpPKtkC6yDNHkcO44LlWdbt4WE5aDGKHNF2hynGneePAN4tl7yitgtnvnPrRv0n4EU
wUGUUCTrOKQBCFnUJ6f5LZY+avjqzOsJnszgNpmqZp/7nFxvZ9N/ajZPGT3mD+03uzGuC8Cu
/RPfpHMx3R7Swg3+NNZZTztTBb9qXJJdl3+aU4j+Dm25XxvbdUf53YV+4Epv8e+NDeqz922L
WE9+yAh1pkgx39xiP55RnZf859f3dJhTkfQN1jyYa6761Qew9SXPfOxNT1LCj9b2TqcY7Acf
cd8/GPina6PrX6974hfZ+OJFPujvufz4Ss7xo+so72vheenbb9rD7nXpi097rXu/hBtntvhP
D+Lnlx/6i+70ZO+qwO3P3cjk33rtz2/DVpqa/uv3cVSWf6lnfvwnTcIXgPc3gO6nfwZ4gAEG
EVD3bw/Ye5X0YMemN5wXCt0HgehHfVfXOvUncLnSXiZYgSg4T/BHaIZHYiOoJyVogbL+N4Kl
0IEe2H90pFkveIHidYI0SAo2eIMX5n/V91M7eE6/54OxB4T7J4RZRoQhiDvXY2O/pn0W+IOj
EIRO+F3UcG8KKHknOINLmIWJBnmdgH14tm0AWFEwaFk1ZndNpQmAdwor6FFqOHpfiGa8ZoX2
p1aeMIemUIcXRVh4yIZYhCdhWGxwqEiiAIg1KGuEaEVHeHaIKIOKuH8F54hMWHitF4XIg4J9
GHE9eIWYGIhYR4eQ2Im5Z4hoh4H6hoY8+AmaSIacmICsSImuKIYGmImnaIq1CIVFiHegSILt
tnngZnpHJYu9+Ii/+G//V4gCdojHR4q76IuV6Itdp3w744X+t/hte3ON7qeFTQhC18iM2Vhq
QsON0diK0UeNyfiHy7iJ5wiMnih0k+hr7aiLkWOGmzCLjZiKIBiMfoUqjUeM7DZol1iNmuaP
HAiQfLIrb4h91tZn7ph9m1iO8rhpFad0BtmL4ViKBMmP0IORtDiP8cVgl3cqcheKE4mQsqeQ
fvg6JPmPzch2XNhFJiaDuNiDNqY5EEmOSEiTJvl1z/hKikcgOfmSO4l/ZkaRxPOTyTSTDVmT
eLdYI3E5SNlCDbiUA9iU4qiUM5d5UgkKghhRzeNdZaeV1IRtLVmG+0ZjbKl1h8WQZOmQDpIm
0VJ1NeGSPjeNBMaQoXh2YhmUU8n+ghOYC2dyeEenlgvYkgUpZnMXmAIok4RZl1TZdoiimHBn
WAIYkeHVhpB5a5K5gJQpmEKpkQyClWKTlIFHmF8CmHIJlKZZmKiZmke5moyZbwc5RLAJlSpY
mcp4mdKgjn7ZmrP5mvHYlFuFSmMZnEPZhT72mLqZnKrlYiu5iMCpgZ+JisJZJMRohX15nN6I
esopZTGZQEk3m5b5nM6od2tonqKoK+l5ca95bG3om25GnbQpaYZXj9RznyA3VR1pawC6nN1D
l85Zm6Fnecz3RXx5QqDpawUqoOvYnPDYnX4lepJYJ3AIURFKoLeCn+JnoWdol8BIbg4aoGE5
oJCGiNj+SZpPqZ/rqaDoqGPjmYiTWYwXuJWzN0Cw2KMkyglluVAO14ATyYgsWnwFKaIwqWoM
N6QEdVo26pIP6qNxiYtUCqEW6aTcyZ5f1xgbmn9HGotBaaQ9KqA+KaMJyp/gx20tqqJIiqRm
Wp4rmqZcio006lYMSoFZGp8GqqPKmYFouqVB2o8m6oyHaYRveqPkSaHyuYGSkpCjqabwWahy
eKiIOqWLOqnZGaMviYyxqZ4peqfmmKfQ2aecKqoHyqh3xqR2aqmZAKX8tF3qJ6YP+KP4g6Nh
OKgeR6oZyaY1JZDLQqW4WqK6aolXmqOB2qlCiqk26YKLOo61Bqq7mqSvyqz+hoqhw6WDvWd3
T9qqluinu2k92Pmt2vozYNegPPqtxwqWHwqSlCqHL6pmdximcPmWqNiu2/eueriBy/qK/rqm
8YcNVreS2Hqp+iqu8YhfCHqhXvqlmqqw7Bqu1cqvY/qpBxurzlqVuZSSEpuvFKuT1tqZTimt
TJeG/OWe9gqoj5iwR4qRDBuvzXqum+BKKouTKjqxGFuxO0qyDVuiNHs66cqnOZuvkxKp5AlY
MEtBPzuzD8uxEduWQLiHgtqV1Vm1AIt6sApUT1akLIunOJqPSwtBTZutTyslQ2t9CkeFJQmu
dyq1+Va2lxq0WZK2ivqxqJBpJeuUYxuiMmu2pjr+W6OHtSmYRGGLf1R4saCqqkB7tifppnjb
pW4Ljl8bs77atsBao/65PXybsVySgZ7ZqJ3LlZXquUCHZCmGonALtiIrtuzot5d7mpl7l7eZ
lpHLumA5nkXbmHtrrp0nkqEUiYVEXYu3u9xJuGzbt6O6taerCrI6NsK7d/8FN31msHkLu6U7
rnH7t3Nbfu8kf3abYHgJVgphdQEQg9c7qtm7sEzLvRqbcuQUvi1XlMM7dtQbhCI3oUBaue0b
u/vJfbO6RSiJJZl5v2CYCvr7r+xLtu7LtavwvFwTvY9EvFmpu6Z7Q463uAsMu8x7slASwLa5
SwbCjZKZv2Sqwb5quR3+LFoAXE+KNTLVyrgO67PGC3EobLIsTAsQfDUvfLNrWMIIfMKsqrzr
C6NO64DfG1bp6LZkOqO8W6VwGq6qWqxx2Ao7LDU9vMQxrKyNS8PONp823LvniUrem06SNX+F
yJLA28UWXMT76sXSuaJqh8RmrLnC6nRbbMUjCqRIu7wKe8FxOgtX3DRea3xHSLhPLLJMOscF
mMS2mKp42qvUKscha8SkNaSDjDSFLMh7nMfb667iQ8VNDAuZPDSbLAuHzMTAOcmWzMWBjMog
/JCtywqiXMlAzMCg7KM5zMkuDG+5y2/2OcuRmXQXS3dUjMmxzHHzym+TG8bOdaPnY6nI3Mv+
tojDsNyuo4zNAFvDfBOwpJzMIEhubOsLoOvMapnHoSbDwvzB1GzHKGrN7FzJ2byk0niQwAvJ
zgvOFNgyxbwL4ViuuvzLflqJPymIpewzOBd+Rgy/y5p19cyS+Him/dvIdRysm8t6/Rx5DW29
HCnQ6SzRuEzR8duJtQrP32zBi1usevuhbVnQZTzS1Vy4sDyKBBieEE136wzSvOzIwbqKJv0K
lEuskLfS9XyKLk3HMO3OCv3TrhDUm0q2iNyh2wyD47zTFZ1+7yzT8cyAT/3Q3ykps2ecIp0b
tSxts1XSWn3SiLupe2hkcXm4uljULewbB90tFo3WdV2fR/vUY2v+jAnrlR091rpR1wuDe1md
17lopn0toVJMjXKtw8BB2DMiywv90s3c0OwYanCdkI9Nrt4MuLwh2ehiaQhcsBnN0Gttq4vN
fqz8keE5hk68G6LdLpl10d932r8LqVm7v6tty517yzi8kJHtZWd92MxQpayMj4rd2FXtlp8N
yLAt28RNj7Va1v7crX0cuqn2fMFMfsY5hb8q3QRxxol6t8+A3OCazsutecH83Uwlu+JNZoKb
1c2gtGv5jQPd3a0d1h0N3pgb33nmy+ct15wtiuuNwgW+k/4ti8tMG7MtL9oYsQP+unFt4Lzd
2BWu4O/9v3eyxpDx4A9z1zV7SMULDSz+zdrIabU7i+EJrtwbHttV/Bgg3jbFXbc7V+ITzpUF
rreYjbEo3tsxvst/GhozjphtOr83XsE5HtEVzuO+naXuTcYcPsagUeQXI+LIwxjopZcqdOKA
PXAH/tvmvLeOesRSjhpWvi3lLbT+mZgGbOIEnsuQud9Q/qii+5VMbed47soyPt0lchZ9t5hw
TuE7vtteWc5fDp+hOeTNN3xo7ucCjLbHZLtLbsxyztgzKNUQ7eQciq8d/tyka1GjkeYTU2pV
GTaUftzJyt5K+9Zgbd+JLtH468oLXrik/g23TgsJDbnV0McZLo08693YS+Zg1uAViNvq/OGQ
vs+0zCjnK7r+3tDirRjsoRuwSAvtFF7mnR3kfpHrI4PlvA4O0o6livym+jvrlV3rjn4a3u5o
4P6e4u7h6z3nHm3havSO207lDWHdlU1qMY3rzd2i1uvaik7sOO3rw7PuKTFW7X5oOcjP45Di
9z7Q5r7bmWaGL6vTE99c4AXq9MDwyw61EC8OEs+R6n3fd17uAg2nPBrNzcfI+wDy463EIx8O
Jd+XJ79lRI3ir43pCf/yQK8PMi/fbqXFEd/KdMrpJU/tY37i6h700S30Qs4P5L3mw3r0qreu
TFmdqe1kdsLRPt8/Ch/1+TD0AR6Qth1ADf87IYmvOc/1XB33Be/0Gz72Lx7zU8/+aP8e7f2u
8vS+3RZ+7sBe91B/91IP87kl4OLe95f+oDvv9cPu1y0f0GRv+GWf95Om+LjO+G/840uv8pvu
4nTq8hyP+R9v+viwkZsvrXG80i3mrcLc+fT+9KWP+Jdv+4W2991w89EMxt3H9IKv4vj8x+zH
55SMD9UU8Jmv+9zA+3/K6V2vyln7mJO/8bzqwfZg9nJ2Dc6O7bvP+ejp6aQC/MPuurR//agv
D9r/Z5rP901Yy78P+t+Zm8rP5HLqwPCw/onP/Nvg/Nk6/oAAADAwKChIWIhoeFi4uIiomOg4
eThgedlIaRjJ+HjJqZkZOkpaanqKmqqqCYraugobKzv+S2sYUIube3r7quv7S6nYiwvJyhhZ
3ErYW7wpWooMPamcOgx8jU0sa53d7R17+y2OyjtuHvu53N28ma7uLDn9/Nh5XE1PCkrtet7v
78gtVMB/BLOFK+ivHMKFwT5hYwfQHrx3ldJ5wiSKoqlo+eR5bIhxIMORqkR+JIny18GU3RSy
fJnL3atk9GTWtJjIXbxTHEfpC9YRXi2ZJmHSKorPqNJVK5fqCuDQqdRqOPHhpCkP4rtlOnkm
FeixK1ih2uatqjqVbMm0bDU1bQsOrtyxAKP+3FmzXrxoJnvSbWc3qN6hg7fJtbl2Lty3isk1
fmxW4l2NzrQ24ntPLdCIewX+R0bXCS0/xUg1QzbK+PSo1KqdMrsc1uzWvBk/by58MjRsn18J
1wZN2nDrpayHOypu/OVriBUxRrQMnapdopZiMz9ZeqzfxHOz206OEDl48eBJLi+8XTLgqBK9
4sUqTOd1zriFpx/dXXj5lOSN99+P0Hnz3DdbZXpRBlJ176nHWYF/ebfZfe7lBxyAJP3XGoYW
+iPgXbg52JmBG8HGEVdVwdcRYkcdqJ90Cv7jHYQb+qLhaTXOOM5yaCFz1W6SWTTQbCX6OBGR
t32XmIQjVgiLaPidJRB7OF5z42NVTvmQlCLWFRKL0XUW5GWBhWiVkdj5YhmT3LUIpZq5YakS
nAD+XCmnOCryltGOJLaXInoMljnfV8poyWV1aTbJJqKJPnlknU/JSaejdiKp2yDUybfniNeh
2KCZ9H0WZqX1LanoLBBaIxqqJREq6ZyQtqoUgsaASpuQzVGz6Z9FBkrWa6QeGqVzlM5qarHS
vPkpo7AKEml4ATxry7OsNbvsNU7GRiyYx1zLYI+5kpmtr8cSGaqlV7HaW6nq/mVauqRWy6xi
0M55y7zzTkItvGhmJq62/gL1JYjbctsru3QBi6e7Ack4bLYGf1caw8nlm1A40NqbGsX6luVe
v7aai0mH64XsJ7gXOdRvwjaVa9rCt67ZpmftzpywvhoTdHG99x4n7c3+G8e8pMeZYmZMwGJ2
afKZDh/M4rsfvhtxody4jKy748tbcq5ZUJ2t3461qLc7O
9Hbt999s3eklyMLKyra/GiUDt6Ca0/c4K2N+OpPVeA/Ydu6qAg0v599kbTHXoYve1uCAUkd7
rZnqtvZvuRueOm8EWs1j4JgP2OPV7d55rcT+vUo8XB5PF5isfn28uq5HpgwWiIN2Of32ZSqL
u+mX2z4q+8v6jlr44o99PXXIrTYoIpynXLf+tKKRa2zxq8+Q6DcZ+92ue1J7naP4BxMM+u8X
J0KeZkjnPNWRrXnpAeGtkJa/ssmGevPL2wf55Dai0Cor02jd5sC3wamgboKIQ2EPSWZATgHw
eTSM3cDYM0A9FQmC1skM8x62Q+8NR4MsoWIOtUG35yywcigDlPqURkQtGpGAd9tSwRq2wk71
xYsJjCLYrIgSOF6xWFnMCySUWLqV/ShVtguT+xj4MsHVT1QxbBtW8nZIZEXReliS40gcOUdE
1dF0ZDsjnwpojyRa0EBjLFwghXSp7qEMbhFco+KgSCsp2giHkTTKDimnPAMWTG0nxIslM8fJ
1ImFe/bDVNOM96b+RF5NSbhq38YguRBktpInNpSdsERkPtfRkoyaTCDAtthCW5aOTA+spRmF
cqkang6PinxG4nrHymWihGBN8iACL+nLTK7Pcl75Y4POKUghyjKbq9tlL4N5yQrSsIjsnJEy
C3JQdTqNMML8pSz/GM1vjkqM4dwVDNWjzy9+CXPwmZ1FT7kXfOYSjftJ6D9MqlBYrkiY2dPO
AinjrTAqjmof06IImedPMgJ0QS+UBPfGFE2BXTCdKV0nB3/4oswJ0aeGe2ZoyAfEji0vln9a
KiEjysboTfVxJgplQXGE0n6Etai4XEtDlXpEQzUzhLd8YpDSds2bVvJ4kmsiFxWUSK7+5lSA
E53SWM3xV7LSbFUyLOSszAfCBtaPrz6RDwW7GjLEbMWrkfXoE2d5k8SVMKCUJGl5Ats5wa6T
kUA7lUjl2dnqFfZcMqwrT08IxL3uSY/gKqU0a9nNxf4mPl8FEGh/J9oZTc6s1zwfP/GpWDby
EZOZYC76sPnRbdr0qnaFbfPAStTgJme4zHQqU2dpw7yWsbrGlOsduyhCvUrJk2Zk3YnSyo7r
GjS72h0Od2GGWl1+krNPAxjhygffeujRudB9KvzCmyfNApibgazTb73x4PrGBH/ryi9Sa9pL
2l5vqpjlqTC8BNfa1s2TKgKWrUp8UThF2CASBk9FOWbHd27+lb9pJG+H09fcbQrQsT8d5+GK
uUQdj/C0KqZvi/3HiZXF8584Xa8hiazaxU1kdrSFbEfX2lxSHu10VE3aK+db5CO3cnf69akv
J/tMVVVTspqtFS/lGq4VVhOjuuqtb40sZrCJ0ofThO8iC9veJa9vNy09o5X7KugUa9m96J3p
UMOcZ+0qMZQjHBeQ7+ldIUvNuFUOCmOnG0vKHlCVclkxNkwd6fxMb5pfzvFANxznb0Xozbvb
MVSTLEGkBhjW2IV0qoP7YngG23EC5bXZvNlHANZaxOJMqkQ7K9MNoRoY0/716KDh5Cln2rrH
pnBTyTky/xaUJmN82gA3aaFq08j+2oI1IRgP+BEsQ9uJzS5KlJ09UyPy8tz13nZJ8czuSL51
tX0FZ8F9o4vsLfu/2760tx8+HoAH/Irek6+rSR1t3zCjg711+LxVCuZGTrzdQ3G3dfENDIxj
S3toTDJpY4jyO/t65OpUOWFf7qaUH3yQCx1sz0PuV5or1M7esHmDOfhVcBvW2D+XtsSFDnX8
PiZGE94Gzqf49MVMLOKNAYW6i22LxvrbhYsK+9anKDxKfN0x3zt7cI7D9Ukdpc0DKAeBLZ4Q
tINO7XHPUN+LJ6W178JOV++3kOxO56RdSO+f47vb/f541QjeFJMvIOLTu/O8Z0hnGQOe5z8P
+tCLfvT+pC+96U+P+tSrfvWsb73rXw/72Mt+9rSvve1vj/ueoV1nUe+97wty98QL9ch9+73x
j5+jSmPes+osPvKfD/19BR/OkdZ99K+P/exrf/vc7773vw/+8It//OQvv/nPj/70q3/97Fe/
9dnyfrWnvfH8mT++0h7/C9mfZ/fKP0P8vwj+5zn/t38BKDzAU38DaID9V4DJ1IBdgzEAWGrB
A38UKH9794A4Y4H3BzrO90gbyDMdCIIOuHccSH/0ElgeaIAimEEjGC0seIIEWIIhGIMKOIEz
6BQq+IKfw3v1h4M76DUxSIJCiIFNYYMn5YJFGC1HiIQ/qIRZ44NE2HhB+BL+OqiEO2gjnFeB
KOgWWuiEzsKFF2h/TFgxYciBB5iBZdiAOUOBZChWXniGQeiGgAWHIcg1UOiDa1iHczgVd/h1
ftiFPSiFCEWF/Ed/aVgxgmiInsOHnVOIC8iDSUiIigiJjPhXgCh/vGeFkyiFmDiIN/iJVSSJ
zteI47CJxVeKdPiFKpiKLTGKbYiH+reKI9iKBvGKMFiLtjiLJZiLLfiFvniIB0GKoDWMQEiJ
cbSBDIiLl5iMwkiLKdiMxjiFzMiLzriMUaiMV9iLqIGIosiA3xg8EqiB4LiE4diNW4N/6UiO
CeiM67iC7FiO8QiBxKiO7WiPVViP8hiKv7eN+DjzMZO3GgAZkP4hkKHQj+2HkAmpkAvJkA3p
kA8JkREpkRNJkRVpkReJkRmpkRvJkR3pkR8JkiEpkiNJkiVpklckjuRwjifJkhOHgE24gNMC
ggfZkjUZXDRJC3xziyeIkzbpkwplifmIgt/oKpvIk29hfR7Ykz/JlHOEh30DlbB4jV34jtLo
KleJhU2plalmiSKoe10Zib+ok1z4khZDg1uJlmIGlo84gEkpiUcJgzF5lmlJl9q1lj3Ilkgp
jmMZlsFognUJmIKFgKiIhvfXljIJjuZYlUMZmI15fKxYkI4pmc0nlJNpmZeJmZmpmZvJmZ0p
foEAADs=
--------------050008020502040502000001
Content-Type: image/gif;
  name="ProtYield_energy_usryield.gif"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
  filename="ProtYield_energy_usryield.gif"
R0lGODlhgALgAcYAAP///wAAAKCgoP8AAADAAACA/8AA/wDu7sBAAO7uACAgwP/AIACAQKCA
/4BAAP+A/wDAYADAwABggMBggACAAED/gDBggIBgAEBAQECAAAAAgIBgEIBgYIBggAAAwAAA
/wBgAOOwwEDAgGCgwGDAAGDAoIAAAIAAgGAggGBgYCAgICBAQCBAgGCAIGCAYGCAgICAQCCA
IICAgKCgoKDQ4MAgIACAgMBgAIDA4MBgwMCAAMCAYP9AAP9AQIDA//+AYP+AgMCgAMDAwMD/
wP8AAP8A//+AoMDAoP9gYAD/AP+AAP+gAIDg4KDg4KD/IMAAAMAAwKAgIKAg/4AgAIAgIIBA
IIBAgIBgwIBg/4CAAMDAAP+AQP+gQP+gYP+gcP/AwP//AP//gP//wAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAACH5BAEAAAAALAAAAACAAuABAAf+gACCg4SFhoeIiYqLjI2O
j5CRkpOUlZaXmJmam5ydnp+goaKjpKWmp6ipqqusra6vsLGys7S1tre4ubq7vL2+v8DBwsPE
xcbHyMnKy8zNzs/Q0dLT1NXW19jZ2tvc3d7f4OHi4+Tl5ufo6err7O3u7/Dx8vP09fb3+Pn6
+/z9/v8AAwocSLCgwYMDAyhEyLChQ08KI0ZaiCmiREEUH2rcyJERxYyNQFYSiTFAx5MoT340
adFkyZUXAVh8SXJQzYwradJ0KTPmzJRAg45rWbJozp6Ech41RNQmT6UuF0I1GpWn0KtYr4nc
WrWoU6pemVr96lUq0alJm2Zdy/YZ10L+R0GirYlULFmzcLu+bcu377K9d/WyFBw27yGchJHO
HVzYr+PHuwA7lYs47liyky82/Wm58mXIoEOLHk26tOnTqFOrXs26tevXsGPLnk27tu3buHPr
3s27t+/fwIMLH068uPHjyJOXo0xXufPnkwNDn06dcdnP1LMbR9wTe8vv4MOLH0++vPnz6NOr
X8++vfv38OPLn0+/vv37+PPrB/+rcl2mbjmD3TEDFlOggQIG6At3/8GlIDMHEhOhMBNSmOCF
vTBIV4XBcNjfg3+BqIyHH/LynU6IkLgghhCyGGIzKmYoTYza1QgPjSvaqKMmPjU344875qaW
Iz5CgldjuAz+0AmOMgaJ21IhMZmWXawo6YiVjWCpiJSROfmkYD3C5NlOW16mYVViPhXmgFou
0qYibx4GpJe1bZYYlIulGFNYefbZ1SEDBCrooIQWaiihIc1J52yS1YXnnQOSdGZ3EqGVGV1x
IpIpoIlGw+WijjVq2ZSK/XmYmWqCReqRmCWyqSGvFhIrRoqCCpuof+KVp55UOmrdkcD+yuWs
E9Vqq2uNUpqWmqPyeimaYl6HorK2fJqLtcdm66mx2nYbDra3gOvtuARyS+652Ihbi7rottsk
NOy6Ky8u8cpS77z42mtuvvyWu22/AI+4b8AE5whvwQgPcy8sCyfsMCcNuxLxwxT+WzIxKxdX
rLGRA2/scSoZqxLyxyQD+G/JKIPcccosQ7xyyzBfMjIqM8escc2m4GzzwzqT0vPOCP880p4O
coIo0CULPQmUSX1CLNJBS2hq0Z08DTXBSkvC9KWbWH01vkN2+JPJVX9NctZDp+i02R+jbeTU
TXviNdtgG0h03GXTffPLksytt7tuL7n23xQH7rLchBfOdyR+J06u4Zuw27jj3kLO4+CUF2x5
JpJnHvXJRnuu+eKQTC66rZtXhPnp/KYu8+qs1w1614IcHfu5rluMydFamn57jbmPpEnvhPj+
e3bBU2It8YMYf/x0yS89fPHUP/846XBWX7v1189eCfP+23NfOfauag+A8+InF73W0xt6fvrd
rl9sJuC/D3+28nP8/aDmo39/cfl7xLDsF74C/m9RASTS7gjIQP8dUDgJ7JQllFQ/Bz4QOBH0
yAIreEEEkg9WDTRfB3eUwUUMkIMjdFIJt7RB2zEwhTZaYSJOuL36wRB4H5SV/SgowhsiL4fU
sxIPm+dDHTVMPBW5m002WMMeFhE6RwyPzOC2xEsMcYgGfCL0fCYni1GRVlbcYQi1+MNRfIZL
YxOLErMkRhuS8Thh688XZcLE97nxjeqTGpLAuD9EUbBQeHROxqKSxD3SMYxCHGMgBXkKs6Bx
jYecYBuduEgANpJWEaOhHSn+WcnhXIwlfBRFjACZyDt20pONZJUov9fETTaPf6ckziBDuUpK
uA+LroylLIFYylb6UpcQpBkoS4GjXubymMDEYCqtY0ZWupJ3L0ymbz55yEw684ogtKA0TUNN
QnLRlq18Ey632ZtPxvETxQxnNsmpzGVOLJ3PXCc7p8lLdepQm/MMTc2sWTpYYvOe+Swnxog5
CWy675cBfZLPpPhNxk0ymuNMqG2iyJ+Glu6h4oymRGlzxC42s28YledGb7NPgoLUmBnF50j5
IkO1ndSeQVwpSWlGy1CMMqQAPR8sZdoacxpScC+NZ06zyNPVULOKH3UoSkWq0aKi5qg1RWdB
d+j+wl8KcadOJY0502jTqSaSiK8UVP+y+lRhFilyXjVgplBI1tK01KOPONRXQQhWorYVNCW1
aFw1utax3lWrA9XrlfiqKb/+VTRRpNRZOZdWujpWkYfF60IVy0+lcqp4pORkZNsC1UjW0rKP
NWxTN8vZwCb1onHqa11HS9q1BPBEcN1ragu7WpW2NiGrANeGGqvDy6r1to+B6mKjNEPeUhKQ
qwVuaXOGtynOkKvZm+1lcWlb5e7jnBCpqHOL6yZ/Ena6kLWua4kJ20JylxHEky5d/4hV8Wal
s5w7UIXSK6v2hrS67r2u3Uw42IOy8Zb5ZcssacpGInptqfgNMD58SmD+9MZ0r/FMsILtAV+T
OtjAF43whLFSYcHSFqFuUqeENzyPrULycBcGcfY0TOKgmHi4ryuwisvH4han5K1kCzGGIWxQ
GwNlkNCVqoyRqeMe+xglL66sjmdMWyMfuSPdjCpEBhvWQFGZh/59skOi7NmuDtlvCNYylC8Z
5Cl/uW+ZHbGY1TEyGMc4xaxt8kPX/BAcU23JcQavMem8ZbNe8swOrTGfGfJi5sIps90Va5gH
3RDhGhrOT2OvopnM6IA42sJ4jrQ8I1pp3OZMu6CYr2Z7aOSqdpofSXaz9AAN59odKrmn7oed
m7tiu/q2tuFFr31jvRy+mTK0Ifx1rXmdjln+IzXTVB6rsGlM7GL7etSwpi+sh9xscxhbyo/V
dG+tOtVqo+PaXWZ2nokqbVfv+tbeJge455vlD+OaeQdO9znW3d9AK7t68ZZ3rz2M6ejaW6di
PeZcqa1vcCQ2t/XOcPngrfCCf6uhIRP13GJFXVu72+HfoGh5PxviiXcc3w3HuDcyCWrlQVLi
/+YUwyEs8ow32Ly7/W/KmTrwVrd8G6+NUIRMjWeVg5zlN+cGkFV9qpOXGdrmPnfFxz3toE8D
u+iEOpHm+BkXevzQAf8qzy3udGp0Vkpbw3ZGQw7sgVO860L/M7anbkg2AbvnZf/tsNGuFXcS
XY2REjfX0b1e7zb+/e90954oT4ziuJN920Ldu5oDPwt6a8rvQK+vopFr63Yz/mCC5/bhFy53
vV8eiO+i+cz9nWvEfx70JuI84PlO+mUv/vQS29fSI99xRHse9i6Cker3/nbVux73T5d96ef+
8eGvHvj+wrzPj2/64v8e+crvN78FPnrfj/r10H+5GSO+e6Yzf8fPz77uH632dVr+4s6/vvhz
72XtP5j3zW+9+tffok8fvfBD9T78AXx+/dM/9fZHeJdTWMvWewTYeej3fwiSedR3Jee2afOn
gMgAZOEWat2HX/+UdMiFfRJogfbXIO1XZQHnf6BVbsbXgSVCUN50WvmHgWK0YwiIggr+Q2Zr
h1aAYntrY3aQxYEyiH+Dd382uHzwx0o6yFY9aCHuRH4tKAoZ2IBDeITVInwFKElFiG8PCIX0
YlbvdIFM+IKah4ULSF5AOIBCSIL9hFWzB3BWBoYAiGlKtoShcFVZl2s8yIar5obSB4NmSHt6
OHZ22CXbV3IeWIa2xXNp+IR/2AoaJ3VBKHqrcIh7mIjlx35fmAqQWIeSqEB2l4eVOAqvFm2Z
eC1294bvF4nghGvxF4qxsFU8gl1uN4Whg4pIp4oiI0ytqEEH+H20U1t+RItRGIB3tyzN8Yqz
uItyhyWY6IsuJYYCaEJUx1WQ6Ilq+ImmqIyt2IyR4R2+d4X+eXOMumiN05ddYyhA2vh4I2iJ
mAV5yQiOtPaD2PhcbXeD56gKDDdX68iONRiCJjeMulCPGjiP+Oh+SjiIt+CPBHSP+JhknCgL
BomMAfkKhbaQsdCQ3/iQZhZ9tkCRiGiRQBWIjEiGBQly9siRpuWO46g6uaCR1UiSyuNnA0kL
KtkmW8eSKLlQ3JeRiJaB6kWTF2mS76g7KSmSBbSTPOmDUbdx+lgLMcl6RVmTIqKUQnmQTNmU
QImRUAlWOvmPVNmRSVmBPZmRUemQWLmVRsmVx9aVs7CUY0mWIBlqDMWCVylXVgZNK8mW+ZhE
gihkSaJylBc+CJmJjreX2SZndun+lJTIkLemWn8piYEZkoMJXoX5Zk8Jk+nIP4oZmZJ5mBOZ
mISpgZgpQZOZlpwpjzu1mEfYmGD5mOtlgJ95Z5oJCzt5manYmq45fjOJClVlSjVXlzxJgVtY
kVVyXORGm7gohiBIkKYpbrJZjJ7ZlFymZMnpbssJnLxph8+JadE5ms3ph5tHkgppUdmpmv/I
nXzIkkP3mxsZnKC4m40Tnk7XmO4pWjBFnfF5cwOmW+XYnI45nOzZbVTJYMGoRstIZFfJn8MJ
Ulv5dfEVc5S2met5oLeZnr6ooM61IRbBaZT5oBA1i/VZcB8ZOR9aJvGoh/u5oX0Jix0qcueJ
l2tkFRj+KprTWIqSB5AJaJ4TmJ8EmqHdF6M5ypoWeZ7YwqA9ipikV5mTRnA0+Z2n9aJEuqMi
lmzOOYkE1X+w+YCc5mRwl6RSCp7VGZyztWdZaqPACJ1dipt9+aRYx43sqKQcx6S9cKVzFqbe
qShuygtwCqbEp6Wp5JU9Wae7cKcNSp20yIoBqj/TuIbCkJtxmqdiqoJ36ZR+2o/Wp23/SYN8
WniRGpS5GKkpinFsGoKZKpjLx6lkuaLYWab0yHlU6pdqOqh0iqrouKlMB4uuOqan2qm3V0Nn
Kp29OYq3egybgqXDhqvVlldLCqu4aaTCqpxV0qreZqxtiqynMHaltmu0Cjv+QfepyCmtpiBO
acZJu5qsdBeR3xSh/YiGWhmu85mq42osxHplD/aiTciuaAetSfmuBSZtVwp5jxh49kqQEvqn
v4WnvcitNmefNomUyOmsvlCFBCqWyYivCbWIJ5mZ0OCwcOqXriCxAdVROYaWzYCx3+WFrsKw
UPqeJQmyzCCy6rWqsMqx+YSaFzuwuqqOQ6hSMDtPpPiSy5BlpaY9u1mjYdSu33ST02B2BQs+
SDuPOOuvJmW00qB1cqlIUkujTUu0HwW10RC08wmxMhqwcuqpzFixVRm1FndFcomoD2s06uqh
Y/uTesKP1MC1LOa1BkaX3dh1OTeiR1t5OYmVScv+nOWpojc6oFt7tjh1XzE4PVhLIUYHt8BA
tz+rsUslqGGbbiGqC0K6tYaYuIvWOyZ7uW5LIJcRKAGgts8wK5ImV2hqgvTjtApDddWgup7b
un5krgfbcgl0dODGo7Jqj9D0uSVLo7Mptgx4DcE7lGGZuEJbvA4HpBKpDZU7Rj9LtzcLu28b
vdkwvRQ5vc0ruM3WYXDpDdy7vADGqD46utr7ld9QvmMpvOjrvOqbsIUqQOSwunPpT9+aeN9r
uaemrXoZDt5rvuuavmD7v6gXsuuatpZJslMpgvWqhTzLDRsoZ/B7e7j7vL4ajtL7c6uJv22L
uKdHoeO7DXR5mfBrdbn+qrfldy+CcrrE2w1cK7w6OJs5K1EAPGU3DAoz3LokS57gK2+mykU7
7DSmRsM4KMKfR29FjDgrVr1qxa/yq7sIN6XncHZ1+8N2tKtNrLNqV78K1MWxSGNQHGzgCnsk
vEpibIwW7MN9uMIoS7+QKz1rzLhkrLhm/Hd1LE0kl7mvs8cTZKU+7LBKfHniq8bzgFKB+8Y8
CsjAxMTx8In6WorpFcKJhrpiNsTN5MjTysAmaISSVGk5bDQwHMPsQLt9qMKvK8otbMX3gMrc
trSY3J1PlsZGDMP4QLtpq8XOw8nWM8ou48ueeIWT3ISh2r82psmITBDFzMu2FLoKBskCkbz+
m0TIxsVnA3ZJwqyU++u6oMVoDKbNBfG3aNvARNhph9xV20yk5byGwmquxxzNW+pF/LjODpqV
dZvB8RxgwExcx0kr9gyb4TS1S0csLkti/ewRshvQr+C++OvBxQfOhXsY0JzIMAWmK+dgDE05
fuxF4nheBuHQO3i9GczPzAgx8sUQIo1CmrbP4sVPQapzKn3RlWi9MbXRmQPTIFpmvZuaIIzR
EE2aFf3SbwvGoDnOeAvUa6nRBltUfdzRLYkQGd1LnjzLNy3Rx/sPU13TJM1aJX1XMgsQ5/vQ
LC3UkyaTUuxeYT3NWgmHGvu7I9vUfJzAkYxu5/vWo2qtQUxWax3+ELGJrksNgXaFokOtS30t
1kxZ1pvqs8yJ0zxD1/BAlLLowPEqgjV8sk4N2e8g2QgouV4oh5OHu47tMIet1YlthYI8kvyb
b3yt2e7A2eHVueS8uPG7UaXtD7AtbNQclqwruhPr2u2A1gBZgN5cSvt7YV8dS7fdD6ydZ8U9
yFA62gCz3PzQ3LpMjTR8yUcqU9SdD/oMbUl8X9/tv4vU3Q5xR7PdW1hKqXL9QObdEKAtXesd
uCiK2geNR+890/m6qIeKZRwqWtLdPVaJEvCMpqN6gqpM24YN3A+xaLCSZm2bhtdKRvl9FQ4O
gcQ90oBrylpU4ULBfwv3arrc3xq+4Fn+3RcX7tbOza8nmkweHhQgPFoSbnMxPuEw9OJAkeKU
7ERcvMgH7N4MLtWkdH4z3pngR+IBfiw4nhJ4yoV6bKDPlNzAWtiPSOUHseQEHqg1bXUE7cbY
wLFJnoKh+RhN/uDnGN8nfNVZbOZWbYlW3t5OzFJBDt9aPtmMTF8IxsUFOrSJ6hcR5CM9Pbt1
TtsPbcYV7Nk6yufBEJ1SPufGCdKiUeZCPdjMe7s+zuYI6T9g7oAcbsdDISH/fJahQY2cjrqT
q7zvi7iLGd+6NofD0NzSuG969FxznOVxWOmoOK9f2wrXauOlfrWYfevfQrbVUiCBPrdd6L1i
2dbB6n+NPqz+DeeCPEzLefvsyfDnxk4a1m6Ou9ztJG7ZHI639BPe0WXJp8i2oj3pjLvt124M
gJ4cBpXUfkWts7rXyAZ00o6t6KtN3pwNOcfTyHHXQDvvHzyjoKyf0Z7wmB7K+r7C/L7r6eLo
qLGsUsmLvZ2x8mnvvX7Gip63rcbugxngWB4bFL/spYdlLY54zX7u+ldBvGPu9804EfqluyPi
cJ6FJ94bhb5trE5YZc6dK++fDaTXoJjH3G7ToWzQZei7+02i2jDysBGNJh9/1WpfvxZp3uVf
3lqw9j2XkDY8fvcqjA3bT/ztaZfzvCH1T36obB/YCj70VrudBuhGOt72rNrpzDr+bjS/8PnK
9Wc/4MCR4J4Z9Fvs98A7tc0n3G1OnoI/nj162dWY0cQ39hmG8iAP6miv8xFopJRL8CqvrG8n
r6FPVZLM34FN92qK5kgKqDxm+hGf+Wm/+Wk6gqBrO6hfwUyV+HyZ9XgMmeL+3wZc+1bf8gb+
92Me+L1t8Mit4Tu/oU4YRO0FxKudxSkv/Pnr6ger9Nm02/st6XUH+6CiyK4OsUfciWZfVxf/
w7dE6o6P16435Le59zGu/Ov/5sl3/N7i/aqt60oNUJItsoAAIDgAQEhYiIh4eCg46Nj4yMiY
CFnZOHlpObC5+Vhp6KlomfnJ2Xlpijm6ytrq+gobCxD+IFtre1tLi7vL2+v7CxwsLJzaqVos
6biImkrpjKzqHFlY/KwInXg8aLpdjewaHf3cLEr5/dlt7B063O4++x6Pqytfb3+Pn++7LJ5N
Ckqt2jVuATmVY3YqErSEoibxY+fPYTqJsASKo7iuoDFWEilO0wfyFb2Q+UaSPIkypcpYiwii
00gOoMyB3FoK5KhM2UKXD81lnAaQHbmPEBXWnHgU09CD1qStBGnyabuoUqtaverupiZPDmNG
zPmv3NKwy0jlxOjvZdCvWs2GdQspHEywhrx2jIsUqz6qennx7Qs4sGC8FZvVNazzm0WX2Qgq
JnyW8GGvDXcalIuQ4cVRazP+9ZS0ziPkwe7+koZl+rTq1ST7vRyN0THdtE2NXg51F23PvGQb
L8TpdidYdBbx3oVblLVf5bZSM38OPZjrsdaEA018O7gmbNUdY+sYOiNlyzy9cYc50zP6skhd
R2/l/H389/Trt6JuPK7st5K/ZsbNmVPVbbWVVl0dNWBtAuaX4H9zebbUWO7Zt8p80FlIYYYa
gsPYa7Rl9pt+5A03EYkR/qTgZLc9Fpxs7DHlIXkv3rehLBgyd2ONOm44YVPaiLWgdg1y1laA
u6X4YXoAxphkkJhFFqSHO1Y4JTxVXonlOAZpRFRvMBL4EX69zCTacSrWhCJ23zkZoIM9Rlll
jqz+yZllnaf9SKSBebmHloP7eLelTdbN9g9jZorIEJew8QeOnZXQqRqkjk5a1TlwrrLYlsAV
mdV3HbZJl5i/lSUmcO1pauqkkpK2KqWuroQnL7GWguql7bxJ6ylkJkdWqYUFihxHn2LZqmDF
vopsSKD52qiUHDJLzC8yJtqfaH9ml2wsxwK2bbbeyjOrLTENSyNIuD6rK0unUltLhN+iFue7
8uoVbkrnAnPvr+Sul+8r9c7bSLd6CQxwwbj8e1K/u0B7sLjgIoorwxkSfBXFBl+MrqAK35rh
vRsLK3GqdVpcFckYn8zcx1h5/DC+lJr8FMwoyUwSzTVHt5TNw4g5Usj+7V47rM5QRRdA0X8J
vVeNSJektI42+4zpbAPQAspnyi0dTNGCaD0K1vd4XQ/YYTdNtnJWT90QoauJ3QvXXFvCdjxx
l+Z02RrO3S5iilAN5IysXkiL0RUaTXjhhh+OeOKKL854444/Dnnkkk9OeeWWX4555ppvznnn
nn8OOuiAB443yqafbvCuaQOJui9ul9567LJTqnrVas/enC5v4857775HfbZ/KsdO+O/GH498
8sovz3zzzj8PffTST0999dZfj3322m/Pfffefw++9sVfLTgr428Ne8DlD/72+RWv3/X67pcM
P9z1z5I+/rvb//r8Md8PifP5TyoDVF//ADj+mNfNSXd8UaD+BsbAwaEPfxN8XwXtN0EH0u+C
j4qgAbGiwQ5mMH8V3J/6RshBAnoQgxQMIAldEUJWrVCEFCQdBFN4whrOUIU4LKFJTKiSGOaw
fQhMiRB9qD8gBnGHQ7Qht5iIxByS74FrI10DrShFq2gNgVvcnxKNiMX4UbGAYKQiBsvXRRCG
8YxpbOEGuYjFwt0QjmZ04wKdGKk28s+NXwQjHjuIxyL60Y4u/CMhVyI4JSbSiX08ySLZR8hG
huSR8bPhEZdoyA8SEXA9DMwlNSjJmkGxhZf8XydjGMq9jFKBhlPjKXeYypKsMoKlnNksLxhL
Vx7SWCvcJB/zB0r+erDyhaQUpjFRuMsy/rKJNNRiL4+pw07aEpfQbGMu8RFMZiYzj9dUIRG/
yUAy2hKcBgynIMdpTHKir5tfk587TchOewjwnekEJj3LiU8IqjOJWfRePD1JTL8EdB4DJeiG
/gnQ8Cl0oQxtqEMfCtGISnSiFK2oRS+K0YxqdKMc7ahHPwrSkIp0pCQtqUlPitKUqnSlLG2p
686JGpi6dKYnbaUjf1jEECKUpjzN6E6BsbtPzvCnPS3qRL04zwOmEx6lDGpUxqfTghp1qr2T
YzGvOkxkftGpNHTb1ppJ1bCCFKkjLB5Zo7lNri4Sql91oVjf+tGz6hGHUB2lHbMKSJxDwnWv
G5UrIw0JTzJyFa2+dCtfD2vRVjqwrvx76tHoaU7DrlOqiK2s8lBJWctq9qJJzexmPwva0Ip2
tKQtrWlPS59AAAA7
--------------050008020502040502000001
Content-Type: image/gif;
  name="ProtYield_energy_16.gif"
Content-Transfer-Encoding: base64
Content-Disposition: attachment;
  filename="ProtYield_energy_16.gif"
R0lGODlhgALgAcYAAP///wAAAKCgoP8AAADAAACA/8AA/wDu7sBAAO7uACAgwP/AIACAQKCA
/4BAAP+A/wDAYADAwABggMBggACAAED/gDBggIBgAEBAQECAAAAAgIBgEIBgYIBggAAAwAAA
/wBgAOOwwEDAgGCgwGDAAGDAoIAAAIAAgGAggGBgYCAgICBAQCBAgGCAIGCAYGCAgICAQCCA
IICAgKCgoKDQ4MAgIACAgMBgAIDA4MBgwMCAAMCAYP9AAP9AQIDA//+AYP+AgMCgAMDAwMD/
wP8AAP8A//+AoMDAoP9gYAD/AP+AAP+gAIDg4KDg4KD/IMAAAMAAwKAgIKAg/4AgAIAgIIBA
IIBAgIBgwIBg/4CAAMDAAP+AQP+gQP+gYP+gcP/AwP//AP//gP//wAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
AAAAAAAAAAAAAAAAAAAAAAAAACH5BAEAAAAALAAAAACAAuABAAf+gACCg4SFhoeIiYqLjI2O
j5CRkpOUlZaXmJmam5ydnp+goaKjpKWmp6ipqqusra6vsLGys7S1tre4ubq7vL2+v8DBwsPE
xcbHyMnKy8zNzs/Q0dLT1NXW19jZ2tvc3d7f4OHi4+Tl5ufo6err7O3u7/Dx8vP09fb3+Pn6
+/z9/v8AAwocSLCgwYMHAyhEyLChQ08KI0ZaiCmiREEUH2rcyJERxYyNQFYSiTFAx5MoT340
adFkyZUXAVh8SXJQzYwradJ0KTPmzJRAg45rWbJozp6Ech41RNQmT6UuF0I1GpWn0KtYr4nc
WrWoU6pemVr96lUq0alJm2Zdy/YZ10L+R0GirYlULFmzcLu+bcu377K9d/WyFBw27yGchJHO
HVzYr+PHuwA7lYs47liyky82/Wm58mXIoEOLHk26tOnTqFOrXs26tevXsGPLnk27tu3buHPr
3s27t+/fwIMLH068uPHjyJOHo+xZufPnaZ8mpgu9OnDGjZda3z7cs+DPPVuKH0++vPnz6NOr
X8++vfv38OPLn0+/vv37+PPr389//y/EVGHHlDPgFVOggc0cSIyCwjDY4H/SgVXXgAlWaOFf
BF6IoS8AKtbVYRky42AwI/4X4oYi9iKeTiyCqGEyJZqY4ovIxMhhNjZyp+M6OR7T445AhjfZ
Nj8iGKRsajn+Qt1EAi45ywCdFLngkUh+qKSUadnFCpSOcNmIl4pg+SCVsEHlk2Y/hVdVTIdd
1iFMazbXoiFgLlKnIne6qBWZZXI2nZV14TVhXoAtFih2gtI1wKKMNuroo5A6GhKOfL4m2aGF
GYoISW8ONhNamSn6SJ6IkAoXpZW2dqllWXrYWKuwglqWgK/SOeqtHqGa6mqrfoiXpm1qiamh
wFZiKiVikrgra5eqGV1msybCHJrU+hontLYkC4y2y3arDLc3eivuN+CqOO653JTLi7rotqsL
u++6Ky818OZS77z4vnLvLfvm628q/dYS8L8EjzLwLAcXrHCUui7ssLnYJPzwxMj+NkzxxbRI
DIvGGHccpsUeh7wKx66QLPLJJrOS8sker6yKyyxfDDMqM8f8cM2m4GyzwjqT0vPO//4sitBA
40s0KEcX7W7SEHUiqdIsM82wJ8dCTbHUnLBbtdU3g5zJ1lwvjPUmWoeNsteYgG02wWNrUvba
LaN9idpw59t2Jm/XLbPcltCtt7xj33Ttpp/4/Xe7WLP5EqZ6cmL44edKTauz2iVVOOQTB27V
r55ummbamBecZGQReji5TZeHLjaEoSIKXt6qs826XZWjTnXsPKs4OpzSpo570Hwb+7vsEaf+
9PCRB/9lo4KA+TjyO95dkSbOE/I89NxJf4m21Q9yPfb+1mlvCffWlw/+uOKPRL35zZ+PvvJd
sg/A9+47l37FX8tPf/3J3T8J+ZCaH/+65T9JkM978hsgnwo4EWMdT4AIVGCqGAgJLHGpe+2T
YKUo+AgLZvCDENQgmTiopLmBEIMiPBIJJ9W3EyYwhdGDH55cGEEYBmmFuWphCHe4Pxv+BoeL
8OD8HthDH/YGiB/TIZQu+EIjVgeJvVMiBFHoxPDJMBFMZGINq2hFg5FHZSYUoBZByEXo9Ks8
YJTiGHdYxuecsXE0C+MSadhGM/oMjqcQ4hx5WMcnXrFU7VtjEftIGygSrm+PEmMiCYmchEWl
FXqcYhMZSZyDmYVkkRQkJZX+Y0mMyCSNwhMjHTdpnE4+cmSTeNoeqUjK7uRRUC8TXhZPyLxW
CseRngRlJAKoyUHaUjWGxOOtZslHX/4SNY5kCSRlGUgXGvOYprEkLAHGTEXWco/Q/GEeP8nN
WFJilWxsZja1mTNuYrKa2CzfM8cpGml+Lo7flGSe0slO3STTSUNDZzhFWU/eBJNCqZTnIejZ
z9uck0nRkQ51YpRIYhaCoAWtzUE72CnGARRXcwwgAtcZUb+c8YsdnFCxLLdLcUIUoh1Fks9A
WsHPcO5178SiSfeJ0pQyK6ZNu+iVCLU4uuSImCfdp03LpMudYsaipyopP4PK0aGuRXNHRSpJ
IQH+VJoK1amqgqeNVrQiZx2SqjOd5xRriVVevbJWpfhpWAfKx7JmtZy28yZYl2pVVrr1NJ2M
KzWVusQHbtR8fr0raNypuG3yta3eIytiBTuae5pMrfwko/4myVjH/DOpVNWoqVDY1MompKiG
neudNjtZz7YTtHAVra0AucXOmpYgH1UTPpEmiTVKtrWUfS1bYttVuWJ0tdZraG51m5W85tK3
8QsnaXFL3MdMdK+qfahMmdtcj/4xg6Od7iir2xaOPba2kl3udrn71NDO9hOQzW6piEjevhi3
m9D9rXQBqcn2djet40HtQ5ln26uCc4v2xUrA8qtf4PaXtBcMbICD8t7+ZX4pgs5T7F/puuDi
yjBSYOVlhS0cWgczwkv9nWFVN3wVaW4suZIc5kxJLJQGF3i+KY7fiFkMFBej8sHYzfCKaZwS
wp43px+GsI4pzGOU+Pi5M8yxiolc5I5cdqp2EvKS/9tkJ9MMp2lFcYhlOuMqb+TIL2bflum7
Yy9rxJTHjW+ShyhhLgdykWZuCJrhC88g3zbIXY6znLeJ5Tvi+KpYFC6g9TwQjf0YyFG+c5Tz
TOjPFu/PzwPnmBstkHuWTMsBZTSlYfvKwnZ4oBrF0zUFOuhN+8O434V0clUZY1MXBNUeTnSp
w1tmVxe606PLsqr/LGQtKtjW+Djyod2GaV7+Nw/DAAb2Pp6s1zWrTb1KVjSoGaVsdTA7zc5G
MXPtKutqowOq2PKqMNkq7Rd2j9tr9vY5ElfR2mFbu88Gbqsf51p1+8hc7QYUlLNt7LaiW7v2
Lge4TUfwNnm6hvGG8bzbzNqAL6fP2VIowX3aJTgDnIbV81u9HT6MgY+02W7GVWkzPleOD0XX
2ytdol4n45KzeVHRjmy/Te6N2I5vM2daaMtFzlqS85zm5PIznW+MZ5eT2+faBnrNV0pg5IY8
6atF+syVro0BN13NT5+6lMUJdapXvc6xFrHRo67OsXv9ujX6MMMvTvYJd/3se8J12Lsd6F/3
0uxwtwaYiR5caif+O93T9vsq1z7rvLulw0g2N97/ns6qbdzw/EJ8gaGt9ds2nu6Qr8benZ5s
jSda0OV+fOYzduVhZ43tnk+3Qwct+tHL4tpDT6zgF69KOJMq1K6PBuwLdPmff36x5M49veTW
+7f3HPjyFr7uiY98wP/+34VXvjF2n/XUyxj0DZc+NKhP5r+zPcnFz772T6T36i/+6M0f7vin
z3zoK9zO4Q/++mmUT32Z3/fgT7/353/vnLmM9/onf/AXgNHHfzKCcpdGZrj3fdMVf8lngGmX
WgnIVu6nfoBFgK0HgQjoRRBXfxRogfvXd/z1cn73gBpoJPjVgaEAgBUYgmQHfRl4gh7+WH7d
V26yB3M45oDvJ4NTclap9oEuqEYiuEgxyIPo1WmDAnbHdjydtXrhV4RGODVZdkqcx3UFeFhS
14JRiAtzlnhLOIKg4ISLBYVbSGz+p4K0VYNXqFpZCIJlmC2l94O2ooXgJXNP+IYd135uWIeX
d26Eh4eRJ3eTR4BC6HZrCIikd4YHJ3RzuIe7xF5U9GuIiDCSx3fSRYd8NWp2RYaTWEGVWIV2
GAoI1kSc2IklxIG9hXWGyFEKpoOlaIo5NDQspYpWCIWuCIv2oodBeDmbiItcGIdeWIumgGyW
54uBmIKLaDD3N4yltYPGeGJnSDY5VyAs6IhhNF6H+IwzuFL+bvMnX9WIu7gJnMVq2kiJimh6
huEqSQhyq4gKnIWN5RhmWZNrBoQmQmJwIqGDokiCuPeK5UiPuhMg6/huEyaJTjNy1hiPaYiM
6NgqK/eN+1WC7jiE9aWQ9neOWPJxBGmQo3Bu82aR0AiKFCUn4xYLHplFFgeSSvhpLPkkF1iR
KmmJsoiGR4gLJwmPMbmBM9mQeJMLN0mIObmQNHgLP4mJQWmGqEiTiFYLRZmQR7k9DJmM20gL
TRmOT0k2pTeBtlCVtCSRV4mVKRiSTyJoVKZcXymUK3h1K+mSvTZe/uiLVpeKaykLXGmWZ7mU
w2eTL+mWd4mX00A0ddlXf9iX/8P+iOyok2yZYxU5aYRpQIZJkIhJl3uJWIzZmJ6YlEophVv5
hZoFhm+Jix+lli0pmXOYkpVpmacYd0Qpf7dng6jJQqq5mfLWmtn4mvuWl6s5m8dnm+ODdqcQ
YZJCm58Ji9yXmwonnEPIm7A5lLJ5nOvFhMqZmszJlKy5m07ZmMXZnM04X5QXnVEUdxxZCn7F
bQfmnbGIX7EngcMJcMjpjOZ5m2l5OurJC+RJa5nllU/Zhbq0ntnXngnJn5Cnn3wHoLp5g+VZ
h325efFFoM5Jkd1pfEHpWLHGoNtZa4dIoWeXnRjKlweqVAmqVRM6mK3ASmXpmiaYn2c1kAhY
UyapeB3+en45+V7PxaKwQKKthl0iuqFUJ6P7WZvi6aL2aaI6qnQKCnY0+go2anEtOKRAJ6E9
iqEJJpE1lZIMeJQa6qOk4HhsxGoJ96G4NqNY2pGLRifY53xXWaQddqSuoKVitXpmiqISuEwL
SJeEx1RW+KZWOpojMY3fKXPCYKd+WqUR+qUZ6Y2No6ZbSaWSZidUOqgYmXLpSHGG+Kd1daeC
GqML4if3KBYRgah62aaWKn4RmpkVoUyQOnEsN6nBAKiIyqQ0504qSlESIqn86KrJWXaBKqpn
ymeQKavqSKu5+gujGKwnmqflRIV7SpIX5anGCY6tip2E+qTGMKzPSphOOqD+Ybqmbjanlyii
2niltlqsi7pru/qo2Bqu7mlN3gpiZ4mm84mulDWudpateBiMjMis2kmWhceuqxCeDmevyqiq
wKBe49qKVuk7hueuK0qvW0Kmwal//rqPmedjDhaxiSqlnCmY+PlXpQivEZWdJiqsznSjZdex
owewU+mxDViylgeGDTux8oiWDEuaJCdImqitMIuZUrmCt5oMd/ei+HqQObuTAPmYzvCzQUqs
Pzq0PKtTeroMSGuXoRqROCiOJxuzNSkNUQtt3HqdeEqkvrmqIxulOCh1XxtKARqZMgm1SkpX
79iljnO1fvZ/1GCzs/e2VGuUldekuka30+BrkXL+h2RbtQeLf143YL36tMqwZWLYjlIbt0wb
n0UbsFp7ZyjJX7aXtGnTqHzLdJM7lUdruWvluOQotGn7aJUrVBqWcVzqtZf6qmE7sKLrtrgq
rzNbuOoGss/AuKM7uvp4jQkbu8IaWFXlc7bLuQgavNawKAFAuKErYoELcwXrb5B7urEJDVpK
uwUpT7W3sR6qvNfbDJJYvBRJarabt/MKvtOZDfqKt9pLo3Cbd46FsthwvLVrvqwHoZ0blVhb
Dfa7veR7trjrbTYmkuxLaiMrShomwCEbcAVMi9vwv76LwLoqgHCnsEbLDRpLbcHZvlMbjirb
TxhMueBwvvcbwK/bwPb+prvdYMLR65n5y6jeu7/IOHct/LDdN73Zi765N8Kgq8G4CoQb3LWK
FsIpxaNrqw2tm8P4O7jxKn1I/JvM663XMKVNLIzpasQd5cMtNSenosW3Q61X3HkpvKNxKI2u
Qzhg7DQMp8O1+rjpKr/TJxHu9klrXL1MjMIkWMRQnMQI9VKeMwA8ubtj+r7F9MZxTKSk2pto
5auXpKImccdWK3aG3IYVrL4kXKifom+oI8nrQ8lu7LgWHLlII5p7KlWdbA5sOsGWPMrWS7SD
rCV9FsnkML5j3IdPrHwPHEee/KOXe7dBPFazp8uU0suigGE/icVvZrEO67x6dq2xZMzu+ML+
rQyUe1tlXFw4zUvFqtxtjcde+eNqUUxN0vyb3ky9wOnMMIrNnyjF23wPbBq9wthD5ew+2XyE
9dyRdXrC8ww64nzGYJTPNcrPoRhpzMxi1CfQI3qz2vvBA/jPiptlCk2nOIzL8XS7umVi9jfR
Cw3Dvyy9ocTRw7PL2yTSI2pSD+umHHmaRUbS5WTSL1uWGYXS66rCJHbPTQPT/RpWyEy9anfQ
C4bTDKPTqtBlGQXOYMPSPLZ73EwPRh2KW5d/EM2NYjkQT32Hs+bEjRaXflsQVz2Grummz9y3
VS0QX/2OWe3QTYa46TmX//DRII3Vx7LAlMbVi4yUBEGELLttOYj+0Z5l1zsrswHhkRxbq8SI
fkQNOSwMEITdTGid1B8M1Iy12PywwE782M3M0Jxp038tvPEwuNmF2YhtqXpbWZTtD2LloD79
gsTrum512v2Q2pvN16NNxiDc1OwE25XtyqL9gohszZw9Trq9D7LdoDE3mRkL3H7dR8OtD8W9
nbwL1R3sr4m9Oqib17wNWHaHxTN9sNWdO9dt1dkNYK29xwJFh98tOp4ND89NXZl93Octzwy8
xev9DnVS3qF3gdxp2TmI27bU3PlgfTEM37OkaQVaVgCODwIez8RY4Baq2uk9LwleDxELg6Xb
lqDNaxEOOPX91qTI0ATLZJA9wP/d4aj+PduXWGbCBZPsmbH+zUUTfhDxxmjIVoFoHdwwbuLE
naP4q4A4Cc44WeLhDRnTG3gc/Nsp7t4bPkEGwqcQSeQ9/oEsDo6jpNVC3iCG6rSOwd/01dNi
R9fKXUXio5Ft3RdFDoRBbt7pTLhLfkM1whgH97lYMcRWfOOiSudhDkNy/i4wAcmkceZUjnyK
usQ4rudN7pCpSuQNxa1Y3XZ7CeZXTiJjUcewF7os3ejc6dvqqs7Y++Jb4umOhuV+onOjUaLL
aN4uDuJM5r/D0Oa5OOSPYervDWGB+5HkG2KSzcYzbLi74Oq/COtbTrKq5+h3R7tzup6DBK++
fozrG+vC3uL++z3MNH3kY4bsnyy2VAnq5Kd5f/7ssx7iMj3dF87DOOvP2P7TnI7Hy6HjDnHY
y4OxClzlez2sA43qdVe2rb7O1y4OMR4UTejY8l7Y/rXcybvghaPtZVzo+g7sy1fqEtvddj7u
wYzuz4Terp3wyZu+x9yz7L4uopHrz/nCho25Suru721MkL5oyIu24rjyjIfwtL7ruNnsv1Hg
F45BB1qtfm3xJH7NLK/xF49b1d3vs2HyYE3gtShh/EruLmfwsx5QYUh7eBy0MxK+wOHgBx5C
Ih+osi1eBd/0GO/K4Zx0IF+dQ9/xqIH1xq31Up65Bhb2FXrOBgrajufyqbTSaO7+QCu+dAzf
G1Y+ci5L0KS4f98zjhtro4cc8krNhzP3zTIf6FSPIjTvG2NU3LSZ5IN3+IRfcXfr9hAedTgM
7Zs7zHPdtg/6YHieLmh/GpVfnUcn7mxuVXzM8Tgq7Ztv2yZs73//vacJqnk746FPJKv/59B5
q5e/6SC9ivI9+AaKekw/WnuvhqVN25gXaouf+WWPgpN/9UEI/aQP1l3P9GTs9cVI98BMrJje
/OBH3XkP36CstDPP7cix9Cls/XOP85k9akluwW1a0WMsxIAwIDgAAEBYiIh4mMhYuJj42KgI
OTjIeBgpCWnoyJmpCRoqOkpaKhpgmqq6ypqI2gobKzv+S1trS1qZq6kbisk5OYnJ+3v5+5hJ
mHvc+Jnc+ezs2Tk8TVnpeE3cK/nJzG2Y7Wsa7dx9e37+ir5eq87+Dh8vP0/rK2wpDUwd/NyP
vU9MGT5ywZSBE3jJILRNwHYZRAiqmTSI1gT9uwaQnsZTGztqcucxpMiRJEXZi6jtpERx/qpp
a2lsE8uT+cp9Q+Zv2SiJCQfeTAkzqMKD+EpuBGlUHtKkTJs6nYUQ2VCVFS3pJGrVWzFPD/Vh
xfhymq6oX82FPZuSl01FQ6GNhWj26a2lcufWvYs3L8p8P70mE0juLdi0RY0BtrhQ59piMy8i
DspM7S6GXRV2U8wwpt51dDf+w+rsObTojfceyyxrkSZBwGIz+m19+DFVtmNfA6VtmiXa0qkd
ZyWcO6fX0e2IzwJtPLlyqK75Au9dU3NMcQBVZ6aZ+CzB7F+1Us68lXHtmuFuL2u7fBXy9B/Z
u3/PKi5ue5V9Ww1H7e/Arg2xG2ZN3kP8gdfWZX19hx1Zt4EHHynrNVjIgxBOuJx8MF3FVXmz
YbgIb9v149+GVXXokn4AftMQaohFQ9l4Li2UImQUMiJhgzXOiKNn5qAH1GGnTSRZf5HBpQ9+
lskUFX9KxpZNWSrO19t5TwbEY44A3PgellZu2VRzugX54o+LOQbigdddqFWBUIYZIoPefPkc
b0j+NklmjFxGeKcree4JX5XCScUkoLGFiZV0hULJYkBzgqVgThr+uZeKWVmI1pZasncpn5p2
WYqUb1YjH4k9RkTRZNEJaRt3woV3KqFueproNptmuhytm946UiRqPfrpRDsN11x8JZ75WqAu
FjRorK+OWCmre9qaHLS4TqsRT9ygZyJ0shoabKeC+WQsgi0WBmtVPQWXorKmaiotce1SCy86
pVJ6rpwmdRcPgN1i5qeP1NH5K7LdNpvju6IZHG/CtljbCk6kkktafOX5SnBgAKcCJsHzscux
wh4ndYyfAbvZi8i30PtwoCjBWc/EHzeCsGcxv0zzOChKvF9huERcy7H+2zAZS34109jx0EbL
wzBJKEMlr8gnnuyd0TPrNfXRVtupMT1LwzKwsDa/E7LOKW9da9FXn91yUV3n2yfG65S6Ctnp
VY0X3Wjf7ZTccqGs985M82l3XYGHNHhHhRuuHDWHL3wxnnXKG3SVi2s0uVIBXA5K5fNoHg/n
nU/oOTyHr73Xaqis9aHMNr6CuSShv/M6ZxTGjg7ttRtn3QCnR6fbZra3wzpdv9sF4fC2GE8L
8smPpqQhu68mo+CrX9k6zJdfj3322m/Pfffefw9++OKPT3755p+Pfvrqr89+++6/D3/88s9P
//qro1I93vrvz3+n04E6rP7JonX5E6ABD6j+P+hRLHUITAUBldfACEqQU7n7X9Ym6IrrYXCD
HOygBz8IwhCKcIQkLKEJT4jCFKpwhSxsoQtfCMMYynCGNJwg9tylwczlMELKu6EOq+dDpwTR
dTscYlKMSKMdIqKARkFiBh/oxJFEkYdAVCJTpkg9KFpROQ88WPA+8sUsPqWLROQh9czYFDJa
z4xqbGIY13jGJTKRJG3MIBshWEc53hGNV3xjEvkIyFr5kWqD1OMZ8Ye8PKIRc4qUYiEXCZI5
iqSRasSiRyj5xewJ8ZFxZCQEOyk8P0rSOJ78ZBm3eEd1jHKSiAylGO2YxlaC8ZWWNJwsT9nF
VV7ylmvUYCMJx8v+P/qylrZ8ZS/zp0svIlKZcbTeMgPpxmc6s5mopKMnwShNaDoym088JCdZ
yU05cjOZlLsmEZf5S2CGk4rUNKV6vikXTPKRnEfhJBnpWcxm/lGf/LSmNkGZxWqy8p91xKfl
CJpJeO4SofN0pyrSOcY3VrGhEaUoJPvZR4sCFKL11GhBjVdJVSb0n8D06Egx6kiTXrR4Ar1L
EF+aQ+1tsoo0DR4x61lTcYrUoJsrok+ZyFPL5ZSKO03kT4uq03geNakkdWFQ8/JUqDr0oVOl
KktBV9UaanWrXO2qV78K1rCKdaxkLatZz4rWtKp1rWxtq1vfCte4ynWudK2rXe+K17yu6nWv
Xb1pKPzK18CmVaYlfWJn2hhVwSoWrIk9njvkqc/GLnayWi2gJolK1J1eCaCzNKRhV9pUyopW
rJfN5UnNidrDPjaSYWTdPkcLW7NaNqE3nO0eURpZWRLWtbCMrW9Jy9rbepapv0TmST/b298q
t6/B9aZwoXnZMhoytcMN7XKv20KZ3jOmSEHiKGHKXeMiF7vkletHy4ve9P51qeptr3vfC9/4
yne+9K2vpgIBADs=
--------------050008020502040502000001
Content-Type: text/plain;
  name="FC_estim.inp"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
  filename="FC_estim.inp"
#define PROTON62
#define PROTON60
TITLE
Faraday Cup - backscattered particles and secondaries
GLOBAL                                          1.        1.
DEFAULTS 
NEW-DEFA
#if PROTON62
BEAM          -0.062                           1.0       0.0      -1.0PROTON
#elif PROTON60
BEAM            -.06                           1.0       0.0      -1.0PROTON
#endif
BEAMPOS          0.0       0.0     -10.0
GEOBEGIN 
COMBNAME
     0    0
* Black body
SPH blkbody    0.0 0.0 0.0 100000.0
* VOID
SPH void       0.0 0.0 0.0 10000.0
* Cubic Cu target block
RPP target     -2.0 2.0 -2.0 2.0 0.0 4.0
END
* Black hole
BLKBODY      5 +blkbody -void
* Void around
Vacuum       5 +void -target
* Target
CuBlock      5 +target
END
GEOEND
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
*
* Material assignment
ASSIGNMA    BLCKHOLE   BLKBODY
ASSIGNMA      VACUUM    Vacuum
ASSIGNMA      COPPER   CuBlock
* Physics and transport
* --------------------------------------
*
* Enalbling the MSC option
* -----------------------------
MULSOPT                                                   1. 
1.GLOBHAD
*
* Sefl-shielding of copper bolck
* ---------------------------------
LOW-MAT       COPPER
*
* Electromagnetic Fluka - production cut - minimum  for e- / e+ and photons
* --------------------------------------------------------------------
EMFCUT     -0.000001 0.0000001              VACUUM  @LASTMAT 
PROD-CUT
*
* Setting the physics of electromagnetic interactions
* ---------------------------------------------------------
*
EMFCUT     -0.000001 0.0000001       70.    Vacuum  @LASTREG
*
* Threshold of production of secondary protons
* ---------------------------------------------------
PART-THR      -1.E-7    PROTON                           0.0
*
* Scoring
* ------------
*
* Energy deposition in Copper block
* ---------------------------------------
USRBIN           10.    ENERGY      -21.       2.0       2.0       4.0EneDep
USRBIN          -2.0      -2.0       0.0      200.      200.      200.&
* Charged Hadrons fluence
USRBIN           10.  HAD-CHAR      -22.        5.        5.        6.HadFLU
USRBIN           -5.       -5.       -5.      200.      200.      200.&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -31.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062       0.0      100.  0.214094       0.0       1.0&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -32.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062       0.0      100.  0.841787   0.24094       1.0&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -33.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062       0.0      100.  1.840302  0.841787       1.0&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -34.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062 0.0000001      100.  3.141593  1.840302       1.0&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -35.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062       0.0      100.  4.656977  3.141593       1.0&
* Back scattered proton fluence
* ---------------------------------
USRBDX          101.    PROTON      -36.   CuBlock    Vacuum 
16.00ProtBckS
USRBDX         0.062       0.0      100.  6.283185  4.656977       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -41.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  0.214094       0.0       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -42.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  0.841787  0.214094       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -43.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  1.840302  0.841787       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -44.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  3.141593  1.840302       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -45.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  4.656977  3.141593       1.0&
* Electron flucence - 'back to vacuum'
* ----------------------------------------
USRBDX          101.  ELECTRON      -46.   CuBlock    Vacuum 
16.0ElecBckS
USRBDX         0.062  0.000001      100.  6.283185  4.656977       1.0&
* Flucence (E, angular distribution) for eletrons and positons
* -------------------------------------------------------------------
USRBDX          101.     E+&E-      -47.   CuBlock    Vacuum 
16.0ElePoBck
USRBDX         0.062  0.000001      100.  0.841787       0.0       1.0&
* Flucence (E, angular distribution) for eletrons and positons
* -------------------------------------------------------------------
USRBDX          101.     E+&E-      -48.   CuBlock    Vacuum 
16.0ElePoBck
USRBDX         0.062  0.000001      100.  3.141593  0.841787       1.0&
* Flucence (E, angular distribution) for eletrons and positons
* -------------------------------------------------------------------
USRBDX          101.     E+&E-      -49.   CuBlock    Vacuum 
16.0ElePoBck
USRBDX         0.062  0.000001      100.  6.283185  3.141593       1.0&
* Charged particle fluence (back of copper block)
USRBDX          101.  HAD-CHAR      -51.   CuBlock    Vacuum 
16.0CharFlue
USRBDX         0.062       0.0      100.  0.841787       0.0       1.0&
* Charged particle fluence (back of copper block)
USRBDX          101.  HAD-CHAR      -52.   CuBlock    Vacuum 
16.0CharFlue
USRBDX         0.062       0.0      100.  3.141593  0.841787       1.0&
* Charged particle fluence (back of copper block)
USRBDX          101.  HAD-CHAR      -53.   CuBlock    Vacuum 
16.0CharFlue
USRBDX         0.062       0.0      100.  6.283185  3.141593       1.0&
*
* Linear energy scoring - 200 energy bins , 1 angular bin
* ---------------------------------------------------------------
USRBDX          101.    PROTON      -90.   CuBlock    Vacuum 
16.0OneABin
USRBDX         0.063       0.0      200.  6.283185       0.0        1.&
* Proton double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD        124.    PROTON      -92.   CuBlock    Vacuum 
1.0YielProt
USRYIELD       180.0       90.      10.0     0.062       0.0        3.&
* Prpoton double differ fluence - energy 200, angle 1 bin - corresponds 
to usrbdx
* --------------------------------------------------------------
USRYIELD       2401.    PROTON      -91.   CuBlock    Vacuum 
1.0YilProt
USRYIELD       0.062       0.0      200.      180.       90.        3.&
* Electron double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD        124.  ELECTRON      -61.   CuBlock    Vacuum       1.0YielEl
USRYIELD       180.0      90.0      10.0     0.062  0.000001        3.&
* El / Pos double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD        124.  POSITRON      -62.   CuBlock    Vacuum 
1.0YiPoAng
USRYIELD       180.0       90.      10.0     0.062  0.000001        3.&
USRYIELD       2401.  POSITRON      -62.   CuBlock    Vacuum       1.0YiPoEn
USRYIELD       0.062  0.000001      200.      180.       90.        3.&
* Charged Hadr double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD        124.  HAD-CHAR      -71.   CuBlock    Vacuum 
1.0YielChA
USRYIELD       180.0       90.      10.0     0.062       0.0        3.&
* Charged Hadr double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD       2401.  HAD-CHAR      -72.   CuBlock    Vacuum 
1.0YielChE
USRYIELD       0.062       0.0     200.0     180.0      90.0        3.&
* Charged Hadr double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD        124.  HAD-CHAR      -73.   CuBlock    Vacuum 
1.0YielChr2
USRYIELD       180.0      90.0      10.0     0.062       0.0        3.&
* Charged Hadr double differencial fluence (E,angle) WRT beam direction.
* ------------------------------------------------------------------------
USRYIELD       2401.  HAD-CHAR      -73.   CuBlock    Vacuum 
1.0YielCh2E
USRYIELD       0.062       0.0     200.0     180.0      90.0        3.&
* Helium fluence
AUXSCORE    USRYIELD  4-HELIUM            YielChr2
* Helium fluence
AUXSCORE    USRYIELD  4-HELIUM            YielCh2E
RANDOMIZ         1.0
START       1000000.                                     0.0
STOP
--------------050008020502040502000001--
Received on Mon Jul 04 2011 - 22:44:30 CEST

This archive was generated by hypermail 2.2.0 : Mon Jul 04 2011 - 22:44:30 CEST