anchor.diff:     real(rk),allocatable :: acoef_RAM(:),abscoef_ram(:),nu_ram(:),intens_omp(:,:),gamma_ram(:)
      integer(ik),allocatable :: indexi_RAM(:),indexf_RAM(:)
      integer(ik),allocatable :: ileveli_RAM(:),ilevelf_RAM(:),gamma_idx_RAM(:)
  -   integer(ik) :: nswap_,nswap,iswap,iswap_,iomp
  +   integer(ik) :: nswap_,nswap,iswap,iswap_,iomp,ichunk
      real(rk),allocatable :: energies_vib(:)
      integer(ik),allocatable :: ivib_state(:)
      !
  
positive.code:    logical :: eof
     character(len=cl) :: w
     !
     integer(ik) :: imol,iostat_
     real(rk)    :: gf,gi,mem_t,temp_gamma_n
     character(55) ch_q,ch_broad
     !
  
negative.code:          energyi = energies(ileveli)
           tranfreq = energyf - energyi
           !
           lambda = 1e7/tranfreq ! wavelength in nm 
           !
           wl2ind = int(log(lambda)/ratiolog+0.5_rk)
           !
           !g_l*fij = c2 * Aji * g_u / (c*wn)**2
           !gf = 1.34738e+21*acoef*g_u / (c*tranfreq)**2
           !
           iloggf = int(log(gf)*(1.0_rk/(log(10.0_rk)*0.001_rk)),kind=2) + 2**14
           iener  = int(log(energyi)*(1.0_rk/(log(10.0_rk)*0.001_rk)),kind=2) + 2**14
           !
           jf = jrot(ilevelf)
           ji = jrot(ileveli)
  