      subroutine begin
*
*  initialize variables for central boundary condition
*
      implicit double precision(a-h,o-z)

      common/shells/ sa(400),ra(400),ba(400),pa(400),ta(400),
     1 ea(400),xca(400),fca(400),s(400),r(400),b(400),
     2 p(400),t(400),xc(400),e(400),sk(400),
     3 rk(400),bk(400),pk(400),tk(400)
      common/contrl/ds,g,sm,wc,it,nite,ja,jb,j,k,l
      common/temp/s1,r1,s2,r2,b2,p2,t2,e2,xc2,xo2,fca2,f2,q2,w2,c
      common/fin/ifin

      j=2 
      k=0 
      l=2 
      s1=-10.
      r1=-10.
      s2=s(1)
      r2=r(1)
      b2=b(1)
      p2=p(1)
      t2=t(1)
      ifin = 0
 
      return
      end 

************************************************************************

      block data
*
*  block data:  initiallizes some physical constants, and also the
*  constants used by the conductive opacity regimes.  the numbers
*  have been taken from itoh et al (1983; 1984). 
*  Also added Itoh and company Bremsstrahlung neutrino coefficients.
*  See Itoh etal 1989, ApJ, 339, 354. + erratum 360, 741.
*
      implicit double precision(a-h,l,o-z)

      common/c/az(3),z(3),taux,ts,ps,sm,rs,err,te,gs,rho,ol1x,rl1x,el1x
     1,atgx,rtgx,tl,pl,xmass,acu,am,bk,an0,pi
*
* the following are the fitting coefficients quoted by itoh 
* "cdats" ==> carbon, solid state
* "cdatl" ==> carbon, liquid metal state
* "odats" ==> oxygen, solid state
* "odatl" ==> oxygen, liquid metal state
*  and so on.
*
      common /hydatl/ alhy(4),blhy(3),clhy(3),dlhy(4),elhy(3),flhy(3) 
 
      common /hedatl/ alhe(4),blhe(3),clhe(3),dlhe(4),elhe(3),flhe(3) 
      common/hedats/ sahe(5),sbhe(3),sehe(5),sfhe(3),sihe(5),sjhe(3), 
     1          sphe(5), sqhe(3), alphahe(4), betahe(4), sche,
     2          sdhe, sghe, shhe, skhe, slhe, srhe, sshe
 
      common /cdatl/ alc(4), blc(3), clc(3), dlc(4), elc(3), flc(3)
      common/cdats/ sac(5), sbc(3), sec(5), sfc(3), sic(5), sjc(3),
     1          spc(5), sqc(3), alphac(4), betac(4), scc,
     2          sdc, sgc, shc, skc, slc, src, ssc 
 
      common /odatl/ alo(4), blo(3), clo(3), dlo(4), elo(3), flo(3)
      common/odats/ sao(5), sbo(3), seo(5), sfo(3), sio(5), sjo(3),
     1          spo(5), sqo(3), alphao(4), betao(4), sco,
     2          sdo, sgo, sho, sko, slo, sro, sso 
