anchor.diff:            call read_line(eof) ; if (eof) exit
             !
             call readu(w)
  +          !
  +          if_species_defined = .true.
             !
             do while (trim(w)/="".and.trim(w)/="END")
               !
  
positive.code:   type(selectT),save :: upper(filtermax),lower(filtermax)
    type(QNT),save :: QN
    !
    integer(ik)    :: Nspecies = 0, Nfilters = 0
    type(speciesT),save :: species(nspecies_max)
    type(HitranErrorT),target,save :: HITRAN_E(HITRAN_max_ierr),HITRAN_S(HITRAN_max_ierr)
    type(HitranErrorT),target,save :: HITRAN_Air,HITRAN_Self,HITRAN_N,HITRAN_Delta
  
negative.code:           !
            allocate(cooling(npoints),stat=info)
            call ArrayStart('cooling',info,size(cooling),kind(cooling))
            cooling = 0
            !
          endif
          !
          pf = 0
          !
        endif
        j0 = 0
  