#!/usr/bin/env python3

from baseband import vdif
import argparse, sys, math
import numpy as np

parser = argparse.ArgumentParser()
parser.add_argument('-n', '--nframe', '-nfra,e', help="Number of frames to process", type=int, default=100)
#parser.add_argument('-p', '--plot', '-plot', help="Plot results", action="store_true")
#parser.add_argument('-r', '--robust', '-robust', help="Robust Fit", action="store_true")
parser.add_argument('vdiffile')

args = parser.parse_args()

np.set_printoptions(precision=3)

f = vdif.open(args.vdiffile, 'rb')

for i in range(args.nframe):
    frame = f.read_frame()

    data = frame.data

    print("{:3d}: ".format(i), "Mean=", data.mean(axis=0), "  StdDev=",data.std(axis=0))


f.close()