*
* the following are the Bremsstrahlung neutrino rate fitting 
* coefficients quoted by Itoh. Common'd blocks are:
* Helium: soldath, phondath, liqdath, quandath
* Carbon: soldatc, phondatc, liqdatc, quandatc
* Oxygen: soldato, phondato, liqdato, quandato
*   Iron: soldatf, phondatf, liqdatf, quandatf 
*
      common/soldath/ snah(5), snbh(3), sneh(5), snfh(3), snih(5), 
     1          snjh(3), snph(5), snqh(3), alphasnh(4), 
     2          betasnh(4), snch, sndh, sngh, snhh, snkh, 
     3          snlh, snrh, snsh
      common/phondath/ pnah(5), pnbh(3), pnih(5), pnjh(3), 
     1          alphapnh(4), betapnh(4), pnch, pndh, pnkh, pnlh
      common/liqdath/ lnah(6), lnbh(4), lneh(6), lnfh(4), lnih(6), 
     1          lnjh(4), lnph(6), lnqh(4), alphalnh(4), 
     2          betalnh(4), lnch, lndh, lngh, lnhh, lnkh, 
     3          lnlh, lnrh, lnsh
      common/quandath/ qnah(6), qnbh(4), qneh(6), qnfh(4), 
     1          alphaqnh(5), qnch, qndh, qngh, qnhh
 
      common/soldatc/ snac(5), snbc(3), snec(5), snfc(3), snic(5), 
     1          snjc(3), snpc(5), snqc(3), alphasnc(4), 
     2          betasnc(4), sncc, sndc, sngc, snhc, snkc, 
     3          snlc, snrc, snsc 
      common/phondatc/ pnac(5), pnbc(3), pnic(5), pnjc(3), 
     1          alphapnc(4), betapnc(4), pncc, pndc, pnkc, pnlc
      common/liqdatc/ lnac(6), lnbc(4), lnec(6), lnfc(4), lnic(6), 
     1          lnjc(4), lnpc(6), lnqc(4), alphalnc(4), 
     2          betalnc(4), lncc, lndc, lngc, lnhc, lnkc, 
     3          lnlc, lnrc, lnsc 
      common/quandatc/ qnac(6), qnbc(4), qnec(6), qnfc(4), 
     1          alphaqnc(5), qncc, qndc, qngc, qnhc
 
      common/soldato/ snao(5), snbo(3), sneo(5), snfo(3), snio(5), 
     1          snjo(3), snpo(5), snqo(3), alphasno(4), 
     2          betasno(4), snco, sndo, sngo, snho, snko, 
     3          snlo, snro, snso
      common/phondato/ pnao(5), pnbo(3), pnio(5), pnjo(3), 
     1          alphapno(4), betapno(4), pnco, pndo, pnko, pnlo
      common/liqdato/ lnao(6), lnbo(4), lneo(6), lnfo(4), lnio(6), 
     1          lnjo(4), lnpo(6), lnqo(4), alphalno(4), 
     2          betalno(4), lnco, lndo, lngo, lnho, lnko, 
     3          lnlo, lnro, lnso
      common/quandato/ qnao(6), qnbo(4), qneo(6), qnfo(4), 
     1          alphaqno(5), qnco, qndo, qngo, qnho

      common/sfodatf/ snaf(5), snbf(3), snef(5), snff(3), snif(5), 
     1          snjf(3), snpf(5), snqf(3), alphasnf(4), 
     2          betasnf(4), sncf, sndf, sngf, snhf, snkf, 
     3          snlf, snrf, snsf
      common/phondatf/ pnaf(5), pnbf(3), pnif(5), pnjf(3), 
     1          alphapnf(4), betapnf(4), pncf, pndf, pnkf, pnlf
      common/liqdatf/ lnaf(6), lnbf(4), lnef(6), lnff(4), lnif(6), 
     1          lnjf(4), lnpf(6), lnqf(4), alphalnf(4), 
     2          betalnf(4), lncf, lndf, lngf, lnhf, lnkf, 
     3          lnlf, lnrf, lnsf
      common/quandatf/ qnaf(6), qnbf(4), qnef(6), qnff(4), 
     1          alphaqnf(5), qncf, qndf, qngf, qnhf
 
      data an0,bk,pi/6.022169e23,1.380622e-16,3.1415926535e0/
*
*  variables used in condlhy
*
      data alhy / 0.34542, -.28157, 0.09184, -.03734 /
      data blhy / -.61919, 0.40004, -.16585 /
      data clhy / 0.35742, -.41151, 0.21552 /
      data dlhy / 0.21512, -.10843, -.00596, -.00950 /
      data elhy / -.36667, 0.14040, -.04588 /
      data flhy / 0.10493, -.09537, 0.04682 /
*
*  variables used in condlhe
*
      data alhe / 0.62199, -.16110, 0.15574, -.02893 /
      data blhe / -.65222, 0.48601, -.18266 /
      data clhe / 0.36580, -.52176, 0.26240 /
      data dlhe / 0.36090, 0.02576, 0.05061, 0.00015 /
      data elhe / -.40559, 0.15316, -.04058 /
      data flhe / 0.12140, -.11621, 0.04939 /
