anchor.diff:            halfwidth=dpwcoef*tranfreq
           endif 
           !
  +        !   half width for Doppler profiling
  +        if (proftype(1:5)=='VOIGT') then
  +          halfwidth0=dpwcoef*tranfreq
  +        endif 
  +        !
           !   if transition frequency is out of selected range
  -        if (tranfreq>freqr+500.0*halfwidth.or.tranfreq<freql-500.0*halfwidth) cycle
  +        if (tranfreq>freqr+ioffset*halfwidth.or.tranfreq<freql-ioffset*halfwidth) cycle
           !
           x0 = sqrt(ln2)/halfwidth*dfreq*0.5_rk
           !
  
positive.code:    !
    end subroutine intensity
    !
  end module spectrum
  
  
  
  
negative.code:           !
            call readf(enermax)
            !
         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 ("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
          !
  