anchor.diff:     do i = 1,nfiles
        !
        open(unit=tunit,file=trim(intfilename(i)),action='read',status='old')
  -     do
  +     loop_tran: do
           !
           if (histogram) then 
              !
  
positive.code:           filter  = .true.
            !
            call read_line(eof) ; if (eof) exit
            call readu(w)
            !
            do while (trim(w)/="".and.trim(w)/="END")
              !
              select case(w)
              !
              case('LOWER')
                !
                call readi(lower%i) ; lower%i = lower%i - 4
                call reada(lower%mask)
                !
              case('UPPER')
                !
                call readi(upper%i) ; upper%i = upper%i - 4
                call reada(upper%mask)
                !
              case default
                !
                call report ("Unrecognized unit name "//trim(w),.true.)
                !
              end select
              !
              call read_line(eof) ; if (eof) exit
              call readu(w)
              !
            enddo
            !
            if (trim(w)/="".and.trim(w)/="END") then
               !
  
negative.code:     character(len=cl) :: name         ! Broadener number of quadrature points 
    end type speciesT
    !
    type(selectT) :: upper(filtermax),lower(filtermax)
    !
    integer(ik)    :: Nspecies = 0, Nfilters = 0
    type(speciesT) :: species(nspecies_max)
    !
    logical :: partfunc_do = .true., filter = .false., histogram = .false.
  