if (!strcmp(argv[ap], "-1")) {
			ConstantFit(data, ndata, 1.0);
		if (!strcmp(argv[ap], "-A")) {
			fprintf(stderr, "Computing Modified Allan Variance into %s\n", argv[ap+1]);
		if (!strcmp(argv[ap], "-a")) {
			fprintf(stderr, "Computing Allan Variance into %s\n", argv[ap+1]);
		if (!strcmp(argv[ap], "-c")) {
			ConstantFit(data, ndata, 0.0);
		if (!strcmp(argv[ap], "-C")) {
			ConstantFit(data, ndata, strtod(argv[ap+1], 0));
		}
		/* Fit a frequency model of the form y = A * ln(t + 1) */
		if (!strcmp(argv[ap], "-e")) {
			ExpFit(data, ndata);
		if (!strcmp(argv[ap], "-F")) {
			if (ap + 2 >= argc)
				errx(1,"Missing arguments for -F");
			FFT(data, ndata, atoi(argv[ap + 1]), argv[ap + 2]);
		if (!strcmp(argv[ap], "-J")) {
			AutoJump(data, &ndata, to);
		if (!strcmp(argv[ap], "-j")) {
			Jump(data, &ndata, fm, to);
		}
		if (!strcmp(argv[ap], "-l")) {
			fprintf(stderr, "Computing Least Square Fit\n");
		 Correct for a linear frequency drift.

		 The argument is who many samples to span for each
		 entry into the average sum.  The model only
		 converges of this value is chosen near the flicker floor.
		if (!strcmp(argv[ap], "-L")) {
			fm = strtoul(argv[ap+1], 0, 0);
			fprintf(stderr, "Computing Linear Frequency Drift, using %d segments of %d samples\n", fm, ndata / fm);
			LinFreq(data, ndata, ndata / fm);
			ap+=2;
			continue;
		}
		/* Omega, the period of the input signal.  For instance if you measure two
		 * 5 Mhz signals, this is 200nsec.  This is used to bridge the gaps in an
		 * input sequence.
		 */
		if (!strcmp(argv[ap], "-O")) {
			if (ap + 1 == argc)
				errx(1,"Missing argument for -O");
			omega = strtod(argv[ap+1], 0);
			ap+=2;
			continue;
		}
		if (!strcmp(argv[ap], "-o")) {
			fm = strtoul(argv[ap+1], 0, 0);
			OutLyers(data, ndata, fm);
			ap += 2;
			continue;
		}
		if (!strcmp(argv[ap], "-p")) {
			fprintf(stderr, "Computing Parabolic Fit\n");
			ParabolicFit(data, ndata);
			ap++;
			continue;
		}
		if (!strcmp(argv[ap], "-r")) {
			if (ap + 1 == argc)
				errx(1,"Missing argument for -r");
			fprintf(stderr, "Reading TS data from %s\n", argv[ap+1]);
			TsReadData(argv[ap+1], &data, &ndata);
			ap += 2;
			continue;
		}
		if (!strcmp(argv[ap], "-R")) {
			fprintf(stderr, "Reading data from %s\n", argv[ap+1]);
		}
		if (!strcmp(argv[ap], "-S")) {
			StraightLine(data, ndata, a);
		}
		if (!strcmp(argv[ap], "-s")) {
			StraightFit(data, ndata);
		}
		if (!strcmp(argv[ap], "-t")) {
			tau = strtod(argv[ap+1], 0);
		}
		if (!strcmp(argv[ap], "-w")) {
			WriteData(argv[ap+1], data, ndata, to);