anchor.diff:       !
        !Combine all the gammas
        do i=0,JmaxAll
  -      do j= max(0,i-1),min(JmaxAll,i+1) 
  -       call fast_voigt%generate_indices(gamma_comb(i,i-j),gamma_idx(i,i-j),JmaxAll)
  +      do j= max(0,i-2),min(JmaxAll,i+2) 
  +       call fast_voigt%generate_indices(gamma_comb(j,i-j),gamma_idx(j,i-j),JmaxAll)
         enddo
        enddo  
        !   
  
positive.code:    gamma_comb = -1
     do i=0,JmaxAll
       do j= max(0,i-2),min(JmaxAll,i+2)
        gamma_comb(i,i-j) = get_Voigt_gamma_val(Nspecies,real(i,rk),real(j,rk))
       enddo
     enddo
     !
  
negative.code:                 tranfreq = nu_ram(iswap)
                  halfwidth = gamma_ram(iswap)
                  !
                  call do_Voigt(tranfreq,abscoef,use_resolving_power,freq,halfwidth,dpwcoef,offset,freql,intens_omp(:,iomp))
                  !
                enddo
  