4 SUBROUTINE fhad(IPRMOD,IPRO,PLAB,ELAB,CX,CY,CZ,
5 * ithkk,itta,ieline,irejfh)
6 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
29 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
117 COMMON /hadthr/ ehadth,inthad
120 COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin), czrh
121 +(maxfin),elrh(maxfin),plrh(maxfin),irh
133 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
134 +iibar(210),k1(210),k2(210)
137 COMMON /projk/ iprojk
139 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
147 ELSEIF(inthad.EQ.2)
THEN
159 pptar=
sqrt(phkk(1,ithkk)**2 + phkk(2,ithkk)**2 + phkk(3,ithkk)**2)
160 pptarm=
sqrt(abs(phkk(4,ithkk)-amtar)*(phkk(4,ithkk)+amtar))
162 bgx=phkk(1,ithkk)/amtar
163 bgy=phkk(2,ithkk)/amtar
164 bgz=phkk(3,ithkk)/amtar
165 IF(pptar.GT.
zero)
THEN
172 gam=phkk(4,ithkk)/amtar
177 etes = elab + phkk(4,ithkk)
178 pxtes= pxpro + phkk(1,ithkk)
179 pytes= pypro + phkk(2,ithkk)
180 pztes= pzpro + phkk(3,ithkk)
182 CALL
daltra(gam,-bgx,-bgy,-bgz, pxpro,pypro,pzpro,elab, pprof,
183 +pxprof,pyprof,pzprof,eprof)
185 WRITE(6,
'(2A/A,4E12.5/A,5E12.5/A,4E12.5/A,4E12.5/A,4I4)')
186 &
' FHAD: projectile after LT into target',
187 &
' nucleon rest system',
188 &
' GAM,BGX,BGY,BGZ ', gam,bgx,bgy,bgz,
189 &
' PHKK(target):', (phkk(ik,ithkk),ik=1,5),
190 &
' PX/Y/Z/PRO, ELAB:', pxpro,pypro,pzpro,elab,
191 &
' Proj. after LT:', pxprof,pyprof,pzprof,eprof,
192 &
' IPRO, ITTA, ITHKK, IELINE:', ipro,itta,ithkk,ieline
194 IF(eprof.LE.aam(ipro))
THEN
195 WRITE(6,
'(2A/A,5E12.5/A,4E12.5/A,4E12.5/A,4I4)')
196 &
' FHAD: inconsistent projectile after LT into target',
197 &
' nucleon rest system',
198 &
' PHKK(target):', (phkk(ik,ithkk),ik=1,5),
199 &
' PX/Y/Z/PRO, ELAB:', pxpro,pypro,pzpro,elab,
200 &
' Proj. after LT:', pxprof,pyprof,pzprof,eprof,
201 &
' IPRO, ITTA, ITHKK, IELINE:', ipro,itta,ithkk,ieline
208 CALL
daltra(gam,-bgx,-bgy,-bgz, phkk(1,ithkk),phkk(2,ithkk),phkk
209 + (3,ithkk),phkk(4,ithkk), pptar,pxtarf,pytarf,pztarf,etarf)
211 WRITE(6,
'(A)')
'FHAD: TARGET MOM. BEFORE/AFTER LORENTZ TRANSF.'
212 WRITE(6,
'(3X,A,5(1PE12.4))')
' PHKK(1-4,ITHKK)', (phkk
214 WRITE(6,
'(3X,A,5(1PE12.4))')
' PXTARF,PYTARF,PZTARF,ETARF',
215 + pxtarf,pytarf,pztarf,etarf
216 WRITE(6,
'(A)')
'FHAD: PROJ. MOM. BEFORE/AFTER LORENTZ TRANSF.'
217 WRITE(6,
'(3X,A,5(1PE12.4))')
' PXPRO,PYPRO,PZPRO,ELAB', pxpro,
219 WRITE(6,
'(3X,A,5(1PE12.4))')
' PXPROF,PYPROF,PZPROF,EPROF',
220 + pxprof,pyprof,pzprof,eprof
232 CALL
dhadri(ipro,pprof,eprof,cxf,cyf,czf,itta)
234 IF (ipaupr.GT.2)
WRITE(6,1000)ipro,pprof,eprof,cxf,cyf,czf,itta
235 1000
FORMAT (
' FHAD IPRO,PPROF,EPROF,CXF,CYF,CZF,ITTA',i5,5f10.2,i5)
236 ELSEIF(iielin.EQ.1)
THEN
237 CALL
elhain(ipro,pprof,eprof,cxf,cyf,czf,itta,irej)
249 cxrh(2)=phkk(1,ithkk)/pptar
250 cyrh(2)=phkk(2,ithkk)/pptar
251 czrh(2)=phkk(3,ithkk)/pptar
252 elrh(2)=phkk(4,ithkk)
261 crsum=cxrh(iii)**2 + cyrh(iii)**2 + czrh(iii)**2
262 IF(abs(crsum-1.0).GT.1
e-4)
THEN
263 WRITE(6,
'(A,I3,1PE12.4)')
264 +
' FHAD: INCORRECT NORM. OF DIRECTION COSINES - III,CRSUM',
267 cxrh(iii)=cxrh(iii)/rcrsum
268 cyrh(iii)=cyrh(iii)/rcrsum
269 czrh(iii)=czrh(iii)/rcrsum
278 WRITE(6,
'(A,I3,6(1PE12.4))')
279 +
' FHAD: ITRH(I), AMI,ELR,PPS,PSX,PSY,PSZ', itrh(iii),ami,elrh
280 + (iii),pps,psx,psy,psz
282 CALL
daltra(gam,bgx,bgy,bgz, psx,psy,psz,elrh(iii), ppps,ppsx,
283 + ppsy,ppsz,elrh(iii))
284 IF(elrh(iii).LT.(ppps-1
d-4))
THEN
285 WRITE(6,
'(2A/3I3,6(1PE12.4))')
286 +
' FHAD: INCONSISTENT KINEMATICS AFTER ALTRA: ',
287 +
' IIELIN,III,ITRH(III),ELRH(III),PPPS,PSX,PSY,PSZ,AMI',iielin,
288 + iii,itrh(iii),elrh(iii),ppps,psx,psy,psz,ami
289 elrh(iii)=
sqrt(ppps**2 + ami**2)
290 WRITE(6,
'(A,1PE12.4)')
' CORRECTED ENERGY ELRH:',elrh(iii)
291 WRITE(6,
'(2A/2I5,5(1PE12.4))')
292 +
' FHAD: 4-MOM. OF TARGET NUCLEON -',
293 +
' ITHKK, IDHKK , PHKK(1-4)', ithkk, idhkk(ithkk), (phkk
300 etes = etes - elrh(iii)
306 IF(abs(etes).GT.0.001d0)
THEN
309 WRITE(6,
'(A,I5)')
' FHAD: TEST OF E-P CONSERVATION IELINE=',
311 WRITE(6,
'(3X,A,5(1PE12.4))')
' ETES,PXTES,PYTES,PZTES', etes,
313 WRITE(6,1000)ipro,pprof,eprof,cxf,cyf,czf,itta
317 pps=
sqrt((elrh(iii)-ami)*(elrh(iii)+ami))
322 WRITE(6,
'(A,I3,6(1PE12.4))')
323 +
' FHAD: ITRH(I), AMI,ELRH,PPS,PSX,PSY,PSZ', itrh(iii),ami,
324 + elrh(iii),pps,psx,psy,psz
335 SUBROUTINE elhain(IP,PLA,ELAB,CX,CY,CZ,IT,IREJ)
343 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
345 parameter(lout=6,llook=9)
346 parameter(two=2.0d0,
one=1.0d0,ohalf=0.5d0,
zero=0.0d0,
349 parameter(ennthr = 3.5d0)
350 parameter(plowh=0.01d0,phih=9.0d0,
351 & blowb=0.05d0,bhib=0.2d0,
352 & blowm=0.1d0, bhim=2.0d0)
355 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),
356 & iich(210),iibar(210),k1(210),k2(210)
359 COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin),
360 & czrh(maxfin),elrh(maxfin),plrh(maxfin),irh
366 plab =
sqrt( (elab-aam(ip))*(elab+aam(ip)) )
371 s = amp2+amt2+two*elab*aam(it)
373 ecmp = ohalf*ecm+(amp2-amt2)/(two*ecm)
374 pcm =
sqrt( (ecmp-aam(ip))*(ecmp+aam(ip)) )
377 IF ( ((ip.EQ.1).OR.(ip.EQ.8)).AND.
378 & ((it.EQ.1).OR.(it.EQ.8)).AND.(
ekin.LT.ennthr) )
THEN
385 IF (ip.EQ.8) kproj = 1
388 t = two*pcm**2*(ctcms-
one)
393 tmax = two*two*pcm**2
394 rr = (plab-plowh)/(phih-plowh)
395 IF (iibar(ip).NE.0)
THEN
396 tslope = blowb+rr*(bhib-blowb)
398 tslope = blowm+rr*(bhim-blowm)
400 fmax =
exp(-tslope*tmax)-
one
407 elrh(2) = (two*amt2-
t)/(two*aam(it))
408 plrh(2) =
sqrt( (elrh(2)-aam(it))*(elrh(2)+aam(it)) )
410 elrh(1) = elab+aam(it)-elrh(2)
411 plrh(1) =
sqrt( (elrh(1)-aam(ip))*(elrh(1)+aam(ip)) )
413 ctlabp = (
t-two*amp2+two*elab*elrh(1))/(two*plab*plrh(1))
415 CALL
dsfecf(splabp,cplabp)
418 CALL
drtran(cx,cy,cz,ctlabp,stlabp,splabp,cplabp,
419 & cxrh(1),cyrh(1),czrh(1))
421 pllabt = plab-ctlabp*plrh(1)
422 ctlabt = pllabt/plrh(2)
426 CALL
drtran(cx,cy,cz,ctlabt,stlabt,-splabp,-cplabp,
427 & cxrh(2),cyrh(2),czrh(2))
449 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
477 dimension sikpp(8),sikmp(8),siapp(8),
p(8)
478 DATA p /0.3d0,0.4d0,0.5d0,0.6d0,0.8d0,1.d0,1.5d0,2.d0/
479 DATA sikpp / 12.0d0, 12.5d0, 13.0d0, 13.0d0, 12.7d0,
480 + 12.0d0, 10.2d0, 6.82d0/
481 DATA sikmp / 42.0d0, 33.0d0, 21.0d0, 16.0d0, 19.0d0,
482 + 22.0d0, 9.0d0, 7.5d0 /
483 DATA siapp / 73.0d0, 70.0d0, 62.0d0, 53.0d0, 48.0d0,
484 + 43.0d0, 38.0d0, 33.0d0/
495 ELSEIF(iproj.EQ.14)
THEN
497 ELSEIF(iproj.EQ.1)
THEN
499 ELSEIF(iproj.EQ.8)
THEN
501 ELSEIF(iproj.EQ.15)
THEN
503 ELSEIF(iproj.EQ.16)
THEN
505 ELSEIF(iproj.EQ.24)
THEN
507 ELSEIF(iproj.EQ.25)
THEN
514 ELSEIF(ippr.EQ.17)
THEN
516 ELSEIF(ippr.EQ.18)
THEN
518 ELSEIF(ippr.EQ.24)
THEN
520 ELSEIF(ippr.EQ.25)
THEN
522 ELSEIF(ippr.GE.20.AND.ippr.LE.22)
THEN
528 IF(ippr.EQ.15.OR.ippr.EQ.16.OR.ippr.EQ.2)
THEN
535 IF(ppoo.LE.
p(jk))
THEN
549 ELSEIF(ippr.EQ.16)
THEN
553 ELSEIF(ippr.EQ.2)
THEN
559 +
' LOGICAL ERROR IN SIHNEL - EXECUTION STOPPED'
562 siel=s1 + (s2-s1)*(ppoo-
p(kk))/(
p(k)-
p(kk)+1
d-7)
565 ELSEIF(ppoo.LE.10.0d0)
THEN
573 ELSEIF(ippr.EQ.16)
THEN
580 ELSEIF(ippr.EQ.2)
THEN
590 siel=a1 + a2*ppoo**an + a3*alp**2 + a4*alp
599 CALL
dsige(ippr,aa,ppoo,siel,zlel)
607 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
669 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
670 +iibar(210),k1(210),k2(210)
675 dimension sig(13,9),seg(16,9),
p(16),
a(9),itt(39)
676 dimension plab(19),sito(19,4),alp(3),bet(3)
677 dimension rea(9,9),stot(9)
678 SAVE a,
p, sig, seg, itt, sito, plab, alp, bet, stot, rea
679 DATA a/9.d0,12.d0,27.d0,47.9d0,55.9d0,63.5d0,112.4d0,
681 DATA p/.13d0,.19d0,.25d0,.3d0,.4d0,.5d0,.6d0,.8d0,1.d0,
682 &1.5d0,2.d0,3.d0,4.d0,5.d0,6.d0,10.d0/
683 DATA sig/ 485.d0,223.d0,112.d0,82.d0,66.d0,78.d0,96.d0,102.d0,
684 &100.d0,98.d0,95.d0,90.d0,79.d0,
685 (680.d0,348.d0,175.d0,103.d0,84.d0,87.d0,106.d0,112.d0,111.d0,
686 &108.d0,107.d0,105.d0,101.d0,
687 (1200.d0,738.d0,387.d0,196.d0,191.d0,200.d0,248.d0,264.d0,
688 &264.d0,257.d0,252.d0,247.d0,228.d0,
689 (1658.d0,1110.d0,635.d0,364.d0,332.d0,356.d0,404.d0,408.d0,
690 &407.d0,404.d0,398.d0,396.d0,384.d0,
691 (1730.d0,1270.d0,725.d0,400.d0,375.d0,412.d0,495.d0,505.d0,
692 &495.d0,492.d0,487.d0,485.d0,475.d0,
693 (1875.d0,1470.d0,835.d0,480.d0,450.d0,450.d0,535.d0,
694 &580.d0,555.d0,540.d0,535.d0,530.d0,
695 (525.d0,2040.d0,2160.d0,1335.d0,850.d0,740.d0,760.d0,880.d0,
696 &905.d0,860.d0,840.d0,820.d0,815.d0,800.d0,
697 (2340.d0,2980.d0,2270.d0,1450.d0,1230.d0,1230.d0,
698 &1380.d0,1420.d0,1410.d0,1380.d0,1360.d0,
699 (1350.d0,1320.d0,2680.d0,3220.d0,2530.d0,1630.d0,1420.d0,1450.d0,
700 &1570.d0,1600.d0,1590.d0,1575.d0,1560.d0,1550.d0,1540.d0/
701 DATA seg/24.d0,128.d0,249.d0,256.d0,202.d0,124.d0,73.d0,60.d0,
702 &64.d0,69.d0,62.d0,50.d0,44.d0,42.d0,42.d0,41.d0,21.d0,156.d0,
703 (273.d0,280.d0,220.d0,212.d0,94.d0,80.d0,82.d0,85.d0,80.d0,73.d0,
704 (69.d0,67.d0,66.d0,64.d0,56.d0,296.d0,560.d0,574.d0,467.d0,350.d0,
705 &235.d0,210.d0,210.d0,200.d0,190.d0,183.d0,176.d0,170.d0,165.d0,
706 (155.d0,100.d0,500.d0,895.d0,880.d0,690.d0,520.d0,378.d0,
707 (355.d0,384.d0,373.d0,352.d0,320.d0,300.d0,288.d0,280.d0,
708 &262.d0,75.d0,500.d0,965.d0,990.d0,775.d0,525.d0,410.d0,410.d0,
709 (433.d0,440.d0,425.d0,395.d0,374.d0,355.d0,340.d0,303.d0,125.d0,
710 (570.d0,1025.d0,1100.d0,825.d0,575.d0,418.d0,458.d0,500.d0,
711 &480.d0,460.d0,440.d0,422.d0,400.d0,384.d0,355.d0,300.d0,880.d0,
712 (1480.d0,1550.d0,1380.d0,940.d0,710.d0,720.d0,810.d0,760.d0,
713 (740.d0,700.d0,665.d0,645.d0,620.d0,570.d0,550.d0,1475.d0,
714 &2250.d0,2350.d0,1850.d0,1500.d0,
715 (1120.d0,1210.d0,1480.d0,1440.d0,1400.d0,1320.d0,
716 &1250.d0,1210.d0,1170.d0,1065.d0,540.d0,
717 (1300.d0,2220.d0,2560.d0,1980.d0,1650.d0,1160.d0,
718 &1360.d0,1600.d0,1560.d0,1510.d0,1410.d0,
719 (1350.d0,1300.d0,1270.d0,1200.d0/
721 DATA itt/1,7,0,0,0,0,0,2,8,0,0,9,3,4,6,5,2,8,9,1,1,2,3,9,10,
722 & 3,0,0,0,0,7,2,7,2,8,1,7,1,7/
723 DATA plab/.3d0,.4d0,.5d0,.6d0,.7d0,.8d0,.9d0,1.d0,1.1d0,
724 &1.2d0,1.3d0,1.4d0,1.5d0,2.d0,3.d0,4.d0,
726 DATA sito/66.8d0,63.6d0,40.35d0,31.25d0,31.1d0,
727 *35.1d0,36.7d0,44.15d0,38.3d0,33.25d0,
728 *29.75d0,29.3d0,29.95d0,26.55d0,24.6d0,22.95d0,
729 *22.75d0,22.95d0,21.55d0,
730 *12.5d0,14.1d0,13.5d0,12.75d0,12.85d0,13.9d0,15.6d0,
731 *17.25d0,18.9d0,19.5d0,18.95d0,18.85d0,
732 *18.45d0,18.2d0,17.5d0,17.7d0,17.5d0,17.25d0,17.4d0,
733 *39.65d0,38.75d0,26.9d0,22.d0,22.d0,24.5d0,26.15d0,30.7d0,28.6d0,
734 &26.4d0,24.35d0,24.1d0,24.2d0
735 (,22.4d0,21.05d0,20.3d0,20.1d0,20.1d0,19.5d0,
736 (280.d0,199.7d0,171.1d0,154.3d0,140.d0,130.d0,116.8d0,117.4d0,
737 &111.6d0,109.d0,106.5d0,
738 (102.8d0,100.d0,90.2d0,76.7d0,68.d0,62.8d0,60.7d0,56.d0/
739 DATA alp/0.823d0,0.843d0,0.630d0/
740 DATA bet/1.26d0,1.31d0,0.90d0/
741 DATA stot /15.d0,20.d0,30.d0,40.d0,60.d0,80.d0,
742 &100.d0,150.d0,200.d0/
743 DATA rea / .20d0,.23d0,.27d0,.30d0,.35d0,.40d0,.47d0,.55d0,.60d0,
744 2 .22d0,.26d0,.31d0,.35d0,.40d0,.45d0,.51d0,.59d0,.63d0,
745 3 .24d0,.29d0,.36d0,.42d0,.50d0,.56d0,.60d0,.66d0,.68d0,
746 4 .26d0,.32d0,.42d0,.49d0,.58d0,.63d0,.66d0,.71d0,.72d0,
747 5 .27d0,.33d0,.44d0,.51d0,.61d0,.65d0,.68d0,.72d0,.74d0,
748 6 .28d0,.35d0,.46d0,.53d0,.63d0,.66d0,.69d0,.73d0,.745d0,
749 7 .35d0,.42d0,.53d0,.62d0,.69d0,.72d0,.74d0,.77d0,.78d0,
750 8 .42d0,.51d0,.62d0,.69d0,.75d0,.77d0,.79d0,.81d0,.82d0,
751 9 .44d0,.53d0,.64d0,.70d0,.76d0,.78d0,.80d0,.81d0,.82d0 /
770 IF((aa.LT.1.5).AND. (ipr.EQ.1.OR.ipr.EQ.8.OR.ipr.EQ.13.OR.ipr.EQ.
772 eke=
sqrt(po**2+aam(ipr)**2) - aam(ipr)
773 CALL
dsihae(ipr,eke,po,aa,sel)
784 IF(it.EQ.8.AND.po.LT.20.0d0)
THEN
789 eke=
sqrt(po**2+aam(it)**2) - aam(it)
790 CALL
dsihae(it,eke,po,aa,sel)
791 IF(ipol.EQ.1) goto 240
802 IF((iit.EQ.7).OR.(iit.EQ.8)) goto 250
803 IF (po.GT.20.d0) goto 250
804 IF(po.LE.10.d0) goto 10
808 IF(iit.LE.4) go to 40
815 IF(po.LE.plab(k)) go to 30
827 IF(po.LE.
p(k)) go to 60
832 IF(aa.LE.
a(i)) go to 70
846 go to(100,100,110,110,120,130,140,150,130,120),iit
850 aloga=
log(
a(j+1)/
a(j))
852 si1=sig(k,j)* aaa **(
log(sig(k,j+1)/sig(k,j))/aloga)
855 si2=sig(kk,j)* aaa **(
log(sig(kk,j+1)/sig(kk,j))/aloga)
858 si=si1+(po-
p(k))*(si2-si1)/(
p(kk)-
p(k))
865 aloga=
log(
a(j+1)/
a(j))
867 si1=seg(k,j)* aaa **(
log(seg(k,j+1)/seg(k,j))/aloga)
870 si2=seg(kk,j)* aaa **(
log(seg(kk,j+1)/seg(kk,j))/aloga)
871 si=si1+(po-
p(k))*(si2-si1)/(
p(kk)-
p(k))
903 si=(sito(k,is)-sikk)*(po-pkk)/(plab(k)-pkk)+sikk
907 si=bet(ia)*si*aa**alp(ia)
909 IF(iv.NE.24) go to 190
913 190
IF(iv.EQ.25) iv=16
915 200 CALL
sihnin(iv,1,po,sinel)
918 IF(ipol.EQ.1) goto 240
926 IF(ipol.EQ.1) goto 240
929 IF(sel.LT.1.
d-15) sel=1.
d-15
930 zl=10000.d0*aa/(6.022*sel)
951 IF((it.EQ.2).OR.(it.EQ.9)) it1=1
958 IF(aa.GT.
a(ia)) goto 260
968 IF(sto.GT.stot(is)) goto 280
978 r1=rr + da2*(rea(js,ja+1)-rr)/da1
980 r2=rr + da2*(rea(js+1,ja+1)-rr)/da1
981 ract=r1 + (sto-stot(js))*(r2-r1)/(stot(js+1)-stot(js))
984 CALL
sihnin(it,1,po,sinel1)
986 IF(ipol.EQ.1) goto 300
991 sel=sel + (sel1-sel)*(poo-10.)/10.
999 1000
FORMAT(
' WARNING AT CALL SIGEL ',i5)
1005 SUBROUTINE dsihae(KPROJ,EKIN,PLAB,ANUC,SIGELA)
1006 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1032 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1033 +iibar(210),k1(210),k2(210)
1039 dimension ekihn(nen),ekiha(nea),amass(nnaa)
1040 dimension sepimp(nen),sepipp(nen),sepp(nen),senp(nen)
1041 dimension sena(nea,nnaa),sepa(nea,nnaa)
1088 parameter( kalgnm = 2 )
1089 parameter( anglgb = 5.0
d-16 )
1090 parameter( anglsq = 2.5
d-31 )
1091 parameter( axcssv = 0.2
d+16 )
1092 parameter( andrfl = 1.0
d-38 )
1093 parameter( avrflw = 1.0
d+38 )
1094 parameter( ainfnt = 1.0
d+30 )
1095 parameter( azrzrz = 1.0
d-30 )
1096 parameter( einfnt = +69.07755278982137
d+00 )
1097 parameter( ezrzrz = -69.07755278982137
d+00 )
1098 parameter( onemns = 0.999999999999999
d+00 )
1099 parameter( onepls = 1.000000000000001
d+00 )
1100 parameter( csnnrm = 2.0
d-15 )
1101 parameter( dmxtrn = 1.0
d+08 )
1133 parameter( zerzer = 0.
d+00 )
1134 parameter( oneone = 1.
d+00 )
1135 parameter( twotwo = 2.
d+00 )
1136 parameter( thrthr = 3.
d+00 )
1137 parameter( foufou = 4.
d+00 )
1138 parameter( fivfiv = 5.
d+00 )
1139 parameter( sixsix = 6.
d+00 )
1140 parameter( sevsev = 7.
d+00 )
1141 parameter( eigeig = 8.
d+00 )
1142 parameter( aninen = 9.
d+00 )
1143 parameter( tenten = 10.
d+00 )
1144 parameter( hlfhlf = 0.5
d+00 )
1145 parameter( onethi = oneone / thrthr )
1146 parameter( twothi = twotwo / thrthr )
1147 parameter( pipipi = 3.1415926535897932270
d+00 )
1148 parameter( eneper = 2.7182818284590452354
d+00 )
1149 parameter( sqrent = 1.6487212707001281468
d+00 )
1191 parameter( clight = 2.99792458
d+10 )
1192 parameter( avogad = 6.0221367
d+23 )
1193 parameter( amelgr = 9.1093897
d-28 )
1194 parameter( plckbr = 1.05457266
d-27 )
1195 parameter( elccgs = 4.8032068
d-10 )
1196 parameter( elcmks = 1.60217733
d-19 )
1197 parameter( amugrm = 1.6605402
d-24 )
1198 parameter( ammumu = 0.113428913
d+00 )
1213 parameter( alpfsc = 7.2973530791728595
d-03 )
1214 parameter( fscto2 = 5.3251361962113614
d-05 )
1215 parameter( fscto3 = 3.8859399018437826
d-07 )
1216 parameter( fscto4 = 2.8357075508200407
d-09 )
1217 parameter( plabrc = 0.197327053
d+00 )
1218 parameter( amelct = 0.51099906
d-03 )
1219 parameter( amugev = 0.93149432
d+00 )
1220 parameter( ammuon = 0.105658389
d+00 )
1221 parameter( rclsel = 2.8179409183694872
d-13 )
1222 parameter( gevmev = 1.0
d+03 )
1223 parameter( emvgev = 1.0
d-03 )
1224 parameter( raddeg = 180.
d+00 / pipipi )
1225 parameter( degrad = pipipi / 180.
d+00 )
1252 parameter( lunin = 5 )
1253 parameter( lunout = 6 )
1254 parameter( lunerr = 66 )
1255 parameter( lunber = 14 )
1256 parameter( lunech = 8 )
1257 parameter( lunflu = 86 )
1258 parameter( lungeo = 16 )
1259 parameter( lunpgs = 12 )
1260 parameter( lunran = 2 )
1261 parameter( lunxsc = 81 )
1262 parameter( lunrdb = 1 )
1291 parameter( mxxrgn = 500 )
1292 parameter( mxxmdf = 56 )
1293 parameter( mxxmde = 50 )
1294 parameter( mfstck = 1000 )
1295 parameter( mestck = 100 )
1296 parameter( nallwp = 39 )
1297 parameter( mpdpdx = 8 )
1298 parameter( icomax = 180 )
1299 parameter( nstbis = 304 )
1300 parameter( idmaxp = 210 )
1305 SAVE ekihn, ekiha, amass, sepimp, sepipp, sepp, senp, sena, sepa,
1308 sigmlw(
e) = 3.
d+03 * pipipi
1309 & / ( 1.206
d+03 *
e + ( -1.86
d+00
1310 & + 0.09415
d+03 *
e + 0.0001306
d+06 *
e**2
1311 & )**2 ) + 1.
d+03 * pipipi / ( 1.206
d+03 *
e
1312 & + ( 0.4223
d+00 + 0.13
d+03 *
e )**2 )
1313 sigmpn(betapr) = 34.10
d+00 / betapr**2 - 82.2
d+00 / betapr
1315 sigmpp(betapr) = 10.63
d+00 / betapr**2 - 29.92
d+00 / betapr
1321 & 0.00d0, 0.02d0, 0.04d0, 0.06d0, 0.08d0, 0.10d0, 0.12d0, 0.14d0,
1322 & 0.16d0, 0.18d0, 0.20d0, 0.22d0, 0.24d0, 0.26d0, 0.28d0, 0.30d0,
1323 & 0.32d0, 0.34d0, 0.36d0, 0.38d0, 0.40d0, 0.42d0, 0.44d0, 0.46d0,
1324 & 0.48d0, 0.50d0, 0.52d0, 0.54d0, 0.56d0, 0.58d0, 0.60d0, 0.62d0,
1325 & 0.64d0, 0.66d0, 0.68d0, 0.70d0, 0.72d0, 0.74d0, 0.76d0, 0.78d0,
1326 & 0.80d0, 0.82d0, 0.84d0, 0.86d0, 0.88d0, 0.90d0, 0.92d0, 0.94d0,
1327 & 0.96d0, 0.98d0, 1.00d0, 1.02d0, 1.04d0, 1.06d0, 1.08d0, 1.10d0,
1328 & 1.12d0, 1.14d0, 1.16d0, 1.18d0, 1.20d0, 1.22d0, 1.24d0, 1.26d0,
1329 & 1.28d0, 1.30d0, 1.32d0, 1.34d0, 1.36d0, 1.38d0, 1.40d0, 1.42d0,
1330 & 1.44d0, 1.46d0, 1.48d0, 1.50d0, 1.52d0, 1.54d0, 1.56d0, 1.58d0,
1331 & 1.60d0, 1.62d0, 1.64d0, 1.66d0, 1.68d0, 1.70d0, 1.72d0, 1.74d0,
1332 & 1.76d0, 1.78d0, 1.80d0, 1.82d0, 1.84d0, 1.86d0, 1.88d0, 1.90d0,
1333 & 1.92d0, 1.94d0, 1.96d0, 1.98d0, 2.00d0, 2.5d0, 3.0d0, 3.5d0,
1336 & 0.015d0, 0.02d0, 0.025d0, 0.03d0, 0.04d0, 0.05d0, 0.06d0,
1337 & 0.08d0, 0.10d0, 0.125d0, 0.15d0, 0.175d0, 0.20d0, 0.225d0,
1338 & 0.25d0, 0.3d0, 0.4d0, 0.6d0, 1.0d0, 2.0d0, 5.0d0,
1341 & 4.d0, 9.d0, 12.d0, 27.d0, 47.9d0, 55.9d0, 63.5d0, 112.4d0,
1346 DATA (sepimp(ie),ie=1,50) /
1347 * 1.250
d+00, 1.500
d+00, 1.750
d+00, 2.450
d+00, 3.800
d+00,
1348 * 6.000
d+00, 9.700
d+00, 1.500
d+01, 2.140
d+01, 2.310
d+01,
1349 * 2.295
d+01, 2.070
d+01, 1.795
d+01, 1.550
d+01, 1.360
d+01,
1350 * 1.230
d+01, 1.130
d+01, 1.070
d+01, 1.050
d+01, 1.070
d+01,
1351 * 1.120
d+01, 1.175
d+01, 1.235
d+01, 1.300
d+01, 1.400
d+01,
1352 * 1.500
d+01, 1.600
d+01, 1.700
d+01, 1.835
d+01, 1.970
d+01,
1353 * 2.050
d+01, 1.915
d+01, 1.770
d+01, 1.650
d+01, 1.570
d+01,
1354 * 1.520
d+01, 1.510
d+01, 1.525
d+01, 1.550
d+01, 1.600
d+01,
1355 * 1.685
d+01, 1.800
d+01, 2.000
d+01, 2.230
d+01, 2.475
d+01,
1356 * 2.635
d+01, 2.510
d+01, 2.300
d+01, 2.140
d+01, 2.000
d+01/
1357 DATA (sepimp(ie),ie=51,106) /
1358 * 1.870
d+01, 1.750
d+01, 1.670
d+01, 1.585
d+01, 1.505
d+01,
1359 * 1.440
d+01, 1.395
d+01, 1.340
d+01, 1.299
d+01, 1.260
d+01,
1360 * 1.215
d+01, 1.175
d+01, 1.140
d+01, 1.099
d+01, 1.060
d+01,
1361 * 1.040
d+01, 1.010
d+01, 9.990
d+00, 9.900
d+00, 9.750
d+00,
1362 * 9.600
d+00, 9.550
d+00, 9.450
d+00, 9.350
d+00, 9.250
d+00,
1363 * 9.250
d+00, 9.350
d+00, 9.650
d+00, 9.850
d+00, 1.000
d+01,
1364 * 1.015
d+01, 1.030
d+01, 1.060
d+01, 1.080
d+01, 1.095
d+01,
1365 * 1.100
d+01, 1.095
d+01, 1.090
d+01, 1.070
d+01, 1.035
d+01,
1366 * 1.000
d+01, 9.600
d+00, 9.050
d+00, 8.550
d+00, 8.200
d+00,
1367 * 8.000
d+00, 7.850
d+00, 7.800
d+00, 7.750
d+00, 7.700
d+00,
1369 * 7.600
d+00, 7.240
d+00, 6.770
d+00, 5.840
d+00, 4.570
d+00/
1375 DATA (sepipp(ie),ie=1,50) /
1376 * 1.800
d+00, 4.000
d+00, 9.900
d+00, 2.170
d+01, 4.000
d+01,
1377 * 6.580
d+01, 9.680
d+01, 1.392
d+02, 1.800
d+02, 2.000
d+02,
1378 * 1.655
d+02, 1.420
d+02, 1.225
d+02, 1.032
d+02, 8.400
d+01,
1379 * 6.725
d+01, 5.510
d+01, 4.725
d+01, 4.130
d+01, 3.690
d+01,
1380 * 3.230
d+01, 2.885
d+01, 2.600
d+01, 2.300
d+01, 2.090
d+01,
1381 * 1.875
d+01, 1.675
d+01, 1.500
d+01, 1.340
d+01, 1.200
d+01,
1382 * 1.100
d+01, 9.980
d+00, 9.200
d+00, 8.600
d+00, 8.200
d+00,
1383 * 8.100
d+00, 8.100
d+00, 8.250
d+00, 8.500
d+00, 8.750
d+00,
1384 * 9.000
d+00, 9.400
d+00, 9.750
d+00, 1.000
d+01, 1.030
d+01,
1385 * 1.075
d+01, 1.130
d+01, 1.200
d+01, 1.275
d+01, 1.330
d+01/
1386 DATA (sepipp(ie),ie=51,106) /
1387 * 1.350
d+01, 1.335
d+01, 1.330
d+01, 1.330
d+01, 1.345
d+01,
1388 * 1.355
d+01, 1.380
d+01, 1.400
d+01, 1.460
d+01, 1.500
d+01,
1389 * 1.555
d+01, 1.625
d+01, 1.700
d+01, 1.800
d+01, 1.875
d+01,
1390 * 1.920
d+01, 1.925
d+01, 1.890
d+01, 1.830
d+01, 1.790
d+01,
1391 * 1.725
d+01, 1.690
d+01, 1.640
d+01, 1.600
d+01, 1.550
d+01,
1392 * 1.505
d+01, 1.475
d+01, 1.430
d+01, 1.400
d+01, 1.365
d+01,
1393 * 1.335
d+01, 1.300
d+01, 1.280
d+01, 1.250
d+01, 1.225
d+01,
1394 * 1.205
d+01, 1.195
d+01, 1.175
d+01, 1.150
d+01, 1.135
d+01,
1395 * 1.105
d+01, 1.095
d+01, 1.080
d+01, 1.060
d+01, 1.030
d+01,
1396 * 1.020
d+01, 1.005
d+01, 9.900
d+00, 9.800
d+00, 9.700
d+00,
1398 * 7.350
d+00, 7.200
d+00, 7.000
d+00, 5.800
d+00, 4.800
d+00/
1402 DATA (sepp(ie),ie=1,50) /
1403 * 6.750
d+02, 1.550
d+02, 6.750
d+01, 4.420
d+01, 3.230
d+01,
1404 * 2.800
d+01, 2.520
d+01, 2.370
d+01, 2.300
d+01, 2.275
d+01,
1405 * 2.260
d+01, 2.260
d+01, 2.260
d+01, 2.260
d+01, 2.270
d+01,
1406 * 2.280
d+01, 2.295
d+01, 2.300
d+01, 2.310
d+01, 2.330
d+01,
1407 * 2.350
d+01, 2.380
d+01, 2.395
d+01, 2.420
d+01, 2.460
d+01,
1408 * 2.485
d+01, 2.500
d+01, 2.530
d+01, 2.565
d+01, 2.600
d+01,
1409 * 2.620
d+01, 2.640
d+01, 2.660
d+01, 2.675
d+01, 2.690
d+01,
1410 * 2.700
d+01, 2.705
d+01, 2.710
d+01, 2.715
d+01, 2.720
d+01,
1411 * 2.725
d+01, 2.725
d+01, 2.720
d+01, 2.715
d+01, 2.710
d+01,
1412 * 2.700
d+01, 2.695
d+01, 2.680
d+01, 2.670
d+01, 2.660
d+01/
1413 DATA (sepp(ie),ie=51,106) /
1414 * 2.640
d+01, 2.625
d+01, 2.605
d+01, 2.590
d+01, 2.570
d+01,
1415 * 2.545
d+01, 2.525
d+01, 2.500
d+01, 2.480
d+01, 2.470
d+01,
1416 * 2.450
d+01, 2.430
d+01, 2.410
d+01, 2.395
d+01, 2.370
d+01,
1417 * 2.360
d+01, 2.340
d+01, 2.325
d+01, 2.305
d+01, 2.290
d+01,
1418 * 2.275
d+01, 2.270
d+01, 2.260
d+01, 2.250
d+01, 2.230
d+01,
1419 * 2.225
d+01, 2.210
d+01, 2.200
d+01, 2.195
d+01, 2.190
d+01,
1420 * 2.175
d+01, 2.165
d+01, 2.150
d+01, 2.140
d+01, 2.125
d+01,
1421 * 2.120
d+01, 2.105
d+01, 2.100
d+01, 2.090
d+01, 2.075
d+01,
1422 * 2.065
d+01, 2.055
d+01, 2.045
d+01, 2.030
d+01, 2.020
d+01,
1423 * 2.005
d+01, 2.000
d+01, 1.995
d+01, 1.980
d+01, 1.975
d+01,
1425 * 17.15
d+00, 14.45
d+00, 13.00
d+00, 11.50
d+00, 10.50
d+00/
1429 DATA (senp(ie),ie=1,50) /
1430 * 1.965
d+03, 4.750
d+02, 2.200
d+02, 1.300
d+02, 9.180
d+01,
1431 * 7.300
d+01, 6.030
d+01, 5.180
d+01, 4.680
d+01, 4.320
d+01,
1432 * 4.080
d+01, 3.910
d+01, 3.760
d+01, 3.650
d+01, 3.550
d+01,
1433 * 3.480
d+01, 3.415
d+01, 3.370
d+01, 3.325
d+01, 3.290
d+01,
1434 * 3.275
d+01, 3.250
d+01, 3.255
d+01, 3.275
d+01, 3.285
d+01,
1435 * 3.275
d+01, 3.220
d+01, 3.150
d+01, 3.075
d+01, 2.990
d+01,
1436 * 2.875
d+01, 2.775
d+01, 2.695
d+01, 2.630
d+01, 2.590
d+01,
1437 * 2.565
d+01, 2.560
d+01, 2.560
d+01, 2.560
d+01, 2.565
d+01,
1438 * 2.570
d+01, 2.575
d+01, 2.578
d+01, 2.580
d+01, 2.585
d+01,
1439 * 2.580
d+01, 2.575
d+01, 2.560
d+01, 2.540
d+01, 2.505
d+01/
1440 DATA (senp(ie),ie=51,106) /
1441 * 2.470
d+01, 2.425
d+01, 2.375
d+01, 2.315
d+01, 2.275
d+01,
1442 * 2.230
d+01, 2.200
d+01, 2.175
d+01, 2.155
d+01, 2.145
d+01,
1443 * 2.130
d+01, 2.125
d+01, 2.115
d+01, 2.105
d+01, 2.100
d+01,
1444 * 2.095
d+01, 2.090
d+01, 2.080
d+01, 2.070
d+01, 2.060
d+01,
1445 * 2.050
d+01, 2.045
d+01, 2.040
d+01, 2.030
d+01, 2.025
d+01,
1446 * 2.020
d+01, 2.015
d+01, 2.010
d+01, 2.005
d+01, 2.002
d+01,
1447 * 2.000
d+01, 1.999
d+01, 1.990
d+01, 1.985
d+01, 1.975
d+01,
1448 * 1.970
d+01, 1.965
d+01, 1.960
d+01, 1.950
d+01, 1.945
d+01,
1449 * 1.940
d+01, 1.925
d+01, 1.920
d+01, 1.915
d+01, 1.910
d+01,
1450 * 1.900
d+01, 1.898
d+01, 1.895
d+01, 1.890
d+01, 1.880
d+01,
1452 * 17.00
d+00, 14.40
d+00, 12.00
d+00, 11.00
d+00, 10.00
d+00/
1457 DATA (sena(ie,1),ie=1,nea) /
1458 * 5.103
d-01, 5.157
d-01, 5.103
d-01, 4.777
d-01, 4.072
d-01,
1459 * 3.420
d-01, 2.714
d-01, 1.683
d-01, 6.700
d-02, 6.100
d-02,
1460 * 5.800
d-02, 4.900
d-02, 3.800
d-02, 3.300
d-02, 3.000
d-02,
1461 * 2.400
d-02, 2.300
d-02, 2.900
d-02, 3.600
d-02, 4.100
d-02,
1462 * 4.000
d-02, 3.700
d-02, 3.400
d-02/
1465 DATA (sena(ie,2),ie=1,nea) /
1466 * 8.762
d-01, 8.856
d-01, 8.762
d-01, 8.203
d-01, 6.991
d-01,
1467 * 5.873
d-01, 4.661
d-01, 2.890
d-01, 1.401
d-01, 1.305
d-01,
1468 * 1.238
d-01, 1.069
d-01, 8.495
d-02, 7.480
d-02, 6.750
d-02,
1469 * 5.565
d-02, 5.230
d-02, 6.470
d-02, 7.765
d-02, 8.722
d-02,
1470 * 8.440
d-02, 7.821
d-02, 7.259
d-02/
1473 DATA (sena(ie,3),ie=1,nea) /
1474 * 9.200
d-01, 9.500
d-01, 9.400
d-01, 8.800
d-01, 7.500
d-01,
1475 * 6.100
d-01, 5.000
d-01, 3.700
d-01, 1.820
d-01, 1.710
d-01,
1476 * 1.620
d-01, 1.410
d-01, 1.130
d-01, 1.000
d-01, 9.000
d-02,
1477 * 7.500
d-02, 7.000
d-02, 8.600
d-02, 1.020
d-01, 1.140
d-01,
1478 * 1.100
d-01, 1.020
d-01, 9.500
d-02/
1481 DATA (sena(ie,4),ie=1,nea) /
1482 * 1.090
d+00, 1.180
d+00, 1.240
d+00, 1.280
d+00, 1.260
d+00,
1483 * 1.160
d+00, 9.300
d-01, 6.300
d-01, 3.580
d-01, 3.450
d-01,
1484 * 3.350
d-01, 2.990
d-01, 2.480
d-01, 2.220
d-01, 2.020
d-01,
1485 * 1.730
d-01, 1.610
d-01, 1.920
d-01, 2.200
d-01, 2.420
d-01,
1486 * 2.370
d-01, 2.220
d-01, 2.060
d-01/
1489 DATA (sena(ie,5),ie=1,nea) /
1490 * 1.029
d+00, 9.469
d-01, 1.091
d+00, 1.284
d+00, 1.591
d+00,
1491 * 1.691
d+00, 1.258
d+00, 9.241
d-01, 5.620
d-01, 5.493
d-01,
1492 * 5.375
d-01, 4.907
d-01, 4.182
d-01, 3.800
d-01, 3.484
d-01,
1493 * 3.038
d-01, 2.823
d-01, 3.307
d-01, 3.720
d-01, 4.040
d-01,
1494 * 3.959
d-01, 3.743
d-01, 3.517
d-01/
1497 DATA (sena(ie,6),ie=1,nea) /
1498 * 1.178
d+00, 9.793
d-01, 1.090
d+00, 1.271
d+00, 1.650
d+00,
1499 * 1.799
d+00, 1.339
d+00, 1.009
d+00, 6.223
d-01, 6.132
d-01,
1500 * 6.042
d-01, 5.572
d-01, 4.812
d-01, 4.402
d-01, 4.053
d-01,
1501 * 3.554
d-01, 3.304
d-01, 3.814
d-01, 4.244
d-01, 4.603
d-01,
1502 * 4.523
d-01, 4.293
d-01, 4.053
d-01/
1505 DATA (sena(ie,7),ie=1,nea) /
1506 * 1.386
d+00, 1.050
d+00, 1.134
d+00, 1.302
d+00, 1.722
d+00,
1507 * 1.922
d+00, 1.449
d+00, 1.103
d+00, 6.762
d-01, 6.686
d-01,
1508 * 6.602
d-01, 6.131
d-01, 5.344
d-01, 4.912
d-01, 4.541
d-01,
1509 * 4.004
d-01, 3.728
d-01, 4.273
d-01, 4.725
d-01, 5.103
d-01,
1510 * 5.022
d-01, 4.781
d-01, 4.524
d-01/
1513 DATA (sena(ie,8),ie=1,nea) /
1514 * 2.029
d+00, 1.537
d+00, 1.660
d+00, 1.906
d+00, 2.520
d+00,
1515 * 2.812
d+00, 2.121
d+00, 1.614
d+00, 1.014
d+00, 1.012
d+00,
1516 * 1.006
d+00, 9.557
d-01, 8.607
d-01, 8.038
d-01, 7.541
d-01,
1517 * 6.775
d-01, 6.334
d-01, 7.080
d-01, 7.669
d-01, 8.156
d-01,
1518 * 8.074
d-01, 7.769
d-01, 7.404
d-01/
1521 DATA (sena(ie,9),ie=1,nea) /
1522 * 3.050
d+00, 2.310
d+00, 2.495
d+00, 2.865
d+00, 3.789
d+00,
1523 * 4.228
d+00, 3.188
d+00, 2.426
d+00, 1.536
d+00, 1.538
d+00,
1524 * 1.536
d+00, 1.488
d+00, 1.384
d+00, 1.317
d+00, 1.256
d+00,
1525 * 1.153
d+00, 1.089
d+00, 1.185
d+00, 1.255
d+00, 1.315
d+00,
1526 * 1.307
d+00, 1.269
d+00, 1.224
d+00/
1529 DATA (sena(ie,10),ie=1,nea) /
1530 * 3.346
d+00, 2.535
d+00, 2.738
d+00, 3.143
d+00, 4.157
d+00,
1531 * 4.639
d+00, 3.498
d+00, 2.662
d+00, 1.685
d+00, 1.687
d+00,
1532 * 1.685
d+00, 1.632
d+00, 1.518
d+00, 1.445
d+00, 1.378
d+00,
1533 * 1.265
d+00, 1.194
d+00, 1.300
d+00, 1.377
d+00, 1.443
d+00,
1534 * 1.434
d+00, 1.392
d+00, 1.343
d+00/
1541 DATA (sepa(ie,1),ie=1,nea) /
1543 * 5.103
d-01, 5.157
d-01, 5.103
d-01, 4.777
d-01, 4.072
d-01,
1544 * 3.420
d-01, 2.714
d-01, 1.683
d-01, 6.700
d-02, 6.100
d-02,
1545 * 5.800
d-02, 4.900
d-02, 3.800
d-02, 3.300
d-02, 3.000
d-02,
1546 * 2.400
d-02, 2.300
d-02, 2.900
d-02, 3.600
d-02, 4.100
d-02,
1547 * 4.000
d-02, 3.700
d-02, 3.400
d-02/
1550 DATA (sepa(ie,2),ie=1,nea) /
1552 * 8.762
d-01, 8.856
d-01, 8.762
d-01, 8.203
d-01, 6.991
d-01,
1553 * 5.873
d-01, 4.661
d-01, 2.890
d-01, 1.401
d-01, 1.305
d-01,
1554 * 1.238
d-01, 1.069
d-01, 8.495
d-02, 7.480
d-02, 6.750
d-02,
1555 * 5.565
d-02, 5.230
d-02, 6.470
d-02, 7.765
d-02, 8.722
d-02,
1556 * 8.440
d-02, 7.821
d-02, 7.259
d-02/
1559 DATA (sepa(ie,3),ie=1,nea) /
1561 * 9.200
d-01, 9.500
d-01, 9.400
d-01, 8.800
d-01, 7.500
d-01,
1562 * 6.100
d-01, 5.000
d-01, 3.700
d-01, 1.820
d-01, 1.710
d-01,
1563 * 1.620
d-01, 1.410
d-01, 1.130
d-01, 1.000
d-01, 9.000
d-02,
1564 * 7.500
d-02, 7.000
d-02, 8.600
d-02, 1.020
d-01, 1.140
d-01,
1565 * 1.100
d-01, 1.020
d-01, 9.500
d-02/
1568 DATA (sepa(ie,4),ie=1,nea) /
1570 * 1.090
d+00, 1.180
d+00, 1.240
d+00, 1.280
d+00, 1.260
d+00,
1571 * 1.160
d+00, 9.300
d-01, 6.300
d-01, 3.650
d-01, 3.540
d-01,
1572 * 3.420
d-01, 3.060
d-01, 2.530
d-01, 2.260
d-01, 2.040
d-01,
1573 * 1.750
d-01, 1.610
d-01, 1.900
d-01, 2.200
d-01, 2.430
d-01,
1574 * 2.370
d-01, 2.220
d-01, 2.070
d-01/
1577 DATA (sepa(ie,5),ie=1,nea) /
1579 * 1.029
d+00, 9.469
d-01, 1.091
d+00, 1.284
d+00, 1.591
d+00,
1580 * 1.691
d+00, 1.258
d+00, 9.241
d-01, 5.828
d-01, 5.726
d-01,
1581 * 5.594
d-01, 5.100
d-01, 4.310
d-01, 3.897
d-01, 3.561
d-01,
1582 * 3.084
d-01, 2.829
d-01, 3.262
d-01, 3.714
d-01, 4.066
d-01,
1583 * 3.985
d-01, 3.764
d-01, 3.517
d-01/
1586 DATA (sepa(ie,6),ie=1,nea) /
1588 * 1.178
d+00, 9.793
d-01, 1.090
d+00, 1.271
d+00, 1.650
d+00,
1589 * 1.799
d+00, 1.339
d+00, 1.009
d+00, 6.383
d-01, 6.313
d-01,
1590 * 6.212
d-01, 5.732
d-01, 4.913
d-01, 4.483
d-01, 4.113
d-01,
1591 * 3.594
d-01, 3.304
d-01, 3.764
d-01, 4.243
d-01, 4.623
d-01,
1592 * 4.543
d-01, 4.313
d-01, 4.053
d-01/
1595 DATA (sepa(ie,7),ie=1,nea) /
1597 * 1.386
d+00, 1.050
d+00, 1.134
d+00, 1.302
d+00, 1.722
d+00,
1598 * 1.922
d+00, 1.449
d+00, 1.103
d+00, 6.950
d-01, 6.895
d-01,
1599 * 6.803
d-01, 6.322
d-01, 5.471
d-01, 5.014
d-01, 4.619
d-01,
1600 * 4.048
d-01, 3.728
d-01, 4.211
d-01, 4.722
d-01, 5.135
d-01,
1601 * 5.051
d-01, 4.804
d-01, 4.527
d-01/
1604 DATA (sepa(ie,8),ie=1,nea) /
1606 * 2.029
d+00, 1.537
d+00, 1.660
d+00, 1.906
d+00, 2.520
d+00,
1607 * 2.812
d+00, 2.121
d+00, 1.614
d+00, 1.045
d+00, 1.043
d+00,
1608 * 1.036
d+00, 9.718
d-01, 8.822
d-01, 8.211
d-01, 7.679
d-01,
1609 * 6.828
d-01, 6.325
d-01, 6.951
d-01, 7.647
d-01, 8.232
d-01,
1610 * 8.138
d-01, 7.935
d-01, 7.415
d-01/
1613 DATA (sepa(ie,9),ie=1,nea) /
1615 * 3.050
d+00, 2.310
d+00, 2.495
d+00, 2.865
d+00, 3.789
d+00,
1616 * 4.228
d+00, 3.188
d+00, 2.426
d+00, 1.589
d+00, 1.584
d+00,
1617 * 1.577
d+00, 1.528
d+00, 1.417
d+00, 1.345
d+00, 1.277
d+00,
1618 * 1.159
d+00, 1.086
d+00, 1.159
d+00, 1.252
d+00, 1.331
d+00,
1619 * 1.320
d+00, 1.278
d+00, 1.256
d+00/
1622 DATA (sepa(ie,10),ie=1,nea) /
1624 * 3.346
d+00, 2.535
d+00, 2.738
d+00, 3.143
d+00, 4.157
d+00,
1625 * 4.639
d+00, 3.498
d+00, 2.662
d+00, 1.743
d+00, 1.738
d+00,
1626 * 1.730
d+00, 1.676
d+00, 1.554
d+00, 1.475
d+00, 1.401
d+00,
1627 * 1.271
d+00, 1.191
d+00, 1.271
d+00, 1.373
d+00, 1.460
d+00,
1628 * 1.448
d+00, 1.402
d+00, 1.378
d+00/
1630 DATA relto / 0.175
d+00, 6*0.
d+00, 0.175
d+00, 4*0.
d+00, 0.14
d+00,
1638 IF(anuc.LT.1.5d0)
THEN
1642 IF(
ekin.GT.20.0d0)
THEN
1643 sigela=relto(kproj)*
dshpto(kproj,plab)
1645 ELSEIF(
ekin.GT.10.0d0)
THEN
1648 ek2=
sqrt(po2**2+aam(kproj)**2) - aam(kproj)
1649 sel2=relto(kproj)*
dshpto(kproj,po2)
1652 ELSE IF (
ekin .LT. 0.06
d+00 )
THEN
1653 IF ( kproj .EQ. 1 )
THEN
1654 IF (
ekin .LT. 0.02
d+00 )
THEN
1655 sigela = onethi * sigmlw(
ekin)
1657 betapr = plab / (
ekin + aam(kproj) )
1658 sigela = sigmpp(betapr)
1661 ELSE IF ( kproj .EQ. 8 )
THEN
1662 IF (
ekin .LT. 0.04
d+00 )
THEN
1663 sigela = sigmlw(
ekin)
1665 betapr = plab / (
ekin + aam(kproj) )
1666 sigela = sigmpn(betapr)
1674 IF(ek1.LT.ekihn(ie))
THEN
1679 ddee=ekihn(je2) - ekihn(je1)
1693 ELSEIF(kproj.EQ.8)
THEN
1697 ELSEIF(kproj.EQ.13)
THEN
1701 ELSEIF(kproj.EQ.14)
THEN
1710 sigela=s1 + (s2-s1)*(ek1-ekihn(je1))/ddee
1715 sigela=sel1 + (sel2-sel1)*(
ekin-ek1)/(ek2-ek1)
1724 IF(
ekin.LT.ekiha(ie))
THEN
1729 IF(
ekin.EQ.ekiha(nea))
THEN
1737 IF(anuc.LT.amass(ia))
THEN
1742 IF(anuc.EQ.amass(nnaa))
THEN
1750 70
IF (je) 190,150,80
1751 80 temp1=anuc/amass(ja)
1752 temp2=
log(amass(ja+1)/amass(ja))
1757 power=
log(sena(ke,ja+1)/slow)/temp2
1760 power=
log(sepa(ke,ja+1)/slow)/temp2
1762 tsig(i)=slow*temp1**power
1766 100 sigela=tsig(1) + (
ekin-ekiha(je))*(tsig(2)-tsig(1)) /(ekiha(je+1)
1774 temp1= (anuc/amass(ja)) **0.66667d0
1775 120
IF (je) 200,170,130
1776 130
IF(kproj.EQ.8)
THEN
1777 tsig(1) = sena(je,ja) * temp1
1778 tsig(2) = sena(je+1,ja) *temp1
1780 tsig(1) = sepa(je,ja) * temp1
1781 tsig(2) = sepa(je+1,ja) *temp1
1787 temp1= (anuc/amass(ja))**.66667
1792 160 temp1=anuc/amass(ja)
1793 temp2=
log(amass(ja+1)/amass(ja))
1796 power=
log(sena(je,ja+1)/slow)/temp2
1799 power=
log(sepa(je,ja+1)/slow)/temp2
1801 sigela=slow*temp1**power
1806 180
IF(kproj.EQ.8)
THEN
1807 sigela=sena(je,ja)*temp1
1809 sigela=sepa(je,ja)*temp1
1825 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1843 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1844 +iibar(210),k1(210),k2(210)
1848 umo2=amit2 + aam(1)**2 + 2.*aam(1)*(po+0.5*amit2/po)
1856 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1880 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1881 +iibar(210),k1(210),k2(210)
1888 CALL
sihnel(kpro,1,poo,siel)
1889 CALL
sihnin(kpro,1,poo,siin)
1897 IF (poo.LT.2100.)
THEN
1898 xshpto=45.6+219.*poo**(-4.23)+0.41*
log(poo)**2-3.41*
log(poo)
1900 ELSEIF(poo.LE.432000.)
THEN
1901 xshpto=41.1+77.2*poo**(-0.68)+0.293*
log(poo)**2-1.82*
log(poo)
1905 umo2=amit2 + aam(1)**2 + 2.*aam(1)*(po+0.5*amit2/po)
1907 IF (umo.LT.1800.)
THEN
1908 xshpto=61.65+(umo-900.)*6.6/900.
1910 ELSEIF(umo.LT.16000.)
THEN
1911 xshpto=68.3+(umo-1800.)*25./14200.
1914 xshpto=93.3+(umo-16000.)*12.5/24000.
1925 IF(poo.LT.10.0)
THEN
1927 goto(10,20,260,260,260,260,260,30,40,260,260,50,60,70,80,90,
1928 + 100,20,50, 10,10,10,110,120,130), itt
1970 IF(
rndm(v).GE.0.5) goto 80
2048 xshpto=a1 + a2*poo**a3 + a4*alp**2 + a5*alp
2053 ELSEIF(ipio.EQ.2)
THEN
2061 umo2=amit2 + aam(1)**2 + 2.*aam(1)*(po+0.5*amit2/po)
2064 goto(150,150,260,260,260,260,260,160,160,260,260,210,170,170,190,
2065 +190,240,240,210,240,240,240,250,220,230), itt
2072 IF(itt.EQ.1) goto 270
2084 IF(itt.EQ.8) goto 270
2091 IF(umo.LT.47.0) goto 180
2100 IF(itt.EQ.14) goto 270
2107 IF(umo.LT.110.) goto 200
2116 IF(itt.EQ.16) goto 270
2124 IF(
rndm(v).LT.0.5) itt=16
2166 xshpto=a1 + a2*(
log(umo2/a3))** 2+ a4/umo2 + a5*umo2**a6
2177 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2235 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2237 dimension seeg(4),pee(4),sigkp(20),sigkm(20),sigap(20)
2238 dimension seg(20),sig(20),segp(20),
p(20)
2240 DATA p/0.3d0,0.4d0,0.5d0,0.6d0,0.8d0,1.d0,
2241 +1.5d0,2.d0,3.d0,4.d0,5.d0,6.d0,10.d0,
2242 *20.d0,50.d0,100.d0,200.d0,400.d0,1000.d0,10000.d0/
2244 DATA seeg/0.1d0,16.d0,35.d0,42.d0/
2245 DATA pee /0.13d0,0.19d0,0.25d0,0.30d0/
2247 DATA (sig(ie),ie=1,20) / 3*0.0001d0,0.1d0,1.d0,4.d0,
2248 +24.5d0,25.d0,27.2d0,27.8d0,
2249 +28.5d0,29.2d0,29.7d0, 30.5d0,31.5d0,
2250 +31.7d0,32.1d0,32.9d0,34.5d0,41.2d0/
2253 DATA (seg(ie),ie=1,20) / 42.d0,19.d0,16.1d0,17.d0,
2254 +22.7d0,32.5d0,24.6d0,26.2d0,
2255 +25.0d0,23.7d0,23.0d0,22.5d0,22.d0,
2256 + 21.2d0,20.8d0,20.7d0,21.d0,21.9d0,23.8d0,28.4d0/
2259 DATA (segp(ie),ie=1,20) / 0.1d0,0.1d0,0.1d0,0.1d0,
2260 +11.d0,12.5d0,22.d0,19.2d0,21.5d0,
2261 +21.4d0,20.8d0,20.6d0,20.2d0, 19.8d0,19.9d0,
2262 +20.d0,20.5d0,21.5d0,23.5d0,32.2d0/
2265 DATA sigap/ 164.d0,126.d0,114.d0,98.d0,86.d0,
2266 +72.4d0,59.d0,57.d0,53.d0,52.d0,48.d0,45.5d0,
2267 +43.5d0, 40.4d0,36.5d0,35.2d0,34.5d0,34.5d0,
2271 DATA sigkm/ 38.d0,43.d0,23.d0,18.5d0,20.d0,
2272 +29.d0,25.d0,23.d0,22.5d0,21.d0,20.5d0,20.d0,
2273 +19.2d0, 18.5d0,17.8d0,17.8d0,18.3d0,19.2d0,
2277 DATA sigkp/ 0.001d0,0.001d0,0.001d0,0.001d0,
2278 +0.2d0,4.5d0,8.9d0,11.6d0,12.2d0,13.4d0,
2279 +13.6d0, 13.7d0,13.7d0,14.9d0,15.9d0,16.5d0,
2280 +17.4d0,18.6d0,20.9d0,28.8d0/
2286 IF(itar.NE.1.AND.itar.NE.8)
THEN
2287 IF(ipri.GE.1)
WRITE(6,
'(A/A,2I5,2(1PE15.6))')
2288 +
' WRONG CALL OF SIHNIN/ITAR',
2289 +
' IPROJ,ITAR,PO,SIIN :', iproj,itar,po,siin
2294 IF(iproj.GE.3.AND.iproj.LE.7)
RETURN
2297 IF(iproj.EQ.23) iipp=13
2301 IF(iitt.EQ.8) iitt=1
2303 IF(iproj.EQ.19.OR.iproj.EQ.12)
THEN
2306 IF(rnd.GT.0.5) iipp=25
2314 IF(po.LE.
p(jk))
THEN
2327 IF(iipp.EQ.1.OR.iipp.EQ. 8.OR.iipp.EQ.17.OR.(iipp.GE.20.and
2328 + .iipp.LE.22))
THEN
2334 ELSEIF(k.EQ.21)
THEN
2341 ELSEIF(iipp.EQ.14)
THEN
2346 IF(po.LE.pee(jk))
THEN
2362 siin=si1 + (po-pee(kkk))*(si2-si1)/(pee(kk)-pee(kkk))
2364 ELSEIF(k.EQ.21)
THEN
2372 ELSEIF(iipp.EQ.13)
THEN
2377 ELSEIF(k.EQ.21)
THEN
2384 ELSEIF(iipp.EQ.16.OR.iipp.EQ.25)
THEN
2388 ELSEIF(k.EQ.21)
THEN
2395 ELSEIF(iipp.EQ.15.OR.iipp.EQ.24)
THEN
2399 ELSEIF(k.EQ.21)
THEN
2407 ELSEIF(iipp.EQ.2.OR.iipp.EQ.9.OR.iipp.EQ.18)
THEN
2411 ELSEIF(k.EQ.21)
THEN
2423 60 siin=si1 + (po-
p(k))*(si1-si2)/(
p(k)-
p(kk))
2434 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2451 COMMON /dgamre/ redu,amo,amm(15 )
2452 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
2453 common/dreac/umo(296),plabf(296),siin(296),wk(5184),
2454 *nrk(2,268),nure(30,2)
2455 COMMON /dablti/ amh(110),gah(110),tauh(110),ichh(110),ibarh(110),
2457 COMMON /dspli/nzk(460,3),wt(460)
2458 COMMON /dmetls/ cxs(149),cys(149),
2459 +czs(149),els(149),pls(149),
2461 COMMON /drun/ runtes,eftes
2473 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
2474 +iibar(210),k1(210),k2(210)
2477 parameter(maxfin=10)
2478 COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin), czrh
2479 +(maxfin),elrh(maxfin),plrh(maxfin),irh
2481 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2483 dimension itprf(110)
2486 DATA itprf/-1,-1,5*1,-1,-1,1,1,1,-1,-1,-1,-1,6*1,-1,-1,-1,85*1/
2488 IF (
n.LE.0.OR.
n.GE.111)
n=1
2489 IF (itprf(
n ).GT.0 .OR. itta.GT.8)
THEN
2498 IF (abs(plab-5.).LT.4.99999d0) go to 20
2499 IF(ipri.GE.1)
WRITE (6,1010) plab
2501 1010
FORMAT (
' PROJECTILE HADRON MOMENTUM OUTSIDE OF THE
2502 + ALLOWED REGION, PLAB=',1e15.5)
2505 umodat=
n*1.11111d0+itta*2.19291d0
2506 IF(umodat.NE.umoda) CALL
dcalum(
n,itta)
2513 IF (lowp.GT.20) go to 280
2515 IF (nnn.EQ.
n) go to 50
2524 IF(itta.GT.1) ire=nure(
n,2)
2529 CALL
dsigin(ire,plab,
n,ie,amt ,amn,ecm,si,itta)
2531 IF (amh(1).NE.0.9383d0) ianth=1
2532 IF (ianth.GE.0) si=1.
2540 IF (si.LT.1.
d-6) go to 280
2541 IF (
n.LE.nstab) go to 60
2543 IF (runtes.LT.20.d0)
WRITE(6,1020)
n
2544 1020
FORMAT(3h
n=,i10,30h the proektile is
a resonance )
2545 IF(ibarh(
n).EQ.1)
n=8
2546 IF(ibarh(
n).EQ.-1)
n=9
2549 IF (imach.GT.10) go to 280
2553 ecmn=(ecm**2+amn2-amt2)/(2.*ecm )
2554 IF(ecmn.LE.amn) ecmn=amn
2555 pcmn=
sqrt(ecmn**2-amn2)
2558 IF (ianth.GE.0) ecm=2.1
2572 idwk=ieii(ire+1)-iiei
2580 humo=2.*umo(iiei+idwk)-umo(iiei+idwk-1)
2581 IF (humo.LT.ecm) ecm=humo
2598 iwk=iiwk+(ik-1)*idwk+ie-iiei
2607 IF (plab.LT.plabf(iiei+2)) wdk=0.
2608 wico=wok*1.23459876d0+wdk*1.735218469d0
2609 IF (wico.EQ.wicor) go to 70
2610 IF (umo(iiei+idwk).LT.hecm) wdk=0.
2616 eklim=-thresh(iiki+ik)
2618 delim=umo(ielim)+eklim
2620 dete=(ecm-(ecmo-eklim)*.5)*2.
2621 IF (delim*delim-dete*dete) 90,90,80
2626 wkk=wok-wdk*dec/(decc+1.
d-9)
2632 IF (vv.GT.wkk) go to 70
2647 IF (i1001.GT.50) go to 60
2649 IF (it2*ams.GT.it2*ecm) go to 110
2652 IF (ianth.GE.0) ecm=elab+amt+0.00001d0
2655 IF (it2.GT.0) go to 120
2673 IF(ww.LT. 0.5d0) go to 130
2688 IF(ib1.EQ.ibn) go to 140
2698 CALL
dtwopa(ecm1,ecm2,pcm1,pcm2,cod1,cod2,cof1,cof2,sif1,sif2,
2699 *it1,it2,ecm,ecmn,pcmn,
n,am1,am2)
2707 CALL
dtrafo(gam,bgam,cx,cy,cz,cod1,cof1,sif1,pcm1,ecm1,pls(ist),
2708 *cxs(ist),cys(ist),czs(ist),els(ist))
2712 CALL
dtrafo(gam,bgam,cx,cy,cz,cod2,cof2,sif2,
2713 *pcm2,ecm2,pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
2719 IF(its(ist).GT.nstab) go to 160
2733 IF(ist.GE.1) go to 150
2752 IF (vv.GT.wt(iik)) go to 180
2763 IF (it2-1.LT.0) go to 240
2770 IF (ieco.LE.10) go to 200
2772 IF(iatmpt.GT.3) go to 280
2775 IF (i310.GT.50) go to 170
2776 IF (ams.GT.eco) go to 190
2781 IF (redu.LT.0.d0) go to 30
2784 IF(it3.EQ.0) go to 220
2787 CALL
dthrep(eco,ecm1,ecm2,ecm3,pcm1,pcm2,pcm3,cod1,cof1,sif1,
2788 *cod2,cof2,sif2,cod3,cof3,sif3,am1,am2,am3)
2790 220 CALL
dtwopd(eco,ecm1,ecm2,pcm1,pcm2,cod1,cof1,sif1,cod2,cof2,sif2,
2796 IF (redu.GT.0.d0) go to 240
2798 IF (itwthc.GT.100) go to 30
2799 IF (itwth) 220,220,210
2802 IF (it2-1.LT.0) go to 250
2809 CALL
dtrafo(gam,bgam,rx,ry,rz,cod1,cof1,sif1,pcm1,ecm1,
2810 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
2813 CALL
dtrafo(gam,bgam,rx,ry,rz,cod2,cof2,sif2,pcm2,ecm2,
2814 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
2815 IF (it3.LE.0) go to 250
2818 CALL
dtrafo(gam,bgam,rx,ry,rz,cod3,cof3,sif3,pcm3,ecm3,
2819 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
2844 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2846 DATA runt/100.d0,100.d0/
2850 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2862 IF(
n.EQ.2*nn) go to 10
2869 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2872 COMMON /dslope/ ambmbb(75)
2873 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
2876 DATA ikii/ 0, 15, 41, 67, 82, 93, 110, 133, 148, 159, 172, 183,
2877 & 207, 224, 241, 252, 268 /
2878 DATA ieii/ 0, 21, 46, 71, 92, 109, 126, 143, 160, 173, 186, 199,
2879 & 220, 241, 262, 279, 296 /
2880 DATA irii/ 0, 315, 965, 1615, 1930, 2117, 2406, 2797, 3052, 3195,
2881 & 3364, 3507, 4011, 4368, 4725, 4912, 5184/
2889 DATA ambmbb/ 0.8d0, 0.85d0, 0.9d0, 0.95d0, 1.d0,
2890 & 1.05d0, 1.1d0, 1.15d0, 1.2d0, 1.25d0,
2891 & 1.3d0, 1.35d0, 1.4d0, 1.45d0, 1.5d0,
2892 & 1.55d0, 1.6d0, 1.65d0, 1.7d0, 1.75d0,
2893 & 1.8d0, 1.85d0, 1.9d0, 1.95d0, 2.d0,
2894 & 15.6d0, 14.95d0, 14.3d0, 13.65d0, 13.d0,
2895 & 12.35d0, 11.7d0, 10.85d0, 10.d0, 9.15d0,
2896 & 8.3d0, 7.8d0, 7.3d0, 7.25d0, 7.2d0,
2897 & 6.95d0, 6.7d0, 6.6d0, 6.5d0, 6.3d0,
2898 & 6.1d0, 5.85d0, 5.6d0, 5.35d0, 5.1d0,
2899 & 15.d0, 15.d0, 15.d0, 15.d0, 15.d0, 15.d0, 15.d0,
2900 & 14.2d0, 13.4d0, 12.6d0,
2901 & 11.8d0, 11.2d0, 10.6d0, 9.8d0, 9.d0,
2902 & 8.25d0, 7.5d0, 6.25d0, 5.d0, 4.5d0, 5*4.d0 /
2908 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2911 common/dablti/amm(110),ga(110),
tau(110),ich(110)
2912 *,ibar(110),k1(110),k2(110)
2913 COMMON /dgamre/ redu,amo,am(15 )
2914 dimension gasuni(14)
2916 *-1.d0,-.98d0,-.95d0,-.87d0,-.72d0,-.48d0,
2917 *-.17d0,.17d0,.48d0,.72d0,.87d0,.95d0,.98d0,1.d0/
2923 IF (it.LE.0) go to 30
2924 IF (it.LE.nstab) go to 20
2925 dgauni=gauno*gaunon/(io-1.)
2932 IF (vv.GT.v1) go to 10
2933 uniga=dgauni*(i-2.+(vv-vo+1.
e-16)/(v1-vo)-(io-1.)*.5)
2934 dam=ga(it)*uniga/gauno
2947 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2951 common/dablti/am(110),ga(110),
tau(110),ich(110)
2952 *,ibar(110),k1(110),k2(110)
2953 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
2954 common/dspli/nzk(460,3),wt(460)
2955 COMMON /dreac/umo( 296),plabf( 296),siin( 296),wk( 5184),
2956 *nrk(2, 268),nure(30,2)
2957 ire=nure(
n,itta/8+1)
2966 elab=
sqrt(am12+am22+2.*
sqrt(plab2+am12)*am2)
2973 IF(nrk(2,ik).GT.0) go to 30
2982 IF(
in.GT.0)ams=ams+am(
in)
2984 IF(
in.GT.0) ams=ams+am(
in)
2985 IF (ams.LT.amss) amss=ams
2987 IF(umoo.LT.amss) umoo=amss
2994 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2995 common/dablti/am(110),ga(110),
tau(110),ich(110)
2996 *,ibar(110),k1(110),k2(110)
2997 common/dspli/ nzk(460,3),wt(460)
2999 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
3001 COMMON /dreac/umo( 296),plabf( 296),siin( 296),wk( 5184),
3002 *nrk(2, 268),nure(30,2)
3004 equivalence(wk(1),si(1))
3015 iee=ieii(ire+1)-ieii(ire)
3016 ike=ikii(ire+1)-ikii(ire)
3024 iwk=iwko+iee*(ik-1)+ie
3025 IF(nrk(2,iika+ik).EQ.0) sinorc=1.0
3026 sis=sis+si(iwk)*sinorc
3030 IF (sis.GE.1.
d-12) go to 20
3036 iwk=iwko+iee*(ik-1)+ie
3037 IF(nrk(2,iika+ik).EQ.0) sinorc=1.0
3038 sio=sio+si(iwk)*sinorc/sis
3042 iwk=iwko+iee*(ik-1)+ie
3047 inrk1=nrk(1,iiki+ik)
3048 IF (inrk1.GT.0) am111=am(inrk1)
3050 inrk2=nrk(2,iiki+ik)
3051 IF (inrk2.GT.0) am222=am(inrk2)
3052 thresh(iiki+ik)=am111 +am222
3053 IF (inrk2-1.GE.0) go to 60
3057 DO 50 inrk1=inrkk,inrko
3061 IF (inzk1.LE.0.OR.inzk1.GT.110) go to 50
3062 IF (inzk2.LE.0.OR.inzk2.GT.110) go to 50
3063 IF (inzk3.LE.0.OR.inzk3.GT.110) go to 50
3066 ams=am(inzk1)+am(inzk2)
3067 IF (inzk3-1.GE.0) ams=ams+am(inzk3)
3068 IF (amss.GT.ams) amss=ams
3071 IF (ams.LT.umo(ieo)) ams=umo(ieo)
3083 IF (ik2.GT.460)ik2=460
3099 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3100 CHARACTER*8 anamc,zknamc,aname,zkname,rkname,anamz,zknamz
3102 common/dpar/anamc(210),amc(210),gac(210),tauc(210),ichc(210),
3103 *ibarc(210),k1c(210),k2c(210)
3106 parameter(idmax9=602)
3108 common/ddecac/ zknamc(idmax9),wtc(idmax9),nzkc(idmax9,3)
3111 COMMON /dnams/ aname(110),zkname(460),rkname(268)
3112 COMMON /dablti/am(110),ga(110),
tau(110),ich(110)
3113 *,ibar(110),k1(110),k2(110)
3114 COMMON /dspli/nzk(460,3),wt(460)
3115 COMMON /daddhp/ amz(16),gaz(16),tauz(16),ichz(16),ibarz(16),k1z
3116 +(16),k2z(16),wtz(153),ii22, nzkz(153,3)
3118 common/daddhn/anamz(16),zknamz(153)
3122 IF (iretur.GT.1)
RETURN
3142 nzk( i,1)=nzk
c(i, 1)
3143 nzk( i,2)=nzk
c(i, 2)
3144 nzk( i,3)=nzk
c(i, 3)
3161 nzk( l,3)=nzk
z(i, 3)
3162 nzk( l,2)=nzk
z(i, 2)
3163 nzk( l,1)=nzk
z(i, 1)
3169 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3173 COMMON /drun/ runtes,eftes
3174 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
3175 COMMON /dreac/umo( 296),plabf( 296),siin( 296),wk( 5184),
3176 *nrk(2, 268),nure(30,2)
3180 IF (pl.LT.0.) go to 30
3183 IF (pl.LE.plabf(i)) go to 60
3186 IF ( eftes.GT.40.d0) go to 20
3194 IF (-pl.LE.umo(i)) go to 60
3197 IF ( eftes.GT.40.d0) go to 50
3205 1000
FORMAT(14h plab or -ecm=,e12.4,27h is out of considered range ,
3209 SUBROUTINE dsigin(IRE ,PLAB,N,IE ,AMT ,AMN,ECM ,SI ,ITAR)
3210 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3211 common/dablti/am(110),ga(110),
tau(110),ich(110)
3212 *,ibar(110),k1(110),k2(110)
3213 COMMON /dredve/ thresh(268), irii(17),ikii(17),ieii(17)
3214 COMMON /dreac/umo( 296),plabf( 296),siin( 296),wk( 5184),
3215 *nrk(2, 268),nure(30,2)
3217 IF (ie.LE.ieii(ire)) ie=ie+1
3222 ecm=
sqrt(amn2+amt2+2.*amt*
sqrt(amn2+plab**2))
3232 IF (ecm.GT.ecmo) wdk=0.
3235 delim=umo(ielim)+eklim
3237 dete=(ecm-(ecmo-eklim)*.5)*2.
3238 IF (delim*delim-dete*dete) 20,20,10
3243 wkk=wok-wdk*dec/(decc+1.
d-9)
3244 IF (wkk.LT.0.) wkk=0.
3246 IF (-eklim.GT.ecm) si=1.
d-14
3250 SUBROUTINE dtchoi(T,P,PP,E,EE,I,II,N,AM1,AM2)
3251 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3256 COMMON /dablti/ am(110),ga(110),
tau(110),ich(110), ibar(110),k1
3258 COMMON /dslope/ sm(25),bbm(25),bbb(25)
3261 IF (i.GT.30.AND.ii.GT.30) go to 20
3264 IF (i.LE.30) go to 10
3272 IF (ama.LE.amb) go to 30
3287 tmin=-2.*(
e*ee-
p*pp)+am(
n)**2+am1 **2
3288 tmax=-2.*(
e*ee+
p*pp)+am(
n)**2+am1 **2
3298 IF (abs(tma).GT.120.d0) go to 70
3301 an=(1./bm)*(
exp(tmi)-etma)
3304 t=(1./bm)*
log(etma+
r*an*bm)
3308 SUBROUTINE dtwopa(E1,E2,P1,P2,COD1,COD2,COF1,COF2,SIF1,SIF2,
3309 1it1,it2,umoo,ecm,
p,
n,am1,am2)
3310 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3319 COMMON /dablti/ am(110),ga(110),
tau(110),ich(110), ibar(110),k1
3325 e1=((umoo-amb)*(umoo+amb) + ama2)/(2.*umoo)
3327 IF (e1.LT.ama*1.00001d0) e1=ama*1.00001d0
3328 amte=(e1-ama)*(e1+ama)
3338 CALL
dcosi(cof1,sif1)
3342 CALL
dtchoi(tr,
p,
p1,ecm,e1,it1,it2,
n,am1,am2)
3343 cod1=(tr-ama2-am(
n)*am(
n)+2.*ecm*e1)/(2.*
p*
p1+1.
e-18)
3344 IF (cod1.GT.0.9999999d0) cod1=0.9999999d0
3350 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3356 IF(
xy.GT.1.) go to 10
3359 IF(
rndm(
c).LT.0.5d0) go to 20
3366 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3368 IF(is.NE.0) go to 10
3382 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3392 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3393 CHARACTER*8 zknam5,zknam6,anamz,zknam4
3394 COMMON /daddhp/ amz(16),gaz(16),tauz(16),ichz(16),ibarz(16),
3395 & k1z(16),k2z(16),wtz(153),ii22,nzk1(153),nzk2(153),nzk3(153)
3396 COMMON /daddhn/ anamz(16),zknam4(9),zknam5(90),zknam6(54)
3398 DATA amz/ 3*2.2d0, 0.9576d0, 3*1.887d0, 2.4d0, 2.03d0, 2*1.44d0,
3401 DATA gaz/ 3*.2d0, .1d0, 4*.2d0, .18d0, 2*.2d0, 2*.15d0, 3*0.d0 /
3403 DATA tauz / 16*0.d0 /
3405 DATA ichz/ 0, 1, 3*0, 1, -1, 0, 1, -1, 0, 0, 1 , 3*0 /
3407 DATA ibarz/ 2, 7*0, 1, -1, -1, 1, 1, 3*0 /
3410 DATA k1z/ 308,310,313,317,322,365,393,421,425,434,440,446,449,
3414 DATA k2z/ 309,312,316,321,364,392,420,424,433,439,445,448,451,
3417 DATA wtz/ .17d0, .83d0, 2*.33d0, .34d0, .17d0, 2*.33d0, .17d0,
3418 & .01d0, .13d0, .36d0, .27d0, .23d0, .0014d0, .0029d0, .0014d0,
3419 & .0029d0, 4*.0007d0, .0517d0, .0718d0, .0144d0, .0431d0, .0359d0,
3420 & .0718d0, .0014d0, .0273d0, .0014d0, .0431d0, 2*.0129d0, .0259d0,
3421 & .0517d0, .0359d0, .0014d0, 2*.0144d0, .0129d0, .0014d0, .0259d0,
3422 & .0359d0, .0072d0, .0474d0, .0948d0, .0259d0, .0072d0, .0144d0,
3423 & .0287d0, .0431d0, .0144d0, .0287d0, .0474d0, .0144d0, .0075d0,
3424 & .0057d0, .0019d0, .0038d0, .0095d0, 2*.0014d0, .0191d0, .0572d0,
3425 & .1430d0, 2*.0029d0, 5*.0477d0, .0019d0, .0191d0, .0686d0,.0172d0,
3426 & .0095d0, .1888d0, .0172d0, .0191d0, .0381d0, 2*.0571d0, .0190d0,
3427 & .0057d0, .0019d0, .0038d0, .0095d0, .0014d0, .0014d0, .0191d0,
3428 & .0572d0, .1430d0, 2*.0029d0, 5*.0477d0, .0019d0, .0191d0,.0686d0,
3429 & .0172d0, .0095d0, .1888d0, .0172d0, .0191d0, .0381d0, 2*.0571d0,
3430 & .0190d0, 4*.25d0, 2*.2d0, .12d0, .1d0, .07d0, .07d0, .14d0,
3431 & 2*.05d0, .0d0, .3334d0, .2083d0, 2*.125d0, .2083d0, .0d0, .125d0,
3432 & .2083d0, .3334d0, .2083d0, .125d0, .3d0, .05d0, .65d0, .3d0,
3433 & .05d0, .65d0, 9*1.d0 /
3435 DATA nzk1/ 8, 1, 2, 9, 1, 2, 9, 2, 9, 7, 13, 31, 15, 24, 23, 13,
3436 & 23, 13, 2*23, 14, 13, 23, 31, 98, 2*33, 32, 23, 14, 13, 35, 2*23,
3437 & 14, 13, 33, 23, 98, 31, 23, 14, 13, 35, 2*33, 32, 23, 35, 33, 32,
3438 & 98, 5*35, 4*13, 23, 13, 98, 32, 33, 23, 13, 23, 13, 14, 13, 32,
3439 & 13, 98, 23, 13, 2*32, 13, 33, 32, 98, 2*35, 4*14, 23, 14, 98,
3440 & 2*34, 23, 14, 23, 2*14, 13, 34, 14, 98, 23, 14, 2*34, 14, 33, 32,
3441 & 98, 2*35, 104, 61, 105, 62, 1, 17, 21, 17, 22, 2*21, 22, 21, 2,
3442 & 67, 68, 69, 2, 2*9, 68, 69, 70, 2, 9, 2*24, 15, 2*25, 16, 9*0/
3443 DATA nzk2/ 2*8, 1, 8, 9, 2*8, 2*1, 7, 14, 13, 16, 25, 23, 14, 23,
3444 & 14, 31, 33, 32, 34, 35, 31, 23, 31, 33, 34, 31, 32, 34, 31, 33,
3445 & 32, 2*33, 35, 31, 33, 31, 33, 32, 34, 35, 31, 33, 34, 35, 31,
3446 & 4*33, 32, 3*35, 2*23, 13, 31, 32, 33, 13, 31, 32, 2*31, 32, 33,
3447 & 32, 32, 35, 31, 2*32, 33, 31, 33, 35, 33, 3*32, 35, 2*23, 14,
3448 & 31, 34, 33, 14, 31, 33, 2*31, 34, 32, 33, 34, 35, 31, 2*34, 33,
3449 & 31, 33, 35, 33, 2*34, 33, 35, 1, 2, 8, 9, 25, 13, 35, 2*32, 33,
3450 & 31, 13, 23, 31, 13, 23, 14, 79, 80, 31, 13, 23, 14, 78, 79, 8,
3451 & 1, 8, 1, 8, 1, 9*0 /
3452 DATA nzk3/ 23, 14, 2*13, 23, 13, 2*23, 14, 0, 7, 14, 4*0, 2*23,
3453 & 10*0, 33, 2*31, 0, 33, 34, 32, 34, 0, 35, 0, 31, 3*35, 0, 3*31,
3454 & 35, 31, 33, 34, 31, 33, 34, 31, 33, 35, 0, 23, 14, 6*0, 32, 3*33,
3455 & 32, 34, 0, 35, 0, 2*35, 2*31, 35, 32, 34, 31, 33, 32, 0, 23, 13,
3456 & 6*0, 34, 2*33, 34, 33, 34, 0, 35, 0,2*35, 2*31, 35, 2*34, 31,
3457 & 2*34, 25*0, 23, 2*14, 23, 2*13, 9*0 /
3459 DATA anamz /
'NNPI',
'ANPPI',
'ANNPI',
' ETS ',
' PAP ',
' PAN ',
3460 &
'APN',
'DEO ',
'S+2030',
'AN*-14',
'AN*014',
'KONPI ',
'AKOPPI',
3463 DATA zknam4/
'NNPI0',
'PNPI-',
'APPPI+',
'ANNPI+',
'ANPPI0',
'APNPI+',
3464 &
'ANNPI0',
'APPPI0',
'ANPPI-'/
3465 DATA zknam5/
' GAGA ',
'P+P-GA',
'ETP+P-',
'K+K- ',
'K0AK0 ',
3466 &
' POPO ',
' P+P- ',
'POPOPO',
'P+P0P-',
'P0ET ',
'&0R0 ',
'P-R+ ',
3467 &
'P+R- ',
'POOM ',
' ETET ',
'ETSP0 ',
'R0ET ',
' R0R0 ',
'R+R- ',
3468 &
'P0ETR0',
'P-ETR+',
'P+ETR-',
' OMET ',
'P0R0R0',
'P0R+R-',
'P-R+R0',
3469 &
'P+R-R0',
'R0OM ',
'P0ETOM',
'ETSR0 ',
'ETETET',
'P0R0OM',
'P-R+OM',
3470 &
'P+R-OM',
'OMOM ',
'R0ETET',
'R0R0ET',
'R+R-ET',
'P0OMOM',
'OMETET',
3471 &
'R0R0R0',
'R+R0R-',
'ETSRET',
'OMR0R0',
'OMR+R-',
'OMOMET',
'OMOMR0',
3473 &
' P+PO ',
'P+POPO',
'P+P+P-',
'P+ET ',
'P0R+ ',
'P+R0 ',
'ETSP+ ',
3474 &
'R+ET ',
' R0R+ ',
'POETR+',
'P+ETR0',
'POR+R-',
'P+R0R0',
'P-R+R+',
3475 &
'P+R-R+',
'R+OM ',
'P+ETOM',
'ETSR+ ',
'POR+OM',
'P+R0OM',
'R+ETET',
3476 &
'R+R0ET',
'P+OMOM',
'R0R0R+',
'R+R+R-',
'ETSR+E',
'OMR+R0',
'OMOMR+',
3477 &
'P-PO ',
'P-POPO',
'P-P-P+',
'P-ET ',
'POR- ',
'P-R0 ',
'ETSP- ',
3478 &
'R-ET ',
'R-R0 ',
'POETR-',
'P-ETR0',
'POR-R0',
'P-R+R-',
'P-R0R0'/
3479 DATA zknam6/
'P+R-R-',
'R-OM ',
'P-ETOM',
'ETSR- ',
'POR-OM',
'P-R0OM',
3480 &
'R-ETET',
'R-R0ET',
'P-OMOM',
'R0R0R-',
'R+R-R-',
'ETSR-E',
'OMR0R-',
3481 &
'OMOMR-',
'PAN-14',
'APN+14',
'NAN014',
'ANN014',
'PAKO ',
'LPI+ ',
3482 &
'SI+OM',
'LAMRO+',
'SI0RO+',
'SI+RO0',
'SI+ETA',
'SI0PI+',
'SI+PI0',
3483 &
'APETA ',
'AN=P+ ',
'AN-PO ',
'ANOPO ',
'APRHOO',
'ANRHO-',
'ANETA ',
3484 &
'AN-P+ ',
'AN0PO ',
'AN+P- ',
'APRHO+',
'ANRHO0',
3485 &
'KONPIO',
'KOPPI-',
'K+NPI-',
'AKOPPO',
'AKONP+',
'K-PPI+',
3500 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3523 COMMON /dreac/ umo( 296), plabf( 296), siin( 296),
3524 + wk(5184), nrk(2,268), nure(30,2)
3527 & umopi(92), umokc(68), umop(39), umon(63), umok0(34),
3528 & plapi(92), plakc(68), plap(39), plan(63), plak0(34),
3529 & spikp1(315), spikpu(278), spikpv(372),
3530 & spikpw(278), spikpx(372), spikp4(315),
3531 & spikp5(187), spikp6(289),
3532 & skmpel(102), spikp7(289), skmnel(68), spikp8(187),
3533 & spikp9(143), spikp0(169), spkpv(143),
3534 & sappel(105), spikpe(399), sapnel(84), spikpz(273),
3535 & sanpel(84) , spikpf(273),
3536 & spkp15(187), spkp16(272),
3537 & nrkpi(164), nrkkc(132), nrkp(70), nrkn(116), nrkk0(54),
3540 dimension nrklin(532)
3541 equivalence(nrk(1,1), nrklin(1))
3542 equivalence( umo( 1), umopi(1)), ( umo( 93), umokc(1))
3543 equivalence( umo(161), umop(1)), ( umo(200), umon(1))
3544 equivalence( umo(263), umok0(1))
3545 equivalence( plabf( 1), plapi(1)), ( plabf( 93), plakc(1))
3546 equivalence( plabf(161), plap(1)), ( plabf(200), plan(1))
3547 equivalence( plabf(263), plak0(1))
3548 equivalence( wk( 1), spikp1(1)), ( wk( 316), spikpu(1))
3549 equivalence( wk( 594), spikpv(1)), ( wk( 966), spikpw(1))
3550 equivalence( wk(1244), spikpx(1)), ( wk(1616), spikp4(1))
3551 equivalence( wk(1931), spikp5(1)), ( wk(2118), spikp6(1))
3552 equivalence( wk(2407), skmpel(1)), ( wk(2509), spikp7(1))
3553 equivalence( wk(2798), skmnel(1)), ( wk(2866), spikp8(1))
3554 equivalence( wk(3053), spikp9(1)), ( wk(3196), spikp0(1))
3555 equivalence( wk(3365), spkpv(1)), ( wk(3508), sappel(1))
3556 equivalence( wk(3613), spikpe(1)), ( wk(4012), sapnel(1))
3557 equivalence( wk(4096), spikpz(1)), ( wk(4369), sanpel(1))
3558 equivalence( wk(4453), spikpf(1)), ( wk(4726), spkp15(1))
3559 equivalence( wk(4913), spkp16(1))
3560 equivalence(nrk(1,1), nrklin(1))
3561 equivalence(nrklin( 1), nrkpi(1)), (nrklin( 165), nrkkc(1))
3562 equivalence(nrklin( 297), nrkp(1)), (nrklin( 367), nrkn(1))
3563 equivalence(nrklin( 483), nrkk0(1))
3564 equivalence(nure(1,1), nureln(1))
3568 DATA plapi / 0.d0, .3d0, .5d0, .6d0, .7d0, .8d0, .9d0, .95d0,1.d0,
3569 & 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0, 2.d0, 2.3d0, 2.5d0, 2.8d0,
3570 & 3.d0, 3.5d0, 4.d0, 0.d0, .285d0, .4d0, .45d0, .5d0, .6d0, .7d0,
3571 & .75d0, .8d0, .85d0, .9d0, 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0,
3572 & 1.8d0, 2.d0, 2.3d0, 2.5d0, 2.8d0, 3.d0, 3.5d0, 4.d0, 4.5d0, 0.d0,
3573 & .285d0, .4d0, .45d0, .5d0, .6d0, .7d0, .75d0, .8d0, .85d0, .9d0,
3574 & 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0, 2.d0, 2.3d0, 2.5d0,
3575 & 2.8d0, 3.d0, 3.5d0, 4.d0, 4.5d0, 0.d0, .3d0, .5d0, .6d0, .7d0,
3576 & .8d0, .9d0, .95d0, 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0,
3577 & 2.d0, 2.3d0, 2.5d0, 2.8d0, 3.d0, 3.5d0, 4.d0 /
3579 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3580 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3581 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
3582 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3583 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3584 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
3585 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3586 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3587 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
3588 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3589 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3590 & 3.51d0, 3.84d0, 4.16d0, 4.49d0/
3592 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3593 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3594 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
3595 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
3596 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
3597 & 3.51d0, 3.84d0, 4.16d0, 4.49d0/
3600 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3601 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
3602 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3603 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
3604 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3605 & 3.43d0, 3.75d0, 4.07d0, 4.43d0 /
3608 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
3609 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3610 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
3611 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
3612 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3613 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
3614 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
3615 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
3616 & 3.43d0, 3.75d0, 4.07d0, 4.43d0 /
3617 DATA siin / 296*0.d0 /
3618 DATA umopi/ 1.08d0,1.233d0,1.302d0,1.369d0,1.496d0,
3619 & 1.557d0,1.615d0,1.6435d0,
3620 & 1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,2.071d0,2.159d0,
3621 & 2.286d0,2.366d0,2.482d0,2.56d0,
3623 & 1.08d0,1.222d0,1.302d0,1.3365d0,1.369d0,1.434d0,
3624 & 1.496d0,1.527d0,1.557d0,
3625 & 1.586d0,1.615d0,1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,
3626 & 2.071d0,2.159d0,2.286d0,2.366d0,
3627 & 2.482d0,2.560d0,2.735d0,2.90d0,3.06d0,
3628 & 1.08d0,1.222d0,1.302d0,1.3365d0,1.369d0,1.434d0,
3629 & 1.496d0,1.527d0,1.557d0,
3630 & 1.586d0,1.615d0,1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,
3631 & 2.071d0,2.159d0,2.286d0,2.366d0,
3632 & 2.482d0,2.560d0,2.735d0,2.90d0,3.06d0,
3633 & 1.08d0,1.233d0,1.302d0,1.369d0,1.496d0,
3634 & 1.557d0,1.615d0,1.6435d0,
3635 & 1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,2.071d0,2.159d0,
3636 & 2.286d0,2.366d0,2.482d0,2.56d0,
3639 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3640 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3642 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3643 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3645 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3646 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3648 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3649 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3652 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3653 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3655 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
3656 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
3660 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3662 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3664 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3668 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
3669 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3671 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
3672 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3674 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
3675 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
3678 DATA nrkpi / 13, 1, 15, 21, 81, 0, 13, 54, 23, 53, 13, 63, 13, 58,
3679 & 23, 57, 13, 65, 1, 32, 53, 31, 54, 32, 53, 33, 53, 35, 63, 32,
3680 & 13, 8, 23, 1, 17, 15, 21, 24, 22, 15, 82, 0, 61, 0, 13, 55, 23,
3681 & 54, 14, 53, 13, 64, 23, 63, 13, 59, 23, 58, 14, 57, 13, 66, 23,
3682 & 65, 1, 31, 8, 32, 1, 33, 1, 35, 54, 31, 55, 32, 54, 33, 53, 34,
3683 & 54, 35, 14, 1, 23, 8, 17, 24, 20, 15, 22, 24, 83, 0, 62, 0, 14,
3684 & 54, 23, 55, 13, 56, 14, 63, 23, 64, 14, 58, 23, 59, 13, 60, 14,
3685 & 65, 23, 66, 8, 31, 1, 34, 8, 33, 8, 35, 55, 31, 54, 34, 55, 33,
3686 & 56, 32, 55, 35, 14, 8, 24, 20, 84, 0, 14, 55, 23, 56, 14, 64, 14,
3687 & 59, 23, 60, 14, 66, 8, 34, 56, 31, 55, 34, 56, 33, 56, 35, 64,34/
3688 DATA nrkkc/ 15, 1, 89, 0, 24, 53, 15, 54, 1, 36, 1, 40, 1, 44, 36,
3689 & 63, 15, 63, 45, 53, 44, 54, 15, 8, 24, 1, 91, 0, 24, 54, 15, 55,
3690 & 8, 36, 1, 37, 8, 40, 1, 41, 8, 44, 1, 45, 36, 64, 37, 63, 15, 64,
3691 & 24, 63, 45, 54, 44, 55, 16, 1, 25, 8, 17, 23, 21, 14, 20,
3692 & 13, 22, 23, 90, 0, 38, 1, 39, 8, 16, 54, 25, 55, 1, 42, 8, 43,
3693 & 16, 63, 25, 64, 39, 64, 38, 63, 46, 54, 47, 55, 8, 47, 1, 46, 52,
3694 & 0, 51, 0, 16, 8, 17, 14, 20, 23, 22, 14, 92, 0, 8, 38, 16, 55,
3695 & 25, 56, 8, 42, 16, 64, 38, 64, 46, 55, 47, 56, 8, 46, 94, 0 /
3699 DATA nrkk0 / 24, 8, 106, 0, 15, 56, 24, 55, 37, 8, 41, 8, 45, 8,
3700 & 37, 64, 24, 64, 44, 56, 45, 55, 25, 1, 17, 13, 22, 13, 21, 23,
3701 & 107, 0, 39, 1, 25, 54, 16, 53, 43, 1, 25, 63, 39, 63, 47, 54, 46,
3702 & 53, 47, 1, 103, 0, 93, 0/
3704 DATA nrkp / 1, 1, 85, 0, 8, 53, 1, 54, 1, 63, 8, 57, 1, 58, 2*54,
3705 & 53, 55, 63, 54, 64, 53, 1, 8, 86, 0, 8, 54, 1, 55, 8, 63, 1, 64,
3706 & 8, 58, 1, 59, 64, 54, 63, 55, 54, 55, 53, 56, 77, 0, 2*8, 95, 0,
3707 & 8, 55, 1, 56, 8, 64, 8, 59, 1, 60, 2*55, 54, 56, 64, 55, 63, 56 /
3709 DATA nrkn/ 1, 2, 17, 18, 15, 16, 8, 9, 13, 14, 99, 0, 87, 0, 1,
3710 & 68, 8, 69, 2, 54, 9, 55, 102, 0, 2, 63, 9, 64, 1, 75, 8, 76, 53,
3711 & 67, 54, 68, 55, 69, 56, 70, 63, 68, 64, 69, 75, 54, 76, 55, 2, 8,
3712 & 18, 20, 16, 24, 14, 23, 101, 0, 88, 0, 2, 55, 9, 56, 1, 67, 8,
3713 & 68, 2, 64, 8, 75, 2, 59, 8, 72, 68, 55, 67, 54, 69, 56, 1, 9, 18,
3714 & 21, 15, 25, 13, 23, 100, 0, 96, 0, 2, 53, 9, 54, 1, 69, 8, 70, 1,
3715 & 76, 9, 63, 1, 73, 9, 58, 55, 70, 53, 68, 54, 69 /
3717 DATA spikp1/ 0.d0, 300.d0, 40.d0, 20.d0, 13.d0,8.5d0,8.d0, 9.5d0,
3718 & 12.d0,14.d0,15.5d0,20.d0,17.d0,13.d0,10.d0,9.d0,8.5d0,8.d0,7.8d0,
3719 & 7.3d0, 6.7d0, 9*0.d0,.23d0,.35d0,.7d0,.52d0,.4d0,.3d0,.2d0,.15d0,
3720 & .13d0, .11d0, .09d0, .07d0, 0.d0, .033d0,.8d0,1.35d0,1.35d0,.5d0,
3721 & 15*0.d0, 3*0.d0,.00d0,0.80d0,2.2d0,3.6d0,4.6d0,4.7d0,3.5d0,2.4d0,
3722 &1.8d0,1.4d0,.75d0,.47d0,.25d0,.13d0,.08d0,6*0.d0,0.d0,1.2d0,3.3d0,
3723 & 5.4d0,6.9d0,7.3d0,5.3d0,3.6d0,2.7d0,2.2d0,1.1d0,.73d0,.4d0,.22d0,
3724 & .12d0,9*0.d0,.0d0,0.d0,2.0d0,4.4d0,6.8d0,9.9d0,7.9d0,6.0d0,3.8d0,
3725 &2.5d0,2.d0,1.4d0,1.d0,.6d0,.35d0,10*0.d0,.25d0,.55d0,.75d0,1.25d0,
3726 & 1.9d0,2.d0,1.8d0,1.5d0,1.25d0,1.d0,.8d0,6*0.d0,4*0.d0,.4d0,.85d0,
3727 & 1.1d0, 1.85d0, 2.8d0, 3.d0,2.7d0,2.2d0,1.85d0,1.5d0,1.2d0,6*0.d0,
3728 & 6*0.d0, .5d0, 1.2d0, 1.7d0, 3.4d0, 5.2d0, 6.4d0, 6.1d0, 5.6d0,
3729 & 5.2d0, 6*0.d0, 2*0.d0, .0d0, 1.d0, 3.3d0, 5.2d0, 4.45d0, 3.6d0,
3730 & 2.75d0, 1.9d0, 1.65d0, 1.3d0, .95d0, .6d0, .45d0, 6*0.d0, 3*0.d0,
3731 & .0d0, .45d0, 1.4d0, 1.5d0, 1.1d0, .85d0, .5d0, .3d0, .2d0, .15d0,
3732 & 8*0.d0, 5*0.d0, .0d0, .0d0, .6d0, .8d0, .95d0, .8d0, .7d0, .6d0,
3733 & .5d0, .4d0, 6*0.d0, 5*0.d0, .0d0, .00d0, .85d0, 1.2d0, 1.4d0,
3734 & 1.2d0, 1.05d0, .9d0, .7d0, .55d0, 6*0.d0, 5*0.d0, .0d0, .00d0,
3735 & 1.d0, 1.5d0, 3.5d0, 4.15d0, 3.7d0, 2.7d0, 2.3d0, 1.75d0, 6*0.d0,
3736 & 10*0.d0, .5d0, 2.0d0, 3.3d0, 5.4d0, 7.d0 /
3738 DATA spikpu/ 0.d0, 25.d0, 13.d0, 11.d0, 10.5d0, 14.d0, 20.d0,
3739 & 20.d0, 16.d0, 14.d0, 19.d0, 28.d0, 17.5d0, 13.5d0, 12.d0, 10.5d0,
3740 & 10.d0, 10.d0, 9.5d0, 9.d0, 8.d0, 7.5d0, 7.d0, 6.5d0, 6.d0, 0.d0,
3741 & 48.d0, 19.d0, 15.d0, 11.5d0, 10.d0, 8.d0, 6.5d0, 5.5d0, 4.8d0,
3742 & 4.2d0, 7.5d0, 3.4d0, 2.5d0, 2.5d0, 2.1d0, 1.4d0, 1.d0, .8d0,
3743 & .6d0, .46d0, .3d0, .2d0, .15d0, .13d0, 11*0.d0, .95d0, .65d0,
3744 & .48d0, .35d0, .2d0, .18d0, .17d0, .16d0, .15d0, .1d0, .09d0,
3745 & .065d0, .05d0, .04d0, 12*0.d0, .2d0, .25d0, .25d0, .2d0, .1d0,
3746 & .08d0, .06d0, .045d0, .03d0, .02d0, .01d0, .005d0, .003d0,
3747 & 12*0.d0, .3d0, .24d0, .18d0, .15d0, .13d0, .12d0, .11d0, .1d0,
3748 & .09d0, .08d0, .05d0, .04d0, .03d0, 0.d0, 0.16d0, .7d0, 1.3d0,
3749 & 3.1d0, 4.5d0, 2.d0, 18*0.d0, 3*.0d0, 0.d0, 0.d0, 4.0d0, 11.d0,
3750 & 11.4d0, 10.3d0, 7.5d0, 6.8d0, 4.75d0, 2.5d0, 1.5d0, .9d0, .55d0,
3751 & .35d0, 13*0.d0, .1d0, .34d0, .5d0, .8d0, 1.1d0, 2.25d0, 3.3d0,
3752 & 2.3d0, 1.6d0, .95d0, .45d0, .28d0, .15d0, 10*0.d0, 2*0.d0, .17d0,
3753 & .64d0, 1.d0, 1.5d0, 2.1d0, 4.25d0, 6.2d0, 4.4d0, 3.d0, 1.8d0,
3754 & .9d0, .53d0, .28d0, 10*0.d0, 2*0.d0, .25d0, .82d0,
3755 & 1.3d0, 1.9d0, 2.8d0, 5.5d0 , 8.d0, 5.7d0, 3.9d0, 2.35d0, 1.15d0,
3756 & .69d0, .37d0, 10*0.d0, 7*0.d0, .0d0, .34d0, 1.5d0, 3.47d0,
3757 & 5.87d0, 6.23d0, 4.27d0, 2.6d0, 1.d0, .6d0, .3d0, .15d0, 6*0.d0/
3759 DATA spikpv/ 7*0.d0, .00d0, .16d0, .75d0, 1.73d0, 2.93d0, 3.12d0,
3760 & 2.13d0, 1.3d0, .5d0, .3d0, .15d0, .08d0, 6*0.d0, 10*0.d0, .2d0,
3761 & .6d0, .92d0, 2.4d0, 4.9d0, 6.25d0, 5.25d0, 3.5d0, 2.15d0, 1.4d0,
3762 & 1.d0, .7d0, 13*0.d0, .13d0, .4d0, .62d0, 1.6d0, 3.27d0, 4.17d0,
3763 & 3.5d0, 2.33d0, 1.43d0, .93d0, .66d0, .47d0, 13*0.d0, .07d0, .2d0,
3764 & .31d0, .8d0, 1.63d0, 2.08d0, 1.75d0, 1.17d0, .72d0, .47d0, .34d0,
3765 & .23d0, 17*0.d0, .33d0, 1.d0, 1.8d0, 2.67d0, 5.33d0, 6.d0, 5.53d0,
3766 & 5.d0, 17*0.d0, .17d0, .5d0, .9d0, 1.83d0, 2.67d0, 3.0d0, 2.77d0,
3767 & 2.5d0, 3*0.d0, 3*0.d0, 1.d0, 3.3d0, 2.8d0, 2.5d0, 2.3d0, 1.8d0,
3768 & 1.5d0, 1.1d0, .8d0, .7d0, .55d0, .3d0, 10*0.d0, 9*0.d0, .1d0,
3769 & .4d0, 1.d0, 1.4d0, 2.2d0, 2.5d0, 2.2d0, 1.65d0, 1.35d0, 1.1d0,
3770 & .8d0, .6d0, .4d0, 12*0.d0, .15d0, .6d0, 1.5d0, 2.1d0, 3.3d0,
3771 & 3.8d0, 3.3d0, 2.45d0, 2.05d0, 1.65d0, 1.2d0, .9d0, .6d0, 3*0.d0,
3772 & 9*0.d0, .10d0, .2d0, .5d0, .7d0, 1.3d0, 1.55d0, 1.9d0, 1.8d0,
3773 & 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0, 13*0.d0, .2d0, .5d0, .7d0,
3774 & 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0,
3775 & 17*0.d0, .2d0, .5d0, .85d0, 2.d0, 2.15d0, 2.05d0, 1.75d0, 1.d0,
3776 & 17*0.d0, .13d0, .33d0, .57d0, 1.33d0, 1.43d0, 1.36d0, 1.17d0,
3777 & .67d0, 17*0.d0, .07d0, .17d0, .28d0, .67d0, .72d0, .69d0, .58d0,
3778 & .33d0,17*0.d0,.4d0, .7d0, 1.d0, 1.6d0, 1.8d0, 2.3d0,1.9d0,1.7d0 /
3780 DATA spikpw/ 0.d0, 25.d0, 13.d0, 11.d0, 10.5d0, 14.d0, 2*20.d0,
3781 & 16.d0, 14.d0, 19.d0, 28.d0, 17.5d0, 13.5d0, 12.d0, 10.5d0,
3782 & 2*10.d0, 9.5d0, 9.d0, 8.d0, 7.5d0, 7.d0, 6.5d0, 6.d0, 0.d0,
3783 & 48.d0, 19.d0, 15.d0, 11.5d0, 10.d0, 8.d0, 6.5d0, 5.5d0, 4.8d0,
3784 & 4.2d0, 7.5d0, 3.4d0, 2*2.5d0, 2.1d0, 1.4d0, 1.d0, .8d0, .6d0,
3785 & .46d0, .3d0, .2d0, .15d0, .13d0, 11*0.d0, .95d0, .65d0, .48d0,
3786 & .35d0, .2d0, .18d0, .17d0, .16d0, .15d0, .1d0, .09d0, .065d0,
3787 & .05d0, .04d0, 12*0.d0, .2d0, 2*.25d0, .2d0, .1d0, .08d0, .06d0,
3788 & .045d0, .03d0, .02d0, .01d0, .005d0, .003d0, 12*0.d0, .3d0,
3789 & .24d0, .18d0, .15d0, .13d0, .12d0, .11d0, .1d0, .09d0, .08d0,
3790 & .05d0, .04d0, .03d0, 0.d0, 0.16d0, .7d0, 1.3d0, 3.1d0, 4.5d0,
3791 & 2.d0, 23*0.d0, 4.0d0, 11.d0, 11.4d0, 10.3d0, 7.5d0, 6.8d0,
3792 & 4.75d0, 2.5d0, 1.5d0, .9d0, .55d0, .35d0, 13*0.d0, .1d0, .34d0,
3793 & .5d0, .8d0, 1.1d0, 2.25d0, 3.3d0, 2.3d0, 1.6d0, .95d0, .45d0,
3794 & .28d0, .15d0, 12*0.d0, .17d0, .64d0, 1.d0, 1.5d0, 2.1d0, 4.25d0,
3795 & 6.2d0, 4.4d0, 3.d0, 1.8d0, .9d0, .53d0, .28d0, 12*0.d0, .25d0,
3796 & .82d0, 1.3d0, 1.9d0, 2.8d0, 5.5d0, 8.d0, 5.7d0, 3.9d0, 2.35d0,
3797 & 1.15d0, .69d0, .37d0, 18*0.d0, .34d0, 1.5d0, 3.47d0, 5.87d0,
3798 & 6.23d0, 4.27d0, 2.6d0, 1.d0, .6d0, .3d0, .15d0, 6*0.d0/
3800 DATA spikpx/ 8*0.d0, .16d0, .75d0, 1.73d0, 2.93d0, 3.12d0,
3801 & 2.13d0, 1.3d0, .5d0, .3d0, .15d0, .08d0, 16*0.d0, .2d0, .6d0,
3802 & .92d0, 2.4d0, 4.9d0, 6.25d0, 5.25d0, 3.5d0, 2.15d0, 1.4d0, 1.d0,
3803 & .7d0, 13*0.d0, .13d0, .4d0, .62d0, 1.6d0, 3.27d0, 4.17d0, 3.5d0,
3804 & 2.33d0, 1.43d0, .93d0, .66d0, .47d0, 13*0.d0, .07d0, .2d0, .31d0,
3805 & .8d0, 1.63d0, 2.08d0, 1.75d0, 1.17d0, .72d0, .47d0, .34d0, .23d0,
3806 & 17*0.d0, .33d0, 1.d0, 1.8d0, 2.67d0, 5.33d0, 6.d0, 5.53d0, 5.d0,
3807 & 17*0.d0, .17d0, .5d0, .9d0, 1.83d0, 2.67d0, 3.0d0, 2.77d0, 2.5d0,
3808 & 6*0.d0, 1.d0, 3.3d0, 2.8d0, 2.5d0, 2.3d0, 1.8d0, 1.5d0, 1.1d0,
3809 & .8d0, .7d0, .55d0, .3d0, 19*0.d0, .1d0, .4d0, 1.d0, 1.4d0, 2.2d0,
3810 & 2.5d0, 2.2d0, 1.65d0, 1.35d0, 1.1d0, .8d0, .6d0, .4d0, 12*0.d0,
3811 & .15d0, .6d0, 1.5d0, 2.1d0, 3.3d0, 3.8d0, 3.3d0, 2.45d0, 2.05d0,
3812 & 1.65d0, 1.2d0, .9d0, .6d0, 12*0.d0, .10d0, .2d0, .5d0, .7d0,
3813 & 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0,
3814 & 13*0.d0, .2d0, .5d0, .7d0, 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0,
3815 & 1.35d0, 1.15d0, .95d0, .7d0, 17*0.d0, .2d0, .5d0, .85d0, 2.d0,
3816 & 2.15d0, 2.05d0, 1.75d0, 1.d0, 17*0.d0, .13d0, .33d0, .57d0,
3817 & 1.33d0, 1.43d0, 1.36d0, 1.17d0, .67d0, 17*0.d0, .07d0, .17d0,
3818 & .28d0, .67d0, .72d0, .69d0, .58d0, .33d0, 17*0.d0, .4d0, .7d0,
3819 & 1.d0, 1.6d0, 1.8d0, 2.3d0, 1.9d0, 1.7d0 /
3821 DATA spikp4 / 0.d0, 300.d0, 40.d0, 20.d0, 13.d0, 8.5d0, 8.d0,
3822 & 9.5d0, 12.d0, 14.d0, 15.5d0, 20.d0, 17.d0, 13.d0, 10.d0, 9.d0,
3823 & 8.5d0, 8.d0, 7.8d0, 7.3d0, 6.7d0, 9*0.d0, .23d0, .35d0, .7d0,
3824 & .52d0, .4d0, .3d0, .2d0, .15d0, .13d0, .11d0, .09d0, .07d0, 0.d0,
3825 & .033d0, .8d0, 2*1.35d0, .5d0, 19*0.d0, 0.8d0, 2.2d0, 3.6d0,
3826 & 4.6d0, 4.7d0, 3.5d0, 2.4d0, 1.8d0, 1.4d0, .75d0, .47d0, .25d0,
3827 & .13d0, .08d0, 7*0.d0, 1.2d0, 3.3d0, 5.4d0, 6.9d0, 7.3d0, 5.3d0,
3828 & 3.6d0, 2.7d0, 2.2d0, 1.1d0, .73d0, .4d0, .22d0, .12d0, 11*0.d0,
3829 & 2.0d0, 4.4d0, 6.8d0, 9.9d0, 7.9d0, 6.0d0, 3.8d0, 2.5d0, 2.d0,
3830 & 1.4d0, 1.d0, .6d0, .35d0, 10*0.d0, .25d0, .55d0, .75d0, 1.25d0,
3831 & 1.9d0, 2.d0, 1.8d0, 1.5d0, 1.25d0, 1.d0, .8d0, 10*0.d0, .4d0,
3832 & .85d0, 1.1d0, 1.85d0, 2.8d0, 3.d0, 2.7d0, 2.2d0, 1.85d0, 1.5d0,
3833 & 1.2d0, 12*0.d0, .5d0, 1.2d0, 1.7d0, 3.4d0, 5.2d0, 6.4d0, 6.1d0,
3834 & 5.6d0, 5.2d0, 9*0.d0, 1.d0, 3.3d0, 5.2d0, 4.45d0, 3.6d0, 2.75d0,
3835 & 1.9d0, 1.65d0, 1.3d0, .95d0, .6d0, .45d0, 10*0.d0, .45d0, 1.4d0,
3836 & 1.5d0, 1.1d0, .85d0, .5d0, .3d0, .2d0, .15d0, 15*0.d0, .6d0,
3837 & .8d0, .95d0, .8d0, .7d0, .6d0, .5d0, .4d0, 13*0.d0, .85d0, 1.2d0,
3838 & 1.4d0, 1.2d0, 1.05d0, .9d0, .7d0, .55d0, 13*0.d0, 1.d0, 1.5d0,
3839 & 3.5d0, 4.15d0, 3.7d0, 2.7d0, 2.3d0, 1.75d0, 16*0.d0, .5d0, 2.0d0,
3840 & 3.3d0, 5.4d0, 7.d0 /
3842 DATA spikp5/ 0.d0, 20.d0, 14.d0, 12.d0, 11.5d0, 10.d0, 8.d0,
3843 & 7.d0, 6.d0, 5.5d0, 5.3d0, 5.d0, 4.5d0, 4.4d0, 3.8d0, 3.d0, 2.8d0,
3844 & 0.d0, .5d0, 1.15d0, 2.d0, 1.3d0, .8d0, .45d0, 13*0.d0, 0.9d0,
3845 & 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0, 1.5d0, 1.2d0, .9d0, .6d0,
3846 & .45d0, .21d0, .2d0, 3*0.d0, .9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0,
3847 & 2.d0, 1.7d0, 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
3848 & 4*0.d0, 1.d0, 2.1d0, 2.6d0, 2.3d0, 2.1d0, 1.8d0, 1.7d0, 1.4d0,
3849 & 1.2d0, 1.05d0, .9d0, .66d0, .5d0, 7*0.d0, .3d0, 2*1.d0, .9d0,
3850 & .7d0, .4d0, .3d0, .2d0, 11*0.d0, .1d0, 1.d0, 2.2d0, 3.5d0, 4.2d0,
3851 & 4.55d0, 4.85d0, 4.9d0, 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0,
3852 & 1.71d0, 1.6d0, 6*0.d0, 1.4d0, 3.8d0, 5.d0, 4.7d0, 4.4d0, 4.d0,
3853 & 3.5d0, 2.85d0, 2.35d0, 2.01d0, 1.8d0, 12*0.d0, .1d0, .8d0,2.05d0,
3854 & 3.31d0, 3.5d0, 12*0.d0, .034d0, .2d0, .75d0, 1.04d0, 1.24d0 /
3856 DATA spikp6/ 0.d0, 6.d0, 11.d0, 13.d0, 6.d0, 5.d0, 3.d0, 2.2d0,
3857 & 1.5d0, 1.2d0, 1.d0, .7d0, .6d0, .5d0, .45d0, .35d0, .3d0, 0.d0,
3858 & 6.d0, 11.d0, 13.d0, 6.d0, 5.d0, 3.d0, 2.2d0, 1.5d0, 1.2d0, 1.d0,
3859 & .7d0, .6d0, .5d0, .45d0, .35d0, .3d0, 0.d0, .5d0, 1.3d0, 2.8d0,
3860 & 2.3d0, 1.6d0, .9d0, 13*0.d0, 0.9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0,
3861 & 2.d0, 1.7d0, 1.5d0,1.2d0,.9d0,.6d0,.45d0,.21d0,.2d0,3*0.d0,0.9d0,
3862 & 2.5d0, 3.d0, 2.5d0, 2.3d0,2.d0,1.7d0,1.5d0,1.2d0,.9d0,.6d0,.45d0,
3863 & .21d0, .2d0,4*0.d0,1.d0,2.1d0,2.6d0,2.3d0,2.d0,1.8d0,1.7d0,1.4d0,
3864 & 1.2d0,1.15d0,.9d0,.66d0,.5d0,4*0.d0,1.d0,2.1d0,2.6d0,2.3d0,2.1d0,
3865 & 1.8d0,1.7d0,1.4d0,1.2d0, 1.15d0, .9d0, .66d0, .5d0, 7*0.d0, .3d0,
3866 & 2*1.d0, .9d0, .7d0, .4d0, .35d0, .2d0, 9*0.d0, .3d0, 2*1.d0,.9d0,
3867 & .7d0, .4d0, .35d0, .2d0, 11*0.d0, .1d0, 1.d0, 2.4d0,3.5d0,4.25d0,
3868 & 4.55d0, 4.85d0, 4.9d0, 9*0.d0, .1d0, 1.d0, 2.4d0, 3.5d0, 4.25d0,
3869 & 4.55d0, 4.85d0, 4.9d0, 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0,
3870 & 1.71d0, 1.6d0, 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0, 1.71d0,
3871 & 1.6d0, 6*0.d0, 1.4d0, 3.8d0, 5.d0, 4.7d0,4.4d0,4.d0,3.5d0,2.85d0,
3872 & 2.35d0, 2.01d0, 1.8d0, 6*0.d0, 1.4d0,3.8d0,5.d0,4.7d0,4.4d0,4.d0,
3873 & 3.5d0,2.85d0,2.35d0,2.01d0,1.8d0,12*0.d0,.1d0,.8d0,2.05d0,3.31d0,
3874 & 3.5d0, 12*0.d0, .034d0,.2d0,.75d0,1.04d0,1.24d0 /
3876 DATA skmpel/ 0.d0, 35.d0, 22.d0, 25.d0, 17.d0, 9.d0, 9.5d0, 8.d0,
3877 & 7.d0, 6.5d0, 6.1d0, 5.d0, 4.8d0, 4.6d0, 4.45d0, 4.3d0, 4.2d0,
3878 & 0.d0, 8.d0, 3.5d0, 8.d0, 3.d0, 1.9d0, 1.7d0, 1.d0, .9d0, .8d0,
3879 & .75d0, .5d0, .42d0, .38d0, .34d0, .25d0, .2d0,
3880 & 0.d0, 3.d0, 3.2d0, 3.5d0, 1.5d0, 1.4d0, 1.1d0, .6d0, .5d0,
3881 & .35d0, .28d0, .25d0, .18d0, .12d0, .1d0, .08d0, .04d0,
3882 & 0.d0, 8.5d0, 2.4d0, 1.7d0, 1.3d0, 1.3d0, 1.1d0, .5d0,
3883 & .4d0, .4d0, .35d0, .3d0, .28d0, .2d0, .16d0, .13d0, .11d0,
3884 & 0.d0, 7.d0, 4.8d0, 1.4d0, 1.9d0, .9d0, .4d0, .2d0, .13d0,
3885 & .1d0, .08d0, .06d0, .04d0, .02d0, .015d0, .01d0, .01d0,
3886 & 0.d0, 5.5d0, 1.d0, .8d0, .75d0, .32d0, .2d0, .1d0, .09d0,
3887 & .08d0, .065d0, .05d0, .04d0, .022d0, .017d0, 2*.01d0/
3888 DATA spikp7 / 0.d0, .56d0, 1.46d0, 3.16d0, 2.01d0, 1.28d0, .74d0,
3889 & 14*0.d0, 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0,
3890 & 1.91d0, 1.57d0, 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0, 4*0.d0,
3891 & 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0, 1.91d0, 1.57d0,
3892 & 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0, 3*0.d0, 1.0d0, 3.03d0,
3893 & 3.36d0, 2.8d0, 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0, 1.01d0,
3894 & .67d0, .5d0, .24d0, .23d0, 3*0.d0, 1.0d0, 3.03d0, 3.36d0, 2.8d0,
3895 & 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0, 1.01d0, .67d0, .5d0,
3896 & .24d0, .23d0, 7*0.d0, .34d0, 1.12d0, 1.12d0, 1.01d0, .78d0,
3897 & .45d0, .39d0, .22d0, .07d0, 0.d0, 7*0.d0, .34d0, 1.12d0, 1.12d0,
3898 & 1.01d0, .78d0, .45d0, .39d0, .22d0, .07d0, 0.d0, 6*0.d0, 1.71d0,
3899 & 4.26d0, 5.6d0, 5.57d0, 4.93d0, 4.48d0, 3.92d0, 3.19d0, 2.63d0,
3900 & 2.25d0, 2.d0, 6*0.d0, 1.71d0, 4.26d0, 5.6d0, 5.57d0, 4.93d0,
3901 & 4.48d0, 3.92d0, 3.19d0, 2.63d0, 2.25d0, 2.d0, 10*0.d0, .22d0,
3902 & .8d0, .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0, 10*0.d0, .22d0, .8d0,
3903 & .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0, 13*0.d0, .1d0, .3d0, .7d0,1.d0,
3904 & 13*0.d0, .1d0, .3d0, .7d0, 1.d0, 9*0.d0, .11d0, 1.72d0, 2.69d0,
3905 & 3.92d0, 4.76d0, 5.10d0, 5.44d0, 5.3d0, 9*0.d0, .11d0, 1.72d0,
3906 & 2.69d0, 3.92d0, 4.76d0, 5.1d0, 5.44d0, 5.3d0, 5*0.d0,9.2d0,4.7d0,
3907 & 1.9d0, 10*0.d0, 2.5d0, 15.d0, 21.5d0, 15.3d0, 3.d0, 1.5d0,
3910 DATA skmnel/0.d0, 4.d0, 9.5d0, 20.d0, 13.d0, 9.5d0, 6.d0, 4.4d0,
3911 & 3.d0, 2.4d0, 2.d0, 1.4d0, 1.2d0, 1.d0, .9d0, .7d0, .6d0,
3912 & 0.d0, 4.5d0, 6.d0, 5.d0, 2.5d0, 2.d0, 1.7d0, 2.1d0,
3913 & 1.9d0, .9d0, .5d0, .3d0, .24d0, .2d0, .18d0, .1d0, .09d0,
3914 & 0.d0, 1.8d0, 2.d0, 1.1d0, .9d0, .5d0, .5d0, .4d0, .4d0,
3915 & .2d0, .1d0, .06d0, .05d0, .04d0, .03d0, .02d0, .02d0,
3916 & 0.d0, 1.5d0, 2.d0, .9d0, 1.1d0, .4d0, .6d0, .7d0, .65d0,
3917 & .3d0, .17d0, .1d0, .08d0, .07d0, .06d0, .04d0, .03d0/
3918 DATA spikp8/0.d0, .56d0, 1.29d0, 2.26d0, 1.01d0, .64d0, .37d0,
3919 & 14*0.d0, 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0,
3920 & 1.91d0, 1.57d0, 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0,
3921 & 3*0.d0, 1.d0, 3.03d0, 3.36d0, 2.8d0, 2.58d0, 2.24d0,
3922 & 1.91d0, 1.68d0, 1.35d0, 1.01d0, .67d0, .5d0, .24d0, .23d0,
3923 & 3*0.d0, 1.d0, 3.03d0, 3.36d0, 2.8d0, 2.58d0, 2.24d0,
3924 & 1.91d0, 1.68d0, 1.35d0, 1.01d0, .67d0, .5d0, .24d0, .23d0,
3925 & 7*0.d0, .34d0, 1.12d0, 1.12d0, 1.01d0, .78d0, .45d0,
3926 & .39d0, .22d0, .07d0, 0.d0,
3927 & 6*0.d0, 1.71d0, 4.26d0, 5.6d0, 5.57d0, 4.93d0,
3928 & 4.48d0, 3.92d0, 3.19d0, 2.63d0, 2.25d0, 2.d0,
3929 & 10*0.d0, .22d0, .8d0, .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0,
3930 & 13*0.d0, .1d0, .3d0, .7d0, 1.d0,
3931 & 13*0.d0, .1d0, .3d0, .7d0, 1.d0,
3932 & 9*0.d0, .11d0, 1.72d0, 2.69d0, 3.92d0, 4.76d0,
3933 & 5.10d0, 5.44d0, 5.3d0,
3934 & 4*0.d0, 0.00d0, 9.2d0, 4.7d0, 1.9d0, 9*0.d0/
3936 DATA spikp9/ 0.d0, 24.d0, 25.d0, 27.d0, 23.d0, 21.d0, 20.d0,
3937 & 19.d0, 17.d0, 15.5d0, 14.d0, 13.5d0, 13.d0,
3938 & 0.d0, 3.6d0, 1.7d0, 10*0.d0,
3939 & .0d0, 0.d0, 8.7d0, 17.7d0, 18.8d0, 15.9d0,
3940 & 11.7d0, 8.d0, 6.d0, 5.3d0, 4.5d0, 3.9d0, 3.5d0,
3941 & .0d0, .0d0, 2.8d0, 5.8d0, 6.2d0, 5.1d0, 3.8d0,
3942 & 2.7d0, 2.1d0, 1.8d0, 1.5d0, 1.3d0, 1.1d0,
3943 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0,
3944 & 16.9d0, 16.5d0, 11.d0, 5.5d0, 3.5d0,
3945 & 10*0.d0, 4.3d0, 7.6d0, 9.d0,
3946 & 10*0.d0, 1.7d0, 2.6d0, 3.d0,
3947 & 6*0.d0, .3d0, .6d0, 1.d0, 1.6d0, 1.3d0, .8d0, .6d0,
3948 & 6*0.d0, .7d0, 1.2d0, 1.8d0, 2.5d0, 1.8d0, 1.3d0,
3949 & 1.2d0, 10*0.d0, .6d0, 1.4d0, 1.7d0,
3950 & 10*0.d0, 1.9d0, 4.1d0, 5.2d0/
3952 DATA spikp0/ 0.d0, 24.d0, 25.d0, 27.d0, 23.d0, 21.d0, 20.d0,
3953 & 19.d0, 17.d0, 15.5d0, 14.d0, 13.5d0, 13.d0,
3954 & 0.d0, 1.8d0, .2d0, 12*0.d0,
3955 & 3.2d0, 6.05d0, 9.9d0, 5.1d0,
3956 & 3.8d0, 2.7d0, 1.9d0, 1.5d0, 1.4d0, 1.3d0, 1.1d0,
3957 & 2*.0d0, 3.2d0, 6.05d0, 9.9d0, 5.1d0,
3958 & 3.8d0, 2.7d0, 1.9d0, 1.5d0, 1.4d0, 1.3d0, 1.1d0,
3959 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0,
3960 & 16.4d0, 15.2d0, 11.d0, 5.4d0, 3.5d0,
3961 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0,
3962 & 16.4d0, 15.2d0, 11.d0, 5.4d0, 3.5d0,
3963 & 10*0.d0, .7d0, 5.1d0, 8.d0,
3964 & 10*0.d0, .7d0, 5.1d0, 8.d0,
3965 & 10*.0d0, .3d0, 2.8d0, 4.7d0,
3966 & 10*.0d0, .3d0, 2.8d0, 4.7d0,
3967 & 7*0.d0, 1.2d0, 2.5d0, 3.5d0, 6.d0, 5.3d0, 2.9d0,
3968 & 7*0.d0, 1.7d0, 3.6d0, 5.4d0, 9.d0, 7.6d0, 4.2d0,
3969 & 5*0.d0, 7.7d0, 6.1d0, 2.9d0, 5*0.d0/
3972 DATA spkp v/ 0.d0, 24.d0, 25.d0, 27.d0, 23.d0, 21.d0, 20.d0,
3973 & 19.d0, 17.d0, 15.5d0, 14.d0, 13.5d0, 13.d0,
3974 & 0.d0, 3.6d0, 1.7d0, 12*0.d0,
3975 & 8.7d0, 17.7d0, 18.8d0, 15.9d0,
3976 & 11.7d0, 8.d0, 6.d0, 5.3d0, 4.5d0, 3.9d0, 3.5d0,
3977 & .0d0, .0d0, 2.8d0, 5.8d0, 6.2d0, 5.1d0, 3.8d0,
3978 & 2.7d0, 2.1d0, 1.8d0, 1.5d0, 1.3d0, 1.1d0,
3979 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0, 16.9d0, 16.5d0,
3980 & 11.d0, 5.5d0, 3.5d0,
3981 & 10*0.d0, 4.3d0, 7.6d0, 9.d0,
3982 & 10*0.d0, 1.7d0, 2.6d0, 3.d0,
3983 & 6*0.d0, .3d0, .6d0, 1.d0, 1.6d0, 1.3d0, .8d0, .6d0,
3984 & 6*0.d0, .7d0, 1.2d0, 1.8d0, 2.5d0, 1.8d0, 1.3d0,
3985 & 1.2d0, 10*0.d0, .6d0, 1.4d0, 1.7d0,
3986 & 10*0.d0, 1.9d0, 4.1d0, 5.2d0/
3988 DATA sappel/ 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0,
3989 & 50.d0, 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
3990 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0,
3991 & .05d0, .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0, .85d0,
3992 & 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0, .14d0,
3993 & .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
3994 & 0.d0, 55.d0, 50.d0, 25.d0, 15.d0, 15.d0, 14.d0, 12.d0,
3995 & 10.d0, 7.d0, 6.d0, 4.d0, 3.3d0, 2.8d0, 2.4d0, 2.d0, 1.8d0,
3996 & 1.55d0, 1.3d0, .95d0, .75d0,
3997 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0,
3998 & .25d0, .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
3999 & .01d0, .008d0, .006d0, .005d0/
4000 DATA spikpe/0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4001 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4002 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0,
4003 & 9*0.d0, 2.d0, 2.5d0, .2d0, 19*0.d0, .3d0, 1.4d0, 2.2d0, 1.2d0,
4004 & 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0, .3d0, 10*0.d0, .3d0, 1.4d0,
4005 & 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0, .3d0, 10*0.d0,
4006 & .3d0, 1.4d0, 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0,
4007 & .3d0, 10*0.d0, .3d0, 1.4d0, 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0,
4008 & .6d0, .5d0, .4d0, .3d0, 9*0.d0, .6d0, 2.5d0, 5.d0, 5.2d0, 5.1d0,
4009 & 5.4d0, 5.8d0, 2.8d0, 2.1d0, 1.8d0, 1.6d0, 1.2d0, 13*0.d0, 1.3d0,
4010 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4011 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4012 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4013 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 14*0.d0, .2d0,
4014 & .5d0, 1.1d0, 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0,
4015 & 1.1d0, 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0, 1.1d0,
4016 & 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0, 1.1d0, 1.6d0,
4017 & 1.4d0, 1.1d0, .9d0, 17*0.d0, .3d0, 1.6d0, 2.6d0, 3.6d0, 17*0.d0,
4018 & .3d0, 1.6d0, 2.6d0, 3.6d0, 17*0.d0, .3d0, 1.6d0, 2.6d0,
4019 & 3.6d0, 17*0.d0, .3d0, 1.6d0, 2.6d0, 3.6d0 /
4022 & 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0,
4023 & 50.d0, 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
4024 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0,
4025 & .05d0, .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0,
4026 & .85d0, 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0,
4027 & .14d0, .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
4028 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0,
4029 & .25d0, .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
4030 & .01d0, .008d0, .006d0, .005d0 /
4031 DATA spikpz/ 0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4032 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4033 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0, 9*0.d0, 2.4d0, .2d0,
4034 & 20*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4035 & .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0,
4036 & 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0,
4037 & 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0,
4038 & 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4039 & .7d0, .5d0, .3d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4040 & 7.6d0, 6.d0, 5.d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4041 & 7.6d0, 6.d0, 5.d0, 18*0.d0, 1.d0, 4.9d0, 8.5d0, 18*0.d0, 1.d0,
4042 & 4.9d0, 8.5d0, 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0,
4043 & 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0, 15*0.d0, 1.9d0,
4044 & 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0 /
4050 & 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0, 50.d0,
4051 & 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
4052 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0, .05d0,
4053 & .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0, .85d0,
4054 & 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0, .14d0,
4055 & .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
4056 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0, .25d0,
4057 & .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
4058 & .01d0, .008d0, .006d0, .005d0 /
4059 DATA spikpf/ 0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4060 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4061 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0, 9*0.d0, 2.4d0, .2d0,
4062 & 20*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4063 & .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0,
4064 & 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0,
4065 & 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0,
4066 & 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4067 & .7d0, .5d0, .3d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4068 & 7.6d0, 6.d0, 5.d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4069 & 7.6d0, 6.d0, 5.d0, 18*0.d0, 1.d0, 4.9d0, 8.5d0, 18*0.d0, 1.d0,
4070 & 4.9d0, 8.5d0, 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0,
4071 & 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0, 15*0.d0, 1.9d0,
4072 & 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0 /
4074 DATA spkp15/0.d0, 20.d0, 14.d0, 12.d0, 11.5d0, 10.d0, 8.d0, 7.d0,
4075 & 6.d0, 5.5d0, 5.3d0, 5.d0, 4.5d0, 4.4d0, 3.8d0, 3.d0, 2.8d0,
4076 & 0.d0, .5d0, 1.15d0, 2.d0, 1.3d0, .8d0, .45d0, 10*0.d0,
4077 & 3*0.d0, 0.9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0,
4078 & 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
4079 & 3*0.d0, 0.9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0,
4080 & 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
4081 & 4*0.d0, 1.d0, 2.1d0, 2.6d0, 2.3d0, 2.1d0, 1.8d0, 1.7d0,
4082 & 1.4d0, 1.2d0, 1.05d0, .9d0, .66d0, .5d0,
4083 & 7*0.d0, .3d0, 1.d0, 1.d0, .9d0, .7d0, .4d0, .30d0, .2d0,
4084 & 11*0.d0, .1d0, 1.d0, 2.2d0, 3.5d0, 4.20d0, 4.55d0,
4086 & 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0, 1.71d0, 1.6d0,
4087 & 6*0.d0, 1.4d0, 3.8d0, 5.d0, 4.7d0, 4.4d0, 4.d0, 3.5d0,
4088 & 2.85d0, 2.35d0, 2.01d0, 1.8d0,
4089 & 12*0.d0, .1d0, .8d0, 2.05d0, 3.31d0, 3.5d0,
4090 & 12*0.d0, .034d0, .20d0, .75d0, 1.04d0, 1.24d0 /
4092 DATA spkp16/ 0.d0, 4.d0, 9.5d0, 20.d0, 13.d0, 9.5d0, 6.d0, 4.4d0,
4093 & 3.d0, 2.4d0, 2.d0, 1.4d0, 1.2d0, 1.d0, .9d0, .7d0, .6d0, 0.d0,
4094 & 4.5d0, 6.d0, 5.d0, 2.5d0, 2.d0, 1.7d0, 2.1d0, 1.9d0, .9d0, .5d0,
4095 & .3d0, .24d0, .2d0, .18d0, .1d0, .09d0, 0.d0, 1.8d0, 2.d0, 1.1d0,
4096 & .9d0, .5d0, .5d0, .4d0, .4d0, .2d0, .1d0, .06d0, .05d0, .04d0,
4097 & .03d0, .02d0, .02d0, 0.d0, 1.5d0, 2.d0, .9d0, 1.1d0, .4d0, .6d0,
4098 & .7d0, .65d0, .3d0, .17d0, .1d0, .08d0, .07d0, .06d0, .04d0,
4099 & .03d0, 0.d0, .56d0, 1.29d0, 2.26d0, 1.01d0, .64d0, .37d0,
4100 & 14*0.d0, 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0, 1.91d0,
4101 & 1.57d0, 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0, 3*0.d0, 1.0d0,
4102 & 3.03d0, 3.36d0, 2.8d0, 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0,
4103 & 1.01d0, .67d0, .5d0, .24d0, .23d0, 3*0.d0, 1.0d0, 3.03d0, 3.36d0,
4104 & 2.8d0, 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0, 1.01d0, .67d0,
4105 & .5d0, .24d0, .23d0, 7*0.d0, .34d0, 1.12d0, 1.12d0, 1.01d0, .78d0,
4106 & .45d0, .39d0, .22d0, .07d0, 7*0.d0, 1.71d0, 4.26d0, 5.6d0,5.57d0,
4107 & 4.93d0, 4.48d0, 3.92d0, 3.19d0, 2.63d0, 2.25d0, 2.d0, 10*0.d0,
4108 & .22d0, .8d0, .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0, 13*0.d0, .1d0,
4109 & .3d0, .7d0, 1.d0, 13*0.d0, .1d0, .3d0, .7d0, 1.d0, 9*0.d0, .11d0,
4110 & 1.72d0, 2.69d0, 3.92d0, 4.76d0, 5.10d0, 5.44d0, 5.3d0, 5*0.d0,
4111 & 9.2d0, 4.7d0, 1.9d0, 9*0.d0, .0d0,2.5d0,15.d0,
4112 & 21.5d0, 15.3d0, 3.d0, 1.5d0, 10*0.d0 /
4113 DATA nureln/9, 12, 5*0, 10, 14, 3*0, 1, 3, 5, 7, 6*0, 2, 6, 16,
4114 & 5*0, 10, 13, 5*0, 11, 12, 3*0, 2, 4, 6, 8, 6*0, 3, 15, 7, 5*0 /
4120 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4126 *umopi( 92),umokc( 68),umop( 39),umon( 63),umok0( 34),
4127 *plapi( 92),plakc( 68),plap( 39),plan( 63),plak0( 34),
4129 *spikp1( 315),spikp u(278),spikp v(372),
4130 *spikp
w(278),spikp
x(372),spikp 4(315),
4131 *spikp 5(187),spikp 6(306),
4132 *
s kmpel(102),spikp 7(289),
s kmnel( 68),spikp 8(187),
4133 *spikp 9(143),spikp 0(169),spkp v(143),
4134 *
s appel(105),spikp
e(399),
s apnel( 84),spikp
z(273),
4135 *
s anpel( 84),spikp
f(273),
4136 *spkp15(187),spkp16(255),
4137 *nrkpi( 164),nrkkc( 134),nrkp( 70),nrkn( 116),nrkk0( 52),
4143 DATA plapi / 0.d0, .3d0, .5d0, .6d0, .7d0, .8d0, .9d0, .95d0,1.d0,
4144 & 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0, 2.d0, 2.3d0, 2.5d0, 2.8d0,
4145 & 3.d0, 3.5d0, 4.d0, 0.d0, .285d0, .4d0, .45d0, .5d0, .6d0, .7d0,
4146 & .75d0, .8d0, .85d0, .9d0, 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0,
4147 & 1.8d0, 2.d0, 2.3d0, 2.5d0, 2.8d0, 3.d0, 3.5d0, 4.d0, 4.5d0, 0.d0,
4148 & .285d0, .4d0, .45d0, .5d0, .6d0, .7d0, .75d0, .8d0, .85d0, .9d0,
4149 & 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0, 2.d0, 2.3d0, 2.5d0,
4150 & 2.8d0, 3.d0, 3.5d0, 4.d0, 4.5d0, 0.d0, .3d0, .5d0, .6d0, .7d0,
4151 & .8d0, .9d0, .95d0, 1.d0, 1.15d0, 1.3d0, 1.5d0, 1.6d0, 1.8d0,
4152 & 2.d0, 2.3d0, 2.5d0, 2.8d0, 3.d0, 3.5d0, 4.d0 /
4154 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4155 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4156 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
4157 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4158 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4159 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
4160 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4161 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4162 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
4163 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4164 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4165 & 3.51d0, 3.84d0, 4.16d0, 4.49d0/
4167 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4168 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4169 & 3.51d0, 3.84d0, 4.16d0, 4.49d0,
4170 & 0.d0, .58d0, .8d0, 1.01d0, 1.23d0, 1.45d0, 1.68d0, 1.94d0,
4171 & 2.18d0, 2.42d0, 2.68d0, 2.96d0, 3.24d0,
4172 & 3.51d0, 3.84d0, 4.16d0, 4.49d0/
4175 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4176 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
4177 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4178 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
4179 & 0.d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4180 & 3.43d0, 3.75d0, 4.07d0, 4.43d0 /
4183 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
4184 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4185 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
4186 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
4187 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4188 & 3.43d0, 3.75d0, 4.07d0, 4.43d0,
4189 & 0.d0, 1.
d-3, .1d0, .2d0, .3d0, .4d0, .5d0, .6d0,
4190 & .74d0, 1.06d0, 1.34d0, 1.63d0, 1.92d0, 2.2d0, 2.5d0,2.8d0,3.1d0,
4191 & 3.43d0, 3.75d0, 4.07d0, 4.43d0 /
4192 DATA siin / 296*0.d0 /
4193 DATA umopi/ 1.08d0,1.233d0,1.302d0,1.369d0,1.496d0,
4194 & 1.557d0,1.615d0,1.6435d0,
4195 & 1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,2.071d0,2.159d0,
4196 & 2.286d0,2.366d0,2.482d0,2.56d0,
4198 & 1.08d0,1.222d0,1.302d0,1.3365d0,1.369d0,1.434d0,
4199 & 1.496d0,1.527d0,1.557d0,
4200 & 1.586d0,1.615d0,1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,
4201 & 2.071d0,2.159d0,2.286d0,2.366d0,
4202 & 2.482d0,2.560d0,2.735d0,2.90d0,3.06d0,
4203 & 1.08d0,1.222d0,1.302d0,1.3365d0,1.369d0,1.434d0,
4204 & 1.496d0,1.527d0,1.557d0,
4205 & 1.586d0,1.615d0,1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,
4206 & 2.071d0,2.159d0,2.286d0,2.366d0,
4207 & 2.482d0,2.560d0,2.735d0,2.90d0,3.06d0,
4208 & 1.08d0,1.233d0,1.302d0,1.369d0,1.496d0,
4209 & 1.557d0,1.615d0,1.6435d0,
4210 & 1.672d0,1.753d0,1.831d0,1.930d0,1.978d0,2.071d0,2.159d0,
4211 & 2.286d0,2.366d0,2.482d0,2.56d0,
4214 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4215 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4217 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4218 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4220 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4221 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4223 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4224 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4227 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4228 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4230 & 1.598d0,1.7d0,1.8d0,1.9d0,2.0d0,2.1d0,2.2d0,2.3d0,2.4d0,2.5d0,
4231 & 2.6d0,2.7d0,2.8d0,2.9d0,3.0d0,
4235 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4237 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4239 & 1.88d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4243 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
4244 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4246 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
4247 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4249 & 1.877d0,1.87701d0,1.879d0,1.887d0,1.9d0,1.917d0,1.938d0,1.962d0,
4250 & 2.d0,2.102d0,2.2d0,2.3d0,2.4d0,2.5d0,2.6d0,2.7d0,2.8d0,2.9d0,
4258 *13,54,23,53,13,63,13,58,23,57,13,65,1,32,53,31,54,32,53,33,53,35,
4260 *13,8,23,1,17,15,21,24,22,15,82,0,
4261 *61,0,13,55,23,54,14,53,13,64,
4262 *23,63,13,59,23,58,14,57,13,66,23,65,1,31,8,32,1,33,1,35,54,31,55,
4263 *32,54,33,53,34,54,35,
4264 *14,1,23,8,17,24,20,15,22,24,83,0,
4265 *62,0,14,54,23,55,13,56,14,63,
4266 *23,64,14,58,23,59,13,60,14,65,23,66,8,31,1,34,8,33,8,35,55,31,54,
4267 *34,55,33,56,32,55,35,
4269 *14,55,23,56,14,64,14,59,23,60,14,66,8,34,56,31,55,34,56,33,56,35,
4274 *24,53,15,54,1,36,1,40,1,44,36,63,15,63,45,53,44,54,
4276 *24,54,15,55,8,36,1,37,8,40,1,41,8,44,1,45,36,64,37,63,15,64,24,63,
4278 *16,1,25,8, 17,23,21,14, 20,13,22,23, 90,0,
4279 *38,1,39,8,16,54,25,55,1,42,8,43,16,63,25,64,39,64,38,63,46,54,
4280 *47,55,8,47,1,46,52,0,51,0,
4281 *16,8,17,14,20,23,22,14,92,0,
4282 *8,38,16,55,25,56,8,42,16,64,38,64,46,55,47,56,8,46,94,0
4289 *106,0,15,56,24,55,37,8,41,8,45,8,37,64,24,64,44,56,45,55,
4290 *25,1,17,13, 22,13,21,23,
4291 *107,0,39,1,25,54,16,53,43,1,25,63,39,63,47,54,46,53,47,1,103,0/
4293 DATA nrk
p/1,1,85,0,
4294 f8,53,1,54,1,63,8,57,1,58,54,54,53,55,63,54,64,53,
4296 f8,54,1,55,8,63,1,64,8,58,1,59,64,54,63,55,54,55,53,56,77,0,
4298 *95,0,8,55,1,56,8,64,8,59,1,60,55,55,54,56,64,55,63,56/
4300 DATA nrk
n/ 1,2,17,18,15,16,8,9,13,14,99,0,87,0, 1,68,8,69,2,54,9,
4301 +55,102,0, 2,63,9,64,1,75,8,76,53,67,54,68, 55,69,56,70,63,68,64,
4302 +69,75,54,76,55, 2,8,18,20, 16,24,14,23, 101,0,88,0, 2,55,9,56,1,
4303 +67,8,68,2,64,8,75,2,59,8,72,68,55,67,54,69,56, 1,9,18,21,15,25,13,
4304 +23,100,0, 96,0,2,53,9,54,1,69,8,70,1,76,9,63,1,73,9,58,55,70,53,
4308 DATA spikp1/ 0.d0, 300.d0, 40.d0, 20.d0, 13.d0,8.5d0,8.d0, 9.5d0,
4309 & 12.d0,14.d0,15.5d0,20.d0,17.d0,13.d0,10.d0,9.d0,8.5d0,8.d0,7.8d0,
4310 & 7.3d0, 6.7d0, 9*0.d0,.23d0,.35d0,.7d0,.52d0,.4d0,.3d0,.2d0,.15d0,
4311 & .13d0, .11d0, .09d0, .07d0, 0.d0, .033d0,.8d0,1.35d0,1.35d0,.5d0,
4312 & 15*0.d0, 3*0.d0,.00d0,0.80d0,2.2d0,3.6d0,4.6d0,4.7d0,3.5d0,2.4d0,
4313 &1.8d0,1.4d0,.75d0,.47d0,.25d0,.13d0,.08d0,6*0.d0,0.d0,1.2d0,3.3d0,
4314 & 5.4d0,6.9d0,7.3d0,5.3d0,3.6d0,2.7d0,2.2d0,1.1d0,.73d0,.4d0,.22d0,
4315 & .12d0,9*0.d0,.0d0,0.d0,2.0d0,4.4d0,6.8d0,9.9d0,7.9d0,6.0d0,3.8d0,
4316 &2.5d0,2.d0,1.4d0,1.d0,.6d0,.35d0,10*0.d0,.25d0,.55d0,.75d0,1.25d0,
4317 & 1.9d0,2.d0,1.8d0,1.5d0,1.25d0,1.d0,.8d0,6*0.d0,4*0.d0,.4d0,.85d0,
4318 & 1.1d0, 1.85d0, 2.8d0, 3.d0,2.7d0,2.2d0,1.85d0,1.5d0,1.2d0,6*0.d0,
4319 & 6*0.d0, .5d0, 1.2d0, 1.7d0, 3.4d0, 5.2d0, 6.4d0, 6.1d0, 5.6d0,
4320 & 5.2d0, 6*0.d0, 2*0.d0, .0d0, 1.d0, 3.3d0, 5.2d0, 4.45d0, 3.6d0,
4321 & 2.75d0, 1.9d0, 1.65d0, 1.3d0, .95d0, .6d0, .45d0, 6*0.d0, 3*0.d0,
4322 & .0d0, .45d0, 1.4d0, 1.5d0, 1.1d0, .85d0, .5d0, .3d0, .2d0, .15d0,
4323 & 8*0.d0, 5*0.d0, .0d0, .0d0, .6d0, .8d0, .95d0, .8d0, .7d0, .6d0,
4324 & .5d0, .4d0, 6*0.d0, 5*0.d0, .0d0, .00d0, .85d0, 1.2d0, 1.4d0,
4325 & 1.2d0, 1.05d0, .9d0, .7d0, .55d0, 6*0.d0, 5*0.d0, .0d0, .00d0,
4326 & 1.d0, 1.5d0, 3.5d0, 4.15d0, 3.7d0, 2.7d0, 2.3d0, 1.75d0, 6*0.d0,
4327 & 10*0.d0, .5d0, 2.0d0, 3.3d0, 5.4d0, 7.d0 /
4329 DATA spikpu/ 0.d0, 25.d0, 13.d0, 11.d0, 10.5d0, 14.d0, 20.d0,
4330 & 20.d0, 16.d0, 14.d0, 19.d0, 28.d0, 17.5d0, 13.5d0, 12.d0, 10.5d0,
4331 & 10.d0, 10.d0, 9.5d0, 9.d0, 8.d0, 7.5d0, 7.d0, 6.5d0, 6.d0, 0.d0,
4332 & 48.d0, 19.d0, 15.d0, 11.5d0, 10.d0, 8.d0, 6.5d0, 5.5d0, 4.8d0,
4333 & 4.2d0, 7.5d0, 3.4d0, 2.5d0, 2.5d0, 2.1d0, 1.4d0, 1.d0, .8d0,
4334 & .6d0, .46d0, .3d0, .2d0, .15d0, .13d0, 11*0.d0, .95d0, .65d0,
4335 & .48d0, .35d0, .2d0, .18d0, .17d0, .16d0, .15d0, .1d0, .09d0,
4336 & .065d0, .05d0, .04d0, 12*0.d0, .2d0, .25d0, .25d0, .2d0, .1d0,
4337 & .08d0, .06d0, .045d0, .03d0, .02d0, .01d0, .005d0, .003d0,
4338 & 12*0.d0, .3d0, .24d0, .18d0, .15d0, .13d0, .12d0, .11d0, .1d0,
4339 & .09d0, .08d0, .05d0, .04d0, .03d0, 0.d0, 0.16d0, .7d0, 1.3d0,
4340 & 3.1d0, 4.5d0, 2.d0, 18*0.d0, 3*.0d0, 0.d0, 0.d0, 4.0d0, 11.d0,
4341 & 11.4d0, 10.3d0, 7.5d0, 6.8d0, 4.75d0, 2.5d0, 1.5d0, .9d0, .55d0,
4342 & .35d0, 13*0.d0, .1d0, .34d0, .5d0, .8d0, 1.1d0, 2.25d0, 3.3d0,
4343 & 2.3d0, 1.6d0, .95d0, .45d0, .28d0, .15d0, 10*0.d0, 2*0.d0, .17d0,
4344 & .64d0, 1.d0, 1.5d0, 2.1d0, 4.25d0, 6.2d0, 4.4d0, 3.d0, 1.8d0,
4345 & .9d0, .53d0, .28d0, 10*0.d0, 2*0.d0, .25d0, .82d0,
4346 & 1.3d0, 1.9d0, 2.8d0, 5.5d0 , 8.d0, 5.7d0, 3.9d0, 2.35d0, 1.15d0,
4347 & .69d0, .37d0, 10*0.d0, 7*0.d0, .0d0, .34d0, 1.5d0, 3.47d0,
4348 & 5.87d0, 6.23d0, 4.27d0, 2.6d0, 1.d0, .6d0, .3d0, .15d0, 6*0.d0/
4350 DATA spikpv/ 7*0.d0, .00d0, .16d0, .75d0, 1.73d0, 2.93d0, 3.12d0,
4351 & 2.13d0, 1.3d0, .5d0, .3d0, .15d0, .08d0, 6*0.d0, 10*0.d0, .2d0,
4352 & .6d0, .92d0, 2.4d0, 4.9d0, 6.25d0, 5.25d0, 3.5d0, 2.15d0, 1.4d0,
4353 & 1.d0, .7d0, 13*0.d0, .13d0, .4d0, .62d0, 1.6d0, 3.27d0, 4.17d0,
4354 & 3.5d0, 2.33d0, 1.43d0, .93d0, .66d0, .47d0, 13*0.d0, .07d0, .2d0,
4355 & .31d0, .8d0, 1.63d0, 2.08d0, 1.75d0, 1.17d0, .72d0, .47d0, .34d0,
4356 & .23d0, 17*0.d0, .33d0, 1.d0, 1.8d0, 2.67d0, 5.33d0, 6.d0, 5.53d0,
4357 & 5.d0, 17*0.d0, .17d0, .5d0, .9d0, 1.83d0, 2.67d0, 3.0d0, 2.77d0,
4358 & 2.5d0, 3*0.d0, 3*0.d0, 1.d0, 3.3d0, 2.8d0, 2.5d0, 2.3d0, 1.8d0,
4359 & 1.5d0, 1.1d0, .8d0, .7d0, .55d0, .3d0, 10*0.d0, 9*0.d0, .1d0,
4360 & .4d0, 1.d0, 1.4d0, 2.2d0, 2.5d0, 2.2d0, 1.65d0, 1.35d0, 1.1d0,
4361 & .8d0, .6d0, .4d0, 12*0.d0, .15d0, .6d0, 1.5d0, 2.1d0, 3.3d0,
4362 & 3.8d0, 3.3d0, 2.45d0, 2.05d0, 1.65d0, 1.2d0, .9d0, .6d0, 3*0.d0,
4363 & 9*0.d0, .10d0, .2d0, .5d0, .7d0, 1.3d0, 1.55d0, 1.9d0, 1.8d0,
4364 & 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0, 13*0.d0, .2d0, .5d0, .7d0,
4365 & 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0,
4366 & 17*0.d0, .2d0, .5d0, .85d0, 2.d0, 2.15d0, 2.05d0, 1.75d0, 1.d0,
4367 & 17*0.d0, .13d0, .33d0, .57d0, 1.33d0, 1.43d0, 1.36d0, 1.17d0,
4368 & .67d0, 17*0.d0, .07d0, .17d0, .28d0, .67d0, .72d0, .69d0, .58d0,
4369 & .33d0,17*0.d0,.4d0, .7d0, 1.d0, 1.6d0, 1.8d0, 2.3d0,1.9d0,1.7d0 /
4371 DATA spikpw/ 0.d0, 25.d0, 13.d0, 11.d0, 10.5d0, 14.d0, 2*20.d0,
4372 & 16.d0, 14.d0, 19.d0, 28.d0, 17.5d0, 13.5d0, 12.d0, 10.5d0,
4373 & 2*10.d0, 9.5d0, 9.d0, 8.d0, 7.5d0, 7.d0, 6.5d0, 6.d0, 0.d0,
4374 & 48.d0, 19.d0, 15.d0, 11.5d0, 10.d0, 8.d0, 6.5d0, 5.5d0, 4.8d0,
4375 & 4.2d0, 7.5d0, 3.4d0, 2*2.5d0, 2.1d0, 1.4d0, 1.d0, .8d0, .6d0,
4376 & .46d0, .3d0, .2d0, .15d0, .13d0, 11*0.d0, .95d0, .65d0, .48d0,
4377 & .35d0, .2d0, .18d0, .17d0, .16d0, .15d0, .1d0, .09d0, .065d0,
4378 & .05d0, .04d0, 12*0.d0, .2d0, 2*.25d0, .2d0, .1d0, .08d0, .06d0,
4379 & .045d0, .03d0, .02d0, .01d0, .005d0, .003d0, 12*0.d0, .3d0,
4380 & .24d0, .18d0, .15d0, .13d0, .12d0, .11d0, .1d0, .09d0, .08d0,
4381 & .05d0, .04d0, .03d0, 0.d0, 0.16d0, .7d0, 1.3d0, 3.1d0, 4.5d0,
4382 & 2.d0, 23*0.d0, 4.0d0, 11.d0, 11.4d0, 10.3d0, 7.5d0, 6.8d0,
4383 & 4.75d0, 2.5d0, 1.5d0, .9d0, .55d0, .35d0, 13*0.d0, .1d0, .34d0,
4384 & .5d0, .8d0, 1.1d0, 2.25d0, 3.3d0, 2.3d0, 1.6d0, .95d0, .45d0,
4385 & .28d0, .15d0, 12*0.d0, .17d0, .64d0, 1.d0, 1.5d0, 2.1d0, 4.25d0,
4386 & 6.2d0, 4.4d0, 3.d0, 1.8d0, .9d0, .53d0, .28d0, 12*0.d0, .25d0,
4387 & .82d0, 1.3d0, 1.9d0, 2.8d0, 5.5d0, 8.d0, 5.7d0, 3.9d0, 2.35d0,
4388 & 1.15d0, .69d0, .37d0, 18*0.d0, .34d0, 1.5d0, 3.47d0, 5.87d0,
4389 & 6.23d0, 4.27d0, 2.6d0, 1.d0, .6d0, .3d0, .15d0, 6*0.d0/
4391 DATA spikpx/ 8*0.d0, .16d0, .75d0, 1.73d0, 2.93d0, 3.12d0,
4392 & 2.13d0, 1.3d0, .5d0, .3d0, .15d0, .08d0, 16*0.d0, .2d0, .6d0,
4393 & .92d0, 2.4d0, 4.9d0, 6.25d0, 5.25d0, 3.5d0, 2.15d0, 1.4d0, 1.d0,
4394 & .7d0, 13*0.d0, .13d0, .4d0, .62d0, 1.6d0, 3.27d0, 4.17d0, 3.5d0,
4395 & 2.33d0, 1.43d0, .93d0, .66d0, .47d0, 13*0.d0, .07d0, .2d0, .31d0,
4396 & .8d0, 1.63d0, 2.08d0, 1.75d0, 1.17d0, .72d0, .47d0, .34d0, .23d0,
4397 & 17*0.d0, .33d0, 1.d0, 1.8d0, 2.67d0, 5.33d0, 6.d0, 5.53d0, 5.d0,
4398 & 17*0.d0, .17d0, .5d0, .9d0, 1.83d0, 2.67d0, 3.0d0, 2.77d0, 2.5d0,
4399 & 6*0.d0, 1.d0, 3.3d0, 2.8d0, 2.5d0, 2.3d0, 1.8d0, 1.5d0, 1.1d0,
4400 & .8d0, .7d0, .55d0, .3d0, 19*0.d0, .1d0, .4d0, 1.d0, 1.4d0, 2.2d0,
4401 & 2.5d0, 2.2d0, 1.65d0, 1.35d0, 1.1d0, .8d0, .6d0, .4d0, 12*0.d0,
4402 & .15d0, .6d0, 1.5d0, 2.1d0, 3.3d0, 3.8d0, 3.3d0, 2.45d0, 2.05d0,
4403 & 1.65d0, 1.2d0, .9d0, .6d0, 12*0.d0, .10d0, .2d0, .5d0, .7d0,
4404 & 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0, 1.35d0, 1.15d0, .95d0, .7d0,
4405 & 13*0.d0, .2d0, .5d0, .7d0, 1.3d0, 1.55d0, 1.9d0, 1.8d0, 1.55d0,
4406 & 1.35d0, 1.15d0, .95d0, .7d0, 17*0.d0, .2d0, .5d0, .85d0, 2.d0,
4407 & 2.15d0, 2.05d0, 1.75d0, 1.d0, 17*0.d0, .13d0, .33d0, .57d0,
4408 & 1.33d0, 1.43d0, 1.36d0, 1.17d0, .67d0, 17*0.d0, .07d0, .17d0,
4409 & .28d0, .67d0, .72d0, .69d0, .58d0, .33d0, 17*0.d0, .4d0, .7d0,
4410 & 1.d0, 1.6d0, 1.8d0, 2.3d0, 1.9d0, 1.7d0 /
4412 DATA spikp4 / 0.d0, 300.d0, 40.d0, 20.d0, 13.d0, 8.5d0, 8.d0,
4413 & 9.5d0, 12.d0, 14.d0, 15.5d0, 20.d0, 17.d0, 13.d0, 10.d0, 9.d0,
4414 & 8.5d0, 8.d0, 7.8d0, 7.3d0, 6.7d0, 9*0.d0, .23d0, .35d0, .7d0,
4415 & .52d0, .4d0, .3d0, .2d0, .15d0, .13d0, .11d0, .09d0, .07d0, 0.d0,
4416 & .033d0, .8d0, 2*1.35d0, .5d0, 19*0.d0, 0.8d0, 2.2d0, 3.6d0,
4417 & 4.6d0, 4.7d0, 3.5d0, 2.4d0, 1.8d0, 1.4d0, .75d0, .47d0, .25d0,
4418 & .13d0, .08d0, 7*0.d0, 1.2d0, 3.3d0, 5.4d0, 6.9d0, 7.3d0, 5.3d0,
4419 & 3.6d0, 2.7d0, 2.2d0, 1.1d0, .73d0, .4d0, .22d0, .12d0, 11*0.d0,
4420 & 2.0d0, 4.4d0, 6.8d0, 9.9d0, 7.9d0, 6.0d0, 3.8d0, 2.5d0, 2.d0,
4421 & 1.4d0, 1.d0, .6d0, .35d0, 10*0.d0, .25d0, .55d0, .75d0, 1.25d0,
4422 & 1.9d0, 2.d0, 1.8d0, 1.5d0, 1.25d0, 1.d0, .8d0, 10*0.d0, .4d0,
4423 & .85d0, 1.1d0, 1.85d0, 2.8d0, 3.d0, 2.7d0, 2.2d0, 1.85d0, 1.5d0,
4424 & 1.2d0, 12*0.d0, .5d0, 1.2d0, 1.7d0, 3.4d0, 5.2d0, 6.4d0, 6.1d0,
4425 & 5.6d0, 5.2d0, 9*0.d0, 1.d0, 3.3d0, 5.2d0, 4.45d0, 3.6d0, 2.75d0,
4426 & 1.9d0, 1.65d0, 1.3d0, .95d0, .6d0, .45d0, 10*0.d0, .45d0, 1.4d0,
4427 & 1.5d0, 1.1d0, .85d0, .5d0, .3d0, .2d0, .15d0, 15*0.d0, .6d0,
4428 & .8d0, .95d0, .8d0, .7d0, .6d0, .5d0, .4d0, 13*0.d0, .85d0, 1.2d0,
4429 & 1.4d0, 1.2d0, 1.05d0, .9d0, .7d0, .55d0, 13*0.d0, 1.d0, 1.5d0,
4430 & 3.5d0, 4.15d0, 3.7d0, 2.7d0, 2.3d0, 1.75d0, 16*0.d0, .5d0, 2.0d0,
4431 & 3.3d0, 5.4d0, 7.d0 /
4433 DATA spikp5/ 0.d0, 20.d0, 14.d0, 12.d0, 11.5d0, 10.d0, 8.d0,
4434 & 7.d0, 6.d0, 5.5d0, 5.3d0, 5.d0, 4.5d0, 4.4d0, 3.8d0, 3.d0, 2.8d0,
4435 & 0.d0, .5d0, 1.15d0, 2.d0, 1.3d0, .8d0, .45d0, 13*0.d0, 0.9d0,
4436 & 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0, 1.5d0, 1.2d0, .9d0, .6d0,
4437 & .45d0, .21d0, .2d0, 3*0.d0, .9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0,
4438 & 2.d0, 1.7d0, 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
4439 & 4*0.d0, 1.d0, 2.1d0, 2.6d0, 2.3d0, 2.1d0, 1.8d0, 1.7d0, 1.4d0,
4440 & 1.2d0, 1.05d0, .9d0, .66d0, .5d0, 7*0.d0, .3d0, 2*1.d0, .9d0,
4441 & .7d0, .4d0, .3d0, .2d0, 11*0.d0, .1d0, 1.d0, 2.2d0, 3.5d0, 4.2d0,
4442 & 4.55d0, 4.85d0, 4.9d0, 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0,
4443 & 1.71d0, 1.6d0, 6*0.d0, 1.4d0, 3.8d0, 5.d0, 4.7d0, 4.4d0, 4.d0,
4444 & 3.5d0, 2.85d0, 2.35d0, 2.01d0, 1.8d0, 12*0.d0, .1d0, .8d0,2.05d0,
4445 & 3.31d0, 3.5d0, 12*0.d0, .034d0, .2d0, .75d0, 1.04d0, 1.24d0 /
4453 DATA s pikp6/ 0.,6.,11.,13.,6.,5.,3.,2.2,1.5,1.2,1.,.7,.6,.5,.45,.
4454 +35,.3, 0.,6.,11.,13.,6.,5.,3.,2.2,1.5,1.2,1.,.7,.6,.5,.45,.35,.3,
4455 +0.,.5,1.3,2.8,2.3,1.6,.9,10*0., 3*0.,0.9,2.5,3.,2.5,2.3,2.,1.7,
4456 +1.5,1.2,.9,.6,.45,.21,.2, 3*0.,0.9,2.5,3.,2.5,2.3,2.,1.7,1.5,1.2,.
4457 +9,.6,.45,.21,.2, 4*0.,1.0,2.1,2.6,2.3,2.0,1.8,1.7,1.4,1.2,1.15,.9,
4458 +.66, .5, 4*0.,1.0,2.1,2.6,2.3,2.1,1.8,1.7,1.4,1.2,1.15,.9,.66, .5,
4459 +7*0.,.3,1.,1.,.9,.7,.4,.35,.2,.00,0., 7*0.,.3,1.,1.,.9,.7,.4,.35,.
4460 +2,.00,0., 9*0.,.1,1.,2.4,3.5,4.25,4.55,4.85,4.9, 9*0.,.1,1.,2.4,
4461 +3.5,4.25,4.55,4.85,4.9, 10*0.,.2,.7,1.6,2.5,2.2,1.71,1.6, 10*0.,.
4462 +2,.7,1.6,2.5,2.2,1.71,1.6, 6*0.,1.4,3.8,5.,4.7,4.4,4.,3.5,2.85,
4463 +2.35,2.01,1.8, 6*0.,1.4,3.8,5.,4.7,4.4,4.,3.5,2.85,2.35,2.01,1.8,
4464 +12*0.,.1,.8,2.05,3.31,3.5, 12*0.,.034,.20,.75,1.04,1.24, .0,2.5,
4465 +15.,21.5,15.3,3.,1.5,10*0./
4470 DATA skmpel/ 0.d0, 35.d0, 22.d0, 25.d0, 17.d0, 9.d0, 9.5d0, 8.d0,
4471 & 7.d0, 6.5d0, 6.1d0, 5.d0, 4.8d0, 4.6d0, 4.45d0, 4.3d0, 4.2d0,
4472 & 0.d0, 8.d0, 3.5d0, 8.d0, 3.d0, 1.9d0, 1.7d0, 1.d0, .9d0, .8d0,
4473 & .75d0, .5d0, .42d0, .38d0, .34d0, .25d0, .2d0,
4474 & 0.d0, 3.d0, 3.2d0, 3.5d0, 1.5d0, 1.4d0, 1.1d0, .6d0, .5d0,
4475 & .35d0, .28d0, .25d0, .18d0, .12d0, .1d0, .08d0, .04d0,
4476 & 0.d0, 8.5d0, 2.4d0, 1.7d0, 1.3d0, 1.3d0, 1.1d0, .5d0,
4477 & .4d0, .4d0, .35d0, .3d0, .28d0, .2d0, .16d0, .13d0, .11d0,
4478 & 0.d0, 7.d0, 4.8d0, 1.4d0, 1.9d0, .9d0, .4d0, .2d0, .13d0,
4479 & .1d0, .08d0, .06d0, .04d0, .02d0, .015d0, .01d0, .01d0,
4480 & 0.d0, 5.5d0, 1.d0, .8d0, .75d0, .32d0, .2d0, .1d0, .09d0,
4481 & .08d0, .065d0, .05d0, .04d0, .022d0, .017d0, 2*.01d0/
4482 DATA spikp7 / 0.d0, .56d0, 1.46d0, 3.16d0, 2.01d0, 1.28d0, .74d0,
4483 & 14*0.d0, 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0,
4484 & 1.91d0, 1.57d0, 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0, 4*0.d0,
4485 & 1.13d0, 2.61d0, 2.91d0, 2.58d0, 2.35d0, 2.02d0, 1.91d0, 1.57d0,
4486 & 1.35d0, 1.29d0, 1.01d0, .74d0, .65d0, 3*0.d0, 1.0d0, 3.03d0,
4487 & 3.36d0, 2.8d0, 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0, 1.01d0,
4488 & .67d0, .5d0, .24d0, .23d0, 3*0.d0, 1.0d0, 3.03d0, 3.36d0, 2.8d0,
4489 & 2.58d0, 2.24d0, 1.91d0, 1.68d0, 1.35d0, 1.01d0, .67d0, .5d0,
4490 & .24d0, .23d0, 7*0.d0, .34d0, 1.12d0, 1.12d0, 1.01d0, .78d0,
4491 & .45d0, .39d0, .22d0, .07d0, 0.d0, 7*0.d0, .34d0, 1.12d0, 1.12d0,
4492 & 1.01d0, .78d0, .45d0, .39d0, .22d0, .07d0, 0.d0, 6*0.d0, 1.71d0,
4493 & 4.26d0, 5.6d0, 5.57d0, 4.93d0, 4.48d0, 3.92d0, 3.19d0, 2.63d0,
4494 & 2.25d0, 2.d0, 6*0.d0, 1.71d0, 4.26d0, 5.6d0, 5.57d0, 4.93d0,
4495 & 4.48d0, 3.92d0, 3.19d0, 2.63d0, 2.25d0, 2.d0, 10*0.d0, .22d0,
4496 & .8d0, .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0, 10*0.d0, .22d0, .8d0,
4497 & .75d0, 1.d0, 1.3d0, 1.5d0, 1.3d0, 13*0.d0, .1d0, .3d0, .7d0,1.d0,
4498 & 13*0.d0, .1d0, .3d0, .7d0, 1.d0, 9*0.d0, .11d0, 1.72d0, 2.69d0,
4499 & 3.92d0, 4.76d0, 5.10d0, 5.44d0, 5.3d0, 9*0.d0, .11d0, 1.72d0,
4500 & 2.69d0, 3.92d0, 4.76d0, 5.1d0, 5.44d0, 5.3d0, 5*0.d0,9.2d0,4.7d0,
4501 & 1.9d0, 10*0.d0, 2.5d0, 15.d0, 21.5d0, 15.3d0, 3.d0, 1.5d0,
4511 *0.,4.,9.5,20.,13.,9.5,6.,4.4,3.,2.4,2.,1.4,1.2,1.,.9,.7,.6,
4512 *0.,4.5,6.,5.,2.5,2.,1.7,2.1,1.9,.9,.5,.3,.24,.2,.18,.1,.09,
4513 *0.,1.8,2.,1.1,.9,.5,.5,.4,.4,.2,.1,.06,.05,.04,.03,.02,.02,
4514 *0.,1.5,2.,.9,1.1,.4,.6,.7,.65,.3,.17,.1,.08,.07,.06,.04,.03/
4515 DATA s pikp8/ 0.,.56,1.29,2.26,1.01,.64,.37,10*0., 4*0.,1.13,2.61,
4516 +2.91,2.58,2.35,2.02,1.91,1.57,1.35,1.29,1.01,.74, .65, 3*0.,1.00,
4517 +3.03,3.36,2.8,2.58,2.24,1.91,1.68,1.35,1.01,.67,.5,.24, .23, 3*0.,
4518 +1.00,3.03,3.36,2.8,2.58,2.24,1.91,1.68,1.35,1.01,.67,.5,.24, .23,
4519 +7*0.,.34,1.12,1.12,1.01,.78,.45,.39,.22,.07,0., 6*0.,1.71,4.26,
4520 +5.6,5.57,4.93,4.48,3.92,3.19,2.63,2.25,2., 10*0.,.22,.8,.75,1.,
4521 +1.3,1.5,1.3, 13*0.,.1,.3,.7,1., 13*0.,.1,.3,.7,1., 9*0.,.11,1.72,
4522 +2.69,3.92,4.76,5.10,5.44,5.3, 4*0.,0.00,9.2,4.7,1.9,9*0. /
4532 DATA s pikp9/ 0.,24.,25.,27.,23.,21.,20.,19.,17.,15.5,14.,13.5,
4533 +13., 0.,3.6,1.7, 10*0., .0,0.0,8.7,17.7,18.8,15.9,11.7,8.,6.,5.3,
4534 +4.5,3.9,3.5, .0,.0,2.8,5.8,6.2,5.1,3.8,2.7,2.1,1.8,1.5,1.3,1.1, 4
4535 +*0.,0.0,4.6,10.2,15.1,16.9,16.5,11.,5.5,3.5, 10*0.,4.3,7.6,9., 10
4536 +*0.,1.7,2.6,3., 6*0.,.3,.6,1.,1.6,1.3,.8,.6, 6*0.,.7,1.2,1.8,2.5,
4537 +1.8,1.3,1.2, 10*0.,.6,1.4,1.7, 10*0.,1.9,4.1,5.2/
4542 DATA spikp0/ 0.d0, 24.d0, 25.d0, 27.d0, 23.d0, 21.d0, 20.d0,
4543 & 19.d0, 17.d0, 15.5d0, 14.d0, 13.5d0, 13.d0,
4544 & 0.d0, 1.8d0, .2d0, 12*0.d0,
4545 & 3.2d0, 6.05d0, 9.9d0, 5.1d0,
4546 & 3.8d0, 2.7d0, 1.9d0, 1.5d0, 1.4d0, 1.3d0, 1.1d0,
4547 & 2*.0d0, 3.2d0, 6.05d0, 9.9d0, 5.1d0,
4548 & 3.8d0, 2.7d0, 1.9d0, 1.5d0, 1.4d0, 1.3d0, 1.1d0,
4549 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0,
4550 & 16.4d0, 15.2d0, 11.d0, 5.4d0, 3.5d0,
4551 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0,
4552 & 16.4d0, 15.2d0, 11.d0, 5.4d0, 3.5d0,
4553 & 10*0.d0, .7d0, 5.1d0, 8.d0,
4554 & 10*0.d0, .7d0, 5.1d0, 8.d0,
4555 & 10*.0d0, .3d0, 2.8d0, 4.7d0,
4556 & 10*.0d0, .3d0, 2.8d0, 4.7d0,
4557 & 7*0.d0, 1.2d0, 2.5d0, 3.5d0, 6.d0, 5.3d0, 2.9d0,
4558 & 7*0.d0, 1.7d0, 3.6d0, 5.4d0, 9.d0, 7.6d0, 4.2d0,
4559 & 5*0.d0, 7.7d0, 6.1d0, 2.9d0, 5*0.d0/
4562 DATA spkp v/ 0.d0, 24.d0, 25.d0, 27.d0, 23.d0, 21.d0, 20.d0,
4563 & 19.d0, 17.d0, 15.5d0, 14.d0, 13.5d0, 13.d0,
4564 & 0.d0, 3.6d0, 1.7d0, 12*0.d0,
4565 & 8.7d0, 17.7d0, 18.8d0, 15.9d0,
4566 & 11.7d0, 8.d0, 6.d0, 5.3d0, 4.5d0, 3.9d0, 3.5d0,
4567 & .0d0, .0d0, 2.8d0, 5.8d0, 6.2d0, 5.1d0, 3.8d0,
4568 & 2.7d0, 2.1d0, 1.8d0, 1.5d0, 1.3d0, 1.1d0,
4569 & 5*0.d0, 4.6d0, 10.2d0, 15.1d0, 16.9d0, 16.5d0,
4570 & 11.d0, 5.5d0, 3.5d0,
4571 & 10*0.d0, 4.3d0, 7.6d0, 9.d0,
4572 & 10*0.d0, 1.7d0, 2.6d0, 3.d0,
4573 & 6*0.d0, .3d0, .6d0, 1.d0, 1.6d0, 1.3d0, .8d0, .6d0,
4574 & 6*0.d0, .7d0, 1.2d0, 1.8d0, 2.5d0, 1.8d0, 1.3d0,
4575 & 1.2d0, 10*0.d0, .6d0, 1.4d0, 1.7d0,
4576 & 10*0.d0, 1.9d0, 4.1d0, 5.2d0/
4578 DATA sappel/ 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0,
4579 & 50.d0, 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
4580 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0,
4581 & .05d0, .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0, .85d0,
4582 & 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0, .14d0,
4583 & .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
4584 & 0.d0, 55.d0, 50.d0, 25.d0, 15.d0, 15.d0, 14.d0, 12.d0,
4585 & 10.d0, 7.d0, 6.d0, 4.d0, 3.3d0, 2.8d0, 2.4d0, 2.d0, 1.8d0,
4586 & 1.55d0, 1.3d0, .95d0, .75d0,
4587 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0,
4588 & .25d0, .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
4589 & .01d0, .008d0, .006d0, .005d0/
4590 DATA spikpe/0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4591 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4592 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0,
4593 & 9*0.d0, 2.d0, 2.5d0, .2d0, 19*0.d0, .3d0, 1.4d0, 2.2d0, 1.2d0,
4594 & 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0, .3d0, 10*0.d0, .3d0, 1.4d0,
4595 & 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0, .3d0, 10*0.d0,
4596 & .3d0, 1.4d0, 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0, .6d0, .5d0, .4d0,
4597 & .3d0, 10*0.d0, .3d0, 1.4d0, 2.2d0, 1.2d0, 1.1d0, 1.d0, .8d0,
4598 & .6d0, .5d0, .4d0, .3d0, 9*0.d0, .6d0, 2.5d0, 5.d0, 5.2d0, 5.1d0,
4599 & 5.4d0, 5.8d0, 2.8d0, 2.1d0, 1.8d0, 1.6d0, 1.2d0, 13*0.d0, 1.3d0,
4600 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4601 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4602 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 13*0.d0, 1.3d0,
4603 & 1.5d0, 2.d0, 2.5d0, 2.5d0, 2.3d0, 1.8d0, 1.4d0, 14*0.d0, .2d0,
4604 & .5d0, 1.1d0, 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0,
4605 & 1.1d0, 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0, 1.1d0,
4606 & 1.6d0, 1.4d0, 1.1d0, .9d0, 14*0.d0, .2d0, .5d0, 1.1d0, 1.6d0,
4607 & 1.4d0, 1.1d0, .9d0, 17*0.d0, .3d0, 1.6d0, 2.6d0, 3.6d0, 17*0.d0,
4608 & .3d0, 1.6d0, 2.6d0, 3.6d0, 17*0.d0, .3d0, 1.6d0, 2.6d0,
4609 & 3.6d0, 17*0.d0, .3d0, 1.6d0, 2.6d0, 3.6d0 /
4612 & 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0,
4613 & 50.d0, 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
4614 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0,
4615 & .05d0, .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0,
4616 & .85d0, 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0,
4617 & .14d0, .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
4618 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0,
4619 & .25d0, .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
4620 & .01d0, .008d0, .006d0, .005d0 /
4621 DATA spikpz/ 0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4622 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4623 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0, 9*0.d0, 2.4d0, .2d0,
4624 & 20*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4625 & .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0,
4626 & 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0,
4627 & 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0,
4628 & 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4629 & .7d0, .5d0, .3d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4630 & 7.6d0, 6.d0, 5.d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4631 & 7.6d0, 6.d0, 5.d0, 18*0.d0, 1.d0, 4.9d0, 8.5d0, 18*0.d0, 1.d0,
4632 & 4.9d0, 8.5d0, 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0,
4633 & 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0, 15*0.d0, 1.9d0,
4634 & 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0 /
4640 & 0.d0, 176.d0, 160.d0, 105.d0, 75.d0, 68.d0, 65.d0, 50.d0,
4641 & 50.d0, 43.d0, 42.d0, 40.5d0, 35.d0, 30.d0, 28.d0,
4642 & 25.d0, 22.d0, 21.d0, 20.d0, 18.d0, 17.d0, 11*0.d0, .05d0,
4643 & .15d0, .18d0, .2d0, .2d0, .3d0, .4d0, .6d0, .7d0, .85d0,
4644 & 0.d0, 1.d0, .9d0, .46d0, .3d0, .23d0, .18d0, .16d0, .14d0,
4645 & .1d0, .08d0, .05d0, .02d0, .015d0, 4*.011d0, 3*.005d0,
4646 & 0.d0, 3.3d0, 3.d0, 1.5d0, 1.d0, .7d0, .4d0, .35d0, .4d0, .25d0,
4647 & .18d0, .08d0, .04d0, .03d0, .023d0, .016d0, .014d0,
4648 & .01d0, .008d0, .006d0, .005d0 /
4649 DATA spikpf/ 0.d0, 215.d0, 193.d0, 170.d0, 148.d0, 113.d0, 97.d0,
4650 & 84.d0, 78.d0, 68.d0, 64.d0, 61.d0, 46.d0, 36.d0, 31.3d0, 28.5d0,
4651 & 25.7d0, 22.6d0, 21.4d0, 20.7d0, 19.9d0, 9*0.d0, 2.4d0, .2d0,
4652 & 20*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4653 & .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0,
4654 & 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0, 10*0.d0, 1.8d0, 2.8d0,
4655 & 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0, .7d0, .5d0, .3d0,
4656 & 10*0.d0, 1.8d0, 2.8d0, 3.6d0, 2.3d0, 1.8d0, 1.5d0, 1.3d0, 1.d0,
4657 & .7d0, .5d0, .3d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4658 & 7.6d0, 6.d0, 5.d0, 13*0.d0, 5.2d0, 8.7d0, 11.4d0, 14.d0, 11.9d0,
4659 & 7.6d0, 6.d0, 5.d0, 18*0.d0, 1.d0, 4.9d0, 8.5d0, 18*0.d0, 1.d0,
4660 & 4.9d0, 8.5d0, 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0,
4661 & 15*0.d0, 1.9d0, 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0, 15*0.d0, 1.9d0,
4662 & 2.3d0, 4.d0, 6.5d0, 5.2d0, 3.4d0 /
4664 DATA spkp15/0.d0, 20.d0, 14.d0, 12.d0, 11.5d0, 10.d0, 8.d0, 7.d0,
4665 & 6.d0, 5.5d0, 5.3d0, 5.d0, 4.5d0, 4.4d0, 3.8d0, 3.d0, 2.8d0,
4666 & 0.d0, .5d0, 1.15d0, 2.d0, 1.3d0, .8d0, .45d0, 10*0.d0,
4667 & 3*0.d0, 0.9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0,
4668 & 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
4669 & 3*0.d0, 0.9d0, 2.5d0, 3.d0, 2.5d0, 2.3d0, 2.d0, 1.7d0,
4670 & 1.5d0, 1.2d0, .9d0, .6d0, .45d0, .21d0, .2d0,
4671 & 4*0.d0, 1.d0, 2.1d0, 2.6d0, 2.3d0, 2.1d0, 1.8d0, 1.7d0,
4672 & 1.4d0, 1.2d0, 1.05d0, .9d0, .66d0, .5d0,
4673 & 7*0.d0, .3d0, 1.d0, 1.d0, .9d0, .7d0, .4d0, .30d0, .2d0,
4674 & 11*0.d0, .1d0, 1.d0, 2.2d0, 3.5d0, 4.20d0, 4.55d0,
4676 & 10*0.d0, .2d0, .7d0, 1.6d0, 2.5d0, 2.2d0, 1.71d0, 1.6d0,
4677 & 6*0.d0, 1.4d0, 3.8d0, 5.d0, 4.7d0, 4.4d0, 4.d0, 3.5d0,
4678 & 2.85d0, 2.35d0, 2.01d0, 1.8d0,
4679 & 12*0.d0, .1d0, .8d0, 2.05d0, 3.31d0, 3.5d0,
4680 & 12*0.d0, .034d0, .20d0, .75d0, 1.04d0, 1.24d0 /
4686 *0.,4.,9.5,20.,13.,9.5,6.,4.4,3.,2.4,2.,1.4,1.2,1.,.9,.7,.6,
4687 *0.,4.5,6.,5.,2.5,2.,1.7,2.1,1.9,.9,.5,.3,.24,.2,.18,.1,.09,
4688 *0.,1.8,2.,1.1,.9,.5,.5,.4,.4,.2,.1,.06,.05,.04,.03,.02,.02,
4689 *0.,1.5,2.,.9,1.1,.4,.6,.7,.65,.3,.17,.1,.08,.07,.06,.04,.03,
4690 *0.,.56,1.29,2.26,1.01,.64,.37,10*0.,
4691 *4*0.,1.13,2.61,2.91,2.58,2.35,2.02,1.91,1.57,1.35,1.29,1.01,.74,
4693 *3*0.,1.00,3.03,3.36,2.8,2.58,2.24,1.91,1.68,1.35,1.01,.67,.5,.24,
4695 *3*0.,1.00,3.03,3.36,2.8,2.58,2.24,1.91,1.68,1.35,1.01,.67,.5,.24,
4697 *7*0.,.34,1.12,1.12,1.01,.78,.45,.39,.22,.07,0.,
4698 *6*0.,1.71,4.26,5.6,5.57,4.93,4.48,3.92,3.19,2.63,2.25,2.,
4699 *10*0.,.22,.8,.75,1.,1.3,1.5,1.3,
4700 *13*0.,.1,.3,.7,1., 13*0.,.1,.3,.7,1.,
4701 *9*0.,.11,1.72,2.69,3.92,4.76,5.10,5.44,5.3,
4702 *4*0.,0.00,9.2,4.7,1.9,9*0.
4704 DATA nure/9,12,5*0,10,14,3*0,1,3,5,7,6*0,2,6,16,5*0,
4705 *10,13,5*0,11,12,3*0,2,4,6,8,6*0,3,15,7,5*0/
4713 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4723 dimension dclin(195),dchn(143),dchna(36),dchnb(60)
4724 dimension pdci(60),pdch(55)
4728 DATA (dclin(i),i=1,80) /
4730 * 5.000
d-01, 1.000
d+00, 0.000
d+00, 1.000
d+00, 0.000
d+00,
4731 * 4.993
d-01, 9.881
d-01, 5.963
d-02, 9.851
d-01, 5.945
d-02,
4732 * 4.936
d-01, 8.955
d-01, 5.224
d-01, 8.727
d-01, 5.091
d-01,
4733 * 4.889
d-01, 8.228
d-01, 8.859
d-01, 7.871
d-01, 8.518
d-01,
4734 * 4.874
d-01, 7.580
d-01, 1.210
d+00, 7.207
d-01, 1.117
d+00,
4735 * 4.912
d-01, 6.969
d-01, 1.516
d+00, 6.728
d-01, 1.309
d+00,
4736 * 5.075
d-01, 6.471
d-01, 1.765
d+00, 6.667
d-01, 1.333
d+00,
4737 * 5.383
d-01, 6.054
d-01, 1.973
d+00, 7.059
d-01, 1.176
d+00,
4738 * 5.397
d-01, 5.990
d-01, 2.005
d+00, 7.023
d-01, 1.191
d+00,
4739 * 5.336
d-01, 6.083
d-01, 1.958
d+00, 6.959
d-01, 1.216
d+00,
4740 * 5.317
d-01, 6.075
d-01, 1.962
d+00, 6.897
d-01, 1.241
d+00,
4741 * 5.300
d-01, 6.016
d-01, 1.992
d+00, 6.786
d-01, 1.286
d+00,
4742 * 5.281
d-01, 6.063
d-01, 1.969
d+00, 6.786
d-01, 1.286
d+00,
4743 * 5.280
d-01, 5.960
d-01, 2.020
d+00, 6.667
d-01, 1.333
d+00,
4744 * 5.273
d-01, 5.920
d-01, 2.040
d+00, 6.604
d-01, 1.358
d+00,
4745 * 5.273
d-01, 5.862
d-01, 2.069
d+00, 6.538
d-01, 1.385
d+00/
4746 DATA (dclin(i),i=81,160) /
4748 * 5.223
d-01, 5.980
d-01, 2.814
d+00, 6.538
d-01, 1.385
d+00,
4749 * 5.202
d-01, 5.969
d-01, 2.822
d+00, 6.471
d-01, 1.412
d+00,
4750 * 5.183
d-01, 5.881
d-01, 2.883
d+00, 6.327
d-01, 1.469
d+00,
4751 * 5.159
d-01, 5.866
d-01, 2.894
d+00, 6.250
d-01, 1.500
d+00,
4752 * 5.133
d-01, 5.850
d-01, 2.905
d+00, 6.170
d-01, 1.532
d+00,
4753 * 5.106
d-01, 5.833
d-01, 2.917
d+00, 6.087
d-01, 1.565
d+00,
4754 * 5.084
d-01, 5.801
d-01, 2.939
d+00, 6.000
d-01, 1.600
d+00,
4755 * 5.063
d-01, 5.763
d-01, 2.966
d+00, 5.909
d-01, 1.636
d+00,
4756 * 5.036
d-01, 5.730
d-01, 2.989
d+00, 5.814
d-01, 1.674
d+00,
4757 * 5.014
d-01, 5.683
d-01, 3.022
d+00, 5.714
d-01, 1.714
d+00,
4758 * 4.986
d-01, 5.641
d-01, 3.051
d+00, 5.610
d-01, 1.756
d+00,
4759 * 4.964
d-01, 5.580
d-01, 3.094
d+00, 5.500
d-01, 1.800
d+00,
4760 * 4.936
d-01, 5.573
d-01, 3.099
d+00, 5.431
d-01, 1.827
d+00,
4761 * 4.909
d-01, 5.509
d-01, 3.144
d+00, 5.313
d-01, 1.875
d+00,
4762 * 4.885
d-01, 5.512
d-01, 3.142
d+00, 5.263
d-01, 1.895
d+00,
4763 * 4.857
d-01, 5.437
d-01, 3.194
d+00, 5.135
d-01, 1.946
d+00/
4764 DATA (dclin(i),i=161,195) /
4765 * 4.830
d-01, 5.353
d-01, 3.253
d+00, 5.000
d-01, 2.000
d+00,
4766 * 4.801
d-01, 5.323
d-01, 3.274
d+00, 4.915
d-01, 2.034
d+00,
4767 * 4.770
d-01, 5.228
d-01, 3.341
d+00, 4.767
d-01, 2.093
d+00,
4768 * 4.738
d-01, 5.156
d-01, 3.391
d+00, 4.643
d-01, 2.143
d+00,
4769 * 4.701
d-01, 5.010
d-01, 3.493
d+00, 4.444
d-01, 2.222
d+00,
4770 * 4.672
d-01, 4.990
d-01, 3.507
d+00, 4.375
d-01, 2.250
d+00,
4771 * 4.634
d-01, 4.856
d-01, 3.601
d+00, 4.194
d-01, 2.323
d+00/
4774 * 4.400
d+02, 1.896
d-01, 1.931
d-01, 1.982
d-01, 1.015
d-01,
4775 * 1.029
d-01, 4.180
d-02, 4.228
d-02, 4.282
d-02, 4.350
d-02,
4776 * 2.204
d-02, 2.236
d-02, 5.900
d+02, 1.433
d-01, 1.555
d-01,
4777 * 1.774
d-01, 1.000
d-01, 1.128
d-01, 5.132
d-02, 5.600
d-02,
4778 * 6.158
d-02, 6.796
d-02, 3.660
d-02, 3.820
d-02, 6.500
d+02,
4779 * 1.192
d-01, 1.334
d-01, 1.620
d-01, 9.527
d-02, 1.141
d-01,
4780 * 5.283
d-02, 5.952
d-02, 6.765
d-02, 7.878
d-02, 4.796
d-02,
4781 * 6.957
d-02, 8.000
d+02, 4.872
d-02, 6.694
d-02, 1.152
d-01,
4782 * 9.348
d-02, 1.368
d-01, 6.912
d-02, 7.953
d-02, 9.577
d-02,
4783 * 1.222
d-01, 7.755
d-02, 9.525
d-02, 1.000
d+03, 3.997
d-02,
4784 * 5.456
d-02, 9.804
d-02, 8.084
d-02, 1.208
d-01, 6.520
d-02,
4785 * 8.233
d-02, 1.084
d-01, 1.474
d-01, 9.328
d-02, 1.093
d-01/
4788 * 1.000
d+03, 9.453
d-02, 9.804
d-02, 8.084
d-02, 1.208
d-01,
4789 * 6.520
d-02, 8.233
d-02, 1.084
d-01, 1.474
d-01, 9.328
d-02,
4790 * 1.093
d-01, 1.400
d+03, 1.072
d-01, 7.450
d-02, 6.645
d-02,
4791 * 1.136
d-01, 6.750
d-02, 8.580
d-02, 1.110
d-01, 1.530
d-01,
4792 * 1.010
d-01, 1.350
d-01, 2.170
d+03, 4.004
d-02, 3.013
d-02,
4793 * 2.664
d-02, 5.511
d-02, 4.240
d-02, 7.660
d-02, 1.364
d-01,
4794 * 2.300
d-01, 1.670
d-01, 2.010
d-01, 2.900
d+03, 1.870
d-02,
4795 * 1.804
d-02, 1.320
d-02, 2.970
d-02, 2.860
d-02, 5.160
d-02,
4796 * 1.020
d-01, 2.400
d-01, 2.250
d-01, 3.370
d-01, 4.400
d+03,
4797 * 1.196
d-03, 8.784
d-03, 1.517
d-02, 2.874
d-02, 2.488
d-02,
4798 * 4.464
d-02, 8.330
d-02, 2.008
d-01, 2.360
d-01, 3.567
d-01/
4800 DATA (dchn(i),i=1,90) /
4801 * 4.770
d-01, 4.750
d-01, 4.715
d-01, 4.685
d-01, 4.650
d-01,
4802 * 4.610
d-01, 4.570
d-01, 4.550
d-01, 4.500
d-01, 4.450
d-01,
4803 * 4.405
d-01, 4.350
d-01, 4.300
d-01, 4.250
d-01, 4.200
d-01,
4804 * 4.130
d-01, 4.060
d-01, 4.000
d-01, 3.915
d-01, 3.840
d-01,
4805 * 3.760
d-01, 3.675
d-01, 3.580
d-01, 3.500
d-01, 3.400
d-01,
4806 * 3.300
d-01, 3.200
d-01, 3.100
d-01, 3.000
d-01, 2.900
d-01,
4807 * 2.800
d-01, 2.700
d-01, 2.600
d-01, 2.500
d-01, 2.400
d-01,
4808 * 2.315
d-01, 2.240
d-01, 2.150
d-01, 2.060
d-01, 2.000
d-01,
4809 * 1.915
d-01, 1.850
d-01, 1.780
d-01, 1.720
d-01, 1.660
d-01,
4810 * 1.600
d-01, 1.550
d-01, 1.500
d-01, 1.450
d-01, 1.400
d-01,
4811 * 1.360
d-01, 1.320
d-01, 1.280
d-01, 1.250
d-01, 1.210
d-01,
4812 * 1.180
d-01, 1.150
d-01, 1.120
d-01, 1.100
d-01, 1.070
d-01,
4813 * 1.050
d-01, 1.030
d-01, 1.010
d-01, 9.900
d-02, 9.700
d-02,
4814 * 9.550
d-02, 9.480
d-02, 9.400
d-02, 9.200
d-02, 9.150
d-02,
4815 * 9.100
d-02, 9.000
d-02, 8.990
d-02, 8.900
d-02, 8.850
d-02,
4816 * 8.750
d-02, 8.700
d-02, 8.650
d-02, 8.550
d-02, 8.500
d-02,
4817 * 8.499
d-02, 8.450
d-02, 8.350
d-02, 8.300
d-02, 8.250
d-02,
4818 * 8.150
d-02, 8.100
d-02, 8.030
d-02, 8.000
d-02, 7.990
d-02/
4819 DATA (dchn(i),i=91,143) /
4820 * 7.980
d-02, 7.950
d-02, 7.900
d-02, 7.860
d-02, 7.800
d-02,
4821 * 7.750
d-02, 7.650
d-02, 7.620
d-02, 7.600
d-02, 7.550
d-02,
4822 * 7.530
d-02, 7.500
d-02, 7.499
d-02, 7.498
d-02, 7.480
d-02,
4823 * 7.450
d-02, 7.400
d-02, 7.350
d-02, 7.300
d-02, 7.250
d-02,
4824 * 7.230
d-02, 7.200
d-02, 7.100
d-02, 7.050
d-02, 7.020
d-02,
4825 * 7.000
d-02, 6.999
d-02, 6.995
d-02, 6.993
d-02, 6.991
d-02,
4826 * 6.990
d-02, 6.870
d-02, 6.850
d-02, 6.800
d-02, 6.780
d-02,
4827 * 6.750
d-02, 6.700
d-02, 6.650
d-02, 6.630
d-02, 6.600
d-02,
4828 * 6.550
d-02, 6.525
d-02, 6.510
d-02, 6.500
d-02, 6.499
d-02,
4829 * 6.498
d-02, 6.496
d-02, 6.494
d-02, 6.493
d-02, 6.490
d-02,
4830 * 6.488
d-02, 6.485
d-02, 6.480
d-02/
4833 * 6.300
d+02, 7.810
d-02, 1.421
d-01, 1.979
d-01, 2.479
d-01,
4834 * 3.360
d-01, 5.400
d-01, 7.236
d-01, 1.000
d+00, 1.540
d+03,
4835 * 2.225
d-01, 3.950
d-01, 5.279
d-01, 6.298
d-01, 7.718
d-01,
4836 * 9.405
d-01, 9.835
d-01, 1.000
d+00, 2.560
d+03, 2.625
d-01,
4837 * 4.550
d-01, 5.963
d-01, 7.020
d-01, 8.380
d-01, 9.603
d-01,
4838 * 9.903
d-01, 1.000
d+00, 3.520
d+03, 4.250
d-01, 6.875
d-01,
4839 * 8.363
d-01, 9.163
d-01, 9.828
d-01, 1.000
d+00, 1.000
d+00,
4843 * 6.300
d+02, 3.800
d-02, 7.164
d-02, 1.275
d-01, 2.171
d-01,
4844 * 3.227
d-01, 4.091
d-01, 5.051
d-01, 6.061
d-01, 7.074
d-01,
4845 * 8.434
d-01, 1.000
d+00, 2.040
d+03, 1.200
d-01, 2.115
d-01,
4846 * 3.395
d-01, 5.295
d-01, 7.251
d-01, 8.511
d-01, 9.487
d-01,
4847 * 9.987
d-01, 1.000
d+00, 1.000
d+00, 1.000
d+00, 2.200
d+03,
4848 * 1.344
d-01, 2.324
d-01, 3.754
d-01, 5.674
d-01, 7.624
d-01,
4849 * 8.896
d-01, 9.808
d-01, 1.000
d+00, 1.000
d+00, 1.000
d+00,
4850 * 1.000
d+00, 2.850
d+03, 2.330
d-01, 4.130
d-01, 6.610
d-01,
4851 * 9.010
d-01, 9.970
d-01, 1.000
d+00, 1.000
d+00, 1.000
d+00,
4852 * 1.000
d+00, 1.000
d+00, 1.000
d+00, 3.500
d+03, 3.300
d-01,
4853 * 5.450
d-01, 7.950
d-01, 1.000
d+00, 1.000
d+00, 1.000
d+00,
4854 * 1.000
d+00, 1.000
d+00, 1.000
d+00, 1.000
d+00, 1.000
d+00/
4861 IF (
ekin.GT.3.5d0)
RETURN
4863 IF(kproj.EQ.8) goto 101
4864 IF(kproj.EQ.1) goto 102
4867 &
' INVALID PARTICLE TYPE IN DNUPRE - KPROJ=',kproj,
4868 &
' COS(THETA) = 1D0 RETURNED'
4872 IF (
ekin.GT.0.740d0)goto 1000
4873 IF (
ekin.LT.0.300d0)
THEN
4882 ie=abs(ener/0.020d0)
4883 univ=(ener-float(ie)*0.020d0)/0.020d0
4886 bwfw=(dclin(k+5)-dclin(k))*univ + dclin(k)
4887 IF (
rndm(v).LT.bwfw)
THEN
4895 coef=(dclin(k+5)-dclin(k))*univ + dclin(k)
4907 IF(value2.GT.0.0)
THEN
4908 cst=
max(r1,r2,r3,r4)
4914 cst=-
max(r1,r2,r3,r4,r5)
4918 cst=-
max(r1,r2,r3,r4,r5,r6,r7)
4929 1000 ener=
ekin - 0.66
4934 univ=(ener-float(ie)*0.020d0)/0.020d0
4936 bwfw=(dchn(k+1)-dchn(k))*univ + dchn(k)
4939 IF (rnd.GE.bwfw)
THEN
4941 IF (dchna(k).GT.emev)
THEN
4942 unive=(emev-dchna(k-9))/(dchna(k)-dchna(k-9))
4946 p=(dchna(ii)-dchna(ii-9))*unive + dchna(ii-9)
4951 goto(290,290,290,290,330,340,350,360) i
4960 IF (dchnb(k).GT.emev)
THEN
4961 unive=(emev-dchnb(k-12))/(dchnb(k)-dchnb(k-12))
4965 p=(dchnb(ii)-dchnb(ii-12))*unive + dchnb(ii-12)
4970 goto(120,120,140,150,160,160,180,190,200,210,220) i
4977 120 cst=1.0
d-2*flti-1.0d0
4979 140 cst=2.0
d-2*univ-0.98d0
4981 150 cst=4.0
d-2*univ-0.96d0
4983 160 cst=6.0
d-2*flti-1.16d0
4985 180 cst=8.0
d-2*univ-0.80d0
4987 190 cst=1.0
d-1*univ-0.72d0
4989 200 cst=1.2
d-1*univ-0.62d0
4991 210 cst=2.0
d-1*univ-0.50d0
4993 220 cst=3.0
d-1*(univ-1.0d0)
4996 290 cst=1.0d0-2.5
d-2*flti
4998 330 cst=0.85d0+0.5
d-1*univ
5000 340 cst=0.70d0+1.5
d-1*univ
5002 350 cst=0.50d0+2.0
d-1*univ
5013 IF (
ekin.LE.0.500d0)
THEN
5018 ELSEIF (
ekin.LT.1.0d0)
THEN
5020 IF (pdci(k).GT.emev)
THEN
5021 unive=(emev-pdci(k-12))/(pdci(k)-pdci(k-12))
5026 sum=sum + (pdci(ii)-pdci(ii-12))*unive + pdci(ii-12)
5028 IF (univ.LT.sum)
THEN
5031 goto(55,55,55,60,60,65,65,65,65,70,70) i
5038 IF (pdch(k).GT.emev)
THEN
5039 unive=(emev-pdch(k-11))/(pdch(k)-pdch(k-11))
5044 sum=sum + (pdch(ii)-pdch(ii-11))*unive + pdch(ii-11)
5046 IF (univ.LT.sum)
THEN
5049 goto(50,55,60,60,65,65,65,65,70,70) i
5060 60 cst=0.3d0+0.1d0*flti
5062 65 cst=0.6d0+0.04d0*flti
5064 70 cst=0.78d0+0.02d0*flti
5067 IF (
rndm(v).GT.0.5d0) cst=-cst
5084 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5093 COMMON /paname/ btype(30)
5105 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
5106 +iibar(210),k1(210),k2(210)
5109 parameter(npa=45,ipmax=16)
5110 dimension kproj(ipmax)
5111 dimension
p(npa),ek(npa),pl(180)
5112 dimension sigma(npa,4)
5113 DATA p /0.13, 0.19, 0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0,
5114 +1.1, 1.2, 1.3, 1.4, 1.5, 2.0, 3.0, 4.0, 5.0, 6.0, 10.0, 20., 50.,
5115 +100., 200., 400., 1000., 5000., 10000.,
5116 +20000.,50000.,100000.,200000.,500000.,1000000.,2000000.,
5117 +5000000.,10000000.,20000000.,50000000.,100000000.,
5118 +200000000.,500000000.,1000000000./
5119 DATA kproj / 1, 8, 13, 14, 15, 16, 24, 25, 2, 9, 17, 18, 20, 21,
5123 WRITE(6,
'(2A)')
' TEST OF CROSS SECTIONS ROUTINES APPLIED:',
5124 +
' SHPTOT, SIHNIN, SIHNEL(SIGEL)'
5125 WRITE(6,
'(A,5X,A)')
' PRINTED FOR EACH MOMENTUM/KINETIC ENERGY:',
5126 +
' SIHNEL , SIHNIN / (SIHNEL+SIHNIN), SHPTOT'
5133 ek(ie)=
sqrt(
p(ie)**2+aam(kpro)**2) - aam(kpro)
5135 CALL
sihnel(kpro,1,
p(ie),siel)
5137 CALL
sihnin(kpro,1,
p(ie),siin)
5140 sigma(ie,4)=siin + siel
5141 sigma(ie,3)=
dshpto(kpro,
p(ie))
5144 pl(iie)=log10(
p(ie))
5148 WRITE(6,
'(1H1,A,I3,2A)')
' PROJECTILE IPRO=', kpro,
' - ',
5150 WRITE(6,
'(//A/2A/)')
5151 +
' TABLE OF CROSS SECTION VALUES VERSUS PLAB / EKIN:',
5152 +
' PLAB, EKIN / SIG(EL), SIG(INEL) /',
5153 +
' SIG(TOT)-USED,SIG(EL)+SIG(INEL) - TEST PRINT'
5156 WRITE(6,1000)
p(ie),ek(ie), (sigma(ie,is),is=1,4)
5161 sigma(ie,im)=log10(sigma(ie,im))
5164 +
' DOUBLE-LOG PLOT OF SIGMA(EL/INEL/TOT/SUM EL+INEL) VS PLAB'
5166 CALL
plot(pl,sigma,npoi,4,npa,-1.0d0,0.2d0,0.0d0,0.050d0)
5170 1000
FORMAT(2(1pe10.2),5
x,4(1pe12.3))