pro etarho_sgl,ipx, jpx, csn, eta, hgen, fgen, maxsz, er

;  modified from etarho_vec.pro for single profile

;	purpose:  Compute eta's & rho's 
;
;	THIS VERSION PROCESSES ARRAYS OF FIT PARAMETERS


;   INPUTS:
;   ipx,jpx = first, last wavelength pixels to compute
;   csn      = cosine of field inclination to LOS [npix]
;	eta		 = ratio of line center to continuum absorption [npix]
;	hgen	 = voigt functions [maxsz,3,npix]
;	fgen	 = faraday-voigt functions [maxsz,3,npix]
;   maxsz   = maximum wavelength dimension

;   OUTPUTS:
;   er      = eta/rho parameters


;  size of arrays
;npix = n_elements(csn)
;ftn2 = fltarr(npix) & ftn3 = ftn2 & ftn4 = ftn2 & ftn5 = ftn2

;  Output arrays of eta's & rho's
	er	= fltarr(maxsz,5)

	sp2 = 1.-csn*csn
	ez2 = .5*eta

	for iw = ipx,jpx do begin
;  Combinations that are in both eta's rho's
		ftn2 = ez2*(hgen(iw,1)-.5*(hgen(iw,0)+hgen(iw,2)))
		ftn3 = ez2*(hgen(iw,2)-hgen(iw,0))
		ftn4 = ez2*(fgen(iw,1)-.5*(fgen(iw,0)+fgen(iw,2)))
		ftn5 = ez2*(fgen(iw,2)-fgen(iw,0))

;  Compute eta's & rho's
		er(iw,1) = ftn2*sp2
		er(iw,2) = ftn3*csn
		er(iw,3) = ftn4*sp2
		er(iw,4) = ftn5*csn
		er(iw,0) = er(iw,1)+ez2*(hgen(iw,0)+hgen(iw,2))

	endfor

	return
	end