anchor.diff:    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
     !
  
positive.code:             halfwidth = 0
              !
              if ( trim(species(ispecies)%filename)/="" ) then
                !
                gamma_ = species(ispecies)%gammaQN(Jpp,Jp-Jpp)
                n_     = species(ispecies)%nQN(Jpp,Jp-Jpp)
  
negative.code:            !
             nfiles = 1 
             call reada(intfilename(1))
             !
             ! for J-dependent histograms a J-value is expected in the last column 
             !
             if (histogramJ)   call readi(intJvalue(1))
             !
             cycle 
             !
  