anchor.diff:               Jf = Ji
              endif
              !
  -        elseif (hitran) then 
  +        elseif (hitran_do) then 
              !
              ! HITRAN format
              !
  
positive.code:           if (trim(w(1:3))=="PSE") ioffset = 500
            offset = 25.0_rk
            !
         case ("HWHM","HALFWIDTH")
            !
            call readf(halfwidth)
  
negative.code:     character(len=cl) :: model="const"  ! Broadening model
      real(rk),pointer  :: gammaQN(:,:) ! Voigt parameter gamma as a function of QN
      real(rk),pointer  :: nQN(:,:)     ! Voigt parameter n as a function of QN
      real(rk)          :: delta = 0    ! Pressure shift induced by air, referred to p=1 atm
      !
    end type speciesT
    !
    type HitranT
      integer(ik)  :: molec_id     = 1 ! The HITRAN integer ID for this molecule in all its isotopologue forms
      integer(ik)  :: molec_iso_id = 1 ! Integer ID of a particular Isotopologue, unique only to a given molecule, in order or abundance (1 = most abundant)
    end type HitranT
    !
    type(selectT) :: upper(filtermax),lower(filtermax)
    !
    integer(ik)    :: Nspecies = 0, Nfilters = 0
    type(speciesT) :: species(nspecies_max)
    type(hitranT) :: hitran
    !
    logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., stick_hitran = .false.
    !
    contains
    !
  