anchor.diff:         !
      end select
      !
  -   if (Nspecies>0) then
  +   if (lineprofile_do) then
         write(out,"(10x,'Pressure = ',e18.7)") pressure
         write(out,"(10x,'Voigt parameters:  gamma       n         T0            P0')")
         do i =1,Nspecies
  
positive.code:    !
     nswap = 0
     Noffset = max(nint(offset/dfreq,ik),1)
     !
     if (verbose>=3) print('(/"Transition files ... "/)')
     !
  
negative.code:     !
      implicit none
      !
      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
      type(HitranErrorT),pointer :: HITRAN
  