anchor.diff:                  tranfreq = nu_ram(iswap)
                   halfwidth = gamma_ram(iswap)
                   !
  -                call do_pseud(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,dpwcoef,intens_omp(:,iomp))
  +                call do_Voi_Q(tranfreq,abscoef,dfreq,freq,abciss,weight,halfwidth,offset,freql,dpwcoef,intens_omp(:,iomp))
                   !
                 enddo
                 !
               enddo
               !$omp end parallel do
               !
  -            !
  -        !case ('PSE-R') ! pseudo_Voigt_Rocco_Cruzado_ActaPhysPol_2012.pdf
  -            !
  -            !
  -        !case ('PSE-L') ! pseudo_Voigt_Liu_Lin_JOptSocAmB_2001.pdf
  -            !
  -            !
  -        case ('VOI-Q') ! VOIGT-QUADRATURES
  -            !
  +        case ('VOI-9') ! VOIGT-Alg-916
               !
               !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth) shared(intens_omp) schedule(dynamic)
               do iomp = 1,N_omp_procs
  
positive.code:             enddo
              !$omp end parallel do
              !
          case ('GAUS0')
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  !
                  call do_gauss0(tranfreq,abscoef,dfreq,freq,halfwidth,freql,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
          case ('DOPPL')
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth0) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  halfwidth0=dpwcoef*tranfreq
                  if (halfwidth0<100.0_rk*small_) cycle
                  !
                  call do_gauss(tranfreq,abscoef,dfreq,freq,halfwidth0,offset,freql,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
              !
          case ('DOPP0')
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth0) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  halfwidth0=dpwcoef*tranfreq
                  if (halfwidth0<100.0_rk*small_) cycle
                  !
                  call do_gauss0(tranfreq,abscoef,dfreq,freq,halfwidth0,freql,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
          case ('GAUSS')
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
                !
                do iswap = iomp,nswap,N_omp_procs
                  !
                  abscoef = abscoef_ram(iswap)
                  tranfreq = nu_ram(iswap)
                  !
                  call do_gauss(tranfreq,abscoef,dfreq,freq,halfwidth,offset,freql,intens_omp(:,iomp))
                  !
                enddo
                !
              enddo
              !$omp end parallel do
              !
              !
          case ('LOREN')
              !
              !$omp parallel do private(iomp,iswap,abscoef,tranfreq,halfwidth) shared(intens_omp) schedule(dynamic)
              do iomp = 1,N_omp_procs
  
negative.code:   integer(ik)   :: GNS=1,npoints=1001,nchar=1,nfiles=1,ipartf=0,verbose=2,ioffset = 10,iso=1,imolecule =-1
    real(rk)      :: temp=298.0,partfunc=-1.0,partfunc_ref=-1.0,freql=-small_,freqr= 200000.0,thresh=1.0d-70
    real(rk)      :: voigt_gamma = 0.05, voigt_n = 0.44, offset = 25.0, pressure = 1.0_rk,halfwidth=1e-2,meanmass=1.0,maxtemp=10000.0
    real(rk)      :: lineshift=0
    real(rk)      :: enermax = 1e7, abscoef_thresh = 1.0d-50, abundance = 1.0d0, gf_factor = 1.0d0
    real(rk)      :: S_crit = 1e-29      ! cm/molecule, HITRAN cut-off paramater
    real(rk)      :: nu_crit = 2000.0d0  ! cm-1, HITRAN cut-off paramater
  