C+
C NAME:
C	ThomsonTang
C PURPOSE:
C	(internal use only) Tangential intensity in 1/r^2 density
C CALLING SEQUENCE:
	function ThomsonTang(S)
C CALLS:
C	ElSunDistance, ThomsonBase
C INCLUDE
	include		'sun.h'
C-
	    real	S

	parameter	(denAU = 5.)
	save		ScSun, Elo, U

	ElSun = ElSunDistance(ScSun,S,Elo, SinChi)
	dI  = ThomsonBase(ElSun,SinChi,U, P,dIt,dItr)

	Den = denAU/(ElSun*SUN__RAU)**2		! # el/cm^3 for 1/r^2 density

	ThomsonTang = Den*dIt

	return
C+
C NAME:
C	ThomsonTangMRad
C PURPOSE:
C	(internal use only) Tangential minus radial intensity in 1/r^2 density
C CALLING SEQUENCE:
	entry ThomsonTangMRad(S)
C CALLS:
C	ElSunDistance, ThomsonBase
C-
	ElSun = ElSunDistance(ScSun,S,Elo, SinChi)
	dI  = ThomsonBase(ElSun,SinChi,U, P,dIt,dItr)

	Den = denAU/(ElSun*SUN__RAU)**2		! # el/cm^3 for 1/r^2 density

	ThomsonTangMRad = Den*dItr

	return
C+
C NAME:
C	ThomsonTang3D
C PURPOSE:
C	(internal use only)
C CALLING SEQUENCE:
	entry ThomsonTang3D(S)
C CALLS:
C	ElSunDistance, ThomsonBase, ThomsonLOSDensity
C-
	ElSun = ElSunDistance(ScSun,S,Elo, SinChi)
	dI  = ThomsonBase(ElSun,SinChi,U, P,dIt,dItr)

	Den = ThomsonLOSDensity(S)		! # el/cm^3 for 1/r^2 density

	ThomsonTang = Den*dIt

	return
C+
C NAME:
C	ThomsonTangMRad3D
C PURPOSE:
C	(internal use only) Tangential minus radial intensity
C CALLING SEQUENCE:
	entry ThomsonTangMRad3D(S)
C CALLS:
C	ElSunDistance, ThomsonBase, ThomsonLOSDensity
C-
	ElSun = ElSunDistance(ScSun,S,Elo, SinChi)
	dI  = ThomsonBase(ElSun,SinChi,U, P,dIt,dItr)

	Den = ThomsonLOSDensity(S)		! # el/cm^3 for 1/r^2 density

	ThomsonTangMRad = Den*dItr

	return
C+
C NAME:
C	ThomsonSetupIntegrand
C PURPOSE:
C	(internal use only)
C CALLING SEQUENCE:
	entry ThomsonSetupIntegrand(ScSunIn,EloIn,UIn)
C-
	ThomsonSetupIntegrand = 0.
	ScSun = ScSunIn
	Elo   = EloIn
	U     = UIn

	return
	end