anchor.diff:         !
          Jmax = JmaxAll
          !
  -       allocate(species(i)%gammaQN(0:Jmax,-1:1),stat=info)
  +       allocate(species(i)%gammaQN(0:Jmax,-2:2),stat=info)
          call ArrayStart('gammaQN',info,size(species(i)%gammaQN),kind(species(i)%gammaQN))
          species(i)%gammaQN(:,:) = species(i)%gamma
          !
  -       allocate(species(i)%nQN(0:Jmax,-1:1),stat=info)
  +       allocate(species(i)%nQN(0:Jmax,-2:2),stat=info)
          !         
          call ArrayStart('nQN',info,size(species(i)%nQN),kind(species(i)%nQN))
          species(i)%nQN(:,:) = species(i)%N
  
positive.code:     !
      logical :: eof,if_halfwidth_defined = .false., if_species_defined = .false., if_QN_defined = .false.
      character(len=cl) :: w,v
      integer(ik)   :: i,ifilter,iE,iS,npoints0,ierror,igrid
      type(HitranErrorT),pointer :: HITRAN
      real(rk) :: f_t
  
negative.code:            !
             if (i>nfiles_max) call report ("Too many files, increase nfiles_max"//trim(w),.true.)
             !
             intfilename(i) = trim(vl)
             !
             ! for J-dependent histograms a J-value is expected in the last column
             !
  