float *mapres(MapBeam *mb, Model *mod, float *clnmp, float bmaj, float bmin,
  float bpa, int dosub, int noresid, int dosmth);