C+
C NAME:
C	Thomson3DDensity
C PURPOSE:
C	Implements 1/r^2 density, for use in ThomsonLOSDensity
C CATEGORY:
C CALLING SEQUENCE:
	function Thomson3DDensity(ElLng,ElLat,ElSun)
C INPUTS:
C	ElLng			real	heliocentric (ecliptic) longitude (deg)
C	ElLat			real	heliocentric (ecliptic) latitude (deg)
C	ElSun			real	heliocentric distance (solar radii)
C OUTPUTS:
C	Thomson3DDensity	real	density (electrons/cm^3)
C CALLS:
C INCLUDE:
	include		'sun.h'
C PROCEDURE:
C >	At 1 AU a density of 5 cm^-3 is assumed.
C >	Called by ThomsonLOSDensity, which in turn is called by the
C	ThomsonLOS3D* functions to  calculate line of sight intensities.
C >	Can be replaced by any function calculating the density in a
C	heliocentric coordinate system.
C MODIFICATION HISTORY:
C	JAN-1998, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
C-
	    real	ElLng
	    real	ElLat
	    real	ElSun

	parameter	(denAU = 5.)

	Thomson3DDensity = ElSun*SUN__RAU	! El-Sun distance in AU
	Thomson3DDensity = denAU/(Thomson3DDensity*Thomson3DDensity)

	return
	end