;+
; NAME:
;	ThomsonPDistanceFar
; PURPOSE:
;	Find the distance to the plane of the sky where a single electron
;	would have a given polarization in the limit of large electron-Sun
;	distance
; CALLING SEQUENCE:
	FUNCTION ThomsonPDistanceFar, ElSun, P
; INPUTS:
;	ElSun		array	distance Sun-Electron (in solar radii)
;	P			array	polarization	(-1<=P<=1)
; OUTPUTS:
;	ThomsonPDistance
;				array	distance to the plane of the sky (in solar radii) (>=0)
; INCLUDE:
	@compile_opt.pro		; On error, return to caller
; CALLS:
;	SyncArgs
; PROCEDURE:
; >	The distance to the plane of the sky is given as a positive number.
;	There are two locations on either side of the plane of the sky
;	which match the polarization.
; >	An electron can only produce a positive polarization. If a negative
;	polarization is specified, then D=-1. is returned
; >	An electron can only produce a polarization below a certain maximum
;	(less than one). If a polarization above the maximum value is specified
;	the D=-1. is returned.
; MODIFICATION HISTORY:
;	JUL-1996, Paul Hick (UCSD)
;-

D = 0.
SyncArgs, ElSun, P, D

I = where(P LT 0 OR P GT 1)
IF I[0] NE -1 THEN D[I] = BadValue(D)

I = where(0 LE P AND P LE 1)
IF I[0] ne -1 THEN D[I] = ElSun[I]*sqrt((1-P[I])/(1+P[I]))

RETURN, D  &  END