This simple example/test program computes the basic astronomical
statistics for the Sun at 00h TDB on any given date according to the
The VSOP87 C series.

NOTES:

• This demo version automatically adjusts for Julian or Gregorian
  calendar dates.
  
  
• Dates up to 1582 Oct 04 = Julian Calendar Dates.
  
  
• Dates from 1582 Oct 15 onwards = Gregorian Calendar Dates.
  
  
• BC years are represented here by using negative year values.
  However, there is no calendar year 0 (zero), by definition.
  
  

---

Enter Numerical Calendar Date Elements  (yNum  mNum  dNum)

Initial default startup date = Current date (2026 02 19) at 00:00:00 TT.

The program does not check for errors and simply assumes the calendar
date elements (yNum, mNum, dNum) are valid values.
Valid date elements should return a valid result.

If possible, check the test computations against another VSOP87 C
source to see how closely they match.

EXAMPLE:
October 5th, 1991 date elements would be entered as:   1991 10 05

------------------------------------------------------
VSOP87C Geocentric Astronomical Statistics For the Sun

For Ecliptic and Equinox of Date

yyyy mm dd  =  2026 02 19
JDate  00h  =  2461090.5

------------------------------------
Statistics returned by the function:

VSOP87C_Geocentric_Sun (JDTT)

RAHrs      = 22.161005175083 = Equatorial RA of the Sun in hours
DeclDeg    = -11.351397483625 = Equatorial declination of the Sun in degrees
DistAU     = 0.98851102733468 = True Earth/Sun distance in AUs
AngDiamDeg = 0.53933109160723 = Angular diameter of the Sun in degrees
EclLngDeg  = 330.34120460981 = Ecliptical longitude of the sun in degrees
EclLatDeg  = 4.351740622034E-5 = Ecliptical latitude of the sun in degrees