home *** CD-ROM | disk | FTP | other *** search
/ Usenet 1994 October / usenetsourcesnewsgroupsinfomagicoctober1994disk2.iso / misc / volume12 / ephem4.13 / patch01 next >
Text File  |  1990-04-25  |  1KB  |  31 lines
  1. Newsgroups: comp.sources.misc
  2. subject: v12i015: ephem phase correction
  3. From: downey@cs.umn.edu@dimed1.UUCP
  4. Sender: allbery@uunet.UU.NET (Brandon S. Allbery - comp.sources.misc)
  5.  
  6. Posting-number: Volume 12, Issue 15
  7. Submitted-by: downey@cs.umn.edu@dimed1.UUCP
  8. Archive-name: ephem4.13/patch01
  9.  
  10. This is a correction to ephem related to calculation of phase (percent
  11. illumination). All posted versions of ephem to date (up to and including 4.13)
  12. calculate phase correctly only for objects in the ecliptic. Note that the
  13. problem is particularly evident with comets.
  14.  
  15. The fix is just one line in circum.c. Change the line that reads:
  16.         f = 0.5*(1+cos(lam-lpd0));
  17. to be:
  18.         f = 0.25 * (((rp0+rho0)*(rp0+rho0) - rsn*rsn)/(rp0*rho0));
  19.  
  20. It is a simple enough change that I did not make this an automatic patch of
  21. any sort. The exact line number changes with different versions but the line
  22. is unique and unchanged in all releases.
  23.  
  24. Sorry for the problem. Please contact me if you have any questions.
  25.  
  26. Elwood Downey
  27. downey%dimed1.uucp@cs.umn.edu, or
  28. uunet!rutgers!umn-cs!dimed1!downey
  29.  
  30.  
  31.