*ident upiaea1
*/ groupr - A.Trkov, IAEA
*/          Local updates from NJOY97 that were not implemented
*/ Need more space for the flux calculator (Case: U-238)
*d up89.6
      dimension a(400000)
*d up89.8
      iamax=400000
*/
*ident upijs57
*/ groupr  A.Trkov, 9-March-1999
*/ groupr is caught in an infinite loop for Be-9 from ENDF/B-VI
*/ on a PC with Lahey compiler.
*/ It seems the products inside an "if" statements are
*/ computed differently than outside. The program does not
*/ recognise that ep EQUALS "test" when written as a product.
*d groupr.6054
      test=shade*epn
      if (idis.gt.0.and.ep.lt.test) epn=test
*/
*/ Old proposal from C.Dean
*d groupr.2745,2746
               b(iz+3+li)=(sigz(iz)-sam)*wtf*(1-beta)
*/
*ident upijs59
*/ acer   A.Trkov 20-aug-1999
*/ If pointwise representation in CM system, csn should be used
*/ (pointed out by Harry Wienke).
*d acer.3394,3396
                        if (csn.ge.a(ll).and.csn.le.a(ll+2))
     &                    call terp1(a(ll),a(ll+1),a(ll+2),a(ll+3),
     &                    csn,fmu,lang-10)
*/
*ident upijs61
*/ acer   A.Trkov 3-dec-1999
*/        fix6 may run out of space without warning
*d acer.3247
      dimension a(2000)
*d up3.38
      data namax/2000/
