C+
C NAME:
C	EARTH_WLIMB
C PURPOSE:
C	Returns the heliographic or ecliptic longitude (in degrees) of the
C	solar west limb
C CALLING SEQUENCE:
	function EARTH_WLIMB(iYrIn,Doy)
C INPUTS:
C	iYr		integer		year
C					(if iYr<0 then the ecliptic longitude is returned)
C	Doy		real		day of year, including fraction for time of day
C OUTPUTS:
C	EARTH_WLIMB	real		heliographic longitude of east limb (degrees, 0<=EARTH<360)
C CALLS:
C	SunNewcomb, ECLIPTIC_HELIOGRAPHIC
C SEE ALSO:
C	XMAP_SC_POS
C PROCEDURE:
C	Typically used as external function to e.g. XMAP_SC_POS
C MODIFICATION HISTORY:
C	???-????, Paul Hick (UCSD/CASS; pphick@ucsd.edu)
C-
	    integer	iYrIn
	    real	Doy

	double precision	rLng
	double precision	rLat
	double precision	rDist

	iYr = abs(iYrIn)

	call SunNewcomb(0,iYr,Doy,rLng,rLat,rDist)
	XLNG = amod(sngl(rLng)+270.,360.)

	if (iYrIn .ge. 0) then
	    XLAT = 0.
	    call ECLIPTIC_HELIOGRAPHIC(0,iYr,Doy,XLNG,XLAT)
	end if
	EARTH_WLIMB = XLNG

	return
	end