00001 *!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
00002 * THIS VERSION IS CHANGED FOR TESTS:
00003 * INCLUDES LEPTONIC PART + HADRONIC PART WITH NARROW RESONANCES SUBTRACTED
00004 * NARROW RESONANCES == PHI, J/PSI, Ypsilon
00005 *
00006 *;;;;;;;;;;;;;;;;;;;;;;;;;;;; -*- Mode: Fortran -*- ;;;;;;;;;;;;;;;;;;;;;;;;;;;
00007 *; alphaQED.f ---
00008 *; Author : Fred Jegerlehner
00009 *; Created On : Sat Nov 1 22:49:46 2003
00010 *; Last Modified By: Fred Jegerlehner
00011 *; Last Modified On: Mon Nov 3 02:23:59 2003
00012 *; RCS: $Id: vac_pol_hc1.inc,v 1.1.1.1 2007/11/18 09:51:28 azhemchugov Exp $
00013 *;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
00014 *; Copyright (c) 2003 Fred Jegerlehner
00015 *;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
00016 *;
00017 c PROGRAM alphaQEDn
00018 C Calculating the running QED coupling alpha(E)
00019 C Leptons: 1--loop, 2--loop exact, 3--loop in high energy approximation
00020 C Quarks form routine dhadr5n including effective s--channel shifts in low energy range
00021 C r1....,r2....,r3.... : oneloop, twoloop, threeloop result
00022 C ..real: realpart, ..imag: imaginary part;
00023 C ......l: light fermions only (high energy approximation for 3--loops)
00024 C dallepQED1n=r1real,dallepQED2n=r2real,dallepQED3l=r3reall
00025 C dallepQED3l also returns oneloop and twoloop high energy (light fermion) approximations
00026 c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
00027 c
00028 c H.Czyz: - the constants() are called in subroutine 'input' of PHOKHARA
00029 c where also the vacpol at 's' is calculated for FSR contributions;
00030 c different scales (Q^2, Q_a^2, Q_b^2) contributions are calculated
00031 c in proper places of PHOKHARA. The function dggvap is now complex*16
00032 c and contains also the imaginary part of the vacpol and is used
00033 c as additional factor 1.d0/(1.d0-dggvap(s, 0.d0)) in the amplitudes.
00034 c Few commons were added to 'common.f', previously written separately
00035 c
00036 c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
00037 c
00038 c call constants()
00039 C LEPTONflag=all,had,lep,ele,muo,tau -> iLEP=-3,-2,-1,1,2,3
00040 C Default:
00041 C LEPTONflag='all'
00042 C iLEP = -3 ! for sum of leptons + quarks
00043 C Decomment the folowwing two lines for changing the default
00044 c LEPTONflag='tau'
00045 c iLEP=LFLAG(LEPTONflag)
00046 ************************************************************************
00047 c
00048 c **********************************************************************
00049 c **********************************************************************
00050
00051 subroutine constants()
00052 implicit none
00053 include 'Phokhara/common.inc'
00054 data pi,pi2,ln2/3.141592653589793d0,9.869604401089358d0,
00055 & 0.6931471805599453d0/
00056 data zeta2,zeta3/1.644934066848226d0,1.20205690315959d0/
00057 data zeta5/1.036927755143370d0/
00058 DATA ALP/0.0072973544D0/
00059 data ml /0.51099906d-3,105.658389d-3,1.7771d0/
00060 data small,large,large_3/1.d-4,1.d6,8.d0/
00061 save
00062 adp = alp/pi
00063 adp2 = adp*adp
00064 adp3 = adp*adp2
00065 C switch off/on warnings iwarnings=0/1 default is 1
00066 c H.C. - warnings set to null !!!!!!!!!!!!!!!!!!!!!!!!!!!
00067 iwarnings=0
00068 iwarnings=1
00069 C LEPTONflag=all,had,lep,ele,muo,tau -> iLEP=-3,-2,-1,1,2,3
00070 LEPTONflag='all'
00071 iLEP = -3 ! for sum of leptons + quarks
00072 C Overwrite default switches for replacing exact by asymptotic expansions
00073 C (required to enforce numerical stability)
00074 C Do not change without checking whethe new switches give acceptable results
00075 C z=|s|/m^2 (threshold= 4 m^2),
00076 C low energy expansion up to O(z^3): c_1*z+c_2*z^2
00077 C high energy expansion up to O(z^-3): c_1/z+c_2/z^2 up to logs
00078 C
00079 * small=1d-6 ! default
00080 * large=1d6 ! default
00081 C
00082 C corresponding switch for 3--loop high energy approximation
00083 C where exact result is not available at lower energies
00084 C i.e. below large_3*|s| 3--loop contribution is taken to be zero
00085 C
00086 * large_3=8.d0 ! default
00087 C
00088
00089 st2=0.23d0 ! Reference value for weak mixing parameter
00090 als=0.118d0 ! alpha strong
00091 mtop=175.d0 ! top quark mass
00092
00093 return
00094 end
00095
00096 c ***************************************************************************
00097 c ***************************************************************************
00098
00099 FUNCTION dggvap(sp,erfl)
00100 C -----------------------
00101 C Fermionic contributions to the running of alpha
00102 C Function collecting 1-loop, 2-loop, 3-loop leptonic contributions (e, mu, tau)
00103 C + hadronic contribution from light hadrons (udscb)
00104 C + top quark contribution (t); here at 1-loop (good up to 200 GeV LEP I/II)
00105 C erfl is the errorflag: 0.d0->central value +/- 1.d0 -> upper/lower bound
00106 implicit none
00107 complex*16 dggvap
00108 real*8 sp,erfl
00109 real*8 c0
00110 real*8 mtop2
00111 real*8 dallepQED1,dallepQED2,dallepQED3l,
00112 & dalQED1ferm,res_re,res_im
00113 real*8 r1real,r1imag,r2real,r2imag,imag
00114 real*8 r1reall,r1imagl,r2reall,r2imagl,r3reall,r3imagl
00115 real*8 fac,ee,e,der,errder,deg,errdeg,der_help
00116 real*8 hfun
00117 include 'Phokhara/common.inc'
00118 e=dsqrt(dabs(sp))
00119 res_im=0.d0
00120 if (sp.lt.0.d0) e=-e
00121 if ((LEPTONflag.eq.'all').or.(LEPTONflag.eq.'had')) then
00122 c !!!!!!!!!!!
00123 if((e.gt.0.9d0).and.(e.lt.1.1d0))then
00124 call dhadr5n(0.9d0,st2,der,errder,deg,errdeg)
00125 call dhadr5n(1.1d0,st2,der_help,errder,deg,errdeg)
00126 der = der+(e-0.9d0)/(1.1d0-0.9d0)*(der_help-der)
00127 elseif((e.gt.2.4d0).and.(e.lt.3.5d0))then
00128 call dhadr5n(2.4d0,st2,der,errder,deg,errdeg)
00129 call dhadr5n(3.5d0,st2,der_help,errder,deg,errdeg)
00130 der = der+(e-2.4d0)/(3.5d0-2.4d0)*(der_help-der)
00131 elseif((e.gt.3.5d0).and.(e.lt.4.d0))then
00132 call dhadr5n(3.5d0,st2,der,errder,deg,errdeg)
00133 call dhadr5n(4.0d0,st2,der_help,errder,deg,errdeg)
00134 der = der+(e-3.5d0)/(4.d0-3.5d0)*(der_help-der)
00135 else
00136 call dhadr5n(e,st2,der,errder,deg,errdeg)
00137 endif
00138 c !!!!!!!!!!!
00139 dahadr=der
00140 else
00141 dahadr=0.d0
00142 endif
00143 if ((LEPTONflag.eq.'had').or.(LEPTONflag.eq.'top')) then
00144 dalept=0.d0
00145 else
00146 r1real =dallepQED1(sp,r1imag)
00147 r2real =dallepQED2(sp,r2imag)
00148 c !!!!!!!!!!!!!!!
00149 c the third loop switched of as it is not reliable in low energy region
00150 c r3reall=dallepQED3l(sp,r1reall,
00151 c & r1imagl,r2reall,r2imagl,r3imagl)
00152 r3reall = 0.d0
00153 r3imagl = 0.d0
00154 c!!!!!!!!!!!!!!!
00155 res_re=r1real+r2real+r3reall
00156 res_im=r1imag+r2imag+r3imagl
00157 dalept=res_re
00158 endif
00159 if ((LEPTONflag.eq.'all').or.(LEPTONflag.eq.'top')) then
00160 mtop2=mtop*mtop
00161 daltop=dalQED1ferm(sp,mtop2,imag)*4.d0/3.d0
00162 else
00163 daltop=0.d0
00164 endif
00165 dggvap=dcmplx(dalept+dahadr+erfl*errder,res_im)
00166 RETURN
00167 END
00168 c ***********************************************************************
00169 c ***********************************************************************
00170
00171 function dallepQED1(sp,r1imag)
00172 implicit none
00173 integer i
00174 real *8 dallepQED1,sp,m2,r1real,r1imag,x0,y0,dalQED1ferm,imag
00175 include 'Phokhara/common.inc'
00176 x0 = 0.d0
00177 y0 = 0.d0
00178 if ((LEPTONflag.eq.'had').or.(LEPTONflag.eq.'top')) then
00179 dallepQED1=0.d0
00180 r1imag=0.d0
00181 return
00182 endif
00183 if ((LEPTONflag.eq.'all').or.(LEPTONflag.eq.'lep')) then
00184 do i=1,3
00185 m2=ml(i)**2
00186 x0=x0+dalQED1ferm(sp,m2,imag)
00187 y0=y0+imag
00188 enddo
00189 else
00190 iLEP=LFLAG(LEPTONflag)
00191 m2=ml(iLEP)**2
00192 x0=x0+dalQED1ferm(sp,m2,imag)
00193 y0=y0+imag
00194 endif
00195 r1real=x0
00196 r1imag=y0
00197 dallepQED1=r1real
00198 return
00199 end
00200 C
00201 function dalQED1ferm(sp,m2,imag)
00202 implicit none
00203 integer ini
00204 real *8 dalQED1ferm,dalQED1,sp,m2,imag,x
00205 real *8 s0,m2ds,m2ds2,lnm2ds,z,r,y,y2,y3,epy,emy,
00206 & epysq,emysq,emy2,emy3,emy4,lny,abr1,ddilog,dli3,
00207 & lnm2dsimag,lnyimag,tau,sint,sint2,sint3,cost
00208 real *8 null,one,two,four,eight,half,fourth,fourthird,
00209 & twentyninth,sixteenthird,c1,c2
00210 include 'Phokhara/common.inc'
00211 save
00212 data ini/0/
00213 if (ini.eq.0) goto 2
00214 1 continue
00215 C switches for z=sp/m2 input via Phokhara/common.inc from constants.f
00216 c small=eps
00217 c large=one/eps
00218 s0=four*m2 ! threshold
00219 z=sp/m2
00220 x=s0/sp
00221 imag=null
00222 C LOW ENERGY ASYMPTOTE
00223 if (abs(z).lt.small) then
00224 c write (*,*) ' LOW ENERGY ASYMPTOTE'
00225 dalQED1 = c1*z+c2*z*z
00226 goto 9
00227 endif
00228 C TIME-LIKE BRANCH
00229 if (z.gt.null) then
00230 C HIGH ENERGY ASYMPTOTE (TIME-LIKE)
00231 if (z.gt.large) then
00232 c write (*,*) ' HIGH ENERGY ASYMPTOTE (TIME-LIKE)'
00233 m2ds=m2/sp
00234 lnm2ds=log(m2ds)
00235 lnm2dsimag=null
00236 if (sp.gt.s0) lnm2dsimag=pi
00237 dalQED1 = -fourthird*lnm2ds-twentyninth
00238 & +eight*m2ds*(one+m2ds*(lnm2ds-half))
00239 imag = lnm2dsimag*(-fourthird+eight*m2ds*m2ds)
00240 else if (sp.gt.s0) then
00241 r=sqrt(one-x)
00242 y=(r-one)/(r+one)
00243 C for analytic continuation y = y + i epsilon
00244 y2=y*y
00245 epy=one+y
00246 emy=one-y
00247 epysq=one+y2
00248 emy2=emy*emy
00249 emy3=emy2*emy
00250 lny=log(abs(y))
00251 abr1=epysq-four*y
00252 lnyimag=null
00253 if (sp.gt.s0) lnyimag=pi
00254 dalQED1 = -twentyninth+sixteenthird*y/emy2
00255 & -fourthird*epy*abr1/emy3*lny
00256 imag = lnyimag*(-fourthird*epy*abr1/emy3)
00257 else
00258 sint2=sp/s0
00259 sint =dsqrt(sint2)
00260 sint3=sint2*sint
00261 tau=dasin(sint)
00262 cost=cos(tau)
00263 dalQED1 = -twentyninth + fourthird*
00264 & (tau*cost*(one + two*sint2) - sint)/sint3
00265 endif
00266 C SPACE-LIKE BRANCH
00267 else
00268 if (-z.gt.large) then
00269 C HIGH ENERGY ASYMPTOTE (SPACE-LIKE)
00270 c write (*,*) ' HIGH ENERGY ASYMPTOTE (SPACE-LIKE)'
00271 m2ds=-m2/sp
00272 lnm2ds=log(m2ds)
00273 dalQED1 = -fourthird*lnm2ds-twentyninth
00274 & +eight*m2ds*(one+m2ds*(lnm2ds-half))
00275 imag = null
00276 else
00277 r=sqrt(one-x)
00278 y=(r-one)/(r+one)
00279 C for analytic continuation y = y + i epsilon
00280 y2=y*y
00281 epy=one+y
00282 emy=one-y
00283 epysq=one+y2
00284 emy2=emy*emy
00285 emy3=emy2*emy
00286 lny=log(y)
00287 abr1=epysq-four*y
00288 lnyimag=null
00289 if (sp.gt.s0) lnyimag=pi
00290 dalQED1 = -twentyninth+sixteenthird*y/emy2
00291 & -fourthird*epy*abr1/emy3*lny
00292 imag = null
00293 endif
00294 endif
00295 9 dalQED1ferm=dalQED1*adp*fourth
00296 imag=imag*adp*fourth
00297 RETURN
00298 2 null=0.d0
00299 one=1.d0
00300 two=2.d0
00301 four=4.d0
00302 eight=8.d0
00303 half=0.5d0
00304 fourth=0.25d0
00305 fourthird=4.d0/3.d0
00306 twentyninth=20.d0/9.d0
00307 sixteenthird=16.d0/3.d0
00308 c1=-4.d0/15.d0
00309 c2=-1.d0/35.d0
00310 ini=1
00311 c write (*,*) ' numerical constants initialized'
00312 goto 1
00313 END
00314
00315
00316 integer function LFLAG(Label)
00317 implicit none
00318 character*3 Label
00319 integer i
00320 i=-3
00321 if (Label.eq.'had') i=-2
00322 if (Label.eq.'lep') i=-1
00323 if (Label.eq.'ele') i= 1
00324 if (Label.eq.'muo') i= 2
00325 if (Label.eq.'tau') i= 3
00326 LFLAG=i
00327 return
00328 end
00329
00330 function dallepQED2(sp,r2imag)
00331 C Källen&Saby result 2--loop QEC = 2--loop QCD divided by NC=3 and CF=4/3
00332 implicit none
00333 integer i,nl
00334 real *8 dallepqed2,sp,x1,y1,x1l,r2real,r2imag,
00335 & dalQED2ferm,m2,imag,M1
00336 real *8 null,one,two,half,four,qua
00337 include 'Phokhara/common.inc'
00338 data null,one,two,half,four,qua/0.d0,1.d0,2.d0,.5d0,4.d0,.25d0/
00339 x1 = null
00340 y1 = null
00341 if ((LEPTONflag.eq.'had').or.(LEPTONflag.eq.'top')) then
00342 dallepQED2=0.d0
00343 r2imag=0.d0
00344 return
00345 endif
00346 nl=1
00347 iLEP=LFLAG(LEPTONflag)
00348 if ((LEPTONflag.eq.'all').or.(LEPTONflag.eq.'lep')) nl=3
00349 do i=1,nl
00350 M1=ml(i)
00351 if (nl.eq.1) M1=ml(iLEP)
00352 C
00353 m2=M1**2
00354 x1=x1+dalQED2ferm(sp,m2,imag)
00355 y1=y1+imag
00356 enddo
00357 r2real=x1
00358 r2imag=y1
00359 dallepqed2 = r2real
00360 return
00361 end
00362 C
00363 function dalQED2ferm(sp,m2,imag)
00364 C Exact 2-loop result, originally by Kallen and Sabry 1955,
00365 C here based on a compact form worked out by M.Yu. Kalmykov 2003
00366 implicit none
00367 integer ini
00368 real *8 dalQED2ferm,dalQED2,sp,m2,imag,x
00369 real *8 s0,m2ds,m2ds2,lnm2ds,z,omx,r,y,y2,y3,y4,epy,emy,
00370 & epysq,emysq,emy2,emy3,emy4,lny,abr1,ddilog,dli3,
00371 & lnm2dsimag,lnyimag,sint,sint2,sint3,sint4,sint4i,
00372 & tau,cost,phi,chi,fac1,term2,Clausen2,Clausen3
00373 real *8 null,one,two,three,four,five,six,seven,eight,fourty,
00374 & twentytwo,thirtytwo,fourtyeight,third,onesixteenth,
00375 & tenthird,hundertfourthird,sixteenzeta3,eightthird,
00376 & sixteenthird,thirtytwothird,sixtyfourzeta3,c1,c2,
00377 & fourth,sixteen,fourteenthird,twentysixthird
00378 include 'Phokhara/common.inc'
00379 save
00380 data ini/0/
00381 if (ini.eq.0) goto 2
00382 1 continue
00383 C switches for z=sp/m2 input via Phokhara/common.inc from constants.f
00384 c small=eps
00385 c large=one/eps
00386 s0=four*m2 ! threshold
00387 z=sp/m2
00388 x=s0/sp
00389 omx=one-x
00390 imag=null
00391 lnm2dsimag=null
00392 C LOW ENERGY ASYMPTOTE
00393 if (abs(z).lt.small) then
00394 c write (*,*) ' LOW ENERGY ASYMPTOTE'
00395 dalQED2 = c1*z+c2*z*z
00396 c write (*,*) sp,dalQED2*adp2*onesixteenth
00397 goto 9
00398 endif
00399 C TIME-LIKE BRANCH
00400 if (z.gt.null) then
00401 C HIGH ENERGY ASYMPTOTE (TIME-LIKE)
00402 if (z.gt.large) then
00403 c write (*,*) ' HIGH ENERGY ASYMPTOTE (TIME-LIKE)'
00404 m2ds=m2/sp
00405 m2ds2=m2ds*m2ds
00406 lnm2ds=log(m2ds)
00407 lnm2dsimag=pi
00408 dalQED2 = -four*lnm2ds-tenthird+sixteenzeta3
00409 & +fourtyeight*m2ds*lnm2ds
00410 & -m2ds2*(eightthird+sixtyfourzeta3+fourty*lnm2ds
00411 & -fourtyeight*(lnm2ds**2-lnm2dsimag**2))
00412 imag = lnm2dsimag*(-four+fourtyeight*m2ds
00413 & -m2ds2*(fourty-fourtyeight*two*lnm2ds))
00414 else if (sp.gt.s0) then
00415 r=sqrt(omx)
00416 y=(r-one)/(r+one)
00417 C for analytic continuation y = y + i epsilon
00418 y2=y*y
00419 y3=y*y2
00420 y4=y*y3
00421 epy=one+y
00422 emy=one-y
00423 epysq=one+y2
00424 emysq=one-y2
00425 emy2=emy*emy
00426 emy3=emy2*emy
00427 emy4=emy2*emy2
00428 lny=log(abs(y))
00429 abr1=epysq-four*y
00430 lnyimag=pi
00431 dalQED2 = (-tenthird + hundertfourthird*y/emy2
00432 & + sixteenzeta3 *(one - four*y2/emy4)
00433 & +( + eightthird*y*(two+seven*y-twentytwo*y2+six*y3)*
00434 & (lny**2-lnyimag**2)
00435 & - four*emysq*(epysq-eight*y)*lny
00436 & +(- sixteenthird*(log(emy)
00437 & + two*log(epy))*(lny*(epysq*lny - two*emysq)
00438 & -lnyimag*(epysq*lnyimag))
00439 & + thirtytwothird*(ddilog(y) + two*ddilog(-y))
00440 & *(emysq - two*epysq*lny)
00441 & + thirtytwo*epysq*(dli3(y) + two*dli3(-y)))*abr1)/emy4)
00442 imag = lnyimag*(
00443 & +( + eightthird*y*(two+seven*y-twentytwo*y2+six*y3)*two*lny
00444 & - four*emysq*(epysq-eight*y)
00445 & +(- sixteenthird*(log(emy) + two*log(epy))
00446 & *(two*epysq*lny - two*emysq)
00447 & - two*thirtytwothird*(ddilog(y) + two*ddilog(-y))*epysq
00448 & )*abr1)/emy4)
00449 else
00450 sint2=sp/s0
00451 sint =dsqrt(sint2)
00452 sint3=sint2*sint
00453 sint4=sint3*sint
00454 tau=dasin(sint)
00455 cost=cos(tau)
00456 phi=two*tau
00457 chi=pi-phi
00458 sint4i=one/sint4
00459 fac1=one-fourth*sint4i
00460 term2=cost*(one+two*sint2)/sint3
00461 dalQED2 = sixteen*
00462 & (two*Clausen3(phi)+four*Clausen3(chi)+zeta3)*fac1
00463 & + sixteenthird*(Clausen2(phi)-two*Clausen2(chi))*
00464 & (eight*tau*fac1-term2)
00465 & + thirtytwothird*(log(two*sint)+two*log(two*cost))*
00466 & (two*tau*fac1-term2)*tau -tenthird
00467 & + four*tau*cost*(three+two*sint2)/sint3
00468 & - twentysixthird/sint2 + (fourteenthird/sint4
00469 & + sixteenthird/sint2 - thirtytwo)*tau**2
00470 imag = null
00471 endif
00472 C SPACE-LIKE BRANCH
00473 else
00474 if (-z.gt.large) then
00475 C HIGH ENERGY ASYMPTOTE (SPACE-LIKE)
00476 c write (*,*) ' HIGH ENERGY ASYMPTOTE (SPACE-LIKE)'
00477 m2ds=-m2/sp
00478 m2ds2=m2ds*m2ds
00479 lnm2ds=log(m2ds)
00480 dalQED2 = -four*lnm2ds-tenthird+sixteenzeta3
00481 & +fourtyeight*m2ds*lnm2ds
00482 & -m2ds2*(eightthird+sixtyfourzeta3+fourty*lnm2ds
00483 & -fourtyeight*lnm2ds**2)
00484 imag = null
00485 else
00486 r=sqrt(omx)
00487 y=(r-one)/(r+one)
00488 C for analytic continuation y = y + i epsilon
00489 y2=y*y
00490 y3=y*y2
00491 epy=one+y
00492 emy=one-y
00493 epysq=one+y2
00494 emysq=one-y2
00495 emy2=emy*emy
00496 emy3=emy2*emy
00497 emy4=emy2*emy2
00498 lny=log(y)
00499 abr1=epysq-four*y
00500 dalQED2 = (-tenthird + hundertfourthird*y/emy2
00501 & + sixteenzeta3 *(one - four*y2/emy4)
00502 & +( + eightthird*y*(two+seven*y-twentytwo*y2+six*y3)*lny**2
00503 & - four*emysq*(epysq-eight*y)*lny
00504 & +(- sixteenthird*(log(emy)
00505 & + two*log(epy))*lny*(epysq*lny - two*emysq)
00506 & + thirtytwothird*(ddilog(y) + two*ddilog(-y))
00507 & *(emysq - two*epysq*lny)
00508 & + thirtytwo*epysq*(dli3(y) + two*dli3(-y)))*abr1)/emy4)
00509 imag = null
00510 endif
00511 endif
00512 9 dalQED2ferm=dalQED2*adp2*onesixteenth
00513 imag=imag*adp2*onesixteenth
00514 RETURN
00515 2 null=0.d0
00516 one=1.d0
00517 two=2.d0
00518 three=3.d0
00519 four=4.d0
00520 five=5.d0
00521 six=6.d0
00522 seven=7.d0
00523 eight=8.d0
00524 sixteen=16.d0
00525 fourty=40.d0
00526 twentytwo=22.d0
00527 thirtytwo=32.d0
00528 fourtyeight=48.d0
00529 third=1.d0/3.d0
00530 fourth=0.25d0
00531 onesixteenth=1.d0/16.d0
00532 tenthird=10.d0*third
00533 hundertfourthird=104.d0*third
00534 sixteenzeta3=16.d0*zeta3
00535 eightthird=8.d0*third
00536 fourteenthird=14.d0*third
00537 sixteenthird=16.d0*third
00538 twentysixthird=26.d0*third
00539 thirtytwothird=32.d0*third
00540 sixtyfourzeta3=64.d0*zeta3
00541 c1=-328.d0/81.d0
00542 c2=-449.d0/675.d0
00543 ini=1
00544 c write (*,*) ' numerical constants initialized'
00545 goto 1
00546 END
00547 C
00548 function dallepQED3l(sp,r1real,r1imag,r2real,r2imag,r3imag)
00549 C Warning: only works for large sp >> m_tau**2
00550 c QED VP contribution 1-,2-,3-loop
00551 c r1 1-loop contribution from light lepton
00552 c r2 2-loop contribution from light leptons (Källen/Sabry 55)
00553 c r3 3-loop contribution from light leptons (Steinhauser 98, PLB429(1998)158)
00554 c M1,M2 on-shell masses;
00555 c M1 valence lepton (current-current fermion loop), M2 sea lepton (extra fermion loop)
00556 c Piq2 = adp/4.d0*(Pi_0+adp*Pi_1
00557 c & + adp**2*(PiA_2+Pil+PiF_2+Pih_2))
00558 implicit none
00559 integer i,nl,testprint
00560 real *8 dallepQED3l,sp,q2,M1,M12,M22,sabs,sign,fac1,fac2,fac3
00561 real *8 alpha,Lnqmone,Lnqmtwo,Lnqmone2,M12dq2
00562 real *8 iLnqmone,iLnqmtwo
00563 real *8 r1real,r1imag,r2real,r2imag,r3real,r3imag
00564 real *8 rPi_0,rPi_1,rPi_2,rPiA_2,rPil_2,rPiF_2,rPih_2,rPiq2
00565 real *8 iPi_0,iPi_1,iPi_2,iPiA_2,iPil_2,iPiF_2,iPih_2,iPiq2
00566 real *8 e,x0,y0,x1,y1,z,
00567 & rxa,ixa,rxf,ixf,rxle,rxhe,ipih_2m,ipih_2t,
00568 & ixle,ixhe,rxlm,rxhm,ixlm,ixhm,rxlt,ipil_2e,ipil_2m,
00569 & ixlt,rxht,ixht
00570 include 'Phokhara/common.inc'
00571 C testprint=1: print individual terms to fort.2, fort.3
00572 testprint=0
00573 fac1=-adp /4.d0
00574 fac2=-adp2/4.d0
00575 fac3=-adp3/4.d0
00576 sabs=dabs(sp)
00577 e=sqrt(sabs)
00578 dallepQED3l=0.d0
00579 r1real=0.d0
00580 r1imag=0.d0
00581 r2real=0.d0
00582 r2imag=0.d0
00583 r3real=0.d0
00584 r3imag=0.d0
00585 x0 = 0.d0
00586 y0 = 0.d0
00587 x1 = 0.d0
00588 y1 = 0.d0
00589 rxa = 0.d0
00590 ixa = 0.d0
00591 rxf = 0.d0
00592 ixf = 0.d0
00593 if ((LEPTONflag.eq.'had').or.(LEPTONflag.eq.'top')) then
00594 return
00595 endif
00596 nl=1
00597 iLEP=LFLAG(LEPTONflag)
00598 if ((LEPTONflag.eq.'all').or.(LEPTONflag.eq.'lep')) nl=3
00599 do i=1,nl
00600 M1=ml(i)
00601 if (nl.eq.1) M1=ml(iLEP)
00602 M12=M1**2
00603 z=sp/M12
00604 if (abs(z).lt.large_3) then
00605 if (iwarnings.ne.0) then
00606 write (*,*) ' 3--loop high energy approximation'
00607 write (*,*) ' out of range: result=0.0 at energy= ',e
00608 endif
00609 return
00610 endif
00611 x0=x0+rpi_0(sp,M12,ipi_0)
00612 y0=y0+ipi_0
00613 x1=x1+rpi_1(sp,M12,ipi_1)
00614 y1=y1+ipi_1
00615 rxa=rxa+rpia_2(sp,M12,ipia_2)
00616 ixa=ixa+ipia_2
00617 rxf=rxf+rpif_2(sp,M12,ipif_2)
00618 ixf=ixf+ipif_2
00619 if (testprint.eq.1) then
00620 write (2,*) ' m:',M1,' A:',fac3*rpia_2(sp,M12,ipia_2),
00621 & ' F:',fac3*rpif_2(sp,M12,ipif_2)
00622 endif
00623 enddo
00624 C light-heavy contribution:
00625 C electron
00626 rxle= 0.d0
00627 rxhe= 0.d0
00628 ixle= 0.d0
00629 ixhe= 0.d0
00630 if ((iLEP.eq.-1).or.(iLEP.eq.1)) then
00631 rxhe=rpih_2(sp,ml(2)**2,ipih_2m)+rpih_2(sp,ml(3)**2,ipih_2t)
00632 ixhe=ipih_2m+ipih_2t
00633 if (testprint.eq.1) then
00634 write (2,*) ' ele:',
00635 & ' h-m:',fac3*rpih_2(sp,ml(2)**2,ipih_2m),
00636 & ' h-t:',fac3*rpih_2(sp,ml(3)**2,ipih_2t)
00637 endif
00638 endif
00639 C muon
00640 rxlm= 0.d0
00641 rxhm= 0.d0
00642 ixlm= 0.d0
00643 ixhm= 0.d0
00644 if ((iLEP.eq.-1).or.(iLEP.eq.2)) then
00645 rxlm=rpil_2(sp,ml(2)**2,ml(1)**2,ipil_2)
00646 rxhm=rpih_2(sp,ml(3)**2,ipih_2)
00647 ixlm=ipil_2
00648 ixhm=ipih_2
00649 if (testprint.eq.1) then
00650 write (2,*) ' muo:',
00651 & ' l-e:',fac3*rpil_2(sp,ml(2)**2,ml(1)**2,ipil_2),
00652 & ' h-t:',fac3*rpih_2(sp,ml(3)**2,ipih_2)
00653 endif
00654 endif
00655 C tau
00656 rxlt= 0.d0
00657 ixlt= 0.d0
00658 rxht= 0.d0
00659 ixht= 0.d0
00660 if ((iLEP.eq.-1).or.(iLEP.eq.3)) then
00661 rxlt=rpil_2(sp,ml(3)**2,ml(1)**2,ipil_2e)
00662 & +rpil_2(sp,ml(3)**2,ml(2)**2,ipil_2m)
00663 ixlt=ipil_2e
00664 & +ipil_2m
00665 if (testprint.eq.1) then
00666 write (2,*) ' tau:',
00667 & ' l-e:',fac3*rpil_2(sp,ml(3)**2,ml(1)**2,ipil_2e),
00668 & ' l-m:',fac3*rpil_2(sp,ml(3)**2,ml(2)**2,ipil_2m)
00669 endif
00670 endif
00671 C summing up
00672 rPi_2=rxa+rxf+rxle+rxlm+rxlt+rxhe+rxhm+rxht
00673 iPi_2=ixa+ixf+ixle+ixlm+ixlt+ixhe+ixhm+ixht
00674 if (testprint.eq.1) then
00675 write (3,*) LEPTONflag
00676 write (3,*) fac3*rxa,fac3*rxf
00677 write (3,*) fac3*rxle,fac3*rxlm,fac3*rxlt
00678 write (3,*) fac3*rxhe,fac3*rxhm,fac3*rxht
00679 endif
00680 C
00681 r1real=fac1*x0
00682 r1imag=fac1*y0
00683 r2real=fac2*x1
00684 r2imag=fac2*y1
00685 r3real=fac3*rPi_2
00686 r3imag=fac3*iPi_2
00687 dallepQED3l=r3real
00688 return
00689 end
00690 C
00691 function rpi_0(sp,M12,ipi_0)
00692 implicit none
00693 real *8 rpi_0,ipi_0,sp,sabs,M12,r12,Lnqmone,M12ds
00694 real *8 iLnqmone
00695 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00696 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00697 sabs=abs(sp)
00698 M12ds=M12/sp
00699 Lnqmone=log(sabs/M12)
00700 iLnqmone=0.d0
00701 if (sp.gt.4.d0*M12) iLnqmone=-pi
00702 rPi_0 = 20.d0/9.d0 - 4.d0/3.d0*Lnqmone+8.d0*M12ds
00703 iPi_0 = - 4.d0/3.d0*iLnqmone
00704 return
00705 end
00706
00707 function rpi_1(sp,M12,ipi_1)
00708 implicit none
00709 real *8 rpi_1,ipi_1,sp,sabs,M12,r12,Lnqmone,M12ds
00710 real *8 iLnqmone
00711 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00712 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00713 sabs=abs(sp)
00714 M12ds=M12/sp
00715 Lnqmone=log(sabs/M12)
00716 iLnqmone=0.d0
00717 if (sp.gt.4.d0*M12) iLnqmone=-pi
00718 rPi_1 = 5.d0/6.d0 - 4.d0*zeta3
00719 & - Lnqmone - 12.d0*M12ds*Lnqmone
00720 iPi_1 =
00721 & - iLnqmone - 12.d0*M12ds*iLnqmone
00722 return
00723 end
00724
00725 function rpia_2(sp,M12,ipia_2)
00726 implicit none
00727 real *8 rpia_2,ipia_2,sp,sabs,M12,r12,Lnqmone
00728 real *8 iLnqmone
00729 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00730 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00731 sabs=abs(sp)
00732 r12=sabs/M12
00733 Lnqmone=log(r12)
00734 iLnqmone=0.d0
00735 if (sp.gt.4.d0*M12) iLnqmone=-pi
00736 rPiA_2 = - 121.d0/48.d0 + (- 5.d0 + 8.d0*ln2)*zeta2
00737 & - 99.d0/16.d0*zeta3 + 10.d0*zeta5 + 1.d0/8.d0*Lnqmone
00738 iPiA_2= + 1.d0/8.d0*iLnqmone
00739 return
00740 end
00741
00742 function rpil_2(sp,M12,M22,ipil_2)
00743 implicit none
00744 real *8 rpil_2,ipil_2,sp,sabs,M12,r12,Lnqmone,M22,r22,Lnqmtwo
00745 real *8 iLnqmone,iLnqmtwo
00746 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00747 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00748 sabs=abs(sp)
00749 r12=sabs/M12
00750 r22=sabs/M22
00751 Lnqmone=log(r12)
00752 Lnqmtwo=log(r22)
00753 iLnqmone=0.d0
00754 iLnqmtwo=0.d0
00755 if (sp.gt.4.d0*M12) iLnqmone=-pi
00756 if (sp.gt.4.d0*M22) iLnqmtwo=-pi
00757 rPil_2 = - 116.d0/27.d0 + 4.d0/3.d0*zeta2 + 38.d0/9.d0*zeta3
00758 & + 14.d0/9.d0*Lnqmone + (5.d0/18.d0
00759 & - 4.d0/3.d0*zeta3)*Lnqmtwo
00760 & + 1.d0/6.d0*(Lnqmone**2-iLnqmone**2)
00761 & - 1.d0/3.d0*(Lnqmone*Lnqmtwo-iLnqmone*iLnqmtwo)
00762 iPil_2=
00763 & + 14.d0/9.d0*iLnqmone + (5.d0/18.d0
00764 & - 4.d0/3.d0*zeta3)*iLnqmtwo
00765 & + 1.d0/6.d0*2.d0*Lnqmone*iLnqmone
00766 & - 1.d0/3.d0*(Lnqmone*iLnqmtwo+iLnqmone*Lnqmtwo)
00767 return
00768 end
00769
00770 function rpif_2(sp,M12,ipif_2)
00771 implicit none
00772 real *8 rpif_2,ipif_2,sp,sabs,M12,r12,Lnqmone
00773 real *8 iLnqmone
00774 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00775 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00776 sabs=abs(sp)
00777 r12=sabs/M12
00778 Lnqmone=log(r12)
00779 iLnqmone=0.d0
00780 if (sp.gt.4.d0*M12) iLnqmone=-pi
00781 rPiF_2 = - 307.d0/216.d0 - 8.d0/3.d0*zeta2 + 545.d0/144.d0*zeta3
00782 & + (11.d0/6.d0 - 4.d0/3.d0*zeta3)*Lnqmone
00783 & - 1.d0/6.d0*(Lnqmone**2-iLnqmone**2)
00784 iPiF_2=
00785 & + (11.d0/6.d0 - 4.d0/3.d0*zeta3)*iLnqmone
00786 & - 1.d0/6.d0*2.d0*Lnqmone*iLnqmone
00787 return
00788 end
00789
00790 function rpih_2(sp,M22,ipih_2)
00791 implicit none
00792 real *8 rpih_2,ipih_2,sp,sabs,M22,r22,Lnqmtwo
00793 real *8 iLnqmtwo
00794 real *8 pi,pi2,ln2,zeta2,zeta3,zeta5
00795 common /consts/pi,pi2,ln2,zeta2,zeta3,zeta5
00796 sabs=abs(sp)
00797 r22=sabs/M22
00798 Lnqmtwo=log(r22)
00799 iLnqmtwo=0.d0
00800 if (sp.gt.4.d0*M22) iLnqmtwo=-pi
00801 rPih_2 = - 37.d0/6.d0 + 38.d0/9.d0*zeta3
00802 & + (11.d0/6.d0 - 4.d0/3.d0*zeta3)*Lnqmtwo
00803 & - 1.d0/6.d0*(Lnqmtwo**2-iLnqmtwo**2)
00804 iPih_2 =
00805 & + (11.d0/6.d0 - 4.d0/3.d0*zeta3)*Lnqmtwo
00806 & - 1.d0/6.d0*2.d0*(Lnqmtwo*iLnqmtwo)
00807 return
00808 end
00809 C
00810 function ddilog(x)
00811 C ******************************************************************
00812 C * *
00813 C * program for calculating *
00814 C * the real part of the dilogarithm for real arguments *
00815 C * *
00816 C * functions: ddilog(x) *
00817 C * rli2(x) *
00818 C * clausen2(phi) *
00819 C * cl2(phi) *
00820 C * *
00821 C * F. Jegerlehner, Paul Scherrer Institute *
00822 C * *
00823 C * Version: 25-OCT-1990 *
00824 C * *
00825 C ******************************************************************
00826 implicit none
00827 real *8 pi6,null,one,two,half,qua,x,r
00828 real *8 ddilog,rli2,sparg,omx,clo
00829 common /polylog1/pi6,null,one,two,half,qua
00830 save /polylog1/
00831 data pi6/1.644934066848226d0/
00832 data null,one,two,half,qua/0.d0,1.d0,2.d0,.5d0,.25d0/
00833 ddilog=0.d0
00834 if (x.eq.one) then
00835 ddilog=pi6
00836 return
00837 endif
00838 r=dabs(x)
00839 if (r.le.half) then
00840 ddilog=rli2(x)
00841 return
00842 else if (x.le.-two) then
00843 sparg=one/x
00844 ddilog=-rli2(sparg)-half*dlog(-x)**2-pi6
00845 return
00846 else if (x.ge.two) then
00847 sparg=one/x
00848 ddilog=-rli2(sparg)-half*dlog(x)**2+two*pi6
00849 return
00850 else if (x.lt.-one) then
00851 sparg=one/(one-x)
00852 ddilog=rli2(sparg)-dlog(-x)*dlog(one-x)
00853 & +half*dlog(one-x)**2-pi6
00854 return
00855 else if (x.lt.-half) then
00856 omx =one-x
00857 sparg=one-one/omx
00858 ddilog=-rli2(sparg)-half*dlog(omx)**2
00859 return
00860 else if (x.gt.one) then
00861 clo=dlog(x)
00862 sparg=one-one/x
00863 ddilog=rli2(sparg)+half*clo**2
00864 & -clo*dlog(x-one)+pi6
00865 return
00866 else if (x.gt.half) then
00867 sparg=one-x
00868 ddilog=-rli2(sparg)-dlog(x)*dlog(sparg)+pi6
00869 return
00870 endif
00871 end
00872
00873 function rli2(x)
00874 c Spence function for real arguments of modulus smaller than one
00875 real *8 one,half,qua,x,rli2,b,z,z2
00876 integer ini
00877 dimension b(10)
00878 common /polylog2/one,half,qua,b,ini
00879 save /polylog2/
00880 data ini/0/,one,half,qua/1.d0,.5d0,.25d0/
00881 if (ini.eq.0) goto 2
00882 1 z=-dlog(one-x)
00883 z2=z*z
00884 rli2=z*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*b(10)+b(9))
00885 & +b(8))+b(7))+b(6))+b(5))+b(4))+b(3))+b(2))+b(1))
00886 & +one-z*qua)
00887 return
00888 2 b(1) = 2.7777777777777778d-02
00889 b(2) = -2.7777777777777778d-04
00890 b(3) = 4.7241118669690098d-06
00891 b(4) = -9.1857730746619635d-08
00892 b(5) = 1.8978869988970999d-09
00893 b(6) = -4.0647616451442255d-11
00894 b(7) = 8.9216910204564526d-13
00895 b(8) = -1.9939295860721076d-14
00896 b(9) = 4.5189800296199182d-16
00897 b(10)= -1.0356517612181247d-17
00898 ini=1
00899 goto 1
00900 end
00901
00902 function clausen2(phi)
00903 c Clausens integral for arbitrary real arguments ( defined as the
00904 c imaginary part of the complex Spence function on the unit circle
00905 c z=exp(i*phi) ) (2pi-periodic,odd)
00906 implicit none
00907 real *8 null,one,pi,zpi,phi,phiabs,cl2,clausen2
00908 common /polylog3/null,one,pi,zpi
00909 save /polylog3/
00910 data pi /3.141592653589793d0/,zpi /6.283185307179586d0/
00911 data null,one/0.d0,1.d0/
00912 phi=dmod(phi,zpi)
00913 if (phi.gt.pi) phi=phi-zpi
00914 phiabs=dabs(phi)
00915 if (phi.eq.null) then
00916 clausen2=null
00917 else
00918 clausen2=phi/phiabs*cl2(phiabs)
00919 endif
00920 return
00921 end
00922
00923 function cl2(phi)
00924 c Clausens integral for real arguments 0<phi<pi
00925 implicit none
00926 real *8 one,b,phi,z,z2,cl2,pi,pi2
00927 integer ini
00928 dimension b(15)
00929 common /polylog4/one,pi,pi2,b,ini
00930 save /polylog4/
00931 data pi /3.141592653589793d0/,pi2 /1.570796326794897d0/
00932 data ini/0/,one/1.d0/
00933 if (ini.eq.0) goto 2
00934 1 z=phi
00935 z2=z*z
00936 if (phi.le.pi2) then
00937 if (phi.eq.0d0) then
00938 cl2=0d0
00939 return
00940 endif
00941 cl2= z*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2* b(10)+ b(9))
00942 & + b(8))+ b(7))+ b(6))+ b(5))+ b(4))+ b(3))+ b(2))+ b(1))
00943 & +one-dlog(dabs(z)))
00944 else
00945 cl2= z*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*
00946 & (z2*(z2* b(15)+ b(14))+ b(13))+ b(12))+ b(11))+ b(10))+ b(9))
00947 & + b(8))+ b(7))+ b(6))+ b(5))+ b(4))+ b(3))+ b(2))+ b(1))
00948 & +one-dlog(dabs(z)))
00949 endif
00950 return
00951 2 b(1) = 1.3888888888888889E-02
00952 b(2) = 6.9444444444444444E-05
00953 b(3) = 7.8735197782816830E-07
00954 b(4) = 1.1482216343327454E-08
00955 b(5) = 1.8978869988970999E-10
00956 b(6) = 3.3873013709535213E-12
00957 b(7) = 6.3726364431831804E-14
00958 b(8) = 1.2462059912950673E-15
00959 b(9) = 2.5105444608999546E-17
00960 b(10)= 5.1782588060906234E-19
00961 b(11)= 1.0887357368300848E-20
00962 b(12)= 2.3257441143020872E-22
00963 b(13)= 5.0351952131473897E-24
00964 b(14)= 1.1026499294381215E-25
00965 b(15)= 2.4386585509007344E-27
00966 ini=1
00967 goto 1
00968 end
00969 C
00970 function dli3(x)
00971 C ******************************************************************
00972 C * *
00973 C * program for calculations of *
00974 C * trilogatithms and related functions *
00975 C * *
00976 C * functions: dli3(x) *
00977 C * rli3(x) *
00978 C * rs12(x) *
00979 C * clausen3(phi) *
00980 C * cl3(phi) *
00981 C * *
00982 C * F. Jegerlehner, Paul Scherrer Institute *
00983 C * *
00984 C * Version: 25-OCT-1990 *
00985 C * *
00986 C ******************************************************************
00987 c double precision calculation of the real part of the trilogarithm
00988 implicit none
00989 real *8 zeta3,zeta2,null,one,two,half,qua,i3,i6,x,r
00990 real *8 dli3,rli3,rs12,ddilog,sparg,omx
00991 real *8 clo,cloy,clom
00992 common /polylog5/zeta3,zeta2,null,one,two,half,qua,i3,i6
00993 save /polylog5/
00994 data zeta2,zeta3/1.644934066848226d0,1.20205690315959d0/
00995 data null,one,two,half,qua,i3,i6/0.d0,1.d0,2.d0,.5d0,.25d0
00996 & ,.3333333333333333d0,.1666666666666667d0/
00997 dli3=0.d0
00998 if (x.eq.one) then
00999 dli3=zeta3
01000 return
01001 endif
01002 r=dabs(x)
01003 if (r.le.half) then
01004 dli3=rli3(x)
01005 return
01006 else if (x.le.-two) then
01007 sparg=one/x
01008 clo=dlog(-x)
01009 dli3=rli3(sparg)-i6*clo**3-zeta2*clo
01010 return
01011 else if (x.ge.two) then
01012 sparg=one/x
01013 clo=dlog(x)
01014 c imaginary part ignored
01015 dli3=rli3(sparg)-i6*clo**3+two*zeta2*clo
01016 return
01017 else if (x.lt.-one) then
01018 sparg=one/(one-x)
01019 clo =dlog(-x)
01020 cloy=dlog(sparg)
01021 clom=dlog(one-sparg)
01022 dli3= rs12(sparg)-rli3(sparg)+clom*ddilog(sparg)
01023 & +i6*cloy**3+half*cloy*clom*clo-zeta2*clo
01024 return
01025 else if (x.lt.-half) then
01026 omx =one-x
01027 sparg=one-one/omx
01028 clo=dlog(omx)
01029 dli3= rs12(sparg)-rli3(sparg)-clo*ddilog(sparg)-i6*clo**3
01030 return
01031 else if (x.gt.one) then
01032 clo=dlog(x)
01033 sparg=one-one/x
01034 c imaginary part ignored
01035 dli3=-rs12(sparg)+clo*ddilog(sparg)+i3*clo**3
01036 & -half*clo**2*dlog(x-one)+zeta2*clo+zeta3
01037 return
01038 else if (x.gt.half) then
01039 clo=dlog(x)
01040 sparg=one-x
01041 dli3=-rs12(sparg)-clo*ddilog(sparg)
01042 & -half*clo**2*dlog(sparg)+zeta2*clo+zeta3
01043 return
01044 endif
01045 end
01046
01047 function rli3(x)
01048 c function Li(3) for real arguments of modulus smaller than one
01049 implicit none
01050 real *8 one,zeta3,x,rli3,b,z
01051 integer ini
01052 dimension b(21)
01053 common /polylog6/one,zeta3,b,ini
01054 save /polylog6/
01055 data ini/0/,one/1.d0/
01056 data zeta3/1.20205690315959d0/
01057 if (ini.eq.0) goto 2
01058 1 if (x.eq.one) then
01059 rli3=zeta3
01060 return
01061 endif
01062 z=-dlog(one-x)
01063 rli3=z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*(z*
01064 & (z*b(20)+b(19))+b(18))+b(17))+b(16))+b(15))+b(14))+b(13))
01065 & +b(12))+b(11))+b(10))+b( 9))+b( 8))+b( 7))+b( 6))+b( 5))+b( 4))
01066 & +b(3))+b(2))+b(1))+one)
01067 return
01068 C Coefficients for Li3:
01069 C b(0) = 1.0d0
01070 2 b(1) = - 3.75d-01
01071 b(2) = 7.8703703703703703d-02
01072 b(3) = - 8.6805555555555555d-03
01073 b(4) = 1.2962962962962963d-04
01074 b(5) = 8.1018518518518519d-05
01075 b(6) = - 3.4193571608537595d-06
01076 b(7) = - 1.3286564625850340d-06
01077 b(8) = 8.6608717561098513d-08
01078 b(9) = 2.5260875955320400d-08
01079 b(10)= - 2.1446944683640648d-09
01080 b(11)= - 5.1401106220129789d-10
01081 b(12)= 5.2495821146008294d-11
01082 b(13)= 1.0887754406636318d-11
01083 b(14)= - 1.2779396094493695d-12
01084 b(15)= - 2.3698241773087452d-13
01085 b(16)= 3.1043578879654623d-14
01086 b(17)= 5.2617586299125061d-15
01087 b(18)= - 7.5384795499492654d-16
01088 b(19)= - 1.1862322577752285d-16
01089 b(20)= 1.8316979965491383d-17
01090 b(21)= 2.7068171031837350d-18
01091 ini=1
01092 goto 1
01093 end
01094
01095 function rs12(x)
01096 c function S(1,2) for real arguments of modulus smaller than one
01097 implicit none
01098 real *8 one,two,half,qua,i6,zeta3
01099 real *8 x,rs12,dlog,dfloat,b,z,z2,ir
01100 integer ini,i
01101 common /polylog7/one,two,half,qua,i6,zeta3,b,ini
01102 save /polylog7/
01103 dimension b(10)
01104 data ini/0/,one,two,half,qua/1.d0,2.d0,.5d0,.25d0/
01105 data i6,zeta3/.1666666666666667d0,1.20205690315959d0/
01106 if (ini.eq.0) goto 2
01107 1 if (x.eq.one) then
01108 rs12=zeta3
01109 return
01110 endif
01111 z=-dlog(one-x)
01112 z2=z*z
01113 rs12=z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*b(10)+b(9))
01114 & +b(8))+b(7))+b(6))+b(5))+b(4))+b(3))+b(2))+b(1))
01115 & +half-z*i6)*half
01116 return
01117 2 b(1) = 2.7777777777777778d-02
01118 b(2) = -2.7777777777777778d-04
01119 b(3) = 4.7241118669690098d-06
01120 b(4) = -9.1857730746619635d-08
01121 b(5) = 1.8978869988970999d-09
01122 b(6) = -4.0647616451442255d-11
01123 b(7) = 8.9216910204564526d-13
01124 b(8) = -1.9939295860721076d-14
01125 b(9) = 4.5189800296199182d-16
01126 b(10)= -1.0356517612181247d-17
01127 do i=1,10
01128 ir=dfloat(i)
01129 b(i)=b(i)*(two*ir+one)/(two*ir+two)
01130 enddo
01131 ini=1
01132 goto 1
01133 end
01134
01135 function clausen3(phi)
01136 c Clausens integral for arbitrary real arguments ( defined as the
01137 c real part of the complex trilogarithm (Li3) on the unit circle
01138 c z=exp(i*phi) ) (2pi-periodic,even)
01139 implicit none
01140 real *8 null,one,pi,zpi,zeta3
01141 real *8 phi,phiabs,cl3,clausen3
01142 common /polylog8/null,one,pi,zpi,zeta3
01143 save /polylog8/
01144 data pi /3.141592653589793d0/,zpi /6.283185307179586d0/
01145 data zeta3/1.20205690315959d0/
01146 data null,one/0.d0,1.d0/
01147 phi=dmod(phi,zpi)
01148 if (phi.gt.pi) phi=phi-zpi
01149 phiabs=dabs(phi)
01150 if (phi.eq.null) then
01151 clausen3=zeta3
01152 else
01153 clausen3=cl3(phiabs)
01154 endif
01155 return
01156 end
01157
01158 function cl3(phi)
01159 c Clausens integral of 3rd order for real arguments 0<phi<pi
01160 implicit none
01161 real *8 one,half,threequa,b,phi,z,z2,cl3,pi,pi2,zeta3
01162 integer ini
01163 dimension b(15)
01164 common /polylog9/one,half,threequa,pi,pi2,zeta3,b,ini
01165 save /polylog9/
01166 data zeta3/1.20205690315959d0/
01167 data pi /3.141592653589793d0/,pi2 /1.570796326794897d0/
01168 data ini/0/,one,half,threequa/1.0d0,.50d0,.75d0/
01169 if (ini.eq.0) goto 2
01170 1 z=phi
01171 z2=z*z
01172 if (phi.le.pi2) then
01173 if (phi.eq.0d0) then
01174 cl3=zeta3
01175 return
01176 endif
01177 cl3=-z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2* b(10)+ b(9))
01178 & + b(8))+ b(7))+ b(6))+ b(5))+ b(4))+ b(3))+ b(2))+ b(1))
01179 & +threequa-half*dlog(dabs(z)))+zeta3
01180 else
01181 cl3=-z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*(z2*
01182 & (z2*(z2* b(15)+ b(14))+ b(13))+ b(12))+ b(11))+ b(10))+ b(9))
01183 & + b(8))+ b(7))+ b(6))+ b(5))+ b(4))+ b(3))+ b(2))+ b(1))
01184 & +threequa-half*dlog(dabs(z)))+zeta3
01185 endif
01186 return
01187 2 b(1) = 3.4722222222222222E-03
01188 b(2) = 1.1574074074074074E-05
01189 b(3) = 9.8418997228521037E-08
01190 b(4) = 1.1482216343327454E-09
01191 b(5) = 1.5815724990809166E-11
01192 b(6) = 2.4195009792525152E-13
01193 b(7) = 3.9828977769894878E-15
01194 b(8) = 6.9233666183059291E-17
01195 b(9) = 1.2552722304499773E-18
01196 b(10)= 2.3537540027684652E-20
01197 b(11)= 4.5363989034586869E-22
01198 b(12)= 8.9451696703926431E-24
01199 b(13)= 1.7982840046954963E-25
01200 b(14)= 3.6754997647937384E-27
01201 b(15)= 7.6208079715647952E-29
01202 ini=1
01203 goto 1
01204 end
01205
01206
01207 c;;;;;;;;;;;;;;;;;;;;;;;;;;;;; -*- Mode: Fortran -*- ;;;;;;;;;;;;;;;;;;;;;;;;;
01208 c;; hadr5new.f ---
01209 c;; Author : Fred Jegerlehner
01210 c;; Created On : Sat Jul 12 17:45:19 2003
01211 c;; Last Modified By: Fred Jegerlehner
01212 c;; Last Modified On: Sun Nov 2 00:17:36 2003
01213 c;; RCS: $Id: vac_pol_hc1.inc,v 1.1.1.1 2007/11/18 09:51:28 azhemchugov Exp $
01214 c;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
01215 c;; Copyright (c) 2003 Fred Jegerlehner
01216 c;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
01217 c;;
01218 subroutine hadr5n(e,st2,der,errder,deg,errdeg)
01219 c single precision HADR5 ; double precision DHADR5
01220 c ******************************************************************
01221 c * *
01222 c * subroutine for the evaluation of the light hadron *
01223 c * contributions to Delta_r and Delta_g *
01224 c * using fits to the *
01225 c * QED vacuum polarization from e^+ e^- data *
01226 c * *
01227 c * F. Jegerlehner, DESY, Platanenalle 6, D-15738 Zeuthen *
01228 c * *
01229 c * E-mail: fred.jegerlehner@desy.de *
01230 c * Phone : +49-33762-77-259 *
01231 c * *
01232 c * Reference: F. Jegerlehner, Z. Phys. C32 (1986) 195 *
01233 c * H. Burkhardt et al., Z. Phys. C42 (1989) 497 *
01234 c * S. Eidelman, F. Jegerlehner, Z. Phys. C (1995) *
01235 c * *
01236 c ******************************************************************
01237 c VERSION: 24/02/1995
01238 c
01239 C Notation: E energy ( momentum transfer ): E>0 timelike , E<0 spacelike
01240 C st2 is sin^2(Theta); st2=0.2322 is the reference value
01241 C the routine returns the hadronic contribution of 5 flavors (u,d,s,c,b)
01242 C to DER=Delta_r with hadronic error ERRDER
01243 C and DEG=Delta_g with hadronic error ERRDEG
01244 C The effective value of the fine structure constant alphaQED at energy
01245 C E is alphaQED(E)=alphaQED(0)/(1-Delta_r) ,similarly for the SU(2)
01246 C coupling alphaSU2(E)=alphaSU2(0)/(1-Delta_g), where Delta_r(g) is the
01247 C sum of leptonic, hadronic contributions (top to be added).
01248 C
01249 C
01250 IMPLICIT NONE
01251 INTEGER NA,NB,NC,IJ,I,ini
01252 PARAMETER(NA=882,NB=530,NC=250)
01253 real e,st2,der,errder,deg,errdeg
01254 real dal(2),dg2(2)
01255 REAL ETA(NA),DAT(NA,2),ESA(NB),DAS(NB,2),EMA(NC),DAM(NC,2)
01256 COMMON /DATT/ETA,DAT
01257 COMMON /DATS/ESA,DAS
01258 COMMON /DATM/EMA,DAM
01259 data ini/0/
01260 C initialize data
01261
01262 call dalhad_spacelike
01263
01264 call dalhad_timelike
01265
01266 call dalhad_timelike1
01267
01268 if (e.eq.0.0) then
01269 der=0.0
01270 errder=0.0
01271 deg=0.0
01272 errdeg=0.0
01273 return
01274 else if ((e.lt.ETA(NA)).and.(e.ge.ETA(1))) then
01275 ij=NA
01276 do while (ETA(ij).ge.e)
01277 ij=ij-1
01278 enddo
01279 DO I=1,2
01280 dal(i)=DAT(IJ,I)
01281 & +(DAT(IJ+1,I)-DAT(IJ,I))/(ETA(IJ+1)-ETA(IJ))*(E-ETA(IJ))
01282 ENDDO
01283 der=dal(1)
01284 errder=dal(2)
01285 deg=0.0
01286 errdeg=0.0
01287 return
01288 else if ((e.le.EMA(NC)).and.(e.gt.EMA(1))) then
01289 if ((e.gt.3.0).and.(e.lt.12.0).and.(ini.eq.0)) then
01290 write (*,*) ' ******************************************'
01291 write (*,*) ' * Warning: results may not be reliable *'
01292 write (*,*) ' * close to J/Psi and Upsilon resonances! *'
01293 write (*,*) ' ***************FJ@DESY********************'
01294 ini=1
01295 endif
01296 ij=NC
01297 do while (EMA(ij).ge.e)
01298 ij=ij-1
01299 enddo
01300 DO I=1,2
01301 dal(i)=DAM(IJ,I)
01302 & +(DAM(IJ+1,I)-DAM(IJ,I))/(EMA(IJ+1)-EMA(IJ))*(E-EMA(IJ))
01303 ENDDO
01304 der=dal(1)
01305 errder=dal(2)
01306 deg=0.0
01307 errdeg=0.0
01308 return
01309 else if ((e.le.ESA(NB)).and.(e.gt.0.0)) then
01310 ij=NB
01311 do while (ESA(ij).ge.e)
01312 ij=ij-1
01313 enddo
01314 DO I=1,2
01315 dal(i)=DAS(IJ,I)
01316 & +(DAS(IJ+1,I)-DAS(IJ,I))/(ESA(IJ+1)-ESA(IJ))*(E-ESA(IJ))
01317 ENDDO
01318 der=dal(1)
01319 errder=dal(2)
01320 deg=0.0
01321 errdeg=0.0
01322 return
01323 else
01324 write(*,*) ' out of range! '
01325 endif
01326 return
01327 end
01328
01329 subroutine dhadr5n(de,dst2,dder,derrder,ddeg,derrdeg)
01330 c ******************************************************************
01331 c * F. Jegerlehner, DESY Zeuthen, D-15738 Zeuthen *
01332 c ******************************************************************
01333 c Converts hadr5 to doubleprecision variables dhadr5
01334 c
01335 implicit none
01336 real *8 de,dst2,dder,ddeg,derrder,derrdeg
01337 real se,sst2,sder,sdeg,serrder,serrdeg
01338 se =sngl(de)
01339 sst2=sngl(dst2)
01340 call hadr5n(se,sst2,sder,serrder,sdeg,serrdeg)
01341 dder =dble(sder)
01342 ddeg =dble(sdeg)
01343 derrder=dble(serrder)
01344 derrdeg=dble(serrdeg)
01345 return
01346 end
01347
01348 subroutine dalhad_spacelike
01349 * hadronic contribution of the 5 light quark flavors evaluated using
01350 * e^+e^-data. CMD-2 2001 data have a normalization problem. I added a energy
01351 * independent +1.5 % correction and incresed the syst error by 0.5%.
01352 * Results thus are preliminary! Arrays:
01353 * ETA(I) : - energy E in GeV; the minus sign indicates the spacelike region
01354 * DAT(I,1): Delta alpha^5(-E**2)
01355 * DAT(I,2): error (systematics dominated, should be treated as a systematic error)
01356 C
01357 IMPLICIT NONE
01358 INTEGER NA,I,J
01359 PARAMETER(NA=882)
01360 REAL ETA(NA),DAT(NA,2)
01361 COMMON /DATT/ETA,DAT
01362 C 882
01363 DATA (ETA(I), I=1,100) /
01364 & -1.000E+03,-9.900E+02,-9.800E+02,-9.700E+02,-9.600E+02,
01365 & -9.500E+02,-9.400E+02,-9.300E+02,-9.200E+02,-9.100E+02,
01366 & -9.000E+02,-8.900E+02,-8.800E+02,-8.700E+02,-8.600E+02,
01367 & -8.500E+02,-8.400E+02,-8.300E+02,-8.200E+02,-8.100E+02,
01368 & -8.000E+02,-7.900E+02,-7.800E+02,-7.700E+02,-7.600E+02,
01369 & -7.500E+02,-7.400E+02,-7.300E+02,-7.200E+02,-7.100E+02,
01370 & -7.000E+02,-6.900E+02,-6.800E+02,-6.700E+02,-6.600E+02,
01371 & -6.500E+02,-6.400E+02,-6.300E+02,-6.200E+02,-6.100E+02,
01372 & -6.000E+02,-5.900E+02,-5.800E+02,-5.700E+02,-5.600E+02,
01373 & -5.500E+02,-5.400E+02,-5.300E+02,-5.200E+02,-5.100E+02,
01374 & -5.000E+02,-4.900E+02,-4.800E+02,-4.700E+02,-4.600E+02,
01375 & -4.500E+02,-4.400E+02,-4.300E+02,-4.200E+02,-4.100E+02,
01376 & -4.000E+02,-3.900E+02,-3.800E+02,-3.700E+02,-3.600E+02,
01377 & -3.500E+02,-3.400E+02,-3.300E+02,-3.200E+02,-3.100E+02,
01378 & -3.000E+02,-2.900E+02,-2.800E+02,-2.700E+02,-2.600E+02,
01379 & -2.500E+02,-2.400E+02,-2.300E+02,-2.200E+02,-2.100E+02,
01380 & -2.000E+02,-1.900E+02,-1.800E+02,-1.700E+02,-1.600E+02,
01381 & -1.500E+02,-1.400E+02,-1.300E+02,-1.200E+02,-1.100E+02,
01382 & -1.000E+02,-9.900E+01,-9.800E+01,-9.700E+01,-9.600E+01,
01383 & -9.500E+01,-9.400E+01,-9.300E+01,-9.200E+01,-9.100E+01/
01384
01385 DATA (ETA(I), I=101,200) /
01386 & -9.000E+01,-8.900E+01,-8.800E+01,-8.700E+01,-8.600E+01,
01387 & -8.500E+01,-8.400E+01,-8.300E+01,-8.200E+01,-8.100E+01,
01388 & -8.000E+01,-7.900E+01,-7.800E+01,-7.700E+01,-7.600E+01,
01389 & -7.500E+01,-7.400E+01,-7.300E+01,-7.200E+01,-7.100E+01,
01390 & -7.000E+01,-6.900E+01,-6.800E+01,-6.700E+01,-6.600E+01,
01391 & -6.500E+01,-6.400E+01,-6.300E+01,-6.200E+01,-6.100E+01,
01392 & -6.000E+01,-5.900E+01,-5.800E+01,-5.700E+01,-5.600E+01,
01393 & -5.500E+01,-5.400E+01,-5.300E+01,-5.200E+01,-5.100E+01,
01394 & -5.000E+01,-4.900E+01,-4.800E+01,-4.700E+01,-4.600E+01,
01395 & -4.500E+01,-4.400E+01,-4.300E+01,-4.200E+01,-4.100E+01,
01396 & -4.000E+01,-3.900E+01,-3.800E+01,-3.700E+01,-3.600E+01,
01397 & -3.500E+01,-3.400E+01,-3.300E+01,-3.200E+01,-3.100E+01,
01398 & -3.000E+01,-2.900E+01,-2.800E+01,-2.700E+01,-2.600E+01,
01399 & -2.500E+01,-2.400E+01,-2.300E+01,-2.200E+01,-2.100E+01,
01400 & -2.000E+01,-1.900E+01,-1.800E+01,-1.700E+01,-1.600E+01,
01401 & -1.500E+01,-1.400E+01,-1.300E+01,-1.200E+01,-1.100E+01,
01402 & -1.000E+01,-9.910E+00,-9.820E+00,-9.730E+00,-9.640E+00,
01403 & -9.550E+00,-9.460E+00,-9.370E+00,-9.280E+00,-9.190E+00,
01404 & -9.100E+00,-9.010E+00,-8.920E+00,-8.830E+00,-8.740E+00,
01405 & -8.650E+00,-8.560E+00,-8.470E+00,-8.380E+00,-8.290E+00/
01406
01407 DATA (ETA(I), I=201,300) /
01408 & -8.200E+00,-8.110E+00,-8.020E+00,-7.930E+00,-7.840E+00,
01409 & -7.750E+00,-7.660E+00,-7.570E+00,-7.480E+00,-7.390E+00,
01410 & -7.300E+00,-7.210E+00,-7.120E+00,-7.030E+00,-6.940E+00,
01411 & -6.850E+00,-6.760E+00,-6.670E+00,-6.580E+00,-6.490E+00,
01412 & -6.400E+00,-6.310E+00,-6.220E+00,-6.130E+00,-6.040E+00,
01413 & -5.950E+00,-5.860E+00,-5.770E+00,-5.680E+00,-5.590E+00,
01414 & -5.500E+00,-5.410E+00,-5.320E+00,-5.230E+00,-5.140E+00,
01415 & -5.050E+00,-4.960E+00,-4.870E+00,-4.780E+00,-4.690E+00,
01416 & -4.600E+00,-4.510E+00,-4.420E+00,-4.330E+00,-4.240E+00,
01417 & -4.150E+00,-4.060E+00,-3.970E+00,-3.880E+00,-3.790E+00,
01418 & -3.700E+00,-3.610E+00,-3.520E+00,-3.430E+00,-3.340E+00,
01419 & -3.250E+00,-3.160E+00,-3.070E+00,-2.980E+00,-2.890E+00,
01420 & -2.800E+00,-2.710E+00,-2.620E+00,-2.530E+00,-2.440E+00,
01421 & -2.350E+00,-2.260E+00,-2.170E+00,-2.080E+00,-1.990E+00,
01422 & -1.900E+00,-1.810E+00,-1.720E+00,-1.630E+00,-1.540E+00,
01423 & -1.450E+00,-1.360E+00,-1.270E+00,-1.180E+00,-1.090E+00,
01424 & -1.000E+00,-9.910E-01,-9.820E-01,-9.730E-01,-9.640E-01,
01425 & -9.550E-01,-9.460E-01,-9.370E-01,-9.280E-01,-9.190E-01,
01426 & -9.100E-01,-9.010E-01,-8.920E-01,-8.830E-01,-8.740E-01,
01427 & -8.650E-01,-8.560E-01,-8.470E-01,-8.380E-01,-8.290E-01/
01428
01429 DATA (ETA(I), I=301,400) /
01430 & -8.200E-01,-8.110E-01,-8.020E-01,-7.930E-01,-7.840E-01,
01431 & -7.750E-01,-7.660E-01,-7.570E-01,-7.480E-01,-7.390E-01,
01432 & -7.300E-01,-7.210E-01,-7.120E-01,-7.030E-01,-6.940E-01,
01433 & -6.850E-01,-6.760E-01,-6.670E-01,-6.580E-01,-6.490E-01,
01434 & -6.400E-01,-6.310E-01,-6.220E-01,-6.130E-01,-6.040E-01,
01435 & -5.950E-01,-5.860E-01,-5.770E-01,-5.680E-01,-5.590E-01,
01436 & -5.500E-01,-5.410E-01,-5.320E-01,-5.230E-01,-5.140E-01,
01437 & -5.050E-01,-4.960E-01,-4.870E-01,-4.780E-01,-4.690E-01,
01438 & -4.600E-01,-4.510E-01,-4.420E-01,-4.330E-01,-4.240E-01,
01439 & -4.150E-01,-4.060E-01,-3.970E-01,-3.880E-01,-3.790E-01,
01440 & -3.700E-01,-3.610E-01,-3.520E-01,-3.430E-01,-3.340E-01,
01441 & -3.250E-01,-3.160E-01,-3.070E-01,-2.980E-01,-2.890E-01,
01442 & -2.800E-01,-2.710E-01,-2.620E-01,-2.530E-01,-2.440E-01,
01443 & -2.350E-01,-2.260E-01,-2.170E-01,-2.080E-01,-1.990E-01,
01444 & -1.900E-01,-1.810E-01,-1.720E-01,-1.630E-01,-1.540E-01,
01445 & -1.450E-01,-1.360E-01,-1.270E-01,-1.180E-01,-1.090E-01,
01446 & -1.000E-01,-9.910E-02,-9.820E-02,-9.730E-02,-9.640E-02,
01447 & -9.550E-02,-9.460E-02,-9.370E-02,-9.280E-02,-9.190E-02,
01448 & -9.100E-02,-9.010E-02,-8.920E-02,-8.830E-02,-8.740E-02,
01449 & -8.650E-02,-8.560E-02,-8.470E-02,-8.380E-02,-8.290E-02/
01450
01451 DATA (ETA(I), I=401,500) /
01452 & -8.200E-02,-8.110E-02,-8.020E-02,-7.930E-02,-7.840E-02,
01453 & -7.750E-02,-7.660E-02,-7.570E-02,-7.480E-02,-7.390E-02,
01454 & -7.300E-02,-7.210E-02,-7.120E-02,-7.030E-02,-6.940E-02,
01455 & -6.850E-02,-6.760E-02,-6.670E-02,-6.580E-02,-6.490E-02,
01456 & -6.400E-02,-6.310E-02,-6.220E-02,-6.130E-02,-6.040E-02,
01457 & -5.950E-02,-5.860E-02,-5.770E-02,-5.680E-02,-5.590E-02,
01458 & -5.500E-02,-5.410E-02,-5.320E-02,-5.230E-02,-5.140E-02,
01459 & -5.050E-02,-4.960E-02,-4.870E-02,-4.780E-02,-4.690E-02,
01460 & -4.600E-02,-4.510E-02,-4.420E-02,-4.330E-02,-4.240E-02,
01461 & -4.150E-02,-4.060E-02,-3.970E-02,-3.880E-02,-3.790E-02,
01462 & -3.700E-02,-3.610E-02,-3.520E-02,-3.430E-02,-3.340E-02,
01463 & -3.250E-02,-3.160E-02,-3.070E-02,-2.980E-02,-2.890E-02,
01464 & -2.800E-02,-2.710E-02,-2.620E-02,-2.530E-02,-2.440E-02,
01465 & -2.350E-02,-2.260E-02,-2.170E-02,-2.080E-02,-1.990E-02,
01466 & -1.900E-02,-1.810E-02,-1.720E-02,-1.630E-02,-1.540E-02,
01467 & -1.450E-02,-1.360E-02,-1.270E-02,-1.180E-02,-1.090E-02,
01468 & -1.000E-02,-9.910E-03,-9.820E-03,-9.730E-03,-9.640E-03,
01469 & -9.550E-03,-9.460E-03,-9.370E-03,-9.280E-03,-9.190E-03,
01470 & -9.100E-03,-9.010E-03,-8.920E-03,-8.830E-03,-8.740E-03,
01471 & -8.650E-03,-8.560E-03,-8.470E-03,-8.380E-03,-8.290E-03/
01472
01473 DATA (ETA(I), I=501,600) /
01474 & -8.200E-03,-8.110E-03,-8.020E-03,-7.930E-03,-7.840E-03,
01475 & -7.750E-03,-7.660E-03,-7.570E-03,-7.480E-03,-7.390E-03,
01476 & -7.300E-03,-7.210E-03,-7.120E-03,-7.030E-03,-6.940E-03,
01477 & -6.850E-03,-6.760E-03,-6.670E-03,-6.580E-03,-6.490E-03,
01478 & -6.400E-03,-6.310E-03,-6.220E-03,-6.130E-03,-6.040E-03,
01479 & -5.950E-03,-5.860E-03,-5.770E-03,-5.680E-03,-5.590E-03,
01480 & -5.500E-03,-5.410E-03,-5.320E-03,-5.230E-03,-5.140E-03,
01481 & -5.050E-03,-4.960E-03,-4.870E-03,-4.780E-03,-4.690E-03,
01482 & -4.600E-03,-4.510E-03,-4.420E-03,-4.330E-03,-4.240E-03,
01483 & -4.150E-03,-4.060E-03,-3.970E-03,-3.880E-03,-3.790E-03,
01484 & -3.700E-03,-3.610E-03,-3.520E-03,-3.430E-03,-3.340E-03,
01485 & -3.250E-03,-3.160E-03,-3.070E-03,-2.980E-03,-2.890E-03,
01486 & -2.800E-03,-2.710E-03,-2.620E-03,-2.530E-03,-2.440E-03,
01487 & -2.350E-03,-2.260E-03,-2.170E-03,-2.080E-03,-1.990E-03,
01488 & -1.900E-03,-1.810E-03,-1.720E-03,-1.630E-03,-1.540E-03,
01489 & -1.450E-03,-1.360E-03,-1.270E-03,-1.180E-03,-1.090E-03,
01490 & -1.000E-03,-9.910E-04,-9.820E-04,-9.730E-04,-9.640E-04,
01491 & -9.550E-04,-9.460E-04,-9.370E-04,-9.280E-04,-9.190E-04,
01492 & -9.100E-04,-9.010E-04,-8.920E-04,-8.830E-04,-8.740E-04,
01493 & -8.650E-04,-8.560E-04,-8.470E-04,-8.380E-04,-8.290E-04/
01494
01495 DATA (ETA(I), I=601,700) /
01496 & -8.200E-04,-8.110E-04,-8.020E-04,-7.930E-04,-7.840E-04,
01497 & -7.750E-04,-7.660E-04,-7.570E-04,-7.480E-04,-7.390E-04,
01498 & -7.300E-04,-7.210E-04,-7.120E-04,-7.030E-04,-6.940E-04,
01499 & -6.850E-04,-6.760E-04,-6.670E-04,-6.580E-04,-6.490E-04,
01500 & -6.400E-04,-6.310E-04,-6.220E-04,-6.130E-04,-6.040E-04,
01501 & -5.950E-04,-5.860E-04,-5.770E-04,-5.680E-04,-5.590E-04,
01502 & -5.500E-04,-5.410E-04,-5.320E-04,-5.230E-04,-5.140E-04,
01503 & -5.050E-04,-4.960E-04,-4.870E-04,-4.780E-04,-4.690E-04,
01504 & -4.600E-04,-4.510E-04,-4.420E-04,-4.330E-04,-4.240E-04,
01505 & -4.150E-04,-4.060E-04,-3.970E-04,-3.880E-04,-3.790E-04,
01506 & -3.700E-04,-3.610E-04,-3.520E-04,-3.430E-04,-3.340E-04,
01507 & -3.250E-04,-3.160E-04,-3.070E-04,-2.980E-04,-2.890E-04,
01508 & -2.800E-04,-2.710E-04,-2.620E-04,-2.530E-04,-2.440E-04,
01509 & -2.350E-04,-2.260E-04,-2.170E-04,-2.080E-04,-1.990E-04,
01510 & -1.900E-04,-1.810E-04,-1.720E-04,-1.630E-04,-1.540E-04,
01511 & -1.450E-04,-1.360E-04,-1.270E-04,-1.180E-04,-1.090E-04,
01512 & -1.000E-04,-9.910E-05,-9.820E-05,-9.730E-05,-9.640E-05,
01513 & -9.550E-05,-9.460E-05,-9.370E-05,-9.280E-05,-9.190E-05,
01514 & -9.100E-05,-9.010E-05,-8.920E-05,-8.830E-05,-8.740E-05,
01515 & -8.650E-05,-8.560E-05,-8.470E-05,-8.380E-05,-8.290E-05/
01516
01517 DATA (ETA(I), I=701,800) /
01518 & -8.200E-05,-8.110E-05,-8.020E-05,-7.930E-05,-7.840E-05,
01519 & -7.750E-05,-7.660E-05,-7.570E-05,-7.480E-05,-7.390E-05,
01520 & -7.300E-05,-7.210E-05,-7.120E-05,-7.030E-05,-6.940E-05,
01521 & -6.850E-05,-6.760E-05,-6.670E-05,-6.580E-05,-6.490E-05,
01522 & -6.400E-05,-6.310E-05,-6.220E-05,-6.130E-05,-6.040E-05,
01523 & -5.950E-05,-5.860E-05,-5.770E-05,-5.680E-05,-5.590E-05,
01524 & -5.500E-05,-5.410E-05,-5.320E-05,-5.230E-05,-5.140E-05,
01525 & -5.050E-05,-4.960E-05,-4.870E-05,-4.780E-05,-4.690E-05,
01526 & -4.600E-05,-4.510E-05,-4.420E-05,-4.330E-05,-4.240E-05,
01527 & -4.150E-05,-4.060E-05,-3.970E-05,-3.880E-05,-3.790E-05,
01528 & -3.700E-05,-3.610E-05,-3.520E-05,-3.430E-05,-3.340E-05,
01529 & -3.250E-05,-3.160E-05,-3.070E-05,-2.980E-05,-2.890E-05,
01530 & -2.800E-05,-2.710E-05,-2.620E-05,-2.530E-05,-2.440E-05,
01531 & -2.350E-05,-2.260E-05,-2.170E-05,-2.080E-05,-1.990E-05,
01532 & -1.900E-05,-1.810E-05,-1.720E-05,-1.630E-05,-1.540E-05,
01533 & -1.450E-05,-1.360E-05,-1.270E-05,-1.180E-05,-1.090E-05,
01534 & -1.000E-05,-9.910E-06,-9.820E-06,-9.730E-06,-9.640E-06,
01535 & -9.550E-06,-9.460E-06,-9.370E-06,-9.280E-06,-9.190E-06,
01536 & -9.100E-06,-9.010E-06,-8.920E-06,-8.830E-06,-8.740E-06,
01537 & -8.650E-06,-8.560E-06,-8.470E-06,-8.380E-06,-8.290E-06/
01538
01539 DATA (ETA(I), I=801,882) /
01540 & -8.200E-06,-8.110E-06,-8.020E-06,-7.930E-06,-7.840E-06,
01541 & -7.750E-06,-7.660E-06,-7.570E-06,-7.480E-06,-7.390E-06,
01542 & -7.300E-06,-7.210E-06,-7.120E-06,-7.030E-06,-6.940E-06,
01543 & -6.850E-06,-6.760E-06,-6.670E-06,-6.580E-06,-6.490E-06,
01544 & -6.400E-06,-6.310E-06,-6.220E-06,-6.130E-06,-6.040E-06,
01545 & -5.950E-06,-5.860E-06,-5.770E-06,-5.680E-06,-5.590E-06,
01546 & -5.500E-06,-5.410E-06,-5.320E-06,-5.230E-06,-5.140E-06,
01547 & -5.050E-06,-4.960E-06,-4.870E-06,-4.780E-06,-4.690E-06,
01548 & -4.600E-06,-4.510E-06,-4.420E-06,-4.330E-06,-4.240E-06,
01549 & -4.150E-06,-4.060E-06,-3.970E-06,-3.880E-06,-3.790E-06,
01550 & -3.700E-06,-3.610E-06,-3.520E-06,-3.430E-06,-3.340E-06,
01551 & -3.250E-06,-3.160E-06,-3.070E-06,-2.980E-06,-2.890E-06,
01552 & -2.800E-06,-2.710E-06,-2.620E-06,-2.530E-06,-2.440E-06,
01553 & -2.350E-06,-2.260E-06,-2.170E-06,-2.080E-06,-1.990E-06,
01554 & -1.900E-06,-1.810E-06,-1.720E-06,-1.630E-06,-1.540E-06,
01555 & -1.450E-06,-1.360E-06,-1.270E-06,-1.180E-06,-1.090E-06,
01556 & -1.000E-06, 0.000E+00/
01557
01558 DATA ((DAT(I,J), I=1,100), J=1,1) /
01559 & 4.380E-02, 4.372E-02, 4.365E-02, 4.357E-02, 4.349E-02,
01560 & 4.341E-02, 4.333E-02, 4.325E-02, 4.317E-02, 4.308E-02,
01561 & 4.300E-02, 4.291E-02, 4.283E-02, 4.274E-02, 4.265E-02,
01562 & 4.256E-02, 4.248E-02, 4.239E-02, 4.229E-02, 4.220E-02,
01563 & 4.211E-02, 4.201E-02, 4.192E-02, 4.182E-02, 4.172E-02,
01564 & 4.162E-02, 4.152E-02, 4.142E-02, 4.132E-02, 4.121E-02,
01565 & 4.111E-02, 4.100E-02, 4.089E-02, 4.078E-02, 4.067E-02,
01566 & 4.056E-02, 4.045E-02, 4.033E-02, 4.021E-02, 4.009E-02,
01567 & 3.997E-02, 3.985E-02, 3.972E-02, 3.960E-02, 3.947E-02,
01568 & 3.934E-02, 3.921E-02, 3.907E-02, 3.893E-02, 3.880E-02,
01569 & 3.865E-02, 3.851E-02, 3.836E-02, 3.821E-02, 3.806E-02,
01570 & 3.791E-02, 3.775E-02, 3.759E-02, 3.742E-02, 3.726E-02,
01571 & 3.709E-02, 3.691E-02, 3.673E-02, 3.655E-02, 3.636E-02,
01572 & 3.617E-02, 3.597E-02, 3.577E-02, 3.557E-02, 3.535E-02,
01573 & 3.514E-02, 3.491E-02, 3.468E-02, 3.444E-02, 3.419E-02,
01574 & 3.394E-02, 3.368E-02, 3.340E-02, 3.312E-02, 3.282E-02,
01575 & 3.251E-02, 3.219E-02, 3.185E-02, 3.150E-02, 3.112E-02,
01576 & 3.072E-02, 3.030E-02, 2.985E-02, 2.936E-02, 2.883E-02,
01577 & 2.826E-02, 2.820E-02, 2.813E-02, 2.807E-02, 2.801E-02,
01578 & 2.795E-02, 2.788E-02, 2.782E-02, 2.775E-02, 2.769E-02/
01579
01580 DATA ((DAT(I,J), I=101,200), J=1,1) /
01581 & 2.762E-02, 2.756E-02, 2.749E-02, 2.742E-02, 2.735E-02,
01582 & 2.728E-02, 2.721E-02, 2.714E-02, 2.706E-02, 2.699E-02,
01583 & 2.692E-02, 2.684E-02, 2.677E-02, 2.669E-02, 2.661E-02,
01584 & 2.653E-02, 2.645E-02, 2.637E-02, 2.629E-02, 2.620E-02,
01585 & 2.612E-02, 2.603E-02, 2.595E-02, 2.586E-02, 2.577E-02,
01586 & 2.568E-02, 2.558E-02, 2.549E-02, 2.540E-02, 2.530E-02,
01587 & 2.520E-02, 2.510E-02, 2.500E-02, 2.489E-02, 2.479E-02,
01588 & 2.468E-02, 2.457E-02, 2.446E-02, 2.435E-02, 2.423E-02,
01589 & 2.412E-02, 2.400E-02, 2.387E-02, 2.375E-02, 2.362E-02,
01590 & 2.349E-02, 2.336E-02, 2.322E-02, 2.308E-02, 2.294E-02,
01591 & 2.279E-02, 2.264E-02, 2.249E-02, 2.233E-02, 2.217E-02,
01592 & 2.200E-02, 2.183E-02, 2.165E-02, 2.147E-02, 2.128E-02,
01593 & 2.109E-02, 2.089E-02, 2.068E-02, 2.047E-02, 2.024E-02,
01594 & 2.001E-02, 1.977E-02, 1.952E-02, 1.926E-02, 1.899E-02,
01595 & 1.870E-02, 1.840E-02, 1.808E-02, 1.775E-02, 1.739E-02,
01596 & 1.702E-02, 1.662E-02, 1.619E-02, 1.572E-02, 1.522E-02,
01597 & 1.468E-02, 1.462E-02, 1.457E-02, 1.452E-02, 1.446E-02,
01598 & 1.441E-02, 1.436E-02, 1.430E-02, 1.425E-02, 1.419E-02,
01599 & 1.414E-02, 1.408E-02, 1.402E-02, 1.397E-02, 1.391E-02,
01600 & 1.385E-02, 1.379E-02, 1.373E-02, 1.367E-02, 1.361E-02/
01601
01602 DATA ((DAT(I,J), I=201,300), J=1,1) /
01603 & 1.355E-02, 1.349E-02, 1.342E-02, 1.336E-02, 1.330E-02,
01604 & 1.323E-02, 1.317E-02, 1.310E-02, 1.303E-02, 1.297E-02,
01605 & 1.290E-02, 1.283E-02, 1.276E-02, 1.269E-02, 1.262E-02,
01606 & 1.254E-02, 1.247E-02, 1.240E-02, 1.232E-02, 1.225E-02,
01607 & 1.217E-02, 1.209E-02, 1.201E-02, 1.193E-02, 1.185E-02,
01608 & 1.177E-02, 1.169E-02, 1.160E-02, 1.152E-02, 1.143E-02,
01609 & 1.134E-02, 1.125E-02, 1.116E-02, 1.107E-02, 1.098E-02,
01610 & 1.089E-02, 1.079E-02, 1.069E-02, 1.059E-02, 1.049E-02,
01611 & 1.039E-02, 1.029E-02, 1.018E-02, 1.007E-02, 9.964E-03,
01612 & 9.853E-03, 9.740E-03, 9.625E-03, 9.508E-03, 9.389E-03,
01613 & 9.267E-03, 9.143E-03, 9.017E-03, 8.888E-03, 8.756E-03,
01614 & 8.621E-03, 8.484E-03, 8.344E-03, 8.201E-03, 8.054E-03,
01615 & 7.904E-03, 7.750E-03, 7.592E-03, 7.430E-03, 7.265E-03,
01616 & 7.094E-03, 6.919E-03, 6.738E-03, 6.552E-03, 6.361E-03,
01617 & 6.162E-03, 5.957E-03, 5.744E-03, 5.523E-03, 5.292E-03,
01618 & 5.052E-03, 4.800E-03, 4.535E-03, 4.257E-03, 3.963E-03,
01619 & 3.653E-03, 3.620E-03, 3.588E-03, 3.556E-03, 3.523E-03,
01620 & 3.490E-03, 3.457E-03, 3.424E-03, 3.391E-03, 3.357E-03,
01621 & 3.323E-03, 3.289E-03, 3.255E-03, 3.221E-03, 3.186E-03,
01622 & 3.151E-03, 3.116E-03, 3.081E-03, 3.046E-03, 3.010E-03/
01623
01624 DATA ((DAT(I,J), I=301,400), J=1,1) /
01625 & 2.974E-03, 2.938E-03, 2.902E-03, 2.866E-03, 2.829E-03,
01626 & 2.792E-03, 2.755E-03, 2.718E-03, 2.681E-03, 2.643E-03,
01627 & 2.605E-03, 2.567E-03, 2.529E-03, 2.491E-03, 2.452E-03,
01628 & 2.413E-03, 2.375E-03, 2.335E-03, 2.296E-03, 2.257E-03,
01629 & 2.217E-03, 2.177E-03, 2.137E-03, 2.097E-03, 2.057E-03,
01630 & 2.017E-03, 1.977E-03, 1.936E-03, 1.895E-03, 1.855E-03,
01631 & 1.814E-03, 1.773E-03, 1.732E-03, 1.691E-03, 1.650E-03,
01632 & 1.609E-03, 1.567E-03, 1.526E-03, 1.485E-03, 1.444E-03,
01633 & 1.403E-03, 1.362E-03, 1.321E-03, 1.280E-03, 1.239E-03,
01634 & 1.199E-03, 1.158E-03, 1.118E-03, 1.078E-03, 1.038E-03,
01635 & 9.988E-04, 9.596E-04, 9.207E-04, 8.821E-04, 8.438E-04,
01636 & 8.060E-04, 7.686E-04, 7.316E-04, 6.952E-04, 6.592E-04,
01637 & 6.238E-04, 5.891E-04, 5.549E-04, 5.214E-04, 4.886E-04,
01638 & 4.565E-04, 4.253E-04, 3.948E-04, 3.652E-04, 3.364E-04,
01639 & 3.086E-04, 2.818E-04, 2.560E-04, 2.312E-04, 2.075E-04,
01640 & 1.849E-04, 1.634E-04, 1.432E-04, 1.241E-04, 1.064E-04,
01641 & 8.985E-05, 8.827E-05, 8.670E-05, 8.515E-05, 8.361E-05,
01642 & 8.209E-05, 8.057E-05, 7.907E-05, 7.759E-05, 7.612E-05,
01643 & 7.466E-05, 7.321E-05, 7.178E-05, 7.036E-05, 6.896E-05,
01644 & 6.757E-05, 6.619E-05, 6.482E-05, 6.347E-05, 6.213E-05/
01645
01646 DATA ((DAT(I,J), I=401,500), J=1,1) /
01647 & 6.081E-05, 5.950E-05, 5.821E-05, 5.692E-05, 5.566E-05,
01648 & 5.440E-05, 5.316E-05, 5.193E-05, 5.072E-05, 4.952E-05,
01649 & 4.833E-05, 4.716E-05, 4.600E-05, 4.486E-05, 4.373E-05,
01650 & 4.261E-05, 4.151E-05, 4.042E-05, 3.935E-05, 3.829E-05,
01651 & 3.724E-05, 3.621E-05, 3.519E-05, 3.419E-05, 3.320E-05,
01652 & 3.223E-05, 3.127E-05, 3.032E-05, 2.939E-05, 2.847E-05,
01653 & 2.757E-05, 2.668E-05, 2.580E-05, 2.494E-05, 2.409E-05,
01654 & 2.326E-05, 2.244E-05, 2.164E-05, 2.085E-05, 2.008E-05,
01655 & 1.932E-05, 1.857E-05, 1.784E-05, 1.713E-05, 1.642E-05,
01656 & 1.574E-05, 1.506E-05, 1.441E-05, 1.376E-05, 1.313E-05,
01657 & 1.252E-05, 1.192E-05, 1.133E-05, 1.076E-05, 1.021E-05,
01658 & 9.664E-06, 9.138E-06, 8.626E-06, 8.128E-06, 7.645E-06,
01659 & 7.177E-06, 6.724E-06, 6.286E-06, 5.862E-06, 5.453E-06,
01660 & 5.058E-06, 4.679E-06, 4.314E-06, 3.964E-06, 3.628E-06,
01661 & 3.308E-06, 3.002E-06, 2.711E-06, 2.435E-06, 2.174E-06,
01662 & 1.927E-06, 1.695E-06, 1.479E-06, 1.277E-06, 1.089E-06,
01663 & 9.168E-07, 9.004E-07, 8.841E-07, 8.680E-07, 8.520E-07,
01664 & 8.362E-07, 8.205E-07, 8.049E-07, 7.896E-07, 7.743E-07,
01665 & 7.592E-07, 7.443E-07, 7.295E-07, 7.149E-07, 7.004E-07,
01666 & 6.860E-07, 6.718E-07, 6.578E-07, 6.439E-07, 6.301E-07/
01667
01668 DATA ((DAT(I,J), I=501,600), J=1,1) /
01669 & 6.165E-07, 6.031E-07, 5.897E-07, 5.766E-07, 5.636E-07,
01670 & 5.507E-07, 5.380E-07, 5.254E-07, 5.130E-07, 5.007E-07,
01671 & 4.886E-07, 4.766E-07, 4.648E-07, 4.531E-07, 4.416E-07,
01672 & 4.302E-07, 4.190E-07, 4.079E-07, 3.970E-07, 3.862E-07,
01673 & 3.756E-07, 3.651E-07, 3.547E-07, 3.445E-07, 3.345E-07,
01674 & 3.246E-07, 3.149E-07, 3.053E-07, 2.958E-07, 2.865E-07,
01675 & 2.774E-07, 2.684E-07, 2.595E-07, 2.508E-07, 2.422E-07,
01676 & 2.339E-07, 2.256E-07, 2.175E-07, 2.095E-07, 2.017E-07,
01677 & 1.940E-07, 1.865E-07, 1.791E-07, 1.719E-07, 1.648E-07,
01678 & 1.579E-07, 1.512E-07, 1.445E-07, 1.380E-07, 1.317E-07,
01679 & 1.255E-07, 1.195E-07, 1.136E-07, 1.079E-07, 1.023E-07,
01680 & 9.686E-08, 9.157E-08, 8.642E-08, 8.143E-08, 7.659E-08,
01681 & 7.189E-08, 6.734E-08, 6.295E-08, 5.870E-08, 5.459E-08,
01682 & 5.064E-08, 4.684E-08, 4.318E-08, 3.967E-08, 3.631E-08,
01683 & 3.310E-08, 3.004E-08, 2.713E-08, 2.436E-08, 2.175E-08,
01684 & 1.928E-08, 1.696E-08, 1.479E-08, 1.277E-08, 1.089E-08,
01685 & 9.170E-09, 9.006E-09, 8.843E-09, 8.681E-09, 8.522E-09,
01686 & 8.363E-09, 8.206E-09, 8.051E-09, 7.897E-09, 7.745E-09,
01687 & 7.594E-09, 7.444E-09, 7.296E-09, 7.150E-09, 7.005E-09,
01688 & 6.861E-09, 6.719E-09, 6.579E-09, 6.440E-09, 6.302E-09/
01689
01690 DATA ((DAT(I,J), I=601,700), J=1,1) /
01691 & 6.166E-09, 6.031E-09, 5.898E-09, 5.767E-09, 5.636E-09,
01692 & 5.508E-09, 5.381E-09, 5.255E-09, 5.131E-09, 5.008E-09,
01693 & 4.887E-09, 4.767E-09, 4.649E-09, 4.532E-09, 4.417E-09,
01694 & 4.303E-09, 4.190E-09, 4.080E-09, 3.970E-09, 3.862E-09,
01695 & 3.756E-09, 3.651E-09, 3.548E-09, 3.446E-09, 3.345E-09,
01696 & 3.246E-09, 3.149E-09, 3.053E-09, 2.958E-09, 2.865E-09,
01697 & 2.774E-09, 2.684E-09, 2.595E-09, 2.508E-09, 2.423E-09,
01698 & 2.339E-09, 2.256E-09, 2.175E-09, 2.095E-09, 2.017E-09,
01699 & 1.940E-09, 1.865E-09, 1.791E-09, 1.719E-09, 1.649E-09,
01700 & 1.579E-09, 1.512E-09, 1.445E-09, 1.381E-09, 1.317E-09,
01701 & 1.255E-09, 1.195E-09, 1.136E-09, 1.079E-09, 1.023E-09,
01702 & 9.686E-10, 9.157E-10, 8.643E-10, 8.143E-10, 7.659E-10,
01703 & 7.189E-10, 6.735E-10, 6.295E-10, 5.870E-10, 5.459E-10,
01704 & 5.064E-10, 4.684E-10, 4.318E-10, 3.967E-10, 3.631E-10,
01705 & 3.310E-10, 3.004E-10, 2.713E-10, 2.436E-10, 2.175E-10,
01706 & 1.928E-10, 1.696E-10, 1.479E-10, 1.277E-10, 1.090E-10,
01707 & 9.170E-11, 9.006E-11, 8.843E-11, 8.681E-11, 8.522E-11,
01708 & 8.363E-11, 8.206E-11, 8.051E-11, 7.897E-11, 7.745E-11,
01709 & 7.594E-11, 7.444E-11, 7.296E-11, 7.150E-11, 7.005E-11,
01710 & 6.861E-11, 6.719E-11, 6.579E-11, 6.440E-11, 6.302E-11/
01711
01712 DATA ((DAT(I,J), I=701,800), J=1,1) /
01713 & 6.166E-11, 6.031E-11, 5.898E-11, 5.767E-11, 5.636E-11,
01714 & 5.508E-11, 5.380E-11, 5.255E-11, 5.131E-11, 5.008E-11,
01715 & 4.887E-11, 4.767E-11, 4.649E-11, 4.532E-11, 4.417E-11,
01716 & 4.303E-11, 4.191E-11, 4.080E-11, 3.970E-11, 3.862E-11,
01717 & 3.756E-11, 3.651E-11, 3.548E-11, 3.446E-11, 3.345E-11,
01718 & 3.246E-11, 3.149E-11, 3.053E-11, 2.959E-11, 2.865E-11,
01719 & 2.774E-11, 2.684E-11, 2.595E-11, 2.508E-11, 2.423E-11,
01720 & 2.339E-11, 2.256E-11, 2.175E-11, 2.095E-11, 2.017E-11,
01721 & 1.940E-11, 1.865E-11, 1.791E-11, 1.719E-11, 1.649E-11,
01722 & 1.579E-11, 1.511E-11, 1.445E-11, 1.380E-11, 1.317E-11,
01723 & 1.255E-11, 1.195E-11, 1.136E-11, 1.079E-11, 1.023E-11,
01724 & 9.686E-12, 9.157E-12, 8.643E-12, 8.143E-12, 7.659E-12,
01725 & 7.189E-12, 6.735E-12, 6.295E-12, 5.870E-12, 5.459E-12,
01726 & 5.064E-12, 4.684E-12, 4.318E-12, 3.967E-12, 3.631E-12,
01727 & 3.310E-12, 3.004E-12, 2.713E-12, 2.436E-12, 2.175E-12,
01728 & 1.928E-12, 1.696E-12, 1.479E-12, 1.277E-12, 1.090E-12,
01729 & 9.170E-13, 9.006E-13, 8.843E-13, 8.682E-13, 8.522E-13,
01730 & 8.363E-13, 8.206E-13, 8.051E-13, 7.897E-13, 7.745E-13,
01731 & 7.594E-13, 7.444E-13, 7.296E-13, 7.150E-13, 7.005E-13,
01732 & 6.861E-13, 6.719E-13, 6.579E-13, 6.440E-13, 6.302E-13/
01733
01734 DATA ((DAT(I,J), I=801,882), J=1,1) /
01735 & 6.166E-13, 6.031E-13, 5.898E-13, 5.767E-13, 5.636E-13,
01736 & 5.508E-13, 5.381E-13, 5.255E-13, 5.131E-13, 5.008E-13,
01737 & 4.887E-13, 4.767E-13, 4.649E-13, 4.532E-13, 4.417E-13,
01738 & 4.303E-13, 4.190E-13, 4.080E-13, 3.970E-13, 3.862E-13,
01739 & 3.756E-13, 3.651E-13, 3.548E-13, 3.446E-13, 3.345E-13,
01740 & 3.246E-13, 3.149E-13, 3.053E-13, 2.958E-13, 2.865E-13,
01741 & 2.774E-13, 2.684E-13, 2.595E-13, 2.508E-13, 2.423E-13,
01742 & 2.339E-13, 2.256E-13, 2.175E-13, 2.095E-13, 2.017E-13,
01743 & 1.940E-13, 1.865E-13, 1.792E-13, 1.719E-13, 1.648E-13,
01744 & 1.579E-13, 1.511E-13, 1.445E-13, 1.380E-13, 1.317E-13,
01745 & 1.255E-13, 1.195E-13, 1.136E-13, 1.079E-13, 1.023E-13,
01746 & 9.686E-14, 9.157E-14, 8.643E-14, 8.143E-14, 7.659E-14,
01747 & 7.189E-14, 6.735E-14, 6.295E-14, 5.870E-14, 5.459E-14,
01748 & 5.064E-14, 4.684E-14, 4.318E-14, 3.967E-14, 3.631E-14,
01749 & 3.310E-14, 3.004E-14, 2.713E-14, 2.436E-14, 2.175E-14,
01750 & 1.928E-14, 1.696E-14, 1.479E-14, 1.277E-14, 1.090E-14,
01751 & 9.170E-15, 0.000E+00/
01752
01753 DATA ((DAT(I,J), I=1,100), J=2,2) /
01754 & 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04, 3.330E-04,
01755 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
01756 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
01757 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
01758 & 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04,
01759 & 3.330E-04, 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04,
01760 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
01761 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
01762 & 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04, 3.330E-04,
01763 & 3.320E-04, 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04,
01764 & 3.320E-04, 3.320E-04, 3.330E-04, 3.320E-04, 3.320E-04,
01765 & 3.320E-04, 3.320E-04, 3.320E-04, 3.330E-04, 3.320E-04,
01766 & 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04,
01767 & 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04, 3.320E-04,
01768 & 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04, 3.320E-04,
01769 & 3.330E-04, 3.330E-04, 3.320E-04, 3.320E-04, 3.320E-04,
01770 & 3.320E-04, 3.320E-04, 3.330E-04, 3.320E-04, 3.320E-04,
01771 & 3.320E-04, 3.320E-04, 3.320E-04, 3.320E-04, 3.310E-04,
01772 & 3.310E-04, 3.320E-04, 3.310E-04, 3.310E-04, 3.310E-04,
01773 & 3.310E-04, 3.310E-04, 3.300E-04, 3.310E-04, 3.310E-04/
01774
01775 DATA ((DAT(I,J), I=101,200), J=2,2) /
01776 & 3.310E-04, 3.300E-04, 3.310E-04, 3.310E-04, 3.310E-04,
01777 & 3.310E-04, 3.310E-04, 3.310E-04, 3.300E-04, 3.300E-04,
01778 & 3.310E-04, 3.300E-04, 3.300E-04, 3.300E-04, 3.300E-04,
01779 & 3.300E-04, 3.300E-04, 3.300E-04, 3.300E-04, 3.290E-04,
01780 & 3.300E-04, 3.300E-04, 3.300E-04, 3.290E-04, 3.300E-04,
01781 & 3.300E-04, 3.290E-04, 3.290E-04, 3.290E-04, 3.280E-04,
01782 & 3.290E-04, 3.290E-04, 3.280E-04, 3.290E-04, 3.280E-04,
01783 & 3.280E-04, 3.280E-04, 3.280E-04, 3.280E-04, 3.270E-04,
01784 & 3.270E-04, 3.270E-04, 3.260E-04, 3.260E-04, 3.260E-04,
01785 & 3.250E-04, 3.250E-04, 3.250E-04, 3.250E-04, 3.240E-04,
01786 & 3.240E-04, 3.240E-04, 3.230E-04, 3.230E-04, 3.220E-04,
01787 & 3.210E-04, 3.200E-04, 3.200E-04, 3.200E-04, 3.180E-04,
01788 & 3.170E-04, 3.170E-04, 3.160E-04, 3.140E-04, 3.130E-04,
01789 & 3.120E-04, 3.110E-04, 3.080E-04, 3.060E-04, 3.040E-04,
01790 & 3.010E-04, 2.980E-04, 2.960E-04, 2.920E-04, 2.880E-04,
01791 & 2.830E-04, 2.780E-04, 2.710E-04, 2.630E-04, 2.540E-04,
01792 & 2.440E-04, 2.480E-04, 2.470E-04, 2.450E-04, 2.440E-04,
01793 & 2.430E-04, 2.420E-04, 2.410E-04, 2.390E-04, 2.390E-04,
01794 & 2.370E-04, 2.360E-04, 2.350E-04, 2.330E-04, 2.310E-04,
01795 & 2.310E-04, 2.290E-04, 2.280E-04, 2.260E-04, 2.250E-04/
01796
01797 DATA ((DAT(I,J), I=201,300), J=2,2) /
01798 & 2.240E-04, 2.220E-04, 2.210E-04, 2.180E-04, 2.170E-04,
01799 & 2.160E-04, 2.140E-04, 2.120E-04, 2.110E-04, 2.090E-04,
01800 & 2.080E-04, 2.060E-04, 2.050E-04, 2.030E-04, 2.000E-04,
01801 & 1.990E-04, 1.970E-04, 1.960E-04, 1.940E-04, 1.920E-04,
01802 & 1.900E-04, 1.880E-04, 1.860E-04, 1.840E-04, 1.820E-04,
01803 & 1.800E-04, 1.780E-04, 1.760E-04, 1.740E-04, 1.720E-04,
01804 & 1.700E-04, 1.680E-04, 1.660E-04, 1.630E-04, 1.610E-04,
01805 & 1.590E-04, 1.570E-04, 1.540E-04, 1.520E-04, 1.500E-04,
01806 & 1.470E-04, 1.450E-04, 1.430E-04, 1.409E-04, 1.380E-04,
01807 & 1.355E-04, 1.332E-04, 1.308E-04, 1.283E-04, 1.259E-04,
01808 & 1.233E-04, 1.209E-04, 1.185E-04, 1.160E-04, 1.135E-04,
01809 & 1.110E-04, 1.085E-04, 1.060E-04, 1.035E-04, 1.009E-04,
01810 & 9.850E-05, 9.600E-05, 9.340E-05, 9.100E-05, 8.850E-05,
01811 & 8.600E-05, 8.350E-05, 8.100E-05, 7.840E-05, 7.590E-05,
01812 & 7.340E-05, 7.090E-05, 6.820E-05, 6.560E-05, 6.290E-05,
01813 & 6.010E-05, 5.730E-05, 5.430E-05, 5.120E-05, 4.800E-05,
01814 & 4.470E-05, 4.430E-05, 4.400E-05, 4.360E-05, 4.320E-05,
01815 & 4.280E-05, 4.250E-05, 4.210E-05, 4.180E-05, 4.140E-05,
01816 & 4.100E-05, 4.070E-05, 4.030E-05, 3.990E-05, 3.950E-05,
01817 & 3.910E-05, 3.870E-05, 3.830E-05, 3.800E-05, 3.760E-05/
01818
01819 DATA ((DAT(I,J), I=301,400), J=2,2) /
01820 & 3.720E-05, 3.680E-05, 3.640E-05, 3.600E-05, 3.560E-05,
01821 & 3.520E-05, 3.470E-05, 3.430E-05, 3.390E-05, 3.350E-05,
01822 & 3.310E-05, 3.270E-05, 3.220E-05, 3.180E-05, 3.140E-05,
01823 & 3.100E-05, 3.050E-05, 3.010E-05, 2.960E-05, 2.920E-05,
01824 & 2.870E-05, 2.820E-05, 2.780E-05, 2.730E-05, 2.690E-05,
01825 & 2.640E-05, 2.600E-05, 2.540E-05, 2.490E-05, 2.450E-05,
01826 & 2.400E-05, 2.350E-05, 2.300E-05, 2.260E-05, 2.210E-05,
01827 & 2.160E-05, 2.110E-05, 2.060E-05, 2.000E-05, 1.950E-05,
01828 & 1.910E-05, 1.860E-05, 1.800E-05, 1.750E-05, 1.700E-05,
01829 & 1.650E-05, 1.600E-05, 1.550E-05, 1.500E-05, 1.450E-05,
01830 & 1.397E-05, 1.346E-05, 1.296E-05, 1.244E-05, 1.195E-05,
01831 & 1.145E-05, 1.095E-05, 1.045E-05, 9.960E-06, 9.470E-06,
01832 & 8.990E-06, 8.520E-06, 8.050E-06, 7.580E-06, 7.130E-06,
01833 & 6.680E-06, 6.230E-06, 5.800E-06, 5.380E-06, 4.980E-06,
01834 & 4.570E-06, 4.190E-06, 3.820E-06, 3.450E-06, 3.110E-06,
01835 & 2.770E-06, 2.460E-06, 2.160E-06, 1.870E-06, 1.610E-06,
01836 & 1.359E-06, 1.336E-06, 1.312E-06, 1.289E-06, 1.266E-06,
01837 & 1.243E-06, 1.221E-06, 1.198E-06, 1.175E-06, 1.153E-06,
01838 & 1.131E-06, 1.110E-06, 1.088E-06, 1.067E-06, 1.045E-06,
01839 & 1.025E-06, 1.004E-06, 9.830E-07, 9.620E-07, 9.420E-07/
01840
01841 DATA ((DAT(I,J), I=401,500), J=2,2) /
01842 & 9.230E-07, 9.030E-07, 8.840E-07, 8.640E-07, 8.450E-07,
01843 & 8.260E-07, 8.080E-07, 7.890E-07, 7.710E-07, 7.520E-07,
01844 & 7.340E-07, 7.170E-07, 7.000E-07, 6.820E-07, 6.650E-07,
01845 & 6.480E-07, 6.320E-07, 6.150E-07, 5.990E-07, 5.820E-07,
01846 & 5.660E-07, 5.510E-07, 5.350E-07, 5.200E-07, 5.060E-07,
01847 & 4.910E-07, 4.760E-07, 4.620E-07, 4.470E-07, 4.330E-07,
01848 & 4.200E-07, 4.060E-07, 3.930E-07, 3.800E-07, 3.670E-07,
01849 & 3.540E-07, 3.420E-07, 3.290E-07, 3.180E-07, 3.060E-07,
01850 & 2.950E-07, 2.830E-07, 2.720E-07, 2.610E-07, 2.500E-07,
01851 & 2.400E-07, 2.300E-07, 2.190E-07, 2.100E-07, 2.010E-07,
01852 & 1.910E-07, 1.820E-07, 1.730E-07, 1.640E-07, 1.560E-07,
01853 & 1.474E-07, 1.393E-07, 1.316E-07, 1.240E-07, 1.166E-07,
01854 & 1.095E-07, 1.026E-07, 9.590E-08, 8.950E-08, 8.320E-08,
01855 & 7.720E-08, 7.140E-08, 6.580E-08, 6.050E-08, 5.540E-08,
01856 & 5.050E-08, 4.580E-08, 4.140E-08, 3.710E-08, 3.320E-08,
01857 & 2.940E-08, 2.590E-08, 2.250E-08, 1.950E-08, 1.660E-08,
01858 & 1.400E-08, 1.374E-08, 1.350E-08, 1.325E-08, 1.301E-08,
01859 & 1.277E-08, 1.252E-08, 1.229E-08, 1.205E-08, 1.182E-08,
01860 & 1.160E-08, 1.136E-08, 1.114E-08, 1.092E-08, 1.069E-08,
01861 & 1.047E-08, 1.026E-08, 1.004E-08, 9.830E-09, 9.620E-09/
01862
01863 DATA ((DAT(I,J), I=501,600), J=2,2) /
01864 & 9.410E-09, 9.210E-09, 9.010E-09, 8.800E-09, 8.610E-09,
01865 & 8.400E-09, 8.220E-09, 8.020E-09, 7.830E-09, 7.650E-09,
01866 & 7.460E-09, 7.280E-09, 7.100E-09, 6.920E-09, 6.750E-09,
01867 & 6.570E-09, 6.400E-09, 6.220E-09, 6.060E-09, 5.900E-09,
01868 & 5.730E-09, 5.570E-09, 5.410E-09, 5.260E-09, 5.110E-09,
01869 & 4.960E-09, 4.810E-09, 4.660E-09, 4.520E-09, 4.380E-09,
01870 & 4.230E-09, 4.100E-09, 3.960E-09, 3.830E-09, 3.690E-09,
01871 & 3.580E-09, 3.440E-09, 3.320E-09, 3.200E-09, 3.070E-09,
01872 & 2.960E-09, 2.850E-09, 2.730E-09, 2.620E-09, 2.520E-09,
01873 & 2.410E-09, 2.310E-09, 2.200E-09, 2.100E-09, 2.010E-09,
01874 & 1.910E-09, 1.820E-09, 1.740E-09, 1.650E-09, 1.560E-09,
01875 & 1.479E-09, 1.398E-09, 1.320E-09, 1.244E-09, 1.169E-09,
01876 & 1.098E-09, 1.028E-09, 9.610E-10, 8.960E-10, 8.340E-10,
01877 & 7.740E-10, 7.150E-10, 6.590E-10, 6.060E-10, 5.550E-10,
01878 & 5.050E-10, 4.590E-10, 4.140E-10, 3.720E-10, 3.320E-10,
01879 & 2.950E-10, 2.590E-10, 2.260E-10, 1.950E-10, 1.670E-10,
01880 & 1.400E-10, 1.375E-10, 1.350E-10, 1.326E-10, 1.301E-10,
01881 & 1.276E-10, 1.253E-10, 1.229E-10, 1.205E-10, 1.182E-10,
01882 & 1.160E-10, 1.137E-10, 1.114E-10, 1.091E-10, 1.069E-10,
01883 & 1.048E-10, 1.026E-10, 1.004E-10, 9.840E-11, 9.620E-11/
01884
01885 DATA ((DAT(I,J), I=601,700), J=2,2) /
01886 & 9.420E-11, 9.210E-11, 9.010E-11, 8.800E-11, 8.610E-11,
01887 & 8.410E-11, 8.210E-11, 8.020E-11, 7.830E-11, 7.640E-11,
01888 & 7.460E-11, 7.280E-11, 7.100E-11, 6.920E-11, 6.740E-11,
01889 & 6.570E-11, 6.400E-11, 6.230E-11, 6.060E-11, 5.900E-11,
01890 & 5.730E-11, 5.570E-11, 5.410E-11, 5.260E-11, 5.110E-11,
01891 & 4.960E-11, 4.800E-11, 4.670E-11, 4.520E-11, 4.370E-11,
01892 & 4.230E-11, 4.100E-11, 3.960E-11, 3.830E-11, 3.700E-11,
01893 & 3.570E-11, 3.450E-11, 3.320E-11, 3.200E-11, 3.080E-11,
01894 & 2.970E-11, 2.850E-11, 2.740E-11, 2.630E-11, 2.510E-11,
01895 & 2.410E-11, 2.300E-11, 2.210E-11, 2.110E-11, 2.010E-11,
01896 & 1.920E-11, 1.820E-11, 1.730E-11, 1.640E-11, 1.570E-11,
01897 & 1.479E-11, 1.398E-11, 1.319E-11, 1.243E-11, 1.170E-11,
01898 & 1.098E-11, 1.028E-11, 9.610E-12, 8.960E-12, 8.330E-12,
01899 & 7.730E-12, 7.150E-12, 6.600E-12, 6.060E-12, 5.540E-12,
01900 & 5.060E-12, 4.590E-12, 4.140E-12, 3.720E-12, 3.330E-12,
01901 & 2.940E-12, 2.590E-12, 2.260E-12, 1.950E-12, 1.660E-12,
01902 & 1.400E-12, 1.375E-12, 1.350E-12, 1.326E-12, 1.301E-12,
01903 & 1.277E-12, 1.253E-12, 1.230E-12, 1.205E-12, 1.182E-12,
01904 & 1.160E-12, 1.137E-12, 1.114E-12, 1.091E-12, 1.069E-12,
01905 & 1.047E-12, 1.026E-12, 1.004E-12, 9.840E-13, 9.620E-13/
01906
01907 DATA ((DAT(I,J), I=701,800), J=2,2) /
01908 & 9.410E-13, 9.210E-13, 9.010E-13, 8.800E-13, 8.610E-13,
01909 & 8.410E-13, 8.210E-13, 8.030E-13, 7.830E-13, 7.640E-13,
01910 & 7.460E-13, 7.270E-13, 7.100E-13, 6.920E-13, 6.740E-13,
01911 & 6.570E-13, 6.400E-13, 6.230E-13, 6.060E-13, 5.900E-13,
01912 & 5.730E-13, 5.570E-13, 5.410E-13, 5.260E-13, 5.110E-13,
01913 & 4.960E-13, 4.800E-13, 4.670E-13, 4.520E-13, 4.370E-13,
01914 & 4.230E-13, 4.100E-13, 3.960E-13, 3.830E-13, 3.700E-13,
01915 & 3.570E-13, 3.450E-13, 3.320E-13, 3.200E-13, 3.080E-13,
01916 & 2.970E-13, 2.850E-13, 2.740E-13, 2.630E-13, 2.510E-13,
01917 & 2.410E-13, 2.300E-13, 2.210E-13, 2.110E-13, 2.010E-13,
01918 & 1.920E-13, 1.820E-13, 1.730E-13, 1.640E-13, 1.570E-13,
01919 & 1.479E-13, 1.398E-13, 1.319E-13, 1.243E-13, 1.170E-13,
01920 & 1.098E-13, 1.028E-13, 9.610E-14, 8.960E-14, 8.330E-14,
01921 & 7.730E-14, 7.150E-14, 6.600E-14, 6.060E-14, 5.540E-14,
01922 & 5.060E-14, 4.590E-14, 4.140E-14, 3.720E-14, 3.330E-14,
01923 & 2.940E-14, 2.590E-14, 2.260E-14, 1.950E-14, 1.660E-14,
01924 & 1.400E-14, 1.375E-14, 1.350E-14, 1.326E-14, 1.301E-14,
01925 & 1.277E-14, 1.253E-14, 1.230E-14, 1.205E-14, 1.182E-14,
01926 & 1.160E-14, 1.137E-14, 1.114E-14, 1.091E-14, 1.069E-14,
01927 & 1.047E-14, 1.026E-14, 1.004E-14, 9.840E-15, 9.620E-15/
01928
01929 DATA ((DAT(I,J), I=801,882), J=2,2) /
01930 & 9.410E-15, 9.210E-15, 9.010E-15, 8.800E-15, 8.610E-15,
01931 & 8.410E-15, 8.210E-15, 8.030E-15, 7.830E-15, 7.640E-15,
01932 & 7.460E-15, 7.270E-15, 7.100E-15, 6.920E-15, 6.740E-15,
01933 & 6.570E-15, 6.400E-15, 6.230E-15, 6.060E-15, 5.900E-15,
01934 & 5.730E-15, 5.570E-15, 5.410E-15, 5.260E-15, 5.110E-15,
01935 & 4.960E-15, 4.800E-15, 4.670E-15, 4.520E-15, 4.370E-15,
01936 & 4.230E-15, 4.100E-15, 3.960E-15, 3.830E-15, 3.700E-15,
01937 & 3.570E-15, 3.450E-15, 3.320E-15, 3.200E-15, 3.080E-15,
01938 & 2.970E-15, 2.850E-15, 2.740E-15, 2.630E-15, 2.510E-15,
01939 & 2.410E-15, 2.300E-15, 2.210E-15, 2.110E-15, 2.010E-15,
01940 & 1.920E-15, 1.820E-15, 1.730E-15, 1.640E-15, 1.570E-15,
01941 & 1.479E-15, 1.398E-15, 1.319E-15, 1.243E-15, 1.170E-15,
01942 & 1.098E-15, 1.028E-15, 9.610E-16, 8.960E-16, 8.330E-16,
01943 & 7.730E-16, 7.150E-16, 6.600E-16, 6.060E-16, 5.540E-16,
01944 & 5.060E-16, 4.590E-16, 4.140E-16, 3.720E-16, 3.330E-16,
01945 & 2.940E-16, 2.590E-16, 2.260E-16, 1.950E-16, 1.660E-16,
01946 & 1.400E-16, 0.000E+00/
01947 C
01948 save
01949 RETURN
01950 END
01951
01952 subroutine dalhad_timelike
01953 * hadronic contribution of the 5 light quark flavors evaluated using
01954 * e^+e^-data. CMD-2 2001 data have a normalization problem. I added a energy
01955 * independent +1.5 % correction and incresed the syst error by 0.5%.
01956 * Results thus are preliminary! Arrays:
01957 * ESA(I) : energy E in GeV;
01958 * DAS(I,1): Delta alpha^5(E**2)
01959 * DAS(I,2): error (systematics dominated, should be treated as a systematic error)
01960 C
01961 IMPLICIT NONE
01962 INTEGER NB,I,J
01963 PARAMETER(NB=530)
01964 REAL ESA(NB),DAS(NB,2)
01965 COMMON /DATS/ESA,DAS
01966 C 530
01967 DATA (ESA(I), I=1,100) /
01968 & 1.000E-06, 3.001E-03, 6.001E-03, 9.001E-03, 1.200E-02,
01969 & 1.500E-02, 1.800E-02, 2.100E-02, 2.400E-02, 2.700E-02,
01970 & 3.000E-02, 3.300E-02, 3.600E-02, 3.900E-02, 4.200E-02,
01971 & 4.500E-02, 4.800E-02, 5.100E-02, 5.400E-02, 5.700E-02,
01972 & 6.000E-02, 6.300E-02, 6.600E-02, 6.900E-02, 7.200E-02,
01973 & 7.500E-02, 7.800E-02, 8.100E-02, 8.400E-02, 8.700E-02,
01974 & 9.000E-02, 9.300E-02, 9.600E-02, 9.900E-02, 1.020E-01,
01975 & 1.050E-01, 1.080E-01, 1.110E-01, 1.140E-01, 1.170E-01,
01976 & 1.200E-01, 1.230E-01, 1.260E-01, 1.290E-01, 1.320E-01,
01977 & 1.350E-01, 1.380E-01, 1.410E-01, 1.440E-01, 1.470E-01,
01978 & 1.500E-01, 1.530E-01, 1.560E-01, 1.590E-01, 1.620E-01,
01979 & 1.650E-01, 1.680E-01, 1.710E-01, 1.740E-01, 1.770E-01,
01980 & 1.800E-01, 1.830E-01, 1.860E-01, 1.890E-01, 1.920E-01,
01981 & 1.950E-01, 1.980E-01, 2.010E-01, 2.040E-01, 2.070E-01,
01982 & 2.100E-01, 2.130E-01, 2.160E-01, 2.190E-01, 2.220E-01,
01983 & 2.250E-01, 2.280E-01, 2.310E-01, 2.340E-01, 2.370E-01,
01984 & 2.400E-01, 2.430E-01, 2.460E-01, 2.490E-01, 2.520E-01,
01985 & 2.550E-01, 2.580E-01, 2.610E-01, 2.640E-01, 2.670E-01,
01986 & 2.700E-01, 2.710E-01, 2.730E-01, 2.737E-01, 2.760E-01,
01987 & 2.763E-01, 2.770E-01, 3.200E-01, 3.387E-01, 3.573E-01/
01988
01989 DATA (ESA(I), I=101,200) /
01990 & 3.760E-01, 3.947E-01, 4.133E-01, 4.320E-01, 4.507E-01,
01991 & 4.693E-01, 4.880E-01, 5.067E-01, 5.253E-01, 5.440E-01,
01992 & 5.627E-01, 5.813E-01, 6.000E-01, 6.100E-01, 6.150E-01,
01993 & 6.200E-01, 6.250E-01, 6.300E-01, 6.350E-01, 6.400E-01,
01994 & 6.450E-01, 6.500E-01, 6.550E-01, 6.600E-01, 6.650E-01,
01995 & 6.700E-01, 6.750E-01, 6.800E-01, 6.850E-01, 6.900E-01,
01996 & 6.950E-01, 7.000E-01, 7.050E-01, 7.100E-01, 7.150E-01,
01997 & 7.200E-01, 7.250E-01, 7.300E-01, 7.350E-01, 7.400E-01,
01998 & 7.450E-01, 7.500E-01, 7.550E-01, 7.600E-01, 7.650E-01,
01999 & 7.700E-01, 7.704E-01, 7.708E-01, 7.712E-01, 7.716E-01,
02000 & 7.720E-01, 7.724E-01, 7.728E-01, 7.732E-01, 7.736E-01,
02001 & 7.740E-01, 7.744E-01, 7.748E-01, 7.752E-01, 7.756E-01,
02002 & 7.760E-01, 7.764E-01, 7.768E-01, 7.772E-01, 7.776E-01,
02003 & 7.780E-01, 7.784E-01, 7.788E-01, 7.792E-01, 7.796E-01,
02004 & 7.800E-01, 7.804E-01, 7.808E-01, 7.812E-01, 7.816E-01,
02005 & 7.820E-01, 7.824E-01, 7.828E-01, 7.832E-01, 7.836E-01,
02006 & 7.840E-01, 7.844E-01, 7.848E-01, 7.852E-01, 7.856E-01,
02007 & 7.860E-01, 7.864E-01, 7.868E-01, 7.872E-01, 7.876E-01,
02008 & 7.880E-01, 7.884E-01, 7.888E-01, 7.892E-01, 7.896E-01,
02009 & 7.900E-01, 7.905E-01, 7.910E-01, 7.915E-01, 7.920E-01/
02010
02011 DATA (ESA(I), I=201,300) /
02012 & 7.925E-01, 7.930E-01, 7.935E-01, 7.940E-01, 7.945E-01,
02013 & 7.950E-01, 7.955E-01, 7.960E-01, 7.965E-01, 7.970E-01,
02014 & 7.975E-01, 7.980E-01, 7.985E-01, 7.990E-01, 7.995E-01,
02015 & 8.000E-01, 8.005E-01, 8.010E-01, 8.015E-01, 8.020E-01,
02016 & 8.025E-01, 8.030E-01, 8.035E-01, 8.040E-01, 8.045E-01,
02017 & 8.050E-01, 8.055E-01, 8.060E-01, 8.065E-01, 8.070E-01,
02018 & 8.075E-01, 8.080E-01, 8.085E-01, 8.090E-01, 8.095E-01,
02019 & 8.100E-01, 8.105E-01, 8.110E-01, 8.115E-01, 8.120E-01,
02020 & 8.125E-01, 8.130E-01, 8.135E-01, 8.140E-01, 8.145E-01,
02021 & 8.150E-01, 8.155E-01, 8.160E-01, 8.165E-01, 8.170E-01,
02022 & 8.175E-01, 8.180E-01, 8.185E-01, 8.190E-01, 8.195E-01,
02023 & 8.200E-01, 8.220E-01, 8.240E-01, 8.260E-01, 8.280E-01,
02024 & 8.300E-01, 8.320E-01, 8.340E-01, 8.360E-01, 8.380E-01,
02025 & 8.400E-01, 8.420E-01, 8.440E-01, 8.460E-01, 8.480E-01,
02026 & 8.500E-01, 8.520E-01, 8.540E-01, 8.560E-01, 8.580E-01,
02027 & 8.600E-01, 8.620E-01, 8.640E-01, 8.660E-01, 8.680E-01,
02028 & 8.700E-01, 8.720E-01, 8.740E-01, 8.760E-01, 8.780E-01,
02029 & 8.800E-01, 8.820E-01, 8.840E-01, 8.860E-01, 8.880E-01,
02030 & 8.900E-01, 8.920E-01, 8.940E-01, 8.960E-01, 8.980E-01,
02031 & 9.000E-01, 9.020E-01, 9.040E-01, 9.060E-01, 9.080E-01/
02032
02033 DATA (ESA(I), I=301,400) /
02034 & 9.100E-01, 9.120E-01, 9.140E-01, 9.160E-01, 9.180E-01,
02035 & 9.200E-01, 9.240E-01, 9.280E-01, 9.320E-01, 9.360E-01,
02036 & 9.400E-01, 9.440E-01, 9.480E-01, 9.520E-01, 9.560E-01,
02037 & 9.600E-01, 9.640E-01, 9.680E-01, 9.720E-01, 9.760E-01,
02038 & 9.800E-01, 9.840E-01, 9.880E-01, 9.920E-01, 9.960E-01,
02039 & 1.004E+00, 1.008E+00, 1.012E+00, 1.016E+00, 1.020E+00,
02040 & 1.024E+00, 1.028E+00, 1.032E+00, 1.036E+00, 1.044E+00,
02041 & 1.048E+00, 1.052E+00, 1.056E+00, 1.060E+00, 1.064E+00,
02042 & 1.068E+00, 1.072E+00, 1.076E+00, 1.080E+00, 1.084E+00,
02043 & 1.088E+00, 1.092E+00, 1.096E+00, 1.100E+00, 2.000E+01,
02044 & 2.100E+01, 2.200E+01, 2.300E+01, 2.400E+01, 2.500E+01,
02045 & 2.600E+01, 2.700E+01, 2.800E+01, 2.900E+01, 3.000E+01,
02046 & 3.100E+01, 3.200E+01, 3.300E+01, 3.400E+01, 3.500E+01,
02047 & 3.600E+01, 3.700E+01, 3.800E+01, 3.900E+01, 4.000E+01,
02048 & 4.100E+01, 4.200E+01, 4.300E+01, 4.400E+01, 4.500E+01,
02049 & 4.600E+01, 4.700E+01, 4.800E+01, 4.900E+01, 5.000E+01,
02050 & 5.100E+01, 5.200E+01, 5.300E+01, 5.400E+01, 5.500E+01,
02051 & 5.600E+01, 5.700E+01, 5.800E+01, 5.900E+01, 6.000E+01,
02052 & 6.100E+01, 6.200E+01, 6.300E+01, 6.400E+01, 6.500E+01,
02053 & 6.600E+01, 6.700E+01, 6.800E+01, 6.900E+01, 7.000E+01/
02054
02055 DATA (ESA(I), I=401,500) /
02056 & 7.100E+01, 7.200E+01, 7.300E+01, 7.400E+01, 7.500E+01,
02057 & 7.600E+01, 7.700E+01, 7.800E+01, 7.900E+01, 8.000E+01,
02058 & 8.100E+01, 8.200E+01, 8.300E+01, 8.400E+01, 8.500E+01,
02059 & 8.600E+01, 8.700E+01, 8.800E+01, 8.900E+01, 9.000E+01,
02060 & 9.100E+01, 9.200E+01, 9.300E+01, 9.400E+01, 9.500E+01,
02061 & 9.600E+01, 9.700E+01, 9.800E+01, 9.900E+01, 1.000E+02,
02062 & 1.090E+02, 1.180E+02, 1.270E+02, 1.360E+02, 1.450E+02,
02063 & 1.540E+02, 1.630E+02, 1.720E+02, 1.810E+02, 1.900E+02,
02064 & 1.990E+02, 2.080E+02, 2.170E+02, 2.260E+02, 2.350E+02,
02065 & 2.440E+02, 2.530E+02, 2.620E+02, 2.710E+02, 2.800E+02,
02066 & 2.890E+02, 2.980E+02, 3.070E+02, 3.160E+02, 3.250E+02,
02067 & 3.340E+02, 3.430E+02, 3.520E+02, 3.610E+02, 3.700E+02,
02068 & 3.790E+02, 3.880E+02, 3.970E+02, 4.060E+02, 4.150E+02,
02069 & 4.240E+02, 4.330E+02, 4.420E+02, 4.510E+02, 4.600E+02,
02070 & 4.690E+02, 4.780E+02, 4.870E+02, 4.960E+02, 5.050E+02,
02071 & 5.140E+02, 5.230E+02, 5.320E+02, 5.410E+02, 5.500E+02,
02072 & 5.590E+02, 5.680E+02, 5.770E+02, 5.860E+02, 5.950E+02,
02073 & 6.040E+02, 6.130E+02, 6.220E+02, 6.310E+02, 6.400E+02,
02074 & 6.490E+02, 6.580E+02, 6.670E+02, 6.760E+02, 6.850E+02,
02075 & 6.940E+02, 7.030E+02, 7.120E+02, 7.210E+02, 7.300E+02/
02076
02077 DATA (ESA(I), I=501,530) /
02078 & 7.390E+02, 7.480E+02, 7.570E+02, 7.660E+02, 7.750E+02,
02079 & 7.840E+02, 7.930E+02, 8.020E+02, 8.110E+02, 8.200E+02,
02080 & 8.290E+02, 8.380E+02, 8.470E+02, 8.560E+02, 8.650E+02,
02081 & 8.740E+02, 8.830E+02, 8.920E+02, 9.010E+02, 9.100E+02,
02082 & 9.190E+02, 9.280E+02, 9.370E+02, 9.460E+02, 9.550E+02,
02083 & 9.640E+02, 9.730E+02, 9.820E+02, 9.910E+02, 1.000E+03/
02084
02085 DATA ((DAS(I,J), I=1,100), J=1,1) /
02086 & -9.736E-15,-8.259E-08,-3.303E-07,-7.431E-07,-1.321E-06,
02087 & -2.064E-06,-2.973E-06,-4.048E-06,-5.289E-06,-6.696E-06,
02088 & -8.269E-06,-1.001E-05,-1.192E-05,-1.399E-05,-1.624E-05,
02089 & -1.865E-05,-2.123E-05,-2.398E-05,-2.690E-05,-3.000E-05,
02090 & -3.326E-05,-3.670E-05,-4.031E-05,-4.410E-05,-4.806E-05,
02091 & -5.220E-05,-5.652E-05,-6.101E-05,-6.568E-05,-7.054E-05,
02092 & -7.557E-05,-8.079E-05,-8.619E-05,-9.178E-05,-9.755E-05,
02093 & -1.035E-04,-1.097E-04,-1.160E-04,-1.226E-04,-1.293E-04,
02094 & -1.362E-04,-1.434E-04,-1.507E-04,-1.582E-04,-1.660E-04,
02095 & -1.739E-04,-1.821E-04,-1.904E-04,-1.990E-04,-2.078E-04,
02096 & -2.168E-04,-2.261E-04,-2.356E-04,-2.453E-04,-2.552E-04,
02097 & -2.653E-04,-2.758E-04,-2.864E-04,-2.973E-04,-3.084E-04,
02098 & -3.198E-04,-3.315E-04,-3.434E-04,-3.556E-04,-3.680E-04,
02099 & -3.808E-04,-3.938E-04,-4.071E-04,-4.207E-04,-4.346E-04,
02100 & -4.488E-04,-4.634E-04,-4.782E-04,-4.934E-04,-5.089E-04,
02101 & -5.248E-04,-5.410E-04,-5.577E-04,-5.747E-04,-5.921E-04,
02102 & -6.099E-04,-6.281E-04,-6.468E-04,-6.660E-04,-6.856E-04,
02103 & -7.058E-04,-7.266E-04,-7.479E-04,-7.699E-04,-7.926E-04,
02104 & -8.161E-04,-8.089E-04,-8.406E-04,-8.306E-04,-8.662E-04,
02105 & -8.533E-04,-8.751E-04,-1.262E-03,-1.438E-03,-1.612E-03/
02106
02107 DATA ((DAS(I,J), I=101,200), J=1,1) /
02108 & -1.812E-03,-1.991E-03,-2.223E-03,-2.460E-03,-2.697E-03,
02109 & -3.007E-03,-3.356E-03,-3.682E-03,-4.029E-03,-4.532E-03,
02110 & -5.050E-03,-5.545E-03,-6.280E-03,-6.663E-03,-6.863E-03,
02111 & -7.069E-03,-7.281E-03,-7.497E-03,-7.719E-03,-7.944E-03,
02112 & -8.173E-03,-8.404E-03,-8.636E-03,-8.866E-03,-9.093E-03,
02113 & -9.313E-03,-9.522E-03,-9.715E-03,-9.887E-03,-1.003E-02,
02114 & -1.014E-02,-1.020E-02,-1.020E-02,-1.013E-02,-9.987E-03,
02115 & -9.745E-03,-9.399E-03,-8.941E-03,-8.372E-03,-7.702E-03,
02116 & -6.954E-03,-6.172E-03,-5.429E-03,-4.842E-03,-4.600E-03,
02117 & -4.939E-03,-5.012E-03,-5.092E-03,-5.180E-03,-5.274E-03,
02118 & -5.376E-03,-5.485E-03,-5.597E-03,-5.720E-03,-5.846E-03,
02119 & -5.974E-03,-6.104E-03,-6.233E-03,-6.356E-03,-6.470E-03,
02120 & -6.567E-03,-6.640E-03,-6.677E-03,-6.665E-03,-6.584E-03,
02121 & -6.415E-03,-6.127E-03,-5.688E-03,-5.059E-03,-4.197E-03,
02122 & -3.059E-03,-1.607E-03, 1.835E-04, 2.313E-03, 4.745E-03,
02123 & 7.403E-03, 1.017E-02, 1.291E-02, 1.549E-02, 1.777E-02,
02124 & 1.968E-02, 2.120E-02, 2.232E-02, 2.308E-02, 2.354E-02,
02125 & 2.375E-02, 2.377E-02, 2.365E-02, 2.342E-02, 2.313E-02,
02126 & 2.279E-02, 2.242E-02, 2.203E-02, 2.165E-02, 2.127E-02,
02127 & 2.088E-02, 2.035E-02, 1.983E-02, 1.936E-02, 1.888E-02/
02128
02129 DATA ((DAS(I,J), I=201,300), J=1,1) /
02130 & 1.841E-02, 1.798E-02, 1.756E-02, 1.718E-02, 1.679E-02,
02131 & 1.642E-02, 1.609E-02, 1.575E-02, 1.546E-02, 1.515E-02,
02132 & 1.485E-02, 1.458E-02, 1.431E-02, 1.408E-02, 1.383E-02,
02133 & 1.359E-02, 1.338E-02, 1.316E-02, 1.297E-02, 1.277E-02,
02134 & 1.258E-02, 1.239E-02, 1.223E-02, 1.206E-02, 1.188E-02,
02135 & 1.173E-02, 1.157E-02, 1.144E-02, 1.128E-02, 1.115E-02,
02136 & 1.101E-02, 1.089E-02, 1.075E-02, 1.062E-02, 1.051E-02,
02137 & 1.038E-02, 1.029E-02, 1.017E-02, 1.006E-02, 9.953E-03,
02138 & 9.843E-03, 9.760E-03, 9.655E-03, 9.565E-03, 9.465E-03,
02139 & 9.391E-03, 9.295E-03, 9.214E-03, 9.123E-03, 9.033E-03,
02140 & 8.968E-03, 8.883E-03, 8.810E-03, 8.728E-03, 8.669E-03,
02141 & 8.590E-03, 8.320E-03, 8.070E-03, 7.818E-03, 7.602E-03,
02142 & 7.400E-03, 7.210E-03, 7.022E-03, 6.852E-03, 6.692E-03,
02143 & 6.540E-03, 6.379E-03, 6.242E-03, 6.110E-03, 5.983E-03,
02144 & 5.862E-03, 5.738E-03, 5.625E-03, 5.510E-03, 5.581E-03,
02145 & 5.610E-03, 5.586E-03, 5.661E-03, 5.684E-03, 5.705E-03,
02146 & 5.683E-03, 5.688E-03, 5.646E-03, 5.651E-03, 5.621E-03,
02147 & 5.579E-03, 5.588E-03, 5.622E-03, 5.555E-03, 5.545E-03,
02148 & 5.473E-03, 5.490E-03, 5.456E-03, 5.465E-03, 5.382E-03,
02149 & 5.342E-03, 5.337E-03, 5.254E-03, 5.251E-03, 5.205E-03/
02150
02151 DATA ((DAS(I,J), I=301,400), J=1,1) /
02152 & 5.185E-03, 5.102E-03, 5.086E-03, 5.028E-03, 4.931E-03,
02153 & 4.902E-03, 4.894E-03, 4.725E-03, 4.550E-03, 4.365E-03,
02154 & 4.186E-03, 3.968E-03, 3.764E-03, 3.498E-03, 3.227E-03,
02155 & 2.904E-03, 2.587E-03, 2.163E-03, 1.705E-03, 1.144E-03,
02156 & 5.114E-04,-2.812E-04,-1.269E-03,-2.578E-03,-4.424E-03,
02157 & -1.046E-02,-1.577E-02,-2.603E-02,-4.838E-02, 3.453E-02,
02158 & 5.283E-02, 3.471E-02, 2.630E-02, 2.175E-02, 1.621E-02,
02159 & 1.446E-02, 1.324E-02, 1.230E-02, 1.154E-02, 1.092E-02,
02160 & 1.040E-02, 9.966E-03, 9.588E-03, 9.242E-03, 8.955E-03,
02161 & 8.691E-03, 8.450E-03, 8.227E-03, 8.027E-03, 1.862E-02,
02162 & 1.892E-02, 1.920E-02, 1.947E-02, 1.973E-02, 1.998E-02,
02163 & 2.021E-02, 2.044E-02, 2.066E-02, 2.087E-02, 2.107E-02,
02164 & 2.127E-02, 2.146E-02, 2.164E-02, 2.182E-02, 2.200E-02,
02165 & 2.216E-02, 2.232E-02, 2.248E-02, 2.264E-02, 2.279E-02,
02166 & 2.294E-02, 2.308E-02, 2.322E-02, 2.335E-02, 2.349E-02,
02167 & 2.362E-02, 2.374E-02, 2.387E-02, 2.399E-02, 2.411E-02,
02168 & 2.423E-02, 2.434E-02, 2.445E-02, 2.456E-02, 2.467E-02,
02169 & 2.478E-02, 2.488E-02, 2.498E-02, 2.508E-02, 2.518E-02,
02170 & 2.528E-02, 2.537E-02, 2.547E-02, 2.556E-02, 2.565E-02,
02171 & 2.574E-02, 2.583E-02, 2.592E-02, 2.600E-02, 2.609E-02/
02172
02173 DATA ((DAS(I,J), I=401,500), J=1,1) /
02174 & 2.617E-02, 2.625E-02, 2.633E-02, 2.641E-02, 2.649E-02,
02175 & 2.657E-02, 2.664E-02, 2.672E-02, 2.679E-02, 2.686E-02,
02176 & 2.693E-02, 2.701E-02, 2.708E-02, 2.715E-02, 2.722E-02,
02177 & 2.728E-02, 2.735E-02, 2.742E-02, 2.748E-02, 2.755E-02,
02178 & 2.761E-02, 2.767E-02, 2.774E-02, 2.780E-02, 2.786E-02,
02179 & 2.792E-02, 2.798E-02, 2.804E-02, 2.809E-02, 2.815E-02,
02180 & 2.865E-02, 2.910E-02, 2.951E-02, 2.990E-02, 3.025E-02,
02181 & 3.059E-02, 3.090E-02, 3.118E-02, 3.146E-02, 3.171E-02,
02182 & 3.195E-02, 3.217E-02, 3.238E-02, 3.258E-02, 3.276E-02,
02183 & 3.293E-02, 3.309E-02, 3.323E-02, 3.336E-02, 3.347E-02,
02184 & 3.357E-02, 3.366E-02, 3.372E-02, 3.376E-02, 3.373E-02,
02185 & 3.365E-02, 3.338E-02, 3.303E-02, 3.323E-02, 3.345E-02,
02186 & 3.386E-02, 3.415E-02, 3.444E-02, 3.470E-02, 3.499E-02,
02187 & 3.523E-02, 3.548E-02, 3.571E-02, 3.596E-02, 3.618E-02,
02188 & 3.638E-02, 3.660E-02, 3.679E-02, 3.699E-02, 3.718E-02,
02189 & 3.736E-02, 3.754E-02, 3.771E-02, 3.788E-02, 3.804E-02,
02190 & 3.821E-02, 3.836E-02, 3.852E-02, 3.867E-02, 3.882E-02,
02191 & 3.896E-02, 3.909E-02, 3.924E-02, 3.937E-02, 3.951E-02,
02192 & 3.964E-02, 3.976E-02, 3.989E-02, 4.001E-02, 4.013E-02,
02193 & 4.025E-02, 4.037E-02, 4.048E-02, 4.060E-02, 4.071E-02/
02194
02195 DATA ((DAS(I,J), I=501,530), J=1,1) /
02196 & 4.082E-02, 4.092E-02, 4.103E-02, 4.114E-02, 4.124E-02,
02197 & 4.134E-02, 4.145E-02, 4.154E-02, 4.164E-02, 4.174E-02,
02198 & 4.183E-02, 4.192E-02, 4.202E-02, 4.211E-02, 4.219E-02,
02199 & 4.229E-02, 4.238E-02, 4.247E-02, 4.255E-02, 4.264E-02,
02200 & 4.272E-02, 4.280E-02, 4.289E-02, 4.297E-02, 4.304E-02,
02201 & 4.313E-02, 4.321E-02, 4.328E-02, 4.336E-02, 4.344E-02/
02202
02203 DATA ((DAS(I,J), I=1,100), J=2,2) /
02204 & 1.427E-16, 1.261E-09, 5.040E-09, 1.135E-08, 2.020E-08,
02205 & 3.150E-08, 4.540E-08, 6.180E-08, 8.080E-08, 1.023E-07,
02206 & 1.264E-07, 1.530E-07, 1.820E-07, 2.140E-07, 2.490E-07,
02207 & 2.850E-07, 3.250E-07, 3.670E-07, 4.120E-07, 4.590E-07,
02208 & 5.100E-07, 5.630E-07, 6.180E-07, 6.770E-07, 7.370E-07,
02209 & 8.010E-07, 8.680E-07, 9.380E-07, 1.010E-06, 1.084E-06,
02210 & 1.163E-06, 1.244E-06, 1.327E-06, 1.414E-06, 1.505E-06,
02211 & 1.600E-06, 1.700E-06, 1.790E-06, 1.890E-06, 2.000E-06,
02212 & 2.100E-06, 2.220E-06, 2.330E-06, 2.460E-06, 2.580E-06,
02213 & 2.700E-06, 2.830E-06, 2.960E-06, 3.100E-06, 3.240E-06,
02214 & 3.390E-06, 3.540E-06, 3.690E-06, 3.840E-06, 4.000E-06,
02215 & 4.170E-06, 4.330E-06, 4.510E-06, 4.690E-06, 4.870E-06,
02216 & 5.050E-06, 5.250E-06, 5.440E-06, 5.640E-06, 5.850E-06,
02217 & 6.050E-06, 6.270E-06, 6.490E-06, 6.720E-06, 6.950E-06,
02218 & 7.200E-06, 7.440E-06, 7.690E-06, 7.950E-06, 8.210E-06,
02219 & 8.490E-06, 8.770E-06, 9.050E-06, 9.350E-06, 9.650E-06,
02220 & 9.970E-06, 1.029E-05, 1.061E-05, 1.096E-05, 1.131E-05,
02221 & 1.166E-05, 1.204E-05, 1.242E-05, 1.282E-05, 1.323E-05,
02222 & 1.366E-05, 1.670E-05, 1.410E-05, 1.715E-05, 1.456E-05,
02223 & 1.761E-05, 1.471E-05, 2.950E-05, 3.200E-05, 3.550E-05/
02224
02225 DATA ((DAS(I,J), I=101,200), J=2,2) /
02226 & 3.970E-05, 4.330E-05, 4.820E-05, 5.320E-05, 5.800E-05,
02227 & 6.470E-05, 7.230E-05, 7.900E-05, 8.600E-05, 9.730E-05,
02228 & 1.083E-04, 1.183E-04, 1.348E-04, 1.430E-04, 1.472E-04,
02229 & 1.516E-04, 1.560E-04, 1.605E-04, 1.650E-04, 1.697E-04,
02230 & 1.743E-04, 1.788E-04, 1.833E-04, 1.875E-04, 1.915E-04,
02231 & 1.951E-04, 1.982E-04, 2.008E-04, 2.021E-04, 2.030E-04,
02232 & 2.030E-04, 2.000E-04, 1.960E-04, 1.890E-04, 1.797E-04,
02233 & 1.685E-04, 1.541E-04, 1.381E-04, 1.224E-04, 1.114E-04,
02234 & 1.121E-04, 1.303E-04, 1.656E-04, 2.154E-04, 2.797E-04,
02235 & 2.610E-04, 2.677E-04, 2.748E-04, 2.821E-04, 2.898E-04,
02236 & 2.977E-04, 3.060E-04, 3.145E-04, 3.235E-04, 3.328E-04,
02237 & 3.424E-04, 3.521E-04, 3.622E-04, 3.724E-04, 3.827E-04,
02238 & 3.929E-04, 4.027E-04, 4.120E-04, 4.204E-04, 4.274E-04,
02239 & 4.325E-04, 4.348E-04, 4.334E-04, 4.273E-04, 4.152E-04,
02240 & 3.958E-04, 3.682E-04, 3.319E-04, 2.878E-04, 2.391E-04,
02241 & 1.937E-04, 1.660E-04, 1.730E-04, 2.100E-04, 2.590E-04,
02242 & 3.090E-04, 3.520E-04, 3.880E-04, 4.150E-04, 4.330E-04,
02243 & 4.450E-04, 4.510E-04, 4.530E-04, 4.500E-04, 4.450E-04,
02244 & 4.380E-04, 4.300E-04, 4.220E-04, 4.130E-04, 4.030E-04,
02245 & 2.106E-03, 2.063E-03, 2.021E-03, 1.980E-03, 1.940E-03/
02246
02247 DATA ((DAS(I,J), I=201,300), J=2,2) /
02248 & 1.902E-03, 1.863E-03, 1.825E-03, 1.788E-03, 1.753E-03,
02249 & 1.719E-03, 1.685E-03, 1.651E-03, 1.618E-03, 1.586E-03,
02250 & 1.554E-03, 1.524E-03, 1.493E-03, 1.464E-03, 1.434E-03,
02251 & 1.406E-03, 1.377E-03, 1.350E-03, 1.322E-03, 1.296E-03,
02252 & 1.268E-03, 1.243E-03, 1.217E-03, 1.192E-03, 1.167E-03,
02253 & 1.142E-03, 1.118E-03, 1.094E-03, 1.071E-03, 1.048E-03,
02254 & 1.025E-03, 1.003E-03, 9.807E-04, 9.593E-04, 9.375E-04,
02255 & 9.167E-04, 8.949E-04, 8.743E-04, 8.540E-04, 8.337E-04,
02256 & 8.138E-04, 7.939E-04, 7.745E-04, 7.552E-04, 7.363E-04,
02257 & 7.173E-04, 6.988E-04, 6.804E-04, 6.623E-04, 6.444E-04,
02258 & 6.266E-04, 6.091E-04, 5.917E-04, 5.746E-04, 5.575E-04,
02259 & 5.408E-04, 4.756E-04, 4.133E-04, 3.543E-04, 2.983E-04,
02260 & 2.459E-04, 1.983E-04, 1.578E-04, 1.350E-04, 1.640E-04,
02261 & 2.013E-04, 2.425E-04, 2.851E-04, 3.281E-04, 3.708E-04,
02262 & 4.131E-04, 4.551E-04, 4.964E-04, 5.281E-04, 5.300E-04,
02263 & 5.320E-04, 5.332E-04, 5.356E-04, 5.370E-04, 5.382E-04,
02264 & 5.391E-04, 5.402E-04, 5.412E-04, 5.425E-04, 5.439E-04,
02265 & 5.445E-04, 5.455E-04, 5.462E-04, 5.465E-04, 5.472E-04,
02266 & 5.474E-04, 5.482E-04, 5.486E-04, 5.497E-04, 5.500E-04,
02267 & 5.508E-04, 5.516E-04, 5.516E-04, 5.524E-04, 5.532E-04/
02268
02269 DATA ((DAS(I,J), I=301,400), J=2,2) /
02270 & 5.548E-04, 5.545E-04, 5.550E-04, 5.555E-04, 5.557E-04,
02271 & 5.569E-04, 5.575E-04, 5.586E-04, 5.600E-04, 5.612E-04,
02272 & 5.633E-04, 5.654E-04, 5.674E-04, 5.682E-04, 5.692E-04,
02273 & 5.690E-04, 5.694E-04, 5.701E-04, 5.712E-04, 5.735E-04,
02274 & 5.774E-04, 5.824E-04, 5.908E-04, 6.027E-04, 6.230E-04,
02275 & 7.080E-04, 8.030E-04, 1.020E-03, 1.559E-03, 9.470E-04,
02276 & 1.380E-03, 9.520E-04, 7.760E-04, 6.930E-04, 6.110E-04,
02277 & 5.890E-04, 5.750E-04, 5.640E-04, 5.570E-04, 5.500E-04,
02278 & 5.456E-04, 5.426E-04, 5.401E-04, 5.376E-04, 5.355E-04,
02279 & 5.342E-04, 5.326E-04, 5.307E-04, 5.299E-04, 3.760E-04,
02280 & 3.700E-04, 3.670E-04, 3.640E-04, 3.600E-04, 3.580E-04,
02281 & 3.550E-04, 3.540E-04, 3.520E-04, 3.510E-04, 3.500E-04,
02282 & 3.480E-04, 3.470E-04, 3.470E-04, 3.460E-04, 3.450E-04,
02283 & 3.440E-04, 3.430E-04, 3.430E-04, 3.420E-04, 3.420E-04,
02284 & 3.420E-04, 3.410E-04, 3.400E-04, 3.400E-04, 3.400E-04,
02285 & 3.390E-04, 3.390E-04, 3.390E-04, 3.390E-04, 3.390E-04,
02286 & 3.390E-04, 3.380E-04, 3.380E-04, 3.370E-04, 3.370E-04,
02287 & 3.370E-04, 3.370E-04, 3.370E-04, 3.370E-04, 3.370E-04,
02288 & 3.370E-04, 3.360E-04, 3.360E-04, 3.360E-04, 3.360E-04,
02289 & 3.360E-04, 3.360E-04, 3.360E-04, 3.350E-04, 3.350E-04/
02290
02291 DATA ((DAS(I,J), I=401,500), J=2,2) /
02292 & 3.350E-04, 3.360E-04, 3.350E-04, 3.350E-04, 3.360E-04,
02293 & 3.350E-04, 3.350E-04, 3.350E-04, 3.350E-04, 3.350E-04,
02294 & 3.350E-04, 3.350E-04, 3.350E-04, 3.340E-04, 3.350E-04,
02295 & 3.350E-04, 3.350E-04, 3.350E-04, 3.340E-04, 3.350E-04,
02296 & 3.340E-04, 3.350E-04, 3.340E-04, 3.350E-04, 3.350E-04,
02297 & 3.340E-04, 3.350E-04, 3.350E-04, 3.340E-04, 3.340E-04,
02298 & 3.340E-04, 3.340E-04, 3.330E-04, 3.330E-04, 3.340E-04,
02299 & 3.340E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02300 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02301 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02302 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02303 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02304 & 3.320E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02305 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02306 & 3.330E-04, 3.320E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02307 & 3.330E-04, 3.320E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02308 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02309 & 3.330E-04, 3.320E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02310 & 3.330E-04, 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04,
02311 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04/
02312
02313 DATA ((DAS(I,J), I=501,530), J=2,2) /
02314 & 3.330E-04, 3.330E-04, 3.330E-04, 3.320E-04, 3.330E-04,
02315 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02316 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02317 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02318 & 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04, 3.330E-04,
02319 & 3.330E-04, 3.320E-04, 3.330E-04, 3.330E-04, 3.330E-04/
02320 C
02321 save
02322 RETURN
02323 END
02324
02325 subroutine dalhad_timelike1
02326 * hadronic contribution of the 5 light quark flavors evaluated using
02327 * e^+e^-data. CMD-2 2001 data have a normalization problem. I added a energy
02328 * independent +1.5 % correction and incresed the syst error by 0.5%.
02329 * Results thus are preliminary! Arrays:
02330 * ESA(I) : energy E in GeV;
02331 * DAS(I,1): Delta alpha^5(E**2)
02332 * DAS(I,2): error (systematics dominated, should be treated as a systematic error)
02333 C
02334 IMPLICIT NONE
02335 INTEGER NC,I,J
02336 PARAMETER(NC=250)
02337 REAL EMA(NC),DAM(NC,2)
02338 COMMON /DATM/EMA,DAM
02339 C 250
02340 DATA (EMA(I), I=1,100) /
02341 & 1.600E+00, 1.700E+00, 1.800E+00, 1.900E+00, 2.000E+00,
02342 & 2.100E+00, 2.200E+00, 2.300E+00, 2.400E+00, 2.500E+00,
02343 & 2.600E+00, 2.620E+00, 2.640E+00, 2.660E+00, 2.680E+00,
02344 & 2.700E+00, 2.720E+00, 2.740E+00, 2.760E+00, 2.780E+00,
02345 & 2.800E+00, 2.820E+00, 2.840E+00, 2.860E+00, 2.880E+00,
02346 & 2.900E+00, 2.920E+00, 2.940E+00, 2.960E+00, 2.980E+00,
02347 & 3.000E+00, 3.020E+00, 3.040E+00, 3.060E+00, 3.080E+00,
02348 & 3.100E+00, 3.108E+00, 3.115E+00, 3.122E+00, 3.130E+00,
02349 & 3.138E+00, 3.145E+00, 3.152E+00, 3.160E+00, 3.168E+00,
02350 & 3.175E+00, 3.182E+00, 3.190E+00, 3.197E+00, 3.205E+00,
02351 & 3.213E+00, 3.220E+00, 3.227E+00, 3.235E+00, 3.243E+00,
02352 & 3.250E+00, 3.257E+00, 3.265E+00, 3.273E+00, 3.280E+00,
02353 & 3.287E+00, 3.295E+00, 3.303E+00, 3.310E+00, 3.318E+00,
02354 & 3.325E+00, 3.332E+00, 3.340E+00, 3.348E+00, 3.355E+00,
02355 & 3.362E+00, 3.370E+00, 3.378E+00, 3.385E+00, 3.392E+00,
02356 & 3.400E+00, 3.410E+00, 3.420E+00, 3.430E+00, 3.440E+00,
02357 & 3.450E+00, 3.460E+00, 3.470E+00, 3.480E+00, 3.490E+00,
02358 & 3.500E+00, 3.510E+00, 3.520E+00, 3.530E+00, 3.540E+00,
02359 & 3.550E+00, 3.560E+00, 3.570E+00, 3.580E+00, 3.590E+00,
02360 & 3.620E+00, 3.626E+00, 3.631E+00, 3.637E+00, 3.642E+00/
02361
02362 DATA (EMA(I), I=101,200) /
02363 & 3.648E+00, 3.654E+00, 3.659E+00, 3.665E+00, 3.670E+00,
02364 & 3.676E+00, 3.682E+00, 3.687E+00, 3.690E+00, 3.693E+00,
02365 & 3.695E+00, 3.698E+00, 3.699E+00, 3.700E+00, 3.704E+00,
02366 & 3.704E+00, 3.708E+00, 3.710E+00, 3.710E+00, 3.712E+00,
02367 & 3.715E+00, 3.716E+00, 3.720E+00, 3.720E+00, 3.721E+00,
02368 & 3.725E+00, 3.726E+00, 3.729E+00, 3.730E+00, 3.732E+00,
02369 & 3.733E+00, 3.737E+00, 3.740E+00, 3.742E+00, 3.743E+00,
02370 & 3.746E+00, 3.749E+00, 3.750E+00, 3.750E+00, 3.754E+00,
02371 & 3.754E+00, 3.759E+00, 3.760E+00, 3.763E+00, 3.766E+00,
02372 & 3.767E+00, 3.770E+00, 3.771E+00, 3.772E+00, 3.776E+00,
02373 & 3.777E+00, 3.780E+00, 3.782E+00, 3.788E+00, 3.790E+00,
02374 & 3.794E+00, 3.799E+00, 3.800E+00, 3.805E+00, 3.810E+00,
02375 & 3.810E+00, 3.816E+00, 3.820E+00, 3.822E+00, 3.827E+00,
02376 & 3.830E+00, 3.833E+00, 3.838E+00, 3.840E+00, 3.844E+00,
02377 & 3.850E+00, 3.850E+00, 3.855E+00, 3.860E+00, 3.861E+00,
02378 & 3.866E+00, 3.870E+00, 3.872E+00, 3.878E+00, 3.880E+00,
02379 & 3.883E+00, 3.889E+00, 3.890E+00, 3.894E+00, 3.900E+00,
02380 & 4.000E+00, 4.005E+00, 4.010E+00, 4.015E+00, 4.020E+00,
02381 & 4.025E+00, 4.030E+00, 4.035E+00, 4.040E+00, 4.045E+00,
02382 & 4.050E+00, 4.055E+00, 4.060E+00, 4.065E+00, 4.070E+00/
02383
02384 DATA (EMA(I), I=201,250) /
02385 & 4.075E+00, 4.080E+00, 4.085E+00, 4.090E+00, 4.095E+00,
02386 & 4.100E+00, 4.105E+00, 4.110E+00, 4.115E+00, 4.120E+00,
02387 & 4.125E+00, 4.130E+00, 4.135E+00, 4.140E+00, 4.145E+00,
02388 & 4.150E+00, 4.155E+00, 4.160E+00, 4.165E+00, 4.170E+00,
02389 & 4.175E+00, 4.180E+00, 4.185E+00, 4.190E+00, 4.195E+00,
02390 & 4.200E+00, 5.500E+00, 6.500E+00, 7.500E+00, 8.500E+00,
02391 & 1.050E+01, 1.100E+01, 1.150E+01, 1.350E+01, 1.362E+01,
02392 & 1.375E+01, 1.388E+01, 1.400E+01, 1.450E+01, 1.500E+01,
02393 & 1.550E+01, 1.600E+01, 1.650E+01, 1.700E+01, 1.750E+01,
02394 & 1.800E+01, 1.850E+01, 1.900E+01, 1.950E+01, 2.000E+01/
02395
02396 DATA ((DAM(I,J), I=1,100), J=1,1) /
02397 & 4.012E-03, 4.261E-03, 4.427E-03, 4.415E-03, 4.318E-03,
02398 & 4.324E-03, 4.407E-03, 4.464E-03, 4.464E-03, 4.368E-03,
02399 & 4.156E-03, 4.090E-03, 4.013E-03, 3.930E-03, 3.818E-03,
02400 & 3.700E-03, 3.561E-03, 3.397E-03, 3.224E-03, 3.014E-03,
02401 & 2.788E-03, 2.519E-03, 2.202E-03, 1.840E-03, 1.398E-03,
02402 & 8.780E-04, 2.395E-04,-5.553E-04,-1.572E-03,-2.925E-03,
02403 & -4.823E-03,-7.339E-03,-1.176E-02,-2.123E-02,-5.422E-02,
02404 & 3.444E-01, 1.010E-01, 6.094E-02, 4.450E-02, 3.562E-02,
02405 & 3.007E-02, 2.630E-02, 2.359E-02, 2.154E-02, 1.994E-02,
02406 & 1.867E-02, 1.763E-02, 1.678E-02, 1.605E-02, 1.535E-02,
02407 & 1.470E-02, 1.412E-02, 1.360E-02, 1.314E-02, 1.273E-02,
02408 & 1.235E-02, 1.200E-02, 1.169E-02, 1.140E-02, 1.113E-02,
02409 & 1.087E-02, 1.064E-02, 1.042E-02, 1.021E-02, 1.001E-02,
02410 & 9.825E-03, 9.647E-03, 9.473E-03, 9.307E-03, 9.147E-03,
02411 & 8.992E-03, 8.841E-03, 8.696E-03, 8.550E-03, 8.415E-03,
02412 & 8.269E-03, 8.090E-03, 7.911E-03, 7.730E-03, 7.538E-03,
02413 & 7.349E-03, 7.192E-03, 7.020E-03, 6.863E-03, 6.704E-03,
02414 & 6.543E-03, 6.503E-03, 6.464E-03, 6.426E-03, 6.398E-03,
02415 & 6.389E-03, 6.459E-03, 6.549E-03, 6.789E-03, 7.459E-03,
02416 & -2.514E-03,-2.685E-03,-3.083E-03,-3.731E-03,-4.668E-03/
02417
02418 DATA ((DAM(I,J), I=101,200), J=1,1) /
02419 & -5.995E-03,-7.868E-03,-1.062E-02,-1.493E-02,-2.246E-02,
02420 & -3.861E-02,-9.647E-02, 3.645E-01, 1.180E-01, 7.257E-02,
02421 & 5.660E-02, 4.302E-02, 4.072E-02, 3.891E-02, 3.246E-02,
02422 & 3.174E-02, 2.737E-02, 2.572E-02, 2.540E-02, 2.392E-02,
02423 & 2.196E-02, 2.140E-02, 1.968E-02, 1.947E-02, 1.935E-02,
02424 & 1.793E-02, 1.740E-02, 1.665E-02, 1.638E-02, 1.588E-02,
02425 & 1.557E-02, 1.462E-02, 1.410E-02, 1.376E-02, 1.349E-02,
02426 & 1.297E-02, 1.249E-02, 1.229E-02, 1.224E-02, 1.160E-02,
02427 & 1.159E-02, 1.110E-02, 1.101E-02, 1.101E-02, 1.133E-02,
02428 & 1.169E-02, 1.252E-02, 1.291E-02, 1.299E-02, 1.395E-02,
02429 & 1.407E-02, 1.417E-02, 1.406E-02, 1.359E-02, 1.341E-02,
02430 & 1.308E-02, 1.261E-02, 1.255E-02, 1.220E-02, 1.187E-02,
02431 & 1.185E-02, 1.154E-02, 1.135E-02, 1.128E-02, 1.104E-02,
02432 & 1.093E-02, 1.083E-02, 1.063E-02, 1.057E-02, 1.045E-02,
02433 & 1.028E-02, 1.028E-02, 1.013E-02, 1.001E-02, 9.988E-03,
02434 & 9.862E-03, 9.775E-03, 9.735E-03, 9.618E-03, 9.569E-03,
02435 & 9.503E-03, 9.402E-03, 9.379E-03, 9.295E-03, 9.199E-03,
02436 & 8.654E-03, 8.639E-03, 8.621E-03, 8.609E-03, 8.601E-03,
02437 & 8.603E-03, 8.609E-03, 8.623E-03, 8.630E-03, 8.648E-03,
02438 & 8.663E-03, 8.680E-03, 8.694E-03, 8.704E-03, 8.719E-03/
02439
02440 DATA ((DAM(I,J), I=201,250), J=1,1) /
02441 & 8.737E-03, 8.747E-03, 8.758E-03, 8.767E-03, 8.772E-03,
02442 & 8.775E-03, 8.782E-03, 8.797E-03, 8.816E-03, 8.851E-03,
02443 & 8.889E-03, 8.928E-03, 8.956E-03, 8.980E-03, 8.991E-03,
02444 & 8.989E-03, 8.979E-03, 8.973E-03, 8.964E-03, 8.954E-03,
02445 & 8.944E-03, 8.944E-03, 8.942E-03, 8.943E-03, 8.941E-03,
02446 & 8.936E-03, 1.150E-02, 1.180E-02, 1.232E-02, 1.315E-02,
02447 & 1.413E-02, 1.550E-02, 1.572E-02, 1.622E-02, 1.628E-02,
02448 & 1.633E-02, 1.638E-02, 1.644E-02, 1.665E-02, 1.686E-02,
02449 & 1.706E-02, 1.725E-02, 1.744E-02, 1.763E-02, 1.780E-02,
02450 & 1.798E-02, 1.815E-02, 1.831E-02, 1.847E-02, 1.862E-02/
02451
02452 DATA ((DAM(I,J), I=1,100), J=2,2) /
02453 & 4.220E-04, 4.386E-04, 4.473E-04, 4.434E-04, 4.609E-04,
02454 & 4.794E-04, 4.711E-04, 4.734E-04, 4.777E-04, 4.855E-04,
02455 & 4.966E-04, 4.965E-04, 4.980E-04, 5.005E-04, 5.039E-04,
02456 & 5.071E-04, 5.091E-04, 5.129E-04, 5.183E-04, 5.254E-04,
02457 & 5.341E-04, 5.406E-04, 5.500E-04, 5.627E-04, 5.792E-04,
02458 & 6.007E-04, 6.263E-04, 6.617E-04, 7.104E-04, 7.803E-04,
02459 & 8.869E-04, 1.358E-03, 2.081E-03, 3.264E-03, 6.100E-03,
02460 & 2.402E-02, 9.518E-03, 6.259E-03, 4.725E-03, 3.764E-03,
02461 & 3.077E-03, 2.549E-03, 2.123E-03, 1.773E-03, 1.483E-03,
02462 & 1.244E-03, 1.055E-03, 9.160E-04, 8.440E-04, 8.120E-04,
02463 & 7.730E-04, 7.390E-04, 7.120E-04, 6.870E-04, 6.660E-04,
02464 & 6.480E-04, 6.330E-04, 6.180E-04, 6.070E-04, 5.980E-04,
02465 & 5.910E-04, 5.850E-04, 5.793E-04, 5.741E-04, 5.705E-04,
02466 & 5.677E-04, 5.663E-04, 5.660E-04, 5.669E-04, 5.678E-04,
02467 & 5.693E-04, 5.715E-04, 5.746E-04, 5.783E-04, 5.831E-04,
02468 & 5.886E-04, 5.918E-04, 5.970E-04, 6.039E-04, 6.127E-04,
02469 & 6.234E-04, 6.345E-04, 6.473E-04, 6.621E-04, 6.793E-04,
02470 & 6.986E-04, 8.741E-04, 1.109E-03, 1.400E-03, 1.752E-03,
02471 & 2.180E-03, 2.711E-03, 3.401E-03, 4.384E-03, 6.076E-03,
02472 & 5.225E-03, 4.546E-03, 4.036E-03, 3.658E-03, 3.403E-03/
02473
02474 DATA ((DAM(I,J), I=101,200), J=2,2) /
02475 & 3.281E-03, 3.323E-03, 3.594E-03, 4.217E-03, 5.497E-03,
02476 & 8.425E-03, 1.913E-02, 6.645E-02, 2.070E-02, 1.226E-02,
02477 & 9.289E-03, 6.767E-03, 6.340E-03, 6.007E-03, 4.826E-03,
02478 & 4.694E-03, 3.904E-03, 3.608E-03, 3.549E-03, 3.287E-03,
02479 & 2.941E-03, 2.844E-03, 2.549E-03, 2.514E-03, 2.494E-03,
02480 & 2.259E-03, 2.174E-03, 2.056E-03, 2.014E-03, 1.938E-03,
02481 & 1.893E-03, 1.760E-03, 1.694E-03, 1.651E-03, 1.619E-03,
02482 & 1.562E-03, 1.512E-03, 1.493E-03, 1.489E-03, 1.431E-03,
02483 & 1.429E-03, 1.374E-03, 1.356E-03, 1.303E-03, 1.246E-03,
02484 & 1.209E-03, 1.163E-03, 1.153E-03, 1.153E-03, 1.159E-03,
02485 & 1.160E-03, 1.153E-03, 1.136E-03, 1.086E-03, 1.067E-03,
02486 & 1.034E-03, 9.880E-04, 9.830E-04, 9.510E-04, 9.200E-04,
02487 & 9.170E-04, 8.900E-04, 8.720E-04, 8.660E-04, 8.450E-04,
02488 & 8.350E-04, 8.260E-04, 8.149E-04, 8.126E-04, 8.076E-04,
02489 & 8.022E-04, 8.018E-04, 7.979E-04, 7.968E-04, 7.964E-04,
02490 & 7.968E-04, 7.976E-04, 7.984E-04, 8.019E-04, 8.039E-04,
02491 & 8.069E-04, 8.138E-04, 8.154E-04, 8.219E-04, 8.320E-04,
02492 & 1.203E-03, 1.243E-03, 1.275E-03, 1.290E-03, 1.275E-03,
02493 & 1.212E-03, 1.094E-03, 9.506E-04, 8.793E-04, 9.563E-04,
02494 & 1.113E-03, 1.252E-03, 1.338E-03, 1.378E-03, 1.388E-03/
02495
02496 DATA ((DAM(I,J), I=201,250), J=2,2) /
02497 & 1.383E-03, 1.372E-03, 1.361E-03, 1.352E-03, 1.346E-03,
02498 & 1.344E-03, 1.345E-03, 1.346E-03, 1.346E-03, 1.341E-03,
02499 & 1.326E-03, 1.297E-03, 1.247E-03, 1.172E-03, 1.074E-03,
02500 & 9.652E-04, 8.729E-04, 8.347E-04, 8.685E-04, 9.535E-04,
02501 & 1.053E-03, 1.140E-03, 1.204E-03, 1.244E-03, 1.264E-03,
02502 & 1.267E-03, 1.066E-03, 1.222E-03, 1.546E-03, 1.694E-03,
02503 & 1.555E-03, 1.751E-03, 1.066E-03, 4.650E-04, 4.610E-04,
02504 & 4.570E-04, 4.530E-04, 4.480E-04, 4.360E-04, 4.250E-04,
02505 & 4.170E-04, 4.090E-04, 4.030E-04, 3.970E-04, 3.930E-04,
02506 & 3.890E-04, 3.850E-04, 3.810E-04, 3.780E-04, 3.760E-04/
02507 C
02508 save
02509 RETURN
02510 END
