anchor.diff:               !
                if (lineprofile_do) then
                  gamma_RAM(iswap)=halfwidth
  +               gamma_idx_RAM(iswap) = 0
                endif
                !
                cycle
  
positive.code:         call ArrayStart('ivib_state',info,size(ivib_state),kind(ivib_state))
          ivib_state = 1
          energies_vib = 0
          ivib1 = QN%vibcol(1)-3
          ivib2 = QN%vibcol(2)-3
        endif
        !
        if (trim(proftype)=='LIFETIME') THEN
          !
          ! allocate the matrix for the sum of the A-coeffs
          !
  
negative.code:   integer(ik) :: N_omp_procs=1
    !
    integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=2,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
    real(rk)      :: S_crit = 1e-29      ! cm/molecule, HITRAN cut-off paramater
    real(rk)      :: nu_crit = 2000.0d0  ! cm-1, HITRAN cut-off paramater
    real(rk)      :: resolving_power  = 1e6,resolving_f ! using resolving_power to set up grid
  