*i acer.3262
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3272
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3278
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3280
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3282
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3284
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3290
      if(nw.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3296
      if(nw+l.gt.namax) call error('fix6','storage exceeded.',' ')
*i acer.3299
      if(nw+l.gt.namax) call error('fix6','storage exceeded.',' ')
*/
*ident upijs62
*/ groupr  A.Trkov, 29-Apr-2002
*/         - Fix inconsistent usage of output weighting flux unit number
*/         - Fix logic when searching for the right flux point
*i groupr.80
c     *            note: weighting flux file is always written binary  *
*d groupr.99
c     *    ninwt   tape unit for flux parameters (binary)              *
*d groupr.272
*d groupr.2290
         ninwt=iabs( ninwt)
         call openz(-ninwt,1)
*d groupr.2294
     &     ehi,sigpot,nflmax,-ninwt,jsigz
*d groupr.2307,2308
         call openz(-ninwt,0)
         write(nsyso,'(/'' ninwt......'',i4)') -ninwt
*d groupr.3037
      if (e.gt.el*(1-small).and.e.lt.en*(1+small)) go to 230
*/
*ident upijs63
*/ purr  A.Trkov, May 2002
*/       Increase scratch space to process Pu-239 from JEF-2.2
*d purr.95
      maxscr=12000
*/
*ident upijs64
*/ groupr  A.Trkov, July 2002
*/         There is an error in Cl-37 in ENDF/B-VI Rev.8
*/         MAT 1731 MF 15 MT 102
*/         The first energy for the gamma yields is zero. It
*/         redefines enext in getsed to zero, which is
*/         interpreted as a flag to do initialisation in
*/         getgyl (called from getff). Since 'yl' is
*/         already reserved, groupr crashes.
*d groupr.8697
            if (elo.gt.0 .and. elo.lt.enext) enext=elo
*/
*ident upiaea2
*/ groupr - A.Trkov, IAEA, Feb-2003
*/          Increase array size for Legendre coefficients
*i groupr.4756
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.4764
      dimension term(mxlg),terml(mxlg)
*i groupr.5209
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.5212
      dimension term(mxlg),x(10),y(10,mxlg)
*d groupr.5329,5330
c     maximum legendre coefficients
      parameter (mxlg=65)
      dimension cnow(*),term(*),p(mxlg)
      dimension x(10),y(10,mxlg),yt(mxlg)
*d groupr.5333
      external f6ddx,f6psp,f6dis,legndr,error
*i groupr.5358
      if(nl.gt.mxlg) call error('f6cm','nl>mxlg',' ')
*d groupr.5586
c     maximum legendre coefficients
      parameter (mxlg=65)
      dimension cnow(*),p(mxlg)
*d groupr.5768
c     maximum legendre coefficients
      parameter (mxlg=65)
      dimension cnow(*),p(mxlg)
*i groupr.5630
         if(nl.gt.mxlg) call error('f6ddx','nl>mxlg',' ')
*i groupr.5992
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.5995
      dimension term(mxlg),p(mxlg),amu(50),fmu(50),qp(8),qw(8)
*i groupr.6130
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.6133,6134
      dimension term1(mxlg),term2(mxlg),p(mxlg)
      dimension qp(8),qw(8)
*i groupr.6733
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.6740
      dimension flo(mxlg),fhi(mxlg)
*i groupr.7314
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.7318
      dimension b(6),alo(mxlg),ahi(mxlg)
*i groupr.7929
c     maximum legendre coefficients
      parameter (mxlg=65)
*d groupr.7942
      dimension fl(mxlg)
*/
*ident upiaea3
*/ acer   - A.Trkov, IAEA, Feb-2003
*/          declare "error" external to avoid conflict with intrinsic function
*i up69.66
      external error
*/
*ident upfndl1
*/ matxsr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ Need more space  (for processing JENDL-3.3)
*/  subroutine band
*d matxsr.1973
      dimension b(30000)
*d matxsr.1978
      maxb=30000
*/  subroutine shufl
*d matxsr.2078
      dimension b(30000)
*d matxsr.2080
      maxb=30000
*/
*ident upfndl2
*/ gaspr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ Need more space (case Ni-058 from JEFF-3.0)
*d gaspr.28
      dimension egas(80000),sgas(5,80000)
*d gaspr.41
      maxg=80000
*ident upfndl3
*/ purr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ to allow PT for multi-isotope materials when
*/ not all the isotopes have unresolved resonance data
*d purr.624
*i purr.651
      if (ier.eq.ner) go to 110
*/
*ident upfndl4
*/ acer - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ Format extension for negative energies (overlap)
*d acer.4947,4948
            write(nsyso,'(''   energy range: '',1p,e11.4,
     &        '' - '',e11.4,'' ev'')') urlo,urhi
*/
*ident upfndl5
*/ groupr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ Needed to process Be-9 from JEFF-3.0 evaluation
*/  MT=875-891 for (n,2n) splitted reaction
*i groupr.1147
         else if (mtd.ge.875.and.mtd.le.884) then
            write(reac,'(''2n0'',i1)') mtd-875
         else if (mtd.ge.885.and.mtd.le.890) then
            write(reac,'(''2n'',i2)') mtd-875
         else if (mtd.eq.891) then
            reac='2nc'
*i groupr.3964
         if (iverf.ge.6.and.mtd.ge.875.and.mtd.le.891) mt=mtd
*i groupr.4331
      if (mfd.ge.31.and.mfd.le.36.and.iverf.ge.6.and.
     &   (mtd.ge.875.and.mtd.lt.891)) go to 400
*i groupr.8014
         if (mth.ge.875.and.mth.lt.891) za2=1
*i groupr.8330
      if (iverf.ge.6.and.mth.ge.875.and.mth.le.890) mt0=875
*/
*/ matxsr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ Needed to process Be-9 from JEFF-3.0 evaluation
*/  MT=875-891 for (n,2n) splitted reaction
*d matxsr.1418
  290 if (mt.gt.891) go to 300
*i matxsr.1433
      if (mt.ge.875.and.mt.lt.885) write(strng,'(''2n0'',i1)') mt-875
      if (mt.ge.885.and.mt.lt.891) write(strng,'(''2n'',i2)') mt-875
      if (mt.eq.891) write(strng,'(''2ncn'')')
*i matxsr.1504
      k016=0
*i matxsr.1652
      if (mt.eq.16) k016=1
      if (k016.eq.1.and.mt.ge.875.and.mt.le.891) go to 310
*/
*ident upfndl6
*/ groupr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ subroutine namer corrected for mt=659, reac=d09
*/ It only affects groupr printout
*d groupr.1124
         else if (mtd.ge.650.and.mtd.le.659) then
*/
*ident wlup2
*/ thermr - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/ More space for new scattering law
*d thermr.101
      dimension a(800000)
*d thermr.131
      namax=800000
*/
*/ subroutine calcem - D. L. Aldama, NDS/IAEA Consultant, Nov-2004
*/  TAB1 records are larger than npage=306
*/  moreio should be used
*/  required by thermal scattering law generated by NDS/IKE
*/
*d thermr.1543
      lmore=loc
      call tab1io(nendf,0,0,a(loc),nb,nw)
      do while (nb.ne.0)
        lmore=lmore+nw
        call moreio(nendf,0,0,a(lmore),nb,nw)
      enddo
*d thermr.1572
         lmore=loc
         call listio(nendf,0,0,a(loc),nb,nw)
         do while (nb.ne.0)
           lmore=lmore+nw
           call moreio(nendf,0,0,a(lmore),nb,nw)
         enddo
*d thermr.1598
         lmore=loc
         call listio(nendf,0,0,a(loc),nb,nw)
         do while (nb.ne.0)
           lmore=lmore+nw
           call moreio(nendf,0,0,a(lmore),nb,nw)
         enddo
*d thermr.1602
      lmore=loc
      call tab1io(nendf,0,0,a(loc),nb,nw)
      do while (nb.ne.0)
        lmore=lmore+nw
        call moreio(nendf,0,0,a(lmore),nb,nw)
      enddo
*/
*/ident up2_thermr
*/ thermr - M. Mattes/IKE ---------------------------------------mm
*/    maximum.energy transfer at higher temperatures (print only)
*i thermr.1622
      if (lat.eq.1) tmax=tmax*tevz/(bk*temp)
*/ ---------------------------------------------------end up2_thermr
*ident up3_thermr
*/   more incident energies necessary for ZrH for an adequate
*/   represention of the neutron scattering xs in comparison
*/   with experimental results (.1 < E_i < 1 eV)
*d thermr.1422,thermr.1423
       dimension egrid(79)
       dimension ubar(79)
*d thermr.1427
       data ngrid/79/
*d thermr.1433,thermr.1436
     &  .030613d0,.042757d0,.056925d0,.081972d0,
     *      .1d0, .111573d0, .12d0, .13d0, .145728d0, .16d0, .17d0,
     &  .184437d0,
     *    .195d0, .2277d0, .24d0, .2510392d0,.2705304d0,.2907501d0,
     &  .3011332d0,.3206421d0,.3576813d0,
     *   .38d0, .4d0, .4170351d0,  .44d0, .48d0, .5032575d0, .53d0,
     &   .56d0, .59d0,  .6249328d0, .65d0, .69d0, .72d0, .75d0,
     *   .7821141d0, .82d0, .89d0, .9506956d0,1.0137432d0,1.1664337d0,
*d thermr.1455,thermr.1458
     &  .030613e0,.042757e0,.056925e0,.081972e0,
     *      .1e0, .111573e0, .12e0, .13e0, .145728e0, .16e0, .17e0,
     &  .184437e0,
     *    .195e0, .2277e0, .24e0, .2510392e0,.2705304e0,.2907501e0,
     &  .3011332e0,.3206421e0,.3576813e0,
     *   .38e0, .4e0, .4170351e0,  .44e0, .48e0, .5032575e0, .53e0,
     &   .56e0, .59e0,  .6249328e0, .65e0, .69e0, .72e0, .75e0,
     *   .7821141e0, .82e0, .89e0, .9506956e0,1.0137432e0,1.1664337e0,
*/
*/ --------------------------------------------------end up3_thermr
*ident up1_leapr
*/
*/ if number(>200) for alpha and/or beta values then
*/           increase array sizes plus working buffers
*/    e.g. beta(200) -> beta(400) for H2O
*d leapr.164
       common/ab/nalpha,nbeta,naint,nbint,alpha(200),beta(400)
*d leapr.173
       common/lstore/a(7500000)
*d leapr.178
       data nbmax,namax/400,200/
*d leapr.189
       maxa=7500000
*d leapr.391
       mscr=4000
*d leapr.420
       common/ab/nalph1,nbeta1,naint,nbint,alph1(200),beta1(400)
*d leapr.426
       dimension maxt(400)
*d leapr.794
       common/ab/nalpha,nbeta,naint,nbint,alpha(200),beta1(400)
*d leapr.1252
       common/ab/nalpha1,nbeta1,naint,nbint,alpha(200),beta1(400)
*d leapr.1852
       common/ab/nalph1,nbeta1,naint,nbint,alpha(200),beta1(400)
*/
*/    IKE evaluation of H(H2O) creates more pairs of alpha
*/    and S(alpha,beta) than npage
*i leapr.3122
               l_mm=1+nw
               do while (nb.ne.0)
                    call moreio(0,nout,nprnt,scr(l_mm),nb,nw)
                    l_mm=l_mm+nw
               enddo
*/ -------------------------------------------------end up1_leapr
*ident up2_leapr
*/     for H(H2O), D(D2O) NS should be equal 1
*d leapr.2976,2977
      if(nss.gt.0) scr(5)=6*(nss+1)
      scr(6)=nss
*/ -------------------------------------------------end up2_leapr
*ident up3_leapr
*/  T-eff calculation
*d leapr.1557
       tempf(itemp)=(tbeta+twt)*tempf(itemp)+tsave
*/ -------------------------------------------------end up3_leapr
*ident up4_leapr
*/ generation of more than 10 temperatures
*d leapr.166,167
      common/te/tempr(20),tempf(20),tempf1(20)
      common/dw/dwpix(20),dwp1(20)
*d leapr.179
      data ntmax/20/
*d leapr.600,601
      common/te/tempr(20),tempf(20),tempf1(20)
      common/dw/dwpix(20),dwp1(20)
*d leapr.797
      common/te/tempr(20),tempf(20),tempf1(20)
*d leapr.1256
      common/te/tempr(20),tempf(20),tempf1(20)
*d leapr.1258
      common/dw/dwpix(20),dwp1(20)
*d leapr.1856
      common/te/tempr(20),tempf(20),tempf1(20)
*d leapr.2651,2652
      common/te/tempr(20),tempf(20),tempf1(20)
      common/dw/dbw(20),dbw1(20)
*/ --------------------------------------------------end up4_leapr
*ident up5_leapr
*/   correct directory
*i leapr.2797
         scr(5)=scr(5)+1
      if(iel.ne.0) scr(5)=scr(5)+1
*/ --------------------------------------------------end up5_leapr
*/
*ident up_acer
*/ id-name of thermal data set with 6 characters else blank
*i acer.358
         if (nch.eq.0) nch=6
*/
*/ more than 64 bins  wt(65)
*d acer.13079
      dimension wt(401)
*/
*/ increase of scratch buffer (?)
*d acer.13098
      ninmax=8000
*/
*/ ------------------------------------------------jan05------end
*/
*/ Kazuaki Kosako, Shimizu Corporation
*/ MF5 and MF6 INT=21-25 laws modifications
*/ in groupr
*/
*ident upshim
*/
*i groupr.548
      ee=0
*d groupr.550
*/
*/*d groupr.5586
*/      dimension cnow(*),p(64)
*i groupr.5668
            if (na.gt.64)  call error('f6ddx',
     &                   'order of legendre expansion > 64  ','  ')
*/*d groupr.5768
*/      dimension cnow(*),p(64)
*i groupr.5811
            if (na.gt.64)  call error('f6dis',
     &                   'order of legendre expansion > 64  ','  ')
*i groupr.8641
      parameter (npmax=15001,npmax0=npmax-1,mxset=200,mxnrf1=50)
      common/mainio/nsysi,nsyso,nsyse,ntty
      common/nalf5a/npf1(mxset),nbtf1(mxnrf1),intf1(mxnrf1)
      common/nalf5b/eintf1(mxset),xdatwk(npmax),xdatf1(npmax,mxset),
     1              ansf1(npmax),ydatwk(npmax),ydatf1(npmax,mxset),
     2              epmax(mxset)
*i groupr.8647
      external f5xlin,f5xint
*i groupr.8750
               nr=nint(c(l+4))
               intchk=0
               do k=1,nr
                  intnow=nint(c(nnt+ic-1+2*k))
                  if (intnow.le.0.or.(intnow.gt.6.and.intnow.le.20).or.
     &                intnow.gt.25) then
                     write(strng,'(a,i6)') 'int=',intnow
                     call error('getsed','interplation error',strng)
                  endif
                  if (intnow.gt.intchk) intchk=intnow
               enddo
*i groupr.8767
               if (intchk.ge.21.and.intchk.le.25) then
                  if (nne.gt.mxset)
     &               call error('getsed','storage exceeded - nne',' ')
                  eintf1(nne)=c(iraw+1)
                  npnow=nint(c(iraw+5))
                  nrf1=nint(c(iraw+4))
                  if (npnow.gt.npmax0)
     &               call error('getsed','storage exceeded - npnow',' ')
                  if (nrf1.gt.mxnrf1)
     &               call error('getsed','storage exceeded - nrf1',' ')
                  intone=0
                  do j=1,nrf1
                     nbtf1(j)=nint(c(iraw+4+2*j))
                     intf1(j)=nint(c(iraw+5+2*j))
                     if (intf1(j).eq.1) intone=1
                  enddo
                  do i=1,npmax
                     xdatwk(i)=0
                     ydatwk(i)=0
                  enddo
                  ipos=iraw+6+2*nrf1
                  do j=1,npnow
                     isw1=ipos+2*(j-1)
                     xdatwk(j)=c(isw1)
                     ydatwk(j)=c(isw1+1)
                  enddo
                  emaxf1=xdatwk(npnow)
                  epmax(nne)=emaxf1
                  do j=1,npnow
                     xdatwk(j)=xdatwk(j)/emaxf1
                     ydatwk(j)=ydatwk(j)*emaxf1
                  enddo
                  if (xdatwk(1).gt.zero.and.intf1(1).eq.1) then
                     if (ydatwk(1).eq.zero) xdatwk(1)=0
                  endif
c
                  npnew=0
                  call f5xlin(nbtf1,intf1,nrf1,xdatwk,ydatwk,npnow,
     &                        xdatf1(1,nne),ydatf1(1,nne),npmax0,npnew,
     &                        nsyso)
                  sumf1=0
                  npf1(nne)=npnew
                  if (intone.eq.1) then
                     do j=2,npf1(nne)-2
                        if (xdatf1(j,nne).eq.xdatf1(j+1,nne)) then
                           s1=xdatf1(j,nne)*0.999995d+0
                           s2=xdatf1(j,nne)*1.000005d+0
                           if (xdatf1(j-1,nne).lt.s1) xdatf1(j,nne)=s1
                           if (xdatf1(j+2,nne).gt.s2) xdatf1(j+1,nne)=s2
                        endif
                     enddo
                  endif
                  do j=2,npf1(nne)
                     sumf1=sumf1+(xdatf1(j,nne)-xdatf1(j-1,nne))*
     &                     (ydatf1(j,nne)+ydatf1(j-1,nne))/2
                  enddo
                  xdatf1(npnew+1,nne)=0
                  ydatf1(npnew+1,nne)=0
                  anorm=1
                  if (sumf1.gt.zero) anorm=1/sumf1
                  do j=1,npf1(nne)
                     ydatf1(j,nne)=anorm*ydatf1(j,nne)
                  enddo
               endif
*i groupr.8782
               if (intchk.ge.21.and.intchk.le.25) then
                  econst=0
               else
*i groupr.8813
               endif
*i groupr.8903
                  if (int.ge.21.and.int.le.25) then
                     call f5xint(ed,sed,eg,ng,nk,matd,mfd,ne,int,nne,
     &                           ikt,elo,ehi,pe,nsyso)
                  else
*i groupr.8907
                  endif
*i groupr.8966
      subroutine f5xlin(nbt,int,n1,xin,yin,n2,xout,yout,npmax,npout,
     1                  nsyso)
c     ******************************************************************
c     convert data to linear-linear form.
c     originality of this routine is subroutine filex in the linear-2000
c     code.
c     ******************************************************************
*if sw
      implicit real*8 (a-h,o-z)
*endif
      external error,terp1,sigfig
      save xcmin,errx,half
      dimension xin(n2),yin(n2),xout(npmax),yout(npmax)
      dimension nbt(n1),int(n1)
*if sw
      data xcmin/1.0d-15/
      data errx/5.0d-4/
      data half/0.5d0/
*else
      data xcmin/1.0e-15/
      data errx/5.0e-4/
      data half/0.5e0/
*endif
c
c-----initialize count of points in core and on scratch.
      n2core=1
      zero=0
c-----define first point in table.
      xout(1)=xin(1)
      yout(1)=yin(1)
c-----set flag to indicate if this is a threshold reaction.
      ithres=0
      if (yout(1).le.xcmin) ithres=1
      n2p1=2
      n2p2=3
c-----set up loop over interpolation ranges.
      nr2=1
      do 400 ir=1,n1
c-----define points in interpolation range and type of interpolation.
      nr1=nr2+1
      nr2=nbt(ir)
      intype=int(ir)
c-----data check
      minus=0
      do i=nr1,nr2
         if (yin(i).le.zero) minus=1
      enddo
c
      if (intype.gt.3.and.minus.eq.1) intype=2
c
c     non-linear interpolation law requires sub-division. set up loop
c     over points in current interpolation region.
c
c     interpolation law branch.
      if (intype-2) 50,80,100
c
c     for histogram replace each energy point that is internal to the
c     interpolation energy range by two energy points and replace each
c     end point of the interpolation energy range by a single energy
c     point (start of interval has already been stored).
c
   50 continue
      do 70 npt=nr1,nr2
      n2core=n2core+2
      if (n2core.gt.npmax) call error('f5xlin','storage exceeded',' ')
      xout(n2core)=xin(npt)
      yout(n2core)=yin(npt)
c-----eliminate duplicate points.
      if (xout(n2core).ne.xout(n2core-2)) go to 60
      xout(n2core-1)=xout(n2core)
      yout(n2core-1)=yout(n2core)
      n2core=n2core-1
      n2p1=n2core+1
      n2p2=n2p1+1
      go to 70
   60 continue
      xout(n2core-1)=xout(n2core)
      yout(n2core-1)=yout(n2core-2)
   70 continue
      n2core=n2core-1
      if (ir.lt.n1) then
         n2core=n2core+1
         if (n2core.gt.npmax)
     &      call error('f5xlin','storage exceeded - n2core',' ')
         xout(n2core)=xin(nr2)
         yout(n2core)=yin(nr2)
      endif
      n2p1=n2core+1
      n2p2=n2p1+1
      go to 400
c
c     for linear-linear interpolation just copy points to output array.
c
   80 continue
      do 90 npt=nr1,nr2
c-----core space for one more point required. if no room in core, thin
      n2core=n2core+1
      if (n2core.gt.npmax)
     &   call error('f5xlin','storage exceeded - n2core',' ')
      xout(n2core)=xin(npt)
      yout(n2core)=yin(npt)
   90 continue
      n2p1=n2core+1
      n2p2=n2p1+1
      go to 400
c
c     non-linear interpolation law requires sub-division. set up loop
c     over points in current interpolation region.
c
  100 continue
      do 300 npt=nr1,nr2
c-----initialize to range not saved
      imsaved=0
c-----define upper energy limit of interval.
      xend=xin(npt)
      yend=yin(npt)
c     if necessary insert max. energy or thermal point
  130 continue
      xtop=xend
      if (xout(n2core).le.zero) go to 150
      if (xtop.gt.2*xout(n2core)) xtop=2*xout(n2core)
      if (xtop.ge.xend) go to 150
c-----interpolate cross section to new end
      call terp1(xout(n2core),yout(n2core),xend,yend,xtop,ytop,intype)
c-----save end and set range saved flag
      imsaved=1
      xsave=xend
      ysave=yend
c-----define end of interval at thermal
      xend=xtop
      yend=ytop
c-----define energy and cross section at two ends of interval.
  150 continue
      xn2p2=xend
      yn2p2=yend
      xn2=xout(n2core)
      yn2=yout(n2core)
c-----do not subdivide if cross section is constant.
      if (yn2.eq.yn2p2) go to 220
c
c     do not sub-divide intervals within which the absolute value of
c     the cross section is less than the minimum cross section of
c     interest, accept near thresholds.
c
c-----if this is a threshold reaction do not use low cross section
c-----cutoff criteria in the vacinity of the threshold.
      if (ithres.eq.0) go to 160
      if (abs(yn2).lt.xcmin) go to 170
      ithres=0
      go to 170
c-----assume convergence if the absolute value of the cross section
c-----at both ends of the interval is less than the minimum cross
c-----section of interest.
  160 continue
      if (abs(yn2).le.xcmin.and.abs(yn2p2).le.xcmin) go to 220
c
c     define energy at middle of interval.
c
c-----intype is now 3, 4, 5 or 6.
  170 continue
      if (intype-4) 190,180,190
c-----linear energy.
  180 continue
      xn2p1=half*(xn2+xn2p2)
      go to 200
  190 continue
      xn2p1=sqrt(xn2*xn2p2)
c-----round midpoint.
  200 continue
      xn2p1=sigfig(xn2p1,7,0)
c
c     small energy interval convergence tests.
c
c-----if energy at middle of interval is less than allowable energy
c-----spacing only keep the two ends of the interval.
      if (xn2p1.le.xn2.or.xn2p1.ge.xn2p2) go to 220
c
c     define cross section at middle of interval by interpolation law
c     and linear-linear interpolation.
c
      call terp1(xn2,yn2,xn2p2,yn2p2,xn2p1,yn2p1,intype)
      call terp1(xn2,yn2,xn2p2,yn2p2,xn2p1,yaprox,2)
c
c-----test for convergence. if convergence keep midpoint of interval
c-----in order to allow accurate backward thinning before output.
      if (abs(yn2p1-yaprox).le.abs(yn2p1*errx)) go to 210
c
c     no convergence. shorten interval.
c
      xn2p2=xn2p1
      yn2p2=yn2p1
      go to 170
c
c     convergence. keep interval midpoint if backward thinning will be
c     performed. otherwise only keep ends of interval.
c
  210 continue
      if (n2p2.gt.npmax)
     &   call error('f5xlin','storage exceeded - n2p2',' ')
c-----keep interval midpoint and endpoint.
      xout(n2p1)=xn2p1
      yout(n2p1)=yn2p1
      xout(n2p2)=xn2p2
      yout(n2p2)=yn2p2
      n2core=n2p2
      n2p1=n2core+1
      n2p2=n2p1+1
      go to 230
c-----core space for one more point required. if no room in core, thin
c-----and unload one page of points to scratch.
  220 continue
      if(n2p1.gt.npmax)
     &  call error('f5xlin','storage exceeded - n2p1',' ')
c-----midpoint of interval is not required. only keep endpoint of
c-----interval.
      xout(n2p1)=xn2p2
      yout(n2p1)=yn2p2
      n2core=n2p1
      n2p1=n2core+1
      n2p2=n2p1+1
c-----if not end of current interval continue. otherwise
c-----move to next interval.
  230 continue
      if (xout(n2core).lt.xend) go to 150
c-----if range saved, restore point and reset flag
      if (imsaved.le.0) go to 300
      imsaved=0
      xend=xsave
      yend=ysave
      go to 130
  300 continue
c
  400 continue
      npout=n2core
      return
      end
c
      subroutine f5xint(ed,sed,eg,ng,nk,matd,mfd,ne,int,nne,ikt,elo,ehi,
     &                  pe,nsyso )
c     ******************************************************************
c     compute secondary energy distribution for all sink groups
c     simultaneously.  laws 1 with int=21-25 is coded.
c     ******************************************************************
*if sw
      implicit real*8 (a-h,o-z)
*endif
      parameter (npmax=15001,npmax0=npmax-1,mxset=200,mxnrf1=50)
      parameter (mxwork=10+2*npmax)
      common/nalf5a/npf1(mxset),nbtf1(mxnrf1),intf1(mxnrf1)
      common/nalf5b/eintf1(mxset),xdatwk(npmax),xdatf1(npmax,mxset),
     1              ansf1(npmax),ydatwk(npmax),ydatf1(npmax,mxset),
     2              epmax(mxset)
      dimension sed(nk,*),eg(*),c(mxwork)
      character strng*60
      external error,terp1,intega
      save ist,jst1,jst2
*if sw
      data small/1.d-10/
      data ebig/1.d8/
      data zero/0.d0/
*else
      data small/1.e-10/
      data ebig/1.e8/
      data zero/0.e0/
*endif
c
      if (nne.ge.ne) nne=ne-1
      nne1=nne+1
      npnow=npf1(nne)
      ist=2
      if (xdatf1(1,nne).ne.zero) then
         write(strng,'(a,i6)') 'xdatf1(1,nne).ne.0  nne=',nne
         call error('f5xint','invalid evaluation data',strng)
      elseif (xdatf1(1,nne1).ne.zero) then
         write(strng,'(a,i6)') 'xdatf1(1,nne1).ne.0  nne1=',nne1
         call error('f5xint','invalid evaluation data',strng)
      endif
c
c     make union relative secondary nuetron energy mesh
c     at requested energy range ( eintf1(nne) & eintf1(nne1) )
c
      do i=1,npmax
         xdatwk(i)=0
         ydatwk(i)=0
      enddo
      do i=1,npnow
         xdatwk(i)=xdatf1(i,nne)
      enddo
c
      do 170 j=2,npf1(nne1)-1
      enow=xdatf1(j,nne1)
      iadd=0
      do i=ist,npnow
         ehinow=xdatwk(i)
         elonow=xdatwk(i-1)
         ratio=1
         if (ehinow.gt.zero) ratio=enow/ehinow-1
         if (abs(ratio).lt.1.000d-5) go to 150
         if (enow.gt.elonow.and.enow.lt.ehinow) then
            iadd=i
            go to 150
         endif
      enddo
  150 continue
      if (iadd.ge.ist) then
         ratio=xdatwk(iadd)/enow
         if (npnow+1.gt.npmax)
     &      call error('f5xint','storage exceeded - npnow+1',' ')
         do k=iadd,npnow
            i=npnow-k+iadd
            xdatwk(i+1)=xdatwk(i)
         enddo
         xdatwk(iadd)=enow
         npnow=npnow+1
         ist=iadd
      endif
  170 continue
c
c     interpolate secondary neutron energy distribution at ed energy
c
      itype=int-20
      np1=npf1(nne)-1
      np2=npf1(nne1)-1
      jst1=1
      jst2=1
      eint1=eintf1(nne)
      eint2=eintf1(nne1)
c
c  *** set maximum energy of secondary neutron
c
      epmax1=epmax(nne)
      epmax2=epmax(nne1)
      epmaxn=0
      call terp1(eint1,epmax1,eint2,epmax2,ed,epmaxn,2)
      sumf1=0
      eold=0
      yold=0
c
      do 200 i=1,npnow
      enow=xdatwk(i)
      ynow1=0
      ynow2=0
      do j=jst1,np1
         jpos=j
         if (enow.ge.xdatf1(j,nne).and.enow.lt.xdatf1(j+1,nne))
     &      go to 180
      enddo
      jpos=np1
  180 continue
      jpos1=jpos+1
      jst1=jpos
      x1=xdatf1(jpos,nne)
      y1=ydatf1(jpos,nne)
      x2=xdatf1(jpos1,nne)
      y2=ydatf1(jpos1,nne)
      call terp1(x1,y1,x2,y2,enow,ynow1,2)
      do j=jst2,np2
         jpos=j
         if (enow.ge.xdatf1(j,nne1).and.enow.lt.xdatf1(j+1,nne1))
     &      go to 190
      enddo
      jpos=np2
  190 continue
      jpos1=jpos+1
      jst2=jpos
      x1=xdatf1(jpos,nne1)
      y1=ydatf1(jpos,nne1)
      x2=xdatf1(jpos1,nne1)
      y2=ydatf1(jpos1,nne1)
      call terp1(x1,y1,x2,y2,enow,ynow2,2)
      ynow=0
      call terp1(eint1,ynow1,eint2,ynow2,ed,ynow,itype)
      ydatwk(i)=ynow
      if (i.gt.1) sumf1=sumf1+(ynow+yold)*(enow-eold)/2
      eold=enow
      yold=ynow
  200 continue
c
c     define final secondary nuetron energy distribution data
c
      anorm=1
      if (sumf1.gt.zero) anorm=1/sumf1
      do i=1,npnow
         xdatwk(i)=xdatwk(i)*epmaxn
         ydatwk(i)=anorm*ydatwk(i)/epmaxn
      enddo
c
c     calculates groupwise secondary nuetron energy distribution data
c
      do ig=1,ng
         ansf1(ig)=0
      enddo
      if (xdatwk(npnow).le.eg(2)) then
         ansf1(1)=1
         factn=1
         go to 320
      endif
      c(1)=0
      c(2)=0
      c(3)=0
      c(4)=0
      c(5)=1
      c(6)=npnow
      c(7)=npnow
      c(8)=2
      isw=8
      do i=1,npnow
         isw=isw+1
         c(isw)=xdatwk(i)
         isw=isw+1
         c(isw)=ydatwk(i)
      enddo
      ip=2
      ir=1
      do ig=1,ng
         e1=eg(ig)
         if (ig.eq.1) e1=0
         e2=eg(ig+1)
         if (ig.eq.ng) e2=ebig
         call intega(ansf1(ig),e1,e2,c(1),ip,ir)
      enddo
      sump=zero
      do ig=1,ng
         sump=sump+ansf1(ig)
      enddo
      factn=1
      if (sump.gt.zero) factn=1/sump
c
c     store final result to sed array
c
  320 continue
      do ig=1,ng
         sed(ikt,ig)=sed(ikt,ig)+ansf1(ig)*pe*factn
      enddo
c
      return
      end
*/
*/  increase the allowed number of legendre terms in h6ddx
*/  to handle the new ENDF/B-VII mo95 evaluation.
*/
*d up15.9
      dimension cnow(*),p(65)
*d up30.10
*i heatr.3298
      data nlmax/65/
*/
*/ Kazuaki Kosako, Shimizu Corporation
*/ kinematic method
*/
*i heatr.55
c     *    kkerma   0/1=total (mt301) is energy balance / kinematic    *
c     *             (default=0)                                        *
*i heatr.91
      common/heat6/kkerma
*d heatr.139
      kkerma=0
      read(nsysi,*) matd,npk,nqa,ntemp,local,iprint,break,kkerma
*i heatr.151
      if (kkerma.eq.1) then
         if (npk.gt.0) then
            mtkk=0
            do i=1,npk
               if (mtk(i).eq.443) mtkk=i
            enddo
            if (mtkk.eq.0) then
               npk=npk+1
               if (kchk.eq.1) then
                  npkk=3*npk+7
               else
                  npkk=npk+3
               endif
               if (npkk.gt.npkmax) call error('heatr',
     &           'requested too many kerma mt-s (6+mt301 allowed).',' ')
               mtk(npk)=443
            endif
         else
            npk=1
            mtk(npk)=443
         endif
      endif
*i heatr.4997
      common/heat6/kkerma
*i heatr.5045
      if (kkerma.eq.1) then
         do i=2,npk
            if (mtp(i).eq.443) mtkk=i
         enddo
      endif
*i heatr.5121
      if (mt.eq.301.and.kkerma.eq.1) a(ibase+i)=c(mtkk)
*i heatr.5146
      if (kkerma.eq.1) write(nsyso,'(/'' total kerma (mt=301) was '',
     &  ''replaced to kinematic kerma (mt=443).'')')
*/
*ident upcad
*/ groupr -- 11oct2002
*/   add the ecco33, ecco1968, tripoli315, xmas172 and Vitamin-J
*/   group structures with 7 significant decimal digits. those
*/   group structures are used in Europe for fast breeder and
*/   thermal reactor neutronics calculations. for compatibility
*/   with Calendf and Apollo
*b groupr.132
c     *     18           xmas NEA-LANL
c     *     All new additional group structure with 7 significant
c     *     decimal digits compatible with CALENDF
c     *     19           ecco  33-group structure
c     *     20           ecco 1968-group structure
c     *     21           tripoli 315-group structure
c     *     22           xmas LWPC 172-group structure
c     *     23           vit-J LWPC 175-group structure
*b groupr.1297
c        19      ecco  33-group structure
c        20      ecco 1968-group structure
c        21      tripoli 315-group structure
c        22      xmas LWPC 172-group structure
c        23      vit-J LWPC 175-group structure
*b groupr.1312
      dimension eg19(34)
      dimension eg20(1969)
      dimension eg20a(95),eg20b(95),eg20c(95),eg20d(95),
     *          eg20e(95),eg20f(95),eg20g(95),eg20h(95),
     *          eg20i(95),eg20j(95),eg20k(95),eg20l(95),
     *          eg20m(95),eg20n(95),eg20o(95),eg20p(95),
     *          eg20q(95),eg20r(95),eg20s(95),eg20t(95),
     *          eg20u(69)
      dimension eg21(316)
      dimension eg21a(95),eg21b(95),eg21c(95),eg21d(31)
      dimension eg22(173)
      dimension eg22a(95),eg22b(78)
      dimension eg23(176)
      dimension eg23a(95),eg23b(81)
      equivalence (eg20a(1),eg20(1)),(eg20b(1),eg20(96)),
     *            (eg20c(1),eg20(191)),(eg20d(1),eg20(286)),
     *            (eg20e(1),eg20(381)),(eg20f(1),eg20(476)),
     *            (eg20g(1),eg20(571)),(eg20h(1),eg20(666)),
     *            (eg20i(1),eg20(761)),(eg20j(1),eg20(856)),
     *            (eg20k(1),eg20(951)),(eg20l(1),eg20(1046)),
     *            (eg20m(1),eg20(1141)),(eg20n(1),eg20(1236)),
     *            (eg20o(1),eg20(1331)),(eg20p(1),eg20(1426)),
     *            (eg20q(1),eg20(1521)),(eg20r(1),eg20(1616)),
     *            (eg20s(1),eg20(1711)),(eg20t(1),eg20(1806)),
     *            (eg20u(1),eg20(1901))
      equivalence (eg21a(1),eg21(1)),(eg21b(1),eg21(96)),
     *            (eg21c(1),eg21(191)),(eg21d(1),eg21(286))
      equivalence (eg22a(1),eg22(1)),(eg22b(1),eg22(96))
      equivalence (eg23a(1),eg23(1)),(eg23b(1),eg23(96))
*b groupr.1470
      data eg19/
     &1.000010d-05,1.000000d-01,5.400000d-01,4.000000d+00,8.315287d+00,
     &1.370959d+01,2.260329d+01,4.016900d+01,6.790405d+01,9.166088d+01,
     &1.486254d+02,3.043248d+02,4.539993d+02,7.485183d+02,1.234098d+03,
     &2.034684d+03,3.354626d+03,5.530844d+03,9.118820d+03,1.503439d+04,
     &2.478752d+04,4.086771d+04,6.737947d+04,1.110900d+05,1.831564d+05,
     &3.019738d+05,4.978707d+05,8.208500d+05,1.353353d+06,2.231302d+06,
     &3.678794d+06,6.065307d+06,1.000000d+07,1.964033d+07/
      data eg20a/
     &1.000010d-05,3.000000d-03,5.000000d-03,6.900000d-03,1.000000d-02,
     &1.500000d-02,2.000000d-02,2.500000d-02,3.000000d-02,3.500000d-02,
     &4.200000d-02,5.000000d-02,5.800000d-02,6.700000d-02,7.700000d-02,
     &8.000000d-02,9.500000d-02,1.000000d-01,1.150000d-01,1.340000d-01,
     &1.400000d-01,1.463700d-01,1.530300d-01,1.600000d-01,1.697100d-01,
     &1.800000d-01,1.890000d-01,1.988100d-01,2.091400d-01,2.200000d-01,
     &2.335800d-01,2.480000d-01,2.635100d-01,2.800000d-01,3.000000d-01,
     &3.145000d-01,3.200000d-01,3.346600d-01,3.500000d-01,3.699300d-01,
     &3.910000d-01,4.000000d-01,4.139900d-01,4.330000d-01,4.496800d-01,
     &4.670100d-01,4.850000d-01,5.000000d-01,5.196200d-01,5.315800d-01,
     &5.400000d-01,5.669600d-01,5.952800d-01,6.250000d-01,6.531500d-01,
     &6.825600d-01,7.050000d-01,7.415500d-01,7.800000d-01,7.900000d-01,
     &8.194500d-01,8.500000d-01,8.600000d-01,8.764250d-01,9.100000d-01,
     &9.300000d-01,9.500000d-01,9.720000d-01,9.860000d-01,9.960000d-01,
     &1.020000d+00,1.035000d+00,1.045000d+00,1.071000d+00,1.080000d+00,
     &1.097000d+00,1.110000d+00,1.123000d+00,1.150000d+00,1.170000d+00,
     &1.202060d+00,1.235000d+00,1.267080d+00,1.300000d+00,1.337500d+00,
     &1.370000d+00,1.404560d+00,1.440000d+00,1.475000d+00,1.500000d+00,
     &1.544340d+00,1.590000d+00,1.629510d+00,1.670000d+00,1.711970d+00/
      data eg20b/
     &1.755000d+00,1.797000d+00,1.840000d+00,1.855390d+00,1.884460d+00,
     &1.930000d+00,1.974490d+00,2.020000d+00,2.059610d+00,2.100000d+00,
     &2.130000d+00,2.185310d+00,2.242050d+00,2.300270d+00,2.360000d+00,
     &2.382370d+00,2.421710d+00,2.485030d+00,2.550000d+00,2.600000d+00,
     &2.659320d+00,2.720000d+00,2.767920d+00,2.837990d+00,2.909830d+00,
     &2.983490d+00,3.059020d+00,3.137330d+00,3.217630d+00,3.300000d+00,
     &3.380750d+00,3.466330d+00,3.554080d+00,3.644050d+00,3.736300d+00,
     &3.830880d+00,3.927860d+00,4.000000d+00,4.129250d+00,4.233782d+00,
     &4.340961d+00,4.450853d+00,4.563526d+00,4.679053d+00,4.797503d+00,
     &4.918953d+00,5.043477d+00,5.085681d+00,5.128239d+00,5.171153d+00,
     &5.214426d+00,5.258061d+00,5.302061d+00,5.346430d+00,5.391169d+00,
     &5.436284d+00,5.481775d+00,5.527647d+00,5.573904d+00,5.620547d+00,
     &5.667581d+00,5.715008d+00,5.762832d+00,5.811056d+00,5.859684d+00,
     &5.908719d+00,5.958164d+00,6.008022d+00,6.058298d+00,6.108995d+00,
     &6.160116d+00,6.211665d+00,6.263645d+00,6.316060d+00,6.368914d+00,
     &6.422210d+00,6.475952d+00,6.530144d+00,6.584789d+00,6.639892d+00,
     &6.695455d+00,6.751484d+00,6.807981d+00,6.864952d+00,6.922399d+00,
     &6.980326d+00,7.038739d+00,7.097640d+00,7.157034d+00,7.216925d+00,
     &7.277317d+00,7.338215d+00,7.399622d+00,7.461544d+00,7.523983d+00/
      data eg20c/
     &7.586945d+00,7.650434d+00,7.714454d+00,7.779009d+00,7.844105d+00,
     &7.909746d+00,7.975936d+00,8.042680d+00,8.109982d+00,8.177848d+00,
     &8.246281d+00,8.315287d+00,8.384871d+00,8.455037d+00,8.525790d+00,
     &8.597135d+00,8.669077d+00,8.741621d+00,8.814772d+00,8.888536d+00,
     &8.962916d+00,9.037919d+00,9.113550d+00,9.189814d+00,9.266715d+00,
     &9.344261d+00,9.422455d+00,9.501303d+00,9.580812d+00,9.660985d+00,
     &9.741830d+00,9.823351d+00,9.905554d+00,9.988446d+00,1.007203d+01,
     &1.015631d+01,1.024130d+01,1.032701d+01,1.041342d+01,1.050056d+01,
     &1.058843d+01,1.067704d+01,1.076639d+01,1.085648d+01,1.094733d+01,
     &1.103894d+01,1.113132d+01,1.122446d+01,1.131839d+01,1.141311d+01,
     &1.150861d+01,1.160492d+01,1.170203d+01,1.179995d+01,1.189870d+01,
     &1.199827d+01,1.209867d+01,1.219991d+01,1.230201d+01,1.240495d+01,
     &1.250876d+01,1.261343d+01,1.271898d+01,1.282542d+01,1.293274d+01,
     &1.304097d+01,1.315010d+01,1.326014d+01,1.337110d+01,1.348299d+01,
     &1.359582d+01,1.370959d+01,1.382431d+01,1.394000d+01,1.405665d+01,
     &1.417428d+01,1.429289d+01,1.441250d+01,1.453310d+01,1.465472d+01,
     &1.477735d+01,1.490101d+01,1.502570d+01,1.515144d+01,1.527823d+01,
     &1.540608d+01,1.553500d+01,1.566500d+01,1.579609d+01,1.592827d+01,
     &1.606156d+01,1.619597d+01,1.633150d+01,1.646816d+01,1.660597d+01/
      data eg20d/
     &1.674493d+01,1.688506d+01,1.702635d+01,1.716883d+01,1.731250d+01,
     &1.745738d+01,1.760346d+01,1.775077d+01,1.789931d+01,1.804910d+01,
     &1.820013d+01,1.835244d+01,1.850601d+01,1.866087d+01,1.881703d+01,
     &1.897449d+01,1.913328d+01,1.929339d+01,1.945484d+01,1.961764d+01,
     &1.978180d+01,1.994734d+01,2.011426d+01,2.028258d+01,2.045231d+01,
     &2.062345d+01,2.079603d+01,2.097006d+01,2.114554d+01,2.132249d+01,
     &2.150092d+01,2.168084d+01,2.186227d+01,2.204522d+01,2.222969d+01,
     &2.241572d+01,2.260329d+01,2.279244d+01,2.298317d+01,2.317550d+01,
     &2.336944d+01,2.356499d+01,2.376219d+01,2.396104d+01,2.416154d+01,
     &2.436373d+01,2.456761d+01,2.477320d+01,2.498050d+01,2.518954d+01,
     &2.540033d+01,2.561289d+01,2.582722d+01,2.604335d+01,2.626128d+01,
     &2.648104d+01,2.670264d+01,2.692609d+01,2.715141d+01,2.737862d+01,
     &2.760773d+01,2.783875d+01,2.807171d+01,2.830662d+01,2.854349d+01,
     &2.878235d+01,2.902320d+01,2.926607d+01,2.951098d+01,2.975793d+01,
     &3.000695d+01,3.025805d+01,3.051126d+01,3.076658d+01,3.102404d+01,
     &3.128365d+01,3.154544d+01,3.180942d+01,3.207560d+01,3.234401d+01,
     &3.261467d+01,3.288760d+01,3.316281d+01,3.344032d+01,3.372015d+01,
     &3.400233d+01,3.428686d+01,3.457378d+01,3.486310d+01,3.515484d+01,
     &3.544902d+01,3.574566d+01,3.604479d+01,3.634642d+01,3.665057d+01/
      data eg20e/
     &3.695727d+01,3.726653d+01,3.757838d+01,3.789285d+01,3.820994d+01,
     &3.852969d+01,3.885211d+01,3.917723d+01,3.950507d+01,3.983565d+01,
     &4.016900d+01,4.050514d+01,4.084410d+01,4.118589d+01,4.153054d+01,
     &4.187807d+01,4.222851d+01,4.258189d+01,4.293822d+01,4.329753d+01,
     &4.365985d+01,4.402521d+01,4.439361d+01,4.476511d+01,4.513971d+01,
     &4.551744d+01,4.589834d+01,4.628243d+01,4.666972d+01,4.706026d+01,
     &4.745407d+01,4.785117d+01,4.825160d+01,4.865538d+01,4.906253d+01,
     &4.947309d+01,4.988709d+01,5.030456d+01,5.072551d+01,5.114999d+01,
     &5.157802d+01,5.200963d+01,5.244486d+01,5.288373d+01,5.332626d+01,
     &5.377251d+01,5.422248d+01,5.467623d+01,5.513376d+01,5.559513d+01,
     &5.606036d+01,5.652948d+01,5.700253d+01,5.747954d+01,5.796053d+01,
     &5.844556d+01,5.893464d+01,5.942781d+01,5.992511d+01,6.042657d+01,
     &6.093223d+01,6.144212d+01,6.195628d+01,6.247474d+01,6.299754d+01,
     &6.352471d+01,6.405630d+01,6.459233d+01,6.513285d+01,6.567789d+01,
     &6.622749d+01,6.678169d+01,6.734053d+01,6.790405d+01,6.847228d+01,
     &6.904527d+01,6.962305d+01,7.020566d+01,7.079316d+01,7.138556d+01,
     &7.198293d+01,7.258529d+01,7.319270d+01,7.380518d+01,7.442280d+01,
     &7.504558d+01,7.567357d+01,7.630682d+01,7.694537d+01,7.758926d+01,
     &7.823854d+01,7.889325d+01,7.955344d+01,8.021915d+01,8.089044d+01/
      data eg20f/
     &8.156734d+01,8.224991d+01,8.293819d+01,8.363223d+01,8.433208d+01,
     &8.503778d+01,8.574939d+01,8.646695d+01,8.719052d+01,8.792015d+01,
     &8.865588d+01,8.939776d+01,9.014586d+01,9.090021d+01,9.166088d+01,
     &9.242791d+01,9.320136d+01,9.398128d+01,9.476773d+01,9.556076d+01,
     &9.636043d+01,9.716679d+01,9.797990d+01,9.879981d+01,9.962658d+01,
     &1.004603d+02,1.013009d+02,1.021486d+02,1.030034d+02,1.038654d+02,
     &1.047345d+02,1.056110d+02,1.064947d+02,1.073859d+02,1.082845d+02,
     &1.091907d+02,1.101044d+02,1.110258d+02,1.119548d+02,1.128917d+02,
     &1.138364d+02,1.147890d+02,1.157496d+02,1.167182d+02,1.176949d+02,
     &1.186798d+02,1.196729d+02,1.206744d+02,1.216842d+02,1.227024d+02,
     &1.237292d+02,1.247646d+02,1.258087d+02,1.268615d+02,1.279231d+02,
     &1.289935d+02,1.300730d+02,1.311615d+02,1.322590d+02,1.333658d+02,
     &1.344818d+02,1.356072d+02,1.367420d+02,1.378862d+02,1.390401d+02,
     &1.402036d+02,1.413768d+02,1.425599d+02,1.437529d+02,1.449558d+02,
     &1.461688d+02,1.473920d+02,1.486254d+02,1.498691d+02,1.511232d+02,
     &1.523879d+02,1.536631d+02,1.549489d+02,1.562456d+02,1.575531d+02,
     &1.588715d+02,1.602010d+02,1.615415d+02,1.628933d+02,1.642565d+02,
     &1.656310d+02,1.670170d+02,1.684146d+02,1.698239d+02,1.712451d+02,
     &1.726781d+02,1.741231d+02,1.755802d+02,1.770494d+02,1.785310d+02/
      data eg20g/
     &1.800250d+02,1.815315d+02,1.830505d+02,1.845823d+02,1.861269d+02,
     &1.876845d+02,1.892551d+02,1.908388d+02,1.924358d+02,1.940461d+02,
     &1.956699d+02,1.973073d+02,1.989584d+02,2.006233d+02,2.023021d+02,
     &2.039950d+02,2.057021d+02,2.074234d+02,2.091592d+02,2.109095d+02,
     &2.126744d+02,2.144541d+02,2.162487d+02,2.180583d+02,2.198830d+02,
     &2.217230d+02,2.235784d+02,2.254494d+02,2.273360d+02,2.292384d+02,
     &2.311567d+02,2.330910d+02,2.350416d+02,2.370084d+02,2.389917d+02,
     &2.409917d+02,2.430083d+02,2.450418d+02,2.470924d+02,2.491601d+02,
     &2.512451d+02,2.533476d+02,2.554676d+02,2.576054d+02,2.597611d+02,
     &2.619348d+02,2.641267d+02,2.663370d+02,2.685657d+02,2.708131d+02,
     &2.730793d+02,2.753645d+02,2.776688d+02,2.799924d+02,2.823354d+02,
     &2.846980d+02,2.870804d+02,2.894827d+02,2.919052d+02,2.943479d+02,
     &2.968110d+02,2.992948d+02,3.017993d+02,3.043248d+02,3.068715d+02,
     &3.094394d+02,3.120288d+02,3.146399d+02,3.172729d+02,3.199279d+02,
     &3.226051d+02,3.253047d+02,3.280269d+02,3.307719d+02,3.335398d+02,
     &3.363309d+02,3.391454d+02,3.419834d+02,3.448452d+02,3.477309d+02,
     &3.506408d+02,3.535750d+02,3.565338d+02,3.595173d+02,3.625258d+02,
     &3.655595d+02,3.686185d+02,3.717032d+02,3.748137d+02,3.779502d+02,
     &3.811129d+02,3.843021d+02,3.875180d+02,3.907608d+02,3.940308d+02/
      data eg20h/
     &3.973281d+02,4.006530d+02,4.040057d+02,4.073865d+02,4.107955d+02,
     &4.142332d+02,4.176995d+02,4.211949d+02,4.247195d+02,4.282736d+02,
     &4.318575d+02,4.354713d+02,4.391154d+02,4.427900d+02,4.464953d+02,
     &4.502317d+02,4.539993d+02,4.577984d+02,4.616294d+02,4.654923d+02,
     &4.693877d+02,4.733156d+02,4.772763d+02,4.812703d+02,4.852976d+02,
     &4.893587d+02,4.934537d+02,4.975830d+02,5.017468d+02,5.059455d+02,
     &5.101793d+02,5.144486d+02,5.187536d+02,5.230946d+02,5.274719d+02,
     &5.318859d+02,5.363368d+02,5.408249d+02,5.453506d+02,5.499142d+02,
     &5.545160d+02,5.591563d+02,5.638354d+02,5.685536d+02,5.733114d+02,
     &5.781089d+02,5.829466d+02,5.878248d+02,5.927438d+02,5.977040d+02,
     &6.027057d+02,6.077492d+02,6.128350d+02,6.179633d+02,6.231345d+02,
     &6.283489d+02,6.336071d+02,6.389092d+02,6.442557d+02,6.496469d+02,
     &6.550832d+02,6.605651d+02,6.660928d+02,6.716668d+02,6.772874d+02,
     &6.829550d+02,6.886701d+02,6.944330d+02,7.002441d+02,7.061038d+02,
     &7.120126d+02,7.179709d+02,7.239790d+02,7.300373d+02,7.361464d+02,
     &7.423066d+02,7.485183d+02,7.547820d+02,7.610981d+02,7.674671d+02,
     &7.738894d+02,7.803654d+02,7.868957d+02,7.934805d+02,8.001205d+02,
     &8.068160d+02,8.135676d+02,8.203756d+02,8.272407d+02,8.341631d+02,
     &8.411435d+02,8.481824d+02,8.552801d+02,8.624372d+02,8.696542d+02/
      data eg20i/
     &8.769316d+02,8.842699d+02,8.916696d+02,8.991312d+02,9.066553d+02,
     &9.142423d+02,9.218928d+02,9.296074d+02,9.373865d+02,9.452307d+02,
     &9.531405d+02,9.611165d+02,9.691593d+02,9.772694d+02,9.854473d+02,
     &9.936937d+02,1.002009d+03,1.010394d+03,1.018849d+03,1.027375d+03,
     &1.035972d+03,1.044641d+03,1.053383d+03,1.062198d+03,1.071087d+03,
     &1.080050d+03,1.089088d+03,1.098201d+03,1.107391d+03,1.116658d+03,
     &1.126002d+03,1.135425d+03,1.144926d+03,1.154507d+03,1.164168d+03,
     &1.173910d+03,1.183734d+03,1.193639d+03,1.203628d+03,1.213700d+03,
     &1.223857d+03,1.234098d+03,1.244425d+03,1.254839d+03,1.265339d+03,
     &1.275928d+03,1.286605d+03,1.297372d+03,1.308228d+03,1.319176d+03,
     &1.330215d+03,1.341346d+03,1.352571d+03,1.363889d+03,1.375303d+03,
     &1.386811d+03,1.398416d+03,1.410118d+03,1.421919d+03,1.433817d+03,
     &1.445816d+03,1.457915d+03,1.470115d+03,1.482417d+03,1.494822d+03,
     &1.507331d+03,1.519944d+03,1.532663d+03,1.545489d+03,1.558422d+03,
     &1.571463d+03,1.584613d+03,1.597874d+03,1.611245d+03,1.624728d+03,
     &1.638324d+03,1.652034d+03,1.665858d+03,1.679798d+03,1.693855d+03,
     &1.708030d+03,1.722323d+03,1.736735d+03,1.751268d+03,1.765923d+03,
     &1.780701d+03,1.795602d+03,1.810628d+03,1.825780d+03,1.841058d+03,
     &1.856464d+03,1.871999d+03,1.887665d+03,1.903461d+03,1.919389d+03/
      data eg20j/
     &1.935451d+03,1.951647d+03,1.967979d+03,1.984447d+03,2.001053d+03,
     &2.017798d+03,2.034684d+03,2.051710d+03,2.068879d+03,2.086192d+03,
     &2.103650d+03,2.121253d+03,2.139004d+03,2.156904d+03,2.174953d+03,
     &2.193153d+03,2.211506d+03,2.230012d+03,2.248673d+03,2.267490d+03,
     &2.286465d+03,2.305599d+03,2.324892d+03,2.344347d+03,2.363965d+03,
     &2.383747d+03,2.403695d+03,2.423809d+03,2.444092d+03,2.464545d+03,
     &2.485168d+03,2.505965d+03,2.526935d+03,2.548081d+03,2.569403d+03,
     &2.590904d+03,2.612586d+03,2.634448d+03,2.656494d+03,2.678723d+03,
     &2.701139d+03,2.723743d+03,2.746536d+03,2.769519d+03,2.792695d+03,
     &2.816065d+03,2.839630d+03,2.863392d+03,2.887354d+03,2.911515d+03,
     &2.935879d+03,2.960447d+03,2.985221d+03,3.010202d+03,3.035391d+03,
     &3.060792d+03,3.086405d+03,3.112233d+03,3.138276d+03,3.164538d+03,
     &3.191019d+03,3.217722d+03,3.244649d+03,3.271800d+03,3.299179d+03,
     &3.326787d+03,3.354626d+03,3.382698d+03,3.411005d+03,3.439549d+03,
     &3.468332d+03,3.497355d+03,3.526622d+03,3.556133d+03,3.585891d+03,
     &3.615898d+03,3.646157d+03,3.676668d+03,3.707435d+03,3.738460d+03,
     &3.769744d+03,3.801290d+03,3.833099d+03,3.865175d+03,3.897520d+03,
     &3.930135d+03,3.963023d+03,3.996186d+03,4.029627d+03,4.063347d+03,
     &4.097350d+03,4.131637d+03,4.166211d+03,4.201075d+03,4.236230d+03/
      data eg20k/
     &4.271679d+03,4.307425d+03,4.343471d+03,4.379817d+03,4.416468d+03,
     &4.453426d+03,4.490693d+03,4.528272d+03,4.566165d+03,4.604375d+03,
     &4.642906d+03,4.681758d+03,4.720936d+03,4.760441d+03,4.800277d+03,
     &4.840447d+03,4.880952d+03,4.921797d+03,4.962983d+03,5.004514d+03,
     &5.046393d+03,5.088622d+03,5.131204d+03,5.174143d+03,5.217441d+03,
     &5.261101d+03,5.305127d+03,5.349521d+03,5.394287d+03,5.439427d+03,
     &5.484945d+03,5.530844d+03,5.577127d+03,5.623797d+03,5.670858d+03,
     &5.718312d+03,5.766164d+03,5.814416d+03,5.863072d+03,5.912135d+03,
     &5.961609d+03,6.011496d+03,6.061802d+03,6.112528d+03,6.163678d+03,
     &6.215257d+03,6.267267d+03,6.319712d+03,6.372597d+03,6.425924d+03,
     &6.479697d+03,6.533920d+03,6.588597d+03,6.643731d+03,6.699327d+03,
     &6.755388d+03,6.811918d+03,6.868921d+03,6.926401d+03,6.984362d+03,
     &7.042809d+03,7.101744d+03,7.161172d+03,7.221098d+03,7.281525d+03,
     &7.342458d+03,7.403901d+03,7.465858d+03,7.528334d+03,7.591332d+03,
     &7.654857d+03,7.718914d+03,7.783507d+03,7.848641d+03,7.914319d+03,
     &7.980548d+03,8.047330d+03,8.114671d+03,8.182576d+03,8.251049d+03,
     &8.320095d+03,8.389719d+03,8.459926d+03,8.530719d+03,8.602106d+03,
     &8.674090d+03,8.746676d+03,8.819869d+03,8.893675d+03,8.968099d+03,
     &9.043145d+03,9.118820d+03,9.195127d+03,9.272074d+03,9.349664d+03/
      data eg20l/
     &9.427903d+03,9.506797d+03,9.586352d+03,9.666572d+03,9.747463d+03,
     &9.829031d+03,9.911282d+03,9.994221d+03,1.007785d+04,1.016219d+04,
     &1.024723d+04,1.033298d+04,1.041944d+04,1.050664d+04,1.059456d+04,
     &1.068321d+04,1.077261d+04,1.086276d+04,1.095366d+04,1.104532d+04,
     &1.113775d+04,1.123095d+04,1.132494d+04,1.141970d+04,1.151527d+04,
     &1.161163d+04,1.170880d+04,1.180678d+04,1.190558d+04,1.200521d+04,
     &1.210567d+04,1.220697d+04,1.230912d+04,1.241212d+04,1.251599d+04,
     &1.262073d+04,1.272634d+04,1.283283d+04,1.294022d+04,1.304851d+04,
     &1.315770d+04,1.326780d+04,1.337883d+04,1.349079d+04,1.360368d+04,
     &1.371752d+04,1.383231d+04,1.394806d+04,1.406478d+04,1.418247d+04,
     &1.430116d+04,1.442083d+04,1.454151d+04,1.466319d+04,1.478590d+04,
     &1.490963d+04,1.503439d+04,1.516020d+04,1.528706d+04,1.541499d+04,
     &1.554398d+04,1.567406d+04,1.580522d+04,1.593748d+04,1.607085d+04,
     &1.620533d+04,1.634094d+04,1.647768d+04,1.661557d+04,1.675461d+04,
     &1.689482d+04,1.703620d+04,1.717876d+04,1.732251d+04,1.746747d+04,
     &1.761364d+04,1.776104d+04,1.790966d+04,1.805953d+04,1.821066d+04,
     &1.836305d+04,1.851671d+04,1.867166d+04,1.882791d+04,1.898547d+04,
     &1.914434d+04,1.930454d+04,1.946608d+04,1.962898d+04,1.979324d+04,
     &1.995887d+04,2.012589d+04,2.029431d+04,2.046413d+04,2.063538d+04/
      data eg20m/
     &2.080806d+04,2.098218d+04,2.115777d+04,2.133482d+04,2.151335d+04,
     &2.169338d+04,2.187491d+04,2.205796d+04,2.224255d+04,2.242868d+04,
     &2.261636d+04,2.280562d+04,2.299646d+04,2.318890d+04,2.338295d+04,
     &2.357862d+04,2.377593d+04,2.397489d+04,2.417552d+04,2.437782d+04,
     &2.458182d+04,2.478752d+04,2.499495d+04,2.520411d+04,2.541502d+04,
     &2.562770d+04,2.584215d+04,2.605841d+04,2.627647d+04,2.649635d+04,
     &2.671808d+04,2.694166d+04,2.700000d+04,2.716711d+04,2.739445d+04,
     &2.762369d+04,2.785485d+04,2.808794d+04,2.832299d+04,2.850000d+04,
     &2.856000d+04,2.879899d+04,2.903999d+04,2.928300d+04,2.952804d+04,
     &2.977514d+04,3.002430d+04,3.027555d+04,3.052890d+04,3.078437d+04,
     &3.104198d+04,3.130174d+04,3.156368d+04,3.182781d+04,3.209415d+04,
     &3.236272d+04,3.263353d+04,3.290662d+04,3.318198d+04,3.345965d+04,
     &3.373965d+04,3.402199d+04,3.430669d+04,3.459377d+04,3.488326d+04,
     &3.517517d+04,3.546952d+04,3.576633d+04,3.606563d+04,3.636743d+04,
     &3.667176d+04,3.697864d+04,3.728808d+04,3.760011d+04,3.791476d+04,
     &3.823203d+04,3.855196d+04,3.887457d+04,3.919988d+04,3.952791d+04,
     &3.985869d+04,4.019223d+04,4.052857d+04,4.086771d+04,4.120970d+04,
     &4.155455d+04,4.190229d+04,4.225293d+04,4.260651d+04,4.296305d+04,
     &4.332257d+04,4.368510d+04,4.405066d+04,4.441928d+04,4.479099d+04/
      data eg20n/
     &4.516581d+04,4.554376d+04,4.592488d+04,4.630919d+04,4.669671d+04,
     &4.708747d+04,4.748151d+04,4.787884d+04,4.827950d+04,4.868351d+04,
     &4.909090d+04,4.950170d+04,4.991594d+04,5.033364d+04,5.075484d+04,
     &5.117957d+04,5.160785d+04,5.203971d+04,5.247518d+04,5.291430d+04,
     &5.335710d+04,5.380360d+04,5.425384d+04,5.470784d+04,5.516564d+04,
     &5.562728d+04,5.609278d+04,5.656217d+04,5.703549d+04,5.751277d+04,
     &5.799405d+04,5.847935d+04,5.896871d+04,5.946217d+04,5.995976d+04,
     &6.046151d+04,6.096747d+04,6.147765d+04,6.199211d+04,6.251086d+04,
     &6.303396d+04,6.356144d+04,6.409333d+04,6.462968d+04,6.517051d+04,
     &6.571586d+04,6.626579d+04,6.682031d+04,6.737947d+04,6.794331d+04,
     &6.851187d+04,6.908519d+04,6.966330d+04,7.024626d+04,7.083409d+04,
     &7.142684d+04,7.202455d+04,7.262726d+04,7.323502d+04,7.384786d+04,
     &7.446583d+04,7.508897d+04,7.571733d+04,7.635094d+04,7.698986d+04,
     &7.763412d+04,7.828378d+04,7.893887d+04,7.950000d+04,7.959944d+04,
     &8.026554d+04,8.093721d+04,8.161451d+04,8.229747d+04,8.250000d+04,
     &8.298615d+04,8.368059d+04,8.438084d+04,8.508695d+04,8.579897d+04,
     &8.651695d+04,8.724094d+04,8.797098d+04,8.870714d+04,8.944945d+04,
     &9.019798d+04,9.095277d+04,9.171388d+04,9.248135d+04,9.325525d+04,
     &9.403563d+04,9.482253d+04,9.561602d+04,9.641615d+04,9.722297d+04/
      data eg20o/
     &9.803655d+04,9.885694d+04,9.968419d+04,1.005184d+05,1.013595d+05,
     &1.022077d+05,1.030630d+05,1.039254d+05,1.047951d+05,1.056720d+05,
     &1.065563d+05,1.074480d+05,1.083471d+05,1.092538d+05,1.101681d+05,
     &1.110900d+05,1.120196d+05,1.129570d+05,1.139022d+05,1.148554d+05,
     &1.158165d+05,1.167857d+05,1.177629d+05,1.187484d+05,1.197421d+05,
     &1.207441d+05,1.217545d+05,1.227734d+05,1.238008d+05,1.248368d+05,
     &1.258814d+05,1.269348d+05,1.279970d+05,1.290681d+05,1.301482d+05,
     &1.312373d+05,1.323355d+05,1.334429d+05,1.345596d+05,1.356856d+05,
     &1.368210d+05,1.379660d+05,1.391205d+05,1.402847d+05,1.414586d+05,
     &1.426423d+05,1.438360d+05,1.450396d+05,1.462533d+05,1.474772d+05,
     &1.487113d+05,1.499558d+05,1.512106d+05,1.524760d+05,1.537519d+05,
     &1.550385d+05,1.563359d+05,1.576442d+05,1.589634d+05,1.602936d+05,
     &1.616349d+05,1.629875d+05,1.643514d+05,1.657268d+05,1.671136d+05,
     &1.685120d+05,1.699221d+05,1.713441d+05,1.727779d+05,1.742237d+05,
     &1.756817d+05,1.771518d+05,1.786342d+05,1.801291d+05,1.816364d+05,
     &1.831564d+05,1.846891d+05,1.862346d+05,1.877930d+05,1.893645d+05,
     &1.909491d+05,1.925470d+05,1.941583d+05,1.957830d+05,1.974214d+05,
     &1.990734d+05,2.007393d+05,2.024191d+05,2.041130d+05,2.058210d+05,
     &2.075434d+05,2.092801d+05,2.110314d+05,2.127974d+05,2.145781d+05/
      data eg20p/
     &2.163737d+05,2.181844d+05,2.200102d+05,2.218512d+05,2.237077d+05,
     &2.255797d+05,2.274674d+05,2.293709d+05,2.312903d+05,2.332258d+05,
     &2.351775d+05,2.371455d+05,2.391299d+05,2.411310d+05,2.431488d+05,
     &2.451835d+05,2.472353d+05,2.493042d+05,2.513904d+05,2.534941d+05,
     &2.556153d+05,2.577544d+05,2.599113d+05,2.620863d+05,2.642794d+05,
     &2.664910d+05,2.687210d+05,2.709697d+05,2.732372d+05,2.755237d+05,
     &2.778293d+05,2.801543d+05,2.824986d+05,2.848626d+05,2.872464d+05,
     &2.896501d+05,2.920740d+05,2.945181d+05,2.969826d+05,2.972000d+05,
     &2.985000d+05,2.994678d+05,3.019738d+05,3.045008d+05,3.070489d+05,
     &3.096183d+05,3.122093d+05,3.148219d+05,3.174564d+05,3.201129d+05,
     &3.227916d+05,3.254928d+05,3.282166d+05,3.309631d+05,3.337327d+05,
     &3.365254d+05,3.393415d+05,3.421812d+05,3.450446d+05,3.479320d+05,
     &3.508435d+05,3.537795d+05,3.567399d+05,3.597252d+05,3.627354d+05,
     &3.657708d+05,3.688317d+05,3.719181d+05,3.750304d+05,3.781687d+05,
     &3.813333d+05,3.845243d+05,3.877421d+05,3.909868d+05,3.942586d+05,
     &3.975578d+05,4.008846d+05,4.042393d+05,4.076220d+05,4.110331d+05,
     &4.144727d+05,4.179410d+05,4.214384d+05,4.249651d+05,4.285213d+05,
     &4.321072d+05,4.357231d+05,4.393693d+05,4.430460d+05,4.467535d+05,
     &4.504920d+05,4.542618d+05,4.580631d+05,4.618963d+05,4.657615d+05/
      data eg20q/
     &4.696591d+05,4.735892d+05,4.775523d+05,4.815485d+05,4.855782d+05,
     &4.896416d+05,4.937390d+05,4.978707d+05,5.020369d+05,5.062381d+05,
     &5.104743d+05,5.147461d+05,5.190535d+05,5.233971d+05,5.277769d+05,
     &5.321934d+05,5.366469d+05,5.411377d+05,5.456660d+05,5.502322d+05,
     &5.548366d+05,5.594796d+05,5.641614d+05,5.688824d+05,5.736429d+05,
     &5.784432d+05,5.832837d+05,5.881647d+05,5.930866d+05,5.980496d+05,
     &6.030542d+05,6.081006d+05,6.131893d+05,6.183206d+05,6.234948d+05,
     &6.287123d+05,6.339734d+05,6.392786d+05,6.446282d+05,6.500225d+05,
     &6.554620d+05,6.609470d+05,6.664779d+05,6.720551d+05,6.776790d+05,
     &6.833499d+05,6.890683d+05,6.948345d+05,7.006490d+05,7.065121d+05,
     &7.124243d+05,7.183860d+05,7.243976d+05,7.304594d+05,7.365720d+05,
     &7.427358d+05,7.489511d+05,7.552184d+05,7.615382d+05,7.679109d+05,
     &7.743369d+05,7.808167d+05,7.873507d+05,7.939393d+05,8.005831d+05,
     &8.072825d+05,8.140380d+05,8.208500d+05,8.277190d+05,8.346455d+05,
     &8.416299d+05,8.486728d+05,8.557746d+05,8.629359d+05,8.701570d+05,
     &8.774387d+05,8.847812d+05,8.921852d+05,8.996511d+05,9.071795d+05,
     &9.147709d+05,9.224259d+05,9.301449d+05,9.379285d+05,9.457772d+05,
     &9.536916d+05,9.616723d+05,9.697197d+05,9.778344d+05,9.860171d+05,
     &9.942682d+05,1.002588d+06,1.010978d+06,1.019438d+06,1.027969d+06/
      data eg20r/
     &1.036571d+06,1.045245d+06,1.053992d+06,1.062812d+06,1.071706d+06,
     &1.080674d+06,1.089717d+06,1.098836d+06,1.108032d+06,1.117304d+06,
     &1.126654d+06,1.136082d+06,1.145588d+06,1.155175d+06,1.164842d+06,
     &1.174589d+06,1.184418d+06,1.194330d+06,1.204324d+06,1.214402d+06,
     &1.224564d+06,1.234812d+06,1.245145d+06,1.255564d+06,1.266071d+06,
     &1.276666d+06,1.287349d+06,1.298122d+06,1.308985d+06,1.319938d+06,
     &1.330984d+06,1.342122d+06,1.353353d+06,1.364678d+06,1.376098d+06,
     &1.387613d+06,1.399225d+06,1.410934d+06,1.422741d+06,1.434646d+06,
     &1.446652d+06,1.458758d+06,1.470965d+06,1.483274d+06,1.495686d+06,
     &1.508202d+06,1.520823d+06,1.533550d+06,1.546383d+06,1.559323d+06,
     &1.572372d+06,1.585530d+06,1.598797d+06,1.612176d+06,1.625667d+06,
     &1.639271d+06,1.652989d+06,1.666821d+06,1.680770d+06,1.694834d+06,
     &1.709017d+06,1.723318d+06,1.737739d+06,1.752281d+06,1.766944d+06,
     &1.781731d+06,1.796640d+06,1.811675d+06,1.826835d+06,1.842122d+06,
     &1.857538d+06,1.873082d+06,1.888756d+06,1.904561d+06,1.920499d+06,
     &1.936570d+06,1.952776d+06,1.969117d+06,1.985595d+06,2.002210d+06,
     &2.018965d+06,2.035860d+06,2.052897d+06,2.070076d+06,2.087398d+06,
     &2.104866d+06,2.122480d+06,2.140241d+06,2.158151d+06,2.176211d+06,
     &2.194421d+06,2.212785d+06,2.231302d+06,2.249973d+06,2.268802d+06/
      data eg20s/
     &2.287787d+06,2.306932d+06,2.326237d+06,2.345703d+06,2.365332d+06,
     &2.385126d+06,2.405085d+06,2.425211d+06,2.445505d+06,2.465970d+06,
     &2.486605d+06,2.507414d+06,2.528396d+06,2.549554d+06,2.570889d+06,
     &2.592403d+06,2.614096d+06,2.635971d+06,2.658030d+06,2.680272d+06,
     &2.702701d+06,2.725318d+06,2.748124d+06,2.771121d+06,2.794310d+06,
     &2.817693d+06,2.841272d+06,2.865048d+06,2.889023d+06,2.913199d+06,
     &2.937577d+06,2.962159d+06,2.986947d+06,3.011942d+06,3.037147d+06,
     &3.062562d+06,3.088190d+06,3.114032d+06,3.140091d+06,3.166368d+06,
     &3.192864d+06,3.219583d+06,3.246525d+06,3.273692d+06,3.301087d+06,
     &3.328711d+06,3.356566d+06,3.384654d+06,3.412978d+06,3.441538d+06,
     &3.470337d+06,3.499377d+06,3.528661d+06,3.558189d+06,3.587965d+06,
     &3.617989d+06,3.648265d+06,3.678794d+06,3.709579d+06,3.740621d+06,
     &3.771924d+06,3.803488d+06,3.835316d+06,3.867410d+06,3.899773d+06,
     &3.932407d+06,3.965314d+06,3.998497d+06,4.031957d+06,4.065697d+06,
     &4.099719d+06,4.134026d+06,4.168620d+06,4.203504d+06,4.238679d+06,
     &4.274149d+06,4.309916d+06,4.345982d+06,4.382350d+06,4.419022d+06,
     &4.456001d+06,4.493290d+06,4.530890d+06,4.568805d+06,4.607038d+06,
     &4.645590d+06,4.684465d+06,4.723666d+06,4.763194d+06,4.803053d+06,
     &4.843246d+06,4.883775d+06,4.924643d+06,4.965853d+06,5.007408d+06/
      data eg20t/
     &5.049311d+06,5.091564d+06,5.134171d+06,5.177135d+06,5.220458d+06,
     &5.264143d+06,5.308195d+06,5.352614d+06,5.397406d+06,5.442572d+06,
     &5.488116d+06,5.534042d+06,5.580351d+06,5.627049d+06,5.674137d+06,
     &5.721619d+06,5.769498d+06,5.817778d+06,5.866462d+06,5.915554d+06,
     &5.965056d+06,6.014972d+06,6.065307d+06,6.116062d+06,6.167242d+06,
     &6.218851d+06,6.270891d+06,6.323367d+06,6.376282d+06,6.429639d+06,
     &6.483443d+06,6.537698d+06,6.592406d+06,6.647573d+06,6.703200d+06,
     &6.759294d+06,6.815857d+06,6.872893d+06,6.930406d+06,6.988401d+06,
     &7.046881d+06,7.105850d+06,7.165313d+06,7.225274d+06,7.285736d+06,
     &7.346704d+06,7.408182d+06,7.470175d+06,7.532687d+06,7.595721d+06,
     &7.659283d+06,7.723377d+06,7.788008d+06,7.853179d+06,7.918896d+06,
     &7.985162d+06,8.051983d+06,8.119363d+06,8.187308d+06,8.255820d+06,
     &8.324906d+06,8.394570d+06,8.464817d+06,8.535652d+06,8.607080d+06,
     &8.679105d+06,8.751733d+06,8.824969d+06,8.898818d+06,8.973284d+06,
     &9.048374d+06,9.124092d+06,9.200444d+06,9.277435d+06,9.355070d+06,
     &9.433354d+06,9.512294d+06,9.591895d+06,9.672161d+06,9.753099d+06,
     &9.834715d+06,9.917013d+06,1.000000d+07,1.008368d+07,1.016806d+07,
     &1.025315d+07,1.033895d+07,1.042547d+07,1.051271d+07,1.060068d+07,
     &1.068939d+07,1.077884d+07,1.086904d+07,1.095999d+07,1.105171d+07/
      data eg20u/
     &1.114419d+07,1.123745d+07,1.133148d+07,1.142631d+07,1.152193d+07,
     &1.161834d+07,1.171557d+07,1.181360d+07,1.191246d+07,1.201215d+07,
     &1.211267d+07,1.221403d+07,1.231624d+07,1.241930d+07,1.252323d+07,
     &1.262802d+07,1.273370d+07,1.284025d+07,1.294770d+07,1.305605d+07,
     &1.316531d+07,1.327548d+07,1.338657d+07,1.349859d+07,1.361155d+07,
     &1.372545d+07,1.384031d+07,1.395612d+07,1.407291d+07,1.419068d+07,
     &1.430943d+07,1.442917d+07,1.454991d+07,1.467167d+07,1.479444d+07,
     &1.491825d+07,1.504309d+07,1.516897d+07,1.529590d+07,1.542390d+07,
     &1.555297d+07,1.568312d+07,1.581436d+07,1.594670d+07,1.608014d+07,
     &1.621470d+07,1.635039d+07,1.648721d+07,1.662518d+07,1.676430d+07,
     &1.690459d+07,1.704605d+07,1.718869d+07,1.733253d+07,1.747757d+07,
     &1.762383d+07,1.777131d+07,1.792002d+07,1.806998d+07,1.822119d+07,
     &1.837367d+07,1.852742d+07,1.868246d+07,1.883880d+07,1.899644d+07,
     &1.915541d+07,1.931570d+07,1.947734d+07,1.964033d+07/
      data eg21a/
     &1.000010d-05,1.100000d-04,3.000000d-03,5.500100d-03,1.000000d-02,
     &1.500000d-02,2.000000d-02,3.000000d-02,3.200000d-02,3.238000d-02,
     &4.300000d-02,5.900100d-02,7.700100d-02,9.500000d-02,1.000000d-01,
     &1.150000d-01,1.340000d-01,1.600000d-01,1.890000d-01,2.200000d-01,
     &2.480000d-01,2.825000d-01,3.145000d-01,3.520000d-01,3.910100d-01,
     &4.139900d-01,4.330000d-01,4.850100d-01,5.315800d-01,5.400100d-01,
     &6.250100d-01,6.825600d-01,7.050000d-01,7.900100d-01,8.600100d-01,
     &8.764200d-01,9.300100d-01,9.860100d-01,1.010000d+00,1.035000d+00,
     &1.070000d+00,1.080000d+00,1.090000d+00,1.110000d+00,1.125400d+00,
     &1.170000d+00,1.235000d+00,1.305000d+00,1.370000d+00,1.440000d+00,
     &1.445000d+00,1.510000d+00,1.590000d+00,1.670000d+00,1.755000d+00,
     &1.840000d+00,1.855400d+00,1.930000d+00,2.020000d+00,2.130000d+00,
     &2.360000d+00,2.372400d+00,2.767900d+00,3.059000d+00,3.380700d+00,
     &3.927900d+00,4.129200d+00,4.470000d+00,4.670000d+00,5.043500d+00,
     &5.623000d+00,6.160100d+00,6.476000d+00,7.079000d+00,7.524000d+00,
     &7.943000d+00,8.315300d+00,8.913000d+00,9.189800d+00,1.000000d+01,
     &1.067700d+01,1.122400d+01,1.259000d+01,1.371000d+01,1.522700d+01,
     &1.674500d+01,1.760300d+01,1.902800d+01,2.045200d+01,2.260300d+01,
     &2.498000d+01,2.791800d+01,2.920300d+01,3.051100d+01,3.388900d+01/
      data eg21b/
     &3.726700d+01,3.981000d+01,4.551700d+01,4.785100d+01,5.012000d+01,
     &5.559500d+01,6.144200d+01,6.310000d+01,6.790400d+01,7.079000d+01,
     &7.889300d+01,8.527700d+01,9.166100d+01,1.013000d+02,1.122000d+02,
     &1.300700d+02,1.367400d+02,1.585000d+02,1.670200d+02,1.778000d+02,
     &2.039900d+02,2.144500d+02,2.430100d+02,2.753600d+02,3.043200d+02,
     &3.535800d+02,3.981000d+02,4.540000d+02,5.144600d+02,5.829500d+02,
     &6.310000d+02,6.772900d+02,7.079000d+02,7.485200d+02,8.482000d+02,
     &9.611200d+02,1.010400d+03,1.116700d+03,1.234100d+03,1.363900d+03,
     &1.507300d+03,1.584600d+03,1.795600d+03,2.034700d+03,2.113000d+03,
     &2.248700d+03,2.371000d+03,2.485200d+03,2.612600d+03,2.661000d+03,
     &2.746500d+03,2.818000d+03,3.035400d+03,3.162000d+03,3.354600d+03,
     &3.548000d+03,3.707400d+03,3.981000d+03,4.307400d+03,4.642900d+03,
     &5.004500d+03,5.530800d+03,6.267300d+03,7.101700d+03,7.465900d+03,
     &8.251000d+03,9.118800d+03,1.007800d+04,1.113800d+04,1.170900d+04,
     &1.272600d+04,1.383200d+04,1.503400d+04,1.585000d+04,1.661600d+04,
     &1.778000d+04,1.930500d+04,1.995000d+04,2.054000d+04,2.113000d+04,
     &2.187500d+04,2.239000d+04,2.304000d+04,2.357900d+04,2.417600d+04,
     &2.441000d+04,2.478800d+04,2.512000d+04,2.585000d+04,2.605800d+04,
     &2.661000d+04,2.700000d+04,2.738000d+04,2.818000d+04,2.850000d+04/
      data eg21c/
     &2.901000d+04,2.985000d+04,3.073000d+04,3.162000d+04,3.182800d+04,
     &3.430700d+04,3.697900d+04,4.086800d+04,4.358900d+04,4.630900d+04,
     &4.939200d+04,5.247500d+04,5.516600d+04,5.656200d+04,6.172500d+04,
     &6.737900d+04,7.200000d+04,7.499000d+04,7.950000d+04,8.229700d+04,
     &8.250000d+04,8.651700d+04,9.803700d+04,1.110900d+05,1.167900d+05,
     &1.227700d+05,1.290700d+05,1.356900d+05,1.426400d+05,1.499600d+05,
     &1.576400d+05,1.657300d+05,1.742200d+05,1.831600d+05,1.925500d+05,
     &2.024200d+05,2.128000d+05,2.237100d+05,2.351800d+05,2.472400d+05,
     &2.732400d+05,2.872500d+05,2.945200d+05,2.972000d+05,2.985000d+05,
     &3.019700d+05,3.337300d+05,3.688300d+05,3.877400d+05,4.076200d+05,
     &4.504900d+05,5.234000d+05,5.502300d+05,5.784400d+05,6.081000d+05,
     &6.392800d+05,6.720600d+05,7.065100d+05,7.427400d+05,7.808200d+05,
     &8.208500d+05,8.629400d+05,9.071800d+05,9.616400d+05,1.002600d+06,
     &1.108000d+06,1.164800d+06,1.224600d+06,1.287300d+06,1.353400d+06,
     &1.422700d+06,1.495700d+06,1.572400d+06,1.653000d+06,1.737700d+06,
     &1.826800d+06,1.920500d+06,2.019000d+06,2.122500d+06,2.231300d+06,
     &2.306900d+06,2.345700d+06,2.365300d+06,2.385200d+06,2.466000d+06,
     &2.592400d+06,2.725300d+06,2.865000d+06,3.011900d+06,3.166400d+06,
     &3.328700d+06,3.678800d+06,4.065700d+06,4.493300d+06,4.723700d+06/
      data eg21d/
     &4.965900d+06,5.220500d+06,5.488100d+06,5.769500d+06,6.065300d+06,
     &6.376300d+06,6.592400d+06,6.703200d+06,7.046900d+06,7.408200d+06,
     &7.788000d+06,8.187300d+06,8.607100d+06,9.048400d+06,9.512300d+06,
     &1.000000d+07,1.051300d+07,1.105200d+07,1.161800d+07,1.221400d+07,
     &1.284000d+07,1.349900d+07,1.384000d+07,1.419100d+07,1.455000d+07,
     &1.491800d+07,1.568300d+07,1.648700d+07,1.690500d+07,1.733300d+07,
     &1.964000d+07/
      data eg22a/
     &1.000010d-05,3.000000d-03,5.000000d-03,6.900000d-03,1.000000d-02,
     &1.500000d-02,2.000000d-02,2.500000d-02,3.000000d-02,3.500000d-02,
     &4.200000d-02,5.000000d-02,5.800000d-02,6.700000d-02,7.700000d-02,
     &8.000000d-02,9.500000d-02,1.000000d-01,1.150000d-01,1.340000d-01,
     &1.400000d-01,1.600000d-01,1.800000d-01,1.890000d-01,2.200000d-01,
     &2.480000d-01,2.800000d-01,3.000000d-01,3.145000d-01,3.200000d-01,
     &3.500000d-01,3.910000d-01,4.000000d-01,4.330000d-01,4.850000d-01,
     &5.000000d-01,5.400000d-01,6.250000d-01,7.050000d-01,7.800000d-01,
     &7.900000d-01,8.500000d-01,8.600000d-01,9.100000d-01,9.300000d-01,
     &9.500000d-01,9.720000d-01,9.860000d-01,9.960000d-01,1.020000d+00,
     &1.035000d+00,1.045000d+00,1.071000d+00,1.097000d+00,1.110000d+00,
     &1.123000d+00,1.150000d+00,1.170000d+00,1.235000d+00,1.300000d+00,
     &1.337500d+00,1.370000d+00,1.440000d+00,1.475000d+00,1.500000d+00,
     &1.590000d+00,1.670000d+00,1.755000d+00,1.840000d+00,1.930000d+00,
     &2.020000d+00,2.100000d+00,2.130000d+00,2.360000d+00,2.550000d+00,
     &2.600000d+00,2.720000d+00,2.767920d+00,3.300000d+00,3.380750d+00,
     &4.000000d+00,4.129250d+00,5.043477d+00,5.346430d+00,6.160116d+00,
     &7.523983d+00,8.315287d+00,9.189814d+00,9.905554d+00,1.122446d+01,
     &1.370959d+01,1.592827d+01,1.945484d+01,2.260329d+01,2.498050d+01/
      data eg22b/
     &2.760773d+01,3.051126d+01,3.372015d+01,3.726653d+01,4.016900d+01,
     &4.551744d+01,4.825160d+01,5.157802d+01,5.559513d+01,6.790405d+01,
     &7.567357d+01,9.166088d+01,1.367420d+02,1.486254d+02,2.039950d+02,
     &3.043248d+02,3.717032d+02,4.539993d+02,6.772874d+02,7.485183d+02,
     &9.142423d+02,1.010394d+03,1.234098d+03,1.433817d+03,1.507331d+03,
     &2.034684d+03,2.248673d+03,3.354626d+03,3.526622d+03,5.004514d+03,
     &5.530844d+03,7.465858d+03,9.118820d+03,1.113775d+04,1.503439d+04,
     &1.661557d+04,2.478752d+04,2.739445d+04,2.928300d+04,3.697864d+04,
     &4.086771d+04,5.516564d+04,6.737947d+04,8.229747d+04,1.110900d+05,
     &1.227734d+05,1.831564d+05,2.472353d+05,2.732372d+05,3.019738d+05,
     &4.076220d+05,4.504920d+05,4.978707d+05,5.502322d+05,6.081006d+05,
     &8.208500d+05,9.071795d+05,1.002588d+06,1.108032d+06,1.224564d+06,
     &1.353353d+06,1.652989d+06,2.018965d+06,2.231302d+06,2.465970d+06,
     &3.011942d+06,3.678794d+06,4.493290d+06,5.488116d+06,6.065307d+06,
     &6.703200d+06,8.187308d+06,1.000000d+07,1.1618343d+07,
     &1.3840307d+07,1.4918247d+07,1.733253d+07,1.964033d+07/
      data eg23a/
     &1.000010d-05,1.000010d-01,4.139940d-01,5.315790d-01,6.825600d-01,
     &8.764250d-01,1.123000d+00,1.440000d+00,1.855390d+00,2.382370d+00,
     &3.059020d+00,3.927860d+00,5.043480d+00,6.475950d+00,8.315290d+00,
     &1.067700d+01,1.370960d+01,1.760350d+01,2.260330d+01,2.902320d+01,
     &3.726650d+01,4.785120d+01,6.144210d+01,7.889320d+01,1.013010d+02,
     &1.300730d+02,1.670170d+02,2.144540d+02,2.753640d+02,3.535750d+02,
     &4.539990d+02,5.829470d+02,7.485180d+02,9.611170d+02,1.234100d+03,
     &1.584610d+03,2.034680d+03,2.248670d+03,2.485170d+03,2.612590d+03,
     &2.746540d+03,3.035390d+03,3.354630d+03,3.707440d+03,4.307420d+03,
     &5.530840d+03,7.101740d+03,9.118820d+03,1.059460d+04,1.170880d+04,
     &1.503440d+04,1.930450d+04,2.187490d+04,2.357860d+04,2.417550d+04,
     &2.478750d+04,2.605840d+04,2.700010d+04,2.850110d+04,3.182780d+04,
     &3.430670d+04,4.086770d+04,4.630920d+04,5.247520d+04,5.656220d+04,
     &6.737950d+04,7.202450d+04,7.949870d+04,8.250340d+04,8.651700d+04,
     &9.803650d+04,1.110900d+05,1.167860d+05,1.227730d+05,1.290680d+05,
     &1.356860d+05,1.426420d+05,1.499560d+05,1.576440d+05,1.657270d+05,
     &1.742240d+05,1.831560d+05,1.925470d+05,2.024190d+05,2.127970d+05,
     &2.237080d+05,2.351770d+05,2.472350d+05,2.732370d+05,2.872460d+05,
     &2.945180d+05,2.972110d+05,2.984910d+05,3.019740d+05,3.337330d+05/
      data eg23b/
     &3.688320d+05,3.877420d+05,4.076220d+05,4.504920d+05,4.978710d+05,
     &5.233970d+05,5.502320d+05,5.784430d+05,6.081010d+05,6.392790d+05,
     &6.720550d+05,7.065120d+05,7.427360d+05,7.808170d+05,8.208500d+05,
     &8.629360d+05,9.071800d+05,9.616720d+05,1.002590d+06,1.108030d+06,
     &1.164840d+06,1.224560d+06,1.287350d+06,1.353350d+06,1.422740d+06,
     &1.495690d+06,1.572370d+06,1.652990d+06,1.737740d+06,1.826840d+06,
     &1.920500d+06,2.018970d+06,2.122480d+06,2.231300d+06,2.306930d+06,
     &2.345700d+06,2.365330d+06,2.385130d+06,2.465970d+06,2.592400d+06,
     &2.725320d+06,2.865050d+06,3.011940d+06,3.166370d+06,3.328710d+06,
     &3.678790d+06,4.065700d+06,4.493290d+06,4.723670d+06,4.965850d+06,
     &5.220460d+06,5.488120d+06,5.769500d+06,6.065310d+06,6.376280d+06,
     &6.592410d+06,6.703200d+06,7.046880d+06,7.408180d+06,7.788010d+06,
     &8.187310d+06,8.607080d+06,9.048370d+06,9.512290d+06,1.000000d+07,
     &1.051270d+07,1.105170d+07,1.161830d+07,1.221400d+07,1.252320d+07,
     &1.284030d+07,1.349860d+07,1.384030d+07,1.419070d+07,1.454990d+07,
     &1.491820d+07,1.568310d+07,1.648720d+07,1.690460d+07,1.733250d+07,
     &1.964030d+07/
*b groupr.1837
c
c     ***ecco  33-group structure
      else if (ign.eq.19) then
         ng=33
         do ig=1,34
            eg(ig)=eg19(ig)
         enddo
c
c     ***ecco 1968-group structure
      else if (ign.eq.20) then
         ng=1968
         do ig=1,1969
            eg(ig)=eg20(ig)
         enddo
c
c     ***tripoli 315-group structure
      else if (ign.eq.21) then
         ng=315
         do ig=1,316
          eg(ig)=eg21(ig)
         enddo
c
c     ***xmas LWPC 172-group structure
      else if (ign.eq.22) then
         ng=172
         do ig=1,173
            eg(ig)=eg22(ig)
         enddo
c
c     ***vit-J LWPC 175-group structure
      else if (ign.eq.23) then
         ng=175
         do ig=1,176
            eg(ig)=eg23(ig)
         enddo
*b groupr.1888
      if (ign.eq.19) write(nsyso,'(/
     &  '' neutron group structure......ecco  33-group'')')
      if (ign.eq.20) write(nsyso,'(/
     &  '' neutron group structure......ecco  1968-group'')')
      if (ign.eq.21) write(nsyso,'(/
     &  '' neutron group structure......tripoli 315-group'')')
      if (ign.eq.22) write(nsyso,'(/
     &  '' neutron group structure......xmas LWPC 172-group'')')
      if (ign.eq.23) write(nsyso,'(/
     &  '' neutron group structure......vit-J LWPC 175-group'')')
*/
*/ increase the size of egn from 641 to 15000 i.e DICE 13193 APOLLO 11276
*/
*d groupr.1643
      data ngmax/15000/
*d groupr.228
      common/groupn/ign,ngn,egn(15000)
*d groupr.772
      common/groupn/ign,ngn,egn(15000)
*d groupr.1303
      common/groupn/ign,ng,eg(15000)
*d groupr.2643
      common/groupn/ig,ngn,egn(15000)
*d groupr.3074
      common/groupn/ign,ngn,egn(15000)
*d groupr.4274
      common/groupn/ign,ngn,egn(15000)
*d groupr.6415
      common/groupn/ign,ngn,egn(15000)
*d groupr.6918
      common/groupn/ign,ngn,egn(15000)
*/
*ident wlup8x
*/ wimsr
*/ Prepare (n,2n) from MT=875-891 to correct absorption
*/ D. L. Aldama, NDS/IAEA Consultant, 2005
*/
*i wimsr.895
      jn2n=0
*i wimsr.984
      if (mth.eq.16) jn2n=16
*i wimsr.985
      if ((mth.ge.875.and.mth.le.891).and.(jn2n.ne.16)) go to 236
*/
*ident ads01
*/ matxsr - D. L. Aldama, NDS/IAEA Consultant, July 2005
*/ Need more space  (for processing ORNL 421 energy group structure)
*/  subroutine vector
*d up88.21
      dimension b(30000)
*d up88.23
      maxb=30000
*ident ads02
*/ matxsr - D. L. Aldama, NDS/IAEA Consultant, July 2005
*/ Remove a problem in case of multi-temperature matxs files processing.
*/ Scratch tapes iref and nscr should be simultaneously forwarded.
*/ The problem seems to affect only multi-temperature runs, if the mf=6
*/ of gendf tape was not ordered by mt numbers.
*/
*i matxsr.1870
            if (iref.ne.0) call contio(iref,0,0,b(1),nb,nw)
*d matxsr.1878
            if (iskip.eq.0) then
              if (iref.ne.0) call tosend(iref,0,0,b(1))
              call tosend(nscr,0,0,b(1))
            endif
*/
*ident ads03
*/ acer - D. L. Aldama, NDS/IAEA Consultant, July 2005
*/ need more storage for processing jeff-3.1 U-238
*/ ADS library
*/
*d up57.76
      max3=4000000
*d up57.78
      common/xsst/xss(4000000),n3
*d up57.80
      common/xsst/xss(4000000),n3
*d up57.82
      common/xsst/xss(4000000),n3
*d up57.84
      common/xsst/xss(4000000),n3
*d up57.86
      common/xsst/xss(4000000),n3
*d up57.88
      common/xsst/xss(4000000),n3
*d up57.90
      common/xsst/xss(4000000),n3
*d up57.92
      common/xsst/xss(4000000),n3
*d up57.94
      common/xsst/xss(4000000),n3
*d up57.96
      common/xsst/xss(4000000),n3
*d up57.98
      common/xsst/xss(4000000),n3
*d up57.100
      common/xsst/xss(4000000),n3
*d up57.102
      common/xsst/xss(4000000),n3
*d up57.104
      common/xsst/xss(4000000),n3
*d up57.106
      common/xsst/xss(4000000),n3
*d up57.108
      common/xsst/xss(4000000),n3
*d up57.110
      common/xsst/xss(4000000),n3
*d up57.112
      common/xsst/xss(4000000),n3
*d up57.114
      common/xsst/xss(4000000),n3
*d up57.116
      common/xsst/xss(4000000),n3
*d up57.118
      common/xsst/xss(4000000),n3
*d up57.120
      common/xsst/xss(4000000),n3
*d up57.122
      common/xsst/xss(4000000),n3
*d up57.124
      common/xsst/xss(4000000),n3
*d up57.126
      common/xsst/xss(4000000),n3
*d up57.128
      common/xsst/xss(4000000),n3
*d up57.130
      common/xsst/xss(4000000),n3
*d up57.132
      common/xsst/xss(4000000),n3
*d up57.134
      common/xsst/xss(4000000),n3
*d up57.136
      common/xsst/xss(4000000),n3
*d up57.138
      common/xsst/xss(4000000),n3
*d up57.140
      common/xsst/xss(4000000),n3
*d up57.142
      common/xsst/xss(4000000),n3
*d up70.95
      common/xsst/xss(4000000),n3
*d up57.144
      common/xsst/xss(4000000),n3
*d up70.222
      common/xsst/xss(4000000),n3
*d up57.146
      common/xsst/xss(4000000),n3
*d up57.148
      common/xsst/xss(4000000),n3
*d up57.150
      common/xsst/xss(4000000),n3
*/