C+ C NAME: C ExtractPositionn8 C PURPOSE: C CALLING SEQUENCE: subroutine ExtractPositionn8(ITP,iYr,Doy8,RVect) C INPUTS: C OUTPUTS: C CALLS: C SunNewcomb8, UlyssesOrbit, HOSOrbit, PlanetOrbit C ECLIPTIC_HELIOGRAPHIC C- real RVect(3),VVect(5) double precision dLngSun,dLatSun,dDisSun double precision Doy8 if (ITP .eq. 3) then ! Earth call SunNewcomb8(0,iYr,Doy8,dLngSun,dLatSun,dDisSun) RVect(1) = dLngSun+180 RVect(2) = -dLatSun/3600 ! arcsec -> degrees RVect(3) = dDisSun else if (ITP .eq. 10) then ! Ulysses Doy = sngl(Doy8) call UlyssesOrbit(iYr,Doy,RVect,VVect) else if (ITP .eq. 11 .or. ! Helios 1 & ITP .eq. 12) then ! Helios 2 Doy = sngl(Doy8) call HOSOrbit(ITP-10,iYr,Doy,RVect(3),RVect(1),VR,VT) RVect(2) = 0 else Doy = sngl(Doy8) call PlanetOrbit(ITP,iYr,Doy,RVect,VVect) end if ! Convert to heliographic coordinates call ECLIPTIC_HELIOGRAPHIC8(0,iYr,Doy8,RVect(1),RVect(2)) return end