* * $Id$ * * $Log$ * Revision 1.1.1.1 1995/10/24 10:19:35 cernlib * Geant * * #include "geant321/pilot.h" *CMZ : 3.21/02 29/03/94 15.41.45 by S.Giani *-- Author : *$ CREATE BLKDT3.FOR *COPY BLKDT3 * *=== blkdt3 ===========================================================* *== * SUBROUTINE FLKDT3 #include "geant321/dblprc.inc" #include "geant321/dimpar.inc" #include "geant321/iounit.inc" * * * reaction channel cross section data * #include "geant321/reac.inc" DIMENSION & UMOPI(92), UMOKC(68), UMOP(39), UMON(63), UMOK0(34), & PLAPI(92), PLAKC(68), PLAP(39), PLAN(63), PLAK0(34), & SPIKP1(315), SPIKPU(278), SPIKPV(372), & SPIKPW(278), SPIKPX(372), SPIKP4(315), & SPIKP5(187), SPIKP6(289), & SKMPEL(102), SPIKP7(289), SKMNEL(68), SPIKP8(187), & SPIKP9(143), SPIKP0(169), SPKPV(143), & SAPPEL(105), SPIKPE(399), SAPNEL(84), SPIKPZ(273), & SANPEL(84) , SPIKPF(273), & SPKP15(187), SPKP16(272), & NRKPI(164), NRKKC(132), NRKP(70), NRKN(116), NRKK0(54), & NURELN(60) * DIMENSION NRKLIN(532) EQUIVALENCE (NRK(1,1), NRKLIN(1)) EQUIVALENCE ( UMO( 1), UMOPI(1)), ( UMO( 93), UMOKC(1)) EQUIVALENCE ( UMO(161), UMOP(1)), ( UMO(200), UMON(1)) EQUIVALENCE ( UMO(263), UMOK0(1)) EQUIVALENCE ( PLABF( 1), PLAPI(1)), ( PLABF( 93), PLAKC(1)) EQUIVALENCE ( PLABF(161), PLAP(1)), ( PLABF(200), PLAN(1)) EQUIVALENCE ( PLABF(263), PLAK0(1)) EQUIVALENCE ( WK( 1), SPIKP1(1)), ( WK( 316), SPIKPU(1)) EQUIVALENCE ( WK( 594), SPIKPV(1)), ( WK( 966), SPIKPW(1)) EQUIVALENCE ( WK(1244), SPIKPX(1)), ( WK(1616), SPIKP4(1)) EQUIVALENCE ( WK(1931), SPIKP5(1)), ( WK(2118), SPIKP6(1)) EQUIVALENCE ( WK(2407), SKMPEL(1)), ( WK(2509), SPIKP7(1)) EQUIVALENCE ( WK(2798), SKMNEL(1)), ( WK(2866), SPIKP8(1)) EQUIVALENCE ( WK(3053), SPIKP9(1)), ( WK(3196), SPIKP0(1)) EQUIVALENCE ( WK(3365), SPKPV(1)), ( WK(3508), SAPPEL(1)) EQUIVALENCE ( WK(3613), SPIKPE(1)), ( WK(4012), SAPNEL(1)) EQUIVALENCE ( WK(4096), SPIKPZ(1)), ( WK(4369), SANPEL(1)) EQUIVALENCE ( WK(4453), SPIKPF(1)), ( WK(4726), SPKP15(1)) EQUIVALENCE ( WK(4913), SPKP16(1)) EQUIVALENCE (NRK(1,1), NRKLIN(1)) EQUIVALENCE (NRKLIN( 1), NRKPI(1)), (NRKLIN( 165), NRKKC(1)) EQUIVALENCE (NRKLIN( 297), NRKP(1)), (NRKLIN( 367), NRKN(1)) EQUIVALENCE (NRKLIN( 483), NRKK0(1)) EQUIVALENCE (NURE(1,1), NURELN(1)) * DIMENSION UMOPI1(92), UMOKC1(68), UMOP1(39), UMON1(63), & UMOK01(34), PLAPI1(92), PLAKC1(68), PLAP1(39), PLAN1(63), & PLAK01(34), SIIN1(296), & SPKP11(315), SPKPU1(278), CPKPV1(372), & SPKPW1(278), SPKPX1(372), SPKP41(315), & SPKP51(187), SPKP61(289), & SKMPL1(102), SPKP71(289), SKMNL1(68), SPKP81(187), & SPKP91(143), SPKP01(169), SPKPV1(143), & SAPPL1(105), SPKPE1(399), SAPNL1(84), SPKPZ1(273), & SANPL1(84) , SPKPF1(273), & CPKP51(187), CPKP61(272), & NRKPI1(164), NRKKC1(132), NRKP1(70), NRKN1(116), & NRKK01(54), NURE1(60) * **** pi- p data * **** pi+ n data * DATA PLAPI1/ 0.D0, .3D0, .5D0, .6D0, .7D0, .8D0, .9D0, .95D0,1.D0, & 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0, 2.D0, 2.3D0, 2.5D0, 2.8D0, & 3.D0, 3.5D0, 4.D0, 0.D0, .285D0, .4D0, .45D0, .5D0, .6D0, .7D0, & .75D0, .8D0, .85D0, .9D0, 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0, & 1.8D0, 2.D0, 2.3D0, 2.5D0, 2.8D0, 3.D0, 3.5D0, 4.D0, 4.5D0, 0.D0, & .285D0, .4D0, .45D0, .5D0, .6D0, .7D0, .75D0, .8D0, .85D0, .9D0, & 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0, 2.D0, 2.3D0, 2.5D0, & 2.8D0, 3.D0, 3.5D0, 4.D0, 4.5D0, 0.D0, .3D0, .5D0, .6D0, .7D0, & .8D0, .9D0, .95D0, 1.D0, 1.15D0, 1.3D0, 1.5D0, 1.6D0, 1.8D0, & 2.D0, 2.3D0, 2.5D0, 2.8D0, 3.D0, 3.5D0, 4.D0 / DATA PLAKC1/ & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0, & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0, & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0, & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0/ DATA PLAK01/ & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0, & 0.D0, .58D0, .8D0, 1.01D0, 1.23D0, 1.45D0, 1.68D0, 1.94D0, & 2.18D0, 2.42D0, 2.68D0, 2.96D0, 3.24D0, & 3.51D0, 3.84D0, 4.16D0, 4.49D0/ * pp pn np nn * DATA PLAP1/ & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0, & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0, & 0.D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0 / * app apn anp ann * DATA PLAN1/ & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0, & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0, & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0, & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0, & 0.D0, 1.D-3, .1D0, .2D0, .3D0, .4D0, .5D0, .6D0, & .74D0, 1.06D0, 1.34D0, 1.63D0, 1.92D0, 2.2D0, 2.5D0,2.8D0,3.1D0, & 3.43D0, 3.75D0, 4.07D0, 4.43D0 / DATA SIIN1/ 296*0.D0 / DATA UMOPI1/ 1.08D0,1.233D0,1.302D0,1.369D0,1.496D0, & 1.557D0,1.615D0,1.6435D0, & 1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,2.071D0,2.159D0, & 2.286D0,2.366D0,2.482D0,2.56D0, & 2.735D0,2.90D0, & 1.08D0,1.222D0,1.302D0,1.3365D0,1.369D0,1.434D0, & 1.496D0,1.527D0,1.557D0, & 1.586D0,1.615D0,1.672D0,1.753D0,1.831D0,1.930D0,1.978D0, & 2.071D0,2.159D0,2.286D0,2.366D0, & 2.482D0,2.560D0,2.735D0,2.90D0,3.06D0, & 1.08D0,1.222D0,1.302D0,1.3365D0,1.369D0,1.434D0, & 1.496D0,1.527D0,1.557D0, & 1.586D0,1.615D0,1.672D0,1.753D0,1.831D0,1.930D0,1.978D0, & 2.071D0,2.159D0,2.286D0,2.366D0, & 2.482D0,2.560D0,2.735D0,2.90D0,3.06D0, & 1.08D0,1.233D0,1.302D0,1.369D0,1.496D0, & 1.557D0,1.615D0,1.6435D0, & 1.672D0,1.753D0,1.831D0,1.930D0,1.978D0,2.071D0,2.159D0, & 2.286D0,2.366D0,2.482D0,2.56D0, & 2.735D0, 2.90D0/ DATA UMOKC1/ 1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0,1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0,1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0,1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0/ DATA UMOK01/ 1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0,1.44D0, & 1.598D0,1.7D0,1.8D0,1.9D0,2.0D0,2.1D0,2.2D0,2.3D0,2.4D0,2.5D0, & 2.6D0,2.7D0,2.8D0,2.9D0,3.0D0, & 3.1D0/ * pp pn np nn * DATA UMOP1/ & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0, & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0, & 1.88D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0/ * app apn anp ann * DATA UMON1 / & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0, & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0, & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0, & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0, & 1.877D0,1.87701D0,1.879D0,1.887D0,1.9D0,1.917D0,1.938D0,1.962D0, & 2.D0,2.102D0,2.2D0,2.3D0,2.4D0,2.5D0,2.6D0,2.7D0,2.8D0,2.9D0, & 3.D0,3.1D0,3.2D0/ **** reaction channel state particles * DATA NRKPI1/ 13, 1, 15, 21, 81, 0, 13, 54, 23, 53, 13, 63, 13, 58, & 23, 57, 13, 65, 1, 32, 53, 31, 54, 32, 53, 33, 53, 35, 63, 32, & 13, 8, 23, 1, 17, 15, 21, 24, 22, 15, 82, 0, 61, 0, 13, 55, 23, & 54, 14, 53, 13, 64, 23, 63, 13, 59, 23, 58, 14, 57, 13, 66, 23, & 65, 1, 31, 8, 32, 1, 33, 1, 35, 54, 31, 55, 32, 54, 33, 53, 34, & 54, 35, 14, 1, 23, 8, 17, 24, 20, 15, 22, 24, 83, 0, 62, 0, 14, & 54, 23, 55, 13, 56, 14, 63, 23, 64, 14, 58, 23, 59, 13, 60, 14, & 65, 23, 66, 8, 31, 1, 34, 8, 33, 8, 35, 55, 31, 54, 34, 55, 33, & 56, 32, 55, 35, 14, 8, 24, 20, 84, 0, 14, 55, 23, 56, 14, 64, 14, & 59, 23, 60, 14, 66, 8, 34, 56, 31, 55, 34, 56, 33, 56, 35, 64,34/ DATA NRKKC1/15, 1, 89, 0, 24, 53, 15, 54, 1, 36, 1, 40, 1, 44, 36, & 63, 15, 63, 45, 53, 44, 54, 15, 8, 24, 1, 91, 0, 24, 54, 15, 55, & 8, 36, 1, 37, 8, 40, 1, 41, 8, 44, 1, 45, 36, 64, 37, 63, 15, 64, & 24, 63, 45, 54, 44, 55, 16, 1, 25, 8, 17, 23, 21, 14, 20, & 13, 22, 23, 90, 0, 38, 1, 39, 8, 16, 54, 25, 55, 1, 42, 8, 43, & 16, 63, 25, 64, 39, 64, 38, 63, 46, 54, 47, 55, 8, 47, 1, 46, 52, & 0, 51, 0, 16, 8, 17, 14, 20, 23, 22, 14, 92, 0, 8, 38, 16, 55, & 25, 56, 8, 42, 16, 64, 38, 64, 46, 55, 47, 56, 8, 46, 94, 0 / * * * k0 p k0 n ak0 p ak/ n * * * DATA NRKK01/ 24, 8, 106, 0, 15, 56, 24, 55, 37, 8, 41, 8, 45, 8, & 37, 64, 24, 64, 44, 56, 45, 55, 25, 1, 17, 13, 22, 13, 21, 23, & 107, 0, 39, 1, 25, 54, 16, 53, 43, 1, 25, 63, 39, 63, 47, 54, 46, & 53, 47, 1, 103, 0, 93, 0/ * pp pn np nn * DATA NRKP1/ 1, 1, 85, 0, 8, 53, 1, 54, 1, 63, 8, 57, 1, 58, 2*54, & 53, 55, 63, 54, 64, 53, 1, 8, 86, 0, 8, 54, 1, 55, 8, 63, 1, 64, & 8, 58, 1, 59, 64, 54, 63, 55, 54, 55, 53, 56, 77, 0, 2*8, 95, 0, & 8, 55, 1, 56, 8, 64, 8, 59, 1, 60, 2*55, 54, 56, 64, 55, 63, 56 / * app apn anp ann * DATA NRKN1/ 1, 2, 17, 18, 15, 16, 8, 9, 13, 14, 99, 0, 87, 0, 1, & 68, 8, 69, 2, 54, 9, 55, 102, 0, 2, 63, 9, 64, 1, 75, 8, 76, 53, & 67, 54, 68, 55, 69, 56, 70, 63, 68, 64, 69, 75, 54, 76, 55, 2, 8, & 18, 20, 16, 24, 14, 23, 101, 0, 88, 0, 2, 55, 9, 56, 1, 67, 8, & 68, 2, 64, 8, 75, 2, 59, 8, 72, 68, 55, 67, 54, 69, 56, 1, 9, 18, & 21, 15, 25, 13, 23, 100, 0, 96, 0, 2, 53, 9, 54, 1, 69, 8, 70, 1, & 76, 9, 63, 1, 73, 9, 58, 55, 70, 53, 68, 54, 69 / **** channel cross section * DATA SPKP11/ 0.D0, 300.D0, 40.D0, 20.D0, 13.D0,8.5D0,8.D0, 9.5D0, & 12.D0,14.D0,15.5D0,20.D0,17.D0,13.D0,10.D0,9.D0,8.5D0,8.D0,7.8D0, & 7.3D0, 6.7D0, 9*0.D0,.23D0,.35D0,.7D0,.52D0,.4D0,.3D0,.2D0,.15D0, & .13D0, .11D0, .09D0, .07D0, 0.D0, .033D0,.8D0,1.35D0,1.35D0,.5D0, & 15*0.D0, 3*0.D0,.00D0,0.80D0,2.2D0,3.6D0,4.6D0,4.7D0,3.5D0,2.4D0, &1.8D0,1.4D0,.75D0,.47D0,.25D0,.13D0,.08D0,6*0.D0,0.D0,1.2D0,3.3D0, & 5.4D0,6.9D0,7.3D0,5.3D0,3.6D0,2.7D0,2.2D0,1.1D0,.73D0,.4D0,.22D0, & .12D0,9*0.D0,.0D0,0.D0,2.0D0,4.4D0,6.8D0,9.9D0,7.9D0,6.0D0,3.8D0, &2.5D0,2.D0,1.4D0,1.D0,.6D0,.35D0,10*0.D0,.25D0,.55D0,.75D0,1.25D0, & 1.9D0,2.D0,1.8D0,1.5D0,1.25D0,1.D0,.8D0,6*0.D0,4*0.D0,.4D0,.85D0, & 1.1D0, 1.85D0, 2.8D0, 3.D0,2.7D0,2.2D0,1.85D0,1.5D0,1.2D0,6*0.D0, & 6*0.D0, .5D0, 1.2D0, 1.7D0, 3.4D0, 5.2D0, 6.4D0, 6.1D0, 5.6D0, & 5.2D0, 6*0.D0, 2*0.D0, .0D0, 1.D0, 3.3D0, 5.2D0, 4.45D0, 3.6D0, & 2.75D0, 1.9D0, 1.65D0, 1.3D0, .95D0, .6D0, .45D0, 6*0.D0, 3*0.D0, & .0D0, .45D0, 1.4D0, 1.5D0, 1.1D0, .85D0, .5D0, .3D0, .2D0, .15D0, & 8*0.D0, 5*0.D0, .0D0, .0D0, .6D0, .8D0, .95D0, .8D0, .7D0, .6D0, & .5D0, .4D0, 6*0.D0, 5*0.D0, .0D0, .00D0, .85D0, 1.2D0, 1.4D0, & 1.2D0, 1.05D0, .9D0, .7D0, .55D0, 6*0.D0, 5*0.D0, .0D0, .00D0, & 1.D0, 1.5D0, 3.5D0, 4.15D0, 3.7D0, 2.7D0, 2.3D0, 1.75D0, 6*0.D0, & 10*0.D0, .5D0, 2.0D0, 3.3D0, 5.4D0, 7.D0 / **** pi+ n data * DATA SPKPU1/ 0.D0, 25.D0, 13.D0, 11.D0, 10.5D0, 14.D0, 20.D0, & 20.D0, 16.D0, 14.D0, 19.D0, 28.D0, 17.5D0, 13.5D0, 12.D0, 10.5D0, & 10.D0, 10.D0, 9.5D0, 9.D0, 8.D0, 7.5D0, 7.D0, 6.5D0, 6.D0, 0.D0, & 48.D0, 19.D0, 15.D0, 11.5D0, 10.D0, 8.D0, 6.5D0, 5.5D0, 4.8D0, & 4.2D0, 7.5D0, 3.4D0, 2.5D0, 2.5D0, 2.1D0, 1.4D0, 1.D0, .8D0, & .6D0, .46D0, .3D0, .2D0, .15D0, .13D0, 11*0.D0, .95D0, .65D0, & .48D0, .35D0, .2D0, .18D0, .17D0, .16D0, .15D0, .1D0, .09D0, & .065D0, .05D0, .04D0, 12*0.D0, .2D0, .25D0, .25D0, .2D0, .1D0, & .08D0, .06D0, .045D0, .03D0, .02D0, .01D0, .005D0, .003D0, & 12*0.D0, .3D0, .24D0, .18D0, .15D0, .13D0, .12D0, .11D0, .1D0, & .09D0, .08D0, .05D0, .04D0, .03D0, 0.D0, 0.16D0, .7D0, 1.3D0, & 3.1D0, 4.5D0, 2.D0, 18*0.D0, 3*.0D0, 0.D0, 0.D0, 4.0D0, 11.D0, & 11.4D0, 10.3D0, 7.5D0, 6.8D0, 4.75D0, 2.5D0, 1.5D0, .9D0, .55D0, & .35D0, 13*0.D0, .1D0, .34D0, .5D0, .8D0, 1.1D0, 2.25D0, 3.3D0, & 2.3D0, 1.6D0, .95D0, .45D0, .28D0, .15D0, 10*0.D0, 2*0.D0, .17D0, & .64D0, 1.D0, 1.5D0, 2.1D0, 4.25D0, 6.2D0, 4.4D0, 3.D0, 1.8D0, & .9D0, .53D0, .28D0, 10*0.D0, 2*0.D0, .25D0, .82D0, & 1.3D0, 1.9D0, 2.8D0, 5.5D0 , 8.D0, 5.7D0, 3.9D0, 2.35D0, 1.15D0, & .69D0, .37D0, 10*0.D0, 7*0.D0, .0D0, .34D0, 1.5D0, 3.47D0, & 5.87D0, 6.23D0, 4.27D0, 2.6D0, 1.D0, .6D0, .3D0, .15D0, 6*0.D0/ * DATA CPKPV1/ 7*0.D0, .00D0, .16D0, .75D0, 1.73D0, 2.93D0, 3.12D0, & 2.13D0, 1.3D0, .5D0, .3D0, .15D0, .08D0, 6*0.D0, 10*0.D0, .2D0, & .6D0, .92D0, 2.4D0, 4.9D0, 6.25D0, 5.25D0, 3.5D0, 2.15D0, 1.4D0, & 1.D0, .7D0, 13*0.D0, .13D0, .4D0, .62D0, 1.6D0, 3.27D0, 4.17D0, & 3.5D0, 2.33D0, 1.43D0, .93D0, .66D0, .47D0, 13*0.D0, .07D0, .2D0, & .31D0, .8D0, 1.63D0, 2.08D0, 1.75D0, 1.17D0, .72D0, .47D0, .34D0, & .23D0, 17*0.D0, .33D0, 1.D0, 1.8D0, 2.67D0, 5.33D0, 6.D0, 5.53D0, & 5.D0, 17*0.D0, .17D0, .5D0, .9D0, 1.83D0, 2.67D0, 3.0D0, 2.77D0, & 2.5D0, 3*0.D0, 3*0.D0, 1.D0, 3.3D0, 2.8D0, 2.5D0, 2.3D0, 1.8D0, & 1.5D0, 1.1D0, .8D0, .7D0, .55D0, .3D0, 10*0.D0, 9*0.D0, .1D0, & .4D0, 1.D0, 1.4D0, 2.2D0, 2.5D0, 2.2D0, 1.65D0, 1.35D0, 1.1D0, & .8D0, .6D0, .4D0, 12*0.D0, .15D0, .6D0, 1.5D0, 2.1D0, 3.3D0, & 3.8D0, 3.3D0, 2.45D0, 2.05D0, 1.65D0, 1.2D0, .9D0, .6D0, 3*0.D0, & 9*0.D0, .10D0, .2D0, .5D0, .7D0, 1.3D0, 1.55D0, 1.9D0, 1.8D0, & 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0, 13*0.D0, .2D0, .5D0, .7D0, & 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0, & 17*0.D0, .2D0, .5D0, .85D0, 2.D0, 2.15D0, 2.05D0, 1.75D0, 1.D0, & 17*0.D0, .13D0, .33D0, .57D0, 1.33D0, 1.43D0, 1.36D0, 1.17D0, & .67D0, 17*0.D0, .07D0, .17D0, .28D0, .67D0, .72D0, .69D0, .58D0, & .33D0,17*0.D0,.4D0, .7D0, 1.D0, 1.6D0, 1.8D0, 2.3D0,1.9D0,1.7D0 / **** pi- p data * DATA SPKPW1/ 0.D0, 25.D0, 13.D0, 11.D0, 10.5D0, 14.D0, 2*20.D0, & 16.D0, 14.D0, 19.D0, 28.D0, 17.5D0, 13.5D0, 12.D0, 10.5D0, & 2*10.D0, 9.5D0, 9.D0, 8.D0, 7.5D0, 7.D0, 6.5D0, 6.D0, 0.D0, & 48.D0, 19.D0, 15.D0, 11.5D0, 10.D0, 8.D0, 6.5D0, 5.5D0, 4.8D0, & 4.2D0, 7.5D0, 3.4D0, 2*2.5D0, 2.1D0, 1.4D0, 1.D0, .8D0, .6D0, & .46D0, .3D0, .2D0, .15D0, .13D0, 11*0.D0, .95D0, .65D0, .48D0, & .35D0, .2D0, .18D0, .17D0, .16D0, .15D0, .1D0, .09D0, .065D0, & .05D0, .04D0, 12*0.D0, .2D0, 2*.25D0, .2D0, .1D0, .08D0, .06D0, & .045D0, .03D0, .02D0, .01D0, .005D0, .003D0, 12*0.D0, .3D0, & .24D0, .18D0, .15D0, .13D0, .12D0, .11D0, .1D0, .09D0, .08D0, & .05D0, .04D0, .03D0, 0.D0, 0.16D0, .7D0, 1.3D0, 3.1D0, 4.5D0, & 2.D0, 23*0.D0, 4.0D0, 11.D0, 11.4D0, 10.3D0, 7.5D0, 6.8D0, & 4.75D0, 2.5D0, 1.5D0, .9D0, .55D0, .35D0, 13*0.D0, .1D0, .34D0, & .5D0, .8D0, 1.1D0, 2.25D0, 3.3D0, 2.3D0, 1.6D0, .95D0, .45D0, & .28D0, .15D0, 12*0.D0, .17D0, .64D0, 1.D0, 1.5D0, 2.1D0, 4.25D0, & 6.2D0, 4.4D0, 3.D0, 1.8D0, .9D0, .53D0, .28D0, 12*0.D0, .25D0, & .82D0, 1.3D0, 1.9D0, 2.8D0, 5.5D0, 8.D0, 5.7D0, 3.9D0, 2.35D0, & 1.15D0, .69D0, .37D0, 18*0.D0, .34D0, 1.5D0, 3.47D0, 5.87D0, & 6.23D0, 4.27D0, 2.6D0, 1.D0, .6D0, .3D0, .15D0, 6*0.D0/ * DATA SPKPX1/ 8*0.D0, .16D0, .75D0, 1.73D0, 2.93D0, 3.12D0, & 2.13D0, 1.3D0, .5D0, .3D0, .15D0, .08D0, 16*0.D0, .2D0, .6D0, & .92D0, 2.4D0, 4.9D0, 6.25D0, 5.25D0, 3.5D0, 2.15D0, 1.4D0, 1.D0, & .7D0, 13*0.D0, .13D0, .4D0, .62D0, 1.6D0, 3.27D0, 4.17D0, 3.5D0, & 2.33D0, 1.43D0, .93D0, .66D0, .47D0, 13*0.D0, .07D0, .2D0, .31D0, & .8D0, 1.63D0, 2.08D0, 1.75D0, 1.17D0, .72D0, .47D0, .34D0, .23D0, & 17*0.D0, .33D0, 1.D0, 1.8D0, 2.67D0, 5.33D0, 6.D0, 5.53D0, 5.D0, & 17*0.D0, .17D0, .5D0, .9D0, 1.83D0, 2.67D0, 3.0D0, 2.77D0, 2.5D0, & 6*0.D0, 1.D0, 3.3D0, 2.8D0, 2.5D0, 2.3D0, 1.8D0, 1.5D0, 1.1D0, & .8D0, .7D0, .55D0, .3D0, 19*0.D0, .1D0, .4D0, 1.D0, 1.4D0, 2.2D0, & 2.5D0, 2.2D0, 1.65D0, 1.35D0, 1.1D0, .8D0, .6D0, .4D0, 12*0.D0, & .15D0, .6D0, 1.5D0, 2.1D0, 3.3D0, 3.8D0, 3.3D0, 2.45D0, 2.05D0, & 1.65D0, 1.2D0, .9D0, .6D0, 12*0.D0, .10D0, .2D0, .5D0, .7D0, & 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0, 1.35D0, 1.15D0, .95D0, .7D0, & 13*0.D0, .2D0, .5D0, .7D0, 1.3D0, 1.55D0, 1.9D0, 1.8D0, 1.55D0, & 1.35D0, 1.15D0, .95D0, .7D0, 17*0.D0, .2D0, .5D0, .85D0, 2.D0, & 2.15D0, 2.05D0, 1.75D0, 1.D0, 17*0.D0, .13D0, .33D0, .57D0, & 1.33D0, 1.43D0, 1.36D0, 1.17D0, .67D0, 17*0.D0, .07D0, .17D0, & .28D0, .67D0, .72D0, .69D0, .58D0, .33D0, 17*0.D0, .4D0, .7D0, & 1.D0, 1.6D0, 1.8D0, 2.3D0, 1.9D0, 1.7D0 / **** pi- n data * DATA SPKP41 / 0.D0, 300.D0, 40.D0, 20.D0, 13.D0, 8.5D0, 8.D0, & 9.5D0, 12.D0, 14.D0, 15.5D0, 20.D0, 17.D0, 13.D0, 10.D0, 9.D0, & 8.5D0, 8.D0, 7.8D0, 7.3D0, 6.7D0, 9*0.D0, .23D0, .35D0, .7D0, & .52D0, .4D0, .3D0, .2D0, .15D0, .13D0, .11D0, .09D0, .07D0, 0.D0, & .033D0, .8D0, 2*1.35D0, .5D0, 19*0.D0, 0.8D0, 2.2D0, 3.6D0, & 4.6D0, 4.7D0, 3.5D0, 2.4D0, 1.8D0, 1.4D0, .75D0, .47D0, .25D0, & .13D0, .08D0, 7*0.D0, 1.2D0, 3.3D0, 5.4D0, 6.9D0, 7.3D0, 5.3D0, & 3.6D0, 2.7D0, 2.2D0, 1.1D0, .73D0, .4D0, .22D0, .12D0, 11*0.D0, & 2.0D0, 4.4D0, 6.8D0, 9.9D0, 7.9D0, 6.0D0, 3.8D0, 2.5D0, 2.D0, & 1.4D0, 1.D0, .6D0, .35D0, 10*0.D0, .25D0, .55D0, .75D0, 1.25D0, & 1.9D0, 2.D0, 1.8D0, 1.5D0, 1.25D0, 1.D0, .8D0, 10*0.D0, .4D0, & .85D0, 1.1D0, 1.85D0, 2.8D0, 3.D0, 2.7D0, 2.2D0, 1.85D0, 1.5D0, & 1.2D0, 12*0.D0, .5D0, 1.2D0, 1.7D0, 3.4D0, 5.2D0, 6.4D0, 6.1D0, & 5.6D0, 5.2D0, 9*0.D0, 1.D0, 3.3D0, 5.2D0, 4.45D0, 3.6D0, 2.75D0, & 1.9D0, 1.65D0, 1.3D0, .95D0, .6D0, .45D0, 10*0.D0, .45D0, 1.4D0, & 1.5D0, 1.1D0, .85D0, .5D0, .3D0, .2D0, .15D0, 15*0.D0, .6D0, & .8D0, .95D0, .8D0, .7D0, .6D0, .5D0, .4D0, 13*0.D0, .85D0, 1.2D0, & 1.4D0, 1.2D0, 1.05D0, .9D0, .7D0, .55D0, 13*0.D0, 1.D0, 1.5D0, & 3.5D0, 4.15D0, 3.7D0, 2.7D0, 2.3D0, 1.75D0, 16*0.D0, .5D0, 2.0D0, & 3.3D0, 5.4D0, 7.D0 / **** k+ p data * DATA SPKP51/ 0.D0, 20.D0, 14.D0, 12.D0, 11.5D0, 10.D0, 8.D0, & 7.D0, 6.D0, 5.5D0, 5.3D0, 5.D0, 4.5D0, 4.4D0, 3.8D0, 3.D0, 2.8D0, & 0.D0, .5D0, 1.15D0, 2.D0, 1.3D0, .8D0, .45D0, 13*0.D0, 0.9D0, & 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0, 1.5D0, 1.2D0, .9D0, .6D0, & .45D0, .21D0, .2D0, 3*0.D0, .9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, & 2.D0, 1.7D0, 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0, & 4*0.D0, 1.D0, 2.1D0, 2.6D0, 2.3D0, 2.1D0, 1.8D0, 1.7D0, 1.4D0, & 1.2D0, 1.05D0, .9D0, .66D0, .5D0, 7*0.D0, .3D0, 2*1.D0, .9D0, & .7D0, .4D0, .3D0, .2D0, 11*0.D0, .1D0, 1.D0, 2.2D0, 3.5D0, 4.2D0, & 4.55D0, 4.85D0, 4.9D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, & 1.71D0, 1.6D0, 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0, 4.4D0, 4.D0, & 3.5D0, 2.85D0, 2.35D0, 2.01D0, 1.8D0, 12*0.D0, .1D0, .8D0,2.05D0, & 3.31D0, 3.5D0, 12*0.D0, .034D0, .2D0, .75D0, 1.04D0, 1.24D0 / **** k+ n data * DATA SPKP61/ 0.D0, 6.D0, 11.D0, 13.D0, 6.D0, 5.D0, 3.D0, 2.2D0, & 1.5D0, 1.2D0, 1.D0, .7D0, .6D0, .5D0, .45D0, .35D0, .3D0, 0.D0, & 6.D0, 11.D0, 13.D0, 6.D0, 5.D0, 3.D0, 2.2D0, 1.5D0, 1.2D0, 1.D0, & .7D0, .6D0, .5D0, .45D0, .35D0, .3D0, 0.D0, .5D0, 1.3D0, 2.8D0, & 2.3D0, 1.6D0, .9D0, 13*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, & 2.D0, 1.7D0, 1.5D0,1.2D0,.9D0,.6D0,.45D0,.21D0,.2D0,3*0.D0,0.9D0, & 2.5D0, 3.D0, 2.5D0, 2.3D0,2.D0,1.7D0,1.5D0,1.2D0,.9D0,.6D0,.45D0, & .21D0, .2D0,4*0.D0,1.D0,2.1D0,2.6D0,2.3D0,2.D0,1.8D0,1.7D0,1.4D0, & 1.2D0,1.15D0,.9D0,.66D0,.5D0,4*0.D0,1.D0,2.1D0,2.6D0,2.3D0,2.1D0, & 1.8D0,1.7D0,1.4D0,1.2D0, 1.15D0, .9D0, .66D0, .5D0, 7*0.D0, .3D0, & 2*1.D0, .9D0, .7D0, .4D0, .35D0, .2D0, 9*0.D0, .3D0, 2*1.D0,.9D0, & .7D0, .4D0, .35D0, .2D0, 11*0.D0, .1D0, 1.D0, 2.4D0,3.5D0,4.25D0, & 4.55D0, 4.85D0, 4.9D0, 9*0.D0, .1D0, 1.D0, 2.4D0, 3.5D0, 4.25D0, & 4.55D0, 4.85D0, 4.9D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, & 1.71D0, 1.6D0, 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, 1.71D0, & 1.6D0, 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0,4.4D0,4.D0,3.5D0,2.85D0, & 2.35D0, 2.01D0, 1.8D0, 6*0.D0, 1.4D0,3.8D0,5.D0,4.7D0,4.4D0,4.D0, & 3.5D0,2.85D0,2.35D0,2.01D0,1.8D0,12*0.D0,.1D0,.8D0,2.05D0,3.31D0, & 3.5D0, 12*0.D0, .034D0,.2D0,.75D0,1.04D0,1.24D0 / **** k- p data * DATA SKMPL1/ 0.D0, 35.D0, 22.D0, 25.D0, 17.D0, 9.D0, 9.5D0, 8.D0, & 7.D0, 6.5D0, 6.1D0, 5.D0, 4.8D0, 4.6D0, 4.45D0, 4.3D0, 4.2D0, & 0.D0, 8.D0, 3.5D0, 8.D0, 3.D0, 1.9D0, 1.7D0, 1.D0, .9D0, .8D0, & .75D0, .5D0, .42D0, .38D0, .34D0, .25D0, .2D0, & 0.D0, 3.D0, 3.2D0, 3.5D0, 1.5D0, 1.4D0, 1.1D0, .6D0, .5D0, & .35D0, .28D0, .25D0, .18D0, .12D0, .1D0, .08D0, .04D0, & 0.D0, 8.5D0, 2.4D0, 1.7D0, 1.3D0, 1.3D0, 1.1D0, .5D0, & .4D0, .4D0, .35D0, .3D0, .28D0, .2D0, .16D0, .13D0, .11D0, & 0.D0, 7.D0, 4.8D0, 1.4D0, 1.9D0, .9D0, .4D0, .2D0, .13D0, & .1D0, .08D0, .06D0, .04D0, .02D0, .015D0, .01D0, .01D0, & 0.D0, 5.5D0, 1.D0, .8D0, .75D0, .32D0, .2D0, .1D0, .09D0, & .08D0, .065D0, .05D0, .04D0, .022D0, .017D0, 2*.01D0/ DATA SPKP71 / 0.D0, .56D0, 1.46D0, 3.16D0, 2.01D0, 1.28D0, .74D0, & 14*0.D0, 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0, & 1.91D0, 1.57D0, 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0, 4*0.D0, & 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0, 1.91D0, 1.57D0, & 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0, 3*0.D0, 1.0D0, 3.03D0, & 3.36D0, 2.8D0, 2.58D0, 2.24D0, 1.91D0, 1.68D0, 1.35D0, 1.01D0, & .67D0, .5D0, .24D0, .23D0, 3*0.D0, 1.0D0, 3.03D0, 3.36D0, 2.8D0, & 2.58D0, 2.24D0, 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0, & .24D0, .23D0, 7*0.D0, .34D0, 1.12D0, 1.12D0, 1.01D0, .78D0, & .45D0, .39D0, .22D0, .07D0, 0.D0, 7*0.D0, .34D0, 1.12D0, 1.12D0, & 1.01D0, .78D0, .45D0, .39D0, .22D0, .07D0, 0.D0, 6*0.D0, 1.71D0, & 4.26D0, 5.6D0, 5.57D0, 4.93D0, 4.48D0, 3.92D0, 3.19D0, 2.63D0, & 2.25D0, 2.D0, 6*0.D0, 1.71D0, 4.26D0, 5.6D0, 5.57D0, 4.93D0, & 4.48D0, 3.92D0, 3.19D0, 2.63D0, 2.25D0, 2.D0, 10*0.D0, .22D0, & .8D0, .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0, 10*0.D0, .22D0, .8D0, & .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0, 13*0.D0, .1D0, .3D0, .7D0,1.D0, & 13*0.D0, .1D0, .3D0, .7D0, 1.D0, 9*0.D0, .11D0, 1.72D0, 2.69D0, & 3.92D0, 4.76D0, 5.10D0, 5.44D0, 5.3D0, 9*0.D0, .11D0, 1.72D0, & 2.69D0, 3.92D0, 4.76D0, 5.1D0, 5.44D0, 5.3D0, 5*0.D0,9.2D0,4.7D0, & 1.9D0, 10*0.D0, 2.5D0, 15.D0, 21.5D0, 15.3D0, 3.D0,1.5D0,10*0.D0/ ***** k- n data * DATA SKMNL1/0.D0, 4.D0, 9.5D0, 20.D0, 13.D0, 9.5D0, 6.D0, 4.4D0, & 3.D0, 2.4D0, 2.D0, 1.4D0, 1.2D0, 1.D0, .9D0, .7D0, .6D0, & 0.D0, 4.5D0, 6.D0, 5.D0, 2.5D0, 2.D0, 1.7D0, 2.1D0, & 1.9D0, .9D0, .5D0, .3D0, .24D0, .2D0, .18D0, .1D0, .09D0, & 0.D0, 1.8D0, 2.D0, 1.1D0, .9D0, .5D0, .5D0, .4D0, .4D0, & .2D0, .1D0, .06D0, .05D0, .04D0, .03D0, .02D0, .02D0, & 0.D0, 1.5D0, 2.D0, .9D0, 1.1D0, .4D0, .6D0, .7D0, .65D0, & .3D0, .17D0, .1D0, .08D0, .07D0, .06D0, .04D0, .03D0/ DATA SPKP81/0.D0, .56D0, 1.29D0, 2.26D0, 1.01D0, .64D0, .37D0, & 14*0.D0, 1.13D0, 2.61D0, 2.91D0, 2.58D0, 2.35D0, 2.02D0, & 1.91D0, 1.57D0, 1.35D0, 1.29D0, 1.01D0, .74D0, .65D0, & 3*0.D0, 1.D0, 3.03D0, 3.36D0, 2.8D0, 2.58D0, 2.24D0, & 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0, .24D0, .23D0, & 3*0.D0, 1.D0, 3.03D0, 3.36D0, 2.8D0, 2.58D0, 2.24D0, & 1.91D0, 1.68D0, 1.35D0, 1.01D0, .67D0, .5D0, .24D0, .23D0, & 7*0.D0, .34D0, 1.12D0, 1.12D0, 1.01D0, .78D0, .45D0, & .39D0, .22D0, .07D0, 0.D0, & 6*0.D0, 1.71D0, 4.26D0, 5.6D0, 5.57D0, 4.93D0, & 4.48D0, 3.92D0, 3.19D0, 2.63D0, 2.25D0, 2.D0, & 10*0.D0, .22D0, .8D0, .75D0, 1.D0, 1.3D0, 1.5D0, 1.3D0, & 13*0.D0, .1D0, .3D0, .7D0, 1.D0, & 13*0.D0, .1D0, .3D0, .7D0, 1.D0, & 9*0.D0, .11D0, 1.72D0, 2.69D0, 3.92D0, 4.76D0, & 5.10D0, 5.44D0, 5.3D0, & 4*0.D0, 0.00D0, 9.2D0, 4.7D0, 1.9D0, 9*0.D0/ ***** p p data * DATA SPKP91/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0, & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0, & 0.D0, 3.6D0, 1.7D0, 10*0.D0, & .0D0, 0.D0, 8.7D0, 17.7D0, 18.8D0, 15.9D0, & 11.7D0, 8.D0, 6.D0, 5.3D0, 4.5D0, 3.9D0, 3.5D0, & .0D0, .0D0, 2.8D0, 5.8D0, 6.2D0, 5.1D0, 3.8D0, & 2.7D0, 2.1D0, 1.8D0, 1.5D0, 1.3D0, 1.1D0, & 5*0.D0, 4.6D0, 10.2D0, 15.1D0, & 16.9D0, 16.5D0, 11.D0, 5.5D0, 3.5D0, & 10*0.D0, 4.3D0, 7.6D0, 9.D0, & 10*0.D0, 1.7D0, 2.6D0, 3.D0, & 6*0.D0, .3D0, .6D0, 1.D0, 1.6D0, 1.3D0, .8D0, .6D0, & 6*0.D0, .7D0, 1.2D0, 1.8D0, 2.5D0, 1.8D0, 1.3D0, & 1.2D0, 10*0.D0, .6D0, 1.4D0, 1.7D0, & 10*0.D0, 1.9D0, 4.1D0, 5.2D0/ ***** p n data * DATA SPKP01/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0, & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0, & 0.D0, 1.8D0, .2D0, 12*0.D0, & 3.2D0, 6.05D0, 9.9D0, 5.1D0, & 3.8D0, 2.7D0, 1.9D0, 1.5D0, 1.4D0, 1.3D0, 1.1D0, & 2*.0D0, 3.2D0, 6.05D0, 9.9D0, 5.1D0, & 3.8D0, 2.7D0, 1.9D0, 1.5D0, 1.4D0, 1.3D0, 1.1D0, & 5*0.D0, 4.6D0, 10.2D0, 15.1D0, & 16.4D0, 15.2D0, 11.D0, 5.4D0, 3.5D0, & 5*0.D0, 4.6D0, 10.2D0, 15.1D0, & 16.4D0, 15.2D0, 11.D0, 5.4D0, 3.5D0, & 10*0.D0, .7D0, 5.1D0, 8.D0, & 10*0.D0, .7D0, 5.1D0, 8.D0, & 10*.0D0, .3D0, 2.8D0, 4.7D0, & 10*.0D0, .3D0, 2.8D0, 4.7D0, & 7*0.D0, 1.2D0, 2.5D0, 3.5D0, 6.D0, 5.3D0, 2.9D0, & 7*0.D0, 1.7D0, 3.6D0, 5.4D0, 9.D0, 7.6D0, 4.2D0, & 5*0.D0, 7.7D0, 6.1D0, 2.9D0, 5*0.D0/ * nn - data * * * DATA SPKPV1/ 0.D0, 24.D0, 25.D0, 27.D0, 23.D0, 21.D0, 20.D0, & 19.D0, 17.D0, 15.5D0, 14.D0, 13.5D0, 13.D0, & 0.D0, 3.6D0, 1.7D0, 12*0.D0, & 8.7D0, 17.7D0, 18.8D0, 15.9D0, & 11.7D0, 8.D0, 6.D0, 5.3D0, 4.5D0, 3.9D0, 3.5D0, & .0D0, .0D0, 2.8D0, 5.8D0, 6.2D0, 5.1D0, 3.8D0, & 2.7D0, 2.1D0, 1.8D0, 1.5D0, 1.3D0, 1.1D0, & 5*0.D0, 4.6D0, 10.2D0, 15.1D0, 16.9D0, 16.5D0, & 11.D0, 5.5D0, 3.5D0, & 10*0.D0, 4.3D0, 7.6D0, 9.D0, & 10*0.D0, 1.7D0, 2.6D0, 3.D0, & 6*0.D0, .3D0, .6D0, 1.D0, 1.6D0, 1.3D0, .8D0, .6D0, & 6*0.D0, .7D0, 1.2D0, 1.8D0, 2.5D0, 1.8D0, 1.3D0, & 1.2D0, 10*0.D0, .6D0, 1.4D0, 1.7D0, & 10*0.D0, 1.9D0, 4.1D0, 5.2D0/ **************** ap - p - data * DATA SAPPL1/ 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0, & 50.D0, 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0, & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0, & .05D0, .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0, .85D0, & 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0, .14D0, & .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0, & 0.D0, 55.D0, 50.D0, 25.D0, 15.D0, 15.D0, 14.D0, 12.D0, & 10.D0, 7.D0, 6.D0, 4.D0, 3.3D0, 2.8D0, 2.4D0, 2.D0, 1.8D0, & 1.55D0, 1.3D0, .95D0, .75D0, & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0, & .25D0, .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0, & .01D0, .008D0, .006D0, .005D0/ DATA SPKPE1/0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0, & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0, & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0, & 9*0.D0, 2.D0, 2.5D0, .2D0, 19*0.D0, .3D0, 1.4D0, 2.2D0, 1.2D0, & 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0, .3D0, 10*0.D0, .3D0, 1.4D0, & 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0, .3D0, 10*0.D0, & .3D0, 1.4D0, 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0, .6D0, .5D0, .4D0, & .3D0, 10*0.D0, .3D0, 1.4D0, 2.2D0, 1.2D0, 1.1D0, 1.D0, .8D0, & .6D0, .5D0, .4D0, .3D0, 9*0.D0, .6D0, 2.5D0, 5.D0, 5.2D0, 5.1D0, & 5.4D0, 5.8D0, 2.8D0, 2.1D0, 1.8D0, 1.6D0, 1.2D0, 13*0.D0, 1.3D0, & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0, & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0, & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 13*0.D0, 1.3D0, & 1.5D0, 2.D0, 2.5D0, 2.5D0, 2.3D0, 1.8D0, 1.4D0, 14*0.D0, .2D0, & .5D0, 1.1D0, 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0, & 1.1D0, 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0, 1.1D0, & 1.6D0, 1.4D0, 1.1D0, .9D0, 14*0.D0, .2D0, .5D0, 1.1D0, 1.6D0, & 1.4D0, 1.1D0, .9D0, 17*0.D0, .3D0, 1.6D0, 2.6D0, 3.6D0, 17*0.D0, & .3D0, 1.6D0, 2.6D0, 3.6D0, 17*0.D0, .3D0, 1.6D0, 2.6D0, & 3.6D0, 17*0.D0, .3D0, 1.6D0, 2.6D0, 3.6D0 / **************** ap - n - data * DATA SAPNL1/ & 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0, & 50.D0, 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0, & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0, & .05D0, .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0, & .85D0, 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0, & .14D0, .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0, & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0, & .25D0, .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0, & .01D0, .008D0, .006D0, .005D0 / DATA SPKPZ1/ 0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0, & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0, & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0, 9*0.D0, 2.4D0, .2D0, & 20*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, & .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, & 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, & 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, & 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, & .7D0, .5D0, .3D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0, & 7.6D0, 6.D0, 5.D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0, & 7.6D0, 6.D0, 5.D0, 18*0.D0, 1.D0, 4.9D0, 8.5D0, 18*0.D0, 1.D0, & 4.9D0, 8.5D0, 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, & 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, 15*0.D0, 1.9D0, & 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0 / * * * * **************** an - p - data * * * DATA SANPL1/ & 0.D0, 176.D0, 160.D0, 105.D0, 75.D0, 68.D0, 65.D0, 50.D0, & 50.D0, 43.D0, 42.D0, 40.5D0, 35.D0, 30.D0, 28.D0, & 25.D0, 22.D0, 21.D0, 20.D0, 18.D0, 17.D0, 11*0.D0, .05D0, & .15D0, .18D0, .2D0, .2D0, .3D0, .4D0, .6D0, .7D0, .85D0, & 0.D0, 1.D0, .9D0, .46D0, .3D0, .23D0, .18D0, .16D0, .14D0, & .1D0, .08D0, .05D0, .02D0, .015D0, 4*.011D0, 3*.005D0, & 0.D0, 3.3D0, 3.D0, 1.5D0, 1.D0, .7D0, .4D0, .35D0, .4D0, .25D0, & .18D0, .08D0, .04D0, .03D0, .023D0, .016D0, .014D0, & .01D0, .008D0, .006D0, .005D0 / DATA SPKPF1/ 0.D0, 215.D0, 193.D0, 170.D0, 148.D0, 113.D0, 97.D0, & 84.D0, 78.D0, 68.D0, 64.D0, 61.D0, 46.D0, 36.D0, 31.3D0, 28.5D0, & 25.7D0, 22.6D0, 21.4D0, 20.7D0, 19.9D0, 9*0.D0, 2.4D0, .2D0, & 20*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, & .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, & 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, 10*0.D0, 1.8D0, 2.8D0, & 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, .7D0, .5D0, .3D0, & 10*0.D0, 1.8D0, 2.8D0, 3.6D0, 2.3D0, 1.8D0, 1.5D0, 1.3D0, 1.D0, & .7D0, .5D0, .3D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0, & 7.6D0, 6.D0, 5.D0, 13*0.D0, 5.2D0, 8.7D0, 11.4D0, 14.D0, 11.9D0, & 7.6D0, 6.D0, 5.D0, 18*0.D0, 1.D0, 4.9D0, 8.5D0, 18*0.D0, 1.D0, & 4.9D0, 8.5D0, 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, & 15*0.D0, 1.9D0, 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0, 15*0.D0, 1.9D0, & 2.3D0, 4.D0, 6.5D0, 5.2D0, 3.4D0 / **** ko - n - data * DATA CPKP51/0.D0, 20.D0, 14.D0, 12.D0, 11.5D0, 10.D0, 8.D0, 7.D0, & 6.D0, 5.5D0, 5.3D0, 5.D0, 4.5D0, 4.4D0, 3.8D0, 3.D0, 2.8D0, & 0.D0, .5D0, 1.15D0, 2.D0, 1.3D0, .8D0, .45D0, 10*0.D0, & 3*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0, & 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0, & 3*0.D0, 0.9D0, 2.5D0, 3.D0, 2.5D0, 2.3D0, 2.D0, 1.7D0, & 1.5D0, 1.2D0, .9D0, .6D0, .45D0, .21D0, .2D0, & 4*0.D0, 1.D0, 2.1D0, 2.6D0, 2.3D0, 2.1D0, 1.8D0, 1.7D0, & 1.4D0, 1.2D0, 1.05D0, .9D0, .66D0, .5D0, & 7*0.D0, .3D0, 1.D0, 1.D0, .9D0, .7D0, .4D0, .30D0, .2D0, & 11*0.D0, .1D0, 1.D0, 2.2D0, 3.5D0, 4.20D0, 4.55D0, & 4.85D0, 4.9D0, & 10*0.D0, .2D0, .7D0, 1.6D0, 2.5D0, 2.2D0, 1.71D0, 1.6D0, & 6*0.D0, 1.4D0, 3.8D0, 5.D0, 4.7D0, 4.4D0, 4.D0, 3.5D0, & 2.85D0, 2.35D0, 2.01D0, 1.8D0, & 12*0.D0, .1D0, .8D0, 2.05D0, 3.31D0, 3.5D0, & 12*0.D0, .034D0, .20D0, .75D0, 1.04D0, 1.24D0 / **** ako - p - data * DATA CPKP61/ 0.D0,4.D0,9.5D0,20.D0,13.D0,9.5D0,6.D0,4.4D0, & 3.D0,2.4D0,2.D0,1.4D0,1.2D0,1.D0,.9D0,.7D0,.6D0,0.D0, & 4.5D0,6.D0,5.D0,2.5D0,2.D0,1.7D0,2.1D0,1.9D0,.9D0,.5D0, & .3D0,.24D0,.2D0,.18D0,.1D0,.09D0,0.D0,1.8D0,2.D0,1.1D0, & .9D0,.5D0,.5D0,.4D0,.4D0,.2D0,.1D0,.06D0,.05D0,.04D0, & .03D0,.02D0,.02D0,0.D0,1.5D0,2.D0,.9D0,1.1D0,.4D0,.6D0, & .7D0,.65D0,.3D0,.17D0,.1D0,.08D0,.07D0,.06D0,.04D0, & .03D0,0.D0,.56D0,1.29D0,2.26D0,1.01D0,.64D0,.37D0, & 14*0.D0,1.13D0,2.61D0,2.91D0,2.58D0,2.35D0,2.02D0,1.91D0, & 1.57D0,1.35D0,1.29D0,1.01D0,.74D0,.65D0,3*0.D0,1.0D0, & 3.03D0,3.36D0,2.8D0,2.58D0,2.24D0,1.91D0,1.68D0,1.35D0, & 1.01D0,.67D0,.5D0,.24D0,.23D0,3*0.D0,1.0D0,3.03D0,3.36D0, & 2.8D0,2.58D0,2.24D0,1.91D0,1.68D0,1.35D0,1.01D0,.67D0, & .5D0,.24D0,.23D0,7*0.D0,.34D0,1.12D0,1.12D0,1.01D0,.78D0, & .45D0,.39D0,.22D0,.07D0,7*0D0,1.71D0,4.26D0,5.6D0,5.57D0,4.93D0, & 4.48D0,3.92D0,3.19D0,2.63D0,2.25D0,2.D0,10*0.D0,.22D0,.8D0, & .75D0,1.D0,1.3D0,1.5D0,1.3D0,13*0.D0,.1D0,.3D0,.7D0,1.D0, & 13*0.D0,.1D0,.3D0,.7D0,1.D0,9*0.D0,.11D0,1.72D0,2.69D0,3.92D0, & 4.76D0,5.10D0,5.44D0,5.3D0,5*0.D0,9.2D0,4.7D0,1.9D0,9*0.D0, & .0D0,2.5D0,15.D0,21.5D0,15.3D0,3.D0,1.5D0,10*0.D0 / DATA NURE1/9, 12, 5*0, 10, 14, 3*0, 1, 3, 5, 7, 6*0, 2, 6, 16, & 5*0, 10, 13, 5*0, 11, 12, 3*0, 2, 4, 6, 8, 6*0, 3, 15, 7, 5*0 / * DO 10 J=1,92 UMOPI(J)=UMOPI1(J) PLAPI(J)=PLAPI1(J) 10 CONTINUE DO 20 J=1,68 UMOKC(J)=UMOKC1(J) PLAKC(J)=PLAKC1(J) 20 CONTINUE DO 30 J=1,39 UMOP(J)=UMOP1(J) PLAP(J)=PLAP1(J) 30 CONTINUE DO 40 J=1,63 UMON(J)=UMON1(J) PLAN(J)=PLAN1(J) 40 CONTINUE DO 50 J=1,34 UMOK0(J)=UMOK01(J) PLAK0(J)=PLAK01(J) 50 CONTINUE DO 60 J=1,296 SIIN(J)=SIIN1(J) 60 CONTINUE DO 70 J=1,315 SPIKP1(J)=SPKP11(J) SPIKP4(J)=SPKP41(J) 70 CONTINUE DO 80 J=1,278 SPIKPU(J)=SPKPU1(J) SPIKPW(J)=SPKPW1(J) 80 CONTINUE DO 90 J=1,372 SPIKPV(J)=CPKPV1(J) SPIKPX(J)=SPKPX1(J) 90 CONTINUE DO 100 J=1,187 SPIKP5(J)=SPKP51(J) SPIKP8(J)=SPKP81(J) SPKP15(J)=CPKP51(J) 100 CONTINUE DO 110 J=1,289 SPIKP6(J)=SPKP61(J) 110 CONTINUE DO 120 J=1,102 SKMPEL(J)=SKMPL1(J) 120 CONTINUE DO 130 J=1,289 SPIKP7(J)=SPKP71(J) 130 CONTINUE DO 140 J=1,68 SKMNEL(J)=SKMNL1(J) 140 CONTINUE DO 150 J=1,143 SPIKP9(J)=SPKP91(J) SPKPV(J)=SPKPV1(J) 150 CONTINUE DO 160 J=1,169 SPIKP0(J)=SPKP01(J) 160 CONTINUE DO 170 J=1,105 SAPPEL(J)=SAPPL1(J) 170 CONTINUE DO 180 J=1,399 SPIKPE(J)=SPKPE1(J) 180 CONTINUE DO 190 J=1,84 SAPNEL(J)=SAPNL1(J) SANPEL(J)=SANPL1(J) 190 CONTINUE DO 200 J=1,273 SPIKPZ(J)=SPKPZ1(J) SPIKPF(J)=SPKPF1(J) 200 CONTINUE DO 210 J=1,272 SPKP16(J)=CPKP61(J) 210 CONTINUE DO 220 J=1,164 NRKPI(J)=NRKPI1(J) 220 CONTINUE DO 230 J=1,132 NRKKC(J)=NRKKC1(J) 230 CONTINUE DO 240 J=1,70 NRKP(J)=NRKP1(J) 240 CONTINUE DO 250 J=1,116 NRKN(J)=NRKN1(J) 250 CONTINUE DO 270 J=1,54 NRKK0(J)=NRKK01(J) 270 CONTINUE DO 280 J=1,60 NURELN(J)=NURE1(J) 280 CONTINUE *= end*block.blkdt3 * END