  /* incorporate unmatched image stars -- after all catalogs searched... */
  for (i = 0; i < Nstars; i++) {
    if (!found[i]) {
      catalog[0].average[Nave].R = stars[N1[i]].R;
      catalog[0].average[Nave].D = stars[N1[i]].D;
      catalog[0].average[Nave].M = 1000.0*stars[N1[i]].M;
      catalog[0].average[Nave].Nm = 1;
      catalog[0].average[Nave].Xp = catalog[0].average[Nave].Xm = 0.0;
      catalog[0].average[Nave].offset = Nmeas;
      catalog[0].measure[Nmeas].dR  = 0.0;
      catalog[0].measure[Nmeas].dD  = 0.0;
      catalog[0].measure[Nmeas].M   = 1000.0*stars[N1[i]].M;
      catalog[0].measure[Nmeas].dM  = stars[N1[i]].dM;
      catalog[0].measure[Nmeas].t  = 3600.0 + stars[N1[i]].Y*0.1;
      catalog[0].measure[Nmeas].average = (0x00ffffff & Nave);
      /* time of star = time of image + offset + Y*rate */
      Nave ++;
      if (Nave == 0x1000000) {
	fprintf (stderr, "too many stars in catalog\n");
	clear_lockfile ();
	exit (0);
      }
      if (Nave == NAVE) {
	NAVE = Nave + 1000;
	REALLOCATE (next, int, NAVE);
	REALLOCATE (catalog[0].average, Average, NAVE);
      }
      Nmeas ++;
      if (Nmeas == NMEAS) {
	NMEAS = Nmeas + 1000;
	REALLOCATE (next, int, NMEAS);
	REALLOCATE (catalog[0].measure, Measure, NMEAS);
      }
    }
  }
