anchor.diff:             do iswap = 1,N_to_RAM
                !
                read(tunit,*,end=118) nu_ram(iswap),abscoef_RAM(iswap)
  +             !
  +             if (lineprofile_do) then
  +               gamma_RAM(iswap)=halfwidth
  +             endif
                !
                cycle
                !
  
positive.code:   integer(ik),parameter :: max_transitions_to_ram = 1000000000
    integer(ik) :: N_omp_procs=1
    !
    integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=3,ioffset = 10,iso=1
    real(rk)      :: temp=298.0,partfunc=-1.0,freql=-small_,freqr= 20000.0,thresh=1.0d-70,halfwidth=1e-2,meanmass=1.0,maxtemp=10000.0
    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, abundance = 1.0d0
  
negative.code:            !
           case ('J','A0')
             !
             read(bunit,*,end=15) ch_broad(1:3),gamma_,n_,J
             species(i)%gammaQN(J,:) = gamma_
             species(i)%nQN(J,:) = n_
             !
           case ('JJ','A1')
             !
             read(bunit,*,end=15) ch_broad(1:3),gamma_,n_,Jpp,Jp
             !
             if (abs(Jp-Jpp)>1) stop 'Jp-Jpp in broadening file'
             !
  