anchor.diff:               !
                if (tranfreq<small_) cycle
                !
  -             abscoef_ram(iswap)=cmcoef*acoef*gf*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
  -             nu_ram(iswap) = tranfreq
  +             iswap_ = iswap_ + 1
  +             !
  +             abscoef_ram(iswap_)=cmcoef*acoef*gf*exp(-beta*energyi)*(1.0_rk-exp(-beta*tranfreq))/(tranfreq**2*partfunc)
  +             !
  +             nu_ram(iswap_) = tranfreq
  +             !
  +             ! estimate the line shape parameter 
  +             !
  +             if (lineprofile_do) then
  +               !
  +               ! The current version of ExoCross does know how to locate J-values in HITRAN-input. 
  +               ! Therefore we can only use the static values of Voigt-parameters together with Nspecies
  +               !
  +               Jf = 1000.0_rk
  +               Ji = 1000.0_rk
  +               if (Nspecies==0) then
  +                 species(1)%gamma = gamma_
  +                 species(1)%n = n_
  +                 species(2)%gamma = gamma_s
  +                 species(2)%n = 0
  +               endif
  +               gamma_RAM(iswap_) = get_Voigt_gamma_n(Nspecies_,Jf,Ji)
  +             endif
                !
                cycle
                !
              119 continue
  -             nswap_ = iswap-1
  +             nswap = iswap_-1
                exit
              enddo
              !
  -           if (nswap_<1) then
  +           nswap = iswap_
  +           !
  +           if (nswap<1) then
                 write(out,"('... Finish swap')") 
                 exit loop_tran
              endif
  -           !
  -           nswap = nswap_
              !
           else ! ExoMol format of the trans-file
              !
  
positive.code:       !
        halfwidth = 0 
        !
        do i=1,Nspecies
          !
          halfwidth =  halfwidth + species(i)%ratio*species(i)%gamma*(species(i)%T0/Temp)**species(i)%N*pressure/species(i)%P0
  
negative.code:   type(HitranErrorT),target :: HITRAN_E(HITRAN_max_ierr),HITRAN_S(HITRAN_max_ierr),HITRAN_Air,HITRAN_Self,HITRAN_N,HITRAN_Delta
    !
    logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., stick_hitran = .false.
    logical :: lineprofile_do = .false.
    logical :: microns = .false.
    logical :: use_resolving_power = .false.  ! using resolving for creating the grid
    !
  