anchor.diff:            !
             call readf(enermax)
             !
  -       case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0','LOREN','LORENTZIAN','LORENTZ','MAX')
  +       case('GAUSSIAN','GAUSS','DOPPL','DOPPLER','RECT','BOX','BIN','STICKS','STICK','GAUS0','DOPP0','LOREN','LORENTZIAN','LORENTZ','MAX','VOIGT','HITRAN')
             !
             proftype = trim(w)
             !
  +          if (trim(w(1:5))=="LOREN") ioffset = 500
  +          if (trim(w(1:5))=="VOIGT") ioffset = 500
  +          !
          case ("HWHM","HALFWIDTH")
             !
             call readf(halfwidth)
             !
  -       case ("MASSES","MASS")
  -          !
  -          call readf(m1)
  -          call readf(m2)
  +       case ("IOFFSET")
  +          !
  +          call readi(ioffset)
  +          !
  +       case ("MASS")
  +          !
  +          call readf(meanmass)
  +          !
  +       case ("VOIGT_GAMMA")
  +          !
  +          call readf(voigt_gamma)
  +          !
  +       case ("VOIGT_N")
  +          !
  +          call readf(voigt_N)
             !
         case default
           !
  
positive.code:           halfwidth=dpwcoef*tranfreq
          endif 
          !
          !   if transition frequency is out of selected range
          if (tranfreq>freqr+500.0*halfwidth.or.tranfreq<freql-500.0*halfwidth) cycle
          !
          x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
          !
  
negative.code:             enddo
              !$omp end parallel do 
              !
          case ('VOIGT')
              !
              ib =  max(nint( ( tranfreq-halfwidth*ioffset-freql)/dfreq )+1,1)
              ie =  min(nint( ( tranfreq+halfwidth*ioffset-freql)/dfreq )+1,npoints)
              !
              !$omp parallel do private(ipoint,tranfreq_i) shared(intens) schedule(dynamic)
              do ipoint=ib,ie
                 !
                 tranfreq_i = freq(ipoint)
                 !
                 intens(ipoint)=intens(ipoint)+voigt_humlicek(tranfreq_i,tranfreq,halfwidth,halfwidth0)*abscoef
                 !
              enddo
              !$omp end parallel do             !
          case ('MAX');
            !
            ipoint =  max(nint( ( tranfreq-freql)/dfreq )+1,1)
  