program xstar1
c
c
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
      common /mtpass/ tgran
      common /linsel/ nlsv(nnnl),nlsvn
      common /prs   / p,p0
      common /temp  / t,to
      common /sigh  / zeta
      common /nmrc  / numrec,npass
      common /tlim  / tinf
      common /phrate/ pirt(nni)
      common /cdpth / dpthc(ncn)
      common /bdpth / dpthb(ncn)
      common /dpttau/ tauth
      common /tau0ln/ tau0(nnnl)
      common /xcol  / xcc(183)
      common /icc   / lichk(nni),lipin
      common /prtop / elnprnt(400),elimdb(2),nstpt,nlnprnt,lstpt(20)
      common /enerc / epi(ncn),dele(ncn),numcon
      common /spectc/ bremsa(ncn),brems(ncn),zrems(ncn),
     &                zremso(ncn),zremsz(ncn)
      common /radius/ delr,r,rl,rmax,rdel,radexp,rscale,rsave
      common /heato / httoto,cltoto,hmctoto
      common /epden / xee,xpx,xnx,xpxo,xnxo,xpxsave,colsave
      common /abion / xii(nnip),xiin(nni),xiip(nni),xiln(nnnl),
     &                xilp(nnnl),xiio(nnip)
      common /abel  / xel(nl),xeln(nni),xelln(nnnl)
      common /avquant/ tinner,taverage,tavhp,tavhpsq,rh2sv,corfac
      common /itdat / enfmxs,ensfrc,critd,crittd,epss,crits,
     $             kmaxs,lppris,nlimd,lpprid,nnmax
      common /comptn/ de(61,61),e(61),sx(61),dez(61),ez(61)
      common /rrrthh/ rrrth(11,11)
      common /spctcb/ zremsb(ncn),zrmsbo(ncn),bremsb(ncn),
     &                brmsab(ncn)
      common /numit / lnerrs,lnerrd,ntmpit,nelit,ntotit,
     &                npttit,nstpnm
      common /ffesc / fesc(nnnl)
      common /ffescb/ fescb(nnnl)
      common /lindat/ elin(nnnl),blin(nnnl), 
     &                sxlin(nnnl),eex(nnnl),cslin(nnnl),
     &      flin1(nnnl),flin2(nnnl),nblin(nnnl),nilin(nnnl),nlin(nni)
      common /pheat / htt(nni),htcmp
      common /pcool / cll(nni),clbr,clcmp
      common /ceemis/ ceem(nnnl)
      common /rcemis/ rcem(nnnl)
      common /heat  / httot,cltot,hmctot
c
      character*72 ktitle
c
      open(unit=9,file='sigtmp.dat',form='unformatted')
      open (unit=55,file='rskirs.dat',status='unknown')
c      open(unit=11,file='../atdat/MANSIG')
c
c      lpri=4
c
c     opening message
      write (6,*) 'xstar version 1.46'
      call pprint(13)
      write (6,*) 'doing initial setup...'
c
c
c     default parameter values
      radexp=0.
      ilstpt = 0
      elimdb(1) = 1.e-5
      elimdb(2) = 1.e+5
      ilndpt = 0
      nlimd = 300
      lffst=0
      nel = nl
      na = 15
      nnnl2 = nnnl
      nni2 = nni
      lthin = 0
      lpri = 0
      lwri = 0
      lfix = 0
      npass = 0
      lbcase = 0
      lnoinwd=0
      nstpt = 1
      lstpt(1) = 9
      t = 1.
      r = 0.
      xpxcol = 1.e+21
      xpx = 1.
      rmax = xpxcol
c      rmax=xpxcol/xpx
      lcdd = 1
      zeta = 0.
      xi = 10.**zeta
      xlum = 1.
      r19 = sqrt(xlum/amax1(1.e-34,xi*xpx))
      r = r19*1.e+19
c      rmax=xpxcol/(amax1(xpx,1.e-34))
      rmax = xpxcol
      lfast = 0
      numrec = 2
      xee = 1.099
c
      call ener
      xlumo = xlum
      xlum = 1.e-30
      tp = 10.
      call ispec
      xlum = xlumo
      tgran=0.1
c
c     main event loop
      lrfrst=1
 200    continue
c
c       read in
        call pprint(13)
        call rread1(lrfrst)
c
c       set up for iterative calculation if necessary
        numreco = numrec
        npass = max0(npass,1)
        rsave = r
        colsave=rmax
        xpxsave=xpx
        if (npass.gt.1) rewind(27)
        if (npass.gt.1) rewind(28)
        nptmp = 0
 500    nptmp = nptmp + 1
        numrectmp = numrec
        numrec = numreco-1
        if ( nptmp.gt.1 ) call shuffl
        if ( nptmp.gt.1 ) call rfnd0
        numrec = numrectmp
c
c       do set-ups
        call pprint(13)
        write (6,*) 'doing model setup, pass=',nptmp,'  out of ',
     &            npass
        call init2
        lmin = 0
c
c
        rdel = 0.
        r = rsave
        xpx=xpxsave
c
        delr = 0.
        lipin = 1
        call ichk
        call lescpe
        dilfac = 0.
        call diden
        n50kev = nbinc(6.e+4)
        ldon=0
c
c       correct rmax if lfix=2
        rmaxo=rmax
        if ( lfix.eq.2 ) 
     $   rmax=xpx*((1.+rmax/xpx/r)**(1./float(max0(1,numrec)))-1.)
     $         *float(numrec)
        rmaxt = rmax/xpx
c
        xpx0 = xpx
        r0 = r
c
c       setup for special quantities
        sumt=0.
        sumhpt=0.
        sumhp=0.
        sumhpsqt=0.
        sumhpsq=0.
        sumhesh=0.
        sumv=0.
        lh2fnd=0
c
c       step thru radius zones
        tinf = 0.31
        write (6,*) 'running ...'
        numrec2=numrec
c        if ((lfix.eq.1).and.(lcdd.ne.1)) numrec2=numrec*2
        if (lfix.eq.1) numrec2=numrec+1
        do 700 jk = 1,numrec2
c
c
          write (6,99001) jk
          if ( jk.eq.1 ) lipin = 1
          lipin = 1
          call ichk
c          lipin = 0
          call moveo
          ectt = 13.6
          emult = 0.1
c         the order of these statements matters!
          if (lfix.ge.1) emult=1.
          if (ldon.ge.1) emult=1.e-10
          if ((numrec.ge.10).and.(jk.le.5).and.(lfix.lt.1)) 
     $         emult=2.5e-5
          if (jk.eq.1) emult = 1.e-20
          call step(ectt,emult)
          if ( lfix.ge.1 ) 
     $        delr = (emult)*rmaxt/float(max0(1,numrec-1))
          rl = r
          r = r + delr
          if ( lfix.eq.2 ) r = r*(1.+delr*emult)
          delr2 = delr
          if ( lfix.eq.2 ) delr2 = r*delr
          rdel = rdel + delr2
          delr = delr2
          r19 = r/1.e+19
c          if ( lmin.eq.1 ) t = min(t,tinf)
          numrectmp = numrec
          numrec = numreco
          if ( nptmp.gt.1 ) call rfnd
          numrec = numrectmp
c
          if (abs(radexp).le.1.e-34) go to 9822
             if (lcdd.eq.1) xpx=xpx*((r+rscale)/(rl+rscale))**radexp
             if (lcdd.eq.0)  p0=p0*((r+rscale)/(rl+rscale))**radexp
 9822        continue
c
c
          call solve(jk)
c
c
c         assign special quantities
          vol=12.56*r19*r19*delr
          sumv=sumv+vol
          if (jk.eq.1) tinner=t
          sumt=sumt+vol*t
          volhp=vol*xii(2)*xel(1)
          sumhpt=sumhpt+volhp*t
          sumhp=sumhp+volhp
          volhpsq=volhp*xnx
          sumhpsqt=sumhpsqt+volhpsq*t
          sumhpsq=sumhpsq+volhpsq
          sumhesh=sumhesh+vol*xii(4)*xel(2)/(xii(2)+1.e-34)/xel(1)
          if ((xii(1).gt.0.999).and.(lh2fnd.eq.0)) rh2sv=xcc(183)
          if ((xii(1).gt.0.999).and.(lh2fnd.eq.0)) lh2fnd=1
c         write (6,*)r,delr,vol,lh2fnd,tinner,volhp,sumhpt,sumhp,
c     $      volhpsq,sumhpsqt,sumhpsq,sumhesh,rh2sv
          sumhrec=0.
          do 3011 ll=2,7
            sumhrec=sumhrec+rrrth(7,ll)
 3011       continue
          sumhrec=sumhrec*(13.6+(0.8617)*t)*(1.602197e-12)*xii(2)
     $            *xpx*xnx
c         write (6,*)'h cooling:',sumhrec,cll(1)
 3093     continue
c
c
c          call ispec9(tauth)
c
          call stpcut
c
          do 650 ll = 1,nstpt
            call pprint(lstpt(ll))
 650        continue
c          write (10,*)zeta,(pirt(mm),mm=115,140)
c
c         write(6,9696)httoto,cltoto
9696      format (1h ,'httot=',1pe10.2,' cltot=',
     $     1pe10.2,'(electron point of view)')
c
          if (npass.gt.1) call savo(28)
          if (lwri.eq.1) call savo(30)
          if ( t.le.(1.1*tinf) ) lmin = lmin+1
          if ( t.le.(1.1*tinf) ) t=3*tinf
          if (lmin.ne.0) lh2fnd=1
c
          if ( jk.eq.numrec ) ldon=ldon+1
          if ( (xcc(183).gt.rmax) .and. (lfix.ne.2) ) ldon=ldon+1
          if (ldon.ge.3) goto 800
          if ( lh2fnd.eq.1) ldon=ldon+1
c
c
c
 700      continue
 800    numreco = jk
c
        write (6,*) 'done'
        call pprint(13)
c
c       compute average quantities
        taverage=sumt/(1.e-34+sumv)
        tavhp=sumhpt/(1.e-34+sumhp)
        tavhpsq=sumhpsqt/(1.e-34+sumhpsq)
        corfac=sumhesh/(1.e-34+sumv)
c

        if ( nptmp.lt.npass ) goto 500
c
        if (lwri.eq.2) call savo(30)
c
        numrec=numreco-1
        if (lfix.eq.1) numrec=numrec-1
        if (lnoinwd.ne.1) go to 901
          do 900 ll = 1,numcon
            tmpe=expo(-dpthc(ll))
            if (lthin.eq.2) tmpe=1.
            zremsz(ll) = zremsz(ll)*tmpe
c     $           +zrems(ll)
 900        continue
 901      continue
c
c
        goto 200
c
99001 format (' ',' step number = ',i4)
c
      end
      function abswis(e)
c
      real abswis,e
      integer i
c
c   returns the photo-electric cross section at energy e (in kev)
c    in units of per 10**21 hydrogen atoms
c    as formulated in morrison and mccammon
c                                      a.e.szymkowiak jan 83
      real ery(14),c0(14),c1(14),c2(14)
      data ery/0.100,0.284,0.400,0.532,0.707,0.867,1.303,
     &     1.840,2.471,3.210,4.038,7.111,8.331,10.000/
      data c0/17.3,34.6,78.1,71.4,95.5,308.9,120.6,141.3,
     &     202.7,342.7,352.2,433.9,629.0,701.2/
      data c1/608.1,267.9,18.8,66.8,145.8,-380.6,169.3,
     &     146.8,104.7,18.7,18.7,-2.4,30.9,25.2/
      data c2/ - 2150.,-476.1,4.3,-51.4,-61.1,294.0,-47.7,
     &     -31.5,-17.0,0.0,0.0,0.75,0.0,0.0/
      if ( e.ge.10.000 ) then
         abswis = 0.001*(c0(14)+c1(14)*10.0)/(e*e*e)
      else
         do 50 i = 1,14
            if ( e.lt.ery(i) ) goto 100
 50      continue
 100     abswis = 0.001*(c0(i)+c1(i)*e+c2(i)*e*e)/(e*e*e)
      endif
      return
      end
      function alphdi(n,j,l,ln,t)
c
      real a,alphdi,b,blin,cslin,delt,dl,ebar,eex,
     &     elin,expo,sxlin,t,twork,x,z12,zf
      integer j,l,ln,n,nblin,nilin,nl,nlin,nni,nnnl
c
c
c  dielectronic recombination : burgess and tworkowski for h-like
c
      parameter (nnnl=3900,nni=168,nl=13)
c
      common /lindat/ elin(nnnl),blin(nnnl),
     &                sxlin(nnnl),eex(nnnl),cslin(nnnl),
     &    flin1(nnnl),flin2(nnnl),nblin(nnnl),nilin(nnnl),nlin(nni)
c
c  burgess and tworkowski
      z12 = j - 12.
      twork = .84 + .5/j**2 + .03*z12/(1.+4.5e-5*z12**3)
      if ( n.ne.j ) twork = 1.
      alphdi = 0.
      dl = delt(n,j,l)
      if ( dl.le.0. ) return
      zf = j - 1.
      b = sqrt(zf)*(zf+1.)**2.5/sqrt(zf*zf+13.4)
      x = eex(ln)/((zf+1.)*13.6)
      a = sqrt(x)/(1.+.105*x+.015*x*x)
      ebar = eex(ln)/(1.+.015*zf**3/(zf+1.)**2)
c
      alphdi = .0030*t**(-1.5)*a*b*flin1(ln)
     &         *twork*dl*expo(-ebar*11590./t)
      return
      end
      subroutine augcmp
c
c
c
c
c     this routine uses the photoionization rate matrices aumputed in
c     augfil, to compute the effective photoionization rate, as describ
c     in weisheit.
c
c
      external bigdat,newdat,rr3,nfllns
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
      common /aphrte/ aprtc(6,7),aprtn(7,8),aprto(8,9),
     &                aprtne(10,11),aprtsi(14,15),aprts(16,17),
     &                aprtfe(26,27)
      common /dirate/ dirth(1),dirthe(2),dirtc(6),dirtn(7),
     &                dirto(8),dirtne(10),dirtmg(12),dirtsi(14),
     &                dirts(16),dirtar(18),dirtca(20),dirtfe(26),
     &                dirtni(28)
      common /rrrate/ rrrth(1),rrrthe(2),rrrtc(6),rrrtn(7),
     &                rrrto(8),rrrtne(10),rrrtmg(12),rrrtsi(14),
     &                rrrts(16),rrrtar(18),rrrtca(20),rrrtfe(26),
     &                rrrtni(28)
      common /tbrate/ tbrth(1),tbrthe(2),tbrtc(6),tbrtn(7),
     &                tbrto(8),tbrtne(10),tbrtmg(12),tbrtsi(14),
     &                tbrts(16),tbrtar(18),tbrtca(20),tbrtfe(26),
     &                tbrtni(28)
      common /hcxrt / hxrh(2,1),hxrhe(2,2),hxrc(2,6),hxrn(2,7),
     &                hxro(2,8),hxrne(2,10),hxrmg(2,12),
     &                hxrsi(2,14),hxrs(2,16),hxrar(2,18),
     &                hxrca(2,20),hxrfe(2,26),hxrni(2,28),
     &                hexrh(2,1),hexrhe(2,2),hexrc(2,6),hexrn(2,7)
     &               ,hexro(2,8),hexrne(2,10),hexrmg(2,12),
     &                hexrsi(2,14),hexrs(2,16),hexrar(2,18),
     &                hexrca(2,20),hexrfe(2,26),hexrni(2,28)
      common /epden /xee,xpx,xnx,xpxo,xnxo,xpxsave,colsave
      common /abion / xii(nnip),xiin(nni),xiip(nni),xiln(nnnl),
     &                xilp(nnnl),xiio(nnip)
      common /abel  / xel(nl),xeln(nni),xelln(nnnl)
      common /phrate/ pirth(1),pirthe(2),pirtc(6),pirtn(7),
     &                pirto(8),pirtne(10),pirtmg(12),pirtsi(14),
     &                pirts(16),pirtar(18),pirtca(20),pirtfe(26),
     &                pirtni(28)
      common /icc   / lichk(nni),lipin
c
      character*72 ktitle
c
      dimension pirt(nni)
c
      equivalence (pirth(1),pirt(1))
c
      if ( nel.ge.3 ) then
c
         xh1 = xel(1)*xii(1)*xpx
         xh2 = xel(1)*xii(2)*xpx
         xhe1 = xel(2)*xii(3)*xpx
         xhe2 = xel(2)*xii(4)*xpx
c
         if ( lpri.gt.2 ) write (6,99001)
c
c     carbon
         pirtc(1) = aprtc(1,2)
         do 50 ik = 2,6
            prod = 1.
            pirtc(ik) = aprtc(ik,ik+1)
            ikm1 = ik - 1
            do 20 jk = 1,ikm1
               kk = ik - jk
               if ( lpri.gt.2 ) write (6,99002) ik,jk,kk,rrrtc(kk)
     &                                ,dirtc(kk),tbrtc(kk),
     &                                 hxrc(1,kk),hexrc(1,kk),
     &                                 pirtc(kk),prod
               temp = ((rrrtc(kk)+dirtc(kk)+tbrtc(kk))*xnx+hxrc(1,kk)
     &                *xh1+hexrc(1,kk)*xhe1)/pirtc(kk)
               if ( prod.lt.1.e+17 ) prod = prod*temp
               pirtc(ik) = pirtc(ik) + aprtc(kk,ik+1)*prod
 20         continue
 50      continue
c
         if ( nel.ge.4 ) then
c
c           nitrogen
            pirtn(1) = aprtn(1,2)
            do 80 ik = 2,7
               prod = 1.
               pirtn(ik) = aprtn(ik,ik+1)
               ikm1 = ik - 1
               do 60 jk = 1,ikm1
                  kk = ik - jk
                  if ( lpri.gt.2 ) write (6,99002) ik,jk,kk,
     &                 rrrtn(kk),dirtn(kk),tbrtn(kk),hxrn(1,kk),
     &                 hexrn(1,kk),pirtn(kk),prod
                  temp = ((rrrtn(kk)+dirtn(kk)+tbrtn(kk))*xnx+hxrn(1,kk)
     &                   *xh1+hexrn(1,kk)*xhe1)/pirtn(kk)
                  if ( prod.lt.1.e+17 ) prod = prod*temp
                  pirtn(ik) = pirtn(ik) + aprtn(kk,ik+1)*prod
 60            continue
 80         continue
c
            if ( nel.ge.5 ) then
c
c              oxygen
c               go to 3009
               pirto(1) = aprto(1,2) + xh2*hxro(2,1)
               do 90 ik = 2,8
                  prod = 1.
                  pirto(ik) = aprto(ik,ik+1)
                  ikm1 = ik - 1
                  do 85 jk = 1,ikm1
                     kk = ik - jk
                     if ( lpri.gt.2 ) write (6,99002) ik,jk,kk,
     &                    rrrto(kk),dirto(kk),tbrto(kk),hxro(1,kk)
     &                   ,hexro(1,kk),pirto(kk),prod
                     temp = ((rrrto(kk)+dirto(kk)+tbrto(kk))
     &                      *xnx+hxro(1,kk)*xh1+hexro(1,kk)*xhe1)
     &                      /pirto(kk)
c                     temp=min(temp,1.e+33/(1.e-34+prod))
                     if ( prod.lt.1.e+12 ) prod = prod*temp
c                     prod=prod*temp
                     pirto(ik) = pirto(ik) + aprto(kk,ik+1)*prod
 85               continue
 90            continue
 3009          continue
c
               if ( nel.ge.6 ) then
c
c     neon
                  pirtne(1) = aprtne(1,2)
                  do 95 ik = 2,10
                     prod = 1.
                     pirtne(ik) = aprtne(ik,ik+1)
                     ikm1 = ik - 1
                     do 92 jk = 1,ikm1
                        kk = ik - jk
                        if ( lpri.gt.2 ) write (6,99002) ik,jk,kk,
     &                       rrrtne(kk),dirtne(kk),tbrtne(kk),
     &                       hxrne(1,kk),hexrne(1,kk),pirtne(kk),
     &                       prod
                        temp = ((rrrtne(kk)+dirtne(kk)+tbrtne(kk))
     &                         *xnx+hxrne(1,kk)*xh1+hexrne(1,kk)*xhe1)
     &                         /pirtne(kk)
                        if ( prod.lt.1.e+17 ) prod = prod*temp
                        pirtne(ik) = pirtne(ik) + aprtne(kk,ik+1)*prod
 92                  continue
 95               continue
c
                  if ( nel.ge.8 ) then
c
c
c     silicon
                     pirtsi(1) = aprtsi(1,2)
                     do 98 ik = 2,14
                        prod = 1.
                        pirtsi(ik) = aprtsi(ik,ik+1)
                        ikm1 = ik - 1
                        do 96 jk = 1,ikm1
                           kk = ik - jk
                           if ( lpri.gt.2 ) write (6,99002) ik,jk,
     &                          kk,rrrtsi(kk),dirtsi(kk),
     &                          tbrtsi(kk),hxrsi(1,kk),hexrsi(1,kk)
     &                         ,pirtsi(kk),prod
                           temp = ((rrrtsi(kk)+dirtsi(kk)+tbrtsi(kk))
     &                            *xnx+hxrsi(1,kk)*xh1+hexrsi(1,kk)
     &                            *xhe1)/pirtsi(kk)
                           if ( prod.lt.1.e+17 ) prod = prod*temp
                           pirtsi(ik) = pirtsi(ik) + aprtsi(kk,ik+1)
     &                                  *prod
 96                     continue
 98                  continue
c
                     if ( nel.ge.9 ) then
c
c
c     sulfur
                        pirts(1) = aprts(1,2)
                        do 102 ik = 2,16
                           prod = 1.
                           pirts(ik) = aprts(ik,ik+1)
                           ikm1 = ik - 1
                           do 100 jk = 1,ikm1
                              kk = ik - jk
                              if ( lpri.gt.2 ) write (6,99002) ik,jk,
     &                             kk,rrrts(kk),dirts(kk),
     &                             tbrts(kk),hxrs(1,kk),hexrs(1,kk)
     &                            ,pirts(kk),prod
                              temp = ((rrrts(kk)+dirts(kk)+tbrts(kk))
     &                               *xnx+hxrs(1,kk)*xh1+hexrs(1,kk)
     &                               *xhe1)/pirts(kk)
                              if ( prod.lt.1.e+17 ) prod = prod*temp
                              pirts(ik) = pirts(ik) + aprts(kk,ik+1)
     &                           *prod
 100                       continue
 102                    continue
c
c
                        if ( nel.ge.12 ) then
c
c                          iron
                           pirtfe(1) = aprtfe(1,2)
                           do 106 ik = 2,26
                              prod = 1.
                              pirtfe(ik) = aprtfe(ik,ik+1)
                              ikm1 = ik - 1
                              do 104 jk = 1,ikm1
                                 kk = ik - jk
                                 if ( lpri.gt.2 ) write (6,99002) ik,
     &                                jk,kk,rrrtfe(kk),dirtfe(kk)
     &                               ,tbrtfe(kk),hxrfe(1,kk),
     &                                hexrfe(1,kk),pirtfe(kk),prod
                                 temp = ((rrrtfe(kk)+dirtfe(kk)+tbrtfe(
     &                                  kk))*xnx+hxrfe(1,kk)
     &                                  *xh1+hexrfe(1,kk)*xhe1)
     &                                  /pirtfe(kk)
                                 if ( prod.lt.1.e+17 ) prod = prod*temp
                                 pirtfe(ik) = pirtfe(ik)
     &                              + aprtfe(kk,ik+1)*prod
 104                          continue
 106                       continue
                        endif
                     endif
                  endif
               endif
            endif
         endif
      endif
