W-8．JASPERコンクリート体系ギャップストリーミング実験 問題識別番号：PSMC-SH-B04 　　　　 　　 　　（黒澤正彦） sourceルーチンのFORTRANソースと入力データ (1)　subroutine　source c c common/jasper/istatus,c(51),d(51,8) dimension c(51),d(51,8) dimension anf(51,8) dimension eng(52) dimension an8(8),wt8(8),ab8(9) c data mtot/8/ data an8/0.0950124,0.281603, 0.458017, 0.617876, 1 0.755404, 0.865631, 0.944574, 0.989401/ data ab8/0.0, 0.1894509, 0.3720543, 0.5412108, 0.6908067, 1 0.8154356, 0.9105941, 0.9728476, 1.0000000/ data wt8/0.1894506, 0.1826034, 0.1691565, 0.1495959, 1 0.1246289, 0.0951585, 0.0622535, 0.0271524/ c data ntot/51/ data eng/14.918, 12.214, 10.00, 8.187, 6.703, 1 5.488, 4.493, 3.679, 3.012, 2.466, 2 2.019, 1.653, 1.353, 1.1080, 0.9072, 3 0.7427, 0.6081, 0.4978, 0.4076, 0.3337, 4 0.2732, 0.2237, 0.1832, 0.1500, 0.12277, 5 8.652d-2, 5.248d-2, 4.086d-2, 3.183d-2, 2.479d-2, 6 1.930d-2, 1.503d-2, 7.102d-3, 4.307d-3, 3.355d-3, 7 2.613d-3, 2.035d-3, 1.585d-3, 1.234d-3, 9.611d-4, 8 4.540d-4, 2.145d-4, 1.013d-4, 4.785d-5, 2.260d-5, 9 1.068d-5, 5.043d-6, 2.382d-6, 1.125d-6, 4.140d-7, a 1.000d-7, 1.000d-9/ c source data was taken from ORNL-5259 data (anf(i,1),i=1,51)/ 1 2.91d+2, 1.55d+3, 4.95d+3, 1.91d+4, 6.39d+4, 2 1.64d+5, 2.88d+5, 4.19d+5, 6.80d+5, 7.64d+5, 3 8.63d+5, 8.23d+5, 7.62d+5, 7.96d+5, 6.49d+5, 4 7.32d+5, 6.64d+5, 4.68d+5, 5.54d+5, 3.67d+5, 5 3.91d+5, 3.29d+5, 1.62d+5, 5.25d+5, 1.93d+5, 6 6.76d+5, 1.53d+5, 7.81d+4, 3.36d+5, 6.99d+4, 7 7.29d+4, 3.08d+5, 2.23d+5, 1.02d+5, 9.94d+4, 8 1.00d+5, 1.01d+5, 9.93d+4, 9.79d+4, 2.91d+5, 9 2.91d+5, 2.88d+5, 2.88d+5, 2.91d+5, 2.91d+5, a 3.54d+5, 3.51d+5, 4.19d+5, 8.34d+5, 1.40d+6, b 9.43d+6/ data (anf(i,2),i=1,51)/ 1 9.65d+2, 4.73d+3, 1.49d+4, 5.27d+4, 1.60d+5, 2 3.78d+5, 6.28d+5, 8.92d+5, 1.39d+6, 1.56d+6, 3 1.68d+6, 1.59d+6, 1.43d+6, 1.53d+6, 1.16d+6, 4 1.37d+6, 1.19d+6, 8.14d+5, 1.06d+6, 6.40d+5, 5 7.73d+5, 5.50d+5, 2.49d+5, 9.31d+5, 3.15d+5, 6 1.41d+6, 2.90d+5, 1.08d+5, 9.36d+5, 2.25d+5, 7 2.04d+5, 7.70d+5, 5.17d+5, 2.46d+5, 2.38d+5, 8 2.37d+5, 2.37d+5, 2.34d+5, 2.30d+5, 6.83d+5, 9 6.82d+5, 6.75d+5, 6.78d+5, 6.82d+5, 6.83d+5, a 8.29d+5, 8.18d+5, 9.80d+5, 1.96d+6, 3.26d+6, b 2.18d+7/ data (anf(i,3),i=1,51)/ 1 2.22d+3, 1.04d+4, 3.05d+4, 9.75d+4, 2.58d+5, 2 5.43d+5, 8.29d+5, 1.14d+6, 1.67d+6, 1.78d+6, 3 1.81d+6, 1.74d+6, 1.61d+6, 1.61d+6, 1.34d+6, 4 1.50d+6, 1.33d+6, 9.04d+5, 1.17d+6, 7.48d+5, 5 8.93d+5, 6.30d+5, 3.00d+5, 1.02d+6, 3.98d+5, 6 1.61d+6, 3.78d+5, 1.27d+5, 1.07d+6, 4.51d+5, 7 3.86d+5, 1.26d+6, 7.69d+5, 3.77d+5, 3.64d+5, 8 3.61d+5, 3.60d+5, 3.51d+5, 3.46d+5, 1.03d+6, 9 1.02d+6, 1.01d+6, 1.01d+6, 1.01d+6, 1.01d+6, a 1.23d+6, 1.22d+6, 1.46d+6, 2.93d+6, 4.93d+6, b 3.42d+7/ data (anf(i,4),i=1,51)/ 1 6.99d+3, 3.09d+4, 8.66d+4, 2.43d+5, 5.44d+5, 2 9.92d+5, 1.34d+6, 1.70d+6, 2.38d+6, 2.35d+6, 3 2.26d+6, 2.13d+6, 1.96d+6, 1.81d+6, 1.63d+6, 4 1.74d+6, 1.53d+6, 1.05d+6, 1.31d+6, 9.16d+5, 5 1.02d+6, 7.53d+5, 3.88d+5, 1.07d+6, 5.39d+5, 6 1.68d+6, 4.84d+5, 1.57d+5, 8.61d+5, 5.52d+5, 7 4.90d+5, 1.60d+6, 9.80d+5, 4.79d+5, 4.63d+5, 8 4.59d+5, 4.56d+5, 4.49d+5, 4.42d+5, 1.32d+6, 9 1.30d+6, 1.29d+6, 1.29d+6, 1.30d+6, 1.30d+6, a 1.58d+6, 1.55d+6, 1.89d+6, 3.79d+6, 6.45d+6, b 4.60d+7/ data (anf(i,5),i=1,51)/ 1 2.13d+4, 9.39d+4, 2.62d+5, 6.76d+5, 1.37d+6, 2 2.23d+6, 2.60d+6, 2.98d+6, 4.15d+6, 3.89d+6, 3 3.43d+6, 3.06d+6, 2.61d+6, 2.28d+6, 2.03d+6, 4 2.15d+6, 1.86d+6, 1.19d+6, 1.51d+6, 1.07d+6, 5 1.18d+6, 8.58d+5, 4.66d+5, 1.14d+6, 6.59d+5, 6 1.85d+6, 5.53d+5, 1.85d+5, 9.93d+5, 5.78d+5, 7 5.21d+5, 1.75d+6, 1.09d+6, 5.33d+5, 5.19d+5, 8 5.17d+5, 5.15d+5, 5.07d+5, 5.00d+5, 1.50d+6, 9 1.48d+6, 1.47d+6, 1.48d+6, 1.50d+6, 1.50d+6, a 1.83d+6, 1.82d+6, 2.17d+6, 4.37d+6, 7.42d+6, b 5.21d+7/ data (anf(i,6),i=1,51)/ 1 3.62d+4, 1.62d+5, 4.51d+5, 1.16d+6, 2.30d+6, 2 3.62d+6, 4.15d+6, 4.55d+6, 6.25d+6, 5.69d+6, 3 4.81d+6, 4.17d+6, 3.38d+6, 2.78d+6, 2.52d+6, 4 2.60d+6, 2.21d+6, 1.35d+6, 1.71d+6, 1.23d+6, 5 1.32d+6, 9.70d+5, 5.47d+5, 1.22d+6, 7.70d+5, 6 1.96d+6, 6.19d+5, 2.16d+5, 9.79d+5, 6.48d+5, 7 5.87d+5, 1.95d+6, 1.21d+6, 5.96d+5, 5.80d+5, 8 5.78d+5, 5.77d+5, 5.66d+5, 5.63d+5, 1.67d+6, 9 1.67d+6, 1.65d+6, 1.67d+6, 1.68d+6, 1.68d+6, a 2.06d+6, 2.06d+6, 2.45d+6, 4.93d+6, 8.40d+6, b 5.92d+7/ data (anf(i,7),i=1,51)/ 1 3.76d+4, 1.71d+5, 4.73d+5, 1.25d+6, 2.58d+6, 2 4.17d+6, 5.00d+6, 5.57d+6, 7.60d+6, 6.89d+6, 3 5.89d+6, 5.08d+6, 4.07d+6, 3.28d+6, 2.94d+6, 4 3.02d+6, 2.51d+6, 1.48d+6, 1.88d+6, 1.36d+6, 5 1.44d+6, 1.06d+6, 6.06d+5, 1.27d+6, 8.50d+5, 6 2.07d+6, 6.61d+5, 2.38d+5, 1.03d+6, 6.66d+5, 7 6.09d+5, 2.04d+6, 1.28d+6, 6.29d+5, 6.13d+5, 8 6.10d+5, 6.09d+5, 6.02d+5, 5.95d+5, 1.78d+6, 9 1.78d+6, 1.78d+6, 1.78d+6, 1.81d+6, 1.81d+6, a 2.20d+6, 2.20d+6, 2.62d+6, 5.26d+6, 9.02d+6, b 6.25d+7/ data (anf(i,8),i=1,51)/ 1 3.70d+4, 1.71d+5, 4.75d+5, 1.28d+6, 2.68d+6, 2 4.40d+6, 5.44d+6, 6.10d+6, 8.27d+6, 7.53d+6, 3 6.47d+6, 5.59d+6, 4.46d+6, 3.56d+6, 3.21d+6, 4 3.28d+6, 2.70d+6, 1.57d+6, 1.97d+6, 1.44d+6, 5 1.53d+6, 1.11d+6, 6.44d+5, 1.31d+6, 8.97d+5, 6 2.13d+6, 6.93d+5, 2.52d+5, 1.04d+6, 6.94d+5, 7 6.81d+5, 2.11d+6, 1.33d+6, 6.55d+5, 6.40d+5, 8 6.38d+5, 6.36d+5, 6.29d+5, 6.22d+5, 1.85d+6, 9 1.85d+6, 1.83d+6, 1.85d+6, 1.86d+6, 1.88d+6, a 2.30d+6, 2.28d+6, 2.73d+6, 5.50d+6, 9.37d+6, b 6.54d+7/ c source data was taken from ORNL-5259 c c ceutrons/cm**n/min/watt/unit weight c data statements can initialize variables in common blocks c data istatus/1/ data socrad/48.575/ c cumulative will be only calculated at the first call if(istatus .eq. 0) go to 1 c unnormalized pdf(l,m) by each energy group l and angle m do 10 m=1,mtot do 10 n=1,ntot c 10 d(n,m)=anf(n,m)*an8(m)*wt8(m)*0.5 10 d(n,m)=anf(n,m)*an8(m)*wt8(m)*0.5*0.85 c factor 0.5 is due to forward direction only c factor 0.85 is for the sphere source cndition (ORNL-TM-5183) c 10 d(n,m)=anf(n,m)*an8(m)*0.5 c need not consider wt8(m) in this calculation : Mr Takemura, in April 17,1997 c summation over direcions and energy groups a=0.0 do 11 n=1,ntot sum=0.0 do 12 m=1,mtot 12 sum=sum+d(n,m) c(n)=sum 11 a=a+sum c cummulative distribution function for d(n,m) and c(n) cum=0.0 do 13 n=1,ntot acum=0.0 do 14 m=1,mtot acum=acum+d(n,m)/c(n) 14 d(n,m)=acum cum=cum+c(n)/a 13 c(n)=cum c write total a write(6,600) a c write pdf(c(m) and d(m,n) for each energy group m and direction n do 15 n=1,ntot 15 write(6,666) n,eng(n),c(n),(d(n,m),m=1,mtot) c set flag when cdf was calculated istatus=0 c format 600 format(' Source Energy and Angular Distributions. Total:', *1pd12.5,/' energy bin(MeV) Energy Angular(1 to 8) cumulative *distribution') 666 format(i5,' ',1pd11.4,0p9f10.7) c 1 continue c c select source energy c gzai = rang() do 20 n=1,ntot if(c(n) .gt. gzai) go to 21 20 continue n=ntot 21 nn=n erg=eng(nn+1)+rang()*(eng(nn)-eng(nn+1)) c c select source point(x,y,z) www1 = 2.0*rang()-1.0 if(abs(www1) .eq. 1.0) then zzz = socrad yyy = 0.0 xxx = 0.0 else c row1 = sqrt(1.0-www1**2) phai = pie*(2.0*rang()-1.0) vvv1 = row1*sin(phai) uuu1 = row1*cos(phai) zzz = socrad*www1 yyy = socrad*vvv1 xxx = socrad*uuu1 endif c c select direction gzai = rang() do 30 m=1,mtot if(d(nn,m) .gt. gzai) go to 31 30 continue m=mtot 31 mm=m eta=ab8(mm)+rang()*(ab8(mm+1)-ab8(mm)) if(eta .eq. 1.0) then uuu=0.0 vvv=0.0 www=www1 else row2=sqrt(1.0-eta**2) theta = pie*(2.0*rang()-1.0) ccc = cos(theta) ddd = sin(theta) if(www1 .eq. 1.0) then uuu = row2*ccc vvv = row2*ddd www = eta*www1 else uuu = row2*(ccc*www1*uuu1 - ddd*vvv1)/row1 + eta*uuu1 vvv = row2*(ccc*www1*vvv1 + ddd*uuu1)/row1 + eta*vvv1 www =-row2*ccc*row1 + eta*www1 endif endif u2v2w2= uuu**2 + vvv**2 + www**2 c wgt = 1.0 tme = 0.0 psc = 0.5 icl = 4 jsu = 0 c uuu = uuu vvv = vvv www = www xxx = xxx yyy = yyy zzz = zzz erg = erg c c call expirx(0,'source','you need a source subroutine.') return end (2)　入力データ例　(体系U-D) JASPER program Gap streaming Experiment II-D 1 2 1.0038e-1 -1 $ core 2 1 5.9643e-2 1 -2 $ alminum 3 2 1.0038e-1 2 -3 $ water 4 1 5.9643e-2 3 -4 $ alminum - source surface 5 999 4.9636e-5 4 -5 100 -101 $ air around core cylinder 6 999 4.9636e-5 4 -6 -100 $ air around core sphere 7 3 8.5853e-2 5 -7 100 -101 $ ss-304 cylinder 8 3 8.5853e-2 6 -8 -100 $ ss-304 sphere 9 4 9.2365e-2 7 -9 100 -101 $ ss-304 + water cylinder 10 4 9.2365e-2 8 -10 -100 $ ss-304 + water sphere 11 5 7.7175e-2 9 -11 100 -101 $ concrete cylinder 12 5 7.7175e-2 10 -11 -100 $ concrete sphere c 13 999 4.9636e-5 -11 101 -102 $ 1.27cm air gap 14 1 5.9643e-2 12 -13 102 -103 $ alminum or carbon steel 15 1 5.9643e-2 14 -15 104 -105 $ alminum or carbon steel 16 1 5.9643e-2 12 -15 103 -104 $ alminum or carbon steel 17 5 7.7175e-2 -11 13 102 -103 $ concrete in the first 18 5 7.7175e-2 -11 15 103 -105 $ concrete in the second 19 999 4.9636e-5 -12 102 -104 $ air in the first collimator 20 999 4.9636e-5 -14 104 -105 $ air in the second collimator 21 999 4.9636e-5 -11 105 -106 $ air in front of the alminum facing 22 999 4.9636e-5 -16 106 -107 $ air 23 1 5.9643e-2 -11 16 106 -107 $ alminum facing 5.08cm c 24 999 4.9636e-5 -11 107 -108 $ air in front of the shield c 25 6 8.5668e-2 -17 $ iron slab in front of the concrete vessel c 26 6 8.5668e-2 -205 109 -110 $ 4.88cm iron holder of cylinder A 27 999 4.9636e-5 -204 205 109 -110 $ air gap out of X3 28 6 8.5668e-2 -203 204 109 -110 $ iron liner of concrete vessel 29 8 8.1809e-2 -19 203 109 -110 $ concrete vessel out of cylinder c 30 8 8.1809e-2 -206 110 -300 $ cylinder A 31 6 8.5668e-2 -205 206 110 -300 $ X3 and iron linner of cylinder A 32 999 4.9636e-5 -204 205 110 -300 $ air gap out of X3 33 6 8.5668e-2 -203 204 110 -300 $ iron linner of concrete vessel 34 8 8.1809e-2 -19 203 110 -300 $ concrete vessel out of cylinder c 35 8 8.1809e-2 -206 300 -301 $ cylinder A 36 6 8.5668e-2 -205 206 300 -301 $ X3 and iron linner of cylinder A 37 999 4.9636e-5 -204 205 300 -301 $ air gap out of X3 38 6 8.5668e-2 -203 204 300 -301 $ iron linner of concrete vessel 39 8 8.1809e-2 -19 203 300 -301 $ concrete vessel out of cylinder c 40 8 8.1809e-2 -206 301 -302 $ cylinder A 41 6 8.5668e-2 -205 206 301 -302 $ X3 and iron linner of cylinder A 42 999 4.9636e-5 -204 205 301 -302 $ air gap out of X3 43 6 8.5668e-2 -203 204 301 -302 $ iron linner of concrete vessel 44 8 8.1809e-2 -19 203 301 -302 $ concrete vessel out of cylinder c 45 8 8.1809e-2 -206 302 -303 $ cylinder A 46 6 8.5668e-2 -205 206 302 -303 $ X3 and iron linner of cylinder A 47 999 4.9636e-5 -204 205 302 -303 $ air gap out of X3 48 6 8.5668e-2 -203 204 302 -303 $ iron linner of concrete vessel 49 8 8.1809e-2 -19 203 302 -303 $ concrete vessel out of cylinder c 50 8 8.1809e-2 -206 303 -304 $ cylinder A 51 6 8.5668e-2 -205 206 303 -304 $ X3 and iron linner of cylinder A 52 999 4.9636e-5 -204 205 303 -304 $ air gap out of X3 53 6 8.5668e-2 -203 204 303 -304 $ iron linner of concrete vessel 54 8 8.1809e-2 -19 203 303 -304 $ concrete vessel out of cylinder c 55 8 8.1809e-2 -206 304 -305 $ cylinder A 56 6 8.5668e-2 -205 206 304 -305 $ X3 and iron linner of cylinder A 57 999 4.9636e-5 -204 205 304 -305 $ air gap out of X3 58 6 8.5668e-2 -203 204 304 -305 $ iron linner of concrete vessel 59 8 8.1809e-2 -19 203 304 -305 $ concrete vessel out of cylinder c 60 8 8.1809e-2 -206 305 -306 $ cylinder A 61 6 8.5668e-2 -205 206 305 -306 $ X3 and iron linner of cylinder A 62 999 4.9636e-5 -204 205 305 -306 $ air gap out of X3 63 6 8.5668e-2 -203 204 305 -306 $ iron linner of concrete vessel 64 8 8.1809e-2 -19 203 305 -306 $ concrete vessel out of cylinder c 65 8 8.1809e-2 -206 306 -307 $ cylinder A 66 6 8.5668e-2 -205 206 306 -307 $ X3 and iron linner of cylinder A 67 999 4.9636e-5 -204 205 306 -307 $ air gap out of X3 68 6 8.5668e-2 -203 204 306 -307 $ iron linner of concrete vessel 69 8 8.1809e-2 -19 203 306 -307 $ concrete vessel out of cylinder c 70 8 8.1809e-2 -206 307 -111 $ cylinder A 71 6 8.5668e-2 -205 206 307 -111 $ X3 and iron linner of cylinder A 72 999 4.9636e-5 -204 205 307 -111 $ air gap out of X3 73 6 8.5668e-2 -203 204 307 -111 $ iron linner of concrete vessel 74 8 8.1809e-2 -19 203 307 -111 $ concrete vessel out of cylinder c 75 8 8.1809e-2 -206 111 -112 $ cylinder A 76 6 8.5668e-2 -205 206 111 -112 $ X3 and iron linner of cylinder A 77 999 4.9636e-5 -204 205 111 -112 $ air gap out of X3 78 6 8.5668e-2 -203 204 111 -112 $ iron linner of concrete vessel 79 6 8.5668e-2 -19 203 111 -112 150 $ iron linner of three fourths vessel 80 8 8.1809e-2 -19 203 111 -112 -150 $ concrete vessel out of cylinder c 81 8 8.1809e-2 -206 112 -113 $ cylinder A 82 6 8.5668e-2 -202 206 112 -113 $ X4 and iron linner of cylinder A 83 999 4.9636e-5 -201 202 112 -113 $ air gap out of X4 84 6 8.5668e-2 -200 201 112 -113 $ iron linner of concrete vessel 85 6 8.5668e-2 -19 200 112 -113 150 $ iron linner of three fourths vessel 86 8 8.1809e-2 -19 200 112 -113 -150 $ concrete vessel out of cylinder c 87 8 8.1809e-2 -206 113 -309 $ cylinder A 88 6 8.5668e-2 -202 206 113 -309 $ X4 and iron linner of cylinder A 89 999 4.9636e-5 -201 202 113 -309 $ air gap out of cylinder X4 90 6 8.5668e-2 -200 201 113 -309 $ X5 and iron linner of concrete vessel 91 8 8.1809e-2 -19 200 113 -309 151 $ concrete vessel out of cylinder 92 6 8.5668e-2 -19 200 113 -309 -151 150 $ concrete vessel out of cylinder 93 8 8.1809e-2 -19 200 113 -309 -150 $ concrete vessel out of cylinder c 94 8 8.1809e-2 -206 309 -310 $ cylinder A 95 6 8.5668e-2 -202 206 309 -310 $ X4 and iron linner of cylinder A 96 999 4.9636e-5 -201 202 309 -310 $ air gap out of cylinder X4 97 6 8.5668e-2 -200 201 309 -310 $ X5 and iron linner of concrete vessel 98 8 8.1809e-2 -19 200 309 -310 151 $ concrete vessel out of cylinder 99 6 8.5668e-2 -19 200 309 -310 -151 150 $ concrete vessel out of cylinder 100 8 8.1809e-2 -19 200 309 -310 -150 $ concrete vessel out of cylinder c 101 8 8.1809e-2 -206 310 -311 $ cylinder A 102 6 8.5668e-2 -202 206 310 -311 $ X4 and iron linner of cylinder A 103 999 4.9636e-5 -201 202 310 -311 $ air gap out of cylinder X4 104 6 8.5668e-2 -200 201 310 -311 $ X5 and iron linner of concrete vessel 105 8 8.1809e-2 -19 200 310 -311 151 $ concrete vessel out of cylinder 106 6 8.5668e-2 -19 200 310 -311 -151 150 $ concrete vessel out of cylinder 107 8 8.1809e-2 -19 20 310 -311 -150 $ concrete vessel out of cylinder c 108 8 8.1809e-2 -206 311 -312 $ cylinder A 109 6 8.5668e-2 -202 206 311 -312 $ X4 and iron linner of cylinder A 110 999 4.9636e-5 -201 202 311 -312 $ air gap out of cylinder X4 111 6 8.5668e-2 -200 201 311 -312 $ X5 and iron linner of concrete vessel 112 8 8.1809e-2 -19 200 311 -312 151 $ concrete vessel out of cylinder 113 6 8.5668e-2 -19 200 311 -312 -151 150 $ concrete vessel out of cylinder 114 8 8.1809e-2 -19 200 311 -312 -150 $ concrete vessel out of cylinder c 115 8 8.1809e-2 -206 312 -313 $ cylinder A 116 6 8.5668e-2 -202 206 312 -313 $ X4 and iron linner of cylinder A 117 999 4.9636e-5 -201 202 312 -313 $ air gap out of cylinder X4 118 6 8.5668e-2 -200 201 312 -313 $ X5 and iron linner of concrete vessel 119 8 8.1809e-2 -19 200 312 -313 151 $ concrete vessel out of cylinder 120 6 8.5668e-2 -19 200 312 -313 -151 150 $ concrete vessel out of cylinder 121 8 8.1809e-2 -19 200 312 -313 -150 $ concrete vessel out of cylinder c 122 8 8.1809e-2 -206 313 -314 $ cylinder A 123 6 8.5668e-2 -202 206 313 -314 $ X4 and iron linner of cylinder A 124 999 4.9636e-5 -201 202 313 -314 $ air gap out of cylinder X4 125 6 8.5668e-2 -200 201 313 -314 $ X5 and iron linner of concrete vessel 126 8 8.1809e-2 -19 200 313 -314 151 $ concrete vessel out of cylinder 127 6 8.5668e-2 -19 200 313 -314 -151 150 $ concrete vessel out of cylinder 128 8 8.1809e-2 -19 200 313 -314 -150 $ concrete vessel out of cylinder c 129 8 8.1809e-2 -206 314 -315 $ cylinder A 130 6 8.5668e-2 -202 206 314 -315 $ X4 and iron linner of cylinder A 131 999 4.9636e-5 -201 202 314 -315 $ air gap out of cylinder X4 132 6 8.5668e-2 -200 201 314 -315 $ X5 and iron linner of concrete vessel 133 8 8.1809e-2 -19 200 314 -315 151 $ concrete vessel out of cylinder 134 6 8.5668e-2 -19 200 314 -315 -151 150 $ concrete vessel out of cylinder 135 8 8.1809e-2 -19 200 314 -315 -150 $ concrete vessel out of cylinder c 136 8 8.1809e-2 -206 315 -316 $ cylinder A 137 6 8.5668e-2 -202 206 315 -316 $ X4 and iron linner of cylinder A 138 999 4.9636e-5 -201 202 315 -316 $ air gap out of cylinder X4 139 6 8.5668e-2 -200 201 315 -316 $ X5 and iron linner of concrete vessel 140 8 8.1809e-2 -19 200 315 -316 151 $ concrete vessel out of cylinder 141 6 8.5668e-2 -19 200 315 -316 -151 150 $ concrete vessel out of cylinder 142 8 8.1809e-2 -19 200 315 -316 -150 $ concrete vessel out of cylinder c 143 8 8.1809e-2 -206 316 -317 $ cylinder A 144 6 8.5668e-2 -202 206 316 -317 $ X4 and iron linner of cylinder A 145 999 4.9636e-5 -201 202 316 -317 $ air gap out of cylinder X4 146 6 8.5668e-2 -200 201 316 -317 $ X5 and iron linner of concrete vessel 147 8 8.1809e-2 -19 200 316 -317 151 $ concrete vessel out of cylinder 148 6 8.5668e-2 -19 200 316 -317 -151 150 $ concrete vessel out of cylinder 149 8 8.1809e-2 -19 200 316 -317 -150 $ concrete vessel out of cylinder c 150 8 8.1809e-2 -206 317 -114 $ cylinder A 151 6 8.5668e-2 -202 206 317 -114 $ X4 and iron linner of cylinder A 152 999 4.9636e-5 -201 202 317 -114 $ air gap out of cylinder X4 153 6 8.5668e-2 -200 201 317 -114 $ X5 and iron linner of concrete vessel 154 8 8.1809e-2 -19 200 317 -114 151 $ concrete vessel out of cylinder 155 6 8.5668e-2 -19 200 317 -114 -151 150 $ concrete vessel out of cylinder 156 8 8.1809e-2 -19 200 317 -114 -150 $ concrete vessel out of cylinder c 157 6 8.5668e-2 -202 114 -115 $ iron facing of cylinder A 158 999 4.9636e-5 -201 202 114 -115 $ air gap out of cylinder A 159 6 8.5668e-2 -200 201 114 -115 $ X5 and iron liner of concrete vessel 160 6 8.5668e-2 -18 200 114 $ iron liner facing of concrete vessel c 161 999 4.9636e-5 -20 115 $ air in front of the concrete vessel c 162 999 4.9636e-5 -11 116 $ detector area c 163 6 8.5668e-2 -18 19 -114 $ iron linner of the concrete vessel 164 9 7.8564e-2 -21 20 $ 61 x 61 x 30.5 concrete block(1) 165 9 7.8564e-2 -22 21 $ 61 x 61 x 30.5 concrete block(2) 166 9 7.8564e-2 -11 22 108 -116 $ 61 x 61 x 30.5 concrete block(3) 167 0 11 $ out world 1 so 38.74 $ source surface 2 so 40.65 $ pressure vessel 3 so 46.68 $ water 4 so 48.58 $ source surface 5 cz 49.85 $ air cylinder around the core 6 so 49.85 $ air sphere aroud the core 7 cz 52.39 $ ss-304 cylinder around the core 8 so 52.39 $ ss-304 sphere around the core 9 cz 72.71 $ around the core 65 v/o ss-304 + 35 v/o water 10 so 72.71 $ around the core 65 v/o ss-304 + 35 v/o water c 11 rpp -198.2 198.2 -198.2 198.2 -75.0 350.0 $ body boundary c 12 cz 38.38 $ radius of collimator steel 13 cz 39.02 $ alminum or carbon steel 14 cz 43.42 $ radius of collimator outlet 15 cz 44.06 $ alminum or carbon steel 16 cz 59.0 $ alminum facing c 17 rpp -76.2 76.2 -76.2 76.2 91.36 93.7 $ iron slab in front of slit 18 rpp -76.2 76.2 -76.2 76.2 93.7 193.7 $ iron linner of concrete vessel 19 rpp -74.93 74.93 -74.93 74.93 93.7 192.43 $ iron lined concrete vessel 20 rpp -76.2 76.2 -76.2 76.2 91.36 213.36 $ dummy box around the shield c c caution ! 93.4 - 91.36 air gap is different from problem II - A 21 rpp -86.2 86.2 -86.2 86.2 91.36 213.36 $ 61 x 61 x 30.5 concrete block(1) 22 rpp -137.2 137.2 -137.2 137.2 91.36 213.36 $ 61 x 61 x 30.5 concrete block(2) c 100 pz 0.0 $ core center line 101 pz 52.39 $ gap in front of the collimator 102 pz 53.66 $ collimator inlet 103 pz 73.06 $ collimatot bending ( concrete ) 104 pz 73.70 $ collimator bending ( alminum of carbon steel ) 105 pz 83.86 $ collimator outlet 106 pz 86.22 $ alminum facing ( core side ) 107 pz 91.3 $ alminum facing ( shield side ) 108 pz 91.36 $ beginning of shield configuration 109 pz 93.7 $ end of iron plate 110 pz 98.58 $ cylinder holder 111 pz 142.43 $ end of iron liner of three fourths part 112 pz 143.7 $ beginning of iron liner of a fourth part 113 pz 144.97 $ end of iron liner of quater part 114 pz 192.43 $ iron liner facing ( shield side) 115 pz 193.7 $ iron liner facing ( measurement side) 116 pz 213.36 $ beginning of measurement area c 150 py -1.27 $ horizontal iron plate 151 py 1.27 $ horizontal iron plate c 200 cz 37.67 $ radius of X5 and iron liner of concrete vessel 201 cz 32.382 $ inner radius of X5 202 cz 30.4 $ radius of X4 203 cz 29.67 $ radius of sleeve cylinder 204 cz 28.4 $ inner radius of sleeve cylinder 205 cz 26.4 $ radius of X3 206 cz 24.04 $ inner radius of cylinder A c 300 pz 103.7 $ +10cm penetrate into concrete vessel 301 pz 108.7 $ +15cm 302 pz 113.7 $ +20cm 303 pz 118.7 $ +25cm 304 pz 123.7 $ +30cm 305 pz 128.7 $ +35cm 306 pz 133.7 $ +40cm 307 pz 138.7 $ +45cm c 308 pz 143.7 $ +50cm replaeced by 111 112 309 pz 148.7 $ +55cm 310 pz 153.7 $ +60cm 311 pz 158.7 $ +65cm 312 pz 163.7 $ +70cm 313 pz 168.7 $ +75cm 314 pz 173.7 $ +80cm 315 pz 178.7 $ +85cm 316 pz 183.7 $ +90cm 317 pz 188.7 $ +95cm penetrate into concrete vessel mode n phys:n 20.0 0.0 c source intensity is 2.60547*10**7(cm-2*W-1*min-1) c = 4.34245*10**5(cm-2*W-1*sec-1) on the reactor core m999 1001.60c 1.45005e-4 $ air 1.2049e-3 g*cm-3 1002.60c 2.17540e-5 6000.60c 1.52792e-4 7014.60c 7.84737e-1 7015.60c 2.91431e-3 8016.60c 2.11771e-1 8017.60c 2.57527e-4 $ air JAERI-M6928 m1 13027.60c 9.90222e-1 $ alminum 2.69811 g*cm-3 26054.60c 5.57239e-4 26056.60c 8.98315e-3 26057.60c 2.09355e-4 26058.60c 2.84032e-5 $ alminum m2 1001.60c 6.66467e-1 $ core (water) 1.0 g*cm-3 1002.60c 1.99816e-4 8016.60c 3.32449e-1 8017.60c 8.84598e-4 $ core (water) m3 6000.60c 2.27016e-3 $ ss-304 7.92 g*cm-3 14000.60c 8.72422e-3 24050.60c 8.29848e-3 24052.60c 1.66420e-1 24053.60c 1.92318e-2 24054.60c 4.87799e-2 25055.60c 1.33824e-2 26054.60c 3.92233e-2 26056.60c 6.32312e-1 26057.60c 1.47362e-2 26058.60c 1.99926e-3 28058.60c 5.94858e-2 28060.60c 2.37029e-2 28061.60c 1.04764e-3 28062.60c 3.39425e-3 28064.60c 8.92936e-4 $ ss-304 7.92 g*cm-3 m4 1001.60c 2.52806e-1 $ 65 %vol ss-304 + 35 %vol water 1002.60c 7.57950e-5 6000.60c 1.53697e-3 8016.60c 1.26461e-1 8017.60c 3.36493e-4 14000.60c 5.85594e-3 15031.60c 2.49239e-4 16000.60c 1.88413e-4 24050.60c 4.99768e-3 24052.60c 1.00225e-1 24053.60c 1.15822e-2 24054.60c 2.93772e-3 25055.60c 8.24789e-3 26054.60c 2.38564e-2 26056.60c 3.84584e-1 26057.60c 8.96284e-3 26058.60c 1.21599e-3 27059.60c 4.72717e-4 28058.60c 3.57361e-2 28060.60c 1.42395e-2 28061.60c 6.29371e-4 28062.60c 2.03909e-3 28064.60c 5.36370e-4 29063.60c 7.77891e-3 29065.60c 3.57726e-3 42000.60c 8.71128e-4 $ 65 %vol ss-304 + 35 %vol water m5 1001.60c 1.15028e-1 $ TSF concrete 2.39 g*cm-3 1002.60c 3.44871e-5 6000.60c 1.03271e-1 8016.60c 5.42771e-1 8017.60c 1.44424e-3 11023.60c 3.53740e-4 12000.60c 1.86588e-2 13027.60c 5.36441e-3 14000.60c 4.97569e-2 19000.60c 3.03206e-2 20000.60c 1.29575e-1 26054.60c 1.94944e-4 26056.60c 3.14266e-3 26057.60c 7.32406e-5 26058.60c 9.93657e-6 $ TSF concrete 2.39 g*cm-3 m6 6000.60c 1.15002e-2 $ iron slab 7.86 g*cm-3 14000.60c 4.91782e-3 25055.60c 1.00574e-2 26054.60c 5.54794e-2 26056.60c 8.94374e-1 26057.60c 2.08436e-2 26058.60c 2.82786e-3 $ iron slab 7.86 g*cm-3 m8 1001.60c 1.46640e-1 $ concrete in the iron lined vessel 2.218 g*cm-3 1002.60c 4.39648e-5 6000.60c 1.13020e-1 8016.60c 5.86397e-1 8017.60c 1.56032e-3 11023.60c 2.67698e-4 12000.60c 6.50298e-3 13027.60c 6.84402e-3 14000.60c 1.77732e-2 15031.60c 1.44972e-4 16000.60c 1.34949e-3 19000.60c 1.95701e-3 20000.60c 1.14450e-1 24050.60c 3.60542e-6 24052.60c 7.23042e-5 24053.60c 8.35560e-6 24054.60c 2.11933e-6 26054.60c 1.68857e-4 26056.60c 2.72211e-3 26057.60c 6.34395e-5 26058.60c 8.60686e-6 $ concrete in the iron lined vessel 2.218 g*cm-3 m9 1001.60c 8.14378e-2 $ 61cm x 61cm x 30.5cm concrete block 2.40 g*cm-3 1002.60c 2.44162e-5 6000.60c 1.28558e-1 8016.60c 5.68722e-1 8017.60c 1.51329e-3 12000.60c 2.77481e-2 13027.60c 7.94257e-3 14000.60c 5.53689e-2 19000.60c 1.41286e-3 20000.60c 1.25885e-1 26054.60c 7.90656e-5 26056.60c 1.27460e-3 26057.60c 2.97050e-5 26058.60c 4.03008e-6 $ 61cm x 61cm x 30.5cm concrete block 2.40 g*cm-3 c c f5:n 0.0 0.0 223.7 0.5 f15z:n 223.7 5.0 0.5 f25z:n 223.7 10.0 0.5 f35z:n 223.7 15.0 0.5 f45z:n 223.7 20.0 0.5 c f55z:n 223.7 22.5 0.5 c f65z:n 223.7 23.0 0.5 c f75z:n 223.7 25.0 0.5 c f85z:n 223.7 26.0 0.5 c f95z:n 223.7 26.5 0.5 c f105z:n 223.7 27.0 0.5 c f115z:n 223.7 27.5 0.5 c f125z:n 223.7 28.0 0.5 c f135z:n 223.7 28.5 0.5 c f145z:n 223.7 29.0 0.5 f155z:n 223.7 30.0 0.5 c f165z:n 223.7 30.5 0.5 f175z:n 223.7 31.0 0.5 c f185z:n 223.7 31.5 0.5 f195z:n 223.7 32.0 0.5 c f205z:n 223.7 32.5 0.5 f215z:n 223.7 33.0 0.5 c f225z:n 223.7 33.5 0.5 f235z:n 223.7 34.0 0.5 c f245z:n 223.7 34.5 0.5 f255z:n 223.7 35.0 0.5 f265z:n 223.7 36.0 0.5 f275z:n 223.7 37.0 0.5 f285z:n 223.7 37.5 0.5 f295z:n 223.7 38.0 0.5 f305z:n 223.7 39.0 0.5 f315z:n 223.7 40.0 0.5 f325z:n 223.7 45.0 0.5 f335z:n 223.7 50.0 0.5 f345z:n 223.7 60.0 0.5 f355z:n 223.7 70.0 0.5 e0 1.1254-6 3.9279-6 1.3710-5 4.7851-5 1.6702-4 4.5400-4 1.2341-3 3.3546-3 9.1188-3 2.4788-2 6.7379-2 1.6573-1 2.7324-1 4.5049-1 7.4274-1 1.2246 2.0190 3.3287 5.4881 1.4798+1 cut:n j 4.1399e-7 c nps 1000000 ctme 375 print wwp:n 5 3 5 0 0 wwe:n 1.0e-5 1.0e-3 0.1 15.0 wwn1:n -1.0000E+00 -1.0000E+00 -1.0000E+00 1.3124E+04 1.1485E+03 1.1485E+03 4.8923E+03 1.5226E+27 1.7600E+05 3.7949E+27 1.7600E+05 3.7949E+27 5.5963E+02 4.5966E+02 1.3904E+03 6.8377E+02 1.9946E+03 6.2642E+03 2.9865E+02 2.2380E+02 2.6077E+02 2.3119E+02 2.9865E+02 2.5447E+02 2.0760E+02 1.1735E+02 9.8574E+00 2.7581E+01 3.4478E+02 4.1080E+02 3.2487E+01 6.8981E+00 2.5095E+01 3.5227E+02 5.1986E+02 3.8999E+01 4.3312E+00 2.7324E+01 2.7986E+02 3.6608E+02 2.2904E+01 2.5562E+00 1.4695E+01 1.7102E+02 1.5954E+02 1.1611E+01 1.5312E+00 1.0051E+01 1.1019E+02 3.6726E+01 5.6068E+00 9.0164E-01 4.8408E+00 4.6811E+01 1.1291E+01 1.6814E+00 4.7307E-01 1.5008E+00 1.2472E+01 4.4875E+00 5.5945E-01 2.3266E-01 4.5027E-01 2.2077E+00 1.8226E+00 2.3268E-01 1.1708E-01 1.2826E-01 4.2493E-01 1.1917E+00 1.2080E-01 7.1025E-02 5.2132E-02 1.1538E-01 1.5595E+00 1.1969E-01 6.3911E-02 3.7757E-02 5.6396E-02 8.2403E-02 1.3202E+00 5.3399E-02 1.6477E-02 4.6166E-02 2.0940E+00 4.5594E+00 1.2706E+00 6.0765E-02 1.9627E-02 4.3441E-02 1.7766E+00 1.1550E+00 2.9866E+00 1.1933E+00 6.3874E-02 1.8173E-02 4.8839E-02 2.3344E+00 1.0109E+00 2.5999E+00 1.2028E+00 5.7922E-02 1.3223E-02 4.0008E-02 1.9459E+00 5.0455E+00 2.7842E+00 9.2291E-01 3.8815E-02 9.3072E-03 3.1039E-02 1.6129E+00 2.0282E+00 1.5352E+00 3.2823E-01 2.5705E-02 6.6319E-03 2.2544E-02 1.0623E+00 1.0235E+00 1.2212E+00 2.5458E-01 1.7385E-02 4.3871E-03 1.4749E-02 4.2225E-01 4.1654E-01 4.7873E-01 7.9044E-02 1.0337E-02 3.0071E-03 7.9857E-03 1.2195E-01 2.4637E-01 8.4058E-02 1.3002E-02 4.3913E-03 1.9672E-03 4.1189E-03 2.4536E-02 1.4102E-02 1.8507E-02 2.7609E-03 1.4979E-03 1.1110E-03 1.5535E-03 4.6058E-03 3.3736E-03 4.2852E-03 7.5600E-04 5.9413E-04 8.4550E-04 5.5905E-04 1.2472E-03 9.6682E-04 1.1558E-03 3.8074E-04 2.1385E-03 3.8992E-04 5.8879E-04 5.9743E-03 7.7102E-03 7.0676E+02 9.5125E+01 4.2897E+01 1.7354E+02 -1.0000E+00 wwn2:n -1.0000E+00 -1.0000E+00 -1.0000E+00 7.0466E+02 3.7035E+02 2.0831E+06 8.7721E+02 2.5555E+06 1.3586E+04 6.4689E+06 1.3586E+04 6.4689E+06 1.9961E+02 1.6133E+02 1.3543E+02 1.6990E+02 4.3692E+02 5.0620E+02 9.2537E+01 6.1325E+01 7.2364E+01 6.4728E+01 2.0310E+05 7.1163E+01 5.1939E+01 2.1123E+01 3.4093E+00 9.0103E+00 1.0638E+02 4.6253E+01 9.4058E+00 2.3088E+00 7.2799E+00 9.6614E+01 3.3222E+01 7.7123E+00 1.4987E+00 7.6939E+00 7.1019E+01 2.3167E+01 5.0678E+00 8.9178E-01 4.1667E+00 4.0486E+01 1.2210E+01 2.5850E+00 5.0241E-01 2.1365E+00 1.5361E+01 4.6446E+00 1.1041E+00 2.6951E-01 7.9013E-01 4.8363E+00 1.5879E+00 3.7226E-01 1.4035E-01 3.0791E-01 1.4863E+00 6.9176E-01 1.4157E-01 7.4167E-02 1.1910E-01 4.3602E-01 3.3653E-01 6.1076E-02 3.9126E-02 4.4698E-02 1.3158E-01 2.1955E-01 3.8092E-02 2.5373E-02 2.0607E-02 4.8750E-02 2.1851E-01 3.6328E-02 2.3590E-02 1.6467E-02 2.7556E-02 3.8264E-02 2.2667E-01 2.2317E-02 9.3247E-03 2.0272E-02 2.8345E-01 4.8703E-01 1.8729E-01 2.2472E-02 1.0657E-02 1.9396E-02 3.0628E-01 1.2033E-01 4.7252E-01 1.7989E-01 2.1906E-02 9.4092E-03 1.8785E-02 3.0026E-01 3.1717E-01 4.2092E-01 1.6178E-01 1.8117E-02 6.9509E-03 1.5763E-02 3.0545E-01 2.1946E-01 3.4604E-01 1.1543E-01 1.4036E-02 5.1109E-03 1.2867E-02 1.9851E-01 6.5475E-02 2.2510E-01 5.7336E-02 8.7992E-03 3.4370E-03 8.4262E-03 1.0437E-01 4.8437E-02 1.0344E-01 2.2974E-02 5.6668E-03 2.2993E-03 5.1335E-03 4.3067E-02 1.7244E-02 3.6294E-02 8.5945E-03 3.0052E-03 1.5108E-03 2.8611E-03 1.4505E-02 9.1821E-03 1.2681E-02 2.7400E-03 1.4368E-03 9.4095E-04 1.4912E-03 4.7423E-03 4.4202E-03 4.4220E-03 8.1037E-04 6.2065E-04 5.7764E-04 6.6140E-04 1.4695E-03 1.4400E-03 1.3783E-03 3.5707E-04 3.2378E-04 4.7522E-04 3.3031E-04 5.7724E-04 4.4283E-04 5.9016E-04 2.8375E-04 9.5267E-04 2.5303E-04 4.2792E-04 3.9080E-03 4.9014E-03 4.0091E+01 7.4418E+00 3.6023E+01 3.6023E+01 -1.0000E+00 wwn3:n -1.0000E+00 -1.0000E+00 -1.0000E+00 7.7462E+01 4.5347E+01 1.0618E+03 1.4678E+02 1.2397E+03 9.1267E+02 4.1320E+04 9.1267E+02 4.1320E+04 2.2774E+01 2.6738E+01 6.0082E+01 3.2615E+01 7.9749E+01 1.2582E+02 1.2750E+01 1.0037E+01 1.1395E+01 1.0411E+01 1.7299E+03 1.1194E+01 9.6057E+00 4.2469E+00 1.3101E+00 1.8115E+00 1.5529E+01 5.3030E+00 1.2840E+00 8.4158E-01 1.3600E+00 1.2240E+01 4.7502E+00 7.4167E-01 5.1304E-01 9.4260E-01 7.1365E+00 2.4758E+00 4.2662E-01 2.8405E-01 6.7419E-01 4.3345E+00 1.3529E+00 2.6387E-01 1.5899E-01 3.3827E-01 2.0149E+00 6.2272E-01 1.3378E-01 8.1315E-02 1.7129E-01 9.2943E-01 3.1556E-01 7.9466E-02 4.6006E-02 9.6350E-02 3.8410E-01 1.7123E-01 4.9402E-02 2.6747E-02 4.5047E-02 1.4413E-01 9.5661E-02 2.8439E-02 1.7028E-02 1.9936E-02 5.6638E-02 7.1720E-02 1.8687E-02 1.2470E-02 1.1235E-02 2.5421E-02 6.7810E-02 1.6334E-02 1.1101E-02 8.7187E-03 1.8495E-02 1.9363E-02 6.4196E-02 9.5544E-03 7.0734E-03 1.1660E-02 9.8752E-02 1.4606E-01 4.8861E-02 8.9379E-03 7.0135E-03 1.0378E-02 8.3084E-02 6.6169E-02 9.6361E-02 3.6941E-02 7.3194E-03 5.1514E-03 8.5011E-03 6.5265E-02 4.5255E-02 6.8382E-02 2.5709E-02 5.4661E-03 3.6444E-03 6.1080E-03 4.2830E-02 4.1906E-02 4.8818E-02 1.7336E-02 3.6577E-03 2.4956E-03 4.0600E-03 2.5536E-02 3.0823E-02 2.7187E-02 9.4969E-03 2.4705E-03 1.7498E-03 2.5879E-03 1.3337E-02 9.9897E-03 1.4705E-02 4.4503E-03 1.4819E-03 1.1268E-03 1.6001E-03 6.9906E-03 6.0744E-03 7.7490E-03 2.1065E-03 9.3531E-04 7.6997E-04 1.0159E-03 3.8256E-03 3.4875E-03 3.5629E-03 9.7079E-04 5.9158E-04 5.5182E-04 6.3840E-04 1.6511E-03 1.1663E-03 1.5030E-03 4.5539E-04 4.0283E-04 4.3372E-04 4.3144E-04 7.2181E-04 7.3031E-04 7.1938E-04 2.7963E-04 3.3226E-04 4.7059E-04 3.6183E-04 4.2950E-04 4.0476E-04 4.1450E-04 4.3793E-04 9.4988E-04 5.7987E-04 4.6087E-04 3.0416E-03 5.8690E-03 8.4958E-01 1.1392E+00 6.7444E+00 6.7444E+00 -1.0000E+00 wwn4:n -1.0000E+00 -1.0000E+00 -1.0000E+00 5.2657E+00 2.7832E+00 4.8733E+02 1.8664E+01 5.2639E+02 4.2942E+01 1.2031E+03 9.0342E+03 2.1594E+07 1.0668E+00 2.4833E+00 2.8043E+00 1.7324E+00 5.3771E+00 4.7660E+00 6.8538E-01 5.1478E-01 5.0703E-01 4.9506E-01 1.6351E+01 5.0204E-01 5.0000E-01 3.8514E-01 1.4978E-01 2.9336E-01 7.9115E-01 3.9701E-01 2.0586E-01 9.9289E-02 2.0908E-01 4.9447E-01 2.4585E-01 1.2827E-01 6.0161E-02 1.3877E-01 2.9645E-01 1.4901E-01 7.8261E-02 3.6342E-02 9.0115E-02 1.8440E-01 9.0036E-02 4.9165E-02 2.3271E-02 5.9153E-02 1.1536E-01 5.6923E-02 3.0266E-02 1.5038E-02 3.7204E-02 7.3088E-02 3.6117E-02 1.9155E-02 9.8343E-03 2.2451E-02 4.4333E-02 2.4221E-02 1.3364E-02 6.8816E-03 1.2574E-02 2.6810E-02 1.6327E-02 9.1817E-03 5.2387E-03 7.1369E-03 1.5659E-02 1.1808E-02 6.8445E-03 4.3383E-03 4.8090E-03 9.3365E-03 9.4736E-03 5.5916E-03 3.8454E-03 3.9524E-03 7.8714E-03 6.7611E-03 8.6414E-03 3.4340E-03 3.5707E-03 5.7603E-03 1.4777E-02 1.1612E-02 7.9043E-03 3.2573E-03 3.4771E-03 5.0866E-03 1.3447E-02 1.4819E-02 1.0428E-02 6.0356E-03 2.5821E-03 2.5465E-03 3.6691E-03 9.8974E-03 1.0646E-02 8.0520E-03 4.1394E-03 1.8979E-03 1.7808E-03 2.4896E-03 6.6833E-03 4.5939E-03 5.6544E-03 2.8116E-03 1.4134E-03 1.2477E-03 1.6794E-03 4.5323E-03 2.9433E-03 3.7768E-03 1.8689E-03 1.0594E-03 8.5866E-04 1.1474E-03 2.9825E-03 2.7896E-03 2.4814E-03 1.2517E-03 7.8570E-04 6.6787E-04 8.3962E-04 1.9029E-03 1.8325E-03 1.6406E-03 8.0277E-04 5.8052E-04 5.5956E-04 5.9272E-04 1.2411E-03 1.1258E-03 1.0716E-03 5.1969E-04 4.5656E-04 4.7692E-04 4.5662E-04 7.8625E-04 7.1362E-04 6.9503E-04 3.4578E-04 3.9877E-04 4.5192E-04 3.9353E-04 5.2175E-04 5.9413E-04 4.8413E-04 2.7765E-04 3.8739E-04 5.7491E-04 4.5799E-04 4.3611E-04 4.8251E-04 4.4465E-04 8.4733E-04 2.1702E-03 9.7644E-04 9.5118E-04 5.2437E-03 1.5538E-02 3.5100E-01 3.8127E-01 4.9713E-01 7.4084E+00 -1.0000E+00 1