home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 January / usenetsourcesnewsgroupsinfomagicjanuary1994.iso / sources / misc / volume14 / ephem4.13 / patch02 next >
Encoding:
Text File  |  1990-08-03  |  1.2 KB  |  45 lines
  1. Newsgroups: comp.sources.misc
  2. subject: v14i045: patch to ephem, v4.13
  3. From: downey@dimed1.UUCP
  4. Sender: allbery@uunet.UU.NET (Brandon S. Allbery - comp.sources.misc)
  5.  
  6. Posting-number: Volume 14, Issue 45
  7. Submitted-by: downey@dimed1.UUCP
  8. Archive-name: ephem4.13/patch02
  9.  
  10. Here is a patch to ephem, my interactive astronomical ephemeris, v4.13.
  11. It is a very simple change to one file, objx.c.  It fixes a problem that
  12. occurs in certain circumstances computing the location of objects specified
  13. via parabolic orbital elements.o
  14.  
  15. thank you very much.
  16. Elwood Downey
  17. downey@dimed.com
  18.  
  19. *** objx.13    Thu Aug  2 18:56:52 1990
  20. --- objx.new    Thu Aug  2 18:58:54 1990
  21. ***************
  22. *** 241,248
  23.   
  24.           sll = sin(ll);
  25.           cll = cos(ll);
  26. !         *lam = atan(rsn*sll/(rpd-rsn*cll))+lpd;
  27. !         /* *lam = atan(-1*rpd*sll/(rsn-rpd*cll))+lg+PI; */
  28.           range (lam, 2*PI);
  29.           *bet = atan(rpd*spsi*sin(*lam-lpd)/(cpsi*rsn*sll));
  30.   
  31.  
  32. --- 241,250 -----
  33.   
  34.           sll = sin(ll);
  35.           cll = cos(ll);
  36. !         if (rpd < rsn)
  37. !         *lam = atan(-1*rpd*sll/(rsn-rpd*cll))+lg+PI;
  38. !         else
  39. !         *lam = atan(rsn*sll/(rpd-rsn*cll))+lpd;
  40.           range (lam, 2*PI);
  41.           *bet = atan(rpd*spsi*sin(*lam-lpd)/(cpsi*rsn*sll));
  42.   
  43.  
  44.  
  45.