anchor.diff:     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))
  +      gamma_comb(j,i-j) = get_Voigt_gamma_val(Nspecies,real(i,rk),real(j,rk))
        enddo
      enddo
      !
  
positive.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
  
negative.code:         !
          do itemp = 1,npoints
            !
            temp0 = real(itemp,rk)*dtemp
            !
            write(tunit,"(f8.1,1x,f15.4)") temp0,pf(0,itemp)
  