anchor.diff:    integer(ik)    :: Nspecies = 0, Nfilters = 0
     type(speciesT) :: species(nspecies_max)
     !
  -  logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran = .false.
  +  logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran = .false.,  histogramJ = .false.
     !
     contains
     !
  
positive.code:            if (i>nfiles_max) call report ("Too many files, increase nfiles_max"//trim(w),.true.)
             !
             intfilename(i) = trim(v)
             !
             call read_line(eof) ; if (eof) exit
             call reada(w)
  
negative.code:   real(rk)      :: voigt_gamma = 0.05, voigt_n = 0.44, offset = -25.0, pressure = 1.0_rk 
    real(rk)      :: enermax = 1e6, abscoef_thresh = 1.0d-50
    integer(ik)   :: nquad = 20      ! Number of quadrature points 
    !
    character(len=cl) :: specttype="absorption",enrfilename="none",intfilename(nfiles_max),proftype="DOPPL",output="output"
    integer(ik)   :: intJvalue(nfiles_max)
    character(4) a_fmt
    character(9) b_fmt
    !
  