anchor.diff:      character(len=cl) :: name         ! Broadener number of quadrature points 
     end type speciesT
     !
  -  type(selectT) :: upper,lower
  +  type(selectT) :: upper(filtermax),lower(filtermax)
     !
  -  integer(ik)    :: Nspecies = 0
  +  integer(ik)    :: Nspecies = 0, Nfilters = 0
     type(speciesT) :: species(nspecies_max)
     !
     logical :: partfunc_do = .true., filter = .false., histogram = .false.
  
positive.code:         !
          cycle
       20  exit
       enddo
       !
       close(tunit)
       !
  
negative.code:    !
     use  input
     !
     integer(ik) :: info,ipoint,nlevels,i,itemp,enunit,tunit,sunit,j,ilog,ib,ie,j0,ilevelf,ileveli,indexi,indexf,iline,maxitems,kitem,l,nlines,iquad,ifilter
     real(rk)    :: beta,ln2,ln22,dtemp,dfreq,temp0,beta0,intband,dpwcoef,x0,tranfreq,tranfreq_i,abscoef,dfreq_,xp,xm,de,lor,b,lor2,dfreq_2,halfwidth0,dnu_half
     real(rk)    :: cmcoef,emcoef,energy,energyf,energyi,elow,jf,ji,acoef,j0rk,bnorm,f,eta1,eta2,intens1,intens2,Va0,gammaV,a,wg,d
     real(rk)    :: sigma,alpha,gamma,y,x1,x2,voigt_,dx2,xi,L1,L2
  