*
* variables used in condshe
* NOTE: the values stored in a(1), e(1), i(1) and p(1) are 1/2 the
*  values listed in the paper.  this is because they are used in the
*  formulae in the paper as a(1)/2, etc.  (same for carbon and oxygen).
*
      data sahe       / -.01006, 0.05975, -.04290, -.00341, -.00338 / 
      data sbhe       / 0.11689, 0.01560, -.00713 /
      data sche, sdhe / -.16902, 1.4525 /
      data sehe       / -.05144, 0.10844, -.04380, -.00663, -.00657 / 
      data sfhe       / 0.12937, 0.02888, -.00413 /
      data sghe, shhe / -.20406, 1.36576 /
      data sihe       / 0.12429, -.04498, -.06418, -.00931, -.00582 / 
      data sjhe       / 0.11189, 0.00562, -.01176 /
      data skhe, slhe / -.22355, 1.70151 /
      data sphe       / 0.09218, -.00644, -.06552, -.01190, -.00832 / 
      data sqhe       / 0.12575, 0.01598, -.00950 /
      data srhe, sshe / -.25532, 1.90219 /
      data alphahe    / 1.5611, -20.6104, 247.0830, -1003.2600 /
      data betahe     / 1.5874, -21.4772, 255.7040, -1028.9000 /
*
*  variables used in condlc
*
      data alc/  0.9896, -0.1851,  0.1019, -0.0360/
      data blc/ -0.8825,  0.6657, -0.3798/
      data clc/ -0.0915, -1.5848,  1.1882/
      data dlc/  0.4406, -0.0161, -0.0093, -0.0028/
      data elc/ -0.4821,  0.0826, -0.0557/
      data flc/ -0.5193, -0.0830,  0.0147/
*
*  variables used in condsc
*
      data sac /   0.03387,  0.01795, -0.04324, -0.00437, -0.00422/
      data sbc /   0.04865, -0.00473, -0.00865/
      data scc, sdc/-0.09378,  0.94180/ 
      data sec /  -0.00329,  0.08497, -0.06909, -0.00582, -0.00677/
      data sfc /   0.12695,  0.00854, -0.01197/
      data sgc, shc/-0.19961,  1.63031/ 
      data sic /   0.177955,-0.08588, -0.07995, -0.00957, -0.00256/
      data sjc /   0.03417, -0.03395, -0.02519/
      data skc, slc/-0.10127,  2.10721/ 
      data spc /   0.18606, -0.06547, -0.10440, -0.01172, -0.00448/
      data sqc /   0.10304, -0.02601, -0.02922/
      data src, ssc/-0.20326,  2.80479/ 
      data alphac/0.4786,  16.4194, -139.1920, 184.7580/
      data betac/ 0.6146, 11.7296, -89.8319, 32.0399/
*
*  variables used in condlo
*
      data alo/ 1.0779, -0.1838, 0.1059, -.0290/
      data blo/ -.9743,  0.6955, -.3966/
      data clo/ -.1040, -1.7692, 1.3546 /
      data dlo/ 0.4486, -0.0160, -.0014, 0.0039 / 
      data elo/ -.5193,  0.0822, -.0467 /
      data flo/ -.5403, -0.1022, 0.0416 /
*
*  variables used in condso
*
      data sao /     0.02312, 0.02087, -.03684, -.00366, -.00348 /
      data sbo /     0.04177, -.00712, -.00865 /
      data sco, sdo/ -.07178, 0.85073 / 
      data seo /     -.001185, 0.08552, -.07206, -.00554, -.00673 /
      data sfo /     0.12586, 0.00371, -.01360 /
      data sgo, sho/ -.19226, 1.66185 / 
      data sio /      .160905, -.08292, -.07377, -.00577, 0.00154/
      data sjo /     0.03816, -.04206, -.02870 /
      data sko, slo/ -.06064, 2.14532 / 
      data spo /      .19600, -.07598, -.10917, -.00942, -.00142 /
      data sqo /     0.10403, -.03774, -.03483 /
      data sro, sso/ -.17568, 2.99352 / 
      data alphao/    .4682, 18.0197, -171.0090, 313.8890 / 
      data betao/    0.6085, 12.9705, -115.3180, 136.3590 / 
