added nstar routine for simultaneous PSF fitting.  Not yet well-tested