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),5x,4(1pe12.3))
subroutine dtwopa(E1, E2, P1, P2, COD1, COD2, COF1, COF2, SIF1, SIF2, IT1, IT2, UMOO, ECM, P, N, AM1, AM2)
subroutine dhadri(N, PLAB, ELAB, CX, CY, CZ, ITTA)
subroutine dtwopd(UMO, ECM1, ECM2, PCM1, PCM2, COD1, COF1, SIF1, COD2, COF2, SIF2, AM1, AM2)
subroutine dtrafo(GAM, BGAM, CX, CY, CZ, COD, COF, SIF, P, ECM, PL, CXL, CYL, CZL, EL)
G4double ekin(const G4LorentzVector &p) const
subroutine dsige(IT, AA, POO, SEL, ZL)
double precision function dshnto(KPROJ, KTARG, UMO)
subroutine sihnin(IPROJ, ITAR, PO, SIIN)
subroutine dsigin(IRE, PLAB, N, IE, AMT, AMN, ECM, SI, ITAR)
subroutine dtchoi(T, P, PP, E, EE, I, II, N, AM1, AM2)
subroutine plot(X, Y, N, M, MM, XO, DX, YO, DY)
double precision function dshpto(IT, PO)
subroutine fhad(IPRMOD, IPRO, PLAB, ELAB, CX, CY, CZ, ITHKK, ITTA, IELINE, IREJFH)
subroutine daltra(GA, BGX, BGY, BGZ, PCX, PCY, PCZ, EC, P, PX, PY, PZ, E)
double precision function rexp(W)
double precision function rndm(RDUMMY)
subroutine elhain(IP, PLA, ELAB, CX, CY, CZ, IT, IREJ)
subroutine dcosi(SFE, CFE)
subroutine dcalum(N, ITTA)
double precision function damg(IT)
subroutine dsfecf(SFE, CFE)
double precision function xshpto(IT, PO)
static c2_log_p< float_type > & log()
make a *new object
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
static c2_sqrt_p< float_type > & sqrt()
make a *new object
subroutine drtran(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
subroutine dsihae(KPROJ, EKIN, PLAB, ANUC, SIGELA)
subroutine dgauss(X, A, S)
subroutine dthrep(UMO, ECM1, ECM2, ECM3, PCM1, PCM2, PCM3, COD1, COF1, SIF1, COD2, COF2, SIF2, COD3, COF3, SIF3, AM1, AM2, AM3)
subroutine tsamcs(KPROJ, EKIN, CST)
subroutine sihnel(IPROJ, ITAR, POO, SIEL)
static c2_exp_p< float_type > & exp()
make a *new object