c
c
c
c
      return
99001 format (' ',' in augcmp ')
99002 format (' ',3i4,7e12.4)
      end
      subroutine augfil
c
c
c
c
 
c     this routine fills the photoinization rates matrices, in aphrate.
c
c
      external bigdat,newdat,rr3,nfllns
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
      common /phheat/ pihth(1),pihthe(2),pihtc(6),pihtn(7),
     &                pihto(8),pihtne(10),pihtmg(12),pihtsi(14),
     &                pihts(16),pihtar(18),pihtca(20),pihtfe(26),
     &                pihtni(28)
      common /phrate/ pirth(1),pirthe(2),pirtc(6),pirtn(7),
     &                pirto(8),pirtne(10),pirtmg(12),pirtsi(14),
     &                pirts(16),pirtar(18),pirtca(20),pirtfe(26),
     &                pirtni(28)
      common /abion / xii(nnip),xih(1),xihe(2),xic(6),xin(7),
     &                xio(8),xine(10),ximg(12),xisi(14),xis(16)
     &               ,xiar(18),xica(20),xife(26),xini(28),
     &                xiip(nni),xiln(nnnl),xilp(nnnl),xiio(nnip)
      common /pqrate/ piqth(1),piqthe(2),piqtc(6),piqtn(7),
     &                piqto(8),piqtne(10),piqtmg(12),piqtsi(14),
     &                piqts(16),piqtar(18),piqtca(20),piqtfe(26),
     &                piqtni(28)
      common /pphht / phkh(1),phkhe(2),phkc(6),phlsc(4),phlpc(2)
     &               ,phkn(7),phlsn(5),phlpn(3),phko(8),
     &                phlso(6),phlpo(4),phkne(10),phlsne(8),
     &                phlpne(6),phkmg(12),phlsmg(10),phlpmg(8),
     &                phmsmg(2),phksi(14),phlssi(12),phlpsi(10),
     &                phmssi(4),phmpsi(2),phks(16),phlss(14),
     &                phlps(12),phmss(6),phmps(4),phkar(18),
     &                phlsar(16),phlpar(14),phmsar(8),phmpar(6),
     &                phkca(20),phlsca(18),phlpca(16),phmsca(10),
     &                phmpca(8),phmdca(2),phkfe(26),phlsfe(24),
     &                phlpfe(22),phmsfe(16),phmpfe(14),phmdfe(8),
     &                phnfe(2),phkni(28),phlsni(26),phlpni(24),
     &                phmsni(18),phmpni(16),phmdni(10),phnni(4)
      common /ppqrt / pqkh(1),pqkhe(2),pqkc(6),pqlsc(4),pqlpc(2)
     &               ,pqkn(7),pqlsn(5),pqlpn(3),pqko(8),
     &                pqlso(6),pqlpo(4),pqkne(10),pqlsne(8),
     &                pqlpne(6),pqkmg(12),pqlsmg(10),pqlpmg(8),
     &                pqmsmg(2),pqksi(14),pqlssi(12),pqlpsi(10),
     &                pqmssi(4),pqmpsi(2),pqks(16),pqlss(14),
     &                pqlps(12),pqmss(6),pqmps(4),pqkar(18),
     &                pqlsar(16),pqlpar(14),pqmsar(8),pqmpar(6),
     &                pqkca(20),pqlsca(18),pqlpca(16),pqmsca(10),
     &                pqmpca(8),pqmdca(2),pqkfe(26),pqlsfe(24),
     &                pqlpfe(22),pqmsfe(16),pqmpfe(14),pqmdfe(8),
     &                pqnfe(2),pqkni(28),pqlsni(26),pqlpni(24),
     &                pqmsni(18),pqmpni(16),pqmdni(10),pqnni(4)
      common /ppirt / prkh(1),prkhe(2),prkc(6),prlsc(4),prlpc(2)
     &               ,prkn(7),prlsn(5),prlpn(3),prko(8),
     &                prlso(6),prlpo(4),prkne(10),prlsne(8),
     &                prlpne(6),prkmg(12),prlsmg(10),prlpmg(8),
     &                prmsmg(2),prksi(14),prlssi(12),prlpsi(10),
     &                prmssi(4),prmpsi(2),prks(16),prlss(14),
     &                prlps(12),prmss(6),prmps(4),prkar(18),
     &                prlsar(16),prlpar(14),prmsar(8),prmpar(6),
     &                prkca(20),prlsca(18),prlpca(16),prmsca(10),
     &                prmpca(8),prmdca(2),prkfe(26),prlsfe(24),
     &                prlpfe(22),prmsfe(16),prmpfe(14),prmdfe(8),
     &                prnfe(2),prkni(28),prlsni(26),prlpni(24),
     &                prmsni(18),prmpni(16),prmdni(10),prnni(4)
      common /ethph / etkh(1),etkhe(2),etkc(6),etlsc(4),etlpc(2)
     &               ,etkn(7),etlsn(5),etlpn(3),etko(8),
     &                etlso(6),etlpo(4),etkne(10),etlsne(8),
     &                etlpne(6),etkmg(12),etlsmg(10),etlpmg(8),
     &                etmsmg(2),etksi(14),etlssi(12),etlpsi(10),
     &                etmssi(4),etmpsi(2),etks(16),etlss(14),
     &                etlps(12),etmss(6),etmps(4),etkar(18),
     &                etlsar(16),etlpar(14),etmsar(8),etmpar(6),
     &                etkca(20),etlsca(18),etlpca(16),etmsca(10),
     &                etmpca(8),etmdca(2),etkfe(26),etlsfe(24),
     &                etlpfe(22),etmsfe(16),etmpfe(14),etmdfe(8),
     &                etnfe(2),etkni(28),etlsni(26),etlpni(24),
     &                etmsni(18),etmpni(16),etmdni(10),etnni(4)
      common /spectc/ bremsa(ncn),brems(ncn),zrems(ncn),
     &                zremso(ncn),zremsz(ncn)
      common /aphrte/ aprtc(6,7),aprtn(7,8),aprto(8,9),
     &                aprtne(10,11),aprtsi(14,15),aprts(16,17),
     &                aprtfe(26,27)
      common /agdata/ amkc(6,6),amkn(7,7),amko(8,8),amkne(10,10),
     &                amkmg(12,12),amlsmg(10,10),amlpmg(8,8),
     &                amksi(14,14),amlssi(12,12),amlpsi(10,10),
     &                amks(16,16),amlss(14,14),amlps(12,12),
     &                amkfe(26,26),amlsfe(24,24),amlpfe(22,22)
      common /flpass/ fflmh(15,1),fflmhe(15,2),fflmc(15,6),
     &                fflmn(15,7),fflmo(15,8),fflmne(15,10),
     &                fflmmg(15,12),fflmsi(15,14),fflms(15,16),
     &                fflmfe(15,26)
      common /icc   / lichk(nni),lipin
c
      character*72 ktitle
c
      dimension agzro(1464)
      dimension flmzro(102)
      dimension flemh(1),flemc(6),flemn(7),flemo(8),flemne(10),
     &          flemsi(14),flems(16),flemfe(26)
      dimension aheat(26),aion(26)
      dimension nflnk(26),nflnls(26),nflnlp(26),nflnms(26),
     &          nflnmp(26),nflnmd(26)
      dimension piht(nni),piqt(nni)
      dimension esumhe(2),esumc(6),esumn(7),esumo(8),esumne(10)
     &         ,esumsi(14),esums(16),esumfe(26)
c
      equivalence (piht(1),pihth(1)),(piqt(1),piqth(1)),
     &             (agzro(1),aprtc(1,1)),(flmzro(1),flemh(1))
c
c
c
      data ergsev/1.602192e-12/
      data nflnk/3, 2*1,7,7,5*6,2*12,10,5*12,6*8,2*7/
      data nflnls/3*0,1,1,5*5,2*10,9,5*11,6*7,2*6/
      data nflnlp/10*0,2*11,11,5*10,6*6,2*5/
      data nflnms/12*0,1,5*9,6*5,2*4/
      data nflnmp/18*0,6*4,2*3/
      data nflnmd/24*0,2*2/
c     energy sumdata
      data esumhe/79.,54.6/
      data esumc/1030.,1019.,994.4,946.5,882.,490./
      data esumn/1486.,1472.,1442.,1394.,1317.,1219.,667./
      data esumo/2044.,2030.,1995.,1940.,1863.,1749.,1611.,
     &     871.4/
      data esumne/3512.,3490.,3449.,3386.,3289.,3162.,
     &     3004.,2797.,2558.,1362./
      data esumsi/7888.,7880.,7863.,7830.,7785.,7618.,
     &     7413.,7166.,6863.,6512.,6110.,5634.,5111.,
     &     2673./
      data esums/10858.,10848.,10824.,10789.,10742.,10670.,
     &     10581.,10300.,9972.,9593.,9146.,8641.,8077.,
     &     7425.,6718.,3494./
      data esumfe/34628.,34620.,34604.,34573.,34518.,34443.,
     &     34344.,34219.,34068.,33833.,33571.,33281.,
     &     32950.,32589.,32197.,31740.,31250.,29985.,
     &     28627.,27171.,25589.,23900.,22101.,20151.,
     &     18106.,9278./
c
c
c     zero the photoionization rate matrices.
      do 100 jk = 1,1464
         agzro(jk) = 0.
 100  continue
      do 200 jk = 1,102
         flmzro(jk) = 0.
 200  continue
c
      if ( lpri.gt.2 ) write (6,99001)
c
      if ( nel.ge.3 ) then
c
c     carbon
         call aughlp(aprtc,amkc,aheat,aion,prkc,esumc,etkc,flemc,1,6,6)
         do 250 ik = 1,4
            nfl = nflnk(6-ik)
            if ( lpri.gt.2 ) write (6,99002) ik,nfl,flemc(ik)
            if ( (nfl.ne.0) .and. (ik.ne.6) ) fflmc(nfl,ik+1)
     &           = flemc(ik)
            piqtc(ik) = pqkc(ik) + aion(ik)
            pihtc(ik) = phkc(ik) + aheat(ik)
 250     continue
c     l shell
         do 300 ik = 1,4
            nfl = nflnls(6-ik)
            if ( (nfl.ne.0) .and. (ik.ne.6) ) fflmc(nfl,ik+1)
     &           = prlsc(ik)
            aprtc(ik,ik+1) = aprtc(ik,ik+1) + prlsc(ik)
            piqtc(ik) = piqtc(ik) + pqlsc(ik)
            pihtc(ik) = pihtc(ik) + phlsc(ik)
 300     continue
         do 350 ik = 1,2
            nfl = nflnlp(6-ik)
            if ( (nfl.ne.0) .and. (ik.ne.6) ) fflmc(nfl,ik+1)
     &           = aprtc(ik,ik+1) + prlpc(2)
            aprtc(ik,ik+1) = aprtc(ik,ik+1) + prlpc(ik)
            piqtc(ik) = piqtc(ik) + pqlpc(ik)
            pihtc(ik) = pihtc(ik) + phlpc(ik)
 350     continue
c
c
         if ( nel.ge.4 ) then
c
c     nitrogen
c     k shell
            call aughlp(aprtn,amkn,aheat,aion,prkn,esumn,etkn,flemn,1,7,
     &                  7)
            do 360 ik = 1,5
               nfl = nflnk(7-ik)
               if ( (nfl.ne.0) .and. (ik.ne.7) ) fflmn(nfl,ik+1)
     &              = flemn(ik)
               piqtn(ik) = pqkn(ik) + aion(ik)
               pihtn(ik) = phkn(ik) + aheat(ik)
 360        continue
c     l shell
            do 380 ik = 1,5
               nfl = nflnls(7-ik)
               if ( (nfl.ne.0) .and. (ik.ne.7) ) fflmn(nfl,ik+1)
     &              = prlsn(ik)
               aprtn(ik,ik+1) = aprtn(ik,ik+1) + prlsn(ik)
               piqtn(ik) = piqtn(ik) + pqlsn(ik)
               pihtn(ik) = pihtn(ik) + phlsn(ik)
 380        continue
            do 400 ik = 1,3
               nfl = nflnlp(7-ik)
               if ( (nfl.ne.0) .and. (ik.ne.7) ) fflmn(nfl,ik+1)
     &              = prlpn(ik)
               aprtn(ik,ik+1) = aprtn(ik,ik+1) + prlpn(ik)
               piqtn(ik) = piqtn(ik) + pqlpn(ik)
               pihtn(ik) = pihtn(ik) + phlpn(ik)
 400        continue
c
            if ( nel.ge.5 ) then
c
c     oxygen
c     k shell
               call aughlp(aprto,amko,aheat,aion,prko,esumo,etko,flemo,
     &                     1,8,8)
               do 410 ik = 1,6
                  nfl = nflnk(8-ik)
                  if ( (nfl.ne.0) .and. (ik.ne.8) ) fflmo(nfl,ik+1)
     &                 = flemo(ik)
                  piqto(ik) = pqko(ik) + aion(ik)
                  pihto(ik) = phko(ik) + aheat(ik)
 410           continue
c     l shell
               do 420 ik = 1,6
                  nfl = nflnls(8-ik)
                  if ( (nfl.ne.0) .and. (ik.ne.8) ) fflmo(nfl,ik+1)
     &                 = prlso(ik)
                  aprto(ik,ik+1) = aprto(ik,ik+1) + prlso(ik)
                  piqto(ik) = piqto(ik) + pqlso(ik)
                  pihto(ik) = pihto(ik) + phlso(ik)
 420           continue
               do 430 ik = 1,4
                  nfl = nflnlp(8-ik)
                  if ( (nfl.ne.0) .and. (ik.ne.8) ) fflmo(nfl,ik+1)
     &                 = prlpo(ik)
                  aprto(ik,ik+1) = aprto(ik,ik+1) + prlpo(ik)
                  piqto(ik) = piqto(ik) + pqlpo(ik)
                  pihto(ik) = pihto(ik) + phlpo(ik)
 430           continue
c
c
               if ( nel.ge.6 ) then
c     neon
c     k shell
                  call aughlp(aprtne,amkne,aheat,aion,prkne,esumne,
     &                        etkne,flemne,1,10,10)
                  do 435 ik = 1,8
                     nfl = nflnk(10-ik)
                     if ( (nfl.ne.0) .and. (ik.ne.10) ) fflmne(nfl,ik+1)
     &                    = flemne(ik)
                     piqtne(ik) = pqkne(ik) + aion(ik)
                     pihtne(ik) = phkne(ik) + aheat(ik)
 435              continue
c     l shell
                  do 440 ik = 1,8
                     nfl = nflnls(10-ik)
                     if ( (nfl.ne.0) .and. (ik.ne.10) ) fflmne(nfl,ik+1)
     &                    = prlsne(ik)
                     aprtne(ik,ik+1) = aprtne(ik,ik+1) + prlsne(ik)
                     piqtne(ik) = piqtne(ik) + pqlsne(ik)
                     pihtne(ik) = pihtne(ik) + phlsne(ik)
 440              continue
                  do 445 ik = 1,6
                     nfl = nflnlp(10-ik)
                     if ( (nfl.ne.0) .and. (ik.ne.10) ) fflmne(nfl,ik+1)
     &                    = prlpne(ik)
                     aprtne(ik,ik+1) = aprtne(ik,ik+1) + prlpne(ik)
                     piqtne(ik) = piqtne(ik) + pqlpne(ik)
                     pihtne(ik) = pihtne(ik) + phlpne(ik)
 445              continue
c
                  if ( nel.ge.7 ) then
c     silicon
c     k shell
                     call aughlp(aprtsi,amksi,aheat,aion,prksi,esumsi,
     &                           etksi,flemsi,1,14,14)
c                     write (6,*)'si k shell rates:'
                     do 446 ik = 1,12
                        nfl = nflnk(14-ik)
                        if ( (nfl.ne.0) .and. (ik.ne.14) )
     &                       fflmsi(nfl,ik+1) = flemsi(ik)
c                        write (6,*)ik,nfl,flemsi(ik)
                        piqtsi(ik) = pqksi(ik) + aion(ik)
                        pihtsi(ik) = phksi(ik) + aheat(ik)
 446                 continue
c     ls shell
                     call aughlp(aprtsi,amlssi,aheat,aion,prlssi,esumsi,
     &                           etlssi,flemsi,1,12,14)
                     do 448 ik = 1,10
                        nfl = nflnls(14-ik)
                        if ( (nfl.ne.0) .and. (ik.ne.14) )
     &                       fflmsi(nfl,ik+1) = flemsi(ik)
                        pihtsi(ik) = pihtsi(ik) + phlssi(ik) + aheat(ik)
                        piqtsi(ik) = piqtsi(ik) + pqlssi(ik) + aion(ik)
 448                 continue
c     lp shell
                     call aughlp(aprtsi,amlpsi,aheat,aion,prlpsi,esumsi,
     &                           etlpsi,flemsi,1,10,14)
                     do 450 ik = 1,4
                        nfl = nflnlp(14-ik)
                        if ( (nfl.ne.0) .and. (ik.ne.14) )
     &                       fflmsi(nfl,ik+1) = flemsi(ik)
                        pihtsi(ik) = pihtsi(ik) + phlpsi(ik) + aheat(ik)
                        piqtsi(ik) = piqtsi(ik) + pqlpsi(ik) + aion(ik)
 450                 continue
c     ms shell
                     do 452 ik = 1,4
                        nfl = nflnms(14-ik)
                        if ( (nfl.ne.0) .and. (ik.ne.14) )
     &                       fflmsi(nfl,ik+1) = prmssi(ik)
                        aprtsi(ik,ik+1) = aprtsi(ik,ik+1) + prmssi(ik)
 452                 continue
                     do 454 ik = 1,4
                        piqtsi(ik) = piqtsi(ik) + pqmssi(ik)
                        pihtsi(ik) = pihtsi(ik) + phmssi(ik)
 454                 continue
c     mp shell
                     do 456 ik = 1,2
                        nfl = nflnmp(14-ik)
                        if ( (nfl.ne.0) .and. (ik.ne.14) )
     &                       fflmsi(nfl,ik+1) = prmpsi(ik)
                        aprtsi(ik,ik+1) = aprtsi(ik,ik+1) + prmpsi(ik)
 456                 continue
                     do 458 ik = 1,2
                        pihtsi(ik) = pihtsi(ik) + phmpsi(ik)
                        piqtsi(ik) = piqtsi(ik) + pqmpsi(ik)
 458                 continue
c
                     if ( nel.ge.8 ) then
c
c     sulfur
c     k shell
                        call aughlp(aprts,amks,aheat,aion,prks,esums,
     &                              etks,flems,1,16,16)
                        do 460 ik = 1,14
                           nfl = nflnk(16-ik)
                           if ( (nfl.ne.0) .and. (ik.ne.16) )
     &                          fflms(nfl,ik+1) = flems(ik)
                           piqts(ik) = pqks(ik) + aion(ik)
                           pihts(ik) = phks(ik) + aheat(ik)
 460                    continue
c     ls shell
                        call aughlp(aprts,amlss,aheat,aion,prlss,esums,
     &                              etlss,flems,1,14,16)
                        do 462 ik = 1,12
                           nfl = nflnls(16-ik)
                           if ( (nfl.ne.0) .and. (ik.ne.16) )
     &                          fflms(nfl,ik+1) = flems(ik)
                           piqts(ik) = piqts(ik) + pqlss(ik) + aion(ik)
                           pihts(ik) = pihts(ik) + phlss(ik) + aheat(ik)
 462                    continue
c     lp shell
                        call aughlp(aprts,amlps,aheat,aion,prlps,esums,
     &                              etlps,flems,1,12,16)
                        do 464 ik = 1,6
                           nfl = nflnlp(16-ik)
                           if ( (nfl.ne.0) .and. (ik.ne.16) )
     &                          fflms(nfl,ik+1) = flems(ik)
                           piqts(ik) = piqts(ik) + pqlps(ik) + aion(ik)
                           pihts(ik) = pihts(ik) + phlps(ik) + aheat(ik)
 464                    continue
c     ms shell
                        do 466 ik = 1,6
                           nfl = nflnms(16-ik)
                           if ( (nfl.ne.0) .and. (ik.ne.16) )
     &                          fflms(nfl,ik+1) = prmss(ik)
                           aprts(ik,ik+1) = aprts(ik,ik+1) + prmss(ik)
 466                    continue
                        do 468 ik = 1,6
                           piqts(ik) = piqts(ik) + pqmss(ik)
                           pihts(ik) = pihts(ik) + phmss(ik)
 468                    continue
c     mp shell
                        do 470 ik = 1,4
                           nfl = nflnmp(16-ik)
                           if ( (nfl.ne.0) .and. (ik.ne.16) )
     &                          fflms(nfl,ik+1) = prmps(ik)
                           aprts(ik,ik+1) = aprts(ik,ik+1) + prmps(ik)
 470                    continue
                        do 472 ik = 1,4
                           piqts(ik) = piqts(ik) + pqmps(ik)
                           pihts(ik) = pihts(ik) + phmps(ik)
 472                    continue
c
c
                        if ( nel.ge.12) then
