message: fix gslib bug for InterpolateGeneral when ncomp>1 and some points were not found
diff: @@ -837,7 +837,7 @@ void FindPointsGSLIB::InterpolateGeneral(const GridFunction &field_in,
  sarray_transfer(struct send_pt, sendpt, proc, 1, cr);
  sdpt = (struct send_pt *)sendpt->ptr;
  - for (int index = 0; index < nptorig; index++)
  + for (int index = 0; index < sendpt->n; index++)
  {
  int idx = sdpt->index + j*nptorig;
  field_out(idx) = sdpt->ival;
  