*
* variables used in bremsolh
*
      data snah       / -.01148, 0.01601, -.00433, 0.00015, -.00034 / 
      data snbh       / 0.01558, 0.00191, -.00055 /
      data snch, sndh / -.01694, 0.10649 /
      data sneh       / -.01827, 0.02395, -.00448, -.00033, -.00088 / 
      data snfh       / 0.01730, 0.00402, -.00005 /
      data sngh, snhh / -.02222, 0.13969 /
      data snih       / -.00323, 0.00440, -.00110, 0.00001, -.00007 / 
      data snjh       / 0.00294, 0.00059, -.00018 /
      data snkh, snlh / -.00337, 0.02116 /
      data snph       / -.00469, 0.00610, -.00114, -.00010, -.00018 / 
      data snqh       / 0.00320, 0.00107, -.00008 /
      data snrh, snsh / -.00442, 0.02775 /
      data alphasnh   / 1.6449, -23.2588, 272.1670, -1074.7000 /
      data betasnh    / 1.6443, -23.2414, 272.0080, -1074.2500 /
      data pnah       / -.00687, 0.00957, -.00204, -.00005, -.00003 / 
      data pnbh       / 0.00661, 0.00135, -.00035 /
      data pnch, pndh / -.00811, 0.05098 /
      data pnih       / -.00169, 0.00231, -.00047, -.00003, 0.00000 / 
      data pnjh       / 0.00111, 0.00042, -.00010 /
      data pnkh, pnlh / -.00161, 0.01013 /
      data alphapnh   / -.1394, 7.0680, -115.5940, 619.9170 /
      data betapnh    / -.1394, 7.0664, -115.5800, 619.8790 /
*
* variables used in bremsolc
*
      data snac       / 0.01838, -.01066, -.00458, -.00177, -.00138 / 
      data snbc       / -.00244, -.00206, -.00037 /
      data sncc, sndc / -.01093, 0.12431 /
      data snec       / 0.02360, -.01353, -.00619, -.00211, -.00176 / 
      data snfc       / 0.00456, -.00174, -.00031 /
      data sngc, snhc / -.02259, 0.20343 /
      data snic       / 0.00053, -.00048, -.00022, 0.00019, -.00001 / 
      data snjc       / 0.00658, -.00180, 0.00036 /
      data snkc, snlc / -.00398, 0.02499 /
      data snpc       / -.00024, 0.00063, -.00064, 0.00030, -.00006 / 
      data snqc       / 0.01013, -.00247, 0.00052 /
      data snrc, snsc / -.00650, 0.04087 /
      data alphasnc   / 0.6252, 10.6819, -70.6879, -44.3349 /
      data betasnc    / 0.6307, 10.4966, -68.7973, -50.0581 /
      data pnac       / 0.01116, -.00589, -.00279, -.00073, -.00043 / 
      data pnbc       / -.00095, -.00059, 0.00002 /
      data pncc, pndc / -.00729, 0.06630 /
      data pnic       / 0.00012, 0.00018, -.00028, 0.00012, -.00004 / 
      data pnjc       / 0.00339, -.00082, 0.00015 /
      data pnkc, pnlc / -.00212, 0.01332 /
      data alphapnc   / 0.5481, -20.4731, 223.9220, -534.9400 /
      data betapnc    / 0.5413, -20.2069, 220.7060, -524.1240 /
*
* variables used in bremsolo
*
      data snao       / 0.01616, -.00874, -.00413, -.00190, -.00139 / 
      data snbo       / -.00344, -.00261, -.00070 /
      data snco, sndo / -.00791, 0.13980 /
      data sneo       / 0.02210, -.00883, -.00857, -.00257, -.00214 / 
      data snfo       / 0.00629, -.00210, -.00099 /
      data sngo, snho / -.02610, 0.26993 /
      data snio       / 0.00100, -.00112, -.00003, 0.00014, 0.00001 / 
      data snjo       / 0.00745, -.00209, 0.00044 /
      data snko, snlo / -.00447, 0.02811 /
      data snpo       / -.00056, 0.00110, -.00074, 0.00024, -.00004 / 
      data snqo       / 0.01286, -.00281, 0.00057 /
      data snro, snso / -.00861, 0.05414 /
      data alphasno   / 0.4889, 16.1962, -138.4860, 185.7060 /
      data betasno    / 0.5111, 15.4195, -130.1540, 159.6050 /
      data pnao       / 0.00800, -.00191, -.00330, -.00075, -.00047 / 
      data pnbo       / 0.00088, -.00098, -.00036 /
      data pnco, pndo / -.00776, 0.08995 /
      data pnio       / -.00009, 0.00055, -.00038, 0.00011, -.00003 / 
      data pnjo       / 0.00429, -.00088, 0.00014 /
      data pnko, pnlo / -.00287, 0.01803 /
      data alphapno   / 0.3173, -14.4048, 186.9100, -476.8100 /
      data betapno    / 0.3073, -13.7973, 176.9940, -438.7520 /