c
c                          iron
c                          k shell
                           call aughlp(aprtfe,amkfe,aheat,aion,prkfe,
     &                                 esumfe,etkfe,flemfe,1,26,26)
                           do 474 ik = 1,25
                              nfl = nflnk(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = flemfe(ik)
                              nfl2 = nflnls(26-ik)
c                             here's the l shell fluorescence fake
                              if ( (nfl2.ne.0) .and. (ik.le.14) )
     &                             fflmfe(nfl2,ik+1) = flemfe(ik)
     &                             *amlsfe(ik,ik)
                              nfl2 = nflnlp(26-ik)
c                             here's the l shell fluorescence fake
                              if ( (nfl2.ne.0) .and. (ik.le.15) )
     &                             fflmfe(nfl2,ik+1) = flemfe(ik)
     &                             *amlpfe(ik,ik)
                              if ( lpri.gt.2 ) write (6,99003) ik,
     &                             nfl,prkfe(ik),flemfe(ik)
                              piqtfe(ik) = pqkfe(ik) + aion(ik)
                              pihtfe(ik) = phkfe(ik) + aheat(ik)
 474                       continue
c                          ls shell
                           call aughlp(aprtfe,amlsfe,aheat,aion,prlsfe,
     &                                 esumfe,etlsfe,flemfe,1,24,26)
                           do 476 ik = 1,22
                              nfl = nflnls(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = flemfe(ik)
                              pihtfe(ik) = pihtfe(ik) + phlsfe(ik)
     &                           + aheat(ik)
                              piqtfe(ik) = piqtfe(ik) + pqlsfe(ik)
     &                           + aion(ik)
 476                       continue
c                          lp shell
                           call aughlp(aprtfe,amlpfe,aheat,aion,prlpfe,
     &                                 esumfe,etlpfe,flemfe,1,22,26)
                           do 478 ik = 1,22
                              nfl = nflnlp(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = flemfe(ik)
                              piqtfe(ik) = piqtfe(ik) + pqlpfe(ik)
     &                           + aion(ik)
                              pihtfe(ik) = pihtfe(ik) + phlpfe(ik)
     &                           + aheat(ik)
 478                       continue
c                          ms shell
                           do 480 ik = 1,16
                              nfl = nflnms(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = prmsfe(ik)
                              aprtfe(ik,ik+1) = aprtfe(ik,ik+1)
     &                           + prmsfe(ik)
 480                       continue
                           do 482 ik = 1,16
                              piqtfe(ik) = piqtfe(ik) + pqmsfe(ik)
                              pihtfe(ik) = pihtfe(ik) + phmsfe(ik)
 482                       continue
c                          mp shell
                           do 484 ik = 1,14
                              nfl = nflnmp(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = prmpfe(ik)
                              aprtfe(ik,ik+1) = aprtfe(ik,ik+1)
     &                           + prmpfe(ik)
 484                       continue
                           do 486 ik = 1,14
                              piqtfe(ik) = piqtfe(ik) + pqmpfe(ik)
                              pihtfe(ik) = pihtfe(ik) + phmpfe(ik)
 486                       continue
c     md shell
                           do 488 ik = 1,8
                              nfl = nflnmd(26-ik)
                              if ( (nfl.ne.0) .and. (ik.ne.26) )
     &                             fflmfe(nfl,ik+1) = prmdfe(ik)
                              aprtfe(ik,ik+1) = aprtfe(ik,ik+1)
     &                           + prmdfe(ik)
 488                       continue
                           do 490 ik = 1,8
                              piqtfe(ik) = piqtfe(ik) + pqmdfe(ik)
                              pihtfe(ik) = pihtfe(ik) + phmdfe(ik)
 490                       continue
c     n shell
                           do 492 ik = 1,2
                              aprtfe(ik,ik+1) = aprtfe(ik,ik+1)
     &                           + prnfe(ik)
 492                       continue
                           do 494 ik = 1,2
                              piqtfe(ik) = piqtfe(ik) + pqnfe(ik)
                              pihtfe(ik) = pihtfe(ik) + phnfe(ik)
 494                       continue
                        endif
                     endif
                  endif
               endif
            endif
         endif
      endif
c
c
      return
99001 format (' ',' in augfil ')
99002 format (' ',2i4,e12.4)
99003 format (' ',2i4,2e12.4)
      end
      subroutine aughlp(apirt,am,aheat,aion,pirt,esum,et,flem,nzmin,
     &                  nzmax,nzmaxo)
c
c
c
c     this routine does the work in computing auger rates
c
      external bigdat,newdat,rr3,nfllns
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
c
      dimension apirt(nzmaxo,100),aheat(28),aion(28),esum(28),
     &          et(28),am(nzmax,100),pirt(28),flem(28)
c
       character*72 ktitle
c
      data ergsev/1.602197e-12/
      data eps/1./
c
      if ( lpri.gt.2 ) write (6,*) 'in aughlp',nzmin,nzmax,nzmaxo
      do 100 ik = nzmin,nzmax
         nzmm1 = nzmax + 1
         ikp1 = ik + 1
         do 50 jk = ikp1,nzmm1
            apirt(ik,jk) = apirt(ik,jk) + am(jk-1,ik)*pirt(ik)
            if ( lpri.gt.2 ) write (6,*) ik,jk,am(jk-1,ik),
     &                              pirt(ik),apirt(ik,jk)
 50      continue
         if ( ik.ne.nzmax ) then
            sum = 0.
            svm = 0.
            etmp = -esum(ik) + et(ik)
            flem(ik) = 0.
            etst=etmp+esum(ik+1)
c            if (etst.gt.eps ) 
            flem(ik) = am(max0(1,ik-1),ik)
     &           *pirt(ik)
            if (lpri.gt.2) write (6,*)ik,esum(ik+1),etmp,etst,am(ik,ik),
     $          flem(ik)
            ikp2 = ik + 2
            do 60 jk = ikp2,nzmm1
               ak = float(jk-ik) - 1.
               empt = (etmp+esum(jk-1))/ak
               empt = amax1(empt,1.e-4)
               svm = svm + am(jk-1,ik)*ak*phi(empt)
               sum = sum + am(jk-1,ik)*ak*eheat(empt)
 60         continue
            aheat(ik) = sum*pirt(ik)*ergsev
            aion(ik) = svm*pirt(ik)
         endif
 100  continue
c
      return
      end
      subroutine autoi
c
c
c
c
c     this routine computes autoionization
c
c
      external bigdat,newdat,rr3,nfllns
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
      common /airate/ airt(nni)
      common /aidata/ cai(nni),eai(nni)
      common /temp  / t,to
      common /ethrsh/ eth(nni)
      common /abion / xii(nnip),xiin(nni),xiip(nni),xiln(nnnl),
     &                xilp(nnnl),xiio(nnip)
      common /icc   / lichk(nni),lipin
c
      character*72 ktitle
c
      data ergsev/1.602197e-12/
c
      if (lpri.gt.2) write (6,*)'in autoi:',t
c
      tsq = sqrt(t)
      ekt = t*(0.861707)
c
      do 100 j = 1,nni
         airt(j) = 0.
         if ( lichk(j).eq.1 ) then
            airt(j) = cai(j)*expo(-eai(j)/ekt)/tsq
            if (lpri.gt.2) write (6,*)j,cai(j),eai(j),airt(j)
         endif
 100  continue
c
      if (lpri.gt.2) write (6,*)'finishing autoi'
c
      return
      end
      subroutine benchmrk(nbmk)
c
      external bigdat,newdat,rr3,nfllns
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
      common tp,xlum,enlum,ecut,vturb,lpri,lwri,llte,
     &   lfast,lcdd,lrc,lfix,lthin,lbcase,nel,na,
     &   nnnl2,nni2,nnip2,lnoinwd,lffst,ktitle
      common /linsel/ nlsv(nnnl),nlsvn
      common /etot  / etotc,etotco,etotl,etotlo
     $               ,etotc1,etotc2,etotc3,etotl1,etotl2
      common /prs   / p,p0
      common /temp  / t,to
      common /sigh  / zeta
      common /nmrc  / numrec,npass
      common /tlim  / tinf
      common /cdpth / dpthc(ncn)
      common /bdpth / dpthb(ncn)
      common /dpttau/ tauth
      common /tau0ln/ tau0(nnnl)
      common /xcol  / xcc(183)
      common /icc   / lichk(nni),lipin
      common /prtop / elnprnt(400),elimdb(2),nstpt,nlnprnt,lstpt(20)
      common /enerc / epi(ncn),dele(ncn),numcon
      common /spectc/ bremsa(ncn),brems(ncn),zrems(ncn),
     &                zremso(ncn),zremsz(ncn)
      common /radius/delr,r,rl,rmax,rdel,radexp,rscale,rsave
      common /heato / httoto,cltoto,hmctoto
      common /pcool / cll(nni),clbr,clcmp
      common /epden /xee,xpx,xnx,xpxo,xnxo,xpxsave,colsave
      common /abion / xii(nnip),xiin(nni),xiip(nni),xiln(nnnl),
     &                xilp(nnnl),xiio(nnip)
      common /abel  / xel(nl),xeln(nni),xelln(nnnl)
      common /itdat / enfmxs,ensfrc,critd,crittd,epss,crits,
     $             kmaxs,lppris,nlimd,lpprid,nnmax
      common /llumin/ elum(nnnl),oelum(nnnl)
      common /llumnb/ elumb(nnnl),oelmb(nnnl)
      common /lindat/ elin(nnnl),blin(nnnl), 
     &                sxlin(nnnl),eex(nnnl),cslin(nnnl),
     &      flin1(nnnl),flin2(nnnl),nblin(nnnl),nilin(nnnl),nlin(nni)
      common /spctcb/ zremsb(ncn),zrmsbo(ncn),bremsb(ncn),
     &                brmsab(ncn)
c
      character*72 ktitle
c
      dimension bmdat(15,50,10),elnsv(2,50,10),eldat(50,10),
     $   nbmdat(10),fbsv(50),sig(50),err(50),sigsm(50),nsm(50),
     $   nlnsv(10),nbmpt(50,10),nlnsv2(10),fluxrl(50),lllsv(50)
c
      character*8 kinam(169)
      character*8 klbm(50,10)
c
c     'Table3MeudonHIIRegion,
c     ',,Meu,Lex,1,2,3,4,5,6,8,9,10,11
       data nlnsv2(1)/13/
c     Table 2: Cool H II region:  Mean,1,2,3,5,6,8,9,10,11,0.,0./
       data nbmdat(1)/10/
       data klbm(14,1),eldat(14,1),bmdat(1,14,1),bmdat(2,14,1),
     $ bmdat(3,14,1),bmdat(4,14,1),bmdat(5,14,1),bmdat(6,14,1),
     $ bmdat(7,14,1),bmdat(8,14,1),bmdat(9,14,1),bmdat(10,14,1),
     $ bmdat(11,14,1),bmdat(12,14,1)
     $ /'L(Hp)',0.,4.93,4.99,4.98,4.93,4.85,
     $  4.83,4.93,4.94,5.01,4.91,0.,0./
       data klbm(2,1),eldat(2,1),bmdat(1,2,1),bmdat(2,2,1),
     $ bmdat(3,2,1),bmdat(4,2,1),bmdat(5,2,1),bmdat(6,2,1),
     $ bmdat(7,2,1),bmdat(8,2,1),bmdat(9,2,1),bmdat(10,2,1),
     $ bmdat(11,2,1),bmdat(12,2,1)
     $ /'[N II]',6584,0.85,0.82,0.91,0.82,0.97,0.82,
     $  0.84,0.84,0.83,0.83,0.,0./
       data klbm(3,1),eldat(3,1),bmdat(1,3,1),bmdat(2,3,1),
     $ bmdat(3,3,1),bmdat(4,3,1),bmdat(5,3,1),bmdat(6,3,1),
     $ bmdat(7,3,1),bmdat(8,3,1),bmdat(9,3,1),bmdat(10,3,1),
     $ bmdat(11,3,1),bmdat(12,3,1)
     $ /'[O III]',3727,1.18,1.11,1.16,1.22,1.32,1.14,
     $  1.21,1.24,1.14,1.11,0.,0./
       data klbm(4,1),eldat(4,1),bmdat(1,4,1),bmdat(2,4,1),
     $ bmdat(3,4,1),bmdat(4,4,1),bmdat(5,4,1),bmdat(6,4,1),
     $ bmdat(7,4,1),bmdat(8,4,1),bmdat(9,4,1),bmdat(10,4,1),
     $ bmdat(11,4,1),bmdat(12,4,1)
     $ /'[Ne II]',1.28E+05,0.31,0.36,0.35,0.29,0.29,0.29,0.29,
     $  0.35,0.29,0.29,0.,0./
       data klbm(5,1),eldat(5,1),bmdat(1,5,1),bmdat(2,5,1),
     $ bmdat(3,5,1),bmdat(4,5,1),bmdat(5,5,1),bmdat(6,5,1),
     $ bmdat(7,5,1),bmdat(8,5,1),bmdat(9,5,1),bmdat(10,5,1),
     $ bmdat(11,5,1),bmdat(12,5,1)
     $ /'[SII]',6720,0.57,0.69,0.64,0.55,0.61,0.52,0.52,
     $  0.6,0.45,0.58,0.,0./
       data klbm(6,1),eldat(6,1),bmdat(1,6,1),bmdat(2,6,1),
     $ bmdat(3,6,1),bmdat(4,6,1),bmdat(5,6,1),bmdat(6,6,1),
     $ bmdat(7,6,1),bmdat(8,6,1),bmdat(9,6,1),bmdat(10,6,1),
     $ bmdat(11,6,1),bmdat(12,6,1)
     $ /'[S III]',1.87E+05,0.32,0.26,0.27,0.36,0.17,0.37,
     $  0.37,0.33,0.4,0.3,0.,0./
       data klbm(7,1),eldat(7,1),bmdat(1,7,1),bmdat(2,7,1),
     $ bmdat(3,7,1),bmdat(4,7,1),bmdat(5,7,1),bmdat(6,7,1),
     $ bmdat(7,7,1),bmdat(8,7,1),bmdat(9,7,1),bmdat(10,7,1),
     $ bmdat(11,7,1),bmdat(12,7,1)
     $ /'[S III]',3.40E+05,0.52,0.43,0.47,0.6,0.27,0.61,
     $  0.62,0.54,0.67,0.51,0.,0./
       data klbm(8,1),eldat(8,1),bmdat(1,8,1),bmdat(2,8,1),
     $ bmdat(3,8,1),bmdat(4,8,1),bmdat(5,8,1),bmdat(6,8,1),
     $ bmdat(7,8,1),bmdat(8,8,1),bmdat(9,8,1),bmdat(10,8,1),
     $ bmdat(11,8,1),bmdat(12,8,1)
     $ /'[S III]',9532,0.55,0.4,0.48,0.55,0.64,0.6,0.56,
     $  0.49,0.62,0.58,0.,0./
       data klbm(9,1),eldat(9,1),bmdat(1,9,1),bmdat(2,9,1),
     $ bmdat(3,9,1),bmdat(4,9,1),bmdat(5,9,1),bmdat(6,9,1),
     $ bmdat(7,9,1),bmdat(8,9,1),bmdat(9,9,1),bmdat(10,9,1),
     $ bmdat(11,9,1),bmdat(12,9,1)
     $ /'L(total)',0.,21.2,20.3,21.3,21.7,20.7,21,21.8,
     $  21.7,22.1,20.6,0.,0./
       data klbm(10,1),eldat(10,1),bmdat(1,10,1),bmdat(2,10,1),
     $ bmdat(3,10,1),bmdat(4,10,1),bmdat(5,10,1),bmdat(6,10,1),
     $ bmdat(7,10,1),bmdat(8,10,1),bmdat(9,10,1),bmdat(10,10,1),
     $ bmdat(11,10,1),bmdat(12,10,1)
     $ /'T(in)',0.,6793,6860,6952,6749,6980,6870,6747,
     $  6230,6912,6838,0.,0./
       data klbm(11,1),eldat(11,1),bmdat(1,11,1),bmdat(2,11,1),
     $ bmdat(3,11,1),bmdat(4,11,1),bmdat(5,11,1),bmdat(6,11,1),
     $ bmdat(7,11,1),bmdat(8,11,1),bmdat(9,11,1),bmdat(10,11,1),
     $ bmdat(11,11,1),bmdat(12,11,1)
     $ /'T(H+)',0.,6744,6690,6740,6742,6950,6660,6742,
     $  6770,6720,6681,0.,0./
       data klbm(12,1),eldat(12,1),bmdat(1,12,1),bmdat(2,12,1),
     $ bmdat(3,12,1),bmdat(4,12,1),bmdat(5,12,1),bmdat(6,12,1),
     $ bmdat(7,12,1),bmdat(8,12,1),bmdat(9,12,1),bmdat(10,12,1),
     $ bmdat(11,12,1),bmdat(12,12,1)
     $ /'<He+>/<H+>',0.,0.054,0,0.041,0.044,0.068,0.048,
     $  0.034,0.055,0.09,0.,0.,0./
       data klbm(13,1),eldat(13,1),bmdat(1,13,1),bmdat(2,13,1),
     $ bmdat(3,13,1),bmdat(4,13,1),bmdat(5,13,1),bmdat(6,13,1),
     $ bmdat(7,13,1),bmdat(8,13,1),bmdat(9,13,1),bmdat(10,13,1),
     $ bmdat(11,13,1),bmdat(12,13,1)
     $ /'R(out) E18',0.,8.96,9,8.93,8.94,9,8.93,9,8.87,9.02,8.97,0.,0./
c    Table 3 Meudon HII Region, Meu,Lex,1,2,3,4,5,6,8,9,10,11/
       data nbmdat(2)/12/
       data nlnsv2(2)/23/
       data klbm(24,2),eldat(24,2),bmdat(1,24,2),bmdat(2,24,2),     
     $ bmdat(3,24,2),bmdat(4,24,2),bmdat(5,24,2),bmdat(6,24,2),     
     $ bmdat(7,24,2),bmdat(8,24,2),bmdat(9,24,2),bmdat(10,24,2),     
     $ bmdat(11,24,2),bmdat(12,24,2)/     
     $ 'H I',4862.,2.06,2.03,1.96,2.06,2.04,1.86,2.02,2.02,
     $2.05,2.10,2.11,2.09/
       data klbm(1,2),eldat(1,2),bmdat(1,1,2),bmdat(2,1,2),
     $ bmdat(3,1,2),bmdat(4,1,2),bmdat(5,1,2),bmdat(6,1,2),
     $ bmdat(7,1,2),bmdat(8,1,2),bmdat(9,1,2),bmdat(10,1,2),
     $ bmdat(11,1,2),bmdat(12,1,2)
     $ /'He I',5876,0.116,0.116,0.125,0.109,
     $ 0.119,0.11,0.101,0.116,0.,0.125,0.115,0.12/
       data klbm(2,2),eldat(2,2),bmdat(1,2,2),bmdat(2,2,2),
     $ bmdat(3,2,2),bmdat(4,2,2),bmdat(5,2,2),bmdat(6,2,2),
     $ bmdat(7,2,2),bmdat(8,2,2),bmdat(9,2,2),bmdat(10,2,2),
     $ bmdat(11,2,2),bmdat(12,2,2)
     $ /'Cll',2326,0.17,0.16,0.07,0.19,0.17,0.16,0.16,0.14,
     $  0.18,0.28,0.12,0.14/
       data klbm(3,2),eldat(3,2),bmdat(1,3,2),bmdat(2,3,2),
     $ bmdat(3,3,2),bmdat(4,3,2),bmdat(5,3,2),bmdat(6,3,2),
     $ bmdat(7,3,2),bmdat(8,3,2),bmdat(9,3,2),bmdat(10,3,2),
     $ bmdat(11,3,2),bmdat(12,3,2)
     $ /'C',1909,0.051,0.06,0.05,0.059,0.059,
     $0.027,0.078,0.065,0.076,0.082,0.077,0.071/
       data klbm(4,2),eldat(4,2),bmdat(1,4,2),bmdat(2,4,2),
     $ bmdat(3,4,2),bmdat(4,4,2),bmdat(5,4,2),bmdat(6,4,2),
     $ bmdat(7,4,2),bmdat(8,4,2),bmdat(9,4,2),bmdat(10,4,2),
     $ bmdat(11,4,2),bmdat(12,4,2)
     $ /'[N II]',1.22E+06,0.,0.031,0.032,0.033,0.,0.,0.,0.036,
     $  0.031,0.03,0.037,0.034/
       data klbm(5,2),eldat(5,2),bmdat(1,5,2),bmdat(2,5,2),
     $ bmdat(3,5,2),bmdat(4,5,2),bmdat(5,5,2),bmdat(6,5,2),
     $ bmdat(7,5,2),bmdat(8,5,2),bmdat(9,5,2),bmdat(10,5,2),
     $ bmdat(11,5,2),bmdat(12,5,2)
     $ /'[N II]',6584,0.73,0.79,0.61,0.88,0.74,0.94,0.87,0.78,
     $  0.73,0.78,0.81,0.75/
       data klbm(6,2),eldat(6,2),bmdat(1,6,2),bmdat(2,6,2),
     $ bmdat(3,6,2),bmdat(4,6,2),bmdat(5,6,2),bmdat(6,6,2),
     $ bmdat(7,6,2),bmdat(8,6,2),bmdat(9,6,2),bmdat(10,6,2),
     $ bmdat(11,6,2),bmdat(12,6,2)
     $ /'[N III]',5.70E+05,0.3,0.27,0.16,0.27,0.29,0.,0.26,
     $  0.3,0.3,0.17,0.27,0.39/
       data klbm(7,2),eldat(7,2),bmdat(1,7,2),bmdat(2,7,2),
     $ bmdat(3,7,2),bmdat(4,7,2),bmdat(5,7,2),bmdat(6,7,2),
     $ bmdat(7,7,2),bmdat(8,7,2),bmdat(9,7,2),bmdat(10,7,2),
     $ bmdat(11,7,2),bmdat(12,7,2)
     $ /'[0 II]',3727,2.01,2.16,1.5,2.19,2.14,2.56,2.3,2.11,
     $  2.26,2.41,2.2,1.95/
       data klbm(8,2),eldat(8,2),bmdat(1,8,2),bmdat(2,8,2),
     $ bmdat(3,8,2),bmdat(4,8,2),bmdat(5,8,2),bmdat(6,8,2),
     $ bmdat(7,8,2),bmdat(8,8,2),bmdat(9,8,2),bmdat(10,8,2),
     $ bmdat(11,8,2),bmdat(12,8,2)
     $ /'[0 III]',5.18E+05,1.1,1.07,1.1,1.04,1.11,1.04,0.99,
     $  1.08,1.08,1.23,1.04,0.97/
       data klbm(9,2),eldat(9,2),bmdat(1,9,2),bmdat(2,9,2),
     $ bmdat(3,9,2),bmdat(4,9,2),bmdat(5,9,2),bmdat(6,9,2),
     $ bmdat(7,9,2),bmdat(8,9,2),bmdat(9,9,2),bmdat(10,9,2),
     $ bmdat(11,9,2),bmdat(12,9,2)
     $ /'[0 III]',8.84E+05,1.2,1.23,1.3,1.07,1.28,0.,1.16,
     $  1.25,1.26,1.42,1.2,1.14/
       data klbm(10,2),eldat(10,2),bmdat(1,10,2),bmdat(2,10,2),
     $ bmdat(3,10,2),bmdat(4,10,2),bmdat(5,10,2),bmdat(6,10,2),
     $ bmdat(7,10,2),bmdat(8,10,2),bmdat(9,10,2),bmdat(10,10,2),
     $ bmdat(11,10,2),bmdat(12,10,2)
     $ /'[0 III]',5007,2.03,2.06,2.3,1.93,1.96,1.47,2.29,2.17,
     $  2.1,2.23,2.22,1.89/
       data klbm(11,2),eldat(11,2),bmdat(1,11,2),bmdat(2,11,2),
     $ bmdat(3,11,2),bmdat(4,11,2),bmdat(5,11,2),bmdat(6,11,2),
     $ bmdat(7,11,2),bmdat(8,11,2),bmdat(9,11,2),bmdat(10,11,2),
     $ bmdat(11,11,2),bmdat(12,11,2)
     $ /'[N III]',1.28E+05,0.21,0.22,0.26,0.23,0.19,0.23,0.22,
     $  0.2,0.2,0.22,0.22,0.2/
       data klbm(12,2),eldat(12,2),bmdat(1,12,2),bmdat(2,12,2),
     $ bmdat(3,12,2),bmdat(4,12,2),bmdat(5,12,2),bmdat(6,12,2),
     $ bmdat(7,12,2),bmdat(8,12,2),bmdat(9,12,2),bmdat(10,12,2),
     $ bmdat(11,12,2),bmdat(12,12,2)
     $ /'[Ne III]',1.55E+05,0.44,0.38,0.37,0.43,0.43,0.47,0.37,
     $  0.42,0.42,0.22,0.34,0.38/
       data klbm(13,2),eldat(13,2),bmdat(1,13,2),bmdat(2,13,2),
     $ bmdat(3,13,2),bmdat(4,13,2),bmdat(5,13,2),bmdat(6,13,2),
     $ bmdat(7,13,2),bmdat(8,13,2),bmdat(9,13,2),bmdat(10,13,2),
     $ bmdat(11,13,2),bmdat(12,13,2)
     $ /'[Ne III]',3869,0.096,0.086,0.085,0.103,
     $  0.086,0.071,0.1,0.079,0.087,0.081,0.087,0.078/
       data klbm(14,2),eldat(14,2),bmdat(1,14,2),bmdat(2,14,2),
     $ bmdat(3,14,2),bmdat(4,14,2),bmdat(5,14,2),bmdat(6,14,2),
     $ bmdat(7,14,2),bmdat(8,14,2),bmdat(9,14,2),bmdat(10,14,2),
     $ bmdat(11,14,2),bmdat(12,14,2)
     $ /'[S II]',6720,0.14,0.2,0.24,0.23,0.16,0.25,0.22,0.17,
     $  0.13,0.21,0.15,0.21/
       data klbm(15,2),eldat(15,2),bmdat(1,15,2),bmdat(2,15,2),
     $ bmdat(3,15,2),bmdat(4,15,2),bmdat(5,15,2),bmdat(6,15,2),
     $ bmdat(7,15,2),bmdat(8,15,2),bmdat(9,15,2),bmdat(10,15,2),
     $ bmdat(11,15,2),bmdat(12,15,2)
     $ /'[S III]',1.87E+05,0.55,0.55,0.56,0.48,0.56,0.53,0.5,0.55,
     $  0.58,0.58,0.58,0.55/
       data klbm(16,2),eldat(16,2),bmdat(1,16,2),bmdat(2,16,2),
     $ bmdat(3,16,2),bmdat(4,16,2),bmdat(5,16,2),bmdat(6,16,2),
     $ bmdat(7,16,2),bmdat(8,16,2),bmdat(9,16,2),bmdat(10,16,2),
     $ bmdat(11,16,2),bmdat(12,16,2)
     $ /'[S III]',3.40E+05,0.93,0.89,0.91,0.82,0.89,0.,0.81,
     $  0.88,0.94,0.92,0.92,0.91/
       data klbm(17,2),eldat(17,2),bmdat(1,17,2),bmdat(2,17,2),
     $ bmdat(3,17,2),bmdat(4,17,2),bmdat(5,17,2),bmdat(6,17,2),
     $ bmdat(7,17,2),bmdat(8,17,2),bmdat(9,17,2),bmdat(10,17,2),
     $ bmdat(11,17,2),bmdat(12,17,2)
     $ /'[S III]',9532,1.25,1.29,1.16,1.27,1.23,1.15,1.48,1.27,
     $  1.3,1.31,1.32,1.46/
       data klbm(18,2),eldat(18,2),bmdat(1,18,2),bmdat(2,18,2),
     $ bmdat(3,18,2),bmdat(4,18,2),bmdat(5,18,2),bmdat(6,18,2),
     $ bmdat(7,18,2),bmdat(8,18,2),bmdat(9,18,2),bmdat(10,18,2),
     $ bmdat(11,18,2),bmdat(12,18,2)
     $ /'[SIV]',1.05E+05,0.39,0.34,0.22,0.37,0.42,0.35,0.36,
     $  0.41,0.33,0.26,0.38,0.27/
       data klbm(19,2),eldat(19,2),bmdat(1,19,2),bmdat(2,19,2),
     $ bmdat(3,19,2),bmdat(4,19,2),bmdat(5,19,2),bmdat(6,19,2),
     $ bmdat(7,19,2),bmdat(8,19,2),bmdat(9,19,2),bmdat(10,19,2),
     $ bmdat(11,19,2),bmdat(12,19,2)
     $ /'L(total)',0.,24.1,24.2,21.7,24.1,24.1,17.4,24.8,24.3,
     $  24.6,26.4,25.5,24.1/
       data klbm(20,2),eldat(20,2),bmdat(1,20,2),bmdat(2,20,2),
     $ bmdat(3,20,2),bmdat(4,20,2),bmdat(5,20,2),bmdat(6,20,2),
     $ bmdat(7,20,2),bmdat(8,20,2),bmdat(9,20,2),bmdat(10,20,2),
     $ bmdat(11,20,2),bmdat(12,20,2)
     $ /'T(in)',0.,7992,7552,7630,7815,7741,
     $  8057,7670,7650,7399,6530,7582,7445/
       data klbm(21,2),eldat(21,2),bmdat(1,21,2),bmdat(2,21,2),
     $ bmdat(3,21,2),bmdat(4,21,2),bmdat(5,21,2),bmdat(6,21,2),
     $ bmdat(7,21,2),bmdat(8,21,2),bmdat(9,21,2),bmdat(10,21,2),
     $ bmdat(11,21,2),bmdat(12,21,2)
     $ /'T(H+)',0.,7378,8034,7880,8064,8047,
     $  7879,8000,8060,8087,8220,8191,7913/
       data klbm(22,2),eldat(22,2),bmdat(1,22,2),bmdat(2,22,2),
     $ bmdat(3,22,2),bmdat(4,22,2),bmdat(5,22,2),bmdat(6,22,2),
     $ bmdat(7,22,2),bmdat(8,22,2),bmdat(9,22,2),bmdat(10,22,2),
     $ bmdat(11,22,2),bmdat(12,22,2)
     $ /'<He+>/<H+>',0.,0,0.77,0,0.71,0.77,
     $  0.69,0.76,0.75,0.83,0.86,0.79,0.77/
       data klbm(23,2),eldat(23,2),bmdat(1,23,2),bmdat(2,23,2),
     $ bmdat(3,23,2),bmdat(4,23,2),bmdat(5,23,2),bmdat(6,23,2),
     $ bmdat(7,23,2),bmdat(8,23,2),bmdat(9,23,2),bmdat(10,23,2),
     $ bmdat(11,23,2),bmdat(12,23,2)
     $ /'R(out)',0.,1.45,1.48,1.43,1.46,1.46,1.61,1.47,1.46,
     $  1.46,1.46,1.49,1.47/
c     Table 4 Blister HII Region   Mean,1,2,3,4,5,6,8,9,10,11,0./
       data nbmdat(3)/11/
       data nlnsv2(3)/22/
       data klbm(1,3),eldat(1,3),bmdat(1,1,3),bmdat(2,1,3),
     $ bmdat(3,1,3),bmdat(4,1,3),bmdat(5,1,3),bmdat(6,1,3),
     $ bmdat(7,1,3),bmdat(8,1,3),bmdat(9,1,3),bmdat(10,1,3),
     $ bmdat(11,1,3),bmdat(12,1,3)
     $ /' I(Hp)',0.,4.62,4.6,4.59,4.81,3.89,
     $  4.69,4.67,4.7,4.85,4.58,4.78,0./
       data klbm(2,3),eldat(2,3),bmdat(1,2,3),bmdat(2,2,3),
     $ bmdat(3,2,3),bmdat(4,2,3),bmdat(5,2,3),bmdat(6,2,3),
     $ bmdat(7,2,3),bmdat(8,2,3),bmdat(9,2,3),bmdat(10,2,3),
     $ bmdat(11,2,3),bmdat(12,2,3)
     $ /'He I',5876,0.12,0.12,0.13,0.11,0.11,
     $  0.12,0.12,0.,0.12,0.11,0.12,0./
       data klbm(3,3),eldat(3,3),bmdat(1,3,3),bmdat(2,3,3),
     $ bmdat(3,3,3),bmdat(4,3,3),bmdat(5,3,3),bmdat(6,3,3),
     $ bmdat(7,3,3),bmdat(8,3,3),bmdat(9,3,3),bmdat(10,3,3),
     $ bmdat(11,3,3),bmdat(12,3,3)
     $ /'Cll',2326,0.18,0.06,0.14,0.2,0.3,
     $  0.1,0.15,0.23,0.35,0.11,0.16,0./
       data klbm(4,3),eldat(4,3),bmdat(1,4,3),bmdat(2,4,3),
     $ bmdat(3,4,3),bmdat(4,4,3),bmdat(5,4,3),bmdat(6,4,3),
     $ bmdat(7,4,3),bmdat(8,4,3),bmdat(9,4,3),bmdat(10,4,3),
     $ bmdat(11,4,3),bmdat(12,4,3)
     $ /'Cll',1335,0.09,0.002,0.17,0.14,0.02,
     $  0.13,0.16,0.,0.01,0.02,0.13,0./
       data klbm(5,3),eldat(5,3),bmdat(1,5,3),bmdat(2,5,3),
     $ bmdat(3,5,3),bmdat(4,5,3),bmdat(5,5,3),bmdat(6,5,3),
     $ bmdat(7,5,3),bmdat(8,5,3),bmdat(9,5,3),bmdat(10,5,3),
     $ bmdat(11,5,3),bmdat(12,5,3)
     $ /'Clil]',1909,0.17,0.13,0.22,0.17,
     $  0.08,0.18,0.15,0.2,0.25,0.14,0.23,0./
       data klbm(6,3),eldat(6,3),bmdat(1,6,3),bmdat(2,6,3),
     $ bmdat(3,6,3),bmdat(4,6,3),bmdat(5,6,3),bmdat(6,6,3),
     $ bmdat(7,6,3),bmdat(8,6,3),bmdat(9,6,3),bmdat(10,6,3),
     $ bmdat(11,6,3),bmdat(12,6,3)
     $ /'[N II]',6584,0.87,0.67,0.58,0.94,
     $  1.48,0.74,0.9,0.87,0.92,0.82,0.83,0./
       data klbm(7,3),eldat(7,3),bmdat(1,7,3),bmdat(2,7,3),
     $ bmdat(3,7,3),bmdat(4,7,3),bmdat(5,7,3),bmdat(6,7,3),
     $ bmdat(7,7,3),bmdat(8,7,3),bmdat(9,7,3),bmdat(10,7,3),
     $ bmdat(11,7,3),bmdat(12,7,3)
     $ /'[Nlil]',5.70E+04,0.031,0.032,0.035,0.033,0.,
     $  0.033,0.032,0.034,0.014,0.032,0.033,0./
       data klbm(8,3),eldat(8,3),bmdat(1,8,3),bmdat(2,8,3),
     $ bmdat(3,8,3),bmdat(4,8,3),bmdat(5,8,3),bmdat(6,8,3),
     $ bmdat(7,8,3),bmdat(8,8,3),bmdat(9,8,3),bmdat(10,8,3),
     $ bmdat(11,8,3),bmdat(12,8,3)
     $ /'[O III]',7330,0.12,0.06,0.1,0.13,0.19,
     $  0.09,0.12,0.14,0.15,0.08,0.1,0./
       data klbm(9,3),eldat(9,3),bmdat(1,9,3),bmdat(2,9,3),
     $ bmdat(3,9,3),bmdat(4,9,3),bmdat(5,9,3),bmdat(6,9,3),
     $ bmdat(7,9,3),bmdat(8,9,3),bmdat(9,9,3),bmdat(10,9,3),
     $ bmdat(11,9,3),bmdat(12,9,3)
     $ /'[O III]',3727,0.88,0.53,0.73,0.98,
     $  1.39,0.69,0.86,1.04,1.04,0.73,0.86,0./
       data klbm(10,3),eldat(10,3),bmdat(1,10,3),bmdat(2,10,3),
     $ bmdat(3,10,3),bmdat(4,10,3),bmdat(5,10,3),bmdat(6,10,3),
     $ bmdat(7,10,3),bmdat(8,10,3),bmdat(9,10,3),bmdat(10,10,3),
     $ bmdat(11,10,3),bmdat(12,10,3)
     $ /'[Olil]',5.18E+04,0.29,0.29,0.31,
     $  0.29,0.26,0.28,0.28,0.28,0.32,0.28,0.27,0./
       data klbm(11,3),eldat(11,3),bmdat(1,11,3),bmdat(2,11,3),
     $ bmdat(3,11,3),bmdat(4,11,3),bmdat(5,11,3),bmdat(6,11,3),
     $ bmdat(7,11,3),bmdat(8,11,3),bmdat(9,11,3),bmdat(10,11,3),
     $ bmdat(11,11,3),bmdat(12,11,3)
     $ /'[Olil]',5007,4.13,4.5,4.74,3.9,
     $  3.28,4.4,3.9,3.96,4.51,4.16,3.98,0./
       data klbm(12,3),eldat(12,3),bmdat(1,12,3),bmdat(2,12,3),
     $ bmdat(3,12,3),bmdat(4,12,3),bmdat(5,12,3),bmdat(6,12,3),
     $ bmdat(7,12,3),bmdat(8,12,3),bmdat(9,12,3),bmdat(10,12,3),
     $ bmdat(11,12,3),bmdat(12,12,3)
     $ /'[N III]',1.28E+04,0.36,0.45,
     $  0.32,0.33,0.44,0.35,0.33,0.35,0.36,0.37,0.35,0./
       data klbm(13,3),eldat(13,3),bmdat(1,13,3),bmdat(2,13,3),
     $ bmdat(3,13,3),bmdat(4,13,3),bmdat(5,13,3),bmdat(6,13,3),
     $ bmdat(7,13,3),bmdat(8,13,3),bmdat(9,13,3),bmdat(10,13,3),
     $ bmdat(11,13,3),bmdat(12,13,3)
     $ /'[Nelil]',1.55E+04,0.98,0.93,
     $  1.24,1.07,1.09,0.96,1.04,1,0.59,0.92,0.97,0./
       data klbm(14,3),eldat(14,3),bmdat(1,14,3),bmdat(2,14,3),
     $ bmdat(3,14,3),bmdat(4,14,3),bmdat(5,14,3),bmdat(6,14,3),
     $ bmdat(7,14,3),bmdat(8,14,3),bmdat(9,14,3),bmdat(10,14,3),
     $ bmdat(11,14,3),bmdat(12,14,3)
     $ /'[Nelil]',3869,0.33,0.33,0.48,0.32,
     $  0.31,0.35,0.26,0.29,35,0.31,0.31,0./
       data klbm(15,3),eldat(15,3),bmdat(1,15,3),bmdat(2,15,3),
     $ bmdat(3,15,3),bmdat(4,15,3),bmdat(5,15,3),bmdat(6,15,3),
     $ bmdat(7,15,3),bmdat(8,15,3),bmdat(9,15,3),bmdat(10,15,3),
     $ bmdat(11,15,3),bmdat(12,15,3)
     $ /'[S III]',1.87E+04,0.35,0.37,0.31,
     $  0.34,0.37,0.31,0.33,0.35,0.37,0.34,0.39,0./
       data klbm(16,3),eldat(16,3),bmdat(1,16,3),bmdat(2,16,3),
     $ bmdat(3,16,3),bmdat(4,16,3),bmdat(5,16,3),bmdat(6,16,3),
     $ bmdat(7,16,3),bmdat(8,16,3),bmdat(9,16,3),bmdat(10,16,3),
     $ bmdat(11,16,3),bmdat(12,16,3)
     $ /'[S III]',9532,1.53,1.52,1.41,1.46,
     $  1.62,1.51,1.42,1.53,1.61,1.42,1.82,0./
       data klbm(17,3),eldat(17,3),bmdat(1,17,3),bmdat(2,17,3),
     $ bmdat(3,17,3),bmdat(4,17,3),bmdat(5,17,3),bmdat(6,17,3),
     $ bmdat(7,17,3),bmdat(8,17,3),bmdat(9,17,3),bmdat(10,17,3),
     $ bmdat(11,17,3),bmdat(12,17,3)
     $ /'[SIV1',1.05E+04,0.46,0.26,0.54,
     $  0.52,0.42,0.51,0.53,0.43,0.36,0.5,0.51,0./
       data klbm(18,3),eldat(18,3),bmdat(1,18,3),bmdat(2,18,3),
     $ bmdat(3,18,3),bmdat(4,18,3),bmdat(5,18,3),bmdat(6,18,3),
     $ bmdat(7,18,3),bmdat(8,18,3),bmdat(9,18,3),bmdat(10,18,3),
     $ bmdat(11,18,3),bmdat(12,18,3)
     $ /'I(total)',0.,50.3,47.1,52.6,52.4,
     $  44.2,50.4,49.4,50.3,54.9,47.4,52.9,0./
       data klbm(19,3),eldat(19,3),bmdat(1,19,3),bmdat(2,19,3),
     $ bmdat(3,19,3),bmdat(4,19,3),bmdat(5,19,3),bmdat(6,19,3),
     $ bmdat(7,19,3),bmdat(8,19,3),bmdat(9,19,3),bmdat(10,19,3),
     $ bmdat(11,19,3),bmdat(12,19,3)
     $ /'T(in)',0.,7989,8300,8206,7582,
     $  0.,8200,8200,7366,7740,8122,8189,0./
       data klbm(20,3),eldat(20,3),bmdat(1,20,3),bmdat(2,20,3),
     $ bmdat(3,20,3),bmdat(4,20,3),bmdat(5,20,3),bmdat(6,20,3),
     $ bmdat(7,20,3),bmdat(8,20,3),bmdat(9,20,3),bmdat(10,20,3),
     $ bmdat(11,20,3),bmdat(12,20,3)
     $ /'T(H+)',0.,8263,8170,8324,8351,
     $  0.,8310,8200,8328,8220,8217,8250,0./
       data klbm(21,3),eldat(21,3),bmdat(1,21,3),bmdat(2,21,3),
     $ bmdat(3,21,3),bmdat(4,21,3),bmdat(5,21,3),bmdat(6,21,3),
     $ bmdat(7,21,3),bmdat(8,21,3),bmdat(9,21,3),bmdat(10,21,3),
     $ bmdat(11,21,3),bmdat(12,21,3)
     $ /'<He+>/<H+>',0.,0.85,0.,0.94,0.78,
     $  0.,0.93,0.79,0.84,0.86,0.85,0.84,0./
       data klbm(22,3),eldat(22,3),bmdat(1,22,3),bmdat(2,22,3),
     $ bmdat(3,22,3),bmdat(4,22,3),bmdat(5,22,3),bmdat(6,22,3),
     $ bmdat(7,22,3),bmdat(8,22,3),bmdat(9,22,3),bmdat(10,22,3),
     $ bmdat(11,22,3),bmdat(12,22,3)
     $ /'AR',0.,0.,2.9,2.88,3.08,0.,2.93,
     $           2.98,3.09,3.1,2.67,3.03,0./
c     Table 6 High Ionization PN,  Mean,1,2,3,4,5,6,10,11,0.,0.,0./
       data nbmdat(5)/9/
       data nlnsv2(5)/22/
       data klbm(1,5),eldat(1,5),bmdat(1,1,5),bmdat(2,1,5),
     $ bmdat(3,1,5),bmdat(4,1,5),bmdat(5,1,5),bmdat(6,1,5),
     $ bmdat(7,1,5),bmdat(8,1,5),bmdat(9,1,5),bmdat(10,1,5),
     $ bmdat(11,1,5),bmdat(12,1,5)
     $ /' L(H:)',0.,5.85,5.67,6.05,5.96,6.02,
     $      5.65,5.74,5.72,6.02,0.,0.,0./
       data klbm(2,5),eldat(2,5),bmdat(1,2,5),bmdat(2,2,5),
     $ bmdat(3,2,5),bmdat(4,2,5),bmdat(5,2,5),bmdat(6,2,5),
     $ bmdat(7,2,5),bmdat(8,2,5),bmdat(9,2,5),bmdat(10,2,5),
     $ bmdat(11,2,5),bmdat(12,2,5)
     $ /'He I',5876,0.12,0.12,0.13,0.13,
     $           0.1,0.12,0.13,0.13,0.13,0.,0.,0./
       data klbm(3,5),eldat(3,5),bmdat(1,3,5),bmdat(2,3,5),
     $ bmdat(3,3,5),bmdat(4,3,5),bmdat(5,3,5),bmdat(6,3,5),
     $ bmdat(7,3,5),bmdat(8,3,5),bmdat(9,3,5),bmdat(10,3,5),
     $ bmdat(11,3,5),bmdat(12,3,5)
     $ /'He II',4686,0.081,0.096,0.08,0.087,
     $  0.039,0.085,0.092,0.09,0.083,0.,0.,0./
       data klbm(4,5),eldat(4,5),bmdat(1,4,5),bmdat(2,4,5),
     $ bmdat(3,4,5),bmdat(4,4,5),bmdat(5,4,5),bmdat(6,4,5),
     $ bmdat(7,4,5),bmdat(8,4,5),bmdat(9,4,5),bmdat(10,4,5),
     $ bmdat(11,4,5),bmdat(12,4,5)
     $ /'C III]',1909,0.83,0.9,0.6,0.6,
     $          0.89,0.99,0.89,1.03,0.74,0.,0.,0./
       data klbm(5,5),eldat(5,5),bmdat(1,5,5),bmdat(2,5,5),
     $ bmdat(3,5,5),bmdat(4,5,5),bmdat(5,5,5),bmdat(6,5,5),
     $ bmdat(7,5,5),bmdat(8,5,5),bmdat(9,5,5),bmdat(10,5,5),
     $ bmdat(11,5,5),bmdat(12,5,5)
     $ /'C IV',1549,0.34,0.24,0.35,0.29,0.45,
     $           0.4,0.37,0.32,0.29,0.,0.,0./
       data klbm(6,5),eldat(6,5),bmdat(1,6,5),bmdat(2,6,5),
     $ bmdat(3,6,5),bmdat(4,6,5),bmdat(5,6,5),bmdat(6,6,5),
     $ bmdat(7,6,5),bmdat(8,6,5),bmdat(9,6,5),bmdat(10,6,5),
     $ bmdat(11,6,5),bmdat(12,6,5)
     $ /'[N II]',6584,0.12,0.12,0.11,0.11,
     $  0.14,0.15,0.12,0.12,0.12,0.,0.,0./
       data klbm(7,5),eldat(7,5),bmdat(1,7,5),bmdat(2,7,5),
     $ bmdat(3,7,5),bmdat(4,7,5),bmdat(5,7,5),bmdat(6,7,5),
     $ bmdat(7,7,5),bmdat(8,7,5),bmdat(9,7,5),bmdat(10,7,5),
     $ bmdat(11,7,5),bmdat(12,7,5)
     $ /'[Nlil]',5.70E+04,0.39,0.27,0.37,
     $            0.,0.,0.4,0.41,0.4,0.48,0.,0.,0./
       data klbm(8,5),eldat(8,5),bmdat(1,8,5),bmdat(2,8,5),
     $ bmdat(3,8,5),bmdat(4,8,5),bmdat(5,8,5),bmdat(6,8,5),
     $ bmdat(7,8,5),bmdat(8,8,5),bmdat(9,8,5),bmdat(10,8,5),
     $ bmdat(11,8,5),bmdat(12,8,5)
     $ /'[O II]',3727,0.29,0.32,0.22,0.24,
     $  0.35,0.35,0.26,0.32,0.27,0.,0.,0./
       data klbm(9,5),eldat(9,5),bmdat(1,9,5),bmdat(2,9,5),
     $ bmdat(3,9,5),bmdat(4,9,5),bmdat(5,9,5),bmdat(6,9,5),
     $ bmdat(7,9,5),bmdat(8,9,5),bmdat(9,9,5),bmdat(10,9,5),
     $ bmdat(11,9,5),bmdat(12,9,5)
     $ /'[O lil]',5007,11.5,12.1,10,
     $  10.1,12.7,12.2,11.7,11.9,11.2,0.,0.,0./
       data klbm(10,5),eldat(10,5),bmdat(1,10,5),bmdat(2,10,5),
     $ bmdat(3,10,5),bmdat(4,10,5),bmdat(5,10,5),bmdat(6,10,5),
     $ bmdat(7,10,5),bmdat(8,10,5),bmdat(9,10,5),bmdat(10,10,5),
     $ bmdat(11,10,5),bmdat(12,10,5)
     $ /'[O lil]',5.20E+04,2.02,2.03,1.88,1.96,2.39,
     $   1.95,2.02,2.02,1.94,0.,0.,0./
       data klbm(11,5),eldat(11,5),bmdat(1,11,5),bmdat(2,11,5),
     $ bmdat(3,11,5),bmdat(4,11,5),bmdat(5,11,5),bmdat(6,11,5),
     $ bmdat(7,11,5),bmdat(8,11,5),bmdat(9,11,5),bmdat(10,11,5),
     $ bmdat(11,11,5),bmdat(12,11,5)
     $ /'[OIV]',2.60E+04,0.79,0.76,0.68,
     $  0.8,1.09,0.71,0.86,0.77,0.67,0.,0.,0./
       data klbm(12,5),eldat(12,5),bmdat(1,12,5),bmdat(2,12,5),
     $ bmdat(3,12,5),bmdat(4,12,5),bmdat(5,12,5),bmdat(6,12,5),
     $ bmdat(7,12,5),bmdat(8,12,5),bmdat(9,12,5),bmdat(10,12,5),
     $ bmdat(11,12,5),bmdat(12,12,5)
     $ /'[Ne lil]',1.55E+04,1.35,1.35,1.3,1.32,
     $  1.55,1.3,1.35,1.34,1.31,0.,0.,0./
       data klbm(13,5),eldat(13,5),bmdat(1,13,5),bmdat(2,13,5),
     $ bmdat(3,13,5),bmdat(4,13,5),bmdat(5,13,5),bmdat(6,13,5),
     $ bmdat(7,13,5),bmdat(8,13,5),bmdat(9,13,5),bmdat(10,13,5),
     $ bmdat(11,13,5),bmdat(12,13,5)
     $ /'[Ne lil]',3869,1.03,1.15,1.02,0.92,
     $      1.02,1.13,0.89,1.11,1,0.,0.,0./
       data klbm(14,5),eldat(14,5),bmdat(1,14,5),bmdat(2,14,5),
     $ bmdat(3,14,5),bmdat(4,14,5),bmdat(5,14,5),bmdat(6,14,5),
     $ bmdat(7,14,5),bmdat(8,14,5),bmdat(9,14,5),bmdat(10,14,5),
     $ bmdat(11,14,5),bmdat(12,14,5)
     $ /'Mg II',2798,0.13,0.34,0.1,0.07,0.14,0.05,0.1,0.1,0.11,0.,0.,0./
       data klbm(15,5),eldat(15,5),bmdat(1,15,5),bmdat(2,15,5),
     $ bmdat(3,15,5),bmdat(4,15,5),bmdat(5,15,5),bmdat(6,15,5),
     $ bmdat(7,15,5),bmdat(8,15,5),bmdat(9,15,5),bmdat(10,15,5),
     $ bmdat(11,15,5),bmdat(12,15,5)
     $ /'Si lil]',1892,0.15,0.1,0.09,
     $        0.1,0.37,0.15,0.13,0.,0.11,0.,0.,0./
       data klbm(16,5),eldat(16,5),bmdat(1,16,5),bmdat(2,16,5),
     $ bmdat(3,16,5),bmdat(4,16,5),bmdat(5,16,5),bmdat(6,16,5),
     $ bmdat(7,16,5),bmdat(8,16,5),bmdat(9,16,5),bmdat(10,16,5),
     $ bmdat(11,16,5),bmdat(12,16,5)
     $ /'[SIII]',1.87E+04,0.34,0.45,
     $  0.24,0.32,0.,0.26,0.28,0.36,0.49,0.,0.,0./
       data klbm(17,5),eldat(17,5),bmdat(1,17,5),bmdat(2,17,5),
     $ bmdat(3,17,5),bmdat(4,17,5),bmdat(5,17,5),bmdat(6,17,5),
     $ bmdat(7,17,5),bmdat(8,17,5),bmdat(9,17,5),bmdat(10,17,5),
     $ bmdat(11,17,5),bmdat(12,17,5)
     $ /'[S lil]',9532,1.13,1.4,0.81,
     $                  0.92,0.,1.02,0.85,1.1,1.77,0.,0.,0./
       data klbm(18,5),eldat(18,5),bmdat(1,18,5),bmdat(2,18,5),
     $ bmdat(3,18,5),bmdat(4,18,5),bmdat(5,18,5),bmdat(6,18,5),
     $ bmdat(7,18,5),bmdat(8,18,5),bmdat(9,18,5),bmdat(10,18,5),
     $ bmdat(11,18,5),bmdat(12,18,5)
     $ /'1S IV1',1.05E+04,2.05,1.66,2.19,
     $  2.21,2.1,2.11,2.35,2.31,1.46,0.,0.,0./
       data klbm(19,5),eldat(19,5),bmdat(1,19,5),bmdat(2,19,5),
     $ bmdat(3,19,5),bmdat(4,19,5),bmdat(5,19,5),bmdat(6,19,5),
     $ bmdat(7,19,5),bmdat(8,19,5),bmdat(9,19,5),bmdat(10,19,5),
     $ bmdat(11,19,5),bmdat(12,19,5)
     $ /'L(total)',0.,133,133,122,120,140,132,131,134,134,0.,0.,0./
       data klbm(20,5),eldat(20,5),bmdat(1,20,5),bmdat(2,20,5),
     $ bmdat(3,20,5),bmdat(4,20,5),bmdat(5,20,5),bmdat(6,20,5),
     $ bmdat(7,20,5),bmdat(8,20,5),bmdat(9,20,5),bmdat(10,20,5),
     $ bmdat(11,20,5),bmdat(12,20,5)
     $ /'T(in)',0.,1.48,1.45,1.48,1.42,1.83,1.4,1.45,1.36,1.44,0.,0.,0./
       data klbm(21,5),eldat(21,5),bmdat(1,21,5),bmdat(2,21,5),
     $ bmdat(3,21,5),bmdat(4,21,5),bmdat(5,21,5),bmdat(6,21,5),
     $ bmdat(7,21,5),bmdat(8,21,5),bmdat(9,21,5),bmdat(10,21,5),
     $ bmdat(11,21,5),bmdat(12,21,5)
     $ /'T(H+)',0.,1.05,1.07,1.01,
     $               1.01,1.03,1.14,1.05,1.06,1.03,0.,0.,0./
       data klbm(22,5),eldat(22,5),bmdat(1,22,5),bmdat(2,22,5),
     $ bmdat(3,22,5),bmdat(4,22,5),bmdat(5,22,5),bmdat(6,22,5),
     $ bmdat(7,22,5),bmdat(8,22,5),bmdat(9,22,5),bmdat(10,22,5),
     $ bmdat(11,22,5),bmdat(12,22,5)
     $ /'<He+>/<H+>',0.,0.92,0.,0.92,0.92,0.,
     $         0.92,0.91,0.92,0.92,0.,0.,0./
c     Table 5 Meudon Planetary Nebula,  
c             Meu,Lex,1,2,3,4,5,6,10,11,0.,0./
       data nbmdat(4)/10/
       data nlnsv2(4)/39/
       data klbm(1,4),eldat(1,4),bmdat(1,1,4),bmdat(2,1,4),
     $ bmdat(3,1,4),bmdat(4,1,4),bmdat(5,1,4),bmdat(6,1,4),
     $ bmdat(7,1,4),bmdat(8,1,4),bmdat(9,1,4),bmdat(10,1,4),
     $ bmdat(11,1,4),bmdat(12,1,4)
     $ /' L(Hp)',0.,2.6,2.53,2.06,2.63,
     $            2.68,2.35,2.73,2.68,2.3,2.8,0.,0./
       data klbm(2,4),eldat(2,4),bmdat(1,2,4),bmdat(2,2,4),
     $ bmdat(3,2,4),bmdat(4,2,4),bmdat(5,2,4),bmdat(6,2,4),
     $ bmdat(7,2,4),bmdat(8,2,4),bmdat(9,2,4),bmdat(10,2,4),
     $ bmdat(11,2,4),bmdat(12,2,4)
     $ /'He I',5876,0.11,0.09,0.09,0.11,
     $  0.1,0.05,0.1,0.11,0.09,0.11,0.,0./
       data klbm(3,4),eldat(3,4),bmdat(1,3,4),bmdat(2,3,4),
     $ bmdat(3,3,4),bmdat(4,3,4),bmdat(5,3,4),bmdat(6,3,4),
     $ bmdat(7,3,4),bmdat(8,3,4),bmdat(9,3,4),bmdat(10,3,4),
     $ bmdat(11,3,4),bmdat(12,3,4)
     $ /'He II',4686,0.33,0.41,0.4,0.32,
     $  0.33,0.81,0.35,0.32,0.43,0.34,0.,0./
       data klbm(4,4),eldat(4,4),bmdat(1,4,4),bmdat(2,4,4),
     $ bmdat(3,4,4),bmdat(4,4,4),bmdat(5,4,4),bmdat(6,4,4),
     $ bmdat(7,4,4),bmdat(8,4,4),bmdat(9,4,4),bmdat(10,4,4),
     $ bmdat(11,4,4),bmdat(12,4,4)
     $ /'C II]',2326,0.38,0.27,0.2,0.33,
     $  0.43,0.12,0.27,0.3,0.22,0.32,0.,0./
       data klbm(5,4),eldat(5,4),bmdat(1,5,4),bmdat(2,5,4),
     $ bmdat(3,5,4),bmdat(4,5,4),bmdat(5,5,4),bmdat(6,5,4),
     $ bmdat(7,5,4),bmdat(8,5,4),bmdat(9,5,4),bmdat(10,5,4),
     $ bmdat(11,5,4),bmdat(12,5,4)
     $ /'C III]',1909,1.7,2.14,2.4,1.82,
     $  1.66,2.92,1.72,1.87,3.14,1.63,0.,0./
       data klbm(6,4),eldat(6,4),bmdat(1,6,4),bmdat(2,6,4),
     $ bmdat(3,6,4),bmdat(4,6,4),bmdat(5,6,4),bmdat(6,6,4),
     $ bmdat(7,6,4),bmdat(8,6,4),bmdat(9,6,4),bmdat(10,6,4),
     $ bmdat(11,6,4),bmdat(12,6,4)
     $ /'C IV',1549,1.64,2.51,2.6,2.44,
     $  2.05,1.44,2.66,2.18,4.74,1.94,0.,0./
       data klbm(7,4),eldat(7,4),bmdat(1,7,4),bmdat(2,7,4),
     $ bmdat(3,7,4),bmdat(4,7,4),bmdat(5,7,4),bmdat(6,7,4),
     $ bmdat(7,7,4),bmdat(8,7,4),bmdat(9,7,4),bmdat(10,7,4),
     $ bmdat(11,7,4),bmdat(12,7,4)
     $ /'[N II]',6584,1.44,1.49,1.43,1.59,
     $  1.45,1.69,1.47,1.44,1.47,1.38,0.,0./
       data klbm(8,4),eldat(8,4),bmdat(1,8,4),bmdat(2,8,4),
     $ bmdat(3,8,4),bmdat(4,8,4),bmdat(5,8,4),bmdat(6,8,4),
     $ bmdat(7,8,4),bmdat(8,8,4),bmdat(9,8,4),bmdat(10,8,4),
     $ bmdat(11,8,4),bmdat(12,8,4)
     $ /'N III]',1749,0.11,0.12,0.16,0.13,
     $  0.13,0.01,0.11,0.13,0.16,0.16,0.,0./
       data klbm(9,4),eldat(9,4),bmdat(1,9,4),bmdat(2,9,4),
     $ bmdat(3,9,4),bmdat(4,9,4),bmdat(5,9,4),bmdat(6,9,4),
     $ bmdat(7,9,4),bmdat(8,9,4),bmdat(9,9,4),bmdat(10,9,4),
     $ bmdat(11,9,4),bmdat(12,9,4)
     $ /'[N III',5.70E+04,0.,0.13,0.11,
     $  0.12,0.13,0.,0.13,0.13,0.13,0.14,0.,0./
       data klbm(10,4),eldat(10,4),bmdat(1,10,4),bmdat(2,10,4),
     $ bmdat(3,10,4),bmdat(4,10,4),bmdat(5,10,4),bmdat(6,10,4),
     $ bmdat(7,10,4),bmdat(8,10,4),bmdat(9,10,4),bmdat(10,10,4),
     $ bmdat(11,10,4),bmdat(12,10,4)
     $ /'N IV]',1487,0.12,0.2,0.22,
     $  0.2,0.15,0.2,0.21,0.19,0.26,0.17,0.,0./
       data klbm(11,4),eldat(11,4),bmdat(1,11,4),bmdat(2,11,4),
     $ bmdat(3,11,4),bmdat(4,11,4),bmdat(5,11,4),bmdat(6,11,4),
     $ bmdat(7,11,4),bmdat(8,11,4),bmdat(9,11,4),bmdat(10,11,4),
     $ bmdat(11,11,4),bmdat(12,11,4)
     $ /'N V',1240,0.09,0.21,0.17,0.18,
     $  0.12,0.34,0.23,0.15,0.38,0.15,0.,0./
       data klbm(12,4),eldat(12,4),bmdat(1,12,4),bmdat(2,12,4),
     $ bmdat(3,12,4),bmdat(4,12,4),bmdat(5,12,4),bmdat(6,12,4),
     $ bmdat(7,12,4),bmdat(8,12,4),bmdat(9,12,4),bmdat(10,12,4),
     $ bmdat(11,12,4),bmdat(12,12,4)
     $ /'[O I]',6300,0.15,0.14,0.17,0.15,
     $  0.12,0.15,0.14,0.14,0.16,0.13,0.,0./
       data klbm(13,4),eldat(13,4),bmdat(1,13,4),bmdat(2,13,4),
     $ bmdat(3,13,4),bmdat(4,13,4),bmdat(5,13,4),bmdat(6,13,4),
     $ bmdat(7,13,4),bmdat(8,13,4),bmdat(9,13,4),bmdat(10,13,4),
     $ bmdat(11,13,4),bmdat(12,13,4)
     $ /'[O II]',3727,2.23,2.28,2.35,
     $  2.23,2.27,0.,2.31,2.18,2.5,2.14,0.,0./
       data klbm(14,4),eldat(14,4),bmdat(1,14,4),bmdat(2,14,4),
     $ bmdat(3,14,4),bmdat(4,14,4),bmdat(5,14,4),bmdat(6,14,4),
     $ bmdat(7,14,4),bmdat(8,14,4),bmdat(9,14,4),bmdat(10,14,4),
     $ bmdat(11,14,4),bmdat(12,14,4)
     $ /'[O lil]',5007,20.9,20.7,21.8,21.1,
     $  21.4,19.8,19.4,21.1,20.2,20.9,0.,0./
       data klbm(15,4),eldat(15,4),bmdat(1,15,4),bmdat(2,15,4),
     $ bmdat(3,15,4),bmdat(4,15,4),bmdat(5,15,4),bmdat(6,15,4),
     $ bmdat(7,15,4),bmdat(8,15,4),bmdat(9,15,4),bmdat(10,15,4),
     $ bmdat(11,15,4),bmdat(12,15,4)
     $ /'[O lil]',4363,0.16,0.16,0.18,0.16,
     $  0.16,0.19,0.14,0.16,0.16,0.15,0.,0./
       data klbm(16,4),eldat(16,4),bmdat(1,16,4),bmdat(2,16,4),
     $ bmdat(3,16,4),bmdat(4,16,4),bmdat(5,16,4),bmdat(6,16,4),
     $ bmdat(7,16,4),bmdat(8,16,4),bmdat(9,16,4),bmdat(10,16,4),
     $ bmdat(11,16,4),bmdat(12,16,4)
     $ /'[Olil]',5.20E+04,1.43,1.34,
     $  1.39,1.42,1.44,0.96,1.4,1.46,1.26,1.41,0.,0./
       data klbm(17,4),eldat(17,4),bmdat(1,17,4),bmdat(2,17,4),
     $ bmdat(3,17,4),bmdat(4,17,4),bmdat(5,17,4),bmdat(6,17,4),
     $ bmdat(7,17,4),bmdat(8,17,4),bmdat(9,17,4),bmdat(10,17,4),
     $ bmdat(11,17,4),bmdat(12,17,4)
     $ /'[OIV]',2.60E+04,3.62,3.92,3.9,
     $  3.52,3.98,4.48,3.32,3.86,5.01,3.33,0.,0./
       data klbm(18,4),eldat(18,4),bmdat(1,18,4),bmdat(2,18,4),
     $ bmdat(3,18,4),bmdat(4,18,4),bmdat(5,18,4),bmdat(6,18,4),
     $ bmdat(7,18,4),bmdat(8,18,4),bmdat(9,18,4),bmdat(10,18,4),
     $ bmdat(11,18,4),bmdat(12,18,4)
     $ /'O IV]',1403,0.13,0.27,0.36,0.2,
     $  0.23,0.18,0.26,0.33,0.41,0.15,0.,0./
       data klbm(19,4),eldat(19,4),bmdat(1,19,4),bmdat(2,19,4),
     $ bmdat(3,19,4),bmdat(4,19,4),bmdat(5,19,4),bmdat(6,19,4),
     $ bmdat(7,19,4),bmdat(8,19,4),bmdat(9,19,4),bmdat(10,19,4),
     $ bmdat(11,19,4),bmdat(12,19,4)
     $ /'O V]',1218,0.09,0.24,0.,0.2,0.11,0.35,0.29,0.19,0.33,0.,0.,0./
       data klbm(20,4),eldat(20,4),bmdat(1,20,4),bmdat(2,20,4),
     $ bmdat(3,20,4),bmdat(4,20,4),bmdat(5,20,4),bmdat(6,20,4),
     $ bmdat(7,20,4),bmdat(8,20,4),bmdat(9,20,4),bmdat(10,20,4),
     $ bmdat(11,20,4),bmdat(12,20,4)
     $ /'[Ne lil]',1.55E+04,2.51,2.49,2.67,
     $  2.75,2.76,0.72,2.8,2.81,2.71,2.74,0.,0./
       data klbm(21,4),eldat(21,4),bmdat(1,21,4),bmdat(2,21,4),
     $ bmdat(3,21,4),bmdat(4,21,4),bmdat(5,21,4),bmdat(6,21,4),
     $ bmdat(7,21,4),bmdat(8,21,4),bmdat(9,21,4),bmdat(10,21,4),
     $ bmdat(11,21,4),bmdat(12,21,4)
     $ /'[Ne lil]',3869,2.59,2.63,3.2,3.33,
     $  2.27,0.88,2.74,2.44,3.35,2.86,0.,0./
       data klbm(22,4),eldat(22,4),bmdat(1,22,4),bmdat(2,22,4),
     $ bmdat(3,22,4),bmdat(4,22,4),bmdat(5,22,4),bmdat(6,22,4),
     $ bmdat(7,22,4),bmdat(8,22,4),bmdat(9,22,4),bmdat(10,22,4),
     $ bmdat(11,22,4),bmdat(12,22,4)
     $ /'Ne IV]',2423,0.56,0.95,1.05,0.72,
     $  0.74,1.64,0.91,0.74,1.19,0.63,0.,0./
       data klbm(23,4),eldat(23,4),bmdat(1,23,4),bmdat(2,23,4),
     $ bmdat(3,23,4),bmdat(4,23,4),bmdat(5,23,4),bmdat(6,23,4),
     $ bmdat(7,23,4),bmdat(8,23,4),bmdat(9,23,4),bmdat(10,23,4),
     $ bmdat(11,23,4),bmdat(12,23,4)
     $ /'[Ne V]',3426,0.73,0.9,0.79,0.74,0.6,
     $  2.29,0.73,0.61,0.81,0.63,0.,0./
       data klbm(24,4),eldat(24,4),bmdat(1,24,4),bmdat(2,24,4),
     $ bmdat(3,24,4),bmdat(4,24,4),bmdat(5,24,4),bmdat(6,24,4),
     $ bmdat(7,24,4),bmdat(8,24,4),bmdat(9,24,4),bmdat(10,24,4),
     $ bmdat(11,24,4),bmdat(12,24,4)
     $ /'[Ne V]',2.42E+04,1.67,0.88,1.2,0.94,
     $  0.76,1.16,0.81,0.99,0.25,0.95,0.,0./
       data klbm(25,4),eldat(25,4),bmdat(1,25,4),bmdat(2,25,4),
     $ bmdat(3,25,4),bmdat(4,25,4),bmdat(5,25,4),bmdat(6,25,4),
     $ bmdat(7,25,4),bmdat(8,25,4),bmdat(9,25,4),bmdat(10,25,4),
     $ bmdat(11,25,4),bmdat(12,25,4)
     $ /'Mg II',2798,1.48,1.56,2.5,2.33,
     $  1.6,0.63,1.22,1.17,1.15,1.92,0.,0./
       data klbm(26,4),eldat(26,4),bmdat(1,26,4),bmdat(2,26,4),
     $ bmdat(3,26,4),bmdat(4,26,4),bmdat(5,26,4),bmdat(6,26,4),
     $ bmdat(7,26,4),bmdat(8,26,4),bmdat(9,26,4),bmdat(10,26,4),
     $ bmdat(11,26,4),bmdat(12,26,4)
     $ /'[Mg IV]',4.50E+03,0.09,0.12,0.11,
     $  0.12,0.13,0.,0.,0.12,0.14,0.,0.,0./
       data klbm(27,4),eldat(27,4),bmdat(1,27,4),bmdat(2,27,4),
     $ bmdat(3,27,4),bmdat(4,27,4),bmdat(5,27,4),bmdat(6,27,4),
     $ bmdat(7,27,4),bmdat(8,27,4),bmdat(9,27,4),bmdat(10,27,4),
     $ bmdat(11,27,4),bmdat(12,27,4)
     $ /'[Si II]',3.48E+04,0.13,0.18,0.14,0.16,
     $  0.26,0.,0.19,0.17,0.15,0.16,0.,0./
       data klbm(28,4),eldat(28,4),bmdat(1,28,4),bmdat(2,28,4),
     $ bmdat(3,28,4),bmdat(4,28,4),bmdat(5,28,4),bmdat(6,28,4),
     $ bmdat(7,28,4),bmdat(8,28,4),bmdat(9,28,4),bmdat(10,28,4),
     $ bmdat(11,28,4),bmdat(12,28,4)
     $ /'Si II]',2335,0.11,0.23,0.23,0.15,0.,
     $           0.53,0.16,0.16,0.,0.15,0.,0./
       data klbm(29,4),eldat(29,4),bmdat(1,29,4),bmdat(2,29,4),
     $ bmdat(3,29,4),bmdat(4,29,4),bmdat(5,29,4),bmdat(6,29,4),
     $ bmdat(7,29,4),bmdat(8,29,4),bmdat(9,29,4),bmdat(10,29,4),
     $ bmdat(11,29,4),bmdat(12,29,4)
     $ /'Si III]',1892,0.2,0.68,0.79,0.39,
     $               0.32,1.95,0.46,0.45,0.,0.4,0.,0./
       data klbm(30,4),eldat(30,4),bmdat(1,30,4),bmdat(2,30,4),
     $ bmdat(3,30,4),bmdat(4,30,4),bmdat(5,30,4),bmdat(6,30,4),
     $ bmdat(7,30,4),bmdat(8,30,4),bmdat(9,30,4),bmdat(10,30,4),
     $ bmdat(11,30,4),bmdat(12,30,4)
     $ /'Si IV',1397,0.15,0.15,0.1,0.2,
     $  0.15,0.03,0.21,0.17,0.,0.16,0.,0./
       data klbm(31,4),eldat(31,4),bmdat(1,31,4),bmdat(2,31,4),
     $ bmdat(3,31,4),bmdat(4,31,4),bmdat(5,31,4),bmdat(6,31,4),
     $ bmdat(7,31,4),bmdat(8,31,4),bmdat(9,31,4),bmdat(10,31,4),
     $ bmdat(11,31,4),bmdat(12,31,4)
     $ /'[S II]',6720,0.39,0.33,0.24,0.21,
     $  0.45,0.08,0.33,0.43,0.41,0.51,0.,0./
       data klbm(32,4),eldat(32,4),bmdat(1,32,4),bmdat(2,32,4),
     $ bmdat(3,32,4),bmdat(4,32,4),bmdat(5,32,4),bmdat(6,32,4),
     $ bmdat(7,32,4),bmdat(8,32,4),bmdat(9,32,4),bmdat(10,32,4),
     $ bmdat(11,32,4),bmdat(12,32,4)
     $ /'[S 111]',1.87E+04,0.49,0.53,0.6,
     $  0.48,0.49,0.,0.46,0.49,0.6,0.55,0.,0./
       data klbm(33,4),eldat(33,4),bmdat(1,33,4),bmdat(2,33,4),
     $ bmdat(3,33,4),bmdat(4,33,4),bmdat(5,33,4),bmdat(6,33,4),
     $ bmdat(7,33,4),bmdat(8,33,4),bmdat(9,33,4),bmdat(10,33,4),
     $ bmdat(11,33,4),bmdat(12,33,4)
     $ /'[S lil]',9532,2.09,1.91,2.31,2.04,1.89,
     $  0.36,2.05,1.87,2.34,2.42,0.,0./
       data klbm(34,4),eldat(34,4),bmdat(1,34,4),bmdat(2,34,4),
     $ bmdat(3,34,4),bmdat(4,34,4),bmdat(5,34,4),bmdat(6,34,4),
     $ bmdat(7,34,4),bmdat(8,34,4),bmdat(9,34,4),bmdat(10,34,4),
     $ bmdat(11,34,4),bmdat(12,34,4)
     $ /'[S IV]',1.05E+04,1.92,1.84,1.58,1.92,
     $  2.21,0.93,1.81,1.98,2.36,1.94,0.,0./
       data klbm(35,4),eldat(35,4),bmdat(1,35,4),bmdat(2,35,4),
     $ bmdat(3,35,4),bmdat(4,35,4),bmdat(5,35,4),bmdat(6,35,4),
     $ bmdat(7,35,4),bmdat(8,35,4),bmdat(9,35,4),bmdat(10,35,4),
     $ bmdat(11,35,4),bmdat(12,35,4)
     $ /'L(total)',0.,129,132,114,139,136,105,135,136,130,142,0.,0./
       data klbm(36,4),eldat(36,4),bmdat(1,36,4),bmdat(2,36,4),
     $ bmdat(3,36,4),bmdat(4,36,4),bmdat(5,36,4),bmdat(6,36,4),
     $ bmdat(7,36,4),bmdat(8,36,4),bmdat(9,36,4),bmdat(10,36,4),
     $ bmdat(11,36,4),bmdat(12,36,4)
     $ /'T(in)',0.,0.,1.8,0.,1.83,1.78,1.63,1.84,1.78,1.95,1.81,0.,0./
       data klbm(37,4),eldat(37,4),bmdat(1,37,4),bmdat(2,37,4),
     $ bmdat(3,37,4),bmdat(4,37,4),bmdat(5,37,4),bmdat(6,37,4),
     $ bmdat(7,37,4),bmdat(8,37,4),bmdat(9,37,4),bmdat(10,37,4),
     $ bmdat(11,37,4),bmdat(12,37,4)
     $ /'T(H+)',0.,0.,1.26,0.,1.22,1.21,1.32,1.35,1.21,1.29,1.2,0.,0./
       data klbm(38,4),eldat(38,4),bmdat(1,38,4),bmdat(2,38,4),
     $ bmdat(3,38,4),bmdat(4,38,4),bmdat(5,38,4),bmdat(6,38,4),
     $ bmdat(7,38,4),bmdat(8,38,4),bmdat(9,38,4),bmdat(10,38,4),
     $ bmdat(11,38,4),bmdat(12,38,4)
     $ /'<He+>/<H+>',0.,0.,0.69,0.,0.74,
     $       0.74,0.,0.71,0.71,0.6,0.68,0.,0./
       data klbm(39,4),eldat(39,4),bmdat(1,39,4),bmdat(2,39,4),
     $ bmdat(3,39,4),bmdat(4,39,4),bmdat(5,39,4),bmdat(6,39,4),
     $ bmdat(7,39,4),bmdat(8,39,4),bmdat(9,39,4),bmdat(10,39,4),
     $ bmdat(11,39,4),bmdat(12,39,4)
     $ /'R(out) E17',0.,0.,4.02,0.,4.04,4.04,0.,
     $         4.07,4.07,3.83,4.08,0.,0./
c     Table 7 Low Ionization PN, Mean,1,2,3,5,6,10,11,0.,0.,0.,0./
       data nbmdat(6)/8/ 
       data nlnsv2(6)/17/
       data klbm(1,6),eldat(1,6),bmdat(1,1,6),bmdat(2,1,6),
     $ bmdat(3,1,6),bmdat(4,1,6),bmdat(5,1,6),bmdat(6,1,6),
     $ bmdat(7,1,6),bmdat(8,1,6),bmdat(9,1,6),bmdat(10,1,6),
     $ bmdat(11,1,6),bmdat(12,1,6)
     $ /'L(H A)',0.,5.41,5.2,5.56,5.52,
     $                 5.35,5.41,5.38,5.44,0.,0.,0.,0./
       data klbm(2,6),eldat(2,6),bmdat(1,2,6),bmdat(2,2,6),
     $ bmdat(3,2,6),bmdat(4,2,6),bmdat(5,2,6),bmdat(6,2,6),
     $ bmdat(7,2,6),bmdat(8,2,6),bmdat(9,2,6),bmdat(10,2,6),
     $ bmdat(11,2,6),bmdat(12,2,6)
     $ /'He I',5876,0.13,0.12,0.14,0.12,0.14,0.12,0.12,0.15,0.,0.,0.,0./
       data klbm(3,6),eldat(3,6),bmdat(1,3,6),bmdat(2,3,6),
     $ bmdat(3,3,6),bmdat(4,3,6),bmdat(5,3,6),bmdat(6,3,6),
     $ bmdat(7,3,6),bmdat(8,3,6),bmdat(9,3,6),bmdat(10,3,6),
     $ bmdat(11,3,6),bmdat(12,3,6)
     $ /'He II',4686,0.088,0.095,0.082,
     $  0.088,0.088,0.088,0.091,0.086,0.,0.,0.,0./
       data klbm(4,6),eldat(4,6),bmdat(1,4,6),bmdat(2,4,6),
     $ bmdat(3,4,6),bmdat(4,4,6),bmdat(5,4,6),bmdat(6,4,6),
     $ bmdat(7,4,6),bmdat(8,4,6),bmdat(9,4,6),bmdat(10,4,6),
     $ bmdat(11,4,6),bmdat(12,4,6)
     $ /'C III]',1909,1.17,1.53,0.78,0.81,1.25,1.41,1.43,1,0.,0.,0.,0./
       data klbm(5,6),eldat(5,6),bmdat(1,5,6),bmdat(2,5,6),
     $ bmdat(3,5,6),bmdat(4,5,6),bmdat(5,5,6),bmdat(6,5,6),
     $ bmdat(7,5,6),bmdat(8,5,6),bmdat(9,5,6),bmdat(10,5,6),
     $ bmdat(11,5,6),bmdat(12,5,6)
     $ /'C IV',1549,1.38,1.2,1.34,1.31,1.54,1.43,1.51,1.32,0.,0.,0.,0./
       data klbm(6,6),eldat(6,6),bmdat(1,6,6),bmdat(2,6,6),
     $ bmdat(3,6,6),bmdat(4,6,6),bmdat(5,6,6),bmdat(6,6,6),
     $ bmdat(7,6,6),bmdat(8,6,6),bmdat(9,6,6),bmdat(10,6,6),
     $ bmdat(11,6,6),bmdat(12,6,6)
     $ /'[O III]',5007,14.3,16,12.7,13.1,
     $             15.,14.5,14.6,14.3,0.,0.,0.,0./
       data klbm(7,6),eldat(7,6),bmdat(1,7,6),bmdat(2,7,6),
     $ bmdat(3,7,6),bmdat(4,7,6),bmdat(5,7,6),bmdat(6,7,6),
     $ bmdat(7,7,6),bmdat(8,7,6),bmdat(9,7,6),bmdat(10,7,6),
     $ bmdat(11,7,6),bmdat(12,7,6)
     $ /'[OIII]',5.20E+04,0.26,0.28,0.26,
     $  0.26,0.26,0.27,0.,0.26,0.,0.,0.,0./
       data klbm(8,6),eldat(8,6),bmdat(1,8,6),bmdat(2,8,6),
     $ bmdat(3,8,6),bmdat(4,8,6),bmdat(5,8,6),bmdat(6,8,6),
     $ bmdat(7,8,6),bmdat(8,8,6),bmdat(9,8,6),bmdat(10,8,6),
     $ bmdat(11,8,6),bmdat(12,8,6)
     $ /'[O IV]',2.60E+04,0.22,0.24,0.21,
     $  0.23,0.21,0.21,0.22,0.22,0.,0.,0.,0./
       data klbm(9,6),eldat(9,6),bmdat(1,9,6),bmdat(2,9,6),
     $ bmdat(3,9,6),bmdat(4,9,6),bmdat(5,9,6),bmdat(6,9,6),
     $ bmdat(7,9,6),bmdat(8,9,6),bmdat(9,9,6),bmdat(10,9,6),
     $ bmdat(11,9,6),bmdat(12,9,6)
     $ /'[Ne III]',1.55E+04,1.11,1.14,1.09,
     $  1.1,1.1,1.12,1.12,1.1,0.,0.,0.,0./
       data klbm(10,6),eldat(10,6),bmdat(1,10,6),bmdat(2,10,6),
     $ bmdat(3,10,6),bmdat(4,10,6),bmdat(5,10,6),bmdat(6,10,6),
     $ bmdat(7,10,6),bmdat(8,10,6),bmdat(9,10,6),bmdat(10,10,6),
     $ bmdat(11,10,6),bmdat(12,10,6)
     $ /'[Ne III]',3869,1.44,1.66,1.39,1.27,
     $  1.48,1.44,1.45,1.38,0.,0.,0.,0./
       data klbm(11,6),eldat(11,6),bmdat(1,11,6),bmdat(2,11,6),
     $ bmdat(3,11,6),bmdat(4,11,6),bmdat(5,11,6),bmdat(6,11,6),
     $ bmdat(7,11,6),bmdat(8,11,6),bmdat(9,11,6),bmdat(10,11,6),
     $ bmdat(11,11,6),bmdat(12,11,6)
     $ /'Ne IV]',2423,0.1,0.1,0.11,0.11,0.1,0.08,0.08,0.1,0.,0.,0.,0./
       data klbm(12,6),eldat(12,6),bmdat(1,12,6),bmdat(2,12,6),
     $ bmdat(3,12,6),bmdat(4,12,6),bmdat(5,12,6),bmdat(6,12,6),
     $ bmdat(7,12,6),bmdat(8,12,6),bmdat(9,12,6),bmdat(10,12,6),
     $ bmdat(11,12,6),bmdat(12,12,6)
     $ /'[S III]',9532,0.52,0.65,0.23,0.37,
     $      0.26,0.42,0.76,0.99,0.,0.,0.,0./
       data klbm(13,6),eldat(13,6),bmdat(1,13,6),bmdat(2,13,6),
     $ bmdat(3,13,6),bmdat(4,13,6),bmdat(5,13,6),bmdat(6,13,6),
     $ bmdat(7,13,6),bmdat(8,13,6),bmdat(9,13,6),bmdat(10,13,6),
     $ bmdat(11,13,6),bmdat(12,13,6)
     $ /'[SII]',1.05E+04,1.43,1.2,
     $  1.32,1.57,1.32,1.84,1.38,1.37,0.,0.,0.,0./
       data klbm(14,6),eldat(14,6),bmdat(1,14,6),bmdat(2,14,6),
     $ bmdat(3,14,6),bmdat(4,14,6),bmdat(5,14,6),bmdat(6,14,6),
     $ bmdat(7,14,6),bmdat(8,14,6),bmdat(9,14,6),bmdat(10,14,6),
     $ bmdat(11,14,6),bmdat(12,14,6)
     $ /'L(total)',0.,120,126,109,112,122,124,122,121,0.,0.,0.,0./
       data klbm(15,6),eldat(15,6),bmdat(1,15,6),bmdat(2,15,6),
     $ bmdat(3,15,6),bmdat(4,15,6),bmdat(5,15,6),bmdat(6,15,6),
     $ bmdat(7,15,6),bmdat(8,15,6),bmdat(9,15,6),bmdat(10,15,6),
     $ bmdat(11,15,6),bmdat(12,15,6)
     $ /'T(in)',0.,1.78,0.,1.83,1.79,1.76,1.73,1.76,1.81,0.,0.,0.,0./
       data klbm(16,6),eldat(16,6),bmdat(1,16,6),bmdat(2,16,6),
     $ bmdat(3,16,6),bmdat(4,16,6),bmdat(5,16,6),bmdat(6,16,6),
     $ bmdat(7,16,6),bmdat(8,16,6),bmdat(9,16,6),bmdat(10,16,6),
     $ bmdat(11,16,6),bmdat(12,16,6)
     $ /'T(H+)',0.,1.16,0.,1.11,1.11,1.28,1.14,1.2,1.14,0.,0.,0.,0./
       data klbm(17,6),eldat(17,6),bmdat(1,17,6),bmdat(2,17,6),
     $ bmdat(3,17,6),bmdat(4,17,6),bmdat(5,17,6),bmdat(6,17,6),
     $ bmdat(7,17,6),bmdat(8,17,6),bmdat(9,17,6),bmdat(10,17,6),
     $ bmdat(11,17,6),bmdat(12,17,6)
     $ /'<He.>/cH.>',0.,0.91,0.,0.91,0.91,
     $      0.91,0.91,0.9,0.9,0.,0.,0.,0./
c     Table 8 NLR Cloud,Mean,1,2,4,5,6,11,
       data nbmdat(7)/7/
       data nlnsv2(7)/29/
      data klbm(30,7),eldat(30,7),bmdat(1,30,7),bmdat(2,30,7),
     $ bmdat(3,30,7),bmdat(4,30,7),bmdat(5,30,7),bmdat(6,30,7),
     $ bmdat(7,30,7),bmdat(8,30,7),bmdat(9,30,7),bmdat(10,30,7),
     $ bmdat(11,30,7),bmdat(12,30,7)
     $ /'6',0.,0.,0.,0.,1.06,1.37,1.43,1.34,0.,0.,0.,0.,0./
       data klbm(1,7),eldat(1,7),bmdat(1,1,7),bmdat(2,1,7),
     $ bmdat(3,1,7),bmdat(4,1,7),bmdat(5,1,7),bmdat(6,1,7),
     $ bmdat(7,1,7),bmdat(8,1,7),bmdat(9,1,7),bmdat(10,1,7),
     $ bmdat(11,1,7),bmdat(12,1,7)
     $ /'Lyoc',1216,34.2,38.3,32.1,37,32.4,31.5,34.2,0.,0.,0.,0.,0./
       data klbm(2,7),eldat(2,7),bmdat(1,2,7),bmdat(2,2,7),
     $ bmdat(3,2,7),bmdat(4,2,7),bmdat(5,2,7),bmdat(6,2,7),
     $ bmdat(7,2,7),bmdat(8,2,7),bmdat(9,2,7),bmdat(10,2,7),
     $ bmdat(11,2,7),bmdat(12,2,7)
     $ /'He I',5876,0.12,0.11,0.13,0.14,0.12,0.13,0.13,0.,0.,0.,0.,0./
       data klbm(3,7),eldat(3,7),bmdat(1,3,7),bmdat(2,3,7),
     $ bmdat(3,3,7),bmdat(4,3,7),bmdat(5,3,7),bmdat(6,3,7),
     $ bmdat(7,3,7),bmdat(8,3,7),bmdat(9,3,7),bmdat(10,3,7),
     $ bmdat(11,3,7),bmdat(12,3,7)
     $ /'He II',4686,0.24,0.25,0.25,0.,
     $               0.25,0.23,0.24,0.,0.,0.,0.,0./
       data klbm(4,7),eldat(4,7),bmdat(1,4,7),bmdat(2,4,7),
     $ bmdat(3,4,7),bmdat(4,4,7),bmdat(5,4,7),bmdat(6,4,7),
     $ bmdat(7,4,7),bmdat(8,4,7),bmdat(9,4,7),bmdat(10,4,7),
     $ bmdat(11,4,7),bmdat(12,4,7)
     $ /'He II',1640,1.6,1.6,1.74,1.49,1.53,1.56,1.67,0.,0.,0.,0.,0./
       data klbm(5,7),eldat(5,7),bmdat(1,5,7),bmdat(2,5,7),
     $ bmdat(3,5,7),bmdat(4,5,7),bmdat(5,5,7),bmdat(6,5,7),
     $ bmdat(7,5,7),bmdat(8,5,7),bmdat(9,5,7),bmdat(10,5,7),
     $ bmdat(11,5,7),bmdat(12,5,7)
     $ /'CIII]',1909,2.82,2.9,2.99,2.45,2.87,2.83,2.9,0.,0.,0.,0.,0./
       data klbm(6,7),eldat(6,7),bmdat(1,6,7),bmdat(2,6,7),
     $ bmdat(3,6,7),bmdat(4,6,7),bmdat(5,6,7),bmdat(6,6,7),
     $ bmdat(7,6,7),bmdat(8,6,7),bmdat(9,6,7),bmdat(10,6,7),
     $ bmdat(11,6,7),bmdat(12,6,7)
     $ /'CIV',1549,3.18,2.7,3.85,2.28,3.69,3.17,3.36,0.,0.,0.,0.,0./
       data klbm(7,7),eldat(7,7),bmdat(1,7,7),bmdat(2,7,7),
     $ bmdat(3,7,7),bmdat(4,7,7),bmdat(5,7,7),bmdat(6,7,7),
     $ bmdat(7,7,7),bmdat(8,7,7),bmdat(9,7,7),bmdat(10,7,7),
     $ bmdat(11,7,7),bmdat(12,7,7)
     $ /'[NII]',6584,2.33,1.4,3.2,1.21,3.1,2.67,2.4,0.,0.,0.,0.,0./
       data klbm(8,7),eldat(8,7),bmdat(1,8,7),bmdat(2,8,7),
     $ bmdat(3,8,7),bmdat(4,8,7),bmdat(5,8,7),bmdat(6,8,7),
     $ bmdat(7,8,7),bmdat(8,8,7),bmdat(9,8,7),bmdat(10,8,7),
     $ bmdat(11,8,7),bmdat(12,8,7)
     $ /'NIII]',1749,0.19,0.24,0.24,0.01,0.22,0.22,0.22,0.,0.,0.,0.,0./
       data klbm(9,7),eldat(9,7),bmdat(1,9,7),bmdat(2,9,7),
     $ bmdat(3,9,7),bmdat(4,9,7),bmdat(5,9,7),bmdat(6,9,7),
     $ bmdat(7,9,7),bmdat(8,9,7),bmdat(9,9,7),bmdat(10,9,7),
     $ bmdat(11,9,7),bmdat(12,9,7)
     $ /'NIV]',1487,0.2,0.2,0.23,0.12,0.22,0.21,0.21,0.,0.,0.,0.,0./
       data klbm(10,7),eldat(10,7),bmdat(1,10,7),bmdat(2,10,7),
     $ bmdat(3,10,7),bmdat(4,10,7),bmdat(5,10,7),bmdat(6,10,7),
     $ bmdat(7,10,7),bmdat(8,10,7),bmdat(9,10,7),bmdat(10,10,7),
     $ bmdat(11,10,7),bmdat(12,10,7)
     $ /'[OI]',6300,1.61,2.2,1.61,1.41,1.67,1.31,1.46,0.,0.,0.,0.,0./
       data klbm(11,7),eldat(11,7),bmdat(1,11,7),bmdat(2,11,7),
     $ bmdat(3,11,7),bmdat(4,11,7),bmdat(5,11,7),bmdat(6,11,7),
     $ bmdat(7,11,7),bmdat(8,11,7),bmdat(9,11,7),bmdat(10,11,7),
     $ bmdat(11,11,7),bmdat(12,11,7)
     $ /'[O II]',6.30E+04,1.12,0.25,1.13,0.,0.,1.44,1.64,0.,0.,0.,0.,0./
       data klbm(12,7),eldat(12,7),bmdat(1,12,7),bmdat(2,12,7),
     $ bmdat(3,12,7),bmdat(4,12,7),bmdat(5,12,7),bmdat(6,12,7),
     $ bmdat(7,12,7),bmdat(8,12,7),bmdat(9,12,7),bmdat(10,12,7),
     $ bmdat(11,12,7),bmdat(12,12,7)
     $ /'[OIII]',3727,1.72,1.6,1.44,3.18,1.58,1.3,1.2,0.,0.,0.,0.,0./
       data klbm(13,7),eldat(13,7),bmdat(1,13,7),bmdat(2,13,7),
     $ bmdat(3,13,7),bmdat(4,13,7),bmdat(5,13,7),bmdat(6,13,7),
     $ bmdat(7,13,7),bmdat(8,13,7),bmdat(9,13,7),bmdat(10,13,7),
     $ bmdat(11,13,7),bmdat(12,13,7)
     $ /'OIII]',1663,0.56,0.35,0.63,0.,0.61,0.57,0.63,0.,0.,0.,0.,0./
       data klbm(14,7),eldat(14,7),bmdat(1,14,7),bmdat(2,14,7),
     $ bmdat(3,14,7),bmdat(4,14,7),bmdat(5,14,7),bmdat(6,14,7),
     $ bmdat(7,14,7),bmdat(8,14,7),bmdat(9,14,7),bmdat(10,14,7),
     $ bmdat(11,14,7),bmdat(12,14,7)
     $ /'[OIII]',5007,33.1,31.4,34.5,31.1,33,32.8,36,0.,0.,0.,0.,0./
       data klbm(15,7),eldat(15,7),bmdat(1,15,7),bmdat(2,15,7),
     $ bmdat(3,15,7),bmdat(4,15,7),bmdat(5,15,7),bmdat(6,15,7),
     $ bmdat(7,15,7),bmdat(8,15,7),bmdat(9,15,7),bmdat(10,15,7),
     $ bmdat(11,15,7),bmdat(12,15,7)
     $ /'[OIII]',4363,0.32,0.3,0.34,0.,
     $             0.31,0.3,0.33,0.,0.,0.,0.,0./
       data klbm(16,7),eldat(16,7),bmdat(1,16,7),bmdat(2,16,7),
     $ bmdat(3,16,7),bmdat(4,16,7),bmdat(5,16,7),bmdat(6,16,7),
     $ bmdat(7,16,7),bmdat(8,16,7),bmdat(9,16,7),bmdat(10,16,7),
     $ bmdat(11,16,7),bmdat(12,16,7)
     $ /'OIV',1403,0.36,0.49,0.3,0.,0.36,0.42,0.25,0.,0.,0.,0.,0./
       data klbm(17,7),eldat(17,7),bmdat(1,17,7),bmdat(2,17,7),
     $ bmdat(3,17,7),bmdat(4,17,7),bmdat(5,17,7),bmdat(6,17,7),
     $ bmdat(7,17,7),bmdat(8,17,7),bmdat(9,17,7),bmdat(10,17,7),
     $ bmdat(11,17,7),bmdat(12,17,7)
     $ /'[NeIII]',1.55E+04,1.89,1.5,2.01,0.,
     $        1.94,2.05,1.95,0.,0.,0.,0.,0./
       data klbm(18,7),eldat(18,7),bmdat(1,18,7),bmdat(2,18,7),
     $ bmdat(3,18,7),bmdat(4,18,7),bmdat(5,18,7),bmdat(6,18,7),
     $ bmdat(7,18,7),bmdat(8,18,7),bmdat(9,18,7),bmdat(10,18,7),
     $ bmdat(11,18,7),bmdat(12,18,7)
     $ /'[Ne III]',3869,1.91,1.9,2.51,0.84,
     $              2.16,1.72,2.34,0.,0.,0.,0.,0./
       data klbm(19,7),eldat(19,7),bmdat(1,19,7),bmdat(2,19,7),
     $ bmdat(3,19,7),bmdat(4,19,7),bmdat(5,19,7),bmdat(6,19,7),
     $ bmdat(7,19,7),bmdat(8,19,7),bmdat(9,19,7),bmdat(10,19,7),
     $ bmdat(11,19,7),bmdat(12,19,7)
     $ /'[Ne IV]',2423,0.44,0.52,0.42,
     $            0.,0.47,0.41,0.38,0.,0.,0.,0.,0./
       data klbm(20,7),eldat(20,7),bmdat(1,20,7),bmdat(2,20,7),
     $ bmdat(3,20,7),bmdat(4,20,7),bmdat(5,20,7),bmdat(6,20,7),
     $ bmdat(7,20,7),bmdat(8,20,7),bmdat(9,20,7),bmdat(10,20,7),
     $ bmdat(11,20,7),bmdat(12,20,7)
     $ /'[NeV]',3426,0.52,0.59,0.55,0.,0 53,0.44,0.5,0.,0.,0.,0.,0./
       data klbm(21,7),eldat(21,7),bmdat(1,21,7),bmdat(2,21,7),
     $ bmdat(3,21,7),bmdat(4,21,7),bmdat(5,21,7),bmdat(6,21,7),
     $ bmdat(7,21,7),bmdat(8,21,7),bmdat(9,21,7),bmdat(10,21,7),
     $ bmdat(11,21,7),bmdat(12,21,7)
     $ /'MgII',2798,1.78,3.5,1.72,1.48,1.23,1.12,1.61,0.,0.,0.,0.,0./
       data klbm(22,7),eldat(22,7),bmdat(1,22,7),bmdat(2,22,7),
     $ bmdat(3,22,7),bmdat(4,22,7),bmdat(5,22,7),bmdat(6,22,7),
     $ bmdat(7,22,7),bmdat(8,22,7),bmdat(9,22,7),bmdat(10,22,7),
     $ bmdat(11,22,7),bmdat(12,22,7)
     $ /'[SiII]',3.48E+04,0.9,1,0.96,0.,
     $           1.07,0.96,0.52,0.,0.,0.,0.,0./
       data klbm(23,7),eldat(23,7),bmdat(1,23,7),bmdat(2,23,7),
     $ bmdat(3,23,7),bmdat(4,23,7),bmdat(5,23,7),bmdat(6,23,7),
     $ bmdat(7,23,7),bmdat(8,23,7),bmdat(9,23,7),bmdat(10,23,7),
     $ bmdat(11,23,7),bmdat(12,23,7)
     $ /'[SII]',6720,1.33,2.4,1.01,1.58,0.93,0.99,1.1,0.,0.,0.,0.,0./
       data klbm(24,7),eldat(24,7),bmdat(1,24,7),bmdat(2,24,7),
     $ bmdat(3,24,7),bmdat(4,24,7),bmdat(5,24,7),bmdat(6,24,7),
     $ bmdat(7,24,7),bmdat(8,24,7),bmdat(9,24,7),bmdat(10,24,7),
     $ bmdat(11,24,7),bmdat(12,24,7)
     $ /'[S III]',9532,1.88,1.6,2.15,1.73,
     $               2.06,1.67,2.08,0.,0.,0.,0.,0./
       data klbm(25,7),eldat(25,7),bmdat(1,25,7),bmdat(2,25,7),
     $ bmdat(3,25,7),bmdat(4,25,7),bmdat(5,25,7),bmdat(6,25,7),
     $ bmdat(7,25,7),bmdat(8,25,7),bmdat(9,25,7),bmdat(10,25,7),
     $ bmdat(11,25,7),bmdat(12,25,7)
     $ /'[SIII',1.87E+04,0.49,0.36,0.61,0.,
     $             0.57,0.52,0.37,0.,0.,0.,0.,0./
       data klbm(26,7),eldat(26,7),bmdat(1,26,7),bmdat(2,26,7),
     $ bmdat(3,26,7),bmdat(4,26,7),bmdat(5,26,7),bmdat(6,26,7),
     $ bmdat(7,26,7),bmdat(8,26,7),bmdat(9,26,7),bmdat(10,26,7),
     $ bmdat(11,26,7),bmdat(12,26,7)
     $ /'[SIII]',1.05E+04,1.05,0.86,
     $  1.24,1.23,0.82,0.94,1.22,0.,0.,0.,0.,0./
       data klbm(27,7),eldat(27,7),bmdat(1,27,7),bmdat(2,27,7),
     $ bmdat(3,27,7),bmdat(4,27,7),bmdat(5,27,7),bmdat(6,27,7),
     $ bmdat(7,27,7),bmdat(8,27,7),bmdat(9,27,7),bmdat(10,27,7),
     $ bmdat(11,27,7),bmdat(12,27,7)
     $ /'I(total) E0',0.,125,131,128,92,128,131,133,0.,0.,0.,0.,0./
       data klbm(28,7),eldat(28,7),bmdat(1,28,7),bmdat(2,28,7),
     $ bmdat(3,28,7),bmdat(4,28,7),bmdat(5,28,7),bmdat(6,28,7),
     $ bmdat(7,28,7),bmdat(8,28,7),bmdat(9,28,7),bmdat(10,28,7),
     $ bmdat(11,28,7),bmdat(12,28,7)
     $ /'T(in) E4',0.,1.7,1.71,1.7,0.,1.72,1.68,1.68,0.,0.,0.,0.,0./
       data klbm(29,7),eldat(29,7),bmdat(1,29,7),bmdat(2,29,7),
     $ bmdat(3,29,7),bmdat(4,29,7),bmdat(5,29,7),bmdat(6,29,7),
     $ bmdat(7,29,7),bmdat(8,29,7),bmdat(9,29,7),bmdat(10,29,7),
     $ bmdat(11,29,7),bmdat(12,29,7)
     $ /'T(H+) E4',0.,1.17,0.,1.24,1.12,
     $         1.06,1.2,1.23,0.,0.,0.,0.,0./
      data nlnsv(2)/34/
      data (elnsv(1,ml,2),elnsv(2,ml,2),nbmpt(ml,2),ml=1,17)/
     $    4858., 4863.,24,
     $    1215., 1217.,0,
     $    5875., 5877.,1,
     $    2325., 2327.,2,
     $    1334., 1336.,0,
     $    1908., 1910.,3,
     $    6545., 6550.,5,
     $    6580., 6590.,5,
     $    56.5e+4, 57.5e+4, 6,
     $    1218025., 1218028.,4,
     $    2053380., 2053390.,0,
     $    6295., 6305.,0,
     $    6360., 6370.,0,
     $    7320., 7325.,0,
     $    7330., 7335.,0,
     $    3725., 3730.,7,
     $    517500., 518500.,8/
      data (elnsv(1,ml,2),elnsv(2,ml,2),nbmpt(ml,2),ml=18,34)/
     $    883300., 883400.,9,
     $    4950., 4965.,10,
     $    5000., 5010.,10,
     $    4360., 4365.,10,
     $    258500., 259500.,0,
     $    127500., 128500.,11,
     $    154500., 157500.,12,
     $    3865., 3870.,13,
     $    3965., 3970.,13,
     $    6713., 6720.,14,
     $    6728., 6733.,14,
     $    4065., 4078.,0,
     $    180000., 197500.,15,
     $    335000., 345000.,16,
     $    9530., 9535.,17,
     $    9065., 9072.,17,
     $    104500., 105500.,18/
      data nlnsv(1)/32/
      data (elnsv(1,ml,1),elnsv(2,ml,1),nbmpt(ml,1),ml=1,16)/
     $    4858., 4863.,14,
     $    1215., 1217.,0,
     $    5875., 5877.,0,
     $    2325., 2327.,0,
     $    1334., 1336.,0,
     $    1908., 1910.,0,
     $    6545., 6550.,2,
     $    6580., 6590.,2,
     $    127500., 128500.,4,
     $    2053380., 2053390.,0,
     $    6295., 6305.,0,
     $    6360., 6370.,0,
     $    7320., 7325.,0,
     $    7330., 7335.,0,
     $    3725., 3730.,3,
     $    517500., 518500.,0/
      data (elnsv(1,ml,1),elnsv(2,ml,1),nbmpt(ml,1),ml=17,32)/
     $    883300., 883400.,0,
     $    4950., 4965.,0,
     $    5000., 5010.,0,
     $    4360., 4365.,0,
     $    258500., 259500.,0,
     $    127500., 128500.,0,
     $    154500., 155500.,0,
     $    3865., 3870.,0,
     $    3965., 3970.,0,
     $    6713., 6720.,5,
     $    6728., 6733.,5,
     $    3.3e+5, 3.4e+5,7,
     $    180000., 197500.,6,
     $    9530., 9535.,8,
     $    9065., 9072.,8,
     $    104500., 105500.,0/
      data nlnsv(3)/32/
      data (elnsv(1,ml,3),elnsv(2,ml,3),nbmpt(ml,3),ml=1,16)/
     $    4858., 4863.,1,
     $    1215., 1217.,0,
     $    5875., 5877.,2,
     $    2325., 2327.,3,
     $    1334., 1336.,4,
     $    1908., 1910.,5,
     $    6545., 6550.,6,
     $    6580., 6590.,6,
     $    1218025., 1218028.,12,
     $    2053380., 2053390.,0,
     $    6295., 6305.,0,
     $    6360., 6370.,0,
     $    7320., 7325.,8,
     $    7330., 7335.,8,
     $    3725., 3730.,9,
     $    517500., 518500.,10/
      data (elnsv(1,ml,3),elnsv(2,ml,3),nbmpt(ml,3),ml=17,32)/
     $    883300., 883400.,0,
     $    4950., 4965.,11,
     $    5000., 5010.,11,
     $    4360., 4365.,11,
     $    258500., 259500.,0,
     $    127500., 128500.,12,
     $    154500., 155500.,13,
     $    3865., 3870.,14,
     $    3965., 3970.,14,
     $    6713., 6720.,0,
     $    6728., 6733.,0,
     $    4065., 4078.,0,
     $    180000., 197500.,15,
     $    9530., 9535.,16,
     $    9065., 9072.,16,
     $    104500., 105500.,17/
      data nlnsv(4)/45/
      data (elnsv(1,ml,4),elnsv(2,ml,4),nbmpt(ml,4),ml=1,14)/
     $   4861., 4863.,1,
     $   6562., 6564., 0,
     $   1215., 1217., 0,
     $   5875., 5877., 2,
     $   4680., 4690., 3,
     $   2325., 2327., 4,
     $   1908., 1910., 5,
     $   1545., 1555., 6,
     $   5195., 5205., 0,
     $   6545., 6550., 7,
     $   6580., 6590., 7,
     $   1745., 1755., 8,
     $   1485., 1490., 10,
     $   1235., 1245.,11/
      data (elnsv(1,ml,4),elnsv(2,ml,4),nbmpt(ml,4),ml=15,30)/
     $   6295., 6305., 12,
     $   3725., 3730., 13,
     $   4950., 4965., 14,
     $   5000., 5010., 14,
     $   517500., 518500.,16, 
     $   4360., 4365., 15,
     $   258500., 259500., 17,
     $   1400., 1408., 18,
     $   1215., 1220., 19,
     $   127500., 128500.,0, 
     $   154500., 155500., 20,
     $   3865., 3870., 21,
     $   3965., 3970., 21,
     $   2420., 2428., 22,
     $   241500., 242500.,0, 
     $   3420., 3430., 23/
      data (elnsv(1,ml,4),elnsv(2,ml,4),nbmpt(ml,4),ml=31,45)/
     $   3340., 3350., 0,
     $   2795., 2800., 25,
     $   44950., 45500., 26,
     $   55950., 56050., 0,
     $   347950., 348050., 27,
     $   2330., 2340., 28,
     $   1890., 1896., 29,
     $   1390., 1405., 30,
     $   6713., 6720., 31,
     $   6728., 6733., 31,
     $   4065., 4078., 0,
     $   186950., 187050.,32, 
     $   9530., 9535., 33,
     $   9065., 9072., 33,
     $   104500., 105500.,34/
      data nlnsv(5)/45/
      data (elnsv(1,ml,5),elnsv(2,ml,5),nbmpt(ml,5),ml=1,14)/
     $   4861., 4863.,1,
     $   6562., 6564., 0,
     $   1215., 1217., 0,
     $   5875., 5877., 2,
     $   4680., 4690., 3,
     $   2325., 2327., 0,
     $   1908., 1910., 4,
     $   1545., 1555., 5,
     $   5195., 5205., 0,
     $   6545., 6550., 6,
     $   6580., 6590., 6,
     $   1745., 1755., 0,
     $   1485., 1490., 0,
     $   1235., 1245.,0/
      data (elnsv(1,ml,5),elnsv(2,ml,5),nbmpt(ml,5),ml=15,30)/
     $   6295., 6305., 0,
     $   3725., 3730., 7,
     $   4950., 4965., 8,
     $   5000., 5010., 8,
     $   517500., 518500.,9, 
     $   4360., 4365., 8,
     $   258500., 259500., 10,
     $   1400., 1408., 0,
     $   1215., 1220., 0,
     $   127500., 128500.,0, 
     $   154500., 155500., 11,
     $   3865., 3870., 12,
     $   3965., 3970., 12,
     $   2420., 2428., 0,
     $   241500., 242500.,0, 
     $   3420., 3430., 0 /
      data (elnsv(1,ml,5),elnsv(2,ml,5),nbmpt(ml,5),ml=31,45)/
     $   3340., 3350., 0,
     $   2795., 2800., 13,
     $   44950., 45500., 0,
     $   55950., 56050., 0,
     $   347950., 348050., 0,
     $   2330., 2340., 0,
     $   1890., 1896., 14,
     $   1390., 1405., 0,
     $   6713., 6720., 0,
     $   6728., 6733., 0,
     $   4065., 4078., 0,
     $   186950., 187050.,15, 
     $   9530., 9535., 16,
     $   9065., 9072., 16,
     $   104500., 105500.,17/
      data nlnsv(6)/45/
      data (elnsv(1,ml,6),elnsv(2,ml,6),nbmpt(ml,6),ml=1,14)/
     $   4861., 4863.,1,
     $   6562., 6564., 0,
     $   1215., 1217., 0,
     $   5875., 5877., 2,
     $   4680., 4690., 3,
     $   2325., 2327., 0,
     $   1908., 1910., 4,
     $   1545., 1555., 5,
     $   5195., 5205., 0,
     $   6545., 6550., 0,
     $   6580., 6590., 0,
     $   1745., 1755., 0,
     $   1485., 1490., 0,
     $   1235., 1245.,0/
      data (elnsv(1,ml,6),elnsv(2,ml,6),nbmpt(ml,6),ml=15,30)/
     $   6295., 6305., 0,
     $   3725., 3730., 0,
     $   4950., 4965., 6,
     $   5000., 5010., 6,
     $   517500., 518500.,7, 
     $   4360., 4365., 6,
     $   258500., 259500., 8,
     $   1400., 1408., 0,
     $   1215., 1220., 0,
     $   127500., 128500.,0, 
     $   154500., 155500., 9,
     $   3865., 3870., 10,
     $   3965., 3970., 10,
     $   2420., 2428., 11,
     $   241500., 242500.,0, 
     $   3420., 3430., 0 /
      data (elnsv(1,ml,6),elnsv(2,ml,6),nbmpt(ml,6),ml=31,45)/
     $   3340., 3350., 0,
     $   2795., 2800., 0,
     $   44950., 45500., 0,
     $   55950., 56050., 0,
     $   347950., 348050., 0,
     $   2330., 2340., 0,
     $   1890., 1896., 0,
     $   1390., 1405., 0,
     $   6713., 6720., 0,
     $   6728., 6733., 0,
     $   4065., 4078., 0,
     $   186950., 187050.,0, 
     $   9530., 9535., 12,
     $   9065., 9072., 12,
     $   104500., 105500.,13/
      data nlnsv(7)/29/
      data (elnsv(1,ml,7),elnsv(2,ml,7),nbmpt(ml,7),ml=1,14)/
     $    4861., 4863., 1,
     $    1215., 1217., 2,
     $    5875., 5877., 3,
     $    4680., 4690., 4,
     $    1635., 1645., 5,
     $    1908., 1910., 6,
     $    1545., 1555., 7,
     $    5195., 5205., 0,
     $    6545., 6550., 8,
     $    6580., 6590., 8,
     $    1745., 1755., 9,
     $    1485., 1490., 10,
     $    6295., 6305., 11,
     $    630000., 635000., 12/
      data (elnsv(1,ml,7),elnsv(2,ml,7),nbmpt(ml,7),ml=15,29)/
     $    3725., 3730., 13,
     $    4950., 4965., 14,
     $    5000., 5010., 14,
     $    4360., 4365., 15,
     $   1400., 1408., 16,
     $   154500., 155500., 17,
     $    3865., 3870., 18,
     $   2420., 2428., 19,
     $   3420., 3430., 20,
     $    2795., 2800., 21,
     $    6713., 6720., 22,
     $    6728., 6733., 22,
     $    9530., 9535., 23,
     $   186950., 187050.,24, 
     $    104500., 105500., 25/
      data kinam(1)/'h i     '/
      data kinam(2),kinam(3)/'he i    ','he ii   '/
      data kinam(4),kinam(5),kinam(6),kinam(7),kinam(8),
     &     kinam(9)/'c i     ','c ii    ','c iii   ','c iv    ',
     &     'c v     ','c vi    '/
      data kinam(10),kinam(11),kinam(12),kinam(13),kinam(14),
     &     kinam(15),kinam(16)/'n i     ','n ii    ','n iii   ',
     &     'n iv    ','n v     ','n vi    ','n vii   '/
      data kinam(17),kinam(18),kinam(19),kinam(20),kinam(21),
     &     kinam(22),kinam(23),kinam(24)/'o i     ','o ii    ',
     &     'o iii   ','o iv    ','o v     ','o vi    ',
     &     'o vii   ','o viii  '/
      data kinam(25),kinam(26),kinam(27),kinam(28),kinam(29),
     &     kinam(30),kinam(31),kinam(32),kinam(33),kinam(34)
     &     /'ne i    ','ne ii   ','ne iii  ','ne iv   ',
     &     'ne v    ','ne vi   ','ne vii  ','ne viii ',
     &     'ne ix   ','ne x    '/
      data kinam(35),kinam(36),kinam(37),kinam(38),kinam(39),
     &     kinam(40),kinam(41),kinam(42),kinam(43),kinam(44),
     &     kinam(45),kinam(46)/'mg i    ','mg ii   ','mg iii  ',
     &     'mg iv   ','mg v    ','mg vi   ','mg vii  ',
     &     'mg viii ','mg ix   ','mg x    ','mg xi   ',
     &     'mg xii  '/
      data kinam(47),kinam(48),kinam(49),kinam(50),kinam(51),
     &     kinam(52),kinam(53),kinam(54),kinam(55),kinam(56),
     &     kinam(57),kinam(58),kinam(59),kinam(60)/'si i    ',
     &     'si ii   ','si iii  ','si iv   ','si v    ',
     &     'si vi   ','si vii  ','si viii ','si ix   ',
     &     'si x    ','si xi   ','si xii  ','si xiii ',
     &     'si xiv  '/
      data kinam(61),kinam(62),kinam(63),kinam(64),kinam(65),
     &     kinam(66),kinam(67),kinam(68),kinam(69),kinam(70),
     &     kinam(71),kinam(72),kinam(73),kinam(74),kinam(75),
     &     kinam(76)/'s i     ','s ii    ','s iii   ','s iv    ',
     &     's v     ','s vi    ','s vii   ','s viii  ',
     &     's ix    ','s x     ','s xi    ','s xii   ',
     &     's xiii  ','s xiv   ','s xv    ','s xvi   '/
      data kinam(77),kinam(78),kinam(79),kinam(80),kinam(81),
     &     kinam(82),kinam(83),kinam(84),kinam(85),kinam(86),
     &     kinam(87),kinam(88),kinam(89),kinam(90),kinam(91),
     &     kinam(92),kinam(93),kinam(94)/'ar i    ','ar ii   ',
     &     'ar iii  ','ar iv   ','ar v    ','ar vi   ',
     &     'ar vii  ','ar viii ','ar ix   ','ar x    ',
     &     'ar xi   ','ar xii  ','ar xiii ','ar xiv  ',
     &     'ar xv   ','ar xvi  ','ar xvii ','ar xviii'/
      data kinam(95),kinam(96),kinam(97),kinam(98),kinam(99),
     &     kinam(100),kinam(101),kinam(102),kinam(103),
     &     kinam(104),kinam(105),kinam(106),kinam(107),
     &     kinam(108),kinam(109),kinam(110),kinam(111),
     &     kinam(112),kinam(113),kinam(114)/'ca i    ',
     &     'ca ii   ','ca iii  ','ca iv   ','ca v    ',
     &     'ca vi   ','ca vii  ','ca viii ','ca ix   ',
     &     'ca x    ','ca xi   ','ca xii  ','ca xiii ',
     &     'ca xiv  ','ca xv   ','ca xvi  ','ca xvii ',
     &     'ca xviii','ca xix  ','ca xx   '/
      data kinam(115),kinam(116),kinam(117),kinam(118),
     &     kinam(119),kinam(120),kinam(121),kinam(122),
     &     kinam(123),kinam(124),kinam(125),kinam(126),
     &     kinam(127),kinam(128),kinam(129),kinam(130),
     &     kinam(131),kinam(132),kinam(133),kinam(134),
     &     kinam(135),kinam(136),kinam(137),kinam(138),
     &     kinam(139),kinam(140)/'fe i    ','fe ii   ',
     &     'fe iii  ','fe iv   ','fe v    ','fe vi   ',
     &     'fe vii  ','fe viii ','fe ix   ','fe x    ',
     &     'fe xi   ','fe xii  ','fe xiii ','fe xiv  ',
     &     'fe xv   ','fe xvi  ','fe xvii ','fe xviii',
     &     'fe xix  ','fe xx   ','fe xxi  ','fe xxii ',
     &     'fe xxiii','fe xxiv ','fe xxv  ','fe xxvi '/
      data kinam(141),kinam(142),kinam(143),kinam(144),
     &     kinam(145),kinam(146),kinam(147),kinam(148),
     &     kinam(149),kinam(150),kinam(151),kinam(152),
     &     kinam(153),kinam(154),kinam(155),kinam(156),
     &     kinam(157),kinam(158),kinam(159),kinam(160),
     &     kinam(161),kinam(162),kinam(163),kinam(164),
     &     kinam(165),kinam(166),kinam(167),kinam(nni),
     &     kinam(169)/'ni i    ','ni ii   ','ni iii  ',
     &     'ni iv   ','ni v    ','ni vi   ','ni vii  ',
     &     'ni viii ','ni ix   ','ni x    ','ni xi   ',
     &     'ni xii  ','ni xiii ','ni xiv  ','ni xv   ',
     &     'ni xvi  ','ni xvii ','ni xviii','ni xix  ',
     &     'ni xx   ','ni xxi  ','ni xxii ','ni xxiii',
     &     'ni xxiv ','ni xxv  ','ni xxvi ','ni xxvii',
     &     'nixxviii','        '/
c
      scfac=1.
      crtt=0.004
      ergsev=1.602197e-12
c
      nbb=nbmk
      lprii=0
      write (6,*)'benchmark number ',nbb
      nlnprnt=nlnsv2(nbb)
      nlnprs=1
      if (nbb.le.2)  nlnprnt=nlnsv2(nbb)+1
      if (nbb.le.2) nlnprs=nlnprnt
      nlnprnt2=nlnsv(nbb)
c
c     first get the current version
      do 820 ll1=1,nlnprnt2
         ll2=ll1
         elnprnt(ll1)=elnsv(1,ll1,nbb)
         elnprnt(ll1+1)=elnsv(2,ll1,nbb)
         if ((elnprnt(ll1).le.1.e-34).or.(elnprnt(ll1+1).le.1.e-34))
     $           go to 820
         ecen=(elnprnt(ll1)+elnprnt(ll1+1))/2.
         ediff=elnprnt(ll1+1)-elnprnt(ll1)
         eww=0.
         ebar=0.
         asym=0.                            
         fluxbs=0.
         fluxbf=0.
         flxmx=0.
         lllsv(ll2)=1
         if (lprii.ne.0) write (6,*)'ll1=',ll1
         do 821 mm=1,nlsvn
            lll=nlsv(mm)
            if (lprii.ne.0)
     $       write (6,*)'ll2,llk,lll=',ll2,ll1,lll
            if (lll.eq.0) go to 821
            if (elin(lll).le.1.e-8) go to 821 
            if ((elin(lll).gt.elnprnt(ll1+1)).or.
     $          (elin(lll).lt.elnprnt(ll1))) go to 821
            ener=12398.54/elin(lll)  
            tmpflxf=elum(lll)
     $          *(1.e+19/r)*(1.e+19/r)*ener*ergsev/12.56       
            fluxbf=fluxbf+tmpflxf
            tmpflux=(elum(lll)+elumb(lll))
     $          *(1.e+19/r)*(1.e+19/r)*ener*ergsev/12.56       
            if (tmpflux.lt.flxmx) go to 3083
              flxmx=tmpflux
              lllsv(ll2)=lll
3083          continue
            fluxbs=fluxbs+tmpflux
            if (ll1.eq.1) scfac=tmpflux
            fluxrl(ll2)=tmpflux/(scfac+1.e-18)
            fluxrl(ll2)=min(fluxrl(ll2),1.e+6)        
            frac=(elum(lll)+elumb(lll))*ener*ergsev/(etotc+1.e-18)
            asym=elum(lll)/amax1(1.e-34,elum(lll)+elumb(lll))
            nbltp=nblin(lll)
            ewtmp=0.
            if ((nbltp.gt.0).and.(nbltp.lt.numcon))
     $       ewtmp=(elum(lll)+elumb(lll))*ener/
     $         amax1(1.e-34,zremsb(nbltp)+zrems(nbltp))
            eww=eww+ewtmp
            ebar=ebar+tmpflux*elin(lll)
821         continue
         if (ll1.eq.1) scfac=fluxbs
         fluxrl(ll2)=fluxbs/(scfac+1.e-18)                                   
         fluxrl(ll2)=min(fluxrl(ll2),1.e+6)   
         if (lprii.ne.0)
     $     write (6,*)'ll2,fluxrl:',ll2,fluxrl(ll2)
         asym=fluxbf/amax1(1.e-34,fluxbs)
         ebar=ebar/amax1(fluxbs,1.e-34)     
         r19=r*1.e-19
         eltmp=fluxbs*12.56*r19*r19 
820      continue                                                      
c
c     put them together and print out
c     step through the benchmark lines
      do 1039 ll1=1,nlnprnt
         ll2=ll1
         fffrl=0.
c        step through the printed lines
         do 1040 ll3=1,nlnprnt2
            if (nbmpt(ll3,nbb).ne.ll2) go to 1040
            fffrl=fluxrl(ll3)+fffrl
            ll3sv=ll3
 1040       continue
         if (fffrl.le.1.e-34) go to 1039
         fdd=1.
         if (nbb.eq.2) fdd=scfac*12.56*r19*r19*10.
         if (nbb.eq.1) fdd=scfac*12.56*r19*r19*100.
         if (nbb.eq.4) fdd=scfac*12.56*r19*r19*1000.
         if (nbb.eq.5) fdd=scfac*12.56*r19*r19*10000.
         if (nbb.eq.6) fdd=scfac*12.56*r19*r19*10000.
         if (nbb.eq.7) fdd=scfac
         if (nbb.eq.3) fdd=scfac
         if (ll1.eq.nlnprs) fffrl=fffrl*fdd
         bmdat(6,ll2,nbb)=fffrl
 1039    continue
c
c     set up averages over the lexington results
      nlnprnt=nlnsv(nbb)
      do 1013 ll1=1,nlnprnt
         ll2=ll1
         fbar=0.
         nbar=0
         if (lprii.eq.1) write (6,*)'line number ',ll1,ll2
         do 1011 ml1=3,nbmdat(nbb)
            if (bmdat(ml1,ll2,nbb).le.1.e-34) go to 1011
c            if (ml1.eq.6) go to 1011
            fbar=fbar+bmdat(ml1,ll2,nbb)
            nbar=nbar+1
            if (lprii.eq.1) 
     $        write (6,*)'table data',ml1,bmdat(ml1,ll2,nbb),fbar
 1011       continue
         fbar=fbar/float(max0(1,nbar))
         if (nbar.le.0) fbar=0.
         fbsv(ll2)=fbar
         sigsum=0.
         nbar=0
         do 1014 ml1=3,nbmdat(nbb)
            if (bmdat(ml1,ll2,nbb).le.1.e-34) go to 1014
c            if (ml1.eq.6) go to 1014
            sigsum=sigsum+(bmdat(ml1,ll2,nbb)-fbar)**2
            nbar=nbar+1
            if (lprii.eq.1) 
     $        write (6,*)'table disp',ml1,bmdat(ml1,ll2,nbb),
     $        fbar,sigsum
 1014       continue
         sig(ll2)=sqrt(sigsum/float(max0(nbar,1)))
         if (nbar.le.0) sig(ll2)=0.
         if (lprii.eq.1) write (6,*)'average and disp:',
     $               ll2,fbsv(ll2),sig(ll2)
 1013    continue
c
c
c
      nlnprnt=nlnsv2(nbb)
      nlnprs=1
      if (nbb.le.2)  nlnprnt=nlnsv2(nbb)+1
      if (nbb.le.2) nlnprs=nlnprnt
      nlnprnt2=nlnsv(nbb)
      if (lprii.ne.0)
     $  write (6,*)'nlnprnt,nlnprnt2:',nlnprnt,nlnprnt2
c     step through the benchmark lines
      do 1019 ll1=1,nlnprnt
         ll2=ll1
         fffrl=0.
c        step through the printed lines
         do 1020 ll3=1,nlnprnt2
            if (nbmpt(ll3,nbb).ne.ll2) go to 1020
            fffrl=fluxrl(ll3)+fffrl
            if (lprii.ne.0)
     $       write (6,*)'stepping through constituents:',
     $       ll3,lllsv(ll3),elin(lllsv(ll3)),kinam(nilin(lllsv(ll3))),
     $              fluxrl(ll3)
            ll3sv=ll3
 1020       continue
         if (lprii.ne.0)
     $    write (6,*)'fffrl=',ll1,fffrl,ll3sv
         if (fffrl.le.1.e-34) go to 1019
         fdd=1.
         if (nbb.eq.2) fdd=scfac*12.56*r19*r19*10.
         if (nbb.eq.1) fdd=scfac*12.56*r19*r19*100.
         if (nbb.eq.4) fdd=scfac*12.56*r19*r19*1000.
         if (nbb.eq.5) fdd=scfac*12.56*r19*r19*10000.
         if (nbb.eq.6) fdd=scfac*12.56*r19*r19*10000.
         if (nbb.eq.7) fdd=scfac
         if (nbb.eq.3) fdd=scfac
         if (ll1.eq.nlnprs) fffrl=fffrl*fdd
         nbmtmp=nbmdat(nbb)+1
         do 1012 ml1=3,nbmtmp
c            if (ml1.eq.6) go to 1012
            if (ml1.eq.nbmtmp) go to 1015
            if (bmdat(ml1,ll2,nbb).le.1.e-34) go to 1012
            err(ml1)=(bmdat(ml1,ll2,nbb)-fbsv(ll2))
     $                 /amax1(1.e-34,sig(ll2))
            if (lprii.ne.0)
     $       write (6,*)'ml1),bmdat,err:',
     $          ml1,bmdat(ml1,ll2,nbb),fbsv(ll2),sig(ll2),err(ml1)
            go to 1016
 1015       continue
            err(ml1)=(fffrl-fbsv(ll2))
     $                 /amax1(1.e-34,sig(ll2))
            if (lprii.ne.0)
     $       write (6,*)'ml1,bmdat,err:',
     $          ml1,fffrl,fbsv(ll2),sig(ll2),err(ml1)
 1016       continue
            sigsm(ml1)=sigsm(ml1)+err(ml1)*err(ml1)
            err(ml1)=min(amax1(err(ml1),-9.99),9.99)
            nsm(ml1)=nsm(ml1)+1
 1012       continue
          write (6,9901) ll1,klbm(ll1,nbb),eldat(ll1,nbb),
     $    err(nbmtmp),fbsv(ll2),sig(ll2),fffrl
 9901     format (1h ,i4,1x,a8,f10.2,
     $      f5.2,1x,3(1pe10.2))
 1019    continue
      write (6,*)(sigsm(mm),mm=1,nbmtmp)
c
      return
      end
      blockdata bigdat
c
      parameter (ncn=5000,nnnl=3900,nnml=4100,nl=13,nni=168,nnip=nni+nl)
c
c
c
c
      common /aidata/ caih(1),caihe(2),caic(6),cain(7),caio(8),
     &                caine(10),caimg(12),caisi(14),cais(16),
     &                caiar(18),caica(20),caife(26),caini(28),
     &                eaih(1),eaihe(2),eaic(6),eain(7),eaio(8),
     &                eaine(10),eaimg(12),eaisi(14),eais(16),
     &                eaiar(18),eaica(20),eaife(26),eaini(28)
      common /enerc / epi(ncn),dele(ncn),numcon
      common /cxdat1/ alhh(2,1),alhhe(2,2),alhc(2,6),alhn(2,7),
     &                alho(2,8),alhne(2,10),alhmg(2,12),
     &                alhsi(2,14),alhs(2,16),alhar(2,18),
     &                alhca(2,20),alhfe(2,26),alhni(2,28),bethh(1)
     &               ,bethhe(2),bethc(6),bethn(7),betho(8),
     &                bethne(10),bethmg(12),bethsi(14),beths(16),
     &                bethar(18),bethca(20),bethfe(26),bethni(28)
     &               ,hkihh(1),hkihhe(2),hkihc(6),hkihn(7),
     &                hkiho(8),hkihne(10),hkihmg(12),hkihsi(14),
     &                hkihs(16),hkihar(18),hkihca(20),hkihfe(26),
     &                hkihni(28)
      common /cxdat2/ aleh(2,1),alehe(2,2),alec(2,6),alen(2,7),
     &                aleo(2,8),alene(2,10),alemg(2,12),
     &                alesi(2,14),ales(2,16),alear(2,18),
     &                aleca(2,20),alefe(2,26),aleni(2,28),beteh(1)
     &               ,betehe(2),betec(6),beten(7),beteo(8),
     &                betene(10),betemg(12),betesi(14),betes(16),
     &                betear(18),beteca(20),betefe(26),beteni(28)
     &               ,hkieh(1),hkiehe(2),hkiec(6),hkien(7),
     &                hkieo(8),hkiene(10),hkiemg(12),hkiesi(14),
     &                hkies(16),hkiear(18),hkieca(20),hkiefe(26),
     &                hkieni(28)
      common /rrdat1/ aradh(1),aradhe(2),aradc(6),aradn(7),
     &                arado(8),aradne(10),aradmg(12),aradsi(14),
     &                arads(16),aradar(18),aradca(20),aradfe(26),
     &                aradni(28),etah(1),etahe(2),etac(6),
     &                etan(7),etao(8),etane(10),etamg(12),
     &                etasi(14),etas(16),etaar(18),etaca(20),
     &                etafe(26),etani(28)
      common /rrdat2/ adih(1),adihe(2),adic(6),adin(7),adio(8),
     &                adine(10),adimg(12),adisi(14),adis(16),
     &                adiar(18),adica(20),adife(26),adini(28),
     &                bdih(1),bdihe(2),bdic(6),bdin(7),bdio(8),
     &                bdine(10),bdimg(12),bdisi(14),bdis(16),
     &                bdiar(18),bdica(20),bdife(26),bdini(28),
     &                t0h(1),t0he(2),t0c(6),t0n(7),t0o(8),
     &                t0ne(10),t0mg(12),t0si(14),t0s(16),
     &                t0ar(18),t0ca(20),t0fe(26),t0ni(28),t1h(1)
     &               ,t1he(2),t1c(6),t1n(7),t1o(8),t1ne(10),
     &                t1mg(12),t1si(14),t1s(16),t1ar(18),
     &                t1ca(20),t1fe(26),t1ni(28)
      common /rrdat3/ edilh(1),edilhe(2),edilc(6),ediln(7),
     &                edilo(8),edilne(10),edilmg(12),edilsi(14),
     &                edils(16),edilar(18),edilca(20),edilfe(26),
     &                edilni(28),aph(1),aphe(2),apc(6),apn(7),
     &                apo(8),apne(10),apmg(12),apsi(14),aps(16)
     &               ,apar(18),apca(20),apfe(26),apni(28),
     $                cdd1(57),cdd2(111),
     &                ndilh(1),ndilhe(2),ndilc(6),ndiln(7),
     &                ndilo(8),ndilne(10),ndilmg(12),ndilsi(14),
     &                ndils(16),ndilar(18),ndilca(20),ndilfe(26),
     &                ndilni(28)
      common /ethrsh/ ethh(1),ethhe(2),ethc(6),ethn(7),etho(8),
     &                ethne(10),ethmg(12),ethsi(14),eths(16),
     &                ethar(18),ethca(20),ethfe(26),ethni(28)
      common /abel  / xeh,xehe,xec,xen,xeo,xene,xemg,
     &                xesi,xes,xear,xeca,xefe,xeni,xeln(nni)
     &               ,xelln(nnnl)
      common /swdata/ stwtrt(nni)
      common /ethph / etkh(1),etkhe(2),etkc(6),etlsc(4),etlpc(2)
     &               ,etkn(7),etlsn(5),etlpn(3),etko(8),
     &                etlso(6),etlpo(4),etkne(10),etlsne(8),
     &                etlpne(6),etkmg(12),etlsmg(10),etlpmg(8),
     &                etmsmg(2),etksi(14),etlssi(12),etlpsi(10),
     &                etmssi(4),etmpsi(2),etks(16),etlss(14),
     &                etlps(12),etmss(6),etmps(4),etkar(18),
     &                etlsar(16),etlpar(14),etmsar(8),etmpar(6),
     &                etkca(20),etlsca(18),etlpca(16),etmsca(10),
     &                etmpca(8),etmdca(2),etkfe(26),etlsfe(24),
     &                etlpfe(22),etmsfe(16),etmpfe(14),etmdfe(8),
     &                etnfe(2),etkni(28),etlsni(26),etlpni(24),
     &                etmsni(18),etmpni(16),etmdni(10),etnni(4)
      common /rchdat/ acor1(15,13),acor2(15,13),nup1(15,13),
     &                nup2(15,13),prfrac(26),ngnd(26),nz(nni),
     &                cor2(26),nnnz(nni)
      common /fe2dat/ ggl(10),ggu(6),eel(10),eeu(6),b(10,6),
     &                e(10,6),nblfe2(10,6),f(10,6),g(10,6)
      common /agdata/ amkc(6,6),amkn(7,7),amko(8,8),amkne(10,10),
     &                amkmg(12,12),amlsmg(10,10),amlpmg(8,8),
     &                amksi(14,14),amlssi(12,12),amlpsi(10,10),
     &                amks(16,16),amlss(14,14),amlps(12,12),
     &                amkf1(26,13),amkf2(26,13),amlsf1(24,12),
     &                amlsf2(24,12),amlpf1(22,11),amlpf2(22,11)
      common /lab   / kllab,kulab,k
      common /itdat / enfmxs,ensfrc,critd,crittd,epss,crits,
     $             kmaxs,lppris,nlimd,lpprid,nnmax
c
      character*8 kllab(10),kulab(6),k(10,6)
c
c      block data
      dimension bln(1000),nj(13),abunj(13),conce(30)
      dimension gndrec(30),power(220),rf(500),tau(nni)
      dimension tauhe(nni),tplus(nni),cnc(12,30),pm(4)
      dimension ptot(12,220),abin(1000),bin(1000),crit(30,6)
      common /bln   / bln,blnmin,blnsyz,nbln
      common /params/ abunj,abund,binmin,binsyz,nbin,nj
      common /result/ conce,gndrec,power,rhy,heneut,heplus,
     &                dne,pcool,pou,pot,re,tu,pm
      common /pt    / rf,tau,tauhe,tplus
      common /com   / cnc,ptot,abin,bin
      common /crt   / crit
      data crit/5*0.,555.,1953.,5400.,0.,22000.,0.,
     &     65000.,0.,140000.,0.,270000.,0.,1.6e6,0.,
     &     2.9e6,5*0.,1.8e7,0.,3.5e7,0.,0.,5*0.,5.8e7,
     &     9.5e8,6.7e9,0.,1.2e11,0.,1.1e12,0.,6.8e12,
     &     0.,3.1e13,0.,1.1e14,0.,3.4e14,5*0.,9.1e15,
     &     0.,2.3e16,0.,0.,5*0.,7.1e8,4.4e9,2.4e10,0.,
     &     3.1e11,0.,2.7e12,0.,1.9e13,0.,6.7e13,0.,
     &     2.1e14,0.,5.8e14,5*0.,6.5e15,0.,1.9e16,0.,
     &     0.,5*0.,5.8e7,1.5e8,1.7e9,0.,4.3e10,0.,
     &     3.9e11,0.,2.3e12,0.,1.0e13,0.,4.5e13,0.,
     &     1.4e14,5*0.,2.6e15,0.,7.5e15,0.,0.,11*0.,
     &     174.,0.,877.,0.,12700.,0.,64000.,0.,230000.,
     &     5*0.,8.5e5,0.,2.7e6,0.,0.,11*0.,3.3e9,0.,
     &     6.1e10,0.,1.2e12,0.,1.4e13,0.,8.8e13,5*0.,
     &     4.2e15,0.,1.0e16,0.,0./
      data rf/500*1.0/
      data abunj/10.93,8.52,7.96,8.82,7.92,7.42,7.52,7.2,
     &     6.9,6.3,7.6,6.3,0./
c
      data nlimd/100/,lpprid/0/,critd/1.e-05/,crittd/1.e-06/
      data nnmax/40/,epss/0.07/,crits/1.e-4/,kmaxs/0/,
     &     enfmxs/1.e-4/,ensfrc/0.5/,lppris/0/
c
      data xeh/1./,xehe/0.1/,xec/3.7e-4/,xen/1.1e-4/,
     &     xeo/6.8e-4/,xene/2.8e-5/,xemg/3.5e-5/,xesi/3.5e-5/,
     &     xes/1.6e-5/,xear/4.5e-06/,xeca/2.1e-6/,xefe/2.5e-5/,
     &     xeni/2.e-6/
c
c
      data numcon/299/
c
c     ionization thresholds
      data ethh/13.598/
      data ethhe/24.59,54.42/
      data ethc/11.260,24.383,47.887,64.492,392.08,489.98/
      data ethn/14.534,29.601,47.448,77.472,97.89,552.06,
     &     667.03/
      data etho/13.618,35.117,54.934,77.413,113.9,138.12,
     &     739.32,871.39/
      data ethne/21.564,40.962,63.45,97.11,126.21,157.93,
     &     207.26,239.09,1195.8,1362.2/
      data ethmg/7.646,15.035,80.143,109.31,141.27,186.51,
     &     224.95,265.92,328.0,367.5,1761.8,1963./
      data ethsi/8.151,16.345,33.492,45.141,166.77,205.08,
     &     246.49,303.16,351.1,401.4,476.1,523.,2438.,
     &     2673./
      data eths/10.36,23.33,34.83,47.30,72.68,88.05,280.01,
     &     328.33,379.1,447.1,504.7,565.,652.,707.,3224.,
     &     3494./
      data ethar/15.759,27.629,40.74,59.81,75.04,91.01,
     &     124.4,143.45,422.6,478.9,539.0,618.,686.,756.,
     &     855.,918.,4121.,4426./
      data ethca/6.113,11.871,50.91,67.15,84.43,108.78,
     &     127.7,147.4,188.7,211.3,591.6,657.,726.,817.,
     &     895.,974.,1087.,1157.,5129.,5470./
      data ethfe/07.87,16.16,30.65,54.8,75.5,100.,128.3,
     &     151.12,235.,262.1,290.4,331.,361.,392.,457.,
     &     489.,1262.,1358.,1456.,1582.,1689.,1799.,
     &     1950.,2045.,8828.,9278./
      data ethni/7.635,18.168,35.17,54.9,75.5,108.,134.,
     &     164.,193.,224.6,321.,352.,384.,430.,464.,
     &     499.,571.,607.,1546.,1648.,1756.,1894.,2010.,
     &     2131.,2295.,2398.,10280.,10790./
c
c
      data etkh(1)/1.360E+01/
      data etkhe(1)/2.459E+01/
      data etkhe(2)/5.442E+01/
      data etlpc(1)/1.126E+01/
      data etlsc(1)/1.939E+01/
      data etkc(1)/2.910E+02/
      data etlpc(2)/2.438E+01/
      data etlsc(2)/3.047E+01/
      data etkc(2)/3.076E+02/
      data etlsc(3)/4.789E+01/
      data etkc(3)/3.289E+02/
      data etlsc(4)/6.449E+01/
      data etkc(4)/3.522E+02/
      data etkc(5)/3.921E+02/
      data etkc(6)/4.900E+02/
      data etlpn(1)/1.453E+01/
      data etlsn(1)/2.541E+01/
      data etkn(1)/4.048E+02/
      data etlpn(2)/2.960E+01/
      data etlsn(2)/3.796E+01/
      data etkn(2)/4.236E+02/
      data etlpn(3)/4.745E+01/
      data etlsn(3)/5.545E+01/
      data etkn(3)/4.473E+02/
      data etlsn(4)/7.747E+01/
      data etkn(4)/4.753E+02/
      data etlsn(5)/9.789E+01/
      data etkn(5)/5.043E+02/
      data etkn(6)/5.521E+02/
      data etkn(7)/6.671E+02/
      data etlpo(1)/1.362E+01/
      data etlso(1)/2.848E+01/
      data etko(1)/5.380E+02/
      data etlpo(2)/3.512E+01/
      data etlso(2)/4.599E+01/
      data etko(2)/5.581E+02/
      data etlpo(3)/5.494E+01/
      data etlso(3)/6.551E+01/
      data etko(3)/5.840E+02/
      data etlpo(4)/7.741E+01/
      data etlso(4)/8.737E+01/
      data etko(4)/6.144E+02/
      data etlso(5)