/* MCHMAG - Change magnitude zero points in a matched catalog. */

#include	"focas1.h"
#include	"match1.h"


mchmag (mcatalog)
char   mcatalog[];
{
    int     i;
    long    pos;
    float   dm[MAXMATCH];
    char    str[20];

    cfd = mcatopen (mcatalog, 2);
    for (i = 0; i < ncats; i++) {
	printf ("Catalog %2d - ", i);
	pos = ftell (cfd);
	rdcathdr (cfd, 1);
	dm[i] = sp.magfst;
	printf ("Magnitude zero point (%6.2f): ", dm[i]);
	if (strlen (gets (str)))
	    sscanf (str, "%hf", &dm[i]);
	dm[i] = dm[i] - sp.magfst;
	sp.magfst += dm[i];
	fseek (cfd, pos, 0);
	wtcathdr (cfd, 1);
    }

    for (;;) {
	pos = ftell (cfd);
	if (rdmcatob (cfd, &mentry, mob))
	    break;
	for (i = 0; i < ncats; i++)
	    if (mob[i].entnum >= 0)
		mob[i].mag += dm[i];
	fseek (cfd, pos, 0);
	wtmcat (cfd, &mentry, mob);
    }

    fclose (cfd);
}