subroutine ECLIPTIC_EQUATOR8(ID,iYr,Doy8,PHI,RLAT)
	
            integer     ID
            integer     iYr
            real        PHI
            real        RLAT
	    
	    real*8      Doy8
            real*8      EPS
            real*8      JD
            real*8      JEpoch

            call Julian8(0,iYr,Doy8,JD,JEpoch)			!Get date in Julian years
            EPS = 23.439291d0-.00013004d0*(JEpoch-2000.d0)	!Angle (ecliptic-equator)

            BETA = sngl(EPS)                              	! .. to equatorial coord
            if (ID .eq. 1) BETA = -BETA				! .. to ecliptic coord
            call Rotate(90.,BETA,-90.,PHI,RLAT)
        return
        end