*
* variables used in bremsolf
*
      data snaf       / 0.02096, -.01768, -.00007, -.00241, -.00080 / 
      data snbf       / -.01705, -.00268, -.00141 /
      data sncf, sndf / 0.00818, 0.13629 /
      data snef       / 0.04628, -.03290, -.00523, -.00539, -.00276 / 
      data snff       / -.02574, -.00630, -.00285 /
      data sngf, snhf / 0.00022, 0.31871 /
      data snif       / 0.00489, -.00653, 0.00171, -.00024, 0.00017 / 
      data snjf       / 0.00869, -.00323, 0.00075 /
      data snkf, snlf / -.00439, 0.02766 /
      data snpf       / 0.00732, -.00957, 0.00222, -.00022, 0.00025 / 
      data snqf       / 0.01867, -.00615, 0.00133 /
      data snrf, snsf / -.01023, 0.06442 /
      data alphasnf   / 0.6798, 12.7527, -140.1800, 268.8290 /
      data betasnf    / 0.7783, 10.2315, -124.2640, 241.3060 /
      data pnaf       / 0.01292, -.00894, -.00097, -.00125, -.00045 / 
      data pnbf       / -.00739, -.00190, -.00062 /
      data pncf, pndf / 0.00167, 0.09950 /
      data pnif       / 0.00240, -.00271, 0.00047, 0.00001, -.00001 / 
      data pnjf       / 0.00604, -.00197, 0.00044 /
      data pnkf, pnlf / -.00320, 0.02012 /
      data alphapnf   / 0.2847, -13.0828, 192.5030, -543.1480 /
      data betapnf    / 0.3221, -13.8640, 201.0700, -573.1160 /
*
* variables used in bremliqh
*
      data lnah       / 0.04518, -.03009, -.00564, -.00544, -.00290,
     1                  -.00224 /
      data lnbh       / -.02148, -.00817, -.00300, -.00170 /
      data lnch, lndh / 0.00671, 0.28130 /
      data lneh       / -.01003, 0.01790, -.00783, -.00021, 0.00024,
     1                  -.00014 / 
      data lnfh       / 0.00538, -.00175, -.00346, -.00031 /
      data lngh, lnhh / -.02199, 0.17300 /
      data lnih       / 0.00096, -.00301, -.00073, 0.00182, 0.00037,
     1                  0.00116 / 
      data lnjh       / 0.01706, -.00753, 0.00066, -.00060 /
      data lnkh, lnlh / -.01021, 0.06417 /
      data lnph       / -.00556, 0.00603, -.00149, 0.00047, 0.00040,
     1                  0.00028 / 
      data lnqh       / 0.00422, -.00009, -.00066, -.00003 /
      data lnrh, lnsh / -.00561, 0.03522 /
      data alphalnh   / -.07980, 0.17057, 1.51980, -.61058 /
      data betalnh    / -.05881, 0.00165, 1.82700, -.76993 /
      data qnah       / -.02260, 0.02634, -.00424, 0.00055, 0.00017, 
     1                  -.00042 / 
      data qnbh       / 0.00550, 0.00105, -.00179, 0.00028 /
      data qnch, qndh / -.02408, 0.28763 /
      data qneh       / -.00539, 0.00617, -.00205, 0.00089, 0.00022,
     1                  0.00032 / 
      data qnfh       / 0.00615, -.00149, -.00055, -.00004 /
      data qngh, qnhh / -.00764, 0.53346 /
      data alphaqnh   / 2.5344, -7.7188, -9.5605, -5.6433, 1.2671 /
*
* variables used in bremliqc
*
      data lnac       / 0.08973, -.05821, -.01089, -.01147, -.00656, 
     1                  -.00519 / 
      data lnbc       / -.04969, -.01584, -.00504, -.00281 /
      data lncc, lndc / 0.00945, 0.34529 /
      data lnec       / 0.03390, -.00944, -.01289, -.00589, -.00404,
     1                  -.00330 / 
      data lnfc       / -.02213, -.01136, -.00467, -.00131 /
      data lngc, lnhc / -.02342, 0.24819 /
      data lnic       / 0.00383, -.00710, -.00028, 0.00232, 0.00044,
     1                  0.00158 / 
      data lnjc       / 0.02300, -.01078, 0.00118, -.00089 /
      data lnkc, lnlc / -.01259, 0.07917 /
      data lnpc       / -.00384, 0.00356, -.00184, 0.00146, 0.00031,
     1                  0.00069 / 
      data lnqc       / 0.01052, -.00354, -.00014, -.00018 /
      data lnrc, lnsc / -.00829, 0.05211 /
      data alphalnc   / -.05483, -.01946, 1.86310, -.78873 /
      data betalnc    / -.06711, 0.06859, 1.74360, -.74498 /
      data qnac       / -.01900, 0.02446, -.00419, -.00030, -.00047, 
     1                  -.00099 / 
      data qnbc       / 0.00093, 0.00060, -.00176, 0.00018 /
      data qncc, qndc / -.02458, 0.25843 /
      data qnec       / -.00690, 0.01046, -.00351, 0.00010, -.00003,
     1                  -.00026 / 
      data qnfc       / 0.00289, -.00082, -.00155, -.00008 /
      data qngc, qnhc / -.01194, 0.34437 /
      data alphaqnc   / 0.9859, 2.2468, -10.6808, 11.4582, -4.0100 /
