def fvoigt(damp, vv):
    """
	Extract spectral data from the Kitt Peak FTS-Spectral-Atlas
	as provided by H. Neckel, Hamburg.

	INPUTS: 
			DAMP: A scalar with the damping parameter
			VV: Wavelength axis usually in Doppler units.

	OUTPUTS: 
			H: Voigt function
			F: Faraday-Voigt function

	NOTES: 
			A rational approximation to the complex error function is used
			after Hui, Armstrong, and Wray(1978, JQSRT 19, 509). H and F are 
			the real and imaginary parts of such function, respectively.
			The procedure is inspired on that in SIR (Ruiz Cobo & del Toro 
			Iniesta 1992, ApJ 398, 385). On its turn, that routine was taken
			from modifications by A. Wittmann (1986) to modifications by S.K.
			Solanki (1985) to an original FORTRAN routine written by J.W. Harvey
			and A. Nordlund.
	"""
    import numpy as np
    A = [122.60793177710433, 214.38238869470644, 181.92853309218154, 
        93.15558045813845, 30.18014219621059, 5.912626209773153, 
        0.564189583562615]
    B = [122.60793177387535, 352.73062511096356, 457.33447878389774, 
        348.7039177194958, 170.35400182109146, 53.99290691294021, 
        10.4798571142604, 1.0]
    z = np.array(damp * np.ones(len(vv)) + -abs(vv) * 1.0j)
    Z = ((((((A[6] * z + A[5]) * z + A[4]) * z + A[3]) * z + A[2]) * z + A[
        1]) * z + A[0]) / (((((((z + B[6]) * z + B[5]) * z + B[4]) * z + B[
        3]) * z + B[2]) * z + B[1]) * z + B[0])
    h = np.real(Z)
    f = np.sign(vv) * np.imag(Z) * 0.5
    return [h, f]
