anchor.diff:         stop
        endif
        !
  +     ichunk = 0
  +     !
        loop_tran: do
           !
           ! finis this if the eof has been reached 
  
positive.code:   type(HitranErrorT),target,save :: HITRAN_E(HITRAN_max_ierr),HITRAN_S(HITRAN_max_ierr)
    type(HitranErrorT),target,save :: HITRAN_Air,HITRAN_Self,HITRAN_N,HITRAN_Delta
    !
    
    integer(ik),allocatable,save  :: gamma_idx(:,:) !Used for indexing the gamma for the fast_voigt
    real(rk),allocatable,save  :: gamma_comb(:,:) !Used for indexing the gamma for the fast_voigt
    logical :: partfunc_do = .true., filter = .false., histogram = .false., hitran_do = .false.,  histogramJ = .false., &
  
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)
  