- switched from fitting fracred to x = ln(f/(1-f)). Output still in f.
- fixed error in dra calculation in split_ra_dec (was miscalculating average dec)
-added ability to pass ra, dec ranges to fit_ra_dec
- added quality parameter output to fit_ra_dec (lnP, n)
- added test routine for fitting resolution variations
- added check on existence of savefile for fit_ra_dec -- generates temporary file rather than overwrite (fool me once....)