*
* variables used in bremliqo
*
      data lnao       / 0.10466, -.06740, -.01293, -.01352, -.00776,
     1                  -.00613 / 
      data lnbo       / -.05950, -.01837, -.00567, -.00310 /
      data lnco, lndo / 0.00952, 0.36029 /
      data lneo       / 0.04652, -.01656, -.01489, -.00778, -.00520,
     1                  -.00418 / 
      data lnfo       / -.03076, -.01390, -.00522, -.00161 /
      data lngo, lnho / -.02513, 0.27480 /
      data lnio       / 0.00476, -.00838, -.00011, 0.00244, 0.00046,
     1                  0.00168 / 
      data lnjo       / 0.02455, -.01167, 0.00132, -.00097 /
      data lnko, lnlo / -.01314, 0.08263 /
      data lnpo       / -.00350, 0.00295, -.00184, 0.00166, 0.00032,
     1                  0.00082 / 
      data lnqo       / 0.01231, -.00445, 0.00002, -.00026 /
      data lnro, lnso / -.00921, 0.05786 /
      data alphalno   / -.06597, 0.06048, 1.74860, -.74308 /
      data betalno    / -.07356, 0.10865, 1.70150, -.73653 /
      data qnao       / -.01792, 0.02395, -.00425, -.00055, -.00065, 
     1                  -.00115 / 
      data qnbo       / -.00039, 0.00043, -.00177, 0.00015 /
      data qnco, qndo / -.02491, 0.25207 /
      data qneo       / -.00866, 0.01344, -.00448, -.00002, -.00009,
     1                  -.00041 /
      data qnfo       / 0.00291, -.00091, -.00200, -.00010 /
      data qngo, qnho / -.01541, 0.32800 /
      data alphaqno   / 0.8978, 2.5443, -10.3549, 10.2665, -3.3537 /
*
* variables used in bremliqf
*
      data lnaf       / 0.17444, -.11076, -.02349, -.02283, -.01250,
     1                  -.00971 / 
      data lnbf       / -.10661, -.02860, -.00785, -.00385 /
      data lncf, lndf / 0.00766, 0.40991 /
      data lnef       / 0.11580, -.05891, -.02531, -.01747, -.01021,
     1                  -.00778 / 
      data lnff       / -.07767, -.02516, -.00711, -.00260 /
      data lngf, lnhf / -.03076, 0.36908 /
      data lnif       / 0.00946, -.01493, -.00125, 0.00262, 0.00055,
     1                  0.00209 / 
      data lnjf       / 0.03034, -.01519, 0.00204, -.00135 /
      data lnkf, lnlf / -.01494, 0.09395 /
      data lnpf       / -.00006, -.00222, -.00104, 0.00225, 0.00037,
     1                  0.00140 / 
      data lnqf       / 0.02013, -.00883, 0.00090, -.00071 /
      data lnrf, lnsf / -.01249, 0.07850 /
      data alphalnf   / -.07608, 0.11559, 1.75730, -.79677 /
      data betalnf    / -.08034, 0.14368, 1.73140, -.79467 /
      data qnaf       / -.01326, 0.02179, -.00484, -.00146, -.00135, 
     1                  -.00175 / 
      data qnbf       / -.00567, -.00040, -.00174, 0.00012 /
      data qncf, qndf / -.02710, 0.23276 /
      data qnef       / -.00430, 0.01443, -.00677, -.00142, -.00119,
     1                  -.00151 /
      data qnff       / -.00396, -.00245, -.00278, -.00023 /
      data qngf, qnhf / -.02288, 0.26787 /
      data alphaqnf   / 1.0379, 1.0195, -5.1242, 3.9369, -0.8701 /

      end 

************************************************************************

