C+
C NAME:
C	ThomsonMidFar
C CALLING SEQUENCE:
	function ThomsonMidFar(S)		! Internal use only
C INCLUDE:
	include		'math.h'		! Input S in units of Observer-Sun distance
C-
	    real	S

	!-------
	! Entry point ThomsonSetupMidFar

	real		EloIn

	!-------

	save		E

	sE  = sin(E)
	cE  = cos(E)
	Y   = S-cE				! Distance to plane of sky
	Y2  = Y*Y
	sE2 = sE*sE
	YE  = Y2+sE2
	ThomsonMidFar = Y*(Y2+3*sE2/5)/(YE*YE)+cE*(1-2*sE2/5)/2+(atan2(Y,sE)-E/2)/sE

	return
C+
C NAME:
C	ThomsonSetupMidFar
C CALLING SEQUENCE:
	entry ThomsonSetupMidFar(EloIn)		! Internal use only
C PROCEDURE:
C	Entry point in href=ThomsonSetupMidFar=
C-
	ThomsonSetupMidFar = 1.
	E = EloIn*MATH__RPD			! Convert to radians

	return
	end