11 real*8 Z, Y1, Y2, B, s, rts, s0, eta1, eta2
12 parameter(z= 35.45, y1=42.53, y2=33.34, b=0.308,
13 * s0=5.38**2, eta1=0.458, eta2=0.545)
15 data ( px(i), i= 1 , np )/
16 1 9.7289
e-02, 0.1485 , 0.1925 , 0.2330 ,
17 2 0.2821 , 0.3558 , 0.4426 , 0.5736 ,
18 3 0.6850 , 0.8180 , 1.004 , 1.319 ,
19 4 1.686 , 2.098 , 2.911 , 4.207 ,
20 5 5.918 , 9.286 , 13.79 , 21.64 ,
21 6 32.15 , 49.08 , 70.95 , 109.8 ,
22 7 144.3 , 192.2 , 274.0 , 418.3 ,
23 8 702.7 , 1133. , 1877. , 2986. ,
26 data ( mb(i), i= 1 , np )/
27 1 445.7 , 380.5 , 321.4 , 283.2 ,
28 2 249.5 , 214.1 , 179.9 , 153.6 ,
29 3 134.6 , 123.0 , 113.7 , 103.9 ,
30 4 95.54 , 86.43 , 75.75 , 66.40 ,
31 5 59.76 , 54.64 , 51.30 , 47.16 ,
32 6 45.22 , 43.82 , 42.92 , 42.03 ,
33 7 41.82 , 41.39 , 41.84 , 42.07 ,
34 8 42.31 , 43.22 , 43.70 , 45.12 ,
37 if(p .gt. 15000.)
then 42 xs =z + b*log(s/s0)**2 + y1*(1./s)**eta1 +y2*(1./s)**eta2
43 elseif( p .gt. 0.1)
then 62 data ( px(i), i= 1 , np )/
63 1 9.8633
e-02, 0.1485 , 0.2299 , 0.3062 ,
64 2 0.4078 , 0.5358 , 0.6850 , 0.8758 ,
65 3 1.166 , 1.575 , 1.986 , 2.505 ,
66 4 3.572 , 4.957 , 6.784 , 9.286 ,
67 5 11.87 , 15.18 , 21.64 , 32.60 ,
68 6 56.27 , 83.59 , 127.6 , 289.4 ,
69 7 486.1 , 936.0 , 1551. , 3241. ,
70 8 6680. , 1.4542
e+04, 3.3434
e+04, 6.4372
e+04,
71 9 1.5845
e+05, 6.2031
e+05, 1.7030
e+06
73 data ( mb(i), i= 1 , np )/
74 1 106.0 , 96.90 , 86.76 , 78.90 ,
75 2 68.79 , 60.62 , 53.69 , 47.81 ,
76 3 43.25 , 37.51 , 32.19 , 26.90 ,
77 4 20.89 , 16.91 , 13.91 , 11.81 ,
78 5 10.80 , 9.772 , 9.078 , 8.568 ,
79 6 7.960 , 7.513 , 7.358 , 7.170 ,
80 7 7.135 , 7.252 , 7.488 , 8.239 ,
81 8 9.064 , 10.13 , 11.21 , 12.07 ,
82 9 13.43 , 15.42 , 17.42
85 if( p .gt. 2000.)
then 87 if(xssave .lt. 0.)
then 91 xs = xs * mb(27)/xssave
92 elseif(p .gt. 0.1)
then 101 real(8),
intent(in)::p
102 real(8),
intent(out)::xs
107 xs =max( txs - exs, 0.
d0)
integer npitbl real *nx parameter(n=101, npitbl=46, nx=n-1) real *8 uconst
dE dx *! Nuc Int sampling table e
subroutine cpbarptotxs(p, xs)
subroutine kpolintplogxyfe(xa, xstep, ya, ystep, nt, m, logxy, x, y, error)
block data cblkEvhnp ! currently usable models data RegMdls ad *special data *Cekaon d0
subroutine cpbarpelaxs(p, xs)
subroutine cpbarpinelaxs(p, xs)