c..common blocks for the ephemeris routines
      integer          iemax
      parameter        (iemax = 20000)

      character*40     zname(iemax)
      common /zcxa/    zname

      integer          nplanets,nmarsat,njupsat,nsatsat,nurasat,
     1                 ncomet,nnaster,
     2                 nsum1,nsum2,nsum3,nsum4,nsum5,nsum6
      common /zcxb/    nplanets,nmarsat,njupsat,nsatsat,nurasat,
     1                 ncomet,nnaster,
     2                 nsum1,nsum2,nsum3,nsum4,nsum5,nsum6

      double precision begjd,endjd
      common /zcxc/    begjd,endjd




c..data structure for the comet orbital elements
c..comet_name     = name of comet
c..comet_epoch    = modified julian date of the elements
c..comet_q        = perihelion distance in au
c..comet_e        = eccentricity of the orbit
c..comet_i        = inclination of the orbit to j2000 ecliptic plane in degrees
c..comet_w        = argument of perihelion to j2000 ecliptic plane in degrees
c..comet_node     = longitude of the ascending node to j2000 in degrees
c..comet_tp_year  = year of perhelion passage
c..comet_tp_month = month of perhelion passage
c..comet_tp_day   = day of perhelion passage
c..comet_tp_hour  = hour of perhelion passage
c..comet_jpl      = orbital solution reference

      integer          ncomet_max
      parameter        (ncomet_max = 400)
      character*38     comet_name(ncomet_max)
      character*10     comet_jpl(ncomet_max)
      integer          comet_epoch(ncomet_max),
     1                 comet_tp_year(ncomet_max),
     2                 comet_tp_month(ncomet_max),
     3                 comet_tp_day(ncomet_max)
      double precision comet_q(ncomet_max),comet_e(ncomet_max),
     1                 comet_i(ncomet_max),comet_w(ncomet_max),
     2                 comet_node(ncomet_max),comet_tp_hour(ncomet_max)

      common /cmeta1/  comet_name,comet_jpl
      common /cmetb1/  comet_epoch,comet_tp_year,comet_tp_month,
     1                 comet_tp_day
      common /cmetc1/  comet_q,comet_e,comet_i,comet_w,
     1                 comet_node,comet_tp_hour





c..data structure for the numbered astroid orbital elements
c..naster_name     = name of numbered asteroid 
c..naster_epoch    = modified julian date of the elements
c..naster_a        = seme-major axis in au
c..naster_e        = eccentricity of the orbit
c..naster_i        = inclination of orbit to j2000 ecliptic plane in degrees
c..naster_w        = argument of perihelion to j2000 ecliptic plane in degrees
c..naster_node     = longitude of the ascending node to j2000 in degrees
c..naster_m        = mean anomaly in degrees
c..naster_jpl      = orbital solution reference

      integer          nnaster_max
      parameter        (nnaster_max = 15000)
      character*38     naster_name(nnaster_max)
      character*10     naster_jpl(nnaster_max)
      integer          naster_epoch(nnaster_max)
      double precision naster_a(nnaster_max),naster_e(nnaster_max),
     1                 naster_i(nnaster_max),naster_w(nnaster_max),
     2                 naster_node(nnaster_max),naster_m(nnaster_max)


      common /casta1/  naster_name,naster_jpl
      common /castb1/  naster_epoch
      common /castc1/  naster_a,naster_e,naster_i,naster_w,
     1                 naster_node,naster_m