anchor.diff:            !
             histogram = .true.
             !
  +        case ("HISTOGRAM-J")
  +          !
  +          histogram = .true.
  +          histogramJ = .true.
  +          !
           case ("HITRAN")
             !
             hitran = .true.
  
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:            !
             read(tunit,*,end=20) tranfreq,abscoef
             !
             ! for pre-computed J-dependent integrated intensities 
             if (histogramJ) then 
               Ji = IntJvalue(i)
               Jf = Ji
             endif
             !
          elseif (hitran) then 
             !
             ! using pre-computed integrated intensities for a given Temperature and bin
  