From 757fb0837d4f51a3b56bd401be14bec288a89b38 Mon Sep 17 00:00:00 2001 From: Artur Glavic Date: Thu, 3 Oct 2019 11:56:33 +0200 Subject: [PATCH] Add new surface source in bunker wall and test it --- e2-estia.i | 25 ++-- geometry/estia_guides.i | 8 +- geometry/estia_surfaces.i | 1 + sources/bunker_wall.i | 244 ++++++++++++++++++++++++++++++++++++++ 4 files changed, 265 insertions(+), 13 deletions(-) create mode 100644 sources/bunker_wall.i diff --git a/e2-estia.i b/e2-estia.i index 2fd9293..3cea4d4 100644 --- a/e2-estia.i +++ b/e2-estia.i @@ -87,6 +87,11 @@ TR101 0.0 0.0 0.0 0.58778525 -0.80901699 0.0 & 0.0 0.0 1.0 1 c Transformation from SKADI model coordinate system here TR102 0 5.0 1.0 0.66913061 -0.74314483 0 0.74314483 0.66913061 0 0 0 1 1 +c Coordinate for surface source within bunker wall +c close to plane 206305 with same surface normal and center around beam +TR103 700.03 -977.88 -0.86 0.587785 -0.809017 0.000000 + 0.809017 0.587785 0.000000 + 0 0 1 read file=geometry/ess_monolith_transformations.i echo c ------------------------------------------------------- c ++++++++++++++++++++++ END ++++++++++++++++++++++++++++ @@ -123,15 +128,15 @@ c read file=tallies/BBG_f4.i echo c read file=tallies/BBG_f5.i echo c read file=tallies/inbunker_f4.i echo c read file=tallies/chopper_pit_doses.i echo -read file=tallies/bunker_feedthrough_f1.i echo -c read file=tallies/bunker_wall_f4.i echo +c read file=tallies/bunker_feedthrough_f1.i echo +read file=tallies/bunker_wall_f4.i echo c read file=tallies/bunker_wall_f5.i echo c read file=tallies/pre_selene1_mesh.i -c read file=tallies/e1_e2_f4.i echo +read file=tallies/e1_e2_f4.i echo c read file=tallies/e1_e2_f5.i echo c read file=tallies/shutter_f4.i echo c read file=tallies/shutter_f5.i echo -c read file=tallies/shutter_wall_f4.i echo +read file=tallies/shutter_wall_f4.i echo c read file=tallies/detector_f5.i echo c read file=tallies/selene1_walls_f4.i echo c read file=tallies/optical_cave_f5.i echo @@ -170,12 +175,12 @@ c --------------- PHYSICS CARDS -------------------------- c ------------------------------------------------------- c dbcn 375642321 j j j j j j j j j j j 15291711 c stop calculation either at XeY histories, 6.5 days or tally 15 error <0.001 -stop nps 1e9 CTME 4e5 f31 0.01 +stop nps 1e6 CTME 4e5 f31 0.01 rand gen=2 seed=29976355590901 $ default 19073486328125 c histp mode n p h / z $ neutron, photon, proton, pi+, pi0 -c imp:n,p,h,/,z 1 276r 0 0 1 589r -imp:n,p,h,/,z 1 276r 0 0 1 487r 0 0 0 1 0 87r 1 9r +imp:n,p,h,/,z 1 276r 0 0 1 589r +c imp:n,p,h,/,z 1 276r 0 0 1 487r 0 0 0 1 0 87r 1 9r cut:h,/,z 1e+08 1e-3 phys:n 3000 0 phys:p 3000 @@ -202,7 +207,9 @@ c * Source for simulation from proton beam (needs different WW and ESplit) c read file=sources/proton.i echo c * Source Term from Valentina for 2m adapted to size c read file=sources/e2_2m.i echo -read file=sources/e2_2m_mod.i echo +c read file=sources/e2_2m_mod.i echo +c * Estimated source within bunker wall to optimize collimator efficiency +read file=sources/bunker_wall.i echo c * McStas source terms for gamma doeses c read file=sources/mcstas_source_vs.i echo c read file=sources/mcstas_source_sample.i echo @@ -210,7 +217,7 @@ c read file=sources/mcstas_source_mf.i echo c read file=sources/mcstas_gammas.i echo c read file=sources/mcstas_gammas_2.i echo c read file=sources/mcstas_source_shutter.i echo -prdmp 1e8 -30 1 2 1e7 +prdmp 1e6 -30 1 2 1e5 print 30 32 40 c 120 $O importance function quality c 128 $O universe map diff --git a/geometry/estia_guides.i b/geometry/estia_guides.i index beacbdb..edd39ca 100644 --- a/geometry/estia_guides.i +++ b/geometry/estia_guides.i @@ -18,10 +18,10 @@ c #(-5102 5501 5502 -5505 5504) c #(5102 -7100 -5531) 206002 4 -1.2E-3 206305 -207205 -206401 -206402 -206403 -206404 #207022 #207024 #207025 #207031 #207032 #207033 -c fill=50 (703.47 -983.22 -0.89 -c 0.58207953 -0.81304008 0.0 -c -0.81310076 -0.58212297 0.0 -c 0.0 0.0 1.0) + fill=50 (703.47 -983.22 -0.89 + 0.58207953 -0.81304008 0.0 + -0.81310076 -0.58212297 0.0 + 0.0 0.0 1.0) 206003 4 -1.2E-3 207205 -207222 -206401 -206402 -206403 -206404 #(207201 -207202 -207203 207204) 206004 4 -1.2E-3 207206 -206405 -206401 -206402 -206403 -206404 208921 diff --git a/geometry/estia_surfaces.i b/geometry/estia_surfaces.i index 8cfc6c4..7924d23 100644 --- a/geometry/estia_surfaces.i +++ b/geometry/estia_surfaces.i @@ -381,6 +381,7 @@ c Segmentation surfaces for bunker feedthrough f1 tally 206702 4 c/x 3.0 0.0 7 206703 4 c/x 3.0 0.0 13 206704 4 c/x 3.0 0.0 30 +c 206705 103 rpp 0 2 -10 10 -10 10 $ to show location of surface cell TR103 c Cu-collimator c Radii 207102 cz 1271.50 diff --git a/sources/bunker_wall.i b/sources/bunker_wall.i new file mode 100644 index 0000000..33e8376 --- /dev/null +++ b/sources/bunker_wall.i @@ -0,0 +1,244 @@ +c source definition based on tally_181.txt, modified for high E flux +cut:n j 0 +c cut all photons that can't produce neutrons +cut:p 1e+08 1.667 +c but adapted to model geometry and actual size of Estia beam extraction +c make sure to supply the correct WWINP file +WWP:n j j 4 j -1 j 1.0 1j 1.0E+8 1j +c the source definition with neutrons/s +sdef tr=103 pos=0 0 0 vec=1 0 0 + dir=d1 rad=fdir d2 erg=fdir d3 + par=n ara=11309.7 wgt=2.0e+8 +c distribution of radial direction +si1 0 + 5.00000E-01 8.66025E-01 9.65926E-01 + 9.96195E-01 9.99391E-01 1.00000E+00 +sp1 0.00000000e+00 + 4.61452024e+07 1.70198257e+08 1.08329950e+08 + 6.38706911e+07 9.10952408e+06 2.26490368e+06 +c directional dependent distribution of position +ds2 q 0 21 + 5.00000E-01 22 + 8.66025E-01 23 + 9.65926E-01 24 + 9.96195E-01 25 + 9.99391E-01 26 + 1.00000E+00 27 +c Energy dependent probability for each direction pair +ds3 q 0 31 + 5.00000E-01 32 + 8.66025E-01 33 + 9.65926E-01 34 + 9.96195E-01 35 + 9.99391E-01 36 + 1.00000E+00 37 +c direction probabilities for each segment +si21 0 1 +sp21 0 1.0000E-20 +c +si22 0 2 7 13 30 122.51 +sp22 0 2.30230E+04 2.07305E+05 5.25726E+05 1.23189E+06 2.10846E+07 +c +si23 0 2 7 13 30 122.51 +sp23 0 6.76418E+04 1.27453E+06 2.04084E+06 6.08084E+06 7.56352E+07 +c +si24 0 2 7 13 30 122.51 +sp24 0 8.63986E+04 7.45148E+05 1.24766E+06 2.35238E+06 4.97334E+07 +c +si25 0 2 7 13 30 122.51 +sp25 0 2.43288E+04 4.55438E+05 5.77578E+05 1.32580E+06 2.95522E+07 +c +si26 0 2 7 13 30 122.51 +sp26 0 1.63116E+05 1.04990E+05 1.76650E+05 4.00354E+05 3.70965E+06 +c +si27 0 2 7 13 30 122.51 +sp27 0 1.68399E+04 4.20802E+04 9.41964E+04 8.99954E+04 8.89340E+05 +c energy dependent probability for each direction +si31 0 1 +sp31 0 1.0000E-20 +C ------------------------------------------------------------------------ +si32 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp32 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 3.14844E+05 1.50026E+05 0.00000E+00 0.00000E+00 0.00000E+00 + 5.04168E+05 0.00000E+00 2.47386E+05 0.00000E+00 0.00000E+00 + 1.19420E+05 0.00000E+00 0.00000E+00 2.83700E+05 5.18008E+05 + 4.50826E+05 1.17780E+05 0.00000E+00 0.00000E+00 3.24702E+05 + 0.00000E+00 3.04862E+05 5.05966E+04 1.27270E+05 3.51030E+05 + 3.12144E+05 2.54320E+05 4.39610E+05 1.89490E+05 3.07168E+05 + 5.67617E+05 5.24070E+05 6.84092E+05 1.08998E+06 1.05406E+06 + 1.56784E+06 2.06740E+06 1.63649E+06 1.28594E+06 9.86652E+05 + 9.00414E+05 8.92366E+05 1.09287E+06 1.32177E+06 1.10393E+06 + 6.95039E+05 1.87933E+05 4.50658E+04 1.73907E+03 0.00000E+00 + 0.00000E+00 0.00000E+00 2.30726E+07 +c +si33 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp33 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 8.96446E+04 + 1.24880E+05 2.83058E+05 2.51954E+05 6.87496E+05 7.15686E+05 + 6.94112E+05 9.17851E+05 1.41215E+05 6.94870E+04 3.30944E+05 + 2.24454E+05 1.15589E+06 3.71718E+05 6.70040E+05 6.24094E+05 + 9.71906E+05 7.78151E+05 3.52868E+05 2.61442E+05 2.02862E+05 + 3.12048E+05 4.24589E+05 4.77852E+05 1.01942E+06 3.71538E+05 + 5.74270E+05 5.68819E+05 7.53604E+05 9.31588E+05 5.57942E+05 + 8.05481E+05 8.61544E+05 8.70714E+05 8.55040E+05 9.96667E+05 + 1.90695E+06 1.83161E+06 2.15864E+06 2.40250E+06 4.03989E+06 + 3.85865E+06 4.51882E+06 4.32260E+06 3.50218E+06 2.75366E+06 + 2.17294E+06 2.62505E+06 4.26290E+06 6.18216E+06 7.27704E+06 + 6.33514E+06 3.45029E+06 1.01333E+06 1.69420E+05 1.41528E+04 + 3.21214E+02 0.00000E+00 8.50991E+07 +c +si34 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp34 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 3.50044E+05 2.11036E+05 7.25122E+05 4.01714E+05 + 2.85638E+05 0.00000E+00 2.89634E+05 5.13630E+05 2.83914E+05 + 7.57076E+05 1.06405E+05 1.20970E+05 3.91716E+05 1.32458E+05 + 0.00000E+00 1.28368E+05 0.00000E+00 3.83210E+05 3.27540E+05 + 4.95046E+05 4.01546E+05 1.45736E+05 3.95288E+05 4.07088E+05 + 2.15020E+05 2.84168E+05 3.56171E+05 2.78774E+05 4.20790E+05 + 5.99099E+05 4.42646E+05 1.72771E+05 4.67428E+05 3.10148E+05 + 6.13976E+05 8.44794E+05 9.57095E+05 1.36761E+06 1.50384E+06 + 1.43597E+06 2.37344E+06 1.74013E+06 1.49268E+06 1.20460E+06 + 8.77910E+05 1.15478E+06 1.96657E+06 4.33218E+06 6.25133E+06 + 6.14042E+06 5.83465E+06 2.37386E+06 7.54202E+05 1.23595E+05 + 1.91279E+04 1.98462E+03 5.41650E+07 +c +si35 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp35 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 5.80822E+04 + 7.82170E+04 0.00000E+00 0.00000E+00 0.00000E+00 1.18801E+05 + 2.03312E+05 9.17614E+04 0.00000E+00 0.00000E+00 3.13736E+05 + 0.00000E+00 2.85096E+04 1.01698E+05 0.00000E+00 0.00000E+00 + 4.15596E+04 1.82898E+05 0.00000E+00 1.31384E+05 4.61576E+04 + 1.54784E+05 3.45572E+04 1.52229E+05 0.00000E+00 1.33941E+04 + 1.05392E+05 1.80908E+05 2.09700E+05 8.77208E+04 3.44191E+05 + 1.95796E+05 2.88313E+05 3.23630E+05 4.55422E+05 6.71276E+05 + 5.39635E+05 8.32866E+05 2.83956E+05 5.40580E+05 4.28812E+05 + 3.37518E+05 3.94562E+05 7.58188E+05 1.78063E+06 3.63320E+06 + 4.66280E+06 7.41933E+06 3.86641E+06 1.47653E+06 2.91329E+05 + 5.72558E+04 1.82814E+04 3.19354E+07 +c +si36 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp36 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 1.14870E+05 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 1.77094E+03 1.95123E+03 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 6.35549E+04 2.26664E+04 0.00000E+00 0.00000E+00 + 2.35992E+04 7.97576E+04 5.61072E+04 3.09262E+04 4.60783E+04 + 5.68676E+04 9.41971E+04 2.47496E+04 1.27729E+04 6.29456E+04 + 3.05540E+04 9.95486E+04 5.93783E+04 2.44680E+05 5.31452E+05 + 5.44316E+05 1.09089E+06 7.82100E+05 3.84272E+05 7.48612E+04 + 1.54286E+04 4.46045E+03 4.55476E+06 +c +si37 0.00000E+00 + 1.00000E-10 1.58000E-10 2.51000E-10 3.98000E-10 6.31000E-10 + 1.00000E-09 1.58000E-09 2.51000E-09 3.98000E-09 6.31000E-09 + 1.00000E-08 1.58000E-08 2.51000E-08 3.98000E-08 6.31000E-08 + 1.00000E-07 1.58000E-07 2.51000E-07 3.98000E-07 6.31000E-07 + 1.00000E-06 1.58000E-06 2.51000E-06 3.98000E-06 6.31000E-06 + 1.00000E-05 1.58000E-05 2.51000E-05 3.98000E-05 6.31000E-05 + 1.00000E-04 1.58000E-04 2.51000E-04 3.98000E-04 6.31000E-04 + 1.00000E-03 1.58000E-03 2.51000E-03 3.98000E-03 6.31000E-03 + 1.00000E-02 1.58000E-02 2.51000E-02 3.98000E-02 6.31000E-02 + 1.00000E-01 1.58000E-01 2.51000E-01 3.98000E-01 6.31000E-01 + 1.00000E+00 1.58000E+00 2.51000E+00 3.98000E+00 6.31000E+00 + 1.00000E+01 1.58000E+01 2.51000E+01 3.98000E+01 6.31000E+01 + 1.00000E+02 1.58000E+02 2.51000E+02 3.98000E+02 6.31000E+02 + 1.00000E+03 1.58000E+03 2.51000E+03 +sp37 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 5.63064E-10 0.00000E+00 0.00000E+00 0.00000E+00 4.99236E-02 + 0.00000E+00 0.00000E+00 1.17927E-04 0.00000E+00 1.87258E+05 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 1.51526E-02 0.00000E+00 0.00000E+00 1.82373E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 1.94199E+03 0.00000E+00 4.92952E+04 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 2.85002E+04 + 0.00000E+00 2.37029E+04 0.00000E+00 1.65305E+04 0.00000E+00 + 1.11586E+04 7.41496E+04 8.46998E+03 5.06353E+03 5.09402E+01 + 9.70764E+01 1.31903E+04 2.21342E+04 3.09494E+04 5.67470E+04 + 9.29058E+04 2.40455E+05 1.69563E+05 7.28883E+04 2.05490E+04 + 6.57676E+03 2.72816E+02 1.13245E+06