53 SUBROUTINE parpol (ISET, X, Q2, UV, DV, QB, ST, GL)
59 COMMON /pepadm/cpdfnam(2,imxpdf),ipdfnam(2,imxpdf),
60 & iplst(10),cunpol,cpol
61 CHARACTER*256 cpdfnam,cunpol,cpol
74 real*8 xb1,
dfint,xb0,q2,uv,
x,st,gl,dv,qb,
xb,qs,
75 &
parton,arrf,xt,xqbf,xuvf,xdvf,xgf,xsf
76 integer n,iq,m,ix,iini,iset,
npart,nx,nq,narg
80 dimension xuvf(nx,nq), xdvf(nx,nq), xqbf(nx,nq),
82 2 qs(nq),
xb(nx), xt(narg), na(narg), arrf(nx+nq)
83 COMMON / intini / iini
84 SAVE xuvf, xdvf, xqbf, xsf, xgf, na, arrf
86 DATA qs / 0.4d0, 0.75d0, 1.0d0, 1.5d0, 2.5d0,
87 1 4.0d0, 6.4d0, 1.0d1, 1.5d1, 2.5d1, 4.0d1, 6.4d1,
88 2 1.0d2, 1.8d2, 3.2d2, 5.8d2, 1.0d3, 1.8d3,
89 3 3.2d3, 5.8d3, 1.0d4 /
90 DATA xb / 1.
d-4, 1.5
d-4, 2.2
d-4, 3.2
d-4, 4.8
d-4, 7.
d-4,
91 2 1.
d-3, 1.5
d-3, 2.2
d-3, 3.2
d-3, 4.8
d-3, 7.
d-3,
92 3 1.
d-2, 1.5
d-2, 2.2
d-2, 3.2
d-2, 5.0
d-2, 7.5
d-2,
93 4 0.1, 0.125, 0.15, 0.175, 0.2, 0.225, 0.25, 0.275,
94 5 0.3, 0.325, 0.35, 0.375, 0.4, 0.45, 0.5, 0.55,
95 6 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 1.0 /
97 IF ( (
x.LT.1.0
d-4) .OR. (
x.GT.1.0d0) )
THEN
99 91
FORMAT (2
x,
'PARTON INTERPOLATION: X OUT OF RANGE')
102 IF ( (q2.LT.0.4d0) .OR. (q2.GT.1.d4) )
THEN
104 92
FORMAT (2
x,
'PARTON INTERPOLATION: Q2 OUT OF RANGE')
118 IF (iini.NE.0) goto 16
119 OPEN(iplst(3),file=cpol,
status=
'OLD')
123 READ(iplst(3),90)
parton(1,
n,m),
parton(2,
n,m),
parton(3,
n,m),
127 90
FORMAT (5(1pe10.3))
135 xuvf(ix,iq) =
parton(1,iq,ix) / (xb1**3 * xb0**0.7)
136 xdvf(ix,iq) =
parton(2,iq,ix) / (xb1**4 * xb0**0.6)
137 xqbf(ix,iq) =
parton(3,iq,ix) / (xb1**7 * xb0**0.3)
138 xsf(ix,iq) =
parton(4,iq,ix) / (xb1**7 * xb0**0.3)
139 xgf(ix,iq) =
parton(5,iq,ix) / (xb1**10 * xb0**0.3)
150 arrf(ix) = dlog(
xb(ix))
153 arrf(nx+iq) = dlog(qs(iq))
159 uv =
dfint(narg,xt,na,arrf,xuvf) * (1.d0-
x)**3 *
x**0.7
160 dv =
dfint(narg,xt,na,arrf,xdvf) * (1.d0-
x)**4 *
x**0.6
161 qb =
dfint(narg,xt,na,arrf,xqbf) * (1.d0-
x)**7 *
x**0.3
162 st =
dfint(narg,xt,na,arrf,xsf) * (1.d0-
x)**7 *
x**0.3
163 gl =
dfint(narg,xt,na,arrf,xgf) * (1.d0-
x)**10 